Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55Linked to original sources
PubMed2005

Membrane protein structure quality in molecular dynamics simulation.

Our goal was to assess the relationship between membrane protein quality, output from protein quality checkers and output from molecular dynamics (MD) simulations. Membrane transport proteins are essential for a wide range of cellular processes. Structural features of integral membrane proteins are still under-explored due to experimental limitations in structure determination. Computational techniques can be used to exploit biochemical and medium resolution structural data, as well as sequence homology to known structures, and enable us to explore the structure-function relationships in several transmembrane proteins. The quality of the models produced is vitally important to obtain reliable predictions. An examination of the relationship between model stability in molecular dynamics (MD) simulations derived from RMSD (root mean squared deviation) and structure quality assessment from various protein quality checkers was undertaken. The results were compared to membrane protein structures, solved at various resolution, by either X-ray or electron diffraction techniques. The checking programs could predict the potential success of MD in making functional conclusions. MD stability was shown to be a good indicator for the quality of structures. The quality was also shown to be dependent on the resolution at which the structures were determined.

Computer Simulation↗
PubMed2004

Sensitivity of molecular dynamics simulations to the choice of the X-ray structure used to model an enzymatic reaction.

A subject of great practical importance that has not received much attention is the question of the sensitivity of molecular dynamics simulations to the initial X-ray structure used to set up the calculation. We have found two cases in which seemingly similar structures lead to quite different results, and in this article we present a detailed analysis of these cases. The first case is acyl-CoA dehydrogenase, and the chief difference of the two structures is attributed to a slight shift in a backbone carbonyl that causes a key residue (the proton-abstracting base) to be in a bad conformation for reaction. The second case is xylose isomerase, and the chief difference of the two structures appears to be the ligand sphere of a Mg2+ metal cofactor that plays an active role in catalysis.

Acyl-CoA Dehydrogenase↗
PubMed1997

Molecular dynamics simulations of six different fully hydrated monomeric conformers of Escherichia coli re-lipopolysaccharide in the presence and absence of Ca2+.

Six previously published conformational models of Escherichia coli Re lipopolysaccharide (ReLPS) were subjected to molecular dynamics simulations using the CHARMM force field. The monomers of ReLPS were completely immersed in a water box. The dynamic behavior of the solvated models in the presence and absence of calcium cations was compared. The structure of the solvent shell was analyzed in terms of radial distribution functions. Diffusion coefficients and mean residence times were analyzed to characterize the dynamic behavior of the solvent. Order parameters and number of gauche defects were used for the description of the dynamics of the acyl chains. The cations are preferentially located between the carboxylate and phosphate groups of the headgroup. Their presence leads to a rigidification of the headgroup structure and alters the conformation of the backbone, thus influencing the structure and flexibility of the hydrophobic region as well. The effect of calcium on the backbone flexibility was measured in terms of glycosidic torsion angles. The six fatty acid chains of each ReLPS monomer adopt a highly ordered micromembrane structure. The packing parameter indicates that aggregation of these ReLPS monomers will lead to lamellar structures. Evaluation of all data enables us to present one conformation, C, which is thought to best represent the average structure of the ReLPS conformers.

Calcium↗
PubMed1999

Study of the stability and unfolding mechanism of BBA1 by molecular dynamics simulations at different temperatures.

BBA1 is a designed protein that has only 23 residues. It is the smallest protein without disulfide bridges that has a well-defined tertiary structure in solution. We have performed unfolding molecular dynamics simulations on BBA1 and some of its mutants at 300, 330, 360, and 400 K to study their kinetic stability as well as the unfolding mechanism of BBA1. It was shown that the unfolding simulations can provide insights into the forces that stabilize the protein. Packing, hydrophobic interactions, and a salt bridge between Asp12 and Lys16 were found to be important to the protein's stability. The unfolding of BBA1 goes through two major steps: (1) disruption of the hydrophobic core and (2) unfolding of the helix. The beta-hairpin remains stable in the unfolding because of the high stability of the type II' turn connecting the two beta-strands.

Amino Acid Sequence↗
PubMed2005

Molecular dynamics simulation of the aggregation of the core-recognition motif of the islet amyloid polypeptide in explicit water.

