Analysis of conformational parameters in nucleic acid fragments. III. Very short chain oligonucleotides. The effect of base stacking

Studies have been made of conformational parameters in single crystal structures of very short chain oligonucleotides consisting of strands with lengths in the range 2-3 bases. Using published data extracted from the Cambridge structural database for 20 such structures, a total of 14 base-pairs were found, of which 10 were hetero-pairs and 4 homo-pairs. Subjecting these to analysis to examine hydrogen bond parameters, propeller twist, buckle and C1'-C1' separation revealed an average propeller twist of 11.6 degrees, with no dependence of this parameter on hydrogen bonding details. In addition, an analysis of base stacking showed there to be no correlation between in-plane geometry and degree of inter-plane overlap.     

Analysis of conformational parameters in nucleic acid fragments. I. Single crystals of nucleosides and nucleotides

A study has been undertaken of conformational parameters in single crystal structures of nucleosides and nucleotides using the techniques of helical conformational analysis. A "quasi-helix" was generated from the geometry of base-paired structures, using published data extracted from the Cambridge structural database. A total of 54 base-pairs were found in these structures, for each of which were calculated hydrogen bond parameters, propeller twist, buckle and C1'-C1' separations. These were analysed according to various classifications. Propeller twists are found to show a wide range of values and the magnitude of twist appears to be unrelated to hydrogen bond parameters or C1'-C1' separation. The values of the buckle parameter vary over a smaller range of values and are unrelated to propeller twist magnitude. There is found to be a greater tendency to form homo-base-pairs among compounds containing adenine bases.     

Analysis of conformational parameters in nucleic acid fragments. IV. Intercalating drug complexes of very short chain oligonucleotides.

Studies have been made of base-pairing conformational parameters in single crystal structures of very short chain oligonucleotide structures complexed with drug molecules, using data extracted from the Cambridge structural database. The planar portion of the drug has a tendency to intercalate between two bases, utilising strong stacking interactions to stabilise the configuration. The effect of the existence of a formative backbone is seen in the high occurrence of standard base-pairing schemes and a consistent C1'-C1' separation, although the choice of compounds studied does tend to emphasise complementary pairing. In addition to the modulation of the general magnitude which is reduced from that in uncomplexed oligonucleotides, there appears to be some correlation of propeller twist value with the presence of planar groups sandwiching a base-pair. The average twist in such sandwiched pairs is lower than in any other group studied to date.     

Structural diversity and isomorphism of hydrogen-bonded base interactions in nucleic acids

The wide structural diversity of RNA results in part from the diversity of non-Watson-Crick interactions between bases. To examine the repertoire of possible hydrogen bond interactions among bases, we computed databases of base-pairs and base-triples by systematically matching all possible hydrogen-bond donors and acceptors between bases and evaluating the geometries of each planar configuration. For base-pairs, we find 53 arrangements having at least two hydrogen bonds, including 23 pairs with protonated bases that have not previously been modeled. A comparison with experimentally observed base-pairs reveals an unexpected G:U pair recently observed in the ribosome. For base-triples, we find 840 arrangements in which the three bases are constrained by a total of at least three hydrogen bonds. Base-triples in particular exhibit a wide range of structural diversity, suggesting how compact or elongated nucleic acid structures may be constructed using different hydrogen-bonding patterns. Base-pair and base-triple conformations were systematically compared to identify structurally isomorphic combinations, and the experimentally observed arrangements within double and triple helices are among the most isomorphic. Unexpectedly, however, other combinations in the database are even more isomorphic, including several in which all-purine arrangements overlap with all-pyrimidine arrangements. These studies highlight some of the combinatoric and geometric versatility of base interactions and help provide a framework for analyzing and modeling isomorphic interactions and potentially for designing novel nucleic acid structures.     

Model studies of interactions between nucleic acids and proteins: hydrogen bonding of amides with nucleic acid bases.

The formation of hydrogen bonded complexes between nucleic acid bases and acetamide has been studied by nuclear magnetic resonance in CDC13 at different temperatures. Pairs of hydrogen bonds are formed when acetamide binds to nucleic acid bases. Thermodynamic parameters have been computed and compared to those obtained for the association of carboxylic acids with nucleic acid bases. The role of hydrogen bonded complexes in the association of proteins with nucleic acids is discussed.     

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.     

