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 469 records · Page 26Linked to original sources

Cooperative fluctuations of unliganded and substrate-bound HIV-1 protease: a structure-based analysis on a variety of conformations from crystallography and molecular dynamics simulations.

The dynamics of HIV-1 protease, both in unliganded and substrate-bound forms have been analyzed by using an analytical method, Gaussian network model (GNM). The method is applied to different conformations accessible to the protein backbone in the native state, observed in crystal structures and snapshots from fully atomistic molecular dynamics (MD) simulation trajectories. The modes of motion obtained from GNM on different conformations of HIV-1 protease are conserved throughout the MD simulations. The flaps and 40's loop of the unliganded HIV-1 protease structure are identified as the most mobile regions. However, in the liganded structure these flaps lose mobility, and terminal regions of the monomers become more flexible. Analysis of the fast modes shows that residues important for stability are in the same regions of all the structures examined. Among these, Gly86 appears to be a key residue for stability. The contribution of residues in the active site region and flaps to the stability is more pronounced in the substrate-bound form than in the unliganded form. The convergence of modes in GNM to similar regions of HIV-1 protease, regardless of the conformation of the protein, supports the robustness of GNM as a potentially useful and predictive tool.

Binding Sites↗

Molecular dynamics study of CaCl2 in methanol.

The effective potential to describe interactions between the calcium ion and methanol molecule was derived from a potential energy surface. The energy surface was constructed from more than 10,000 points obtained from ab initio calculations, performed at the MP2/6-31G(d,p) level using Gaussian98. The derived effective potential and the Palinkas-Hawlicka-Heinzinger methanol model were employed in the molecular dynamics (MD) simulation of 0.25 M CaCl(2) solution in methanol at 298 K. Average distances between ions and the methanol sites--oxygen and carbon atoms--obtained from MD simulation agreed excellently with the X-ray data. The coordination number of the cation, greater as compared with the experimental results, is in accord with its concentration dependence.

Journal Article↗

Hydration dependence of the mass fractal dimension and anomalous diffusion of vibrational energy in proteins.

Vibrational dynamics of proteins and energy flow depend on protein geometry as well as interactions of a protein molecule with the surrounding solvent. We compute the mass fractal dimension D of proteins ranging from 100 to over 10,000 amino acids comparing values for the bare protein with those computed when buried and hydration waters are included in the calculation. Including water in the calculation increases D by about 0.3 to 2.87 on average above D computed for the dehydrated protein. The mass fractal dimension of proteins that are partially unfolded by molecular dynamics (MD) simulation is also computed and found to vary little when the radius of gyration changes within about 10% of that for the Protein Data Bank structure. MD simulations of vibrational energy diffusion in proteins reveal that the exponent characterizing anomalous diffusion of vibrational energy does not change much with hydration, which is seen to be due to an increase in the spectral dimension with hydration by a factor similar to the increase in D.

Computer Simulation↗

MM/PBSA analysis of molecular dynamics simulations of bovine beta-lactoglobulin: free energy gradients in conformational transitions?

The pH-driven opening and closure of beta-lactoglobulin EF loop, acting as a lid and closing the internal cavity of the protein, has been studied by molecular dynamics (MD) simulations and free energy calculations based on molecular mechanics/Poisson-Boltzmann (PB) solvent-accessible surface area (MM/PBSA) methodology. The forms above and below the transition pH differ presumably only in the protonation state of residue Glu89. MM/PBSA calculations are able to reproduce qualitatively the thermodynamics of the transition. The analysis of MD simulations using a combination of MM/PBSA methodology and the colony energy approach is able to highlight the driving forces implied in the transition. The analysis suggests that global rearrangements take place before the equilibrium local conformation is reached. This conclusion may bear general relevance to conformational transitions in all lipocalins and proteins in general.

Animals↗

Water motion in reverse micelles studied by quasielastic neutron scattering and molecular dynamics simulations.

