Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD 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 433 records · Page 24Linked to original sources

Molecular dynamics study of the molecular weight dependence of surface tensions of normal alkanes and methyl methacrylate oligomers.

Surface tensions (gamma) of normal alkanes and methyl methacrylate (MMA) oligomers at various molecular weights in the low molecular weight range were computed using a newly proposed molecular dynamics (MD) simulation strategy which was developed based on the definition of gamma = ( partial differential U/ partial differential sigma)n,V,S. The MD simulations, even with the use of a generic force field, reproduced the experimentally observed molecular weight dependence of gamma (i.e., gamma proportional Mn(-2/3), where Mn is the number-average molecular weight) for both series of oligomers. Analysis of the data reveals that solvent accessible surface area, one of the key input variables used for the calculation of gamma, exhibits an Mn(2/3) (rather than Mn(1)) dependence. The reason for such dependence is that solvent accessible surface area formed by the chainlike small molecules depends, to a larger extent, on their orientations rather than their size. However, this is not the case for high molecular weight molecules as solvent accessible surface area of such surfaces are determined by the orientations of their segments which are determined by the conformations of the molecules. This may explain why surface tension of polymers experimentally exhibits an Mn(-1) dependence. It is inferred that the corresponding molecular weight dependence of the entropy changes associated with molecules in the low and high molecular weight ranges would be different.

Journal Article↗

Importance of discriminator base stacking interactions: molecular dynamics analysis of A73 microhelix(Ala) variants.

Transfer of alanine from Escherichia coli alanyl-tRNA synthetase (AlaRS) to RNA minihelices that mimic the amino acid acceptor stem of tRNA(Ala) has been shown, by analysis of variant minihelix aminoacylation activities, to involve a transition state sensitive to changes in the 'discriminator' base at position 73. Solution NMR has indicated that this single-stranded nucleotide is predominantly stacked onto G1 of the first base pair of the alanine acceptor stem helix. We report the activity of a new variant with the adenine at position 73 substituted by its non-polar isostere 4-methylindole (M). Despite lacking N7, this analog is well tolerated by AlaRS. Molecular dynamics (MD) simulations show that the M substitution improves position 73 base stacking over G1, as measured by a stacking lifetime analysis. Additional MD simulations of wild-type microhelix(Ala) and six variants reveal a positive correlation between N73 base stacking propensity over G1 and aminoacylation activity. For the two DeltaN7 variants simulated we found that the propensity to stack over G1 was similar to the analogous variants that contain N7 and we conclude that the decrease in aminoacylation efficiency observed upon deletion of N7 is likely due to loss of a direct stabilizing interaction with the synthetase.

Acylation↗

EXAFS and molecular dynamics studies of ionic solutions.

This paper focuses on recent advances in the X-ray absorption spectroscopy (XAS) analysis of ionic solutions. The asymmetry of radial distribution functions g(r) associated with the solvent molecules surrounding the ions has to be taken into account to perform a reliable structural analysis. Molecular dynamics (MD) simulations provide reliable g(r)'s which can be used as starting models in the XAS data analysis. The combined MD-XAS investigation reduces meaningfully the indetermination of the structural parameters, especially for coordination numbers and Debye-Waller factors. Double-electron excitation channels can be present in the XAS spectra of ionic solutions and they have to be accounted for in the background extraction. The ability of the XAS technique to probe three-body correlation functions in ionic solutions with the aid of MD g(r1, r2, theta) has been shown. The analysis of the low-k region of the spectra allows the detection of a weak but significant hydrogen structural signal. The XAS technique is especially well suited to determine the detailed shape of the nearest-neighbor peak in the atom-atom pair correlation functions of disordered systems. The information that they contain about the short-range atom-atom pairwise interactions can be very helpful for specifying and properly modifying model potentials used in MD simulations.

Journal Article↗

Molecular dynamics simulations of high-mannose oligosaccharides.

