PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Molecular Dynamics Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 73 records · Page 4Linked to original sources
PubMed2005

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study.

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.

Binding Sites↗
PubMed2006

Gd(III) polyaminocarboxylate chelate: realistic many-body molecular dynamics simulations for molecular imaging applications.

Realistic molecular dynamics simulations of polyaminocarboxylate complexes of gadolinium (III) ion in water are performed, providing coordination numbers and average residence times in quantitative agreement with available experimental data. A theoretical analysis, based on fitting a fluctuating charges model on ab initio data, also indicates that charge transfer between the ion and the ligand is significant.

Letter↗
PubMed1990

Using molecular dynamics simulations on crambin to evaluate the suitability of different continuum dielectric and hydrogen atom models for protein simulations.

Molecular dynamics simulations of enzymes with enough explicit waters of solvation to realistically account for solute-solvent interactions can burden the computational resources required to perform the simulation by more than two orders of magnitude. Since enzyme simulations even with an implicit solvation model can be imposing for a supercomputer, it is important to assess the suitability of different continuum dielectric models for protein simulations. A series of 100-picosecond molecular dynamics simulations were performed on the X-ray crystal structure of the protein crambin to examine how well computed structures, obtained using seven continuum dielectric and two hydrogen atom models, agreed with the X-ray structure. The best level of agreement between computed and experimental structures was obtained using a constant dielectric of 2 and the all-hydrogen model. Continuum dielectric models of 1, 1r, and 2r also led to computed structures in reasonably good agreement with the X-ray structure. In all cases, the all-hydrogen model gave better agreement than the united atom model, although, in one case, the difference was not significant. Dielectric models of 4, 80, and 4r with either hydrogen model yielded significantly poorer fits. It is especially noteworthy that the observed trends did not semiquantitatively converge until about 50 picoseconds into the simulations, suggesting that validation studies for protein calculations based on energy minimizations or short simulations should be viewed with caution.

Computer Simulation↗
PubMed1995

Fluctuation and cross-correlation analysis of protein motions observed in nanosecond molecular dynamics simulations.

Nanosecond molecular dynamics simulations of bovine pancreatic trypsin inhibitor and lysozyme in water are analyzed in terms of backbone atomic positional fluctuations and dynamical cross-correlations. It is found that although the molecular systems are stable, B-factors calculated over a time period as long as 500 ps are not representative for the motions within the proteins. This is especially true for the most mobile residues. On a nanosecond time-scale, the B-factors calculated from the simulations of the proteins in solution are considerably larger than those obtained by structure refinement of the proteins in crystals, based on X-ray data. The time evolution of the atomic fluctuations shows that for large portions of the proteins under study, atomic positional fluctuations are not yet converged after a nanosecond. Cross-correlations do not converge faster than the fluctuations themselves. Most display very erratic behavior if the sampling covers less than about 200 ps. It is also shown that inclusion of mobile atoms into the procedure used to remove rigid-body motion from the simulation can lead to spurious correlations between the motions of the atoms at the surface of the protein.

Animals↗
PubMed2005

What contributions to protein side-chain dynamics are probed by NMR experiments? A molecular dynamics simulation analysis.

Molecular dynamics simulations of the structurally homologous proteins TNfn3 and FNfn10 have been used to investigate the contributions to side-chain dynamics measured by NMR relaxation experiments. The results reproduce the variation in core side-chain dynamics observed by NMR and highlight the relevance of anharmonic motion and transitions between local minima for explaining NMR side-chain order parameters. A method is described for calculating converged order parameters by use of replica exchange molecular dynamics in conjunction with an implicit solvent model. These simulations allow the influence of various factors, such as the flexibility of side-chains and their free volume, on the mobility to be tested by perturbing the system. Deletion mutations are found to have the largest effect on the more densely packed FNfn10. Some counterintuitive effects are seen, such as an increase in order parameters close to deletion mutation sites, but these can be rationalized in terms of direct interactions with the modified side-chains. A statistical analysis of published order parameters supports the conclusions drawn from the simulations.

Amino Acids↗
PubMed2004

Diffusion of hard sphere fluids in disordered media: a molecular dynamics simulation study.