The formation of amyloid fibrils is associated with major human diseases. Nevertheless, the molecular mechanism that directs the nucleation of these fibrils is not fully understood. Here, we used molecular dynamics simulations to study the initial self-assembly stages of the NH2-NFGAIL-COOH peptide, the core-recognition motif of the type II diabetes associated islet amyloid polypeptide. The simulations were performed using multiple replicas of the monomers in explicit water, in a confined box starting from a random distribution of the peptides at T = 300 K and T = 340 K. At both temperatures the formation of unique clusters was observed after a few nanoseconds. Structural analysis of the clusters clearly suggested the formation of "flat" ellipsoid-shaped clusters through a preferred locally parallel alignment of the peptides. The unique assembly is facilitated by a preference for an extended conformation of the peptides and by intermolecular aromatic interactions. Taken together, our results may provide a description of the molecular recognition determinants involved in fibril formation, in terms of the atomic detailed structure of nascent aggregates. These observations may yield information on new ways to control this process for either materials development or drug design.

Amyloid↗
PubMed2003

Insights into the acylation mechanism of class A beta-lactamases from molecular dynamics simulations of the TEM-1 enzyme complexed with benzylpenicillin.

Herein, we present results from molecular dynamics MD simulations ( approximately 1 ns) of the TEM-1 beta-lactamase in aqueous solution. Both the free form of the enzyme and its complex with benzylpenicillin were studied. During the simulation of the free enzyme, the conformation of the Omega loop and the interresidue contacts defining the complex H-bond network in the active site were quite stable. Most interestingly, the water molecule connecting Glu166 and Ser70 does not exchange with bulk solvent, emphasizing its structural and catalytic relevance. In the presence of the substrate, Ser130, Ser235, and Arg244 directly interact with the beta-lactam carboxylate via H-bonds, whereas the Lys234 ammonium group has only an electrostatic influence. These interactions together with other specific contacts result in a very short distance ( approximately 3 A) between the attacking hydroxyl group of Ser70 and the beta-lactam ring carbonyl group, which is a favorable orientation for nucleophilic attack. Our simulations also gave insight into the possible pathways for proton abstraction from the Ser70 hydroxyl group. We propose that either the Glu166 carboxylate-Wat1 or the substrate carboxylate-Ser130 moieties could abstract a proton from the nucleophilic Ser70.

Acylation↗
PubMed2003

Pores formed by the nicotinic receptor m2delta Peptide: a molecular dynamics simulation study.

The M2delta peptide self-assembles to form a pentameric bundle of transmembrane alpha-helices that is a model of the pore-lining region of the nicotinic acetylcholine receptor. Long (>15 ns) molecular dynamics simulations of a model of the M2delta(5) bundle in a POPC bilayer have been used to explore the conformational dynamics of the channel assembly. On the timescale of the simulation, the bundle remains relatively stable, with the polar pore-lining side chains remaining exposed to the lumen of the channel. Fluctuations at the helix termini, and in the helix curvature, result in closing/opening transitions at both mouths of the channel, on a timescale of approximately 10 ns. On average, water within the pore lumen diffuses approximately 4x more slowly than water outside the channel. Examination of pore water trajectories reveals both single-file and path-crossing regimes to occur at different times within the simulation.

Amino Acid Sequence↗
PubMed1996

Molecular dynamics simulations of adipocyte lipid-binding protein: effect of electrostatics and acyl chain unsaturation.

Molecular dynamics (MD) simulations have been performed on adipocyte lipid-binding protein, using the apo and holo forms, bound with stearic and oleic acid. The contribution of electrostatics to protein dynamics and ligand stabilization was assayed by perturbing the electrostatic charge of Arg106 and Arg126 (positive-->neutral) and the fatty acid (132H) headgroup (negative-->neutral). MD simulations for charged holo forms demonstrated significantly greater electrostatic binding energy and a more stabilized hydrogen bond network than simulations performed using neutral forms. Electrostatics, however, appeared to have little effect on fatty acid behavior, e.g., fluctuation of the dihedral head group; number of dihedral transitions within the acyl chain; and change in the end-to-end distance for fatty acid. Instead, fatty acid behavior appeared to be dictated by the presence or absence of an unsaturated bond within the acyl chain. A significantly greater number of transitions were observed during MD simulations in oleic than stearic acid. In addition, significantly greater fluctuation was observed for oleic acid, within the C2 headgroup and C9 and C11 dihedrals (which lie adjacent to the olefin bond of oleic acid). The dynamic behavior of the acyl chain may thereby be more a property of van der Waals contact, and the degree of acyl chain unsaturation, than a function of electrostatics. In the absence of fatty acid, an increase in distance between guanidino carbon centered atoms of Arg126 and Arg106 was observed during MD simulations of the charged apo form. This effect not observed with the neutral apo form or in any of the holo complexes and, presumbably, was a result of repulsion between the negatively charged arginine sidechains. Conserved waters reflected substantially lower mean-square displacement (msd) in all simulations, except the neutral apo form. This suggests that the presence of either charged amino acids or lipid provides increased order for water within the binding pocket. These results provide a dynamic perspective of the interactive nature within the FABP binding pocket regulated in a complex manner by the electrostatics within the binding cavity, acyl chain structure and behavior, and water energetics.

