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 487 records · Page 27Linked to original sources

The dynamics of protein hydration water: a quantitative comparison of molecular dynamics simulations and neutron-scattering experiments.

We present results from an extensive molecular dynamics simulation study of water hydrating the protein Ribonuclease A, at a series of temperatures in cluster, crystal, and powder environments. The dynamics of protein hydration water appear to be very similar in crystal and powder environments at moderate to high hydration levels. Thus, we contend that experiments performed on powder samples are appropriate for discussing hydration water dynamics in native protein environments. Our analysis reveals that simulations performed on cluster models consisting of proteins surrounded by a finite water shell with free boundaries are not appropriate for the study of the solvent dynamics. Detailed comparison to available x-ray diffraction and inelastic neutron-scattering data shows that current generation force fields are capable of accurately reproducing the structural and dynamical observables. On the time scale of tens of picoseconds, at room temperature and high hydration, significant water translational diffusion and rotational motion occur. At low hydration, the water molecules are translationally confined but display appreciable rotational motion. Below the protein dynamical transition temperature, both translational and rotational motions of the water molecules are essentially arrested. Taken together, these results suggest that water translational motion is necessary for the structural relaxation that permits anharmonic and diffusive motions in proteins. Furthermore, it appears that the exchange of protein-water hydrogen bonds by water rotational/librational motion is not sufficient to permit protein structural relaxation. Rather, the complete exchange of protein-bound water molecules by translational displacement seems to be required.

Computer Simulation↗

Molecular dynamics simulations of Cu,Zn superoxide dismutase: effect of temperature on dimer asymmetry.

Molecular dynamics simulations of solvated dimeric Cu,Zn superoxide dismutase have been carried out at four temperatures, namely 200, 225, 250 and 300 K. Analysis of the backbone-to-backbone hydrogen bonds number indicates that the symmetry observed in the two subunits at 200 K is gradually lost by heating the system. The C(alpha) atoms displacement cross-correlation maps confirm that the asymmetric behaviour of the two subunits increases as a function of temperature. The dynamic cross-correlation of the subunits volumes indicates a fast correlation between the two subunits at 300 K, which is delayed upon lowering the simulation temperature. These results indicate that temperature plays an essential role in injecting such an asymmetry; the two subunits being asymmetric and in rapid communication at 300 K, and almost symmetric and in slow communication at lower temperatures.

Journal Article↗

Molecular dynamic simulations of self-assembled alkylthiolate monolayers on an Au(III) surface.

Molecular dynamics simulations incorporating explicit gold atoms in the simulations have been carried out for alkanethiol self-assembled monolayers chemisorbed on the Au(III) surface. The structural properties of the monolayer are evaluated for two force fields: one in which the Au--S--C bond is fixed (FF I), and the other in which it is flexible (FF II). The influence of these force fields on the structural properties of HS(CH2)14CH3 on the structured Au surface is compared at different temperatures. FF I yields greater tilt angles and a smaller film thickness when compared with FF II. Both of the force fields predict that the tilt angles do not follow a monotonic decrease with temperature but show minima around 200 K. Simulations carried out at different chain lengths at 300 K reveal that FF II predicts a greater film thickness than FF I; however, the difference is within 1 A.

Journal Article↗

Molecular dynamics simulations of amphiphilic graft copolymer molecules at a water/air interface.

Fully atomistic molecular dynamics simulations of amphiphilic graft copolymer molecules have been performed at a range of surface concentrations at a water/air interface. These simulations are compared to experimental results from a corresponding system over a similar range of surface concentrations. Neutron reflectivity data calculated from the simulation trajectories agrees well with experimentally acquired profiles. In particular, excellent agreement in neutron reflectivity is found for lower surface concentration simulations. A simulation of a poly(ethylene oxide) (PEO) chain in aqueous solution has also been performed. This simulation allows the conformational behavior of the free PEO chain and those tethered to the interface in the previous simulations to be compared.

Journal Article↗

Molecular dynamics simulation study of hydrogen bonding in aqueous poly(ethylene oxide) solutions.

A molecular dynamics simulation study of hydrogen bonding in poly(ethylene oxide) (PEO)/water solutions was performed. PEO-water and water-water hydrogen bonding manifested complex dependence on both composition and temperature. Strong water clustering in concentrated solutions was seen. Saturation of hydrogen bonding at w(p) approximately equal to 0.5 and a dramatic decrease in PEO-water hydrogen bonding with increasing temperature, consistent with experimentally observed closed-loop phase behavior, were observed. Little tendency toward intermolecular bridging of PEO chains by water molecules was seen.