Motion of water molecules in Aerosol OT [sodium bis(2-ethylhexyl) sulfosuccinate, AOT] reverse micelles with water content w(0) ranging from 1 to 5 has been explored both experimentally through quasielastic neutron scattering (QENS) and with molecular dynamics (MD) simulations. The experiments were performed at the energy resolution of 85 microeV over the momentum transfer (Q) range of 0.36-2.53 A(-1) on samples in which the nonpolar phase (isooctane) and the AOT alkyl chains were deuterated, thereby suppressing their contribution to the QENS signal. QENS results were analyzed via a jump-diffusion/isotropic rotation model, which fits the results reasonably well despite the fact that confinement effects are not explicitly taken into account. This analysis indicates that in reverse micelles with low-water content (w(0)=1 and 2.5) translational diffusion rate is too slow to be detected, while for w(0)=5 the diffusion coefficient is much smaller than for bulk water. Rotational diffusion coefficients obtained from this analysis increase with w(0) and are smaller than for bulk water, but rotational mobility is less drastically reduced than translational mobility. Using the Faeder/Ladanyi model [J. Phys. Chem. B 104, 1033 (2000)] of reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function F(S)(Q,t) for water hydrogens. Comparison of the time Fourier transform of this F(S)(Q,t) with the QENS dynamic structure factor S(Q,omega), shows good agreement between the model and experiment. Separate intermediate scattering functions F(S) (R)(Q,t) and F(S) (CM)(Q,t) were determined for rotational and translational motion. Consistent with the decoupling approximation used in the analysis of QENS data, the product of F(S) (R)(Q,t) and F(S) (CM)(Q,t) is a good approximation to the total F(S)(Q,t). We find that the decay of F(S) (CM)(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior is due to lower water mobility close to the interface and to confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay. However, rotational mobility of O-H bond vectors in the interfacial region remains fairly high due to the lower density of water-water hydrogen bonds in the vicinity of the interface.

Journal Article↗

Hierarchical multiscale simulation of electrokinetic transport in silica nanochannels at the point of zero charge.

Effects of nanoscale confinement and partial charges that stem from quantum calculations are investigated in silica slit channels filled with 1 M KCl at the point of zero charge by using a hierarchical multiscale simulation methodology. Partial charges of both bulk and surface atoms from ab initio quantum calculations that take into account bond polarization and electronegativity are used in molecular dynamics (MD) simulations to obtain ion and water concentration profiles for channel widths of 1.1, 2.1, 2.75, and 4.1 nm. The interfacial electron density profiles of simulations matched well with that of recent X-ray reflectivity experiments. By simulating corresponding channels with no partial charges, it was observed that the partial charges affect the concentration profiles and transport properties such as diffusion coefficients and mobilities up to a distance of about 3 sigma(O)(-)(O) from the surface. Both in uncharged and partially charged cases, oscillations in concentration profiles of K(+) and Cl(-) ions give rise to an electro-osmotic flow in the presence of an external electric field, indicating the presence of an electric double layer at net zero surface charge, contrary to the expectations from classical continuum theory. I-V curves in a channel-bath system using ionic mobilities from MD simulations were significantly different for channels with and without partial charges for channel widths less than 4.1 nm.

Chlorides↗

Molecular dynamics simulations of HPr under hydrostatic pressure.

The histidine-containing protein (HPr) plays an important role in the phosphotransferase system (PTS). The deformations induced on the protein structure at high hydrostatic pressure values (4, 50, 100, 150, and 200 MPa) were previously (H. Kalbitzer, A. Görler, H. Li, P. Dubovskii, A. Hengstenberg, C. Kowolik, H. Yamada, and K. Akasaka, Protein Science 2000, Vol. 9, pp. 693-703) analyzed by NMR experiments: the nonlinear variations of the amide chemical shifts at high pressure values were supposed to arise from induced shifts in the protein conformational equilibrium. Molecular dynamics (MD) simulations are here performed, to analyze the protein internal mobility at 0.1 MPa, and to relate the nonlinear variations of chemical shifts observed at high pressure, to variations in conformational equilibrium. The global features of the protein structure are only slightly modified along the pressure. Nevertheless, the values of the Voronoi residues volumes show that the residues of alpha-helices are more compressed that those belonging to the beta-sheet. The alpha-helices are also displaying the largest internal mobility and deformation in the simulations. The nonlinearity of the 1H chemical shifts, computed from the MD simulation snapshots, is in qualitative agreement with the nonlinearity of the experimentally observed chemical shifts.

Algorithms↗

Molecular dynamics simulations on the inhibition of cyclin-dependent kinases 2 and 5 in the presence of activators.

Interests in CDK2 and CDK5 have stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. Molecular dynamics (MD) simulations have not only become a viable approach to drug design because of advances in computer technology but are increasingly an integral part of drug discovery processes. It is common in MD simulations of inhibitor/CDK complexes to exclude the activator of the CDKs in the structural models to keep computational time tractable. In this paper, we present simulation results of CDK2 and CDK5 with roscovitine using models with and without their activators (cyclinA and p25). While p25 was found to induce slight changes in CDK5, the calculations support that cyclinA leads to significant conformational changes near the active site of CDK2. This suggests that detailed and structure-based inhibitor design targeted at these CDKs should employ activator-included models of the kinases. Comparisons between P/CDK2/cyclinA/roscovitine and CDK5/p25/roscovitine complexes reveal differences in the conformations of the glutamine around the active sites, which may be exploited to find highly selective inhibitors with respect to CDK2 and CDK5.

Cyclin-Dependent Kinase 2↗

Mechanisms of ammonia activation and ammonium ion inhibition of quinoprotein methanol dehydrogenase: a computational approach.

The mechanism of methanol oxidation by quinoprotein methanol dehydrogenase (MDH.PQQ) in combination with methanol (MDH.PQQ.methanol) involves Glu-171--CO2(-) general base removal of the hydroxyl proton of methanol in concert with hydride equivalent transfer to the >C5=O quinone carbon of pyrroloquinoline quinone (PQQ) and rearrangement to hydroquinone (PQQH2) with release of formaldehyde. Molecular dynamics (MD) studies of the structures of MDH.PQQ.methanol in the presence of activator NH3 and inhibitor NH4(+) have been carried out. In the MD structure of MDH.PQQ.methanol.NH3, the hydrated NH3 resides at a distance of approximately 24 A away from methanol and the ortho-quinone portion of PQQ. As such, influence of NH3 on the oxidation reaction is not probable. We find that NH4(+) competes with the substrate by hydrogen-bonding to Glu-171CO2(-) such that the MDH.PQQ.methanol.NH4(+) complex is not reactive. Ammonia readily forms imines with quinone. Imines are present in solution as neutral (>C5=NH) and protonated (>C5=NH2(+)) species. MD simulations establish that the >C5=NH2(+) derivative of MDH.PQQ(NH2(+).methanol structure is unreactive because of the nonproductive means of methanol binding. The structure obtained by the MD simulations with the neutral >C5=NH imine of MDH.PQQ(NH).methanol structure is similar to the reactive MDH.PQQ.methanol complex. This active site geometry allows for catalysis of hydride equivalent transfer to the >C5=NH of PQQ(NH) by concerted Glu-171CO(2)(-) general-base removal of the H-OCH3 proton and Arg-324H+ general-acid proton transfer to the imine nitrogen. Enzyme-bound <C5(H)NH2 derivative of PQQ [PQQ(NH)] and CH(2)O product are formed.

Alcohol Oxidoreductases↗

Refinement of protein cores and protein-peptide interfaces using a potential scaling approach.

Refinement of side chain conformations in protein model structures and at the interface of predicted protein-protein or protein-peptide complexes is an important step during protein structural modelling and docking. A common approach for side chain prediction is to assume a rigid protein main chain for both docking partners and search for an optimal set of side chain rotamers to optimize the steric fit. However, depending on the target-template similarity in the case of comparative protein modelling and on the accuracy of an initially docked complex, the main chain template structure is only an approximation of a realistic target main chain. An inaccurate rigid main chain conformation can in turn interfere with the prediction of side chain conformations. In the present study, a potential scaling approach (PS-MD) during a molecular dynamics (MD) simulation that also allows the inclusion of explicit solvent has been used to predict side chain conformations on semi-flexible protein main chains. The PS-MD method converges much faster to realistic protein-peptide interface structures or protein core structures than standard MD simulations. Depending on the accuracy of the protein main chain, it also gives significantly better results compared with the standard rotamer search method.

Amino Acid Sequence↗

Dynamics of liquid benzene: a cage analysis.

Dynamics of single molecules in liquids, inspected in the picosecond time scale by means of spectroscopic measurements or molecular-dynamics (MD) simulations, reveals a complex behavior which can be addressed as due to local confinement (cage). This work is devoted to the analysis of cage structures in liquid benzene, obtained from MD simulations. According to a paradigm proposed for previous analysis of atomic and molecular liquids [see, for example, A. Polimeno, G. J. Moro, and J. H. Freed, J. Chem. Phys. 102, 8094 (1995)], the istantaneous cage structure is specified by the frame of axes which identifies the molecular configuration at the closest minimum on the potential-energy landscape. In addition, the modeling of the interaction potential between probe molecule and molecular environment, based on symmetry considerations, and its parametrization from the MD trajectories, allows the estimation of the structural parameters which quantify the strength of molecular confinement. Roto-translational dynamics of probe and related cage with respect to a laboratory frame, dynamics of the probe within the cage (vibrations, librations, re-orientational motions), and the restructuring processes of the cage itself are analyzed in terms of selected time self-correlation functions. A time-scale separation between the processes is established. Moreover, by exploiting the evidence of fast vibrational motions of the probe with respect to the cage center, an orientational effective potential is derived to describe the caging in the time scale longer than approximately 0.2 ps.

Journal Article↗

A structural and mechanistic study of the oxidation of methionine residues in hPTH(1-34) via experiments and simulations.

The relationship between the conformational properties of 1-34 human parathyroid hormone [hPTH(1-34)] and the oxidation of its methionine residues, Met8 and Met18, by hydrogen peroxide is analyzed as a function of pH by measuring the rates of oxidation and by performing MD simulations with an explicit representation of water molecules. Between pH 4 and pH 8, both Met8 and Met18 have nearly pH independent rates of oxidation, and Met18 is oxidized at a rate that is 90-100% of that of freeMet and 10-20% faster than that of Met8. We also found that average 2SWCNs calculated from MD simulations correlate well to the rates of oxidation of Met8 and Met18. The use of 2SWCNs is based on the mechanism that we proposed, the water-mediated mechanism, in which water molecules stabilize the transition state via specific interactions, but the transfer of protons (acid-catalyzed mechanism) does not play a role [Chu, J. W., and Trout, B. L. (2004) J. Am. Chem. Soc. 126 (3), 900-908]. Only at very low pH values, pH 1 for the oxidation of freeMet, does the acid-catalyzed oxidation mechanism become important. For the oxidation of Met8 and Met18 in hPTH(1-34), the acid-catalyzed mechanism becomes significant at a higher pH value, pH 2, probably due to the proximity of nearby acidic residues to Met8 (Glu4) and Met18 (Glu22). In this study, we have demonstrated that the chemistry of oxidation and the structure of polypeptides can be correlated via a detailed understanding of the reaction mechanism, appropriate sampling of configurational space, and a suitable choice of a structural property, water coordination number.

Circular Dichroism↗

Investigation of the solution structure of chymotrypsin inhibitor 2 using molecular dynamics: comparison to x-ray crystallographic and NMR data.

The native solution structure and dynamics of chymotrypsin inhibitor 2 (CI2) have been studied using a long (5.3 ns) molecular dynamics (MD) simulation without any imposed restraints. The majority of the experimentally observed spin-spin coupling constants, short- and long-range nuclear Overhauser effect (NOE) cross peaks and the amide hydrogen exchange behavior were reproduced by the MD simulation. This good correspondence suggests that the major structural features of the protein during the simulation are representative of the true protein structure in solution. Two water molecules formed hydrogen bond bridges between beta2 and beta3, in agreement with X-ray crystallographic data and a recent reassessment of the solution structure using time-averaged NMR restraints during MD refinement. The active-site loop of the protein displayed the greatest structural changes and the highest mobility. When this loop region was excluded, the average Calpha r.m.s. deviation of the simulated solution structures from the crystal structure was approximately 1.5 Angstrom from 0.5 to 5.3 ns. There is structural heterogeneity in particular regions of the NMR-derived solution structures, which could be a result of imprecision or true internal motion. A study of the distribution of mobility through the protein allows us to distinguish between these two alternatives. In particular, deviations in the active-site loop appear to be a result of heightened mobility, which is also supported by good correspondence between calculated and experimental S2 N-H order parameters. On the other hand, other ill-defined regions of the NMR-derived structures are well defined in the simulation and are probably the result of a lack of structural restraints (i.e. NOEs), as opposed to reflecting the true mobility.

Computer Simulation↗

The mechanism of cis-trans isomerization of prolyl peptides by cyclophilin.

The mechanism of cis-trans isomerization of prolyl peptides catalyzed by cyclophilin (CyP) was studied computationally via molecular dynamics (MD) simulations of the transition state (TS) and the cis and trans forms of the ground state (GS), when bound to CyP and when free in aqueous solution. The MD simulations include four enzyme-bound species of tetrapeptide (Suc-Ala-XC([double bond]O)-NPro-Phe-pNA; X = Gly, Trp, Ala, and Leu). In water, the prolyl amide bond is favorably planar with the presence of conformers exhibiting +/-20 degrees twist of the C-N dihedral. In the active site a hydrogen bond between the cis-prolyl amide carbonyl O and the backbone amide N-H of Asn102 retains the 20 degrees twist of the C-N dihedral. The TS structure is characterized by a 90 degrees twist of the amide C-N bond and a more favorable interaction with Asn102 due to the shorter distance between Asn102(HN) and the amide carbonyl O. The conformational change of cis --> TS also involves pyramidalization of the amide N, which results in the formation of a hydrogen bond between the amide N and the guanidino group of Arg55. Both Asn102 and Arg55 are held in the same position in CyP.cis-isomer as in CyP.TS. In the ligand-free CyP the Arg55 guanidino group is highly disorganized and Asn102 is displaced 1 A from the position in the ligand-bound CyP. Thus, the organization of Arg55 and Asn102 occurs upon substrate binding. The geometrical complimentarity of the organized enzyme structure to the TS structure is a result of preferential binding of the proline N and the amide carbonyl of the TS compared to that of GS. However, the N-terminal part (Suc-Ala) becomes repositioned in the TS such that two hydrogen bonds disappear, one hydrogen bond appears and two other hydrogen bonds becomes weaker on the conversion of CyP.cis to CyP.TS. During this conversion, total hydrophobic contact between enzyme and the peptide is preserved. Thus, the interaction energies of GS and TS with enzyme are, as a whole, much alike. This does not support the contention that TS is bound more tightly than GS by K(m)/K(TS) = 10(6) in the cis --> trans reaction. Repositioning of the N-terminal part of the peptide on CyP.TS formation becomes more pronounced when the substrate X residue is changed from Gly < Trp < Ala < Leu. We propose that the larger turning of the N-terminus is responsible for the larger value of the experimentally observed Delta S(++) and Delta H(++), which sum up to little change in Delta G(++). The positioning of the Arg55 and the degree of 20 degrees twist of the amide C-N bond are considered as criteria for Near Attack Conformers (NACs) in cis-trans isomerization. NACs account for approximately 30% of the total GS populations of the cis-isomer. Similar NAC populations were observed with four different substrates. This is consistent with the insensitivity of enzymatic activity to the nature of the X residue. Also, the NAC population in CyP.trans-AAPF was comparable to that in CyP.cis-AAPF, in accord with similar experimentally measured rates of the cis --> trans and trans --> cis reaction in CyP. These NACs, found in CyP.cis and CyP.trans, resemble only one of the four possible TS configurations in the water reaction. The identity of this TS structure (syn/exo) is in accord with experimentally determined KIE values in the enzymatic reaction. However, the geometry of the active site was also complementary to another TS structure (anti/exo) that was not detected in the active site by the same KIE measurements, implying that the geometrical fitness of the TS cannot be a single determining factor for enzymatic reactions.

Binding Sites↗

Structures and properties of Newton black films characterized using molecular dynamics simulations.

We used molecular dynamics (MD) simulations to investigate the structures and properties of Newton black films (NBF) for several surfactants: sodium dodecyl sulfate (SDS), cetyltrimethylammonium bromide (C16TAB), and surfactin using film thicknesses up to 10 nm. By calculating the interface formation energy for various packing conditions on the surface pressure-area isotherm, we found that the most probable surface concentration is approximately 42 A(2)/molecule for SDS and C16TAB and approximately 170 A(2)/molecule for surfactin. We then used this most probable concentration of each surfactant to simulate NBF with various film thicknesses. From analyzing the disjoining pressure-film thickness isotherms with the density profiles and the solvation coordination number, we found that the increase of the disjoining pressure during the film thinning was coupled with the change in inner structure of the NBF (i.e., density profile and the solvation of ionic entities). In the range of film thicknesses less than approximately 30 A, the disjoining pressures for the SDS and C16TAB were found to be larger than that of the surfactin. We predicted the Gibbs elasticity (175 dyn/cm for surfactin; 109 dyn/cm for C16TAB; 38 dyn/cm for SDS) required to assess the stability of NBF against surface concentration fluctuations, and the shear modulus (6.5 GPa for the surfactin; 6.1 GPa for the C16TAB; 3.5 GPa for the SDS) and the yield stress (approximately 0.8 GPa for surfactin; approximately 0.8 GPa for C16TAB; approximately 0.4 GPa for the SDS) to assess the mechanical stability against the externally imposed mechanical perturbation.

Algorithms↗

Prediction of binding affinities between the human amphiphysin-1 SH3 domain and its peptide ligands using homology modeling, molecular dynamics and molecular field analysis.

The SH3 domain of the human protein amphiphysin-1, which plays important roles in clathrin-mediated endocytosis, actin function and signaling transduction, can recognize peptide motif PXRPXR (X is any amino acid) with high affinity and specificity. We have constructed a complex structure of the amphiphysin-1 SH3 domain and a high-affinity peptide ligand PLPRRPPRA using homology modeling and molecular docking, which was optimized by molecular dynamics (MD). Three-dimensional quantitative structure-affinity relationship (3D-QSAR) analyses on the 200 peptides with known binding affinities to the amphiphysin-1 SH3 domain was then performed using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMSIA model showed promising predictive power, giving good predictions for about 95% of the peptides in the test set (absolute prediction errors less than 1.0). It was used to validate peptide-SH3 binding structure and provide insight into the structural requirements for binding of peptides to SH3 domains. Finally, MD simulations were performed to analyze the interaction between the SH3 domain and another peptide GFPRRPPPRG that contains with the PXRPXsR (s represents residues with small side chains) motif. MD simulations demonstrated that the binding conformation of GFPRRPPPRG is quite different from that of PLPRRPPRAA especially the four residues at the C terminal, which may explain why the CoMSIA model cannot give good predictions on the peptides of the PXRPXsR motif. Because of its efficiency and predictive power, the 3D-QSAR model can be used as a scoring filter for predicting peptide sequences bound to SH3 domains.

Amino Acid Sequence↗

Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations.

Computational studies of proteins based on empirical force fields represent a powerful tool to obtain structure-function relationships at an atomic level, and are central in current efforts to solve the protein folding problem. The results from studies applying these tools are, however, dependent on the quality of the force fields used. In particular, accurate treatment of the peptide backbone is crucial to achieve representative conformational distributions in simulation studies. To improve the treatment of the peptide backbone, quantum mechanical (QM) and molecular mechanical (MM) calculations were undertaken on the alanine, glycine, and proline dipeptides, and the results from these calculations were combined with molecular dynamics (MD) simulations of proteins in crystal and aqueous environments. QM potential energy maps of the alanine and glycine dipeptides at the LMP2/cc-pVxZ//MP2/6-31G* levels, where x = D, T, and Q, were determined, and are compared to available QM studies on these molecules. The LMP2/cc-pVQZ//MP2/6-31G* energy surfaces for all three dipeptides were then used to improve the MM treatment of the dipeptides. These improvements included additional parameter optimization via Monte Carlo simulated annealing and extension of the potential energy function to contain peptide backbone phi, psi dihedral crossterms or a phi, psi grid-based energy correction term. Simultaneously, MD simulations of up to seven proteins in their crystalline environments were used to validate the force field enhancements. Comparison with QM and crystallographic data showed that an additional optimization of the phi, psi dihedral parameters along with the grid-based energy correction were required to yield significant improvements over the CHARMM22 force field. However, systematic deviations in the treatment of phi and psi in the helical and sheet regions were evident. Accordingly, empirical adjustments were made to the grid-based energy correction for alanine and glycine to account for these systematic differences. These adjustments lead to greater deviations from QM data for the two dipeptides but also yielded improved agreement with experimental crystallographic data. These improvements enhance the quality of the CHARMM force field in treating proteins. This extension of the potential energy function is anticipated to facilitate improved treatment of biological macromolecules via MM approaches in general.

Algorithms↗

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study.

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.

Biophysical Phenomena↗