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 811 records · Page 45Linked to original sources
PubMed2001

Signatures of beta-peptide unfolding in two-dimensional vibrational echo spectroscopy: a simulation study.

An ensemble of exciton Hamiltonians for the amide-I band of the folded and unfolded states of a helical beta-heptapeptide is generated using a molecular dynamics (MD) simulation. The correlated fluctuations of its parameters and their signatures in two-dimensional (2D) vibrational echo spectroscopy are computed. This technique uses infrared pulse sequences to provide ultrafast snapshots of molecular structural fluctuations, in analogy with multidimensional NMR. The present study demonstrates that, by combining a method of calculating the vibrational Hamiltonian from MD snapshots and the nonlinear exciton equations (NEE), it may be possible to simulate realistic multidimensional IR spectra of chemically and biologically interesting systems.

Computer Simulation↗
PubMed2003

Self-consistent-field modeling of hydrated unsaturated lipid bilayers in the liquid-crystal phase and comparison to molecular dynamics simulations.

A molecular-level self-consistent-field (SCF) theory is applied to model the lipid bilayer structures composed of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (18:0/18:1 omega 9cis PC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (18:0/22:6 omega 3cis PC). As compared to earlier attempts to model (saturated) PC membranes several additional features are implemented: (i) A water model is used which correctly leads to low water concentration in the bilayers. (ii) Free volume is allowed for, which is important to obtain bilayers in the fluid state. (iii) A polarization term is included in the segment potentials; this new feature corrects for a minor thermodynamic inconsistency present in (all) earlier results for charged bilayers. (iv) The CH3 groups in the lipid molecules are assumed to have twice the volume of a CH2 group; this leads to stable noninterdigitated bilayers. (v) A cis double bond is simulated by forcing gauche conformations along the sn-2 acyl chain. Results of an all-atom molecular dynamics (MD) simulation, using the collision dynamics method, on the same system are presented. Both SCF and MD prove, in accordance with experimental facts, that acyl unsaturation effectively reduces the length of the chain which counteracts interdigitation. It is also found that the phosphatidylcholine head group is lying almost flat on the membrane surface and the water penetrates into the bilayer upto the glycerol backbone units. From the SCF results it further followed that the free volume is not exactly evenly distributed over the bilayer. There is a small increase in free volume in the center of the bilayer as well as in the glycerol backbone region.

Journal Article↗
PubMed1997

500 picosecond molecular dynamics simulation of amphiphilic polypeptide Ac(LKKL)4 NHEt with 1,2 di-mysristoyl-sn-glycero-3-phosphorylcholine (DMPC) molecules.

Molecular dynamics (MD) simulation of the interaction between amphiphilic polypeptide Ac(LKKL)4NHEt and 4 DMPC (1,2 di-mysristoyl-sn-glycero-3-phosphorylcholine) molecules has been carried out at 310 K for 500 picoseconds (ps) using AMBER 4.0. Interaction energy and a number of conformational parameters are calculated for the subaveraged coordinates, using P-CURVES 3.1 and our MD trajectory analysis program ANALMD. No significant change in DMPC headgroup conformation was observed. However, the mobility of P atoms was found to be restricted. The chains were quite flexible and their flexibility increased towards the ends. They interacted amongst themselves. The polypeptide remained predominantly in alpha-helical conformation. Leu1 and Lys2 at the N terminus and Leu13 to Leu16 at C terminus assumed non helical conformation and were quite flexible. Average interaction energy between the polypeptide and DMPC molecules was found to be -151.828 kcal*mol-1. The main contributory factor was electrostatic interaction of Lys NH3+ groups with the DMPC phosphates. On an average one Lys chain interacted with 1.5 DMPC molecules. Central region of the polypeptide had better contact with DMPC molecules. A model for the fusogenic properties of the polypeptide is presented on the basis of MD results.

Amino Acid Sequence↗
PubMed2001

Molecular simulation of dioleoylphosphatidylcholine lipid bilayers at differing levels of hydration.

The structure and dynamics of the lipid and water components of dioleoylphosphatidylcholine bilayers at various levels of hydration were studied using molecular dynamics (MD) simulations. Equilibration of these systems proceeded by use of a hybrid MD and configurational-bias Monte Carlo technique using one atmosphere of pressure normal to the membrane and a set point for the lateral area derived from experimental Bragg spacings, combined with experimentally derived specific volumes for each of the system components. Membrane surface tensions were observed to be of the order of tens of dyn/cm. The transbilayer molecular fragment peak positions at low hydration were found to agree with experimental neutron and x-ray scattering profiles and previously published simulations. For hydration levels of 5.4, 11.4, and 16 waters/lipid, molecular fragment distributions and order parameters for the headgroup, lipid chains, and water were quantified. Spin-lattice relaxation rates and lateral self-diffusion coefficients of water agreed well with results from experimental nuclear magnetic resonance studies. Relaxation rates of the choline segments and chemical shift anisotropies for the phosphate and carbonyls were computed. Headgroup orientation, as measured by the P-N vector, showed enhanced alignment with the membrane surface at low hydration. The sign of the membrane dipole potential reversed at low hydration, with the membrane interior negative relative to the interlamellar region. Calculation of the number of water molecules in the headgroup hydration shell, as a function of hydration level, supports the hypothesis that the break point in the curve of Bragg spacing versus hydration level near 12 waters/lipid, observed experimentally by Hristova and White (1988. Biophys. J. 74:2419-2433), marks the completion of the first hydration shell.

Anisotropy↗
PubMed2005

Coexistence of multiple conformations in cysteamine monolayers on Au(111).

The structural organization, catalytic function, and electronic properties of cysteamine monolayers on Au(111) have been addressed comprehensively by voltammetry, in situ scanning tunneling microscopy (STM) in anaerobic environment, and a priori molecular dynamics (MD) simulation and STM image simulation. Two sets of voltammetric signals are observed. One peak at -(0.65-0.70) V (SCE) is caused by reductive desorption of cysteamine. The other signal, at -(0.25-0.40) V consists of a peak doublet. The pH dependence of the latter suggests that the origin is catalytic dihydrogen evolution. The doublet feature is indicative of two distinct cysteamine configurations. Cysteamine monolayer formation from initial nucleation to a highly ordered phase has been successfully observed in real time using oxygen-free in situ STM. Random cellular patterns, disordered adlayer formation accompanied by high step edge mobility, and ultimately a highly ordered (square root 3 x 4) R30 degrees lattice are observed sequentially. Pits are formed due to enclosure of the mobile edges during the adsorption process. In the highly ordered cysteamine layer, each unit has two spots with apparent 0.6 A height difference in STM images. The coverage 5.7 +/- 0.1 x 10(-10) mol cm(-2) determined by voltammetry supports that the spots represent two individual cysteamine molecules. A priori MD and density functional simulations hold other clues to the image interpretation and indicate that the NH(3)(+) groups dominate the tunneling contrast. A wide range of interface structures, showing variations in the sulfur binding site and orientation, gauche and trans conformers, and especially hydrogen-bonding interactions, are examined, from which it is concluded that the adsorbate structure is controlled by interactions with the solvent rather than with the substrate.

Computer Simulation↗
PubMed2002

Molecular dynamics simulation studies of induced fit and conformational capture in U1A-RNA binding: do molecular substates code for specificity?

Molecular dynamics (MD) simulations on stem loop 2 of U1 small nuclear RNA and a construct of the U1A protein were carried out to obtain predictions of the structures for the unbound forms in solution and to elucidate dynamical aspects of induced fit upon binding. A crystal structure of the complex between the U1A protein and stem loop 2 RNA and an NMR structure for the uncomplexed form of the U1A protein are available from Oubridge et al. (Nature, 1994, Vol. 372, pp. 432-438) and Avis et al. (Journal of Molecular Biology, 1996, Vol. 257, pp. 398-411), respectively. As a consequence, U1A-RNA binding is a particularly attractive case for investigations of induced fit in protein-nucleic acid complexation. When combined with the available structural data, the results from simulations indicate that structural adaptation of U1A protein and RNA define distinct mechanisms for induced fit. For the protein, the calculations indicate that induced fit upon binding involves a non-native thermodynamic substate in which the structure is preorganized for binding. In contrast, induced fit of the RNA involves a distortion of the native structure in solution to an unstable form. However, the RNA solution structures predicted from simulation show evidence that structures in which groups of bases are favorably oriented for binding the U1A protein are thermally accessible. These results, which quantify with computational modeling recent proposals on induced fit and conformational capture by Leuillot and Varani (Biochemistry, 2001, Vol. 40, pp. 7947-7956) and by Williamson (Nature Structural Biology, 2000, Vol. 7, pp. 834-837) suggest an important role for intrinsic molecular architecture and substates other than the native form in the specificity of protein-RNA interactions.

Computer Simulation↗
PubMed2006

Speciation of the curium(III) ion in aqueous solution: a combined study by quantum chemistry and molecular dynamics simulation.

The structures of aquo complexes of the curium(III) ion have been systematically studied using quantum chemical and molecular dynamics (MD) methods. The first hydration shell of the Cm3+ ion has been calculated using density functional theory (DFT), with and without inclusion of the conductor-like polarizable continuum medium (CPCM) model of solvation. The calculated results indicate that the primary hydration number of Cm3+ is nine, with a Cm-O bond distance of 2.47-2.48 A. The calculated bond distances and the hydration number are in excellent agreement with available experimental data. The inclusion of a complete second hydration shell of Cm3+ has been investigated using both DFT and MD methods. The presence of the second hydration shell has significant effects on the primary coordination sphere, suggesting that the explicit inclusion of second-shell effects is important for understanding the nature of the first shell. The calculated results indicate that 21 water molecules can be coordinated in the second hydration shell of the Cm3+ ion. MD simulations within the hydrated-ion model suggest that the second-shell water molecules exchange with the bulk solvent with a lifetime of 161 ps.

Journal Article↗
PubMed2000

Simulations of the estrogen receptor ligand-binding domain: affinity of natural ligands and xenoestrogens.

We have carried out molecular dynamics (MD) simulations and free energy calculations on the alpha-subtype of the human estrogen receptor ligand-binding domain (ERalpha LBD) complexed with a number of known agonists and putative xenoestrogens. Our dynamical simulations of ligand-receptor complexes underscore the highly structured nature of the complex and offer some interesting insights into the structure-activity relationship (SAR) for these ligands. With traditional thermodynamic integration (TI) calculations, we calculate relative binding free energies for three known agonists, in good agreement with experimental values. The sheer number of possible xenoestrogenic compounds makes an approach using traditional free energy calculations unfeasible. Instead, we have made use of a single-step perturbation methodology that allows the calculation of relative free energies for a large number of related polyaromatic hydrocarbons (PAHs) from a single simulation. Our results show good (maximum deviation 3.3 kJ mol(-1)) agreement with experimental data, suggesting the possibility of large-scale xenoestrogen screening in silico to obtain strongly estrogenic compounds for subsequent experimental testing.

Diethylstilbestrol↗
PubMed2005

Simulations of nanopore formation and phosphatidylserine externalization in lipid membranes subjected to a high-intensity, ultrashort electric pulse.

A combined MD simulator and time dependent Laplace solver are used to analyze the electrically driven phosphatidylserine externalization process in cells. Time dependent details of nanopore formation at cell membranes in response to a high-intensity (100 kV/cm), ultrashort (10 ns) electric pulse are also probed. Our results show that nanosized pores could typically be formed within about 5 ns. These predictions are in very good agreement with recent experimental data. It is also demonstrated that defect formation and PS externalization in membranes should begin on the anode side. Finally, the simulations confirm that PS externalization is a nanopore facilitated event, rather than the result of molecular translocation across the trans-membrane energy barrier.

Computer Simulation↗
PubMed2006

Detailed mechanism for AmtB conducting NH4+/NH3: molecular dynamics simulations.

The mechanism by which the ammonium transporter, AmtB, conducts NH4+/NH3 into the cytoplasm was investigated using conventional molecular dynamics (MD) simulations. These simulations revealed that the neutral molecule, NH3, passes automatically through the channel upon its arrival at the Am2 site and that the function of the channel as a one-way valve for passage of NH3 is not determined by the cytoplasmic exit gate but, rather, by the periplasmic entrance gate of the channel. The NH3, produced by deprotonation of NH4+ at the entrance gate, is spontaneously conveyed to the central region of the channel via a hydrogen-bond network comprising His-168, His-318, Tyr-32, and the NH3 molecule. Finally, the NH3 molecule exits the channel through the exit gate formed by Phe-31, Ile-266, Val-314, and His-318. In addition, Ser-263 is shown to play a critical role in the final stages, acting as a pivoting arm to shunt the NH3 molecule from the cytoplasmic exit gate of the channel out into the cytoplasm. This finding is further supported by another simulation which shows that NH3 fails to be translocated through the channel formed by the Ser-263-Ala mutation. Thus, this study casts new insights on the mechanism of AmtB-mediated passage of NH3 across cellular membranes.

Cation Transport Proteins↗
PubMed1992

QSAR analysis of a set of benzothiopyranoindazole anti-cancer analogs based upon their DNA intercalation properties as determined by molecular dynamics simulation.

The global low-energy intercalation thermodynamic properties for a set of 14 benzothiopyranoindazole (BTPI) analogs, determined by molecular dynamics (MD) simulation, were used as descriptors in a quantitative structure-activity relationship (QSAR) analysis to describe their observed anti-cancer activities. The calculated free-space binding (intercalation) energies by themselves, or in combination with ligand lipophilicity, do not yield significant QSARs. However, the free-space binding energy combined with the aqueous desolvation energy of the intercalation complex leads to a highly significant QSAR. The optimum QSAR corresponds to both the ligand and DNA model being charged, and a molecular dielectric of 40 used in the force field. The desolvation energy is computed using a hydration-shell model. The QSAR suggests that anti-cancer activity can be enhanced by increasing the intercalating free-space binding energy and/or decreasing the loss in the favorable hydration energies of both the ligand and DNA when intercalation occurs.

Antineoplastic Agents↗
PubMed2002

Molecular dynamics simulations of adrenocorticotropin (1-24) peptide in a solvated dodecylphosphocholine (DPC) micelle and in a dimyistoylphosphatidylcholine (DMPC) bilayer.

The structure and interactions of the 1-24 fragment of the adrenocorticotropin hormone, ACTH (1-24), with membrane have been studied by molecular dynamics (MD) simulation in an NPT ensembles in two explicit membrane mimics, a dodecylphosphocholine (DPC) micelle and a dimyristoylphosphatidylcholine (DMPC) bilayer. The starting configuration of the peptide/lipid systems had the 1-10 segment of the peptide lying on the surface of the model membrane, the same as the equilibrated structure (by MD) of ACTH (1-10) in a DPC micelle. The simulations showed that the peptide adopts the surface-binding mode and essentially the same structure in both systems. Thus the results of this work lend support to the assumption that micelles are reasonable mimics for biological membranes for the study of peptide binding. The 1-10 segment is slightly tilted from the parallel orientation to the interface and interacts strongly with the membrane surface while the more polar 11-24 segment shows little tendency to interact with the membrane surface, preferring to reside primarily in the aqueous phase. Furthermore, the 1-10 segment of the peptide binds to the DPC micelle in essentially the same way as ACTH (1-10). Thus the MD results are in excellent agreement with the model of interaction of ACTH (1-24) with membrane derived from NMR experiments. The secondary structure and the hydration of the peptide and the interactions of specific residues with the lipid head groups have also been analyzed.

Adrenocorticotropic Hormone↗
PubMed2004

DNA sequence-dependent deformability--insights from computer simulations.

The article reviews some recent developments in studying DNA sequence-dependent deformability, with emphasis on computer modeling. After a brief outline of available experimental techniques, we proceed to computational methods and focus on atomic-resolution molecular dynamics (MD) simulations. A sequence-dependent local (base-pair step) force field inferred from MD is compared with force fields obtained by other techniques. Various methods for establishing global (flexible-rod) DNA elastic constants are reviewed, including an approach based on atomic resolution MD. The problem of defining the global deformation variables, as well as the question of anisotropy and nonlocal effects, are discussed. As an example, both local and global deformability calculations from atomic-resolution MD of EcoRI dodecamer are presented.

Base Sequence↗
PubMed2003

Characteristic domain motion in the ribosome recycling factor revealed by 15N NMR relaxation experiments and molecular dynamics simulations.

The backbone dynamics of ribosome recycling factor (RRF) from Escherichia coli in water were characterized by (15)N NMR relaxation analysis and molecular dynamics (MD) simulation. RRF is composed of two domains connected by a joint region that consists of two peptide chains, such that the overall structure seems to mimic that of tRNA. MD trajectories indicated that the relative orientation of domains varies on the nanosecond time scale. We analyzed the observed (15)N T(1), T(2), and NOE using an extended model-free spectral density function in which the domain motions with a nanosecond time scale were considered. At 30 degrees C, the order parameters of slow motion () were determined to be approximately 0.9 for domain I and 0.7 for domain II, respectively. These values indicate that domain I is nearly fixed on the molecular diffusion frame, and domain II is wobbling in a cone for which the semi-angle is about 30 degrees.

Escherichia coli↗
PubMed1994

Molecular dynamics simulation of a leucine zipper motif predicted for the integrase of human immunodeficiency virus type 1.

We have used the molecular dynamics (MD) simulation package AMBER4 to search the conformation of a peptide predicted as a leucine zipper motif for the human immunodeficiency virus type 1 integrase protein (HIV IN-LZM). The peptide is composed of 22 amino acid residues and its location is from Val 151 to Leu 172. The searching procedure also includes two known alpha-helices that served as positive controls--namely, a 22-residue GCN4-p1 (LZM) and a 20-residue poly (L-alanine) (PLA). A 21-residue peptide extracted from a cytochrome C crystal (CCC-t) with determined conformation as a beta-turn is also included as a negative control. At the beginning of the search, two starting conformations--namely, the standard right-handed alpha-helix and the fully stretched conformations--are generated for each peptide. Structures generated as standard alpha-helix are equilibrated at room temperature for 90 ps while structures generated as a fully stretched one are equilibrated at 600 K for 120 ps. The CCC-t and PLA helices are nearly destroyed from the beginning of equilibration. However, for both the HIV IN-LZM and the GCN4-p1 LZM structures, there is substantial helicity being retained throughout the entire course of equilibration. Although helix propagation profiles calculated indicate that both peptides possess about the same propensity to form an alpha-helix, the HIV IN-LZM helix appears to be more stable than the GCN4-p1 one as judged by a variety of analyses on both structures generated during the equilibration course. The fact that predicted HIV IN-LZM can exist as an alpha-helix is also supported by the results of high temperature equilibration run on the fully stretched structures generated. In this run, the RMS deviations between the backbone atoms of the structures with the lowest potential energy (PE) identified within every 2 ps and the structure with the lowest PE searched in the same course of simulation are calculated. For both the HIV IN-LZM and the GCN4-p1 LZM, these rms values decrease with the decrease of PE, which indicates that both structures are closer in conformations as their PEs are moved deeper into the PE well.

Amino Acid Sequence↗
PubMed2002

Constant normal pressure, constant surface tension, and constant temperature molecular dynamics simulation of hydrated 1,2-dilignoceroylphosphatidylcholine monolayer.

A constant normal pressure, constant surface tension, and constant temperature (NP(N)gammaT) molecular dynamics (MD) simulation of the liquid condensed phase of a 1,2-dilignoceroylphosphatidylcholine (DLGPC) monolayer has been performed at 293.15 K. A DLGPC molecule has two saturated 24-carbon acyl chains, giving the hydrocarbon core thickness of the monolayer approximately 28 A, which is close to the hydrocarbon core thickness of a membrane of a living system. NP(N)gammaT ensemble was used to reproduce the experimental observations, such as area/lipid, because surface tension is an essential factor in determining the monolayer structure. Data analysis on DLGPC/water monolayer shows that various liquid condensed-phase properties of the monolayer have been well reproduced from the simulation, indicating that surface tension 22.9 mN/M used in the simulation is an appropriate condition for the condensed-phase NP(N)gammaT simulation. The simulation results suggest that this long-chain phospholipid monolayer shares many structural characteristics with typical short-chain 1,2-diacylphosphatidylcholine systems, such as DPPC/water monolayer in the condensed phase and DPPC/water bilayer in the gel phase. Furthermore, it was found that DLGPC/water monolayer has almost completely rotationally disordered acyl chains, which have not been observed so far in short-chain 1,2-diacylphosphatidylcholine/water bilayers. This study indicates the good biological relevance of the DLGPC/water monolayer which might be useful in protein/lipid studies to reveal protein structure and protein/lipid interactions in a membrane environment.

Acylation↗
PubMed2004

Origin of mechanical strength of bovine carbonic anhydrase studied by molecular dynamics simulation.

The forced unfolding process of bovine carbonic anhydrase II (BCA II) was examined at the atomic level by the molecular dynamics (MD) simulation. By force spectroscopy, experimentally obtained force-extension curves (F-E curves) showed a prominent force peak after 50 nm extension. F-E curves obtained from our simulation had three force peaks appearing after extensions of 10-17 nm, 40 nm, and 53 nm, each signifying a brittle fracture of a specific local structure. Upon undergoing the final fracture at 53 nm of extension, the entire molecule became a single flexible chain and was further extended to its full theoretical length, almost as a random coil. This feature of the 53-nm peak strongly suggested its close correspondence to the experimentally observed force peak at approximately 60-nm extension. The 53-nm peak in the molecular dynamics simulation corresponded to the unfolding process of the beta-sheeted core that includes zinc-coordinating histidine residues. These results suggest that the structural change occurring at 50-60 nm in atomic force microscopy experiments corresponded to the destruction of the zinc coordination site.

Animals↗
PubMed2002

Temperature dependence of protein dynamics: computer simulation analysis of neutron scattering properties.

The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.

Biophysical Phenomena↗