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,459 records · Page 81Linked to original sources
PubMedDate not supplied

Analysis of structure-function relationships of neuropeptide Y using molecular dynamics simulations and pharmacological activity and binding measurements.

Studies on the structure-function relationship of neuropeptide Y (NPY) were undertaken using a combination of in vacuo molecular dynamics (MD) simulations and pharmacological receptor binding and biological activity measurements. Following a conformational search of NPY from which a theoretical structure was determined, a study of the structural and dynamic changes in the region of amino acids 25-36 was performed in a variety of NPY fragments and in the NPY free acid. Results revealed an increased structural change as the fragment size was decreased. Also, the mobility appears to be lowest in the full NPY vs the NPY fragments. Pharmacological measurements showed a decreased receptor binding and biological activity as fragment size decreased. Combination of the two approaches suggests a model where conformational maintenance and low configurational entropy of the 25-36 region of NPY favors both receptor binding and biological activity. Furthermore, the possibility of two receptor interaction modes is suggested. Analysis of the NPY structure suggests the direct importance of the amidated C-terminus, Gln34 and His26, an indirect importance of the Tyr1 sidechain as well as the potential importance of an apparent electric 'dipole' in NPY for receptor binding and biological activity.

Amino Acids↗
PubMed2005

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method.

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.

Chemistry, Physical↗
PubMed1991

Mean field stochastic boundary molecular dynamics simulation of a phospholipid in a membrane.

Computer simulations of phospholipid membranes have been carried out by using a combined approach of molecular and stochastic dynamics and a mean field based on the Marcelja model. First, the single-chain mean field simulations of Pastor et al. [(1988) J. Chem. Phys. 89, 1112-1127] were extended to a complete dipalmitoylphosphatidylcholine molecule; a 102-ns Langevin dynamics simulation is presented and compared with experiment. Subsequently, a hexagonally packed seven-lipid array was simulated with Langevin dynamics and a mean field at the boundary and with molecular dynamics (and no mean field) in the center. This hybrid method, mean field stochastic boundary molecular dynamics, reduces bias introduced by the mean field and eliminates the need for periodic boundary conditions. As a result, simulations extending to tens of nanoseconds may be carried out by using a relatively small number of molecules to model the membrane environment. Preliminary results of a 20-ns simulation are reported here. A wide range of motions, including overall reorientation with a nanosecond decay time, is observed in both simulations, and good agreement with NMR, IR, and neutron diffraction data is found.

Computer Simulation↗
PubMed2005

Interfacial folding and membrane insertion of designed peptides studied by molecular dynamics simulations.

The mechanism of interfacial folding and membrane insertion of designed peptides is explored by using an implicit membrane generalized Born model and replica-exchange molecular dynamics. Folding/insertion simulations initiated from fully extended peptide conformations in the aqueous phase, at least 28 A away from the membrane interface, demonstrate a general mechanism for structure formation and insertion (when it occurs). The predominately hydrophobic peptides from the synthetic WALP and TMX series first become localized at the membrane-solvent interface where they form significant helical secondary structure via a helix-turn-helix motif that inserts the central hydrophobic residues into the membrane interior, and then fluctuations occur that provide a persistent helical structure throughout the peptide and it inserts with its N-terminal end moving across the membrane. More specifically, we observed that: (i) the WALP peptides (WALP16, WALP19, and WALP23) spontaneously insert in the membrane as just noted; (ii) TMX-1 also inserts spontaneously after a similar mechanism and forms a transmembrane helix with a population of approximately 50% at 300 K; and (iii) TMX-3 does not insert, but exists in a fluctuating membrane interface-bound form. These findings are in excellent agreement with available experimental data and demonstrate the potential for new implicit solvent/membrane models together with advanced simulation protocols to guide experimental programs in exploring the nature and mechanism of membrane-associated folding and insertion of biologically important peptides.

Amino Acid Motifs↗
PubMed2002

Molecular dynamics simulation for ligand-receptor studies. Carbohydrates interactions in aqueous solutions.