Molecular dynamic simulations are reported for the static and dynamic properties of hard sphere fluids in matrices (or media) composed of quenched hard spheres. The effect of fluid and matrix density, matrix structure, and fluid to matrix sphere size ratio on the static and dynamic properties is studied using discontinuous molecular dynamics. The matrix density has a stronger effect on the self-diffusion coefficient than the fluid density, especially at high matrix densities where the geometric constraints due to the quenched spheres are significant. When the ratio of the size of the fluid spheres to that of the matrix spheres is equal to or greater than one, the diffusion increases as the fluid density is increased, at constant total volume fraction. This trend is however reversed if the ratio is smaller than one. Different methods of generating the matrix have a very strong effect on the dynamic properties even though the static correlations are similar. An analysis of the single-chain structure factor of the particle trajectories shows a change in the particle diffusive behavior at different time scales, suggestive of a hopping mechanism, although normal diffusion is recovered at long times. At high matrix densities, there is considerable heterogeneity in the diffusion of the fluid particles. The simulations demonstrate that the correlations in the matrix play a significant role on the diffusion of fluid spheres. For example, the diffusion constant in matrices constructed by different methods can be an order of magnitude different even though the pair correlation functions are almost identical.

Journal Article↗
PubMed2003

Association lifetimes of hydrophobic amino acid pairs measured directly from molecular dynamics simulations.

Unrestrained molecular dynamics simulations have been used to estimate the lifetimes and equilibrium constants for association of six different pairs of hydrophobic amino acids in aqueous solution. The results indicate that current computational resources are sufficient that the thermodynamics of association of small molecules in water can be routinely measured, thus demonstrating that the development of molecular force fields to reproduce aqueous phase binding data is now an achievable goal.

Computer Simulation↗
PubMed2003

Deformation-mechanism map for nanocrystalline metals by molecular-dynamics simulation.

Molecular-dynamics simulations have recently been used to elucidate the transition with decreasing grain size from a dislocation-based to a grain-boundary-based deformation mechanism in nanocrystalline f.c.c. metals. This transition in the deformation mechanism results in a maximum yield strength at a grain size (the 'strongest size') that depends strongly on the stacking-fault energy, the elastic properties of the metal, and the magnitude of the applied stress. Here, by exploring the role of the stacking-fault energy in this crossover, we elucidate how the size of the extended dislocations nucleated from the grain boundaries affects the mechanical behaviour. Building on the fundamental physics of deformation as exposed by these simulations, we propose a two-dimensional stress-grain size deformation-mechanism map for the mechanical behaviour of nanocrystalline f.c.c. metals at low temperature. The map captures this transition in both the deformation mechanism and the related mechanical behaviour with decreasing grain size, as well as its dependence on the stacking-fault energy, the elastic properties of the material, and the applied stress level.

Computer Simulation↗
PubMed2003

Interactions of the designed antimicrobial peptide MB21 and truncated dermaseptin S3 with lipid bilayers: molecular-dynamics simulations.

Molecular-dynamics simulations covering 30 ns of both a natural and a synthetic antimicrobial peptide in the presence of a zwitterionic lipid bilayer were performed. In both simulations, copies of the peptides were placed in an alpha-helical conformation on either side of the bilayer about 10 A (1 A=0.1 nm) from the interface, with either the hydrophobic or the positively charged face of the helix directed toward the bilayer surface. The degree of peptide-lipid interaction was dependent on the starting configuration: surface binding and subsequent penetration of the bilayer was observed for the hydrophobically oriented peptides, while the charge-oriented peptides demonstrated at most partial surface binding. Aromatic residues near the N-termini of the peptides appear to play an important role in driving peptide-lipid interactions. A correlation between the extent of peptide-lipid interactions and helical stability was observed in the simulations. Insertion of the peptides into the bilayer caused a dramatic increase in the lateral area per lipid and decrease in the bilayer thickness, resulting in substantial disordering of the lipid chains. Results from the simulations are consistent with early stages of proposed mechanisms for the lytic activity of antimicrobial peptides. In addition to these 'free' simulations, 25 ns simulations were carried out with the peptides constrained at three different distances relative to the bilayer interface. The constraint forces are in agreement with the extent of peptide-bilayer insertion observed in the free simulations.

Amino Acid Sequence↗
PubMed1999

Charge pairing of headgroups in phosphatidylcholine membranes: A molecular dynamics simulation study.

Molecular dynamics simulation of the hydrated dimyristoylphosphatidylcholine (DMPC) bilayer membrane in the liquid-crystalline phase was carried out for 5 ns to study the interaction among DMPC headgroups in the membrane/water interface region. The phosphatidylcholine headgroup contains a positively charged choline group and negatively charged phosphate and carbonyl groups, although it is a neutral molecule as a whole. Our previous study (Pasenkiewicz-Gierula, M., Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi. 1997. J. Phys. Chem. 101:3677-3691) showed the formation of water cross-bridges between negatively charged groups in which a water molecule is simultaneously hydrogen bonded to two DMPC molecules. Water bridges link 76% of DMPC molecules in the membrane. In the present study we show that relatively stable charge associations (charge pairs) are formed between the positively and negatively charged groups of two DMPC molecules. Charge pairs link 93% of DMPC molecules in the membrane. Water bridges and charge pairs together form an extended network of interactions among DMPC headgroups linking 98% of all membrane phospholipids. The average lifetimes of DMPC-DMPC associations via charge pairs, water bridges and both, are at least 730, 1400, and over 1500 ps, respectively. However, these associations are dynamic states and they break and re-form several times during their lifetime.