Conformations of several high-mannose-type oligosaccharides that are generated during the biosynthetic degradation of Man9GlcNAc2 to Man5GlcNAc2 have been studied by molecular dynamics (MD). Simulations were performed on NCI-FCRDC's Cray Y-MP 8D/8128 supercomputer using Biosym's CVFF force field for 1000 ps with different initial conformations. The conformations of the two alpha 1,3- and the two alpha 1,6-linkages in each oligomannose were different, suggesting that deriving oligosaccharide conformations based on the conformational preferences of the constituent disaccharide fragments will not always yield correct results. Unlike other oligomannoses, Man9GlcNAc2 appears to take more than one distinct conformation around the core alpha 1,6-linkage. These various conformations may play an important role in determining the processing pathways. Using the data on the preferred conformations of these oligomannoses and the available experimental results, possible pathways for processing Man9GlcNAc2 to Man5GlcNAc2 by alpha 1,2-linkage-specific mannosidases have been proposed. Conformational analysis of Man5GlcNAc2 indicates that the addition of beta 1,2-GlcNAc to the alpha 1,3-linked core mannose, besides serving as a prerequisite for mannosidase II action as suggested earlier, may also prevent the removal of alpha 1,3-mannose. The MD simulations also suggest that the processing of the precursor oligosaccharide during Asn-linked complex and hybrid glycan biosynthesis proceeds in a well-defined pathway involving more than one alpha 1,2-linkage-specific mannosidase. Knowledge of the conformation of the processing intermediates obtained from the present study can be used to design highly specific substrate analogues to inhibit a particular mannosidase, thereby blocking one processing pathway without interfering with the others.

Animals↗

Nucleotide-dependent conformational changes in HisP: molecular dynamics simulations of an ABC transporter nucleotide-binding domain.

ATP-binding cassette (ABC) transporters mediate the movement of molecules across cell membranes in both prokaryotes and eukaryotes. In ABC transporters, solute translocation occurs after ATP is either bound or hydrolyzed at the intracellular nucleotide-binding domains (NBDs). Molecular dynamics (MD) simulations have been employed to study the interactions of nucleotide with NBD. The results of extended (approximately 20 ns) MD simulations of HisP (total simulation time approximately 80 ns), the NBD of the histidine transporter HisQMP2J from Salmonella typhimurium, are presented. Analysis of the MD trajectories reveals conformational changes within HisP that are dependent on the presence of ATP in the binding pocket of the protein, and are sensitive to the presence/absence of Mg ions bound to the ATP. These changes are predominantly confined to the alpha-helical subdomain of HisP. Specifically there is a rotation of three alpha-helices within the subdomain, and a movement of the signature sequence toward the bound nucleotide. In addition, there is considerable conformational flexibility in a conserved glutamine-containing loop, which is situated at the interface between the alpha-helical subdomain and the F1-like subdomain. These results support the mechanism for ATP-induced conformational transitions derived from the crystal structures of other NBDs.

ATP-Binding Cassette Transporters↗

Hydrogen bonding in supercritical tert-butanol assessed by vibrational spectroscopies and molecular-dynamics simulations.

We have investigated the state of aggregation in supercritical tert-butanol (T = 523 K,0.05 < rho < 0.4 g cm(-3)) by means of vibrational spectroscopies (infrared and Raman) and molecular-dynamics (MD) simulations. A quantitative band shape analysis of the spectra associated with the OH stretching mode of tert-butanol has been done using activities computed by ab initio calculations on small clusters. This allows us to determine the degree of hydrogen bonding and populations of oligomers. These latter quantities have been derived from MD simulations and very consistent results are found with experiments. These results show that hydrogen bond still exist in supercritical tert-butanol and that the fluid mainly consists of oligomers smaller than tetramers.

Journal Article↗

Molecular dynamics simulation of hydrated phospholipid bilayers.

Understanding of microscopic behaviour of biological membrane is crucial for designing of molecules to control transport properties of the membranes. Phospholipid-water forms a good model system to study ligand induced structural and dynamical changes in membrane. The review has its main focus on molecular dynamics (MD) simulation of phospholipid bilayers. A brief summary of the current status of structure of phospholipid membranes based on different physico-chemical measurements is given. We discuss here mainly results of MD simulations in the recent years on hydrated phospholipid bilayers and their interaction with ligands. Simulation parameters as: choice of initial system, force fields, protocols for simulation are compared. Main results on: order parameter, head group and chain conformation, water penetration profile, chain tilts, pair-correlation function between atoms of lipid and water, diffusion of ions and ligands are discussed. The review gives application and limitation of MD method for studying lipid water system.

Computer Simulation↗

Development of a force field for Li(2)SiF(6).

