Modeling nucleic acids

Nucleic acids are an important class of biological macromolecules that carry out a variety of cellular roles. For many functions, naturally occurring DNA and RNA molecules need to fold into precise three-dimensional structures. Due to their self-assembling characteristics, nucleic acids have also been widely studied in the field of nanotechnology, and a diverse range of intricate three-dimensional nanostructures have been designed and synthesized. Different physical terms such as base-pairing and stacking interactions, tertiary contacts, electrostatic interactions and entropy all affect nucleic acid folding and structure. Here we review general computational approaches developed to model nucleic acid systems. We focus on four key areas of nucleic acid modeling: molecular representation, potential energy function, degrees of freedom and sampling algorithm. Appropriate choices in each of these key areas in nucleic acid modeling can effectively combine to aid interpretation of experimental data and facilitate prediction of nucleic acid structure.     

Theoretical calculations of base-base interactions in nucleic acids: I. Stacking interactions in free bases.

Stacking interactions in free bases were computed on the basis of molecular association. The results of the calculations were compared with the stacking patterns observed in a few single crystals of nucleic acid components as examples. The following are the conclusions: (i) there can be two types of stacking pattern classified as normal and inverted types for any two interacting bases and both can be energetically favourable (ii) in both the types the stacking interaction is a combined effect of the overlap of the interacting bases and relative positions and orientations of the atomic centres of the two bases (iii) crystal symmetry and H-bonding interaction may influence stacking patterns.     

A free energy analysis of nucleic acid base stacking in aqueous solution.

This paper reports a theoretical study of the free energy contributions to nucleic acid base stacking in aqueous solution. Electrostatic interactions are treated by using the finite difference Poisson-Boltzmann method and nonpolar effects are treated with explicit calculation of van der Waals interactions and/or free energy-surface area relationships. Although for some pairs of bases there is a favorable Coulombic interaction in the stacked conformation, generally the net effect of electrostatic interactions is to oppose stacking. This result is caused by the loss of favorable base-solvent electrostatic interactions, that accompany the partial removal of polar atoms from water in the stacked conformation. Nonpolar interactions, involving the hydrophobic effect and enhancement of van der Waals interactions caused by close-packing, drive stacking. The calculations qualitatively reproduce the experimental dependence of stacking free energy on purine-pyrimidine composition.     

Stacking and hydrogen bond interactions between adenine and gallic acid

We have performed DFT and DFT-SAPT calculations on dimers of gallic acid, the model system for plant polyphenols, and the DNA base adenine. These dimers were selected for this study as they exhibit simultaneously hydrogen bonds and stacking interactions and it allows to quantify the relative values of these interactions. We calculate the relationships between the stability of the complexes, the charge transfer between monomers and the properties of the intermolecular bonds including hydrogen bonds and other bonds that do not involve hydrogen atoms. DFT-SAPT calculations were also performed to obtain the participation of the different types of energy and so the resulting physical effects. The results show that the presence of hydrogen bonds is the main stabilizing factor for dimers: the higher number and strength, the lower the dimer energy. The contribution of stacking to the stabilization is related to the strength and number of bonds between non-hydrogen atoms and quantified by the contribution of the dispersion terms to the interaction energy. Dimers I and II are mainly stabilized due to hydrogen bonds whereas dimer III is mainly stabilized by stacking interactions.     

Specific recognition of single-stranded nucleic acids. Interaction of tryptophan-containing peptides with native, denatured, and ultraviolet-irradiated DNA

The binding of the tripeptide Lys-Trp-Lys to native, denatured, and ultraviolet-irradiated DNAs has been investigated by fluorescence spectroscopy. Two types of complexes are formed which both involve electrostatic interactions. Only one of them involves a stacking of the tryptophyl ring with nucleic acid bases. Quantitative analysis of fluorescence data shows that this stacking interaction is strongly favored in denatured as compared to native DNA. In ultraviolet-irradiated DNA, the peptide Lys-Trp-Lys binds selectively to unpaired regions around thymine dimers. Due to the stacking interaction of the aromatic amino acid with nucleic acid bases, this simple tripeptide is therefore able to discriminate between single-stranded and double-stranded regions in a nucleic acid.     

Solvent interactions in nucleic acid crystal hydrates

A detailed knowledge of structural and energetic aspects of water-nucleic acid interactions is essential for understanding the role of solvent in stabilizing the various helical forms of nucleic acids. In this study, computer simulation techniques have been used to predict structural properties of solvent networks in small nucleic acid crystal hydrates. A detailed comparison of predicted and experimental results on the structure of the solvent networks is presented and includes an analysis of both the local environment and hydrogen bond pattern of each water molecule. A correlation between the environment of each unique water molecule and its energetic properties (such a dipole moment and binding energy) is seen. As in the previous studies on small amino acid hydrate crystals, non-pair additive (cooperative) effects are found to be non-negligible. It is concluded that the potential functions used in this initial study lead to simulated solvent networks in reasonable agreement with experimental data. Thus, it is now feasible to use them in studies of hydration of larger helical fragments of nucleic acids of more direct biological interest.     