Binding Sites↗
PubMed2000

Structure of the Sm binding site from human U4 snRNA derived from a 3 ns PME molecular dynamics simulation.

A molecular dynamics simulation of the Sm binding site from human U4 snRNA was undertaken to determine the conformational flexibility of this region and to identify RNA conformations that were important for binding of the Sm proteins. The RNA was fully-solvated (>9,000 water molecules) and charge neutralized by inclusion of potassium ions. A three nanosecond MD simulation was conducted using AMBER with long-range electrostatic forces considered using the particle mesh Ewald summation method. The initial model of the Sm binding site region had the central and 3' stem-loops that flanked the Sm site co-axial with one another, and with the single-stranded Sm binding site region ([I] conformation). During the course of the trajectory, the axes of the 3' stem-loop, and later the central stem-loop, became roughly orthogonal from their original anti-parallel orientation. As these conformational changes occurred, the snRNA adopted first an [L] conformation, and finally a [U] conformation. The [U] conformation was more stable than either the [I] or [L] conformations, and persisted for the final 1 ns of the trajectory. Analysis of the structure resulting from the MD simulations revealed the bulged nucleotide, U114, and the mismatched Ag91-G110 base pair provided distinctive structural features that may enhance Sm protein binding. Based on the results of the MD simulation and the available experimental data, we proposed a mechanism for the binding of the Sm protein sub-complexes to the snRNA. In this model, the D1/D2 and E/F/G Sm protein sub-complexes first bind the snRNA in the [U] conformation, followed by conformational re-arrangement to the [I] conformation and binding of the D3/B Sm protein sub-complex.

Base Sequence↗
PubMed1990

Dynamic properties of the first enzymatic reaction steps of porcine pancreatic elastase. How rigid is the active site of the native enzyme? Molecular dynamics simulation.

Two molecular dynamics simulations (100 and 50 ps) of native porcine pancreatic elastase i.e., without bound substrate and with the active site hydrated by a dome of water (630 molecules) have been performed. Dynamical properties of the catalytic tetrad have been examined. While relative conformations of the Asp 102, His 57, and Ser 214 are rather stable in time, the side chain of Ser 195 undergoes several conformational changes. No preferences are observed for the formation of a hydrogen bond between the O gamma-H group (Ser 195) and nitrogen N, (His 57). A cluster of ordered water molecules effectively competes with the H-O gamma group (Ser 195) and thereby prevents the formation of this H bond, which is generally agreed to be crucial for catalysis.

Animals↗
PubMed2002

Modern protein force fields behave comparably in molecular dynamics simulations.

Several molecular dynamics simulations were performed on three proteins--bovine apo-calbindin D9K, human interleukin-4 R88Q mutant, and domain IIA of bacillus subtilis glucose permease--with each of the AMBER94, CHARMM22, and OPLS-AA force fields as implemented in CHARMM. Structural and dynamic properties such as solvent-accessible surface area, radius of gyration, deviation from their respective experimental structures, secondary structure, and backbone order parameters are obtained from each of the 2-ns simulations for the purpose of comparing the protein portions of these force fields. For one of the proteins, the interleukin-4 mutant, two independent simulations were performed using the CHARMM22 force field to gauge the sensitivity of some of these properties to the specific trajectory. In general, the force fields tested performed remarkably similarly with differences on the order of those found for the two independent trajectories of interleukin-4 with CHARMM22. When all three proteins are considered together, no force field showed any consistent trend in variations for most of the properties monitored in the study.

Animals↗
PubMed2004

Phase behavior in model homopolymer/CO2 and surfactant/CO2 systems: discontinuous molecular dynamics simulations.