Journal Article↗

Molecular dynamics simulation study on behaviors of liquid 1,2-dichloroethane under external electric fields.

Molecular dynamics simulation was carried out to study the behavior of liquid 1,2-dichloroethane molecules under external electric fields including direct current field, alternating current field and positive-half-period cosin field. The maximum applied field strength was 10(8) V/m, the maximum frequency of the alternating current field and that of the positive-half-period cosine field was 10(12) Hz. The simulation revealed that the field type and field strength act on the population of the molecular configuration. In the strong direct current field, all trans forms converted completely into gauche forms. Order parameter and the correlation of the system torsion angle were also investigated. The results suggested that these two dynamical parameters depended also on the field type and the field strength. The maximum of order parameter was found to be at 0.6 in the strong direct current field.

Journal Article↗

Diffusive dynamics of protein folding studied by molecular dynamics simulations of an off-lattice model.

We report the results of a molecular dynamics study on the kinetic properties of a small off-lattice model of proteins. The model consists of a linear chain of monomers interacting via a number of potentials. These include hydrophobic, bond-angle, and torsion potentials. The ground-state conformation of the studied model is a beta-sheet motif. Molecular dynamics simulations focused on the time evolution of the reaction coordinate measuring the similarity of a given conformation with the native state. Folding time for the studied model is calculated following the diffusive-rate formula of Bryngelson and Wolynes [J. Phys. Chem. 93, 6902 (1989)] by using a computed separately configurational diffusion coefficient. Comparison of the folding time with the mean-first passage time obtained directly from folding simulations shows that the approximation depicting the dynamics of the reaction coordinate in protein folding as a diffusive motion on a free-energy landscape is quantitatively correct for the studied model.

Chemical Phenomena↗

Molecular dynamics simulations of a lipid bilayer and of hexadecane: an investigation of membrane fluidity.

Molecular dynamics simulations of a fluid-phase dipalmitoyl phosphatidylcholine lipid bilayer in water and of neat hexadecane are reported and compared with nuclear magnetic resonance spin-lattice relaxation and quasi-elastic neutron scattering data. On the 100-picosecond time scale of the present simulations, there is effectively no difference in the reorientational dynamics of the carbons in the membrane interior and in pure hexadecane. Given that the calculated fast reorientational correlation times and the "microscopic" lateral diffusion of the lipids show excellent agreement with the experimental results, it is concluded that the apparently high viscosity of the membrane is more closely related to molecular interactions on the surface rather than in the interior.

1,2-Dipalmitoylphosphatidylcholine↗

Molecular dynamics simulation of [Gd(egta)(H(2)O)](-) in aqueous solution: internal motions of the poly(amino carboxylate) and water ligands, and rotational correlation times.

Molecular dynamics simulations of [Gd(egta)(H(2)O)](-) (egta(4-)=3,12-bis(carboxymethyl)-6,9-dioxa-3,12-diazatetradecanedioate(4-)) have been performed without any artificial constraint on the first coordination sphere, such as covalent bonds between the Gd(3+) and the coordination sites. Two new crystallographic structures were determined for this gadolinium chelate and used to start two molecular dynamics simulations. [Gd(egta)(H(2)O)](-) and [Gd(egta)](-) were both observed during the simulations, with a mean volume for the reaction of dissociation [Gd(egta)(H(2)O)](-)-->[Gd(egta)](-)+H(2)O of +7.2 cm(3)mol(-1), which corroborates the previously published experimental value of +10.5 cm(3)mol(-1). Changes in the conformation of the complex with the inversion of several dihedral angles are observed in the simulations independently from the water dissociation. Very fast changes of the third-order rotation axis direction of the Gd(3+) coordination polyhedron (of symmetry D(3h)) are observed during the simulations and are related to the mechanism of electronic relaxation of the complex. Different rotational correlation times (tau(R)) were calculated from the simulations on various observables of the complex. Protons of the inner sphere have different tau(R). The mean tau(R) of the two Gd-HW(HW=hydrogen of water molecule) vectors is 72% lower than tau(R) of the complex, and 75% lower than tau(R) of the vector Gd-OW (OW=oxygen of water molecule). This discrimination of the tumbling rates should be taken into account in future global (17)O NMR, EPR and NMRD (nuclear magnetic relaxation dispersion) data analysis.

