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 1,765 records · Page 98Linked to original sources
PubMed2006

Periodic boundary condition induced breakdown of the equipartition principle and other kinetic effects of finite sample size in classical hard-sphere molecular dynamics simulation.

We examine consequences of the non-Boltzmann nature of probability distributions for one-particle kinetic energy, momentum, and velocity for finite systems of classical hard spheres with constant total energy and nonidentical masses. By comparing two cases, reflecting walls (NVE or microcanonical ensemble) and periodic boundaries (NVEPG or molecular dynamics ensemble), we describe three consequences of the center-of-mass constraint in periodic boundary conditions: the equipartition theorem no longer holds for unequal masses, the ratio of the average relative velocity to the average velocity is increased by a factor of [N/(N-1)]1/2, and the ratio of average collision energy to average kinetic energy is increased by a factor of N/(N-1). Simulations in one, two, and three dimensions confirm the analytic results for arbitrary dimension.

Journal Article↗
PubMed2007

Effects of chain length on the aggregation of model polyglutamine peptides: molecular dynamics simulations.

Aggregation in the brain of polyglutamine-containing proteins is either a cause or an associated symptom of nine hereditary neurodegenerative disorders including Huntington's disease. The molecular level mechanisms by which these proteins aggregate are still unclear. In an effort to shed light on this important phenomenon, we are investigating the aggregation of model polyglutamine peptides using molecular-level computer simulation with a simplified model of polyglutamine that we have developed. This model accounts for the most important types of intra- and inter-molecular interactions-hydrogen bonding and hydrophobic interactions-while allowing the folding process to be simulated in a reasonable time frame. The model is used to examine the folding of isolated polyglutamine peptides 16, 32, and 48 residues long and the folding and aggregation of systems of 24 model polyglutamine peptides 16, 24, 32, 36, 40, and 48 residues long. Although the isolated polyglutamine peptides did form some alpha and beta backbone-backbone hydrogen bonds they did not have as many of these bonds as they would have if they had folded into a complete alpha helix or beta sheet. In one of the simulations on the isolated polyglutamine peptide 48 residues long, we observed a structure that resembles a beta helix. In the multi-chain simulations we observed amorphous aggregates at low temperatures, ordered aggregates with significant beta sheet character at intermediate temperatures, and random coils at high temperatures. We have found that the temperature at which the model peptides undergo the transition from amorphous aggregates to ordered aggregates and the temperature at which the model peptides undergo the transition from ordered aggregates to random coils increase with increasing chain length. Our finding that the stability of the ordered aggregates increases as the peptide chain length increases may help to explain the experimentally observed relation between polyglutamine tract length and aggregation in vitro and disease progression in vivo. We have also observed in our simulations that the optimal temperature for the formation of beta sheets increases with chain length up to 36 glutamine residues but not beyond. Equivalently, at fixed temperature we find a transition from a region dominated by random coils at chain lengths less than 36 to a region dominated by relatively ordered beta sheet structures at chain lengths greater than 36. Our finding of this critical chain length of 36 glutamine residues is interesting because a critical chain length of 37 glutamine residues has been observed experimentally.

Computer Simulation↗
PubMed1993

Solution structure of a conserved DNA sequence from the HIV-1 genome: restrained molecular dynamics simulation with distance and torsion angle restraints derived from two-dimensional NMR spectra.