Carrier Proteins↗
PubMed2002

Molecular dynamics simulations of conformational behavior of linear RGD peptidomimetics and cyclic prodrugs in aqueous and octane solutions.

Conformations available to a class of cyclic prodrugs and corresponding linear RGD peptidomimetics were explored using 1 ns length molecular dynamics simulations performed with the program CHARMM. Water and octane, modeled explicitly, were used as solvents to mimic the change of the environment experienced by the solutes upon partition from water to membrane in the trans-cellular transport process. In water, the linear peptidomimetics tended to populate extended-like structures, characterized by strong favorable interactions with solvent and low intrinsic stability. In these extended conformations the charged termini are able to assume large distances, above 15 A for the longest systems. These linear peptidomimetics have been found to exhibit the highest potency in experimental studies, in accord with the trends experimentally observed for RGD peptides. In contrast, in octane compact conformers of the linear peptidomimetics were favored, with all charged groups aggregated and shielded from solvent, exhibiting high intrinsic stability and weak solute-solvent interactions. Our calculations predict a large unfavorable energy change for transferring the linear systems from water to octane, in agreement with experimental findings that these compounds are not transported via the trans-cellular pathway. The cyclic prodrugs did not exhibit major structural differences between the simulations in water and octane, adopting turn-like conformations in both solvents. The limited response of the cyclic structures to changes in the environment leads to energies of transfer from water to octane that are also unfavorable, but much less so than for the linear molecules. This effect is in accord with the observed enhanced passive trans-cellular transport of the cyclic prodrugs.

Amino Acid Sequence↗
PubMed2006

Layering at an ionic liquid-vapor interface: a molecular dynamics simulation study of [bmim][PF6].

The structure of the planar liquid-vapor interface of a room-temperature ionic liquid, 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6]), is studied using atomistic molecular dynamics simulations. Layering of the ions at the interface is observed as oscillations in the corresponding number density profiles. These oscillations, however, are diminished in amplitude in the electron density profile, due to a near cancellation in the contributions from the anions and the cations. An enhancement by 12% in the electron density at the interface over its value in the bulk liquid is observed, in excellent agreement with X-ray reflectivity experiments. The anions are found to predominantly contribute to this increase in the interfacial electron density. The cations present at the interface are oriented anisotropically. Their butyl chains are observed to be preferentially oriented along the interface normal and to project outside the liquid surface, thus imparting a hydrophobic character. In the densest region of the interface, the imidazolium ring plane is found to lie parallel to the surface normal, in agreement with direct recoil spectroscopy experiments.

Journal Article↗
PubMed2005

Structure and vibrational spectroscopy of salt water/air interfaces: predictions from classical molecular dynamics simulations.

We report the sum frequency generation (SFG) spectra of aqueous sodium iodide interfaces computed with the methodology outlined by Morita and Hynes (J. Phys. Chem. B 2002, 106, 673), which is based on molecular dynamics simulations. The calculated spectra are in qualitative agreement with experiment. Our simulations show that the addition of sodium iodide to water leads to an increase in SFG intensity in the region of 3400 cm(-1), which is correlated with an increase in ordering of hydrogen-bonded water molecules. Depth-resolved orientational distribution functions suggest that the ion double layer orders water molecules that are approximately one water layer below the Gibbs dividing surface. We attribute the increase in SFG intensity to these ordered subsurface water molecules that are present in the aqueous sodium iodide/air interfaces but are absent in the neat water/air interface.

Air↗
PubMed2006

Structure of isolated tryptophyl-glycine dipeptide and tryptophyl-glycyl-glycine tripeptide: ab initio SCC-DFTB-D molecular dynamics simulations and high-level correlated ab initio quantum chemical calculations.

The tryptophyl-glycine (Trp-Gly) and tryptophyl-glycyl-glycine (Trp-Gly-Gly) peptides have been studied by means of molecular dynamic simulations combined with high-level correlated ab initio quantum chemical and statistical thermodynamic calculations. The lowest energy conformers were localized in the free energy surface. The structures of the different Trp-Gly and Trp-Gly-Gly conformers coexisting in the gas phase have been for the first time reported and their scaled theoretical IR spectra unambiguously assigned and compared with previous gas-phase experimental results. Common geometrical features have been systematically observed for the sequence Trp, Trp-Gly, and Trp-Gly-Gly. In addition, the peptide backbone of Trp-Gly-Gly has been compared with that of the previously studied Phe-Gly-Gly (Reha, D. et. al. Chem. Eur. J. 2005, 11, 6803). From the observed systematic structural behavior between these peptide analogues, it is expected that the gas-phase conformers of other similar aromatic small peptides would present equivalent geometries. The DFT methodology failed to describe the potential energy surface of the studied peptides since the London dispersion energy (not covered in DFT) plays a significant role in the stabilization of most stable conformers.

