Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “molecular dynamics”

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 55 records · Page 3Linked to original sources
PubMed2005

Examination of dynamic facilitation in molecular dynamics simulations of glass-forming liquids.

Using data from molecular dynamics computer simulations of the one-component Dzugutov liquid and of BKS silica in metastable equilibrium supercooled states, we examine ideas introduced by Garrahan and Chandler (GC) in their dynamic facilitation (DF) model of the glass transition. Utilizing a recently introduced measure of DF, we find that DF is important for particle motion in both the supercooled Dzugutov liquid and in the BKS silica melt, that mobility propagates continuously, and that this effect becomes increasingly pronounced with decreasing T. We show that, in both systems, dynamic facilitation is strongest on the time scale of the late-beta relaxation, where clusters of highly mobile neighboring particles escaping from their cages are largest and, except for the silicon atoms in BKS silica, stringlike motion is most prominent. By comparing the two systems, we show that the temperature dependence of one measure of DF as the mode-coupling temperature is approached from high temperature is similar, once the temperature dependence of the structural relaxation time in each system is scaled out.

Journal Article↗
PubMed2005

Molecular dynamics integration and molecular vibrational theory. I. New symplectic integrators.

New symplectic integrators have been developed by combining molecular dynamics integration with the standard theory of molecular vibrations to solve the Hamiltonian equations of motion. The presented integrators analytically resolve the internal high-frequency molecular vibrations by introducing a translating and rotating internal coordinate system of a molecule and calculating normal modes of an isolated molecule only. The translation and rotation of a molecule are treated as vibrational motions with the vibrational frequency zero. All types of motion are thus described in terms of the normal coordinates. The method's time reversibility requirement was used to determine the equations of motion for internal coordinate system of a molecule. The calculation of long-range forces is performed numerically within the generalized second-order leap-frog scheme, in the same way as in standard second-order symplectic methods. The new methods for integrating classical equations of motion using normal mode analysis allow us to use a long integration step and are applicable to any system of molecules with one equilibrium configuration.

Journal Article↗
PubMed1994

Computer simulations of cyclic and acyclic cholinergic agonists: conformational search and molecular dynamics simulations.

Molecular dynamics simulations have been performed on aqueous solutions of two chemically similar nicotinic cholinergic agonists in order to compare their structural and dynamical differences. The cyclic 1,1-dimethyl-4-acetylpiperazinium iodide (HPIP) molecule was previously shown to be a strong agonist for nicotinic acetylcholine receptors (McGroddy et al., 1993), while the acyclic N,N,N,N'-tetramethyl-N'-acetylethylenediamine iodide (HTED) derivative is much less potent. These differences were expected to arise from differences in the solution structures and internal dynamics of the two molecules. HPIP was originally thought to be relatively rigid; however, molecular dynamics simulations suggest that the acetyl portion of the molecule undergoes significant ring dynamics on a psec timescale. The less constrained HTED molecule is relatively rigid, with only one transition observed about any of the major dihedrals in four 100 psec simulations, each started from a different conformation. The average structures obtained from the simulations are very similar to the starting minimized structure in each case, except for the HTED simulation where a single rotation about the N-C-C-N(+) backbone occurred. In each case, HTED had three to five more water molecules in its primary solvation shell than HPIP, indicating that differences in the energetics of desolvation before binding may partially explain the increased potency of HPIP as compared to HTED.

Biophysical Phenomena↗
PubMed1999

Effect of G40R mutation on the binding of human SRY protein to DNA: a molecular dynamics view.

Molecular dynamics simulation was conducted to investigate the reason why the mutant G40R of hSRY protein has a low affinity for DNA. Compared with the previous dynamics results of the wild-type hSRY-HMG-DNA complex, the results of molecular dynamics simulation on the mutant G40R hSRY-HMG-DNA system demonstrated that the whole structure of DNA (especially the second strand) had a major deviation away from the short arm of the HMG box. Consequently, the DNA and the mutant protein could not specifically recognize each other, that is, very different, and low-occupancy, direct, and water-mediated hydrogen bonds were detected at the protein-DNA interface, no conformational changes occurred at the loop region around Met9 during the simulation, and residue IIe13 did not intercalate between the bases of A5 and A6. These results indicated that the mutant G40R did not form a specific complex with the DNA target, hence led to complete gonadal dysgenesis. From the simulation, we realized that the residue Gly40 played a critical structural role in the hSRY-DNA recognition. It might be a structural supporting point of DNA binding because of the absence of a side chain. The reason for the difficulty of the mutant G40R to form a complex with DNA might be that the long and positively charged side chain of Arg40 by its bulk and positive charge hindered the DNA's access to the active sites of the protein.

DNA-Binding Proteins↗
PubMed2004

Change in protein flexibility upon complex formation: analysis of Ras-Raf using molecular dynamics and a molecular framework approach.