The three-dimensional solution structure of a trisdecamer DNA duplex sequence, d(AGCTTGCCTTGAG).d(CTCAAGGCAAGCT), from a conserved region of HIV-1 genome's long terminal repeat, has been investigated using NMR spectroscopy and restrained molecular dynamics calculations. Interproton distances derived from two-dimensional nuclear Overhauser enhancement (2D NOE) experiments, using the iterative complete relaxation matrix algorithm MARDIGRAS, and torsion angles for sugar rings, estimated from stimulated fitting of double-quantum-filtered correlation (2QF-COSY) spectra, were obtained [Mujeeb, A., Kerwin, S. M., Egan, W. M., Kenyon, G. L., & James, T. L. (1992) Biochemistry 31, 9325-38]. These structural restraints have now been employed as the basis for structure refinement using restrained molecular dynamics (rMD) to search conformational space for structures consistent with the experimental restraints. Specifically, upper and lower bounds on the restraints were incorporated into the AMBER (version 4.0) total potential energy function of the system, the bounds being used to define the width of a flat-well penalty term in the AMBER force field. Confidence in the time-averaged structure obtained is engendered by convergence to essentially the same structure (root-mean-square deviation approximately 0.9 A) when two quite different DNA models, A-DNA and B-DNA (RMSD approximately 6.5 A), were employed as starting structures and when various initial trajectories were used for the rMD runs. The derived structure is further supported by the total energy calculated, the restraint violation energy, the restraint deviations, and the fit with experimental data. For the latter, the sixth-root residual index indicated a good fit of the determined structure with experimental 2D NOE spectral intensities (R1x < 0.07), and the RMS difference between vicinal proton coupling constants calculated for the derived structure and experimental coupling constants were also in reasonable agreement (JRMS = 0.9 Hz). While the structure of the trisdecamer is basically in the B-DNA family, some structural parameters manifest interesting local variations. The helix parameters indicate that, compared with classical B-DNA, the structure is longitudinally more compressed. Local structural variations at the two TG steps in particular together create bending into the major groove of the duplex. Comparison of the two-CTTG-tetrads in the duplex reveals that they have similar structures, with the TT moieties being almost identical; however, the -CTTG-pur sequence has a larger roll and slide for the -TG- step than for the -CTTG-pyr sequence, in accord with published X-ray crystallographic conclusions.

Base Sequence↗
PubMed2004

Molecular dynamics simulation of the interaction between the complex iron-sulfur flavoprotein glutamate synthase and its substrates.

Glutamate synthase (GltS) is a complex iron-sulfur flavoprotein that catalyzes the reductive transfer of L-glutamine amide group to the C2 carbon of 2-oxoglutarate yielding two molecules of L-glutamate. Molecular dynamics calculations in explicit solvent were carried out to gain insight into the conformational flexibility of GltS and into the role played by the enzyme substrates in regulating the catalytic cycle. We have modelled the free (unliganded) form of Azospirillum brasilense GltS alpha subunit and the structure of the reduced enzyme in complex with the L-glutamine and 2-oxoglutarate substrates starting from the crystallographically determined coordinates of the GltS alpha subunit in complex with L-methionine sulphone and 2-oxoglutarate. The present 4-ns molecular dynamics calculations reveal that the GltS glutaminase site may exist in a catalytically inactive conformation unable to bind glutamine, and in a catalytically competent conformation, which is stabilized by the glutamine substrate. Substrates binding also induce (1) closure of the loop formed by residues 263-271 with partial shielding of the glutaminase site from solvent, and (2) widening of the ammonia tunnel entrance at the glutaminase end to allow for ammonia diffusion toward the synthase site. The Q-loop of glutamate synthase, which acts as an active site lid in other amidotransferases, seems to maintain an open conformation. Finally, binding of L-methionine sulfone, a glutamine analog that mimics the tetrahedral transient species occurring during its hydrolysis, causes a coordinated rigid-body motion of segments of the glutaminase domain that results in the inactive conformation observed in the crystal structure of GltS alpha subunit.

Azospirillum brasilense↗
PubMed2005

Polyoxometalates in solution: Molecular dynamics simulations on the alpha-PW12O40(3-) Keggin anion in aqueous media.

The Keggin anion, PW12O40(3-), is one of the most representative polyoxometalates (POMs). In recent years increasing theoretical work focused on this family of compounds has explained or even predicted some of their properties using quantum mechanics methods. In this report we applied for the first time molecular dynamics (MD) to the title compound to analyze its interactions with water. We used three atomic charge definitions (Mulliken, ChelpG, and formal charges) to carry out MD simulations. The results show that the terminal oxygens of the cluster are invariably most effectively solvated by water because of their prominent position within the framework. On the other hand, bridging oxygens, which are confined in more internal positions, concentrate a smaller portion of the whole solvation. We investigated the hydrogen bonds existing between water and the cluster, confirming that the terminal positions form more contacts with H2O than any other site of the cluster. In the end, the lifetime of such contacts is longer with bridging oxygens, presumably due to their higher atomic charge.