The review deals with the problem of the study of ligand-receptor interactions and the use of Molecular Dynamics (MD) simulation to approach such a problem. After a short review of the fundamentals of MD we describe the medium in which all biology takes place, water. Emphasis is put on the water models appropriate for simulation of macromolecular systems explicitly including the water molecules. We consider the quality of the water model both in terms of simplicity and performance to describe the liquid water properties. Heavy water, although not a biologically viable medium, is considered since many experiments make use of it as a solvent. Sweetness of carbohydrates is considered as an example of the procedure suitable to characterize active sites on the ligands. Consideration is given to the computation of the binding constants through molecular dynamics. The computation of the Free Energy is described and illustrated. The potentiality of MD for studies of ligand-receptor interactions is limited by the computer resources, for even with large computing facilities the need of relatively long simulation times severely restricts the study of large systems. A method is described in which several shells are treated at different levels of approximation, form mechanical response and mean electrical field to quantum mechanics, through stochastic dynamics and atomic classical MD. The review closes with a brief account of the perspectives of the method.

Binding Sites↗
PubMed2005

The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water.

The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280-380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.

Journal Article↗
PubMed2000

Unfolding of hen egg lysozyme by molecular dynamics simulations at 300K: insight into the role of the interdomain interface.

We present the results of two 1.2 ns molecular dynamics (MD) unfolding simulations on hen egg lysozyme in water at 300K, performed using a new procedure called PEDC (Path Exploration With Distance Constraints). This procedure allows exploration of low energy structures as a function of increasing RMSD from the native structure, and offers especially the possibility of extensive exploration of the conformational space during the initial unfolding stages. The two independent MD simulations gave similar chronology of unfolding events: disruption of the active site, kinking of helix C, partial unfolding of the three-stranded beta-sheet to a two-stranded sheet (during which the helices A, B, and D remain to a great extent native), and finally unfolding of the beta-domain and partial unfolding of the alpha-domain in which hydrophobic clusters persist. We show particularly that the loss of hydrophobic contacts between the beta-sheet turn residues Leu55 and Ile56 and the hydrobic patch of the alpha-domain destabilizes the beta-domain and leads to its unfolding, suggesting that the correct embedding of these residues in the alpha-beta interface may constitute the rate limiting step in folding. These results are in accord with experimental observations on the folding/unfolding behavior of hen egg lysozyme at room temperature. They would also explain the loss of stability and the tendency to aggregation observed for the mutant Leu55Thr, and the slow refolding kinetics observed in the analogous amyloidogenic variant of human lysozyme.

Animals↗
PubMed2004

Role of histidine-85 in the catalytic mechanism of thymidine phosphorylase as assessed by targeted molecular dynamics simulations and quantum mechanical calculations.

The structural changes taking place in the enzyme thymidine phosphorylase (TPase, also known as PD-ECGF) that are required to achieve catalytic competence upon binding thymidine and phosphate have been simulated by means of targeted molecular dynamics (tMD). The hinge regions were characterized by structural homology comparisons with pyrimidine nucleoside phosphorylase, whose X-ray structure has been solved both in a closed and in an open form. The rearrangement of residues around the substrate that was observed during the tMD trajectory suggested that His-85 could be playing an important role in the catalytic mechanism. A quantum mechanical study of the reaction in the presence of the most relevant active site residues was then performed at the semiempirical level. The results revealed that His-85 could be involved in the protonation of the pyrimidine base at the O2 position to yield the enol tautomer of the base. To establish the role of this oxygen atom in the reaction, ground states, transition states, and final products were studied using higher level ab initio methods starting from both thymidine and 2-thiothymidine as alternative substrates. Comparison of both transition states showed that replacing the oxygen at position 2 of the pyrimidine base by sulfur should accelerate the reaction rate. Consistent with this result, 2-thiothymidine was shown to be a better substrate for TPase than the natural substrate, thymidine. For simulating the final step of the reaction, tMD simulations were used to study domain opening upon product formation considering both the enol and keto tautomers of thymine. Product release from the enzyme was easiest in the simulation that incorporated the keto tautomer of thymine, suggesting that the enol intermediate spontaneously tautomerizes back to the more energetically stable keto form. These results highlight a previously unreported role for His-85 in the catalytic mechanism of TPase and can have important implications for the design of novel TPase inhibitors.

Amino Acid Sequence↗
PubMed1999

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.

Base Sequence↗
PubMed1993

Molecular dynamics simulations of a winter flounder "antifreeze" polypeptide in aqueous solution.

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute-solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge.

Amino Acid Sequence↗
PubMed2003

Molecular dynamics simulation of the formation, structure, and dynamics of small phospholipid vesicles.

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-)(3) cm s(-)(1), in agreement with experimental measurements.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2004