Propeller Twisting in Single Crystals of Nucleosides

Abstract

Propeller twists were measured for base-paired nucleosides and nucleic acid bases. The occurrence of propeller twist outside a double helical framework is shown, and conclusions drawn on the relative magnitude of the effect in single nucleosides.     

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

Recognition of nucleic acid bases and base-pairs by hydrogen bonding to amino acid side-chains

Sequence-specific protein-nucleic acid recognition is determined, in part, by hydrogen bonding interactions between amino acid side-chains and nucleotide bases. To examine the repertoire of possible interactions, we have calculated geometrically plausible arrangements in which amino acids hydrogen bond to unpaired bases, such as those found in RNA bulges and loops, or to the 53 possible RNA base-pairs. We find 32 possible interactions that involve two or more hydrogen bonds to the six unpaired bases (including protonated A and C), 17 of which have been observed. We find 186 "spanning" interactions to base-pairs in which the amino acid hydrogen bonds to both bases, in principle allowing particular base-pairs to be selectively targeted, and nine of these have been observed. Four calculated interactions span the Watson-Crick pairs and 15 span the G:U wobble pair, including two interesting arrangements with three hydrogen bonds to the Arg guanidinum group that have not yet been observed. The inherent donor-acceptor arrangements of the bases support many possible interactions to Asn (or Gln) and Ser (or Thr or Tyr), few interactions to Asp (or Glu) even though several already have been observed, and interactions to U (or T) only if the base is in an unpaired context, as also observed in several cases. This study highlights how complementary arrangements of donors and acceptors can contribute to base-specific recognition of RNA, predicts interactions not yet observed, and provides tools to analyze proposed contacts or design novel interactions.     

Affinity chromatography of nucleosides and nucleic acid base derivatives with nucleic acid bases or nitrobenzeneboronic acid substituted silicas

Nucleic acid bases such as adenine and uracil, and nitrobenzeneboronic acid substituted silicas were prepared by the reaction of chloromethylbenzene substituted silica with adenine sodium salt and trimethylsilylated uracil, and nitration of benzeneboronic acid substituted silica, respectively. From the results of HPLC of nucleosides and N-ethyl derivatives of nucleic acid bases using modified silicas, hydrophobic base stacking interaction, selective hydrogen bonding interaction between purine and pyrimidine bases, and reversible cyclic boronate ester formation between diols of nucleosides with boronic acid were effective for the separation of nucleic acid related compounds. Moreover, association constants for hydrogen bonding formation of nucleic acid bases were estimated.     

Analysis of local helix geometry in three B-DNA decamers and eight dodecamers

Local variations in B-DNA helix structure are compared among three decamers and eight dodecamers, which contain examples of all ten base-pair step types. All pairwise combinations of helix parameters are compared by linear regression analysis, in a search for internal relationships as well as correlations with base sequence. The primary conclusions are: (1) Three-center hydrogen bonds between base-pairs occur frequently in the major groove at C-C, C-A, A-A and A-C steps, but are less convincing at C-C and C-T steps in the minor groove. The requirements for large base-pair propeller are (1) that the base-pair should be A.T rather than G.C, and (2) that it be involved in a major groove three-center hydrogen bond with the following base-pair. Either condition alone is insufficient. Hence, a large propeller is expected at the leading base-pair of A-A and A-C steps, but not at A-T, T-A, C-A or C-C steps. (2) A systematic and quantitative linkage exists between helix variables twist, rise, cup and roll, of such strength that the rise between base-pairs can hardly be described as an independent variable at all. Two typical patterns of behavior are observed at steps from one base-pair to the next: high twist profile (HTP), characterized by high twist, low rise, positive cup and negative roll, and low twist profile (LTP), marked by low twist, high rise; negative cup and positive roll. Examples of HTP are steps G-C, G-A and Y-C-A-R, where Y is pyrimidine and R is purine. Examples of LTP steps are C-G, G-G, A-G and C-A steps other than Y-C-A-R. (3) The minor groove is especially narrow across the two base-pairs of the following steps: A-T, T-A, A-A and G-A. (4) In general, base step geometry cannot be correlated solely with the bases that define the step in question; the two flanking steps also must be taken into account. Hence, local helix structure must be studied in the context, not of two base-pairs: A-B, but of four: x-A-B-y. Calladine's rules, although too simple in detail, were correct in defining the length of sequence over which a given perturbation is expressed. Whereas ten different two-base steps are possible, allowing for the identity of complementary sequences, there are 136 different four-base steps. Only 33 of these 136 four-base steps are represented in the decamer and dodecamer structures solved to date, and hence it is premature to try to set up detailed structural algorithms. (5) The sugar-phosphate backbone chains of B-DNA place strong limits on sequence-induced structural variation, damping down most variables within four or five base-pairs, and preventing purine-purine anti-anti mismatches from causing bulges in the double helix. Hence, although short-range sequence-induced deformations (or deformability) are observed, long-range deformations propagated down the helix are not to be expected.     