Journal Article↗

Modeling anticancer drug-DNA interactions via mixed QM/MM molecular dynamics simulations.

The development of anticancer drugs started over four decades ago, with the serendipitous discovery of the antitumor activity of cisplatin and its successful use in the treatment of various cancer types. Despite the efforts made in unraveling the mechanism of the action of cisplatin, as well as in the rational design of new anticancer compounds, in many cases detailed structural and mechanistic information is still lacking. Many of these drugs exert their anticancer activity by covalently binding to DNA inducing a distortion or simply impeding replication, thus triggering a cellular response, which eventually leads to cell death. A detailed understanding of the structural and electronic properties of drug-DNA complexes and their mechanism of binding is the key step in elucidating the principles of their anticancer activity. At the theoretical level, the description of covalent drug-DNA complexes requires the use of state-of-the-art computer simulation techniques such as hybrid quantum/classical molecular dynamics simulations. In this review we provide a general overview on: drugs which covalently bind to DNA duplexes, the basic concepts of quantum mechanics/molecular mechanics (QM/MM), molecular dynamics methods and a list of selected applications of these simulations to the study of drug-DNA adducts. Finally, the potential and the limitations of this approach to the study of such systems are critically evaluated.

Antineoplastic Agents↗

Molecular dynamics simulation approaches to K channels: conformational flexibility and physiological function.

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

Animals↗

Molecular dynamics simulations of p-tert-butylcalix[4]arene with small guest molecules.

Classical molecular dynamics simulations were used to study p-tert-butylcalix[4]arene inclusion compounds with xenon, nitrogen, hydrogen, methane, and sulfur dioxide guest molecules. The calixarene units were taken to be rigid and the intermolecular molecular interactions were modeled as a sum of the van der Waals interactions with parameters from the AMBER force field and electrostatic interactions. Simulations of the high-density alpha phase and low-density beta0 phase of p-tert-butylcalix[4]arene were used to test the force field. The predicted densities of the two phases were found to agree with experimental measurements at 173 K to within 5 %. Simulations were performed with guests placed inside the calixarene cages of the beta0 phase. Guest-host ratios of 1:1 to 1:4 were considered. Changes in the unit-cell volume and density of the phases with the addition of guest molecules and the inclusion energies for the guests were determined. Finally, the dynamics of the guest motion inside the cages were characterized by determining the root-mean-square displacements and velocity autocorrelation functions of the xenon and nitrogen guests.

Journal Article↗

Membrane structure of the human immunodeficiency virus gp41 fusion peptide by molecular dynamics simulation. II. The glycine mutants.

In this work, molecular dynamics (MD) simulation of the interaction of three mutants, G3V, G5V and G10V, of the human immunodeficiency virus (HIV) gp41 16-residue fusion peptide (FP) with an explicit palmitoyloleoylphosphatidyl-ethanolamine (POPE) lipid bilayer was performed. The goals of this work are to study the correlation of the fusogenic activity of the FPs with the mode of their interaction with the bilayer and to examine the roles of the many glycine residues in the FP in the fusion process. The results of this work corroborate the main conclusion of our earlier MD work of the WT FP and several mutants with polar substitution. These two studies provide correlation between the mode of insertion and the fusogenic activity of these peptides and support the hypothesis that an oblique insertion of the fusion domain of the viral protein is required for fusogenic activity. Inactive mutants interact with the bilayer by a surface-binding mode. The results of this work, combined with the results of our earlier work, show that, while the secondary structures of the wild-type FP and its mutants do not affect the fusogenic activities, the conformational flexibility appears to be an important factor. The active WT FP and its partially active mutants, G3V and G5V, all have significant conformational transitions at one of the glycine sites. They occur at Gly(5) in FP-wt, at Gly(10) in FP-G5V and at Gly(13) in FP-G3V. Thus, a glycine site in each of these active (or partially active) FPs provides conformational flexibility. On the other hand, the inactive mutants FP-G10V, FP-L9R and FP-V2E do not have any conformational transitions except at either terminus and thus possess no conformational flexibility. Thus, the results of this work support the suggestion that the role of glycine residues in the fusion domain is to provide the necessary conformational flexibility for fusion activity. The glycines also form a "glycine strip" in the FP that locates on one (the less hydrophobic) face of the helix (the "sided helix"). However, whether this "glycine strip" is disrupted or not does not seem to correlate with the retention of fusogenic activities. Finally, although the FLGFL (8-12) motif is absolutely conserved in the HIV fusion domain, a well-structured motif stabilized by hydrogen bonding does not appear to be required for activity. In fact, hydrogen bonding in this motif was found to be missing in FP-G3V and FP-G5V. Both of these mutants are partially active.