Molecular dynamics simulations of acetonitrile/dimethyl sulfoxide liquid mixtures.

Binary liquid mixtures of dimethyl sulfoxide and acetonitrile at the three molar fractions 0.25, 0.50, and 0.75 have been investigated by molecular dynamics computer simulations. Thermodynamic states corresponding to liquid-vapor coexistence at a temperature of 298 K were considered. Intermolecular interactions were described by potential models of the site-site (12-6) Lennard-Jones plus Coulomb type that have been developed for the description of the pure liquids. Dimethyl sulfoxide has been represented by four interactions sites and acetonitrile by a three- as well as a six-site potential model. We have calculated thermodynamic properties and the intermolecular pair distribution functions. Intermolecular interaction energies indicate deviations from the behavior of ideal mixtures. The local mole fraction analysis demonstrates that dimethyl sulfoxide is preferentially solvated by acetonitrile and that the first solvation shell surrounding acetonitrile molecules is significantly enriched by dimethyl sulfoxide. The nonideal behavior in the mixtures is not affected by the choice of the three- or the six-site potential model for acetonitrile. Orientational correlations of dipole vectors within the first solvation shells indicate that the relative molecular orientations found in pure acetonitrile and dimethyl sulfoxide are maintained in the mixtures. Parallel and antiparallel dipole-dipole configurations determine first shell acetonitrile-dimethyl sulfoxide configurations. Dynamical features of the mixtures are discussed in terms of diffusion constants and orientational correlation times as obtained from the time correlation functions for linear velocities and molecular dipole moments, respectively. Computed relaxation times indicate faster reorientational motion for dimethyl sulfoxide if acetonitrile is added. In contrast, the orientational dynamics of acetonitrile becomes stronger correlated upon dilution with dimethyl sulfoxide. The diffusion coefficients for both compounds follow this tendency.

Journal Article↗
PubMed2003

Development of hardware accelerator for molecular dynamics simulations: a computation board that calculates nonbonded interactions in cooperation with fast multipole method.

Evaluation of long-range Coulombic interactions still represents a bottleneck in the molecular dynamics (MD) simulations of biological macromolecules. Despite the advent of sophisticated fast algorithms, such as the fast multipole method (FMM), accurate simulations still demand a great amount of computation time due to the accuracy/speed trade-off inherently involved in these algorithms. Unless higher order multipole expansions, which are extremely expensive to evaluate, are employed, a large amount of the execution time is still spent in directly calculating particle-particle interactions within the nearby region of each particle. To reduce this execution time for pair interactions, we developed a computation unit (board), called MD-Engine II, that calculates nonbonded pairwise interactions using a specially designed hardware. Four custom arithmetic-processors and a processor for memory manipulation ("particle processor") are mounted on the computation board. The arithmetic processors are responsible for calculation of the pair interactions. The particle processor plays a central role in realizing efficient cooperation with the FMM. The results of a series of 50-ps MD simulations of a protein-water system (50,764 atoms) indicated that a more stringent setting of accuracy in FMM computation, compared with those previously reported, was required for accurate simulations over long time periods. Such a level of accuracy was efficiently achieved using the cooperative calculations of the FMM and MD-Engine II. On an Alpha 21264 PC, the FMM computation at a moderate but tolerable level of accuracy was accelerated by a factor of 16.0 using three boards. At a high level of accuracy, the cooperative calculation achieved a 22.7-fold acceleration over the corresponding conventional FMM calculation. In the cooperative calculations of the FMM and MD-Engine II, it was possible to achieve more accurate computation at a comparable execution time by incorporating larger nearby regions.

Algorithms↗
PubMed2006

Liposomal drug transport: a molecular perspective from molecular dynamics simulations in lipid bilayers.

Computational methods to predict drug permeability across biomembranes prior to synthesis are increasingly desirable to minimize the investment in drug design and development. Significant progress in molecular dynamics (MD) simulation methodologies applied to lipid bilayer membranes, for example, is making it possible to move beyond characterization of the membranes themselves to explore various thermodynamic and kinetic processes governing membrane binding and transport. Such methods are also likely to be directly applicable to the design and optimization of liposomal delivery systems. MD simulations are particularly valuable in addressing issues that are difficult to explore in laboratory experiments due to the heterogeneity of lipid bilayer membranes at the molecular level. Insights emerging from MD simulations are contributing to an understanding of which regions within bilayers are most and least favored by solutes at equilibrium as the solute structure is varied, local diffusivities of permeants, and the origin of the amplified selectivity to permeant size imposed by lipid bilayer membranes, particularly as changes in composition increase acyl chain ordering.

