Conformation of dinucleoside phosphates: the conformation encoding hypothesis

The temperature dependence of the spin-lattice relaxation rate of nucleic bases protones and HI' of ApA, ApC, CpA and CpC (D2O, pH 7) were measured. The possible closed conformers of these dinucleoside phosphates (DNP) were computed by atom-atom potential method. On the basis of conformational calculation and experimental data the composition of closed state was determined. Besides the right-handed "canonic" conformers, the "non-canonic" right- and left-handed conformers were shown to be present in the solution of all DNP studied. It is important to note that, "canonic" conformers of DNP studied being equally probable, the possibility of the realization of "non-canonic" conformers is determined by the nucleotide sequence. It may be expected that different nucleotide sequences have unique "non-canonic" conformations. That type of dependence of the spatial organization of polynucleotides on its nucleotide sequence we call "the conformational encoding".     

Hydration of nucleic acid bases studied using novel atom-atom potential functions

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contains terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases - uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecular interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interactions play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by "bridges" formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

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.     

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.     

Hydration of Nucleic Acid Bases Studied Using Novel Atom-Atom Potential Functions

Abstract

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contain terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases—uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecule interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interaction play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by “bridges” formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

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.     

Base-pairing properties of O-methylated bases of nucleic acids. Energetic and steric considerations

Base-pairing properties of O-methylated nucleic-acid bases have been systematically investigated using both semi-empirical quantum-mechanical methods and a second-order perturbation formalism. The energetic, steric and electronic properties of (a) the individual methylated bases, (b) possible base-pairs formed between O-methylated and normal bases, and (c) mini-helices incorporating O-methylated bases were calculated. Two types of base-paired complexes were obtained: Those involving classical linear hydrogen bonds, and those involving bifurcated hydrogen-donor-hydrogen-acceptor interactions. In most complexes the presence of mispairs in the helical structure of nucleic acids is expected to create a local perturbation in the structure of the helix. Even though the most stable planar configurations of the mispairs may deviate markedly from those in the regular double helix, the induced deformations in the structure of the backbone are relatively small. Internal energies and geometries of mispairs are strongly affected by the conformation of the exocyclic group of the methylated bases. Another important contribution to the stability of various base-pairing schemes comes from stacking interactions.     

Interactions between nucleic acid bases. New parameters of potential functions and energy minima

The parameters of atom-atom potential functions suggested by one of the authors in 1979-1986 were slightly changed. The changes were performed to achieve a better agreement with experimental data of interaction energy values in global minima and hydrogen bond lengths. These changes resulted in better accord with experimental values of distances between the layers in DNA monomer crystals and between the base pairs in oligonucleotide duplexes. The refined potential functions were used to calculate the energy of interactions between nucleic acid bases in various mutual positions. The calculations revealed a few types of mutual base arrangements in minima of interaction energy for each pairwise base combination. A new type of minima was found, which correspond to a nearly perpendicular arrangement of base rings and the formation of the intermolecular hydrogen bond.     

Application of Raman spectroscopy to biological macromolecules.

The Raman spectra of biological macromolecules arise from molecular vibrations of either the backbone chains or the side chains. The frequencies of the Raman bands lie in a region between 200 cm-1 and 3000 cm-1. From certain frequencies of the vibrations of the backbone chains one can determine the conformation or secondary structure of a macromolecule. Thus for polypeptides and proteins the frequencies of the Amide I and Amide III vibrations allow one to determine the averge conformation of their backbone chain. In polynucleotides and nucleic acids, the frequency of the phosphate diester stretch of the phosphate furanose chain varies between 814 cm-1 for A conformation and 790 cm-1 for B conformation. Raman spectra of the bases in nucleic acids can be used to determine base stacking and hydrogen bonding interactions. Thus Raman spectroscopy is an important tool for determining the conformation structure of proteins and nucleic acids.     

Global structure of a DNA three-way junction by solution NMR: towards prediction of 3H fold