A force field has been developed for Li(2)SiF(6) for subsequent use in Molecular Dynamics (MD) simulations involving Li(+) and SiF(2-) (6) ions in a polymer electrolyte host. Both ab initio calculations and available empirical data have been used. The force field has been verified in simulations of the crystal structure of Li(2)SiF(6) in two different space groups: P321 and P3(-)m1. The use of MD simulation to assess the correct space group for Li(2)SiF(6) shows that it is probably P321.

Journal Article↗

Structure and properties of the dendron-encapsulated polyoxometalate (C(52)H(60)NO(12))(12)[(Mn(H(2)O))(3)(SbW(9)O(33))(2)], a first generation dendrizyme.

Combining analytical and theoretical methods, we present a detailed study of a heteropolytungstate cluster encapsulated in a shell of dendritically branching surfactants, namely (C(52)H(60)NO(12))(12)[(Mn(H(2)O))(3)(SbW(9)O(33))(2)], 3. This novel surfactant-encapsulated cluster (SEC) self-assembles spontaneously from polyoxometalate-containing solutions treated with a stoichiometric amount of dendrons. Compound 3 exhibits a discrete supramolecular architecture in which a single polyoxometalate anion resides in a compact shell of dendrons. Our approach attempts to combine the catalytic activity of polyoxometalates with the steric properties of tailored dendritic surfactants into size-selective catalytic systems. The structural characterization of the SEC is based on analytical ultracentrifugation (AUC) and small-angle neutron scattering (SANS). The packing arrangement of dendrons at the cluster surface is gleaned from molecular dynamics (MD) simulations, which suggests a highly porous shell structure due to the dynamic formation of internal clefts and cavities. From analysis of the MD trajectory of 3, a theoretical neutron-scattering function is derived that is in good agreement with experimental SANS data. Force field parameters used in MD simulations are partially derived from a quantum mechanical geometry optimization of [(Zn(H(2)O))(3)(SbW(9)O(33))(2)](12)(-), 2b, at the density functional theory (DFT) level. DFT calculations are corroborated by X-ray structure analysis of Na(6)K(6)[(Zn(H(2)O))(3)(SbW(9)O(33))(2)].23H(2)O, which is isostructural with the catalytically active Mn derivative 2a. The combined use of theoretical and analytical methods aims at rapidly prototyping smart catalysts ("dendrizymes"), which are structurally related to naturally occurring metalloproteins.

Journal Article↗

A computational investigation of the different intermediates during organoalkoxysilane hydrolysis.

Using a combination of atomistic molecular dynamics (MD) simulations and density functional theory (DFT) calculations, the four steps of hydrolysis of aminopropyl-, thiolpropyl-, and butyltrimethoxysilane have been studied. Large box MD simulations at constant pressure and temperature yield appropriate pair distribution functions--which allows us to quantify the number of surrounding water molecules--as well as the density of the systems. These densities serve as input for small box DFT calculations, which allow further geometry optimization and calculation of the electronic structure of the systems. The periodic DFT calculations are compared with gas-phase simulations. In all cases, the first step of hydrolysis is exothermic with the extent depending on the type of silane as well as on the number of hydrogen bonds in the initial stage.

Journal Article↗

Insights into the phosphoryl-transfer mechanism of cAMP-dependent protein kinase from quantum chemical calculations and molecular dynamics simulations.

To investigate the molecular details of the phosphoryl-transfer mechanism catalyzed by cAMP-dependent protein kinase, we performed quantum mechanical (QM) calculations on a cluster model of the active site and molecular dynamics (MD) simulations of a ternary complex of the protein with Mg(2)ATP and a 20-residue peptide substrate. Overall, our theoretical results confirm the participation of the conserved aspartic acid, Asp(166), as an acid/base catalyst in the reaction mechanism catalyzed by protein kinases. The MD simulation shows that the contact between the nucleophilic serine side chain and the carboxylate group of Asp(166) is short and dynamically stable, whereas the QM study indicates that an Asp(166)-assisted pathway is structurally and energetically feasible and is in agreement with previous experimental results.

Catalytic Domain↗

Entropic folding pathway of human epidermal growth factor explored by disulfide scrambling and amplified collective motion simulations.