Journal Article↗
PubMedDate not supplied

Solution structure studies of the cobalt complex of a bleomycin functional model bound to d(CGCAATTGCG)2 by two-dimensional nuclear magnetic resonance methods and restrained molecular dynamics simulation.

The interaction between the cobalt(III) complex of a bleomycin functional model (AMPHIS-NET) and the oligonucleotide d(CGCAATTGCG)2 and the structural features of the 1:1 ligand-DNA complex have been determined by high-resolution two-dimensional nuclear magnetic resonance methods and restrained molecular dynamics calculations. The intermolecular nuclear Overhauser effect (NOE) cross-peaks between ligand protons and the DNA minor groove protons suggest that the cobalt(III) complex of AMPHIS-NET binds in the minor groove of DNA at the central AATT site. The NOE connectivities also clearly indicate that the H8 pyridine proton and the H2 imidazole proton in the metal-binding domain interact with the H4' sugar proton of C19 and the H4' sugar proton of A5, respectively, which defines a structure where the metal binding moiety of Co(III).AMPHIS-NET participates in binding to the DNA and extends into the region two base pairs beyond the central AATT site in the minor groove. This binding model is in accord with the consistently observed nondiffusion DNA cleavage in locations two to three residues beyond the end of AT-rich binding sites induced by the corresponding iron(II) complexes of AMPHIS-NET and other AMPHIS-lexitropsin hybrids of the bleomycin functional model compounds.

Antiviral Agents↗
PubMed2003

Nonequilibrium molecular dynamics simulation of gas-mixtures transport in carbon-nanopore membranes.

The numerical simulation of nonequilibrium cotransport of H2 and alkane gas mixtures with very different molecular sizes through a porous carbon membrane structure was implemented. Simulated permeabilities and selectivities in binary diffusive systems (and in one ternary system), at pressure about tens of atmospheres and at operating condition of room temperature or higher, can be predicted from single-gas permeabilities. This suggests that the effect is geometrical and an approximate model of the transport is proposed. It can be used for an estimation of the separation factor of a membrane. Simulations are compared with experimental results of two- and three-component codiffusion and counterdiffusion in a carbon membrane. It is shown that diffusion in a porous molecular network and in a carbon nanotube are completely different. The uniqueness of this work lies in the comparison of simulated, approximate, and experimental results, which enables us to identify the important parameters.

Journal Article↗
PubMed2004

Mode-coupling theory and molecular dynamics simulation for heat conduction in a chain with transverse motions.

We study heat conduction in a 1D chain of particles with longitudinal as well as transverse motions. The particles are connected by 2D harmonic springs together with bending angle interactions. The problem is analyzed by mode-coupling theory and compared with molecular dynamics. We find very good, quantitative agreement for the damping of modes between a full mode-coupling theory and molecular dynamics result, and a simplified mode-coupling theory gives qualitative description of the damping. The theories predict generically that thermal conductance diverges as N(1/3) as the size N increases for systems terminated with heat baths at the ends. The N(2/5) dependence is also observed in molecular dynamics, which we attribute to crossover effect.

Journal Article↗
PubMed1994

Studies on the three-dimensional behaviour of the selectin ligands Lewis(a) and sulphated Lewis(a) using NMR spectroscopy and molecular dynamics simulations.

Sulphated blood group Lewis(a)/Lewis(x) (Le(a)/Le(x)) type sequences, with sulphate at the 3-position of galactose, have emerged as potent ligands for the endothelial adhesion molecule E-selectin and the leukocyte adhesion molecule L-selectin. As a first step in elucidating the molecular basis of the strong interactions with the selectins, we have performed conformational studies of the sulphated Le(a) in comparison with the non-sulphated analogue which is less strongly bound by E-selectin and not at all by L-selectin. Experimental NMR parameters [nuclear Overhauser effects (NOE) and interglycosidic 3JC,H] and theoretical values back-calculated from the minimum energy structures are in excellent agreement for both molecules. Molecular dynamics calculations for SuLe(a) depict only minor torsional fluctuations around the glycosidic linkages over the time course of the 500 ps simulations, leading to the conclusion that the conformation of SuLe(a) approximates to a single-rigid structure, as does the previously investigated Le(a) molecule. Comparison of experimentally and theoretically obtained parameters for SuLe(a) with those for the non-sulphated Le(a) molecule indicate that no detectable changes occur in the three-dimensional structure of the trisaccharide upon sulphation. Thus, the enhanced selectin binding to the sulphated Le(a) is most likely due to favourable electrostatic interactions between the charged sulphate group and corresponding charged groups on the selectin protein.