Three-way junctions (3H) are the simplest and most commonly occurring branched nucleic acids. They consist of three double helical arms (A to C), connected at the junction point, with or without a number of unpaired bases in one or more of the three different strands. Three-way junctions with two unpaired bases in one strand (3HS2) have a high tendency to adopt either of two alternative stacked conformations in which two of the three arms A, B and C are coaxially stacked, i.e. A/B-stacked or A/C-stacked. Empirical stacking rules, which successfully predict for DNA 3HS2 A/B-stacking preference from sequence, have been extended to A/C-stacked conformations. Three novel DNA 3HS2 sequences were designed to test the validity of these extended stacking rules and their conformational behavior was studied by solution NMR. All three show the predicted A/C-stacking preference even in the absence of multivalent cations. The stacking preference for both classes of DNA 3HS2 can thus be predicted from sequence. The high-resolution NMR solution structure for one of the stacked 3HS2 is also reported. It shows a well-defined local and global structure defined by an extensive set of classical NMR restraints and residual dipolar couplings. Analysis of its global conformation and that of other representatives of the 3H family, shows that the relative orientations of the stacked and non-stacked arms, are restricted to narrow regions of conformational space, which can be understood from geometric considerations. Together, these findings open up the possibility of full prediction of 3HS2 conformation (stacking and global fold) directly from sequence.     

Possible configuration of dimers of Gua-Cyt bases. Computation by molecular mechanic models and density functional theory

The energies of interactions between guanine and cytosine in various mutual positions were calculated by the methods of molecular mechanics with refined atom-atom potential functions and the quantum mechanics theory of density functional. Both methods indicate three types of mutual positions of bases in local energy minima. These types correspond to (1) nearly coplanar base positions with intermolecular hydrogen bond formation (base pairing); (2) arrangements of two bases in nearly parallel planes one above another (base stacking); and (3) nearly perpendicular positions of base planes. According to the calculations, the global energy minimum corresponds to the Watson-Crick base pair with three hydrogen bonds. A specific feature of the pair is a transition from many positions of type (2) to positions of type (1) without any energy barrier. This feature is revealed by both methods. Another special feature of this pair is a deviation, for most of mutual base positions, of the amine group atoms from the ring plane, the deviation being more pronounced for Gua. These features are important for understanding the conformational behavior of DNA fragments and the RNA structure.     

Finding and visualizing nucleic acid base stacking

Base stacking is one of the primary factors stabilizing nucleic acid structure. Yet, methods for locating stacking interactions in DNA and RNA are rare and methods for displaying stacking are rarer still. We present here simple, automated procedures to search nucleic acid molecules for base-base and base-oxygen stacking and to display these interactions graphically in a manner that readily conveys both the location and the quality of the interaction. The method makes no a priori assumptions about relative base positions when searching for stacking, nor does it rely on empirical energy functions. This is a distinct advantage for two reasons. First, the relative contributions of the forces stabilizing stacked bases are unknown. Second, the electrostatic and hydrophobic components of base stacking are both poorly defined by existing potential energy functions.     

Electronic spectra, excited state structures and interactions of nucleic acid bases and base assemblies: a review