Effects of 5-[S-(2,4-dinitrophenyl)-thio]-2'-deoxyuridine analog incorporation on the structure and stability of DNA hybrids: implications for the design of nucleic acid probes

Labeled nucleic acid probes are used as diagnostic tools by detecting changes in gene expression upon hybridization to target RNAs or DNAs that are related to specific disease genes. 5-[S-(2, 4-Dinitrophenyl)-thio]-2'-deoxyuridine analog represents an excellent nucleic acid label, containing the DNP group which functions both as a probe and as a precursor for the introduction of a variety of fluorescent groups. This study describes thermal denaturation hybridization experiments with oligonucleotides containing the 5-[S-(2,4-dinitrophenyl)-thio]-2'-deoxyuridine analog. Using molecular modeling techniques, the effects of this analog on the hybrid structure and stability were examined, including (i) analog conformation, (ii) hydrogen bonding, (iii) stacking interactions and (iv) hybrid helical geometry. This analog does not prohibitively affect the hybrid thermal stability and incorporation of the analog does not compromise the structural integrity of the double helix. In particular, the sequence-dependence of the analog effects and the dependence on the modification site relative to the end(s) of the helix were investigated. Findings described here should provide guidelines in the rational design of nucleic acid probes.     

Stereochemistry of nucleic acids and their constituents. X. Solid-state base-stacking patterns in nucleic acid constituents and polynucleotides

The base-stacking patterns in over 70 published crystal structures of nucleic acid constituents and polynucleotides were examined. Several recurring stacking patterns were found. Base stacking in the solid state apparently is very specific, with particular modes of interaction persisting in various crystalline environments. The vertical stacking of purities and pyrimidines in polynucleotides is similar to that observed in crystals of nucleic acid constituents. Only partial base overlap was found in the majority of the structures examined. Usually, the base overlap is accomplished by positioning polar substituents over the ring system of an adjacent base. The stacking interactions are similar to those found in the crystal structures of other polar aromatic compounds, but are considerably different from the ring–ring interactions in nonpolar aromatic compounds. Apparently, dipole-induced dipole forces are largely responsible for solid-state base stacking. It is found that halogen substituents affect base-stacking patterns. In general, the presence of a halogen substituent results in a stacking pattern which permits intimate contact between the halogen atom and adjacent purine or pyrimidine rings. Considering differences in the stacking patterns found for halogenated and nonhalogenated pyrimidines, a model is proposed to account for the mutagenic effects of halogenated pyrimidines.     

Interaction of a tryptophan-containing peptide with chromatin core particles. A fluorescence study

The binding of a tetrapeptide lysyltryptophylglycyllysine to nucleosome core particles has been investigated using UV absorption and fluorescence spectroscopy. Modifications of the absorption spectra and fluorescence quenching of the tryptophyl residue are consistent with stacking between the indole ring and nucleic acid bases. Therefore DNA interactions with histones do not prevent stacking of the tryptophyl residue with nucleic acid bases in the peptide-core particle complexes. The number of peptide binding sites is reduced to half that of naked DNA.     

Nonplanar DNA Base Pairs

Abstract

Three-dimensional structures of a representative set of more than 30 hydrogen-bonded nucleic acids base pairs have been studied by reliable ab initio quantum mechanical methods. We show that many hydrogen-bonded nucleic acid base pairs are intrinsically nonplanar, mainly due to the partial sp³ hybridization of nitrogen atoms of their amino groups and secondary electrostatic interactions. This finding extends the variability of intermolecular interactions of DNA bases in that i) flexibility of the base pairs is larger than has been assumed before, and ii) attractive proton-proton acceptor interactions oriented out of the base pair plane are allowed. For example, all four G…A mismatch base pairs are propeller twisted, and the energy preferences for the nonplanar structures range from less than 0.1 kcal/mol to 1.8 kcal/mol. We predict that nonplanarity of the amino group of guanine in the G(anti)…A(anti) pair of the ApG step of the d(CCAAGATTGG)₂ crystal structure is an important stabilizing factor that improves the energy of this structure by almost 3 kcal/mol. Currently used empirical potentials are not accurate enough to properly cover the interactions associated with amino-group and base-pair nonplanarity.     

Energy transfer in complexes of E. coli single-stranded DNA-binding protein with single-stranded poly-(2-thiouridylic acid)

The complexes of point-mutated Escherichia coli single-stranded DNA-binding protein (Eco SSB) with poly-(2-thiouridylic acid) (poly S2U) have been studied by optical detection of magnetic resonance spectroscopy (ODMR). Previous work has determined that two of four tryptophan (Trp) residues in Eco SSB undergo stacking interactions with nucleic acid bases. Selective photoexcitation of S2U bases was performed and subsequent triplet----triplet energy transfer from S2U to nearby Trp residues in the protein took place. The zero-field splitting (ZFS) parameters and sublevel kinetics were determined for each Trp residue sensitized by S2U. The sublevel lifetimes of the two sensitized residues are similar to those of normal Trp. The ZFS parameters, on the other hand, show a dramatic reduction relative to those of the uncomplexed protein, implying a more polarizable environment for the sensitized Trp residues and/or charge transfer interactions with the S2U bases.     