Glycine↗

Molecular dynamics simulations of the lipid bilayer edge.

Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.

Anisotropy↗

Orientation and conformation of a lipase at an interface studied by molecular dynamics simulations.

Electron density profiles calculated from molecular dynamics trajectories are used to deduce the orientation and conformation of Thermomyces lanuginosa lipase and a mutant adsorbed at an air-water interface. It is demonstrated that the profiles display distinct fine structures, which uniquely characterize enzyme orientation and conformation. The density profiles are, on the nanosecond timescale, determined by the average enzyme conformation. We outline a computational scheme that from a single molecular dynamics trajectory allows for extraction of electron density profiles referring to different orientations of the lipase relative to an implicit interface. Profiles calculated for the inactive and active conformations of the lipase are compared with experimental electron density profiles measured by x-ray reflectivity for the lipase adsorbed at an air-water interface. The experimental profiles contain less fine structural information than the calculated profiles because the resolution of the experiment is limited by the intrinsic surface roughness of water. Least squares fits of the calculated profiles to the experimental profiles provide areas per adsorbed enzyme and suggest that Thermomyces lanuginosa lipase adsorbs to the air-water interface in a semiopen conformation with the lid oriented away from the interface.

Ascomycota↗

Selectivity analysis of 5-(arylthio)-2,4-diaminoquinazolines as inhibitors of Candida albicans dihydrofolate reductase by molecular dynamics simulations.

A series of 5-(arylthio)-2,4-diaminoquinazolines are known as selective inhibitors of dihydrofolate reductase (DHFR) from Candida albicans. We have performed docking and molecular dynamics simulations of these inhibitors with C. albicans and human DHFR to understand the basis for selectivity of these agents. Study was performed on a selected set of 10 compounds with variation in structure and activity. Molecular dynamics simulations were performed at 300 K for 45 ps with equilibration for 10 ps. Trajectory data was analyzed on the basis of hydrogen bond interactions, energy of binding and conformational energy difference. The results indicate that hydrogen bonds formed between the compound and the active site residues are responsible for inhibition and higher potency. The selectivity index, i.e the ratio of I50 against human DHFR to I50 against fungal DHFR, is mainly determined by the conformation adapted by the compounds within the active site of two enzymes. Since the human DHFR active site is rigid, the compound is trapped in a higher energy conformation. This energy difference between the two conformations deltaE mainly governs the selectivity against fungal DHFR. The information generated from this analysis of potency and selectivity should be useful for further work in the area of antifungal research.

Antifungal Agents↗

Role of cholesterol and polyunsaturated chains in lipid-protein interactions: molecular dynamics simulation of rhodopsin in a realistic membrane environment.

We present a 118-ns molecular dynamics simulation of rhodopsin embedded in a bilayer composed of a 2:2:1 mixture of 1-stearoyl-2-docosahexaenoyl-phosphatidylcholine (SDPC), 1-stearoyl-2-docosahexaenoyl-phosphatidylethanolamine (SDPE), and cholesterol, respectively. The simulation demonstrates that the protein breaks the lateral and transverse symmetry of the bilayer. Lipids near the protein preferentially reorient such that their unsaturated chains interact with the protein, while the distribution of cholesterol in the membrane complements the variations in rhodopsin's transverse profile. The latter phenomenon suggests a molecular-level mechanism for the experimental finding that cholesterol stabilizes the native dark-adapted state of rhodopsin without binding directly to the protein.

Cholesterol↗

Molecular dynamics simulations of phospholipases A2.

An extensive molecular dynamics study of phospholipases A2 from pancreatic bovine and Crotalus atrox venom has shown that the well-conserved homologous core of the phospholipases A2, including the so called catalytic network, is very stable during the course of the calculations. The fluctuations which occur are located in segments which have significantly different three-dimensional conformations in the two phospholipases A2 studied, suggesting that a particularly stable core conformation gives rise to a large homologous family of similar three-dimensional structure. The calcium ion, which exhibits a crucial structural role in the monomeric phospholipases A2, appears not to be required to stabilize the C.atrox dimer. Moreover, the behaviour of the dimeric structure during the dynamics raises the question of a possible dissociation of the two subunits into functional monomers.

Animals↗