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 937 records · Page 52Linked to original sources

Molecular dynamics simulations on the Escherichia coli ammonia channel protein AmtB: mechanism of ammonia/ammonium transport.

Molecular dynamics (MD) simulations have been performed at the atomic level to study the ammonium/ammonia transport across the Escherichia coli AmtB membrane protein. Although ammonia primarily exists in the form of NH(4)(+) in aqueous solution, the recent X-ray structure determination of AmtB reveals that the ammonium/ammonia transporter proteins are ammonia-conducting channels rather than ammonium ion transporters [Khademi, S.; et al. Science 2004, 305, 1587; Zheng, L.; et al. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 17090]. Our simulations showed that the entrance of NH(4)(+) into the periplasmic recruitment vestibule requires only 3.1 kcal/mol of energy. This is consistent with the X-ray crystal structure, where one NH(4)(+) is captured in the binding vestibule. In this vestibule, NH(4)(+) loses one water of hydration, but the loss is compensated by a hydrogen bond, first with the backbone carbonyl oxygen of Phe161 then with the hydroxyl group of Ser219, as well as the stabilizing pi-cation interactions with the aromatic rings of Trp148 and Phe107 in the AmtB protein. In the end of this recruitment vestibule, the phenyl ring of Phe107 dynamically switches to an open state. This is correlated with a slight rotation and shifting of the indole ring of Trp148, which eventually creates a slot for the initially buried carboxylate group of Asp160 to become exposed to the bulk solvent. A hydrogen bond wire between NH(4)(+) and the carboxylate group of Asp160 via two water molecules was observed. Thus, Asp160 is most likely the proton acceptor from NH(4)(+). This explains the high conservation of Asp160 in Amt proteins and why the D160A mutant would completely quench the activity of AmtB [Javelle, A.; et al. J. Biol. Chem. 2004, 279, 8530; Marini, A. M.; et al. Curr. Genet. 2006, 49, 364]. Once NH(4)(+) deprotonates, the phenyl ring of Phe215 rotates to open, and the subsequent passage of NH(3) through the channel is straightforward.

Ammonia↗

Molecular dynamics simulations of the effects of ring-saturated thymine lesions on DNA structure.

The effect of thymine lesions produced by radiation or oxidative damage on DNA structure was studied by molecular dynamics simulations of native and damaged DNA. Thymine in position 7 of native dodecamer d(CGCGAATTCGCG)2 was replaced by one of the four thymine lesions 5-hydroxy-5,6-dihydrothymine, 6-hydroxy-5,6-dihydrothymine (thymine photohydrate), 5,6-dihydroxy-5,6-dihydro-thymine (thymine glycol), and 5,6-dihydrothymine. Simulations were performed with Assisted Model Building with Energy Refinement force field. Solvent was represented by a rectangular box of water with periodic boundary conditions applied. A constant temperature and constant volume protocol was used. The observed level of distortions of DNA structure depends on the specific nature of the lesion. The 5,6-dihydrothymine does not cause distinguishable perturbations to DNA. Other lesions produce a dramatic increase in the rise parameter between the lesion and the 5' adjacent adenine. These changes are accompanied by weakening of Watson-Crick hydrogen bonds in the A6-T19 base pair on the 5' side of the lesion. The lesioned bases also show negative values of inclination relative to the helical axis. No changes in the pattern of backbone torsional angles are observed with any of the lesions incorporated into DNA. The structural distortions in DNA correlate well with known biological effects of 5,6-dihydrothymine and thymine glycol on such processes as polymerase action or recognition by repair enzymes.

Base Sequence↗

New scoring functions for virtual screening from molecular dynamics simulations with a quantum-refined force-field (QRFF-MD). Application to cyclin-dependent kinase 2.

A recently introduced new methodology based on ultrashort (50-100 ps) molecular dynamics simulations with a quantum-refined force-field (QRFF-MD) is here evaluated in its ability both to predict protein-ligand binding affinities and to discriminate active compounds from inactive ones. Physically based scoring functions are derived from this approach, and their performance is compared to that of several standard knowledge-based scoring functions. About 40 inhibitors of cyclin-dependent kinase 2 (CDK2) representing a broad chemical diversity were considered. The QRFF-MD method achieves a correlation coefficient, R(2), of 0.55, which is significantly better than that obtained by a number of traditional approaches in virtual screening but only slightly better than that obtained by consensus scoring (R(2) = 0.50). Compounds from the Available Chemical Directory, along with the known active compounds, were docked into the ATP binding site of CDK2 using the program Glide, and the 650 ligands from the top scored poses were considered for a QRFF-MD analysis. Combined with structural information extracted from the simulations, the QRFF-MD methodology results in similar enrichment of known actives compared to consensus scoring. Moreover, a new scoring function is introduced that combines a QRFF-MD based scoring function with consensus scoring, which results in substantial improvement on the enrichment profile.

