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 955 records · Page 53Linked to original sources
PubMed1999

Molecular dynamics simulation in solvent of the bacteriophage 434 cI repressor protein DNA binding domain amino acids (R1-69) in complex with its cognate operator (OR1) DNA sequence.

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the 434 cI Repressor protein DNA binding domain (DBD) amino acids (R1-69) and DNA of operator (OR1) and its flanks consisting of 28 nucleotide base pairs. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the 434 cI repressor DNA recognition helix 3 formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the OR1 DNA major groove halfsites and flanking regions. In addition, hydrophilic amino acids within the loop between helix 3 and helix 4 have strong electrostatic attraction to codon-anticodon nucleotides located within the central nucleotides of the minor groove between the OR1 major groove halfsites. These interactions together induced significant structural changes in the operator DNA manifested by overtwisting of the central nucleotide base pairs and narrowing of the minor groove between the DNA major groove halfsites. Finally, these findings offer a code for site specific DNA recognition by the 434 cI repressor protein.

Amino Acid Sequence↗
PubMed2001

Algorithm for molecular dynamics simulations of spin liquids.

A new symplectic time-reversible algorithm for numerical integration of the equations of motion in magnetic liquids is proposed. It is tested and applied to molecular dynamics simulations of a Heisenberg spin fluid. We show that the algorithm exactly conserves spin lengths and can be used with much larger time steps than those inherent in standard predictor-corrector schemes. The results obtained for time correlation functions demonstrate the evident dynamic interplay between the liquid and magnetic subsystems.

Journal Article↗
PubMed2005

Molecular-dynamics simulation for the characterization of liquid chromatographic stationary phase: effect of temperature.

The influence of temperature on the surface structure of the octadecylsilica (ODS) bonded phase was investigated with a molecular dynamics (MD) simulation. The MD simulation was applied to a molecular model consisting of three parts: amorphous silica base, dimethyloctadecylsilyl ligands and n-hexane as a mobile phase solvent. More detailed information on the effect of temperature was obtained at the low temperature region than that reported in our previous study. The motion of ODS ligands could be estimated by the mean square displacement (MSD) of the terminal carbon atoms of ODS ligands. The gauche fraction in the ODS ligand conformation can also be estimated to obtain the ligand conformation for each simulation condition in detail. It can be seen that an elevated temperature induced the more bent ligand conformation. The trend has a good agreement to that of the results experimentally observed by using various spectroscopic techniques such as nuclear magnetic resonance (NMR), Fourier-transform infrared spectroscopy (FT-IR), and Raman spectroscopy.

Journal Article↗
PubMed1996

Drug-membrane interactions studied by molecular dynamics simulation: size dependence of diffusion.

