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.     

Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.     

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.     

Analysis of conformational parameters in nucleic acid fragments. II. Co-crystal complexes of nucleic acid bases.

Studies have been made of conformational parameters in co-crystal complexes and compounds of nucleic acid bases in which there is the possibility of formation of hetero-base-pairs. Using published data extracted from the Cambridge structural database, a total of 37 base-pairs were found, of which 25 were hetero-pairs and 12 homo-pairs. These base-pairs were subject to analysis to reveal hydrogen bond parameters, propeller twist, buckle and C1'-C1' separation (or a similar parameter if C1' atoms were not present). Hetero-pairs were found to show larger twists than homo-pairs, the magnitude of twist being unrelated to hydrogen bond parameters or buckle value. The propeller twisting is less pronounced in these nucleic acid bases than in nucleosides, but still has a significant magnitude. Propeller twisting in hetero-pairs is found to be larger than in homo-pairs. Hetero-pairs appear to be formed preferentially in competitive situations.     

Theoretical studies on the interaction of proteins and nucleic acid: II. The binding of α-helix to B-DNA

Interactions between B-DNA and homopolymeric α-helices of glycine, alanine, serine, asparagine and aspartic acid have been studied theoretically. The complexation energy has been minimised taking into account the interactions between DNA and the polypeptides as well as the internal energy of the α-helix and the interaction energy of counterions with the complex. The results obtained indicate the important role of strong hydrogen bonds between the peptide side chains and nucleic acid phosphate groups, these bonds being much stronger than specific interactions with the base-pairs. The formation of these structural bonds depends on the size of the α-helix, which in turn determines whether bridging across the major groove is possible. The steric role of the methyl group of thymine in orienting the peptide helix and the role of DNA screening cations in complex stabilization are also significant.     

Geometry and energetics of DNA basepairs and triplets from first principles quantum molecular relaxations.

A first principles model for calculating hydrogen bonding interactions, previously applied to water, is here applied to the more difficult problem of interactions between DNA bases. We first consider the energetics and geometry for the A-T and the G-C basepairs, comparing our results to other calculated results as well as to experiment. Next, we study the interactions of isomorphic DNA base triplet structures, which are important because of their suggested role in the recombination process. We find that energetically the third base in the triplet tends to favor a position along the dyadic axis, where it is hydrogen bonded to both bases in the duplex.     

Self-programmable,Self-assembling Two-dimensional Genetic Matter

Putative two-dimensional coding systems can beconstructed from aqueous solutions of purine andpyrimidine nucleic acid bases evaporated at moderatetemperatures on the surfaces of inorganic solids. Theresultant structures are monolayers which are formedspontaneously by molecular self-assembly and they havebeen observed with molecular resolution by scanningtunnelling microscopy (STM). When formed fromsolutions of a single base, the monolayers of adenineand uracil have crystalline characteristics and theSTM images can be interpreted in terms of thegeometrical placement of planar arranged moleculesthat interact laterally by intermolecular hydrogenbonding. When formed from solutions containing amixture of adenine and uracil, the monolayers haveaperiodic structures. Small crystalline domainswithin these monolayers can be interpreted in terms ofthe single phase configurations of the molecules andthe remaining aperiodic structures can presumably beinterpreted, geometrically, in terms of the 21theoretically possible adenine-adenine, uracil-uraciland adenine-uracil hydrogen bonding interactions. Wepropose that combinatorial arrangements of planararranged purine and pyrimidine bases could provide thenecessary complexity to act as a primitive geneticmechanism and may have relevance to the origin of life.     

Non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs, mostly present in the RNA, often play a prominent role towards maintaining their structural diversity. Higher order structures like base triples are also important in defining and stabilizing the tertiary folded structure of RNA. We have developed a new program BPFIND to analyze different types of canonical and non-canonical base pairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases or sugar O2' only. We considered 104 possible types of base pairs, out of which examples of 87 base pair types are found to occur in the available RNA crystal structures. Analysis indicates that approximately 32.7% base pairs in the functional RNA structures are non-canonical, which include different types of GA and GU Wobble base pairs apart from a wide range of base pair possibilities. We further noticed that more than 10.4% of these base pairs are involved in triplet formation, most of which play important role in maintaining long-range tertiary contacts in the three-dimensional folded structure of RNA. Apart from detection, the program also gives a quantitative estimate of the conformational deformation of detected base pairs in comparison to an ideal planar base pair. This helps us to gain insight into the extent of their structural variations and thus assists in understanding their specific role towards structural and functional diversity.     

Interaction and solvation energies of nonpolar DNA base analogues and their role in polymerase insertion fidelity.

Although DNA polymerase fidelity has been mainly ascribed to Watson-Crick hydrogen bonds, two nonpolar isosteres for thymine (T) and adenine (A)--difluorotoluene (F) and benzimidazole (Z) --effectively mimic their natural counterparts in polymerization experiments with pol I (KF exo-) [JC Morales and ET Kool. Nature Struct Biol, 5, 950-954, 1998]. By ab initio quantum chemical gas phase methods (HF/6-31G* and MP2/6-31G**) and a solvent phase method (CPCM-HF/6-31G**), we find that the A-F interaction energy is 1/3 the A-T interaction energy in the gas phase and unstable in the solvent phase. The F-Z and T-Z interactions are very weak and T-Z is quite unstable in the solvent. Electrostatic solvation energy calculations on F, Z and toluene yield that Z is two times, and F and toluene are five times, less hydrophilic than the natural bases. Of the new "base-pairs" (F-Z, T-Z, and F-A), only F-A formed an A-T-like arrangement in unconstrained optimizations. F-Z and T-Z do not freely form planar arrangements, and constrained optimizations show that large amounts of energy are required to make these pairs fit the exact A-T geometry, suggesting that the polymerase does not require all bases to conform to the exact A-T geometry. We discuss a model for polymerase/nucleotide binding energies and investigate the forces and conformational range involved in the polymerase geometrical selection.     

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.     

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.     

Dissecting the protein-RNA interface: the role of protein surface shapes and RNA secondary structures in protein-RNA recognition

Protein-RNA interactions are essential for many biological processes. However, the structural mechanisms underlying these interactions are not fully understood. Here, we analyzed the protein surface shape (dented, intermediate or protruded) and the RNA base pairing properties (paired or unpaired nucleotides) at the interfaces of 91 protein-RNA complexes derived from the Protein Data Bank. Dented protein surfaces prefer unpaired nucleotides to paired ones at the interface, and hydrogen bonds frequently occur between the protein backbone and RNA bases. In contrast, protruded protein surfaces do not show such a preference, rather, electrostatic interactions initiate the formation of hydrogen bonds between positively charged amino acids and RNA phosphate groups. Interestingly, in many protein-RNA complexes that interact via an RNA loop, an aspartic acid is favored at the interface. Moreover, in most of these complexes, nucleotide bases in the RNA loop are flipped out and form hydrogen bonds with the protein, which suggests that aspartic acid is important for RNA loop recognition through a base-flipping process. This study provides fundamental insights into the role of the shape of the protein surface and RNA secondary structures in mediating protein-RNA interactions.     

Direct evidence for methyl group coordination by carbon-oxygen hydrogen bonds in the lysine methyltransferase SET7/9

SET domain lysine methyltransferases (KMTs) are S-adenosylmethionine (AdoMet)-dependent enzymes that catalyze the site-specific methylation of lysyl residues in histone and non-histone proteins. Based on crystallographic and cofactor binding studies, carbon-oxygen (CH · · · O) hydrogen bonds have been proposed to coordinate the methyl groups of AdoMet and methyllysine within the SET domain active site. However, the presence of these hydrogen bonds has only been inferred due to the uncertainty of hydrogen atom positions in x-ray crystal structures. To experimentally resolve the positions of the methyl hydrogen atoms, we used NMR (1)H chemical shift coupled with quantum mechanics calculations to examine the interactions of the AdoMet methyl group in the active site of the human KMT SET7/9. Our results indicated that at least two of the three hydrogens in the AdoMet methyl group engage in CH · · · O hydrogen bonding. These findings represent direct, quantitative evidence of CH · · · O hydrogen bond formation in the SET domain active site and suggest a role for these interactions in catalysis. Furthermore, thermodynamic analysis of AdoMet binding indicated that these interactions are important for cofactor binding across SET domain enzymes.     

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.     

Interaction of nucleic acids and glycans

Spectrophotometric analysis and dot-hybridization have shown that amylose forms complexes with polypyrimidines (poly dC), while polyuronides form complexes with polypurines (poly dA). In addition, the formation of complexes genomic thymus DNA-hyaluronic acid has been observed. A certain role in the mechanism of NA-polysaccharide interactions can be played by the links between purines and the carboxylic group of hexuronic acid residue, as well as between pyrimidines and the hydroxymethyl group of hexose residue. The quantum-chemical calculations showed that, between nitric bases of DNA and the carboxyl groups of hexuronic acids or the hydroxymethyl group of hexose, hydrogen bonds can be formed the energy of which is comparable with that in the complementary AT and CG pairs. The strength of these bonds is unequal: carboxyl groups form stronger hydrogen bonds with purines and weaker bonds with pyrimidines. The hydroxymethyl group, on the contrary, forms stronger hydrogen bonds with pyrimidines and weaker bonds with purines. The quantum-chemical modeling shows that, in the complementary pairs purin-uronic acid and pyrimidine-hexose, hydrogen bonds are produced that form a binary chain nucleic acid-polysaccharide. The data obtained suggest the existence of template synthesis of GAG polysaccharide fragments with the participation of NA.     

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.     

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.