Discontinuous molecular dynamics simulations are performed on surfactant (HmTn)/solvent systems modeled as a mixture of single-sphere solvent molecules and freely jointed surfactant chains composed of m slightly solvent-philic head spheres (H) and n solvent-philic tail spheres (T), all of the same size. We use a square-well potential to account for the head-head, head-solvent, tail-tail, and tail-solvent interactions and a hard-sphere potential for the head-tail and solvent-solvent interactions. We first simulate homopolymer/supercritical CO2 (scCO2) systems to establish the appropriate interaction parameters for a surfactant/scCO2 system. Next, we simulate surfactant/scCO2 systems and explore the effect of the surfactant volume fraction, packing fraction, and temperature on the phase behavior. The transition from the two-phase region to the one-phase region is located by monitoring the contrast structure factor of the equilibrated surfactant/scCO2 system, and the micelle to unimer transition is located by monitoring the aggregate size distribution of the equilibrated surfactant/scCO2 system. We find a two-phase region, a micelle phase, and a unimer phase with increasing packing fraction at fixed temperature or with increasing temperature at fixed packing fraction. The phase diagram for the surfactant/scCO2 system in the surfactant volume fraction-packing fraction plane and the density dependence of the critical micelle concentration are in qualitative agreement with experimental observations. The phase behavior of a surfactant/scCO2 system can be directly related to the solubilities of the corresponding homopolymers that serve as the head and tail blocks for the surfactant. The influence of surfactant structure (head and tail lengths) on the phase transitions is explored.

Journal Article↗
PubMed2005

Parametric studies of interaction strengths in polymer/CO2 systems: discontinuous molecular dynamics simulations.

Discontinuous molecular dynamics simulations are performed on homopolymer/solvent and surfactant/solvent systems. The homopolymer and surfactant molecules are modeled as freely jointed square-well chains. Solvent molecules are modeled as both hard spheres and square-well spheres. We explore how the various interaction parameters affect the types of phase behavior and micellization observed in the homopolymer/solvent and surfactant/solvent systems. Increasing the packing fraction of homopolymers in both hard-sphere solvents and square-well solvents increases the solvent's ability to dissolve homopolymers only when the segment-solvent interaction strength exceeds a critical value. Although only upper critical solution temperature (UCST) behavior is observed for homopolymers in hard-sphere solvents, both UCST and lower critical solution temperature (LCST) behavior are observed for homopolymers in square-well solvents, depending upon the interaction strengths and chain length. This indicates that it is necessary to account for the solvent-solvent attraction to model LCST behavior in supercritical CO2. Our simulation results on surfactants in hard-sphere solvents show that it is necessary to account for the interactions experienced by both the head and tail blocks in order to capture the essential features of surfactant/supercritical CO2 systems.

Journal Article↗
PubMed2006

Free energy of a trans-membrane pore calculated from atomistic molecular dynamics simulations.

Atomistic molecular dynamics simulations of a lipid bilayer were performed to calculate the free energy of a trans-membrane pore as a function of its radius. The free energy was calculated as a function of a reaction coordinate using a potential of mean constraint force. The pore radius was then calculated from the reaction coordinate using Monte Carlo particle insertions. The main characteristics of the free energy that comes out of the simulations are a quadratic shape for a radius less than about 0.3 nm, a linear shape for larger radii than this, and a rather abrupt change without local minima or maxima between the two regions. In the outer region, a line tension can be calculated, which is consistent with the experimentally measured values. Further, this line tension can be rationalized and understood in terms of the energetic cost for deforming a part of the lipid bilayer into a hydrophilic pore. The region with small radii can be described and understood in terms of statistical mechanics of density fluctuations. In the region of crossover between a quadratic and linear free energy there was some hysteresis associated with filling and evacuation of the pore with water. The metastable prepore state hypothesized to interpret the experiments was not observed in this region.

Computer Simulation↗
PubMed2005

Analyzing the mechanisms of selectivity in biomimetic self-assemblies via IR and NMR spectroscopy of prepolymerization solutions and molecular dynamics simulations.

Molecularly imprinted polymers (MIPs) for 2,4-dichlorophenoxyacetic acid were synthesized via a noncovalent approach with 4-vinylpyridine as functional monomer and ethylene glycol dimethacrylate as cross-linker in a methanol/water mixture. Templated polymers synthesized in this self-assembly approach rely on complex formation between the target analyte and functional monomers in porogenic solution prior to radical polymerization. Consequently, the achievable selectivity is governed by the nature and stability of these complexes. The nature of noncovalent interactions responsible for complex formation during imprinting of the template 2,4-dichlorophenoxyacetic acid (2,4-D) with the functional monomer 4-vinylpyridine has been investigated. Fourier transform infrared and 1H NMR spectroscopies provide the fundamental analytical basis for rationalizing the mechanisms of recognition during the imprinting process probing the governing interactions for selective binding site formation at a molecular level. Molecular modeling studies in explicit solvent (chloroform and water) corroborate the importance of hydrogen bonding in aprotic solvents and of hydrophobic interactions in protic media in agreement with the experimental spectroscopic investigations of prepolymerization solutions. Furthermore, chromatographic studies of the synthesized MIPs provided insight on the importance of size, shape, and functionality during selective 2,4-D rebinding processes confirming the results obtained during the prepolymerization studies.

Journal Article↗