Bioavailability, largely due to membrane permeation, is an important step in the drug delivery process and therefore drug design. In this study, a series of molecular dynamics simulations (totaling 10's of nanoseconds) of small molecules (varying in size, and functional groups) in lipid bilayer membranes were used to elucidate the mechanism of diffusion of drugs within biomembranes. These simulations accurately reproduce many experimentally observed parameters. The simulations also agree with theory that indicates the lipid bilayer has internal structure that influences the diffusion process and that even within the hydrocarbon regions there are distinct regions between which the rate and mechanism of diffusion varies. In particular, in agreement with experiment and theory the small solutes were found to diffuse by a mechanism different from that of large molecules. These variations are linked to the frequency and size of spontaneously arising voids within the bilayer as well as the rate of torsional isomerization of the hydrocarbon chains.

Biological Availability↗
PubMed2003

Why is Leu55-->Pro55 transthyretin variant the most amyloidogenic: insights from molecular dynamics simulations of transthyretin monomers.

Transthyretin (TTR) is one of the known human amyloidogenic proteins. Its native state is a homotetramer with each monomer having a beta-sandwich structure. Strong experimental evidence suggests that TTR dissociates into monomeric intermediates and that the monomers subsequently self-assemble to form amyloid deposits and insoluble fibrils. However, details on the early steps along the pathway of TTR amyloid formation are unclear, although various experimental approaches with resolutions at the molecular or residue level have provided some clues. It is highly likely that the stability and flexibility of monomeric TTR play crucial roles in the early steps of amyloid formation; thereby, it is essential to characterize initial conformational changes of TTR monomers. In this article we probe the possibility that the differences in the monomeric forms of wild-type (WT) TTR and its variants are responsible for differential amyloidogenesis. We begin with the simulations of WT, Val30-->Met (V30M), and Leu55-->Pro (L55P) TTR monomers. Nanosecond time scale molecular dynamics simulations at 300 K were performed using AMBER. The results indicate that the L55P-TTR monomer undergoes substantial structural changes relative to fluctuations observed in the WT and V30M TTR monomers. The observation supports earlier speculation that the L55P mutation may lead to disruption of the beta-sheet structure through the disorder of the "edge strands" that might facilitate amyloidogenesis.

Amyloid↗
PubMed2006

Molecular dynamics simulations of polyelectrolyte-polyampholyte complexes. Effect of solvent quality and salt concentration.

Complexation between polyelectrolyte and polyampholyte chains in poor solvent conditions for the polyelectrolyte backbone has been studied by molecular dynamics simulations. In a poor solvent a polyelectrolyte forms a necklace-like structure consisting of polymeric globules (beads) connected by strings of monomers. The simulation results can be explained by assuming the existence of two different mechanisms leading to the necklace formation. In the case of weak electrostatic interactions, the necklace formation is driven by optimization of short-range monomer-monomer attraction and electrostatic repulsion between charged monomers on the polymer backbone. In the case of strong electrostatic interactions, the necklace structure appears as a result of counterion condensation. While the short-range attractions between monomers are still important, the correlation-induced attraction between condensed counterions and charged monomers and electrostatic repulsion between uncompensated charges provide significant contribution to optimization of the necklace structure. Upon forming a complex with both random and diblock polyampholytes, a polyelectrolyte chain changes its necklace conformation by forming one huge bead. The collapse of the polyelectrolyte chain occurs due to the neutralization of the polyelectrolyte charge by polyampholytes. In the case of the random polyampholyte, the more positively charged sections of the chain mix with negatively charged polyelectrolyte forming the globular bead while more negatively charged chain sections form loops surrounding the collapsed core of the aggregate. In the case of diblock polyampholyte, the positively charged block, a part of the negatively charged block, and a polyelectrolyte chain form a core of the aggregate with a substantial section of the negatively charged block sticking out from the collapsed core of the aggregate. In both cases the core of the aggregate has a layered structure that is characterized by the variations in the excess of concentration of monomers belonging to polyampholyte and polyelectrolyte chains throughout the core radius. These structures appear as a result of optimization of the net electrostatic energy of the complex and short-range attractive interactions between monomers of the polyelectrolyte chain.

Computer Simulation↗
PubMed1991

Conformation of endothelin in aqueous ethylene glycol determined by 1H-NMR and molecular dynamics simulations.

The solution conformation of a 21-residue vasoconstrictor peptide endothelin-1 (ET-1) in water-ethylene glycol has been determined by two-dimensional 1H-NMR spectroscopy and constrained molecular dynamics simulations. The N-terminus (residues 1-4) appears to undergo conformational averaging and no single structure consistent with the NMR constraints could be found for this region. Residues 5-8 form a turn, and residues 9-16 exist in a helical conformation. A flexible 'hinge' between residues 8-9 allows various orientations of the turn relative to the helix. Another 'hinge' at residue 17 connects the extended C-terminus to the bicyclic core region (residues 1-15). Residues important for binding and biological activity form a contiguous surface on one side of the helix, with the two disulfides extending from the other side of the helix.

Computer Simulation↗
PubMed1997

Molecular dynamics simulations of lipid bilayers.

During the past decade, computer simulations of bilayers have moved from the realm of model systems to realistic systems containing tens of phospholipids along with the requisite number of water molecules hydrating the entire molecular assembly. Concomitant with the ability to model larger and larger systems, simulators have also begun to utilize more accurate numerical tools to ensure that the temperature, pressure, simulation timescales, parameter sets and long-range electrostatic interactions of bilayers are correctly accounted for in a typical molecular dynamics simulation. With these tools in hand, work has already begun to define the structure, function and dynamics of bilayer, bilayer/small molecule and bilayer/protein systems. Thus, we have reached an era in which simulators will tackle more and more detailed issues regarding complex bilayer systems.

Anesthetics↗
PubMed2000

Structure prediction of the dimeric neu/ErbB-2 transmembrane domain from multi-nanosecond molecular dynamics simulations.

Dimerization of the neu/ErbB-2 receptor tyrosine kinase is a necessary but not a sufficient step for signaling. Despite the efforts expended to identify the molecular interactions responsible for receptor-receptor contacts and particularly those involving the transmembrane domain, structural details are still unknown. In this work, molecular dynamics simulations of the helical transmembrane domain (TM) of neu and ErbB-2 receptors are used to predict their dimer structure both in the wild and oncogenic forms. A global conformational search method, applied to define the best orientations of parallel helices, showed an energetically favorable configuration with the specific mutation site within the interface, common for both the nontransforming and the transforming neu/ErbB-2 TM dimers. Starting from this configuration, a total of 10 simulations, about 1.4 ns each, performed in vacuum, without any constraints, show that the two helices preferentially wrap in left-handed interactions with a packing angle at about 20 degrees. The resulting structures are nonsymmetric and the hydrogen bond network analysis shows that helices experience pi local distortions that facilitate inter-helix hydrogen bond interactions and may result in a change in the helix packing, leading to a symmetric interface. For the mutated sequences, we show that the Glu side chain interacts directly with its cognate or with carbonyl groups of the facing backbone. We show that the connectivity between interfacial residues conforms to the knobs-into-holes packing mode of transmembrane helices. The dimeric interface described in our models is discussed with respect to mutagenesis studies.

Amino Acid Sequence↗
PubMed2006

Molecular dynamics simulations of surfactant self-organization at a solid-liquid interface.

Self-organization of aqueous surfactants at a planar graphite-like surface is studied by means of coarse-grain molecular dynamics simulations. The nonionic surfactant, n-alkyl poly(ethylene oxide), and water are both represented by coarse-grain models while an implicit representation is used for the graphite surface. The observed morphology of the aggregated surfactants depends on the alkyl chain length. Surfactants with a short chain form a monolayer on the graphite surface with a thickness roughly equal to that of the alkane tail. On the other hand, longer-tail surfactants form continuous hemicylinders on the surface with diameter approximately 5.0 +/- 0.5 nm, in good agreement with experimental AFM data.

Adsorption↗
PubMed2006

Spontaneous fibril formation by polyalanines; discontinuous molecular dynamics simulations.

Fibrillary protein aggregates rich in beta-sheet structure have been implicated in the pathology of several neurodegenerative diseases. In this work, we investigate the formation of fibrils by performing discontinuous molecular dynamics simulations on systems containing 12 to 96 model Ac-KA(14)K-NH(2) peptides using our newly developed off-lattice, implicit-solvent, intermediate-resolution model, PRIME. We find that, at a low concentration, random-coil peptides assemble into alpha-helices at low temperatures. At intermediate concentrations, random-coil peptides assemble into alpha-helices at low temperatures and large beta-sheet structures at high temperatures. At high concentrations, the system forms beta-sheets over a wide range of temperatures. These assemble into fibrils above a critical temperature which decreases with concentration and exceeds the isolated peptide's folding temperature. At very high temperatures and all concentrations, the system is in a random-coil state. All of these results are in good qualitative agreement with those by Blondelle and co-workers on Ac-KA(14)K-NH(2) peptides. The fibrils observed in our simulations mimic the structural characteristics observed in experiments in terms of the number of sheets formed, the values of the intra- and intersheet separations, and the parallel peptide arrangement within each beta-sheet. Finally, we find that when the strength of the hydrophobic interaction between nonpolar side chains is high compared to the strength of hydrogen bonding, amorphous aggregates, rather than fibrillar aggregates, are formed.

Amyloid↗
PubMed2005

Molecular dynamics simulation of self- and mutual diffusion coefficients for confined mixtures.

The self- and mutual diffusion coefficients for binary mixtures of Ar-Kr both in the bulk and in the nanopores were studied by molecular dynamics simulations. The composition dependences and the relationships between the self- and the mutual diffusion coefficients both in the bulk and in the nanopores were further discussed. It was found that the simulation results (D(c.m.)) are close to the calculated ones (D(s)) for the Ar-Kr system. Both self- and mutual diffusion coefficients in nanopores are much lower than that of the bulk, and they ever decrease as the pore width decreases. Nevertheless, the self- and mutual diffusion coefficients increase as the mole fraction of Ar increases, and as expected, increase as the temperature increases. The self-diffusion coefficients of mixtures both in the bulk and in the nanopores are predicted by the Carman model and by the molecular cluster model.

Journal Article↗
PubMed2006

Mixed quantum-classical molecular dynamics simulation of vibrational relaxation of ions in an electrostatic field.

The vibrational relaxation of ions in low-density gases under the action of an electrostatic field is reproduced through a molecular dynamics simulation method. The vibration is treated though quantum mechanics and the remaining degrees of freedom are considered classical. The procedure is tested through comparison against analytic results for a two-dimensional quantum model and by studying energy exchange during binary ion-atom collisions. Finally, the method has been applied successfully to the calculation of the mobility and the vibrational relaxation rate of O2+ in Kr as a function of the mean collision energy using a model interaction potential that reproduces the potential minimum of a previously known ab initio potential surface. The calculation of the steady mean vibrational motion of the ions in (flow) drift tubes seems straightforward, though at the expense of large amounts of computer time.

Journal Article↗
PubMed2001

Molecular dynamics simulations of spin and pure liquids with preservation of all the conservation laws.

A methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the time step, all the conservation laws inherent in the description without breaking the time reversibility. As a result, an implicit second-order algorithm is derived and applied to pure liquids, as well as spin liquids, for which the dynamics is characterized by the conservation of total energy, linear and angular momenta, as well as magnetization and individual spin lengths. It is demonstrated on the basis of Lennard-Jones and Heisenberg fluid models that when such quantities as energy and magnetization must be conserved perfectly, the algorithm turns out to be more efficient than popular decomposition integrators and standard predictor-corrector schemes.

Journal Article↗
PubMed2006

Molecular dynamics simulation suggests possible interaction patterns at early steps of beta2-microglobulin aggregation.

Early events in aggregation of proteins are not easily accessible by experiments. In this work, we perform a 5-ns molecular dynamics simulation of an ensemble of 27 copies of beta(2)-microglobulin in explicit solvent. During the simulation, the formation of intermolecular contacts is observed. The simulation highlights the importance of apical residues and, in particular, of those at the N-terminus end of the molecule. The most frequently found pattern of interaction involves a head-to-head contact arrangement of molecules. Hydrophobic contacts appear to be important for the establishment of long-lived (on the simulation timescale) contacts. Although early events on the pathway to aggregation and fibril formation are not directly related to the end-state of the process, which is reached on a much longer timescale, simulation results are consistent with experimental data and in general with a parallel arrangement of intermolecular beta-strand pairs.

Computer Simulation↗
PubMed2005

Molecular-dynamics simulations of C- and N-terminal peptide derivatives of GCN4-p1 in aqueous solution.

We report the investigation of two 16-residue peptides in aqueous solution by means of molecular-dynamics simulations. The peptides constitute the C- and N-terminal halves of the 33-residue monomer whose dimer constitutes the leucine zipper of the yeast transcriptional activator, denoted GCN4-p1. To examine a hypothesis about coiled-coil formation, in which the C-terminal half contains a helix-formation trigger site absent in the N-terminal half, experimental studies of the two peptides have determined their helix propensities under several conditions of temperature, pH, and salt concentration with circular dichroism. An NMR experiment provides additional evidence. At temperatures of 278 and 325 K and pH 7.5, mixtures of alpha- and pi-helical secondary structure constitute the most probable conformations in both C- and N-terminal halves. A bifurcated salt bridge between Arg25 and Glu22/20 correlates with the structural fluctuations of the C-terminal half. It also exhibits a persistent loop at the N-terminal end involving the side chains of His18 and Glu22, which is reminiscent of helix-capping boxes. Nonreversible unfolding appears to occur abruptly in the Arg25 mutant, suggesting a cooperative event. Analysis does not indicate that the N-terminal half is less stable than the C-terminal half, indicating that 100 ns is too short a period to observe complete unfolding.

Amino Acid Sequence↗
PubMed1994

Molecular dynamics simulations of thrombin inhibitors at the S' subsites of thrombin.

Molecular dynamics were carried out to simulate binding interactions between S' subsites of thrombin and hirudin-based thrombin inhibitors. These inhibitors include three segments: an active-site segment, N alpha-acetyl-(D-Phe)-Pro-Arg-Pro-; a fibrinogen-recognition exo-site segment, hirudin 55-65; and a 13-atom-long linker. These linkers have been reported (Szewczuk et al. (1993) Biochemistry 32, 3396) to influence the binding potency while keeping the same active and exo-site segments. The study found that, by combining different omega-amino acids, the potency could be increased 8-fold or decreased 4-fold compared to the native hirudin linker, -Gln-Ser-His-Asn-Asp-Gly-. Five typical linkers were simulated and compared. Analyzing the trajectory files led to the classification of three different dynamic behaviours for the linkers. The flexible linkers had no influence on the antithrombin activity. Other linkers formed hydrogen bonds with the thrombin S' subsite residues Glu39, leu40, and Gln 151. Formations of some hydrogen bonds enhanced the potency of the inhibitor. In other cases, the hydrogen bonds caused the distortion of the inhibitor conformation while affected the binding potency. Based on these observations, a general binding mode in the S' subsites of thrombin is proposed and potential applications are discussed.

Amino Acid Sequence↗
PubMed2005

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations.

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.

1,2-Dipalmitoylphosphatidylcholine↗