Epidermal growth factor (EGF) regulates cell proliferation and differentiation by binding to the EGF receptor (EGFR) extra-cellular domains. Human EGF is a small, single-chain protein comprising three distinct loops (A, B, and C), which are connected by three disulfide bridges (Cys6-Cys20, Cys14-Cys31, and Cys33-Cys42). These disulfide bridges are essential for structural stability and biological activity. EGF was extensively studied by disulfide scrambling, an experimental technique for the conformational entrapment of intermediate states, which allows us to study the folding pathway of proteins containing disulfide bonds. The experimental results showed that there is a major 2-disulfide intermediate (denoted EGF-II) and that the native disulfide bonding pattern is less prevalent in one of the mutants. In this article, we investigated for the first time the solution conformations of wild-type EGF, EGF-II, and the mutant S9C through extensive molecular dynamics (MD) simulations in water using both the standard MD technique and a recently developed amplified-collective-motion (ACM) sampling method. Compared to standard MD simulations, we achieved a much more enhanced sampling by the ACM simulations, and the structures were sufficiently relaxed to estimate configurational entropies. The simulation results suggest a predominantly entropic folding pathway governed by the disorder of three functional loop regions. Although EGF-II exhibits two native disulfide bonds (Cys14-Cys31 and Cys33- Cys42), its large configurational entropy inhibits a direct transition to the native structure in the folding process. When Ser9 is mutated into Cys, a non-native disulfide bridge Cys9- Cys20 is slightly more favorable than the native Cys6-Cys20 because a less constrained N-terminus affords larger entropy. Isomers that are functionally less active also exhibit a more localized dynamics of the functional loop regions, which may suggest a possible mechanism for the modulation of EGF activity.

Amino Acid Sequence↗

Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments.

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.

1,2-Dipalmitoylphosphatidylcholine↗

Molecular dynamics and circular dichroism studies of human and rat C-peptides.

Proinsulin C-peptide has been recently described as an endogenous peptide hormone, responsible for important physiological functions others than its role in proinsulin processing. Accumulating evidences that C-peptide exerts beneficial effects in the treatment of long term complications of patients with type 1 diabetes mellitus indicate that this molecule may be administered together with insulin in future therapies. Despite its clear pharmacological interest, the secondary and three-dimensional (3D) structures of human C-peptide are still points of controversy. In the present work we report molecular dynamics (MD) simulations of human, rat I and rat II C-peptides. A common experimental strategy applied to all peptides consisted of homology building followed by multinanosecond MD simulations in vacuum and water. Circular dichroism (CD) experiments of each peptide in the absence and presence of 2,2,2-trifluoroethanol (TFE) were performed to support validation of the theoretical models. A multiple sequence alignment of 23 known mammalian C-peptides was constructed to identify significant conserved sites that would be important for the maintenance of secondary and tertiary structures. The analysis of the molecular dynamics trajectories for the human, rat I and rat II molecules have shown quite different general behavior, being the human C-peptide more flexible than the two others. Human and rat C-peptides exhibit very stable turn-like structures at the middle and C-terminal regions, which have been described as potential active sites of C-peptides. Human C-peptide also presented a short alpha-helix throughout the MD, which was not found in the rat molecules. CD data is in very good agreement with the MD results and both methods were able to identify a greater structural stability and potential in rat C-peptides when compared to the human C-peptide. The simulation results are discussed and validated in the light of multiple sequence alignment, recent experimental data from the literature and our own CD experiments.

Amino Acid Sequence↗

Molecular dynamics simulation of ratchet motion in an asymmetric nanochannel.

The persistence of ratchet effects, i.e., nonzero mass flux under a zero-mean time-dependent drive, when many-body interactions are present, is studied via molecular dynamics (MD) simulations of a simple liquid flowing in an asymmetric nanopore. The results show that (i) ratchet effects persist under many-body density correlations induced by the forcing; (ii) two distinct linear responses (flux proportional to the drive amplitude) appear under strong loads. One regime has the same conductivity of linear response theory up to a forcing of about 10 kT, while the second displays a smaller conductivity, the difference in responses is due to geometric effects alone. (iii) Langevin simulations based on a naive mapping of the many-body equilibrium bulk diffusivity, D, onto the damping rate, gamma are also found to yield two distinct linear responses. However, in both regimes, the flux is significantly smaller than the one of MD simulations.

Computer Simulation↗

Effect of N-terminal region of eIF4E and Ser65-phosphorylation of 4E-BP1 on interaction between eIF4E and 4E-BP1 fragment peptide.