Carbohydrate Conformation↗
PubMed2000

Correlation between changes in nuclear magnetic resonance order parameters and conformational entropy: molecular dynamics simulations of native and denatured staphylococcal nuclease.

Recent work has suggested that changes in NMR order parameters may quantitatively reflect changes in the conformational entropy of a protein ensemble. The extent of the mathematical relationship between local entropy changes as seen by NMR order parameters and the full protein entropy change is a complex issue. As a step towards a fuller understanding of this problem, molecular dynamics calculations of both native and denatured staphylococcal nuclease were performed. The N-H bond vector motion, in both explicit and implicit solvent, was analyzed to estimate local and global entropy changes. The calculated N-H bond vector order parameters from simulation agreed on average with experimental values for both native and denatured structures. However, the inverted-U profile of order parameters versus residue number observed experimentally for denatured nuclease was only partially reproduced by simulation of compact denatured structures. Comparisons made across the full set of simulations revealed a correlation between the N-H order parameter-based conformational entropy change and the total quasiharmonic-based conformational entropy change between the native and denatured structures. The calculations showed that about 25% of the total entropy change was reflected by changes in simulated S2 values. This result suggests that NMR-derived order parameters may be used to provide a reasonable estimate of the total conformational entropy change on protein folding.

Magnetic Resonance Spectroscopy↗
PubMed2005

Stability of C60 chains: molecular dynamics simulations.

A linearly aligned structure of three C60 fullerene, interconnected by two benzorods of same size, have been investigated under heat treatment. The overall structure resembles a section of a beaded string. Nine different lengths of benzorods have been considered, and the effect on the thermal stability have been investigated by means of molecular dynamics method. It has been found that the structure is thermally stable up to elevated temperatures, and the linear alignment of the structure is persistent, up to the temperature of decomposition.

Fullerenes↗
PubMed2005

Self-diffusion and transport diffusion of light gases in metal-organic framework materials assessed using molecular dynamics simulations.

Metal-organic framework (MOF) materials pose an interesting alternative to more traditional nanoporous materials for a variety of separation processes. Separation processes involving nanoporous materials can be controlled by either adsorption equilibrium, diffusive transport rates, or a combination of these factors. Adsorption equilibrium has been studied for a variety of gases in MOFs, but almost nothing is currently known about molecular diffusion rates in MOFs. We have used equilibrium molecular dynamics (MD) to probe the self-diffusion and transport diffusion of a number of small gas species in several MOFs as a function of pore loading at room temperature. Specifically, we have studied Ar, CH4, CO2, N2, and H2 diffusion in MOF-5. The diffusion of Ar in MOF-2, MOF-3, and Cu-BTC has been assessed in a similar manner. Our results greatly expand the range of MOFs for which data describing molecular diffusion is available. We discuss the prospects for exploiting molecular transport properties in MOFs in practical separation processes and the future role of MD simulations in screening families of MOFs for these processes.

Journal Article↗
PubMed2006

Insights into photodissociation dynamics of propionyl chloride from ab initio calculations and molecular dynamics simulations.

The potential energy surfaces of isomerization, dissociation, and elimination reactions for CH3CH2COCl in the S0 and S1 states have been mapped with the different ab initio calculations. Mechanistic photodissociation of CH3CH2COCl at 266 nm has been characterized through the computed potential energy surfaces, the optimized surface crossing structure, intrinsic reaction coordinate, and ab initio molecular dynamics calculations. Photoexcitation at 266 nm leads to the CH3CH2COCl molecules in the S1 state. From this state, the C-Cl bond cleavage proceeds in a time scale of picosecond in the gas phase. The barrier to the C-Cl bond cleavage on the S1 surface is significantly increased by effects of the matrix and the internal conversion to the ground state prevails in the condensed phase. The HCl eliminations as a result of internal conversion to the ground state become the dominant channel upon photodissociation of CH3CH2COCl in the argon matrix at 10 K.