Oligopeptides↗
PubMed1993

Conformation of a cyclic opioid peptide analog by NMR and molecular dynamics simulation.

The conformation of a new Ty-c[Lys-Phe-Asp]-NH2 cyclic opioid peptide synthesized by solid phase method, has been determined from two-dimensional NMR and distance geometry followed by restrained molecular dynamics simulation. The conformation of the ring is well-defined, but the exocylic Tyr-1 and Phe-3 side-chain moiety possesses significant orientational freedom.

Amino Acid Sequence↗
PubMed2004

Molecular dynamics simulations of transport and separation of carbon dioxide-alkane mixtures in carbon nanopores.

The configurational-bias Monte Carlo method, which is used for efficient generation of molecular models of n-alkane chains, is combined for the first time with the dual control-volume grand-canonical molecular-dynamics simulation, which has been developed for studying transport of molecules in pores under an external potential gradient, to investigate transport and separation of binary mixtures of n-alkanes, as well as mixtures of CO2 and n-alkanes, in carbon nanopores. The effect of various factors, such as the temperature of the system, the composition of the mixture, and the pore size, on the separation of the mixtures is investigated. We also report the preliminary results of an experimental study of transport and separation of some of the same mixtures in a carbon molecular-sieve membrane with comparable pore sizes. The results indicate that, for the mixtures considered in this paper, even in very small carbon nanopores the energetic effects still play a dominant role in the transport and separation properties of the mixtures, whereas in a real membrane they are dominated by the membrane's morphological characteristics. As a result, for the mixtures considered, a single pore may be a grossly inadequate model of a real membrane, and hence one must resort to three-dimensional molecular pore network models of the membrane.

Journal Article↗
PubMed2002

Dynamic structure of He-Ne mixtures by molecular dynamics simulation: from hydrodynamic to fast and slow sound modes.

Molecular dynamics (MD) results for the dynamic structure of a He(0.77)Ne(0.23) gas mixture at two densities (15.8 and 36.1 nm(-3)) show a clear crossover from hydrodynamic modes to distinct excitations for the two species. The higher density dispension curve neatly shows high- and low-frequency branches setting on with a rather localized transition. The lower density results agree very well with existing neutron scattering data and, in particular, display hydrodynamic behavior up to k approximately 2 nm(-1), in contrast with the conclusions of previous simulation studies. A smooth transition to fast sound is shown to take place for 2<k/nm(-1)<5, where the present MD data fill the existing gap in the experimental results.

Journal Article↗
PubMed2001

Classical molecular interaction potentials: improved setup procedure in molecular dynamics simulations of proteins.

The latest version of the classical molecular interaction potential (CMIP) has the ability to predict the position of crystallographic waters in several proteins with great accuracy. This article analyzes the ability of the CMIP functional to improve the setup procedure of the molecular system in molecular dynamics (MD) simulations of proteins. To this end, the CMIP strategy is used to include both water molecules and counterions in different protein systems. The structural details of the configurations sampled from trajectories obtained using the CMIP setup procedure are compared with those obtained from trajectories derived from a standard equilibration process. The results show that standard MD simulations can lead to artifactual results, which are avoided using the CMIP setup procedure. Because the CMIP is easy to implement at a low computational cost, it can be very useful in obtaining reliable MD trajectories.

Acetylcholinesterase↗
PubMed1995

Molecular dynamics simulations of isolated helices of myoglobin.

The apo form of myoglobin has two non-native stable states that have been experimentally characterized. Investigation of these states has suggested possible folding pathways for myoglobin. We have performed molecular dynamics simulations on solvated isolated helices of myoglobin to investigate the relationship between the intrinsic stabilities of the isolated helices and the structure and folding pathway of apomyoglobin. Analyses of hydrogen bonding and fluctuations from simulations at 298 and 368 K are used to explore the relative stabilities of the helices of myoglobin. The ordering observed is A approximately G approximately H > B > E > F, which mirrors both the experimental equilibrium and kinetic data available for apomyoglobin. The experimental observation that a subdomain comprising helices A, G, and H is an important early intermediate and our result that these helices are the most stable suggest that the intrinsically more stable helices form early in the folding process and that this significantly influences the folding pathway.

Amino Acid Sequence↗
PubMed1993

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations.

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.

Computer Simulation↗