To clarify the contribution of N-terminal region of eukaryotic initiation factor 4E (eIF4E) to the interaction with 4E-BP and to investigate the effect of 4E-BP phosphorylation on the interaction with eIF4E, the interaction profiles of the Ser65-unphosphorylated and phosphorylated peptides (Thr37-Thr70 fragment of 4E-BP1) with full-length and N-terminal 33 residues-deleted eIF4Es were investigated by fluorescence and SPR methods. The effect of N-terminal region of eIF4E on the interaction with 4E-BP1 peptides was shown to be dependent on the interaction state, that is, the steady-state fluorescence and kinetic-state SRP analyses showed the positive and negative contributions of the N-terminal region to the interaction with the peptide, respectively, despite its unphosphorylated or phosphorylated state. The comparison of the association constants of the peptide with those of full-length 4E-BP1 indicated the importance of N-terminal (1-36) and/or C-terminal (71-118) sequence of 4E-BP1 for the interaction, although the MD simulations suggested that the alpha-helical region (Arg56-Cys62) of 4E-BP1 peptide is sufficient for keeping the interaction. The MD simulations also indicated that a charge-dependent rigid hydration shell formed around the phosphate group makes the molecular conformation rigid, and single Ser65 phosphorylation is insufficient for releasing 4E-PB1 peptide from eIF4E.

Adaptor Proteins, Signal Transducing↗

HIV-1 protease molecular dynamics of a wild-type and of the V82F/I84V mutant: possible contributions to drug resistance and a potential new target site for drugs.

The protease from type 1 human immunodeficiency virus (HIV-1) is a critical drug target against which many therapeutically useful inhibitors have been developed; however, the set of viral strains in the population has been shifting to become more drug-resistant. Because indirect effects are contributing to drug resistance, an examination of the dynamic structures of a wild-type and a mutant could be insightful. Consequently, this study examined structural properties sampled during 22 nsec, all atom molecular dynamics (MD) simulations (in explicit water) of both a wild-type and the drug-resistant V82F/I84V mutant of HIV-1 protease. The V82F/I84V mutation significantly decreases the binding affinity of all HIV-1 protease inhibitors currently used clinically. Simulations have shown that the curling of the tips of the active site flaps immediately results in flap opening. In the 22-nsec MD simulations presented here, more frequent and more rapid curling of the mutant's active site flap tips was observed. The mutant protease's flaps also opened farther than the wild-type's flaps did and displayed more flexibility. This suggests that the effect of the mutations on the equilibrium between the semiopen and closed conformations could be one aspect of the mechanism of drug resistance for this mutant. In addition, correlated fluctuations in the active site and periphery were noted that point to a possible binding site for allosteric inhibitors.

Amino Acid Substitution↗

Probing the instabilities in the dynamics of helical fragments from mouse PrPC.

The first step in the formation of the protease resistant form (PrPSc) of prion proteins involves a conformational transition of the monomeric cellular form of PrPC to a more stable aggregation prone state PrPC*. A search of PDBselect and Escherichia coli and yeast genomes shows that the exact pattern of charges in helix 1 (H1) is rare. Among the 23 fragments in PDBselect with the pattern of charges that match H1, 83% are helical. Mapping of the rarely found (in E. coli and yeast genomes) hydrophobicity patterns in helix 2 (H2) to known secondary structures suggests that the PrPC-->PrPC* transition must be accompanied by alterations in conformations in second half of H2. We probe the dynamical instability in H1 and in the combined fragments of H2 and helix 3 (H3) from mPrPC (H2+H3), with intact disulfide bond, using all atom molecular dynamics (MD) simulations totaling 680 ns. In accord with recent experiments, we found that H1 is helical, whereas the double mutant H1[D147A-R151A] is less stable, implying that H1 is stabilized by the (i,i + 4) charged residues. The stability of H1 suggests that it is unlikely to be involved in the PrPC-->PrPC* transition. MD simulations of H2+H3 shows that the second half of H2 (residues 184-194) and parts of H3 (residues 200-204 and 215-223) undergo a transition from alpha-helical conformation to a beta and/or random coil state. Simulations using two force fields (optimized potentials for liquid simulations and CHARMM) give qualitatively similar results. We use the MD results to propose tentative structures for the PrPC* state.

Animals↗