Theoretical studies on protein-nucleic acid interactions. II. Hydrogen bonding of amino acid side chains with bases and base pairs of nucleic acids

With a view to understanding the role of hydrogen bonds in the recognition of nucleic acids by proteins, hydrogen bonding between the bases and base pairs of nucleic acids and the amino acids (Asn, Gln, Asp and Glu, and charged residues Arg⁺, Glu^?, and Asp^?) has been studied by a second-order perturbation theory. Binding energies have been calculated for all possible configurations involving a pair of hydrogen bonds between the base (or base pair) and the amino acid residue. Our results show that the hydrogen bonding in these cases has a large contribution from electrostatic interaction. In general, the charged amino acids, compared to the uncharged ones, form more stable complexes with bases or base pairs. The hydrogen-bond energies are an order of magnitude smaller than the Coulombic interaction energies between basic amino acids (Lys⁺, Arg⁺, and His⁺) and the phosphate groups of nucleic acids. The stabilities of the complexes of amino acids Asn, Gln, Asp, and Glu with bases are in the order: G–X > C–X > A–X U–X or T–X, and G · C–X > A · T(U)–X, where X is one of these amino acid residues. It has been shown that Glu^? and Asp^? can recognize guanine in single-stranded nucleic acids; Arg⁺ can recognize G · C base pairs from A · T base pairs in double-stranded structures.     

Polyethyleneimine derivative for nucleic acid model

Water-soluble polyethyleneimine (PE) derivatives containing nucleic acid bases and hydrophilic amino acids such as homoserine (Hse) and serine were prepared by the activated ester method as nucleic acid models. From spectroscopic measurements, the polymers were found to interact with DNA accompanied by an induction of conformational change. Hypochromicity in UV spectra indicated that a stable polymer complex was formed between poly (A) with PEI-Hse-Ura by complementary hydrogen bonding with equimolar nucleic base units (adenine∶uracil=1∶1). The induced conformation of DNA by the interaction with the polymer containing uracil and homoserine (PEI-Hse-Ura) was concluded to be a super triple helical structure. The formation of the polymer complex, DNA:PEI-Hse-Ura, was found to be affected by the presence of metal ions such as Ca²⁺ and Cu²⁺.     

Simulation of interactions between nucleic acid bases by refined atom-atom potential functions

Energy of interaction between nitrogen bases of nucleic acid has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H . . . O and (or) N-H . . . N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H . . . O or N-H . . . N and one C-H . . . O or C-H . . . N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum-mechanical approach.     

Conformational analysis of canonical 2-deoxyribonucleotides. 1. Pyrimidine nucleotides

The molecular structure and relative stability of north and south conformers of 2'-deoxyribonucleotides containing pyrimidine nucleic acid bases ( 2'-deoxythymidilic (pdT), 2'-deoxycytidilic (pdC) acids and their mono- and dianions) have been obtained and analyzed at the DFT/B3LYP level using the standard 6-31G(d) basis set. We have revealed that, when the nucleobase moiety is incorporated into the nucleotides, it maintains a nonplanar and nonrigid conformation due to out-of-plane deformation of the amino group and pyrimidine ring. It has been demonstrated that an increase of negative charge of the phosphate group results in increase of amino group pyramidalization, discrimination between conformers with syn and anti orientation of base with respect to sugar, strengthening of intramolecular C-H.O hydrogen bonds leading to deformation and fixation of geometry of nucleotides, and weakening of phosphodiester bond. These results allow to make suggestions about sources of twist and buckle deformations of base pairs, mechanisms of repaire of DNA via change of base orientation, and conditions for breakage of the P-O bonds during hydrolysis.     

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.     