Computer Simulation↗

Mechanism of ABC transporters: a molecular dynamics simulation of a well characterized nucleotide-binding subunit.

ATP-binding cassette (ABC) transporters are membrane-bound molecular pumps that form one of the largest of all protein families. Several of them are central to phenomena of biomedical interest, including cystic fibrosis and resistance to chemotherapeutic drugs. ABC transporters share a common architecture comprising two hydrophilic nucleotide-binding domains (NBDs) and two hydrophobic transmembrane domains (TMDs) that form the substrate pathway across the membrane. The conformational changes in the NBDs induced by ATP hydrolysis and the means by which they are transmitted to the TMDs to effect substrate translocation remain largely unknown. We have performed a molecular dynamics simulation of HisP, the well studied NBD of the bacterial histidine permease, to identify hinges and switches of the NBD conformational transitions and subunit-subunit interfaces. This analysis reveals that the TMDs regulate ATP hydrolysis by controlling conformational transitions of the NBD helical domains, and identifies the conformational changes and the crucial TMD:NBD interface, by which the energy of ATP hydrolysis is transmitted to the TMDs. We also define the conformational transitions of the Q-loop, a key element of the NBD mechanism, and identify pathways by which Q-loop switching is coordinated with TMD and NBD conformational changes. We propose a model for the catalytic cycle of ABC transporters that shows how substrate-binding and transport by the TMDs may be coordinated and coupled with ATP binding and hydrolysis in the NBDs.

ATP-Binding Cassette Transporters↗

Molecular dynamics simulations of peptides from BPTI: a closer look at amide-aromatic interactions.

Molecular dynamics (MD) simulations of short peptides in water were performed to establish whether it is possible to reproduce experimental data from chemical shift measurements by nuclear magnetic resonance spectroscopy. Three different tetrapeptides were studied. The first, YTGP (Tyr-Thr-Gly-Pro), shows an electrostatic interaction between the aromatic ring of Tyr and the backbone amide hydrogen atom of Gly. The second, YTAP (Tyr-Thr-Ala-Pro), cannot make such an interaction because of steric hindrance of the Ala side chain and hence does not show a well-defined conformation. The third, FTGP (Phe-Thr-Gly-Pro), is shown to alternate between two different conformations. It is demonstrated that small differences in chemical shift, corresponding to these slightly different conformations, can be quantitatively modeled in MD simulations when using the proper force-field parameters and water model Explicit inclusion of hydrogen atoms o the aromatic rings is essential for a proper description of electrostatic interactions, but the choice of the water model is equally important. We found that a combination of the SPC/E water model and a revised GROMOS87 force field gives close agreement with experiment, while the same and other force fields in combination with SPC or TIP3P water did not reproduce the NMR data at all. Simulations of a longer peptide from bovine pancreatic trypsin inhibitor, containing the YTGP sequence, did show the interaction between the aromatic ring and the amide hydrogen, but not as pronounced as the simulations of shorter periods.

Amides↗

Molecular dynamics simulation of the structure of dimyristoylphosphatidylcholine bilayers with cholesterol, ergosterol, and lanosterol.

Five molecular dynamics computer simulations were performed on different phospholipid:sterol membrane systems in order to study the influence of sterol structure on membrane properties. Three of these simulated bilayer systems were composed of a 1:8 sterol:phospholipid ratio, each of which employed one of the sterol molecules: cholesterol, ergosterol, and lanosterol. The two other simulations were of a bilayer with a 1:1 sterol:phospholipid ratio. These simulations employed cholesterol and lanosterol, respectively, as their sterol components. The observed differences in simulations with cholesterol and lanosterol may have their implication on the form of the phospholipid/sterol phase diagram.

Carbon↗

An efficient molecular dynamics simulation method for calculating the diffusion-influenced reaction rates.

We present a molecular dynamics (MD) simulation method for calculating the diffusion-influenced reaction rates in the limit of low reactant concentrations. To calculate the reaction rate coefficient, we use MD trajectories of a nonreactive equilibrium system that are initiated with a pair of reactant molecules in reactive configuration. Hence reaction systems involving complicated reactant molecules with geometrically restricted reactivities can be treated with comparable efficiency as the simple hard-sphere reaction system. Compared to the similar MD method proposed by Van Beijeren, Dong, and Bocquet [J. Chem. Phys. 114, 6265 (2001)], the present method has a couple of advantages. First, reactions involving more general sink functions can be treated. Second, more accurate results can be obtained when the reaction probability upon collision is less than unity. As an application, we investigate the effects of nondiffusive dynamics and hydrodynamic interaction of reactants on the reaction rate.

Journal Article↗

Phase transitions in electrorheological fluids using molecular dynamics simulations.