Changes in flexibility upon protein-protein complex formation of H-Ras and the Ras-binding domain of C-Raf1 have been investigated using the molecular framework approach FIRST (Floppy Inclusion and Rigid Substructure Topology) and molecular dynamics simulations (MD) of in total approximately 35 ns length. In a computational time of about one second, FIRST identifies flexible and rigid regions in a single, static three-dimensional molecular framework, whose vertices represent protein atoms and whose edges represent covalent and non-covalent (hydrogen bond and hydrophobic) constraints and fixed bond angles within the protein. The two methods show a very good agreement with respect to the identification of changes in flexibility in both binding partners on a local scale. This implies that flexibility can be successfully predicted by identifying which bonds limit motion within a molecule and how they are coupled. In particular, as identified by MD, the beta-sheet in Raf shows considerably more pronounced orientational correlations in the bound state compared to the unbound state. Similarly, FIRST assigns the beta-sheet to the largest rigid cluster of the complex. Interestingly, FIRST allows us to identify that interactions across the interface (but not conformational changes upon complex formation) result in the observed rigidification. Since regions of the beta-sheet of Raf that do not interact directly with Ras become rigidified, this also demonstrates the long-range aspect to rigidity percolation. Possible implications of the change of flexibility of the Ras-binding domain of Raf on the activation of Raf upon complex formation are discussed. Finally, the sensitivity of FIRST results with respect to the representation of non-covalent interactions used as constraints is probed.

Algorithms↗
PubMed2005

Molecular dynamics integration and molecular vibrational theory. II. Simulation of nonlinear molecules.

A series of molecular dynamics (MD) simulations of nonlinear molecules has been performed to test the efficiency of newly introduced semianalytical second-order symplectic time-reversible MD integrators that combine MD and the standard theory of molecular vibrations. The simulation results indicate that for the same level of accuracy, the new algorithms allow significantly longer integration time steps than the standard second-order symplectic leap-frog Verlet method. Since the computation cost per integration step using new MD integrators with longer time steps is approximately the same as for the standard method, a significant speed-up in MD simulation is achieved.

Journal Article↗
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↗
PubMed2005

Molecular dynamics integration and molecular vibrational theory. III. The infrared spectrum of water.

The new symplectic molecular dynamics (MD) integrators presented in the first paper of this series were applied to perform MD simulations of water. The physical properties of a system of flexible TIP3P water molecules computed by the new integrators, such as diffusion coefficients, orientation correlation times, and infrared (IR) spectra, are in good agreement with results obtained by the standard method. The comparison between the new integrators' and the standard method's integration time step sizes indicates that the resulting algorithm allows a 3.0 fs long integration time step as opposed to the standard leap-frog Verlet method, a sixfold simulation speed-up. The accuracy of the method was confirmed, in particular, by computing the IR spectrum of water in which no blueshifting of the stretching normal mode frequencies is observed as occurs with the standard method.

Journal Article↗
PubMed1997

Photoproducts of bacteriorhodopsin mutants: a molecular dynamics study.

Molecular dynamics simulations of wild-type bacteriorhodopsin (bR) and of its D85N, D85T, D212N, and Y57F mutants have been carried out to investigate possible differences in the photoproducts of these proteins. For each mutant, a series of 50 molecular dynamics simulations of the photoisomerization and subsequent relaxation process were completed. The photoproducts can be classified into four distinct classes: 1) 13-cis retinal, with the retinal N-H+ bond oriented toward Asp-96; 2) 13-cis retinal, with the N-H+ oriented toward Asp-85 and hydrogen-bonded to a water molecule; 3) 13,14-di-cis retinal; 4) all-trans retinal. Simulations of wild-type bR and of its Y57F mutant resulted mainly in class 1 and class 2 products; simulations of D85N, D85T, and D212N mutants resulted almost entirely in class 1 products. The results support the suggestion that only class 2 products initiate a functional pump cycle. The formation of class 1 products for the D85N, D85T, and D212N mutants can explain the reversal of proton pumping under illumination by blue and yellow light.

Amino Acid Sequence↗
PubMed2006

Spin-labeled alkylphospholipids in a dipalmitoylphosphatidylcholine bilayer: molecular dynamics simulations.

Molecular dynamics simulation has been performed to investigate the structural properties of perifosine and its synthetic spin-labeled alkylphospholipid analogues. The conformations adopted by these compounds in water and in a dipalmitoylphosphatidylcholine bilayer as a function of the presence and position of the N-oxyl-4',4'-dimethyloxazolidine ring (doxyl group) have been investigated by all-atom molecular dynamics. No predominant conformation was observed in water, but the molecules adopt specific orientations and conformations in the lipid bilayer. As is expected, alkyl chains tend to insert into the hydrophobic core, while charged groups stay at the lipid-water interface. A doxyl group in the middle of the alkyl chain moves up to the interface region, thus preventing adoption of the extended conformation. Compounds with a doxyl group close to the polar head group adopt conformations similar to that of unlabeled perifosine within the first nanoseconds of simulation. When the doxyl group is at the end of alkyl chain, the spin-labeled molecule needs more time to reach equilibrium. These results indicate a considerable effect of the doxyl position within the alkyl chain on its localization in the lipid bilayer and can be extended further to other similar spin probes used in the electron paramagnetic resonance spectroscopy of biological membranes.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2002

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations.

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.

Animals↗
PubMed2006

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation.

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.

Letter↗
PubMed2006

Molecular dynamics of the Shewanella oneidensis response to chromate stress.

Temporal genomic profiling and whole-cell proteomic analyses were performed to characterize the dynamic molecular response of the metal-reducing bacterium Shewanella oneidensis MR-1 to an acute chromate shock. The complex dynamics of cellular processes demand the integration of methodologies that describe biological systems at the levels of regulation, gene and protein expression, and metabolite production. Genomic microarray analysis of the transcriptome dynamics of midexponential phase cells subjected to 1 mm potassium chromate (K(2)CrO(4)) at exposure time intervals of 5, 30, 60, and 90 min revealed 910 genes that were differentially expressed at one or more time points. Strongly induced genes included those encoding components of a TonB1 iron transport system (tonB1-exbB1-exbD1), hemin ATP-binding cassette transporters (hmuTUV), TonB-dependent receptors as well as sulfate transporters (cysP, cysW-2, and cysA-2), and enzymes involved in assimilative sulfur metabolism (cysC, cysN, cysD, cysH, cysI, and cysJ). Transcript levels for genes with annotated functions in DNA repair (lexA, recX, recA, recN, dinP, and umuD), cellular detoxification (so1756, so3585, and so3586), and two-component signal transduction systems (so2426) were also significantly up-regulated (p < 0.05) in Cr(VI)-exposed cells relative to untreated cells. By contrast, genes with functions linked to energy metabolism, particularly electron transport (e.g. so0902-03-04, mtrA, omcA, and omcB), showed dramatic temporal alterations in expression with the majority exhibiting repression. Differential proteomics based on multidimensional HPLC-MS/MS was used to complement the transcriptome data, resulting in comparable induction and repression patterns for a subset of corresponding proteins. In total, expression of 2,370 proteins were confidently verified with 624 (26%) of these annotated as hypothetical or conserved hypothetical proteins. The initial response of S. oneidensis to chromate shock appears to require a combination of different regulatory networks that involve genes with annotated functions in oxidative stress protection, detoxification, protein stress protection, iron and sulfur acquisition, and SOS-controlled DNA repair mechanisms.

Anion Transport Proteins↗
PubMed2005

A new molecular-dynamics based approach for molecular crystal structure search.

A new molecular-dynamics based approach is proposed to search for candidate crystal structures of molecular solids. The procedure is based on the observation of spontaneous transitions between ordered and disordered states in molecular-dynamics simulations of an artificial periodic system with a small unit cell. In such a way only the most stable structures are automatically selected. The method can be applied to the solution of crystal structures from low-quality or very complex diffraction data. Tests are presented for H2O-ice polymorphs.

Journal Article↗
PubMed2003

B lymphocyte activation during cognate interactions with CD4+ T lymphocytes: molecular dynamics and immunologic consequences.

Productive interaction between T and B lymphocytes is required for humoral immune responses to many foreign protein antigens and production of pathogenic antibodies characteristic of several autoimmune conditions. Thus, much attention has been given in recent years to understand the dynamic molecular interactions and signal transduction required for productive T-B interaction. In this review we highlight current knowledge of signaling and biologic responses that occur in B cells during cognate interactions with helper T cells, focusing on the dynamic function of B cell-surface molecules in T-B synapses.

B-Lymphocytes↗
PubMed2004

Intramolecular diffusive motion in alkane monolayers studied by high-resolution quasielastic neutron scattering and molecular dynamics simulations.

Molecular dynamics simulations of a tetracosane (n-C24H50) monolayer adsorbed on a graphite basal-plane surface show that there are diffusive motions associated with the creation and annihilation of gauche defects occurring on a time scale of approximately 0.1-4 ns. We present evidence that these relatively slow motions are observable by high-energy-resolution quasielastic neutron scattering (QNS) thus demonstrating QNS as a technique, complementary to nuclear magnetic resonance, for studying conformational dynamics on a nanosecond time scale in molecular monolayers.

Journal Article↗
PubMed2006

Vibrational dynamics of DNA. III. Molecular dynamics simulations of DNA in water and theoretical calculations of the two-dimensional vibrational spectra.

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dG(n):dC(n) and dA(n):dT(n) double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG(5)C(5) and dA(5)T(5) double helices in D(2)O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP6-31G(*) calculations for the nine basis modes. The simulated 2D IR spectra of dG(n):dC(n) double helix in D(2)O solution are directly compared with closely related experimental results. The 2D IR spectra of dG(n):dC(n) and dA(n):dT(n) are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions.

Base Pairing↗