Thermal difference spectra: a specific signature for nucleic acid structures

We show that nucleic acid structures may be conveniently and inexpensively characterized by their UV thermal difference spectra. A thermal difference spectrum (TDS) is obtained for a nucleic acid by simply recording the ultraviolet absorbance spectra of the unfolded and folded states at temperatures above and below its melting temperature (T_m). The difference between these two spectra is the TDS. The TDS has a specific shape that is unique for each type of nucleic acid structure, a conclusion that is based on a comparison of >900 spectra from 200 different sequences. The shape of the TDS reflects the subtleties of base stacking interactions that occur uniquely within each type of nucleic acid structure. TDS provides a simple, inexpensive and rapid method to obtain structural insight into nucleic acid structures, which is applicable to both DNA and RNA from short oligomers to polynucleotides. TDS complements circular dichroism as a tool for the structural characterization of nucleic acids in solution.     

Stacking interaction of nucleobases: NMR investigations

The stacking interaction between nucleic acid bases has been investigated by the determination of the self-association of 6-methylpurine in various mixtures of water and nonaqueous solvents in order to elucidate the solvent effect. The parameters of stacking association as well as of local solvent-solute interactions have been measured by means of NMR technique. The influences of local hydration and of solvent-solvent interactions on the stacking ability are discussed.Parts I and II see Schimmack 1975     

Prominent stacking interaction with aromatic amino acid by N-quarternization of nucleic acid base: X-ray crystallographic characteristics and biological implications

In order to investigate the mode of interaction between the N-quarternized cytosine base and the aromatic amino acid, the crystal structure of the 3-methyl-cytidine-5'-monophosphate:tryptamine complex was analyzed by X-ray diffraction. The complex crystals were stabilized by extensive hydrogen bond formations in which eight independent water molecules per complex pair participated. A prominent stacking interaction, characterized by a parallel alignment of both rings with a separation distance of ca. 3.4 A, was observed between the cytosine base and the indole ring. Combining the present results with X-ray crystallographic data on the adenine--and guanine--aromatic amino acid interactions, we summarize the structural characteristics observed in the stacking interaction of the N-quarternized nucleic acid base with the aromatic amino acid and discuss their biological implications, especially in connection with the significance of N-protonation of nucleic acid base for selective recognition by protein.     

Beyond nucleic acid base pairs: from triads to heptads

Hydrogen-bonded base pairs are an important determinant of nucleic acid structure and function. However, other interactions such as base-base stacking, base-backbone, and backbone-backbone interactions as well as effects exerted by the solvent and by metal or NH(4)(+) ions also have to be taken into account. In addition, hydrogen-bonded base complexes involving more than two bases can occur. With the rapidly increasing number and structural diversity of nucleic acid structures known at atomic detail higher-order hydrogen-bonded base complexes, base polyads, have attracted much interest. This review provides an overview on the occurrence of base polyads in nucleic acid structures and describes computational studies on these nucleic acid building blocks.     

Disordered nucleiome: Abundance of intrinsic disorder in the DNA‐ and RNA‐binding proteins in 1121 species from Eukaryota,Bacteria and Archaea

Intrinsically disordered proteins (IDPs) are abundant in various proteomes, where they play numerous important roles and complement biological activities of ordered proteins. Among functions assigned to IDPs are interactions with nucleic acids. However, often, such assignments are made based on the guilty‐by‐association principle. The validity of the extension of these correlations to all nucleic acid binding proteins has never been analyzed on a large scale across all domains of life. To fill this gap, we perform a comprehensive computational analysis of the abundance of intrinsic disorder and intrinsically disordered domains in nucleiomes (～548 000 nucleic acid binding proteins) of 1121 species from Archaea, Bacteria and Eukaryota. Nucleiome is a whole complement of proteins involved in interactions with nucleic acids. We show that relative to other proteins in the corresponding proteomes, the DNA‐binding proteins have significantly increased disorder content and are significantly enriched in disordered domains in Eukaryotes but not in Archaea and Bacteria. The RNA‐binding proteins are significantly enriched in the disordered domains in Bacteria, Archaea and Eukaryota, while the overall abundance of disorder in these proteins is significantly increased in Bacteria, Archaea, animals and fungi. The high abundance of disorder in nucleiomes supports the notion that the nucleic acid binding proteins often require intrinsic disorder for their functions and regulation.     

Simulation of non-specific protein-mRNA interactions

Protein–nucleic acid interactions exhibit varying degrees of specificity. Relatively high affinity, sequence-specific interactions, can be studied with structure determination, but lower affinity, non-specific interactions are also of biological importance. We report simulations that predict the population of nucleic acid paths around protein surfaces, and give binding constant differences for changes in the protein scaffold. The method is applied to the non-specific component of interactions between eIF4Es and messenger RNAs that are bound tightly at the cap site. Adding a fragment of eIF4G to the system changes both the population of mRNA paths and the protein–mRNA binding affinity, suggesting a potential role for non-specific interactions in modulating translational properties. Generally, the free energy simulation technique could work in harness with characterized tethering points to extend analysis of nucleic acid conformation, and its modulation by protein scaffolds.