A parametric study of the properties of electrorheological fluids is conducted using molecular dynamics (MD) simulations. The MD model is based on the solution of the Langevin equation for a number of suspended particles. The equations of motion include inertial effects, polarization forces, Stokes' drag, short range repulsion, and Brownian forces. Different polarization forces are considered to include the effect of enhancements at short range due to multipole moments induced by the suspended particles and other effects. The model is used to investigate the structural changes induced by external electric fields and by shear strains imposed on the system. The response times are studied as a function of two characteristic parameters describing the physical status of the system (temperature and external electric field). Finally, the stress-strain characteristics are studied and the yield stress is calculated as a function of the external electric field. The simulated response is compared with experimental findings.

Journal Article↗

Insights into saquinavir resistance in the G48V HIV-1 protease: quantum calculations and molecular dynamic simulations.

The spread of acquired immune deficiency syndrome has increasingly become a great concern owing largely to the failure of chemotherapies. The G48V is considered the key signature residue mutation of HIV-1 protease developing with saquinavir therapy. Molecular dynamics simulations of the wild-type and the G48V HIV-1 protease complexed with saquinavir were carried out to explore structure and interactions of the drug resistance. The molecular dynamics results combined with the quantum-based and molecular mechanics Poisson-Boltzmann surface area calculations indicated a monoprotonation took place on D25, one of the triad active site residues. The inhibitor binding of the triad residues and its interaction energy in the mutant were similar to those in the wild-type. The overall structure of both complexes is almost identical. However, the steric conflict of the substituted valine results in the conformational change of the P2 subsite and the disruption of hydrogen bonding between the -NH of the P2 subsite and the backbone -CO of the mutated residue. The magnitude of interaction energy changes was comparable to the experimental K(i) data. The designing for a new drug should consider a reduction of steric repulsion on P2 to enhance the activity toward this mutant strain.

Amino Acid Substitution↗

Molecular-dynamics simulation of amphiphilic bilayer membranes and wormlike micelles: a multi-scale modelling approach to the design of viscoelastic surfactant solutions.

Bilayer membranes and wormlike micelles have been studied using molecular-dynamics simulations. The structure of the worm is analysed in terms of radial density distribution functions, and mechanical properties such as the elastic modulus are calculated. From an analysis of the fluctuation spectra of the tensionless states, we have calculated bending rigidities. Micelles consisting of coarse-grained (CG) model surfactants are studied in order to map the properties of the atomistic micelle. We optimize the CG model with respect to the structure factor of the atomistic micelle. The mechanical properties thus obtained will be used as input for a mesoscopic model of wormlike micelles where the persistence length is the smallest length-scale.

Colloids↗

Conformational dynamics of the nicotinic acetylcholine receptor channel: a 35-ns molecular dynamics simulation study.

The nicotinic acetylcholine receptor (AChR) is the paradigm of ligand-gated ion channels, integral membrane proteins that mediate fast intercellular communication in response to neurotransmitters. A 35-ns molecular dynamics simulation has been performed to explore the conformational dynamics of the entire membrane-spanning region, including the ion channel pore of the AChR. In the simulation, the 20 transmembrane (TM) segments that comprise the whole TM domain of the receptor were inserted into a large dipalmitoylphosphatidylcholine (DPPC) bilayer. The dynamic behavior of individual TM segments and their corresponding AChR subunit helix bundles was examined in order to assess the contribution of each to the conformational transitions of the whole channel. Asymmetrical and asynchronous motions of the M1-M3 TM segments of each subunit were revealed. In addition, the outermost ring of five M4 TM helices was found to convey the effects exerted by the lipid molecules to the central channel domain. Remarkably, a closed-to-open conformational shift was found to occur in one of the channel ring positions in the time scale of the present simulations, the possible physiological significance of which is discussed.

1,2-Dipalmitoylphosphatidylcholine↗

Probing conformational disorder in neurotensin by two-dimensional solid-state NMR and comparison to molecular dynamics simulations.

An approach is introduced to characterize conformational ensembles of intrinsically unstructured peptides on the atomic level using two-dimensional solid-state NMR data and their combination with molecular dynamics simulations. For neurotensin, a peptide that binds with high affinity to a G-protein coupled receptor, this method permits the investigation of the changes in conformational preferences of a neurotransmitter transferred from a frozen aqueous solution via a lipid model phase to the receptor-bound form. The results speak against a conformational pre-organization of the ligand in detergents in which the receptor has been shown to be functional. Further extensions to the study of protein folding are possible.

Detergents↗

100ps molecular dynamic simulation of d(TATCACC)2.

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage lambda is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H-bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.

Bacteriophage lambda↗

Molecular dynamics simulation of unsaturated lipid bilayers at low hydration: parameterization and comparison with diffraction studies.