Biological Transport↗
PubMed2005

All-atom molecular dynamic simulations and relative NMR spectra study of weak C-H...O contacts in amide-water systems.

Amide-water mixtures are studied by all-atom molecular dynamics (MD) simulations and the relative temperature-dependent NMR experiment. The weak C-H...O contacts are found in the amide-water systems theoretically and experimentally. The statistical results of the average numbers of hydrogen bonds indicate that the methyl groups in amide molecules represent different capabilities in forming the weak C-H...O contacts. The statistics also imply that the C-H...O contacts are more obvious in the amide-rich region than those in the water-rich region. The temperature-dependent NMR spectra are also adopted to investigate the weak C-H...O contacts in the amide-water systems. The relative chemical shifts of the methyl groups are in good agreement with the MD simulations.

Amides↗
PubMed2003

Structural features of human initiation factor 4E, studied by X-ray crystal analyses and molecular dynamics simulations.

The structural features of human eIF4E were investigated by X-ray crystal analyses of its cap analog (m(7)GTP and m(7)GpppA) complexes and molecular dynamics (MD) simulations of cap-free and cap-bound eIF4Es, as well as the cap-bound Ser209-phosphorylated eIF4E. Crystal structure analyses at 2.0 A resolution revealed that the molecule forms a temple-bell-shaped surface of eight antiparallel beta-structures, three alpha-helices and ten loop structures, where the N-terminal region corresponds to the handle of the bell. This concave backbone provides a scaffold for the mRNA cap-recognition pocket consisting of three receiving parts for the 5'-terminal m(7)G base, the triphosphate, and the second nucleotide. The m(7)G base is sandwiched between the two aromatic side-chains of Trp102 and Trp56. The two (m(7)G)NH-O (Glu103 carboxy group) hydrogen bonds stabilize the stacking interaction. The basic residues of Arg157 and Lys162 and water molecules construct a binding pocket for the triphosphate moiety, where a universal hydrogen-bonding network is formed. The flexible C-terminal loop region unobserved in the m(7)GTP complex was clearly observed in the m(7)GpppA complex, as a result of the fixation of this loop by the interaction with the adenosine moiety, indicating the function of this loop as a receiving pocket for the second nucleotide. On the other hand, MD simulation in an aqueous solution system revealed that the cap-binding pocket, especially its C-terminal loop structure, is flexible in the cap-free eIF4E, and the entrance of the cap-binding pocket becomes narrow, although the depth is relatively unchanged. SDS-PAGE analyses showed that this structural instability is highly related to the fast degradation of cap-free eIF4E, compared with cap-bound or 4E-BP/cap-bound eIF4E, indicating the conferment of structural stability of eIF4E by the binary or ternary complex formation. MD simulation of m(7)GpppA-bound Ser209-phosphorylated eIF4E showed that the size of the cap-binding entrance is dependent on the ionization state in the Ser209 phosphorylation, which is associated with the regulatory function through the switching on/off of eIF4E phosphorylation.

Amino Acid Sequence↗
PubMed2005

Molecular modeling and molecular dynamics simulation of the human A2B adenosine receptor. The study of the possible binding modes of the A2B receptor antagonists.

A molecular model of the human A(2B) adenosine receptor containing seven transmembrane alpha helices connected by three intracellular and three extracellular hydrophilic loops had been constructed. A molecular docking of seven structurally diverse xanthine antagonists of the A(2B) receptor was performed, and the differences in their binding modes were investigated. The 1 ns molecular dynamics (MD) simulations of several obtained ligand-receptor complexes inserted into the phospholipid bilayer were carried out. The conformational changes of the A(2B) receptor occurring during MD simulations were explored, and the stable binding modes of the studied antagonists were determined. According to the models presented in this work, the involvement of the His251, Asn282, Ser92, Thr89, and some aromatic residues in ligand recognition was determined. The obtained binding modes of the A(2B) antagonists demonstrate good agreement with the site-directed mutagenesis data.

Adenosine A2 Receptor Antagonists↗
PubMed2005

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF.

Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.

Algorithms↗