Amide protein exchange and surface conformation of the basic pancreatic trypsin inhibitor in solution. Studies with two-dimensional nuclear magnetic resonance

Dickerson and his colleagues have described the structure of the DNA dodecamer C-G-C-G-A-A-T-T-C-G-C-G in the B form at a level that shows clearly several aspects of some base sequence-dependent departures from the ideal, regular helical structure of B-DNA. I argue that the detailed conformation is a consequence of simple steric repulsive forces between purine bases in consecutive base-pairs but on opposite backbones. These repulsions are a consequence of the “propeller twist” of the base-pairs, together with the larger size of the purine bases, and they may occur in either the major or the minor groove. The argument is conducted in terms of the structural mechanics of a deformable elastic system. These repulsive forces between the base-pairs are resisted by stresses in the helical backbones, which may be studied quantitatively via the variation in torsion angles δ along the backbones, at the points where the sugar rings are connected. There is also a correlation between the cross-chain purine repulsions and the perturbations in helical twist angle between successive base-pairs. The work suggests some comments on the proposed “alternating B” form, the Z form and the A form of DNA.     

Base sequence and helix structure variation in B and A DNA

The observed propeller twist in base-pairs of crystalline double-helical DNA oligomers improves the stacking overlap along each individual helix strand. But, as proposed by Calladine, it also leads to clash or steric hindrance between purines at adjacent base-pairs on opposite strands of the helix. This clash can be relieved by: (1) decreasing the local helix twist angle between base-pairs; (2) opening up the roll angle between base-pairs on the side on which the clash occurs; (3) separating purines by sliding base-pairs along their long axes so that the purines are partially pulled out of the stack (leading to equal but opposite alterations in main-chain torsion angle delta at the two ends of the base-pair); and (4) flattening the propeller twist of the offending base-pairs. Simple sum functions, sigma 1 through sigma 4, are defined, by which the expected local variation in helix twist, base roll angle, torsion angle delta and propeller twist may be calculated from base sequence. All four functions are quite successful in predicting the behavior of B DNA. Only the helix twist and base roll functions are applicable to A DNA, and the helix twist function begins to fail for an A helical RNA/DNA hybrid. Within these limits, the sequence-derived sum functions match the observed helix parameter variation quite closely, with correlation coefficients greater than 0.900 in nearly all cases. Implications of this sequence-derived helix parameter variation for repressor-operator interactions are considered.     

Ab-initio quantum mechanical calculations of NMR chemical shifts in nucleic acid constituents. I. The Watson-Crick base pairs

The magnetic shielding constants of the different nuclei of the four nucleic acid bases adenine, uracile, guanine and cytosine are calculated by a non empirical method using a minimal basis set and compared to the available corresponding experimental data. The same calculations carried out for AU and GC pairs give not only the values of the chemical shift variations due to the formation of the pairs but also the relative importance of the three different contributions (geometric, polarization and charge transfer plus exchange) to the total value of delta delta. Their analysis shows the importance of the polarization term. The magnitude of the charge transfer plus exchange term which is obtained for the nuclei belonging to the hydrogen bonding sites indicates that the hydrogen bond length is the major factor in the determination of the magnetic shielding constants of these atoms. On the other hand it appears that the pairing of the bases has a negligible effect on the "geometric" magnetic shielding due to the bases.     

Simulation of interactions in the co-planar nucleic acid base pairs using atom-atom potential functions

Calculations of the energy of nucleic acid base interactions as a function of parameters determining mutual position of two bases in a plane have been performed. Atom-atom potential functions used include terms proportional to the first (electrostatic), sixth (or tenth for the atoms of hydrogen bond) and 12th power of interatomic distance. The calculations have shown the existence of 27 energy minima which correspond to the formation of co-planar pairs with two (or three for G : C pair) almost linear N--H...O and N--H...N hydrogen bonds. The positions of nitrogen bases bound by two hydrogen bonds in every crystal of nucleic acid components, in the complexes of polynucleotides and in tRNA are near to the positions in one of these minima. In addition for every pair there exist energy minima which correspond to the formation of one N--H...O or N--H...N and one C--H...O or C--H...N hydrogen bond. Energy behavior near minima have been investigated. The results of our calculations are in agreement with experimental data and with the calculations which employ quantum mechanical results.