Ab initio molecular dynamics
Dotaz
Zobrazit nápovědu
Ab initio and molecular simulation methods were used in calculations of the neutral individual betulin molecule, and molecular simulations were used to optimize the betulin molecule immersed in various amounts of water. Individual betulin was optimized in different force fields to find the one exhibiting best agreement with ab initio calculations obtained in the Gaussian03 program. Dihedral torsions of active groups of betulin were determined for both procedures, and related calculated structures were compared successfully. The selected force field was used for subsequent optimization of betulin in a water environment, and a conformational search was performed using quench molecular dynamics. The total energies of betulin and its interactions in water bulk were calculated, and the influence of water on betulin structure was investigated.
The dynamic structure and potential energy surface of adenine...thymine and guanine...cytosine base pairs and their methylated analogues interacting with a small number (from 1 to 16 molecules) of organic solvents (methanol, dimethylsulfoxide, and chloroform) were investigated by various theoretical approaches starting from simple empirical methods employing the Cornell et al. force field to highly accurate ab initio quantum chemical calculations (MP2 and particularly CCSD(T) methods). After the simple molecular dynamics simulation, the molecular dynamics in combination with quenching technique was also used. The molecular dynamics simulations presented here have confirmed previous experimental and theoretical results from the bulk solvents showing that, whereas in chloroform the base pairs create hydrogen-bonded structures, in methanol, stacked structures are preferred. While methanol (like water) can stabilize the stacked structures of the base pairs by a higher number of hydrogen bonds than is possible in hydrogen-bonded pairs, the chloroform molecule lacks such a property, and the hydrogen-bonded structures are preferred in this solvent. The large volume of the dimethylsulfoxide molecule is an obstacle for the creation of very stable hydrogen-bonded and stacked systems, and a preference for T-shaped structures, especially for complexes of methylated adenine...thymine base pairs, was observed. These results provide clear evidence that the preference of either the stacked or the hydrogen-bonded structures of the base pairs in the solvent is not determined only by bulk properties or the solvent polarity but rather by specific interactions of the base pair with a small number of the solvent molecules. These conclusions obtained at the empirical level were verified also by high-level ab initio correlated calculations.
- MeSH
- chloroform chemie MeSH
- dimethylsulfoxid chemie MeSH
- financování organizované MeSH
- konformace nukleové kyseliny MeSH
- methanol chemie MeSH
- metylace DNA MeSH
- molekulární modely MeSH
- nukleotidy chemie MeSH
- párování bází MeSH
- počítačová simulace MeSH
- voda chemie MeSH
- vodíková vazba MeSH
- zastoupení bazí MeSH
Abstract Aromatic stacking of nucleic acid bases is one of the key players in determining the structure and dynamics of nucleic acids. The arrangement of nucleic acid bases with extensive overlap of their aromatic rings gave rise to numerous often contradictory suggestions about the physical origins of stacking and the possible role of delocalized electrons in stacked aromatic π systems, leading to some confusion about the issue. The recent advance of computer hardware and software finally allowed the application of state of the art quantum-mechanical approaches with inclusion of electron correlation effects to study aromatic base stacking, now providing an ultimitate qualitative description of the phenomenon. Base stacking is determined by an interplay of the three most commonly encountered molecular interactions: dispersion attraction, electrostatic interaction, and short-range repulsion. Unusual (aromatic- stacking specific) energy contributions were in fact not evidenced and are not necessary to describe stacking. The currently used simple empirical potential form, relying on atom-centered constant point charges and Lennard-Jones van der Waals terms, is entirely able to reproduce the essential features of base stacking. Thus, we can conclude that base stacking is in principle one of the best described interactions in current molecular modeling and it allows to study base stacking in DNA using large-scale classical molecular dynamics simulations. Neglect of cooperativity of stacking appears to be the most serious approximation of the currently used force field form. This review summarizes recent developments in the field. It is written for an audience that is not necessarily expert in computational quantum chemistry and follows up on our previous contribution (Sponer et. al., J. Biomol. Struct. Dyn. 14, 117, (1997)). First, the applied methodology, its accuracy, and the physical nature of base stacking is briefly overviewed, including a comment on the accuracy of other molecular orbital methods and force fields. Then, base stacking is contrasted with hydrogen bonding, the other dominant force in nucleic acid structure. The sequence dependence and cooperativity of base stacking is commented on, and finally a brief introduction into recent progress in large-scale molecular dynamics simulations of nucleic acids is provided. Using four stranded DNA assemblies as an example, we demonstrate the efficacy of current molecular dynamics techniques that utilize refined and verified force fields in the study of stacking in nucleic acid molecules.
The cis-platin binding to the d(CCTGGTCC)*d(GGACCAGG) model DNA octamer was monitored with infrared absorption (IR) and vibrational circular dichroism (VCD) spectroscopies. The spectra were modeled with the aid of density functional computations and a Cartesian coordinate-based transfer of molecular property tensors from smaller DNA fragments. Because of the fragmentation, the tensors could be calculated with a higher precision. Environmental effects, such as the presence of the solvent or the cis-platin ligand, could be included in the modeling. The solvent was modeled by an explicit inclusion of hydrogen-bound water molecules, positions of which were estimated from a molecular dynamics simulation, or by the polarized continuum COSMO model. The B3LYP and BPW91 functionals used for the calculations of the spectral parameters were combined with the relativistic LANL2DZ platinum pseudo-potentials. The simulations reproduced the main IR and VCD DNA spectral features and explained most of the changes observed experimentally upon metal binding. The results confirmed that the influence of the ligand on DNA vibrational properties is quite complex; it originates in the geometry deformation and normal mode coupling pattern changes of the platinated octamer, as well as in local perturbations of the electronic structure and force field of the GC base pairs to which the platinum is bound. Many of the local effects could be accounted for by a point charge used in place of the metal in the GC complex.
- MeSH
- chemické modely MeSH
- cisplatina chemie MeSH
- DNA chemie MeSH
- farmaceutická chemie metody MeSH
- financování organizované MeSH
- fyzikální chemie metody MeSH
- guanidin chemie MeSH
- ionty MeSH
- ligandy MeSH
- magnetická rezonanční spektroskopie MeSH
- molekulární konformace MeSH
- molekulární struktura MeSH
- oligonukleotidy chemie MeSH
- protinádorové látky chemie MeSH
- spektrofotometrie infračervená metody MeSH
- spektrofotometrie metody MeSH
A complete scan of the potential-energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Moller-Plesset second-order perturbation theory (RI-MP2). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T-shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol(-1). The different behaviors of nonmethylated and methylated trimers is also discussed. The high-level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple-zeta basis set predicts 11 uracil trimer interaction energies with a root-mean-square error of 2.3 kcal mol(-1) relative to highly correlated ab initio theoretical calculations.
Altogether 13 keto and enol tautomers of uracil and 13 keto and enol tautomers of thymine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in a water environment. Bulk water was described using the thermodynamic integration method, Conductor-like polarizable continuum model (C-PCM, COSMO) and hybrid model (C-PCM + 1-2 explicit water molecules). The structures of various tautomers were determined at the RI-MP2 level using the TZVPP basis set while relative energies were determined at the CCSD(T) level. The relative free energies at 298 K were based on the relative energies mentioned above and zero-point vibration energies, and temperature dependent enthalpy terms and entropies evaluated at the MP2/6-31G** level. The effect of bulk solvent on the relative stability of uracil and thymine tautomers was studied using molecular dynamics free energy calculations by means of the thermodynamic integration method and self-consistent reaction field. Despite the completely different nature of these methods they provide comparable solvation free energies. Besides theoretical investigation, experimental detection of uracil and thymine tautomers was performed by means of steady-state fluorescence. We conclude that it is impossible to utilize the method used by Suwaiyan and Morsy (M. A. Morsy, A. M. Al-Somali and A. Suwaiyan, J. Phys. Chem. B, 1999, 103(50), 11205) for tautomer detection, even if a very sensitive fluorimeter is used. Theoretical relative energies and free energies for isolated uracil and thymine tautomers support the existence of the canonical form only. The microhydrated environment and bulk solvent stabilize enol forms more than the canonical keto one, but gas phase destabilization of these enol forms is too high. Population of rare enol forms of uracil and thymine in bulk water will thus be very low and canonical structure will also be dominant in this phase.
- MeSH
- biofyzika metody MeSH
- chemické modely MeSH
- fluorescenční spektrometrie metody MeSH
- fluorometrie metody MeSH
- konformace nukleové kyseliny MeSH
- nukleové kyseliny chemie MeSH
- plyny MeSH
- rozpouštědla chemie MeSH
- termodynamika MeSH
- thymin chemie MeSH
- uracil chemie MeSH
- voda chemie MeSH
- Publikační typ
- práce podpořená grantem MeSH