A comprehensive review of recent theoretical and experimental advances in the singlet electronic transitions, excited state structures and dynamics of nucleic acid bases (NABs) and base assemblies are presented. It is well known that NABs absorb ultraviolet radiation, but the absorbed energy is efficiently dissipated in the form of ultrafast internal conversion processes believed to occur in the subpicosecond time scale and, therefore, enabling NABs highly photostable. It is not known how much evolutionary role was played in evolving these molecules and the ultimate selection by nature as genetic materials, but it is well accepted that survival-of-fittest prevails. Recently, significant efforts have been continuously paid to understand the mechanism of electronic excitation deactivation, but universally acceptable mechanism is still elusive. However, recent investigations reveal that electronic excited state geometries of DNA bases are usually nonplanar and this structural nonplanarity may facilitate nonradiative deactivation. Investigation of excited state structures is challenging and, therefore, it is not surprising that despite the impressive theoretical and computational advances, this research area is still hampered by the methodological and computational limitations. Further, stacking has significant influence on the emission properties of molecules. The 2-aminopurine, a fluorescent adenine derivative frequently used in studying DNA dynamics, shows significant attenuations in fluorescence quantum yield when incorporated in the DNA. Theoretical and computational bottlenecks limit a thorough theoretical understanding of effect of stacking interactions on the excited state dynamics of NABs. Despite these limitations the investigations of excited state properties are progressing in the right direction and our better understanding of excited state structure and dynamics of NABs and nucleic acids may help to design preventive strategy for radiation induced illness and photostable materials.     

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.     

Intramolecular stacking association and conformation properties of a ''cap'' structure, m7G5''pppUm, and the related model compounds

The stacking equilibrium quotient of the m7G5'pppUm unit, which occurs as the 5'-terminal "cap" of certain eukaryotic mRNA's, was determined by temperature-dependent difference spectrophotometry as Kstack = 1.82 at 25 degrees and pH 5. In order to evaluate the contribution of different structural modifications to the net stabilization of the cap structures of mRNA, a variety of compounds related to m7G5'pppUm were synthesized and their stacking properties were studied by the same method and compared. The results are summarized as: (1) Introduction of a methyl group into N-7 of G residue results in an increase in base stacking. (2) Methylation at 2'-OH of U residue also stabilizes the stacked structure of G-containing dimers, but it does not influence stacking interaction in m7G-containing dimers. (3) The effect of different types of internucleotide linkages on the order of stacking tendencies is: N5'ppN' greater than N5'pppN' greater than NpN'. UV hypochromicity and CD spectral measurements of the relevant dimers were also conducted, and the hypochromicity values and CD spectra of dimers in their stacked conformation were estimated by making use of the determined Kstack values. The results indicate that, while 2'-O-methylation exerts very little effect on the stacked conformation of the dimers, methylation at N-7 and the nature of the internucleotide linkage strongly influence the stacked conformation, thereby forming unusual left-handed conformations in m7G5'pppU(m), m7G5'ppU(m), and G5'ppU(m).     

Structural Dynamics of Human Telomeric G-Quadruplex Loops Studied by Molecular Dynamics Simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.     

Influence of base stacking and hydrogen bonding on the fluorescence of 2-aminopurine and pyrrolocytosine in nucleic acids

Fluorescent nucleobase analogues are used extensively to probe the structure and dynamics of nucleic acids. The fluorescence of the adenine analogue 2-aminopurine and the cytosine analogue pyrrolocytosine is significantly quenched when the bases are located in regions of double-stranded nucleic acids. To allow more detailed structural information to be obtained from fluorescence studies using these bases, we have studied the excited-state properties of the bases at the CIS and TDB3LYP level in hydrogen-bonded and base-stacked complexes. The results reveal that the first excited state (the fluorescent state) of a hydrogen-bonded complex containing 2-aminopurine and thymine is just the first excited state of 2-aminopurine alone. However, the same cannot be said for structures in which 2-aminopurine is base stacked with other nucleobases. Stacking causes the molecular orbitals involved in the fluorescence transition to spread over more than one base. The predicted rate for the fluorescence transition is reduced, thus reducing the fluorescence quantum yield. The decrease in radiative rate varies with the stacking arrangement (e.g., A- or B-form DNA) and with the identity of the nucleobase with which 2-aminopurine is stacked. Stacking 2-aminopurine between two guanine moieties is shown to significantly decrease the energy gap between the first and second excited states. We do not find reliable evidence for a low-energy charge-transfer state in any of the systems that were studied. In the case of pyrrolocytosine, base stacking was found to reduce the oscillator strength for the fluorescence transition, but very little spreading of molecular orbitals across more than one base was observed.