We report over 30 μs of unrestrained molecular dynamics simulations of six protein-RNA complexes in explicit solvent. We utilize the AMBER ff99bsc0χ(OL3) RNA force field combined with the ff99SB protein force field and its more recent ff12SB version with reparametrized side-chain dihedrals. The simulations show variable behavior, ranging from systems that are essentially stable to systems with progressive deviations from the experimental structure, which we could not stabilize anywhere close to the starting experimental structure. For some systems, microsecond-scale simulations are necessary to achieve stabilization after initial sizable structural perturbations. The results show that simulations of protein-RNA complexes are challenging and every system should be treated individually. The simulations are affected by numerous factors, including properties of the starting structures (the initially high force field potential energy, resolution limits, conformational averaging, crystal packing, etc.), force field imbalances, and real flexibility of the studied systems. These factors, and thus the simulation behavior, differ from system to system. The structural stability of simulated systems does not correlate with the size of buried interaction surface or experimentally determined binding affinities but reflects the type of protein-RNA recognition. Protein-RNA interfaces involving shape-specific recognition of RNA are more stable than those relying on sequence-specific RNA recognition. The differences between the protein force fields are considerably smaller than the uncertainties caused by sampling and starting structures. The ff12SB improves description of the tyrosine side-chain group, which eliminates some problems associated with tyrosine dynamics.
- MeSH
- časové faktory MeSH
- proteiny chemie MeSH
- RNA chemie MeSH
- simulace molekulární dynamiky * MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- proteiny MeSH
- RNA MeSH
A method is proposed to measure global bending in DNA and RNA structures. It relies on a properly defined averaging of base-fixed coordinate frames, computes mean frames of suitably chosen groups of bases and uses these mean frames to evaluate bending. The method is applied to DNA A-tracts, known to induce considerable bend to the double helix. We performed atomistic molecular dynamics simulations of sequences containing the A(4)T(4) and T(4)A(4) tracts, in a single copy and in two copies phased with the helical repeat. Various temperature and salt conditions were investigated. Our simulations indicate bending by roughly 10 degrees per A(4)T(4) tract into the minor groove, and an essentially straight structure containing T(4)A(4), in agreement with electrophoretic mobility data. In contrast, we show that the published NMR structures of analogous sequences containing A(4)T(4) and T(4)A(4) tracts are significantly bent into the minor groove for both sequences, although bending is less pronounced for the T(4)A(4) containing sequence. The bending magnitudes obtained by frame averaging are confirmed by the analysis of superhelices composed of repeated tract monomers.
This review summarizes results concerning molecular interactions of nucleic acid bases as revealed by advanced ab initio quantum chemical (QM) calculations published in last few years. We first explain advantages and limitations of modern QM calculations of nucleobases and provide a brief history of this still rather new field. Then we provide an overview of key electronic properties of standard and selected modified nucleobases, such as their charge distributions, dipole moments, polarizabilities, proton affinities, tautomeric equilibria, and amino group hybridization. Then we continue with hydrogen bonding of nucleobases, by analyzing energetics of standard base pairs, mismatched base pairs, thio-base pairs, and others. After this, the nature of aromatic stacking interactions is explained. Also, nonclassical interactions in nucleic acids such as interstrand bifurcated hydrogen bonds, interstrand close amino group contacts, C [bond] H...O interbase contacts, sugar-base stacking, intrinsically nonplanar base pairs, out-of-plane hydrogen bonds, and amino-acceptor interactions are commented on. Finally, we overview recent calculations on interactions between nucleic acid bases and metal cations. These studies deal with effects of cation binding on the strength of base pairs, analysis of specific differences among cations, such as the difference between zinc and magnesium, the influence of metalation on protonation and tautomeric equlibria of bases, and cation-pi interactions involving nucleobases. In this review, we do not provide methodological details, as these can be found in our preceding reviews. The interrelation between advanced QM approaches and classical molecular dynamics simulations is briefly discussed.
- MeSH
- DNA chemie metabolismus MeSH
- elektrony * MeSH
- kationty metabolismus MeSH
- konformace nukleové kyseliny * MeSH
- molekulární modely MeSH
- RNA chemie metabolismus MeSH
- statická elektřina MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
- Research Support, U.S. Gov't, Non-P.H.S. MeSH
- Research Support, U.S. Gov't, P.H.S. MeSH
- Názvy látek
- DNA MeSH
- kationty MeSH
- RNA MeSH
Protonation of adenine carrying a Pt(II) moiety either at N7, N3, or N1 is possible in solution, but the site of protonation is influenced by the location of the Pt(II) electrophile and to some extent also by the overall charge of the metal entity (+2, +1, 0, -1), hence the other ligands (NH(3), Cl(-), OH(-)) bound to Pt(II). Quantum chemical calculations based on density functional theory (DFT) have been carried out for intrinsic protonation energies of adenine complexes carrying the following Pt(II) species at either of the three ring N atoms: [Pt(NH(3))(3)](2+) (1), trans- [Pt(NH(3))(2)Cl](+) (2a), cis-[Pt(NH(3))(2)Cl](+) (2b), trans-[Pt(NH(3))(2)Cl(2)] (3a), cis-[Pt(NH(3))Cl(2)] (3b), [PtCl(3)](-) (4), trans-[Pt(NH(3))(2)OH](+) (5a), cis-[Pt(NH(3))(2)(OH)](+) (5b), trans-[Pt(NH(3))(OH)(2)] (6a), cis-[Pt(NH(3))(OH)(2)] (6b), and [Pt(OH)(3)](-) (7). The data have been compared with results derived from solution studies (water) and X-ray crystallography, whenever available. The electrostatic effects associated with the charge of the metal entity have the major influence on the calculated intrinsic (gas phase) proton affinities, unlike the condensed phase data. Nevertheless, the relative gas phase trends correlate surprisingly well with condensed phase data; i.e., variation of the pK(a) values measured in solution is consistent with the calculated gas phase protonation energies. In addition to a systematic study of the ring proton affinities, proton transfer processes within the platinated adenine species were often observed when investigating Pt adducts with OH(-) ligands, and they are discussed in more detail. To the best of our knowledge, this is the first study attempting to find a systematic correlation between gas phase and condensed phase data on protonation of metalated nucleobases. The gas phase data provide a very useful complement to the condensed phase and X-ray experiments, showing that the gas phase studies are capable of valuable predictions and contribute to our understanding of the solvent and counterion effects on metal-assisted proton shift processes.
- MeSH
- adeninnukleotidy chemie MeSH
- chemické jevy MeSH
- chemie fyzikální MeSH
- difrakce rentgenového záření MeSH
- elektrolyty chemie MeSH
- kvantová teorie MeSH
- platina chemie MeSH
- protony MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, U.S. Gov't, Non-P.H.S. MeSH
- Research Support, U.S. Gov't, P.H.S. MeSH
- Názvy látek
- adeninnukleotidy MeSH
- elektrolyty MeSH
- platina MeSH
- protony MeSH
Large-scale molecular dynamics (MD) simulations have been utilized to study G-DNA quadruplex molecules containing mixed GCGC and all-guanine GGGG quartet layers. Incorporation of mixed GCGC quartets into G-DNA stems substantially enhances their sequence variability. The mixed quadruplexes form rigid assemblies that require integral monovalent cations for their stabilization. The interaction of cations with the all-guanine quartets is the leading contribution for the stability of the four-stranded assemblies, while the mixed quartets are rather tolerated within the structure. The simulations predict that two cations are preferred to stabilize a four-layer quadruplex stem composed of two GCGC and two all-guanine quartets. The distribution of cations in the structure is influenced by the position of the GCGC quartets within the quadruplex, the presence and arrangement of thymidine loops connecting the guanine/cytosine stretches forming the stems, and the cation type present (Na(+) or K(+)). The simulations identify multiple nanosecond-scale stable arrangements of the thymidine loops present in the molecules investigated. In these thymidine loops, several structured pockets are identified capable of temporarily coordinating cations. However, no stable association of cations to a loop has been observed. The simulations reveal several paths through the thymidine loop regions that can be followed by the cations when exchanging between the central ion channel in the quadruplex stem and the surrounding solvent. We have carried out 20 independent simulations while the length of simulations reaches a total of 90 ns, rendering this study one of the most extensive MD investigations carried out on nucleic acids so far. The trajectories provide a largely converged characterization of the structural dynamics of these four-stranded G-DNA molecules.
- MeSH
- chemické modely * MeSH
- chlorid draselný chemie MeSH
- chlorid sodný chemie MeSH
- cytosin chemie MeSH
- DNA chemie MeSH
- guanin chemie MeSH
- kationty jednomocné MeSH
- konformace nukleové kyseliny * MeSH
- molekulární modely MeSH
- oligonukleotidy chemie MeSH
- počítačová simulace MeSH
- roztoky MeSH
- sekvence nukleotidů MeSH
- termodynamika MeSH
- thymidin chemie MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- chlorid draselný MeSH
- chlorid sodný MeSH
- cytosin MeSH
- DNA MeSH
- guanin MeSH
- kationty jednomocné MeSH
- oligonukleotidy MeSH
- roztoky MeSH
- thymidin MeSH
The ability of the four-stranded guanine (G)-DNA motif to incorporate nonstandard guanine analogue bases 6-oxopurine (inosine, I), 6-thioguanine (tG), and 6-thiopurine (tI) has been investigated using large-scale molecular dynamics simulations. The simulations suggest that a G-DNA stem can incorporate inosines without any marked effect on its structure and dynamics. The all-inosine quadruplex stem d(IIII)(4) shows identical dynamical properties as d(GGGG)(4) on the nanosecond time scale, with both molecular assemblies being stabilized by monovalent cations residing in the channel of the stem. However, simulations carried out in the absence of these cations show dramatic differences in the behavior of d(GGGG)(4) and d(IIII)(4). Whereas vacant d(GGGG)(4) shows large fluctuations but does not disintegrate, vacant d(IIII)(4) is completely disrupted within the first nanosecond. This is a consequence of the lack of the H-bonds involving the N2 amino group that is not present in inosine. This indicates that formation of the inosine quadruplex could involve entirely different intermediate structures than formation of the guanosine quadruplex, and early association of cations in this process appears to be inevitable. In the simulations, the incorporation of 6-thioguanine and 6-thiopurine sharply destabilizes four-stranded G-DNA structures, in close agreement with experimental data. The main reason is the size of the thiogroup leading to considerable steric conflicts and expelling the cations out of the channel of the quadruplex stem. The G-DNA stem can accommodate a single thioguanine base with minor perturbations. Incorporation of a thioguanine quartet layer is associated with a large destabilization of the G-DNA stem whereas the all-thioguanine quadruplex immediately collapses.
- MeSH
- biofyzika MeSH
- biofyzikální jevy MeSH
- DNA chemie MeSH
- inosin chemie MeSH
- iontové kanály chemie MeSH
- konformace nukleové kyseliny * MeSH
- merkaptopurin chemie MeSH
- molekulární modely MeSH
- sodík chemie MeSH
- termodynamika MeSH
- thioguanin chemie MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- DNA MeSH
- inosin MeSH
- iontové kanály MeSH
- merkaptopurin MeSH
- sodík MeSH
- thioguanin 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.
Harmonic elastic constants of 3-11 bp duplex DNA fragments were evaluated using four 5 ns unrestrained molecular dynamics simulation trajectories of 17 bp duplexes with explicit inclusion of solvent and counterions. All simulations were carried out with the Cornell et al. force-field and particle mesh Ewald method for long-range electrostatic interactions. The elastic constants including anisotropic bending and all coupling terms were derived by analyzing the correlations of fluctuations of structural properties along the trajectories. The following sequences have been considered: homopolymer d(ApA)(n) and d(GpG)(n), and alternating d(GPC)(n) and d(APT)(n). The calculated values of elastic constants are in very good overall agreement with experimental values for random sequences. The atomic-resolution molecular dynamics approach, however, reveals a pronounced sequence-dependence of the stretching and torsional rigidity of DNA, while sequence-dependence of the bending rigidity is smaller for the sequences considered. The earlier predicted twist-bend coupling emerged as the most important cross-term for fragments shorter than one helical turn. The calculated hydrodynamic relaxation times suggest that damping of bending motions may play a role in molecular dynamics simulations of long DNA fragments. A comparison of elasticity calculations using global and local helicoidal analyses is reported. The calculations reveal the importance of the fragment length definition. The present work shows that large-scale molecular dynamics simulations represent a unique source of data to study various aspects of DNA elasticity including its sequence-dependence.
- MeSH
- algoritmy MeSH
- anizotropie MeSH
- DNA chemie genetika metabolismus MeSH
- konformace nukleové kyseliny * MeSH
- oligodeoxyribonukleotidy chemie genetika metabolismus MeSH
- párování bází MeSH
- počítačová simulace MeSH
- pohyb těles MeSH
- pružnost MeSH
- rozpouštědla MeSH
- sekvence nukleotidů MeSH
- statická elektřina MeSH
- termodynamika MeSH
- voda metabolismus MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- DNA MeSH
- oligodeoxyribonukleotidy MeSH
- rozpouštědla MeSH
- voda MeSH
Geometries, interaction energies and vibrational frequencies of base pairs, nucleoside pairs and nucleotide pairs were studied by ab initio Hartree-Fock (HF) method using MINI-1 basis set and empirical Cornell et al. force field (AMBER 4.1). A good agreement was found between HF/MINI-1 and AMBER results. In addition, both methods provide reasonable agreement with available high-level ab initio data. Finally, AMBER potential was used to determine the structure, energetics and vibrational frequencies of B-DNA pairs of trinucleotides. Stabilization energies of clusters are lowered when passing from base pairs to nucleoside pairs, nucleotide pairs and to pairs of trinucleotides. The lowest vibrations of base pairs and nucleoside pairs correspond to intermolecular motions of bases, specifically to buckle and propeller motions. In the case of pairs of larger subunits the lowest vibrations are of intramolecular nature (rotation around glycosidic bond, sugar and phosphate vibration). The spectra of these clusters became more complicated and quasi-degenerate. Intermolecular charge transfer between bases in H-bonded and stacked pairs is negligible, while a significant intramolecular charge transfer was observed.
Cation-pi interactions between cytosine and hexahydrated cations have been characterized using ab initio method with inclusion of electron correlation effects, assuming idealized and crystal geometries of the interacting species. Hydrated metal cations can interact with nucleobases in a cation-pi manner. The stabilization energy of such complexes would be large and comparable to the one for cation-pi complex with benzene. Further, polarized water molecules belonging to the hydration shell of the cation are capable to form a strong hydrogen bond interaction with the nitrogen lone electron pair of the amino groups of bases and enforce a pronounced sp3 pyramidalization of the nucleobase amino groups. However, in contrast to the benzene-cation complexes, the cation-pi configurations are highly unstable for a nucleobase since the conventional in plane binding of hydrated cations to the acceptor sites on the nucleobase is strongly preferred. Thus, a cation-pi interaction with a nucleobase can occur only if the position of the cation is locked above the nucleobase plane by another strong interaction. This indeed can occur in biopolymers and may have an effect on the local DNA architecture. Nevertheless, nucleobases have no intrinsic propensity to form cation-pi interactions.
- MeSH
- cytosin chemie metabolismus MeSH
- DNA chemie metabolismus MeSH
- kationty * MeSH
- kovy chemie MeSH
- molekulární modely MeSH
- voda metabolismus MeSH
- zinek MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, U.S. Gov't, Non-P.H.S. MeSH
- Research Support, U.S. Gov't, P.H.S. MeSH
- Názvy látek
- cytosin MeSH
- DNA MeSH
- kationty * MeSH
- kovy MeSH
- voda MeSH
- zinek MeSH
