Dissociation of the Watson-Crick base pairs in vacuum and in aqueous solution: a first-principles molecular dynamics study

The damage of the DNA structure can affect the correct functioning of the cellular processes. This work investigates the required forces to dissociate the Watson-Crick (WC) base pairs AT into A and T, and GC into G and C. The WC base pairs are immersed in water under realistic conditions of temperature, volume, and density that reproduce the main characteristics of a biological system. The simulations are based on first-principles molecular dynamics combined with steering atomic forces. In addition to the force intensities, the charge transfers between the nucleic acid bases, energy variations, and temperature fluctuations in the cleavage moments are reported. With the purpose of evaluating the effects of the aqueous medium, simulations of the WC base pairs in vacuum are included. The results considering the solvated medium are consistent with the experimental measurements, and show the importance of the aqueous solution to regulate the structural modifications of the nucleic acid bases. The investigation contributes with a novel molecular model in molecular simulations, and to better understand the biological processes where the DNA compounds play an active role in life forms.Communicated by Ramaswamy H. Sarma.

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Understanding the role of non-Watson-Crick base pairs in DNA-protein recognition: Structural and energetic aspects using crystallographic database analysis and quantum chemical calculation

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.

Determination of two-photon absorption in nucleobase analogues: a QR-DFT perspective

With the prevalence of fluorescence spectroscopy in biological systems, and the benefits of two-photon absorption techniques, presented here is an assessment of the two-photon accessibility of modern fluorescent nucleobase analogues utilising quadratic response DFT. Due to the complex environment experienced by these nucleobases, the two-photon spectra of each analogue has been assessed in the presence of both [Formula: see text]-stacked and hydrogen-bonding interactions involving the canonical nucleobases. Findings suggest that the [Formula: see text]-stacking environment provides a more significant effect on the spectra of the analogues studies than a hydrogen-bonding environment; analogue structures presenting high two-photon cross-section values for one or more states coincide with polycyclic extensions to preserved canonical base structure, as observed in the qA family of analogues, while analogue structures more closely resembling the structure of the base in question present a much more muted spectra in comparison. Results from this investigation have also allowed for the derivation of a number of design rules for the development of potential, two-photon specific, analogues for future use in both imaging and potential photochemical activation.

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BACKGROUND

Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.

METHODS

We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.

RESULTS

We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.

CONCLUSIONS

The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.

GENERAL SIGNIFICANCE

We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.