A potential energy function for unsaturated hydrocarbons is proposed and is shown to agree well with experiment, using molecular dynamics simulations of a water/octene interface and a dioleoyl phosphatidylcholine (DOPC) bilayer. The simulation results verify most of the assumptions used in interpreting the DOPC experiments, but suggest a few that should be reconsidered. Comparisons with recent results of a simulation of a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer show that disorder is comparable, even though the temperature, hydration level, and surface area/lipid for DOPC are lower. These observations highlight the dramatic effects of unsaturation on bilayer structure.

1,2-Dipalmitoylphosphatidylcholine↗

Hydration of the calcium ion. An EXAFS, large-angle x-ray scattering, and molecular dynamics simulation study.

The structure of the hydrated calcium(II) ion in aqueous solution has been studied by means of extended X-ray absorption fine structure spectroscopy (EXAFS), large-angle X-ray scattering (LAXS), and molecular dynamics (MD) methods. The EXAFS data displayed a broad and asymmetric distribution of the Ca-O bond distances with the centroid at 2.46(2) A. LAXS studies on four aqueous calcium halide solutions (1.5-2 mol dm(-)(3)) gave a mean Ca-O bond distance of 2.46(1) A. This is consistent with a hydration number of 8 determined from correlations between mean distances and coordination numbers from crystal structures. The LAXS studies showed a second coordination sphere with a mean Ca.O(II) distance of 4.58(5) A, and for the hydrated halide ions the distances Cl.O 3.25(1) A, Br.O 3.36(1) A, and I.O 3.61(1) A were obtained. Molecular dynamics simulations of CaCl(2)(aq) were performed using three different Ca(2+)-OH(2) pair potentials. The potential from the GROMOS program gave results in agreement with experiments, i.e., a coordination number of 8 and an average Ca-O distance of 2.46 A, and was used for further comparisons. Theoretical EXAFS oscillations were computed for individual MD snapshots and showed very large variations, though the simulated average spectrum from 2000 snapshots gave satisfactory agreement with the experimental EXAFS spectra. The effect of thermal motions of the coordinated atoms is inherent in the MD simulation method. Thermal disorder parameters evaluated from simulated spatial atom distribution functions of the oxygen atoms coordinated to the calcium ion were in close agreement with those from the current LAXS and EXAFS analyses. The combined results are consistent with a root-mean-square displacement from the mean Ca-O distance of 0.09(2) A in aqueous solution at 300 K.

Body Water↗

Spinodal decomposition in thin films: molecular-dynamics simulations of a binary Lennard-Jones fluid mixture.

We use molecular dynamics (MD) to simulate an unstable homogeneous mixture of binary fluids (AB), confined in a slit pore of width D. The pore walls are assumed to be flat and structureless and attract one component of the mixture (A) with the same strength. The pairwise interactions between the particles are modeled by the Lennard-Jones potential, with symmetric parameters that lead to a miscibility gap in the bulk. In the thin-film geometry, an interesting interplay occurs between surface enrichment and phase separation. We study the evolution of a mixture with equal amounts of A and B, which is rendered unstable by a temperature quench. We find that A-rich surface enrichment layers form quickly during the early stages of the evolution, causing a depletion of A in the inner regions of the film. These surface-directed concentration profiles propagate from the walls towards the center of the film, resulting in a transient layered structure. This layered state breaks up into a columnar state, which is characterized by the lateral coarsening of cylindrical domains. The qualitative features of this process resemble results from previous studies of diffusive Ginzburg-Landau-type models [S. K. Das, S. Puri, J. Horbach, and K. Binder, Phys. Rev. E 72, 061603 (2005)], but quantitative aspects differ markedly. The relation to spinodal decomposition in a strictly two-dimensional geometry is also discussed.

Journal Article↗

Molecular dynamics simulation and essential dynamics study of mutated plastocyanin: structural, dynamical and functional effects of a disulfide bridge insertion at the protein surface.

A molecular dynamics simulation (1.1 ns) at 300 K, of fully hydrated Ile21Cys, Glu25Cys plastocyanin mutant has been performed to investigate the structural, dynamical and functional effects of a disulfide bridge insertion at the surface of the protein. A detailed analysis of the root mean square fluctuations, H-bonding pattern and dynamical cross-correlation map has been performed. An essential dynamics method has also been applied as complementary analysis to identify concerted motions (essential modes), that could be relevant to the electron transfer function. The results have been compared with those previously obtained for wild-type plastocyanin and have revealed that the mutant shows a different pattern of H-bonds, with several interactions lost and a higher flexibility, especially around the electron transfer copper site. The analysis of dynamical cross-correlation map and of essential modes, has shown that the mutant performs different functional concerted motions, which might be related to the binding recognition with its electron transfer partners in comparison with the wild-type protein.

Algorithms↗

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↗