Journal Article↗
PubMed2006

Molecular dynamics simulations of folding processes of a beta-hairpin in an implicit solvent.

Computer simulations of beta-hairpin folding are relatively difficult, especially those based on the explicit water model. This greatly limits the complete analysis and understanding of their folding mechanisms. In this paper, we use the generalized Born/solvent accessible implicit solvent model to simulate the folding processes of a nine-residue beta-hairpin. We find that the beta-hairpin can fold into its native structure very easily, even using the traditional molecular dynamics method. This allows us to extract 21 complete folding events and investigate the folding process sufficiently. Our results show that there exist four most stable states on the free energy landscape of the short peptide, one native state and three intermediates. We find that two of the non-native stable states have almost the same potential energy as the native state but with lower entropy. This suggests that the native state can be stabilized entropically. Furthermore, we find that the folding processes of this peptide have common features: to fold into its native state, the peptide undergoes a continuous collapsing-extending-recollapsing process to adjust the positions of the side chains in order to form the native middle inter-strand hydrogen bonds. The formations of these bonds are the key step of the folding process. Once these bonds are formed, the peptide can fold into the native state quickly.

Computer Simulation↗
PubMed2000

Magnitude of structural changes of the T-cell receptor binding regions determine the strength of T-cell antagonism: molecular dynamics simulations of HLA-DR4 (DRB1*0405) complexed with analogue peptide.

In our model system, we generated T cell clones specific for the HLA-DR4 (DRB1*0405)-index peptide (YWALEAAAD) complex. Based on response patterns of the T cell clones, analogue peptides containing single amino acid substitutions of the index peptide were classified into three types, agonists, antagonists or null peptides (non-agonistic and non-antagonistic peptides). Subtle structural changes induced by the antagonists in the T-cell receptor (TCR) binding regions have already been explained using the root mean square (r.m.s.) deviations from the DR4-index peptide complex in the molecular dynamics (MD) trajectory. In this work, we performed additional MD simulations at 300 K with explicit solvent molecules to reveal the structural character of the HLA-DR4 complexed with the analogue peptides. We examined the r.m.s. deviations of the TCR-binding sites and the exposed areas of the bound peptides. Remarkable differences of the r.m.s. deviations among the DR4-antagonist complexes, together with our previous data, suggest that the magnitude of structural changes of TCR-binding regions would determine the strength of TCR antagonism. The simulations also indicate that TCR could discriminate null peptides from other ligands mainly through the changes of exposed side chains of the bound peptide, rather than the conformational changes of TCR-binding surfaces on HLA molecule.

Binding Sites↗
PubMed2005

Molecular dynamics simulations of alkylsilane stationary-phase order and disorder. 1. Effects of surface coverage and bonding chemistry.

"Shape-selective" polymeric alkylsilane stationary phases are routinely employed over the more common monomeric phases in reversed-phase liquid chromatography (RPLC) to improve the separation of geometric isomers of shape-constrained solutes. We have investigated the molecular dynamics of chromatographic models that represent both monomeric and polymeric stationary phases with alkylsilane surface coverages and bonding chemistries typical of actual materials in an effort to elucidate the molecular-level structural features that control shape-selective separations. The structural characterization of these models is consistent with previous experimental observations of alkyl chain order and disorder: (1) alkyl chain order increases with increased surface coverage; and (2) monomeric and polymeric phases with similar surface coverages yield similar alkyl chain order (although subtle differences exist). In addition, a significant portion of the alkyl chain proximal to the silica surface is disordered (primarily gauche conformations) and the distal end is most ordered. Models that represent shape-selective RPLC phases possess a significant region of distal end chain order with primarily trans dihedral angle conformations. This is consistent with the view that the alkyl chains comprising polymeric stationary phases contain a series of well-defined and rigid voids in which shape-constrained solutes can penetrate and hence be selectively retained.

Journal Article↗