On the potential role of the amino nitrogen atom as a hydrogen bond acceptor in macromolecules

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Unconventional N-H…N Hydrogen Bonds Involving Proline Backbone Nitrogen in Protein Structures

Contrary to DNA double-helical structures, hydrogen bonds (H-bonds) involving nitrogen as the acceptor are not common in protein structures. We systematically searched N-H…N H-bonds in two different sets of protein structures. Data set I consists of neutron diffraction and ultrahigh-resolution x-ray structures (0.9 Å resolution or better) and the hydrogen atom positions in these structures were determined experimentally. Data set II contains structures determined using x-ray diffraction (resolution ≤ 1.8 Å) and the positions of hydrogen atoms were generated using a computational method. We identified 114 and 14,347 potential N-H…N H-bonds from these two data sets, respectively, and 56-66% of these were of the Ni+1-Hi+1…Ni type, with Ni being the proline backbone nitrogen. To further understand the nature of such unusual contacts, we performed quantum chemical calculations on the model compound N-acetyl-L-proline-N-methylamide (Ace-Pro-NMe) with coordinates taken from the experimentally determined structures. A potential energy profile generated by varying the ψ dihedral angle in Ace-Pro-NMe indicates that the conformation with the N-H…N H-bond is the most stable. An analysis of H-bond-forming proline residues reveals that more than 30% of the proline carbonyl groups are also involved in n → π(∗) interactions with the carbonyl carbon of the preceding residue. Natural bond orbital analyses demonstrate that the strength of N-H…N H-bonds is less than half of that observed for a conventional H-bond. This study clearly establishes the H-bonding capability of proline nitrogen and its prevalence in protein structures. We found many proteins with multiple instances of H-bond-forming prolines. With more than 15% of all proline residues participating in N-H…N H-bonds, we suggest a new, to our knowledge, structural role for proline in providing stability to loops and capping regions of secondary structures in proteins.

The exocyclic amino group of adenine in PtII and PdII complexes: a critical comparison of the X-ray crystallographic structural data and gas phase calculations

A detailed computational (DFT level of theory) study regarding the nature of the exocyclic amino group, N6H₂, of the model nucleobase 9-methyladenine (9MeA) and its protonated (9MeAH⁺) and deprotonated forms (9MeA_-H), free and metal-complexed, has been conducted. The metals are Pt^II and Pd^II, bonded to nitrogen-containing co-ligands (NH₃, dien, bpy), with N1, N6, and N7 being the metal-binding sites, individually or in different combinations. The results obtained from gas phase calculations are critically compared with X-ray crystallography data, whenever possible. In the majority of cases, there is good qualitative agreement between calculated and experimentally determined C6-N6 bond lengths, but calculated values always show a trend to larger values, by 0.02-0.08 Å. Both methods indicate, with few exceptions, a high degree of double-bond character of C6-N6, consistent with an essentially sp²-hybridized N6 atom. The shortest values for C6-N6 distances in X-ray crystal structures are around 1.30 Å. Exceptions refer to cases in which DFT calculations suggest the existence of a hydrogen bond with N6H₂ acting as a H bond acceptor, hence a situation with N6 having undergone a substantial hybridization shift toward sp³. Nevertheless, even in these cases the C6-N6 bond (1.392 Å) is still halfway between a typical C-N single bond (1.48 Å) and a typical C=N double bond (1.28 Å). This scenario is, however, not borne out by X-ray crystallographic results, and is attributed to the absence of counter anions and solvent molecules in the calculated structures.

N-H···N Hydrogen Bonds Involving Histidine Imidazole Nitrogen Atoms: A New Structural Role for Histidine Residues in Proteins

The amino acid histidine can play a significant role in the structure and function of proteins. Its various functions include enzyme catalysis, metal binding activity, and involvement in cation-π, π-π, salt-bridge, and other types of noncovalent interactions. Although histidine's imidazole nitrogens (Nδ and Nε) are known to participate in hydrogen bond (HB) interactions as an acceptor or a donor, a systematic study of N-H···N HBs with the Nδ/Nε atom as the acceptor has not been conducted. In this study, we have examined two data sets of ultra-high-resolution (data set I) and very high-resolution (data set II) protein structures and identified 28 and 4017 examples of HBs of the N-H···Nδ/Nε type from both data sets involving histidine imidazole nitrogen as the acceptor. In nearly 70% of them, the main-chain N-H bond is the HB donor, and a majority of the examples are from the N-H group separated by two residues (Ni+2-Hi+2) from histidine. Quantum chemical calculations using model compounds were performed with imidazole and N-methylacetamide, and they assumed conformations from 19 examples from data set I with N-H···Nδ/Nε HBs. Basis set superposition error-corrected interaction energies varied from -5.0 to -6.78 kcal/mol. We also found that the imidazole nitrogen of 9% of histidine residues forming N-H···Nδ/Nε interactions in data set II participate in bifurcated HBs. Natural bond orbital analyses of model compounds indicate that the strength of each HB is mutually influenced by the other. Histidine residues involved in Ni+2-Hi+2···Nδi/Nεi HBs are frequently observed in a specific N-terminal capping position giving rise to a novel helix-capping motif. Along with their predominant occurrence in loop segments, we propose a new structural role for histidines in protein structures.

Connectivity patterns and rotamer states of nucleobases determine acid-base properties of metalated purine quartets

Potentiometric pH titrations and pD dependent (1)H NMR spectroscopy have been applied to study the acidification of the exocyclic amino group of adenine (A) model nucleobases (N9 position blocked by alkyl groups) when carrying trans-a2Pt(II) (with a=NH3 or CH3NH2) entities both at N1 and N7 positions. As demonstrated, in trinuclear complexes containing central A-Pt-A units, it depends on the connectivity pattern of the adenine bases (N7/N7 or N1/N1) and their rotamer states (head-head or head-tail), how large the acidifying effect is. Specifically, a series of trinuclear complexes with (A-N7)-Pt-(N7-A) and (A-N1)-Pt-(N1-A) cross-linking patterns and terminal 9-alkylguanine ligands (9MeGH, 9EtGH) have been analyzed in this respect, and it is shown that, for example, the 9MeA ligands in trans-,trans-,trans-[Pt(NH3)2(N7-9MeA-N1)2{Pt(NH3)2(9EtGH-N7)}2](ClO4)6·6H2O (4a) and trans-,trans-,trans-[Pt(NH3)2(N7-9EtA-N1)2{Pt(CH3NH2)2(9-MeGH-N7)}2](ClO4)6·3H2O (4b) are more acidic, by ca. 1.3 units (first pKa), than the linkage isomer trans-,trans-,trans-[Pt(CH3NH2)2(N1-9MeA-N7)2{Pt(NH3)2(9MeGH-N7)}2](NO3)6·6.25H2O (1b). Overall, acidifications in these types of complexes amount to 7-9 units, bringing the pKa values of such adenine ligands in the best case close to the physiological pH range. Comparison with pKa values of related trinuclear Pt(II) complexes having different co-ligands at the Pt ions, confirms this picture and supports our earlier proposal that the close proximity of the exocyclic amino groups in a head-head arrangement of (A-N7)-Pt-(N7-A), and the stabilization of the resulting N6H(-)⋯H2N6 unit, is key to this difference.

Bifunctional homodimeric triokinase/FMN cyclase: contribution of protein domains to the activities of the human enzyme and molecular dynamics simulation of domain movements

Mammalian triokinase, which phosphorylates exogenous dihydroxyacetone and fructose-derived glyceraldehyde, is neither molecularly identified nor firmly associated to an encoding gene. Human FMN cyclase, which splits FAD and other ribonucleoside diphosphate-X compounds to ribonucleoside monophosphate and cyclic X-phosphodiester, is identical to a DAK-encoded dihydroxyacetone kinase. This bifunctional protein was identified as triokinase. It was modeled as a homodimer of two-domain (K and L) subunits. Active centers lie between K1 and L2 or K2 and L1: dihydroxyacetone binds K and ATP binds L in different subunits too distant (≈ 14 Å) for phosphoryl transfer. FAD docked to the ATP site with ribityl 4'-OH in a possible near-attack conformation for cyclase activity. Reciprocal inhibition between kinase and cyclase reactants confirmed substrate site locations. The differential roles of protein domains were supported by their individual expression: K was inactive, and L displayed cyclase but not kinase activity. The importance of domain mobility for the kinase activity of dimeric triokinase was highlighted by molecular dynamics simulations: ATP approached dihydroxyacetone at distances below 5 Å in near-attack conformation. Based upon structure, docking, and molecular dynamics simulations, relevant residues were mutated to alanine, and kcat and Km were assayed whenever kinase and/or cyclase activity was conserved. The results supported the roles of Thr(112) (hydrogen bonding of ATP adenine to K in the closed active center), His(221) (covalent anchoring of dihydroxyacetone to K), Asp(401) and Asp(403) (metal coordination to L), and Asp(556) (hydrogen bonding of ATP or FAD ribose to L domain). Interestingly, the His(221) point mutant acted specifically as a cyclase without kinase activity.

How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists

In this review primarily written for non-experts we explain basic methodological aspects and interpretation of modern quantum chemical (QM) computations applied to nucleic acids. We introduce current reference QM computations on small model systems consisting of dozens of atoms. Then we comment on recent advance of fast and accurate dispersion-corrected density functional theory methods, which will allow computations of small but complete nucleic acids building blocks in the near future. The qualitative difference between QM and molecular mechanics (MM, force field) computations is discussed. We also explain relation of QM and molecular simulation computations to experiments.

Recognition of asymmetrically dimethylated arginine by TDRD3

Asymmetric dimethylarginine (aDMA) marks are placed on histones and the C-terminal domain (CTD) of RNA Polymerase II (RNAP II) and serve as a signal for recruitment of appropriate transcription and processing factors in coordination with transcription cycle. In contrast to other Tudor domain-containing proteins, Tudor domain-containing protein 3 (TDRD3) associates selectively with the aDMA marks but not with other methylarginine motifs. Here, we report the solution structure of the Tudor domain of TDRD3 bound to the asymmetrically dimethylated CTD. The structure and mutational analysis provide a molecular basis for how TDRD3 recognizes the aDMA mark. The unique aromatic cavity of the TDRD3 Tudor domain with a tyrosine in position 566 creates a selectivity filter for the aDMA residue. Our work contributes to the understanding of substrate selectivity rules of the Tudor aromatic cavity, which is an important structural motif for reading of methylation marks.

Can We Accurately Describe the Structure of Adenine Tracts in B-DNA? Reference Quantum-Chemical Computations Reveal Overstabilization of Stacking by Molecular Mechanics

Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

A duplex DNA model with regular inter-base-pair hydrogen bonds

It is well known that base-pair stacking is the main factor in stabilizing DNA duplex and plays an important role in determining DNA sequence-dependence. What is the dominant force in base-pair stacking? This fundamental biological question remains a challenging problem. Here, based on recent studies about the non-planarity of amino groups on DNA bases, we propose a new duplex DNA model, in which all base amino groups are non-planar and participate in forming regular inter-base-pair hydrogen bonds (IBP H-bonds). This model implies that IBP H-bonds are the dominant force stabilizing base-pair stacking and play a crucial role in determining the geometry and physical properties of sequence-dependent twisted stacking between adjacent base pairs. The model presents a new insight into the link, through regular IBP H-bonds, between base-sequence, fine structure and physical properties at dinucleotide step level, and provides an attractively concise, uniform and quantitative interpretation for various experimentally observed DNA sequence-dependent properties in terms of regular IBP H-bonds. It would provide a new approach to understanding the dynamics and underlying mechanisms of DNA sequence-dependent biological processes, sequence-structure-property relationships, DNA strand separation during replication and transcriptions, etc.

On the many roles of NH3 ligands in mono- and multinuclear complexes of platinum

The role of the NH(3) ligands in the highly successful antitumour agents cisplatin and carboplatin is not fully understood. Suggestions that the ammonia ligands are involved in target recognition through hydrogen bond formation, e.g. with guanine-O6, have been questioned. Here, we review the roles and functions of NH(3) ligands of cis-PtCl(2)(NH(3))(2) and likewise of its trans-isomer in complexes with model nucleobases as well as other N-heterocyclic ligands. Specifically, their roles in hydrogen bonding interactions with nucleobases as well as anions, the influence on acid-base properties of co-ligands, their involvement in condensation reactions, as well as a variety of displacement reactions will be examined. As a result, it can be stated that the ammonia ligands in cis- and trans-Pt(II)(NH(3))(2) entities display additional features to those generally discussed in the last four decades since the discovery of the antitumour activity of cisplatin.

Trans Hoogsteen/sugar edge base pairing in RNA. Structures, energies, and stabilities from quantum chemical calculations

Trans Hoogsteen/sugar edge (H/SE) RNA base pairs form one of the six families of RNA base pairs that utilize the 2'-hydroxyl group of ribose for base pairing and play key roles in stabilizing folded RNA molecules. Here, we provide a detailed quantum chemical characterization of intrinsic structures and interaction energies of this base pair family, along with the evaluation of solvent screening effects by a continuum solvent approach. We report DFT-optimized geometries and MP2 interaction energies for all 10 crystallographically identified members of the family, for a representative set of them, using complete basis set extrapolation. For 6 of the 10 base pairs, we had to apply geometric constraints to keep the geometries relevant to RNA. We confirm that the remaining, hitherto undetected, possible members of this family do not have appropriate steric features required to establish stable base pairing in the trans H/SE fashion. The interaction patterns in the trans H/SE family are highly diverse, with gas-phase interaction energies in the range from -1 to -17 kcal/mol. Except for the C/rC and G/rG trans H/SE base pairs, the interaction energy is roughly evenly distributed between the HF and correlation components. Thus, in the trans H/SE base pairs, the relative importance of electron correlation is noticeably smaller than in the cis WC/SE or cis and trans SE/SE base pairs, but still larger than in canonical base pairs. The trans H/SE A/rG base pair is the intrinsically most stable member of this family. This base pair is also known as the sheared AG base pair and belongs to the most prominent set of RNA base pairs utilized in molecular building blocks of functional RNAs. For all trans H/SE base pairs that we identified, in addition to conventional base pairing, viable alternative structures were stabilized by amino-acceptor interactions. In the QM calculations, these amino-acceptor complexes appear to be equally as stable as those with common H-bonds, and more importantly, the switch to amino-acceptor interaction does not require any significant geometrical rearrangement of the base pairs. Such interactions are worthy of further investigations, as X-ray crystallography cannot unambiguously distinguish between conventional and amino-acceptor interactions involving the 2'-hydroxyl group, formation of such interactions is usually not considered, and molecular modeling force fields do not include such interactions properly as a result of neglect of aminogroup pyramidalization.

Confirmation of a second natural preQ1 aptamer class in Streptococcaceae bacteria

Bioinformatics searches of eubacterial genomes have yielded many riboswitch candidates where the identity of the ligand is not immediately obvious on examination of associated genes. One of these motifs is found exclusively in the family Streptococcaceae within the 5' untranslated regions (UTRs) of genes encoding the hypothetical membrane protein classified as COG4708 or DUF988. While the function of this protein class is unproven, a riboswitch binding the queuosine biosynthetic intermediate pre-queuosine(1) (preQ(1)) has been identified in the 5' UTR of homologous genes in many Firmicute species of bacteria outside of Streptococcaceae. Here we show that a representative of the COG4708 RNA motif from Streptococcus pneumoniae R6 also binds preQ(1). Furthermore, representatives of this RNA have structural and molecular recognition characteristics that are distinct from those of the previously described preQ(1) riboswitch class. PreQ(1) is the second metabolite for which two or more distinct classes of natural aptamers exist, indicating that natural aptamers utilizing different structures to bind the same metabolite may be more common than is currently known. Additionally, the association of preQ(1) binding RNAs with most genes encoding proteins classified as COG4708 strongly suggests that these proteins function as transporters for preQ(1) or another queuosine biosynthetic intermediate.

Tautomeric equilibrium, stability, and hydrogen bonding in 2'-deoxyguanosine monophosphate complexed with Mg2+

The tautomeric equilibrium and hydrogen bonding in nucleotide 2'-deoxyguanosine monophosphate that interacts with hydrated Mg2+ cation (4H2O.Mg[dGMP]) were studied at the MP2/cc-pVDZ//B3LYP/cc-pVDZ and B3LYP/aug-cc-pVTZ//B3LYP/cc-pVDZ levels of theory. The Mg2+ ion forms two inner-shell contacts with the nucleotide, similar to small phosphorylated molecules under physiological conditions. The presence of the phosphate group and the hydrated magnesium cation leads to a change in guanine tautomeric equilibrium of 4H2O.Mg[dGMP] in comparison to free guanine. The influence of the phosphate group and the magnesium cation on tautomeric equilibrium is larger in the anti conformation where the P=O-->Mg and Mg<--N7 coordinate bonds are formed. The canonical oxo form of guanine is more stable (by 6-8 kcal/mol) than the O6-hydroxo form in anti conformation. Thus, the interaction with Mg2+ ion is capable of further suppressing the likelihood of a spontaneous transient formation of the rare tautomer. In the syn conformation of 4H2O.Mg[dGMP], the interaction of the guanine nucleobase with the phosphate group and the magnesium cation is not as strong as in the anti conformation, and the relative stability of guanine tautomers is close to those in free guanine.

Molecular dynamics simulations and their application to four-stranded DNA

This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields. Several methodological approaches complementing the classical explicit solvent molecular dynamics simulations are commented on, including enhanced sampling methods, continuum solvent methods, free energy calculations and gas phase simulations. The successes and pitfalls of recent simulation studies of G-DNA are demonstrated on selected results, including studies of cation interactions and dynamics of G-DNA stems, studies of base substitutions (inosine, thioguanine and mixed tetrads), analysis of possible kinetic intermediates in folding pathway of a G-DNA stem and analysis of loop regions of G-DNA molecules.

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.

Non-Watson-Crick base pairing in RNA. quantum chemical analysis of the cis Watson-Crick/sugar edge base pair family

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/6-31G and RIMP2/aug-cc-pVDZ calculations. Further, a prediction is made for the remaining 3 cis WC/SE base pairs which are yet to be seen in the experiments. The interaction energy calculations point at the key role of the 2'-OH group in stabilizing the sugar-base contact and predict all 16 cis WC/SE base-pairing patterns to be nearly isoenergetic. The perfect correlation of the main geometrical parameters in the gas-phase optimized and X-ray structures shows that the principle of isosteric substitutions in RNA is rooted from the intrinsic structural similarity of the isolated base pairs. The present quantum chemical calculations for the first time analyze base pairs involving the ribose 2'-OH group and unambiguously correlate the structural information known from experiments with the energetics of interactions. The calculations further show that the relative importance and absolute value of the dispersion energy in the cis WC/SE base pairs are enhanced compared to the standard base pairs. This may by an important factor contributing to the strength of such interactions when RNA folds in its polar environment. The calculations further demonstrate that the Cornell et al. force field commonly used in molecular modeling and simulations provides satisfactory performance for this type of RNA interactions.

Sugar edge/sugar edge base pairs in RNA: stabilities and structures from quantum chemical calculations

Cis and trans sugar edge/sugar edge (SE/SE) binding patterns are essential building units of RNAs. For example, SE/SE interactions form the A-minor motifs, the most important tertiary interaction type in functional RNAs. This study provides an in-depth structure and stability analysis for these two base pair families. Gas-phase-optimized geometries are reported for 12 cis and 7 trans SE/SE base pairs and contrasted to their X-ray counterparts. Interaction energies are computed at the RIMP2 level of theory using the density-functional-theory-optimized geometries. There is a good overall agreement between the optimized and X-ray geometries of the cis SE/SE base pairs. In contrast, only three of the seven trans SE/SE binding patterns could be optimized without a significant distortion of the X-ray geometry. Note, however, that many SE/SE base pairs participate in broader networks of interactions; thus it is not surprising to see some of them to deviate from the X-ray geometry in a complete isolation. Computed interaction energies reveal that all 12 known cis SE/SE binding patterns are very stable. Among the trans SE/SE binding patterns, only the rG/rG, rG/rC, and rA/rG base pairs are sufficiently stable in the crystal geometry. Prediction has been made for some structures not yet detected by crystallography, namely, cis rC/rC, rG/rC, rG/rU, and rU/rU and trans rG/rA base pairs. Interestingly, the new cis SE/SE binding patterns are not necessarily isosteric with the remaining 12 members of this family. The trans rG/rA base pair represents a viable option for base pairing in RNA to be identified by future X-ray studies. In a complete lack of structural information, prediction of other unknown members of the trans SE/SE family was not attempted. Analysis of the interaction energies shows a very large electron correlation component of the interaction energy, pointing at the elevated role of dispersion energy as compared to other types of base pairs. This likely is profitable for stabilization of SE/SE binding patterns in polar environments and could be one of the reasons why the A-minor motif is the leading type of tertiary interactions in RNAs.

Principles of RNA base pairing: structures and energies of the trans Watson-Crick/sugar edge base pairs

Due to the presence of the 2'-OH hydroxyl group of ribose, RNA molecules utilize an astonishing variability of base pairing patterns to build up their structures and perform the biological functions. Many of the key RNA base pairing families have no counterparts in DNA. In this study, the trans Watson-Crick/sugar edge (trans WC/SE) RNA base pair family has been characterized using quantum chemical and molecular mechanics calculations. Gas-phase optimized geometries from density functional theory (DFT) calculations and RIMP2 interaction energies are reported for the 10 crystallographically identified trans WC/SE base pairing patterns. Further, stable structures are predicted for all of the remaining six possible members of this family not seen in RNAs so far. Among these novel six base pairs, the computations substantially refine two structures suggested earlier based on simple isosteric considerations. For two additional trans WC/SE base pairs predicted in this study, no arrangement was suggested before. Thus, our study brings a complete set of trans WC/SE base pairing patterns. The present results are also contrasted with calculations reported recently for the cis WC/SE base pair family. The computed base pair sizes are in sound correlation with the X-ray data for all WC/SE pairing patterns including both their cis and trans isomers. This confirms that the isostericity of RNA base pairs, which is one of the key factors determining the RNA sequence conservation patterns, originates in the properties of the isolated base pairs. In contrast to the cis structures, however, the isosteric subgroups of the trans WC/SE family differ not only in their H-bonding patterns and steric dimensions but also in the intrinsic strength of the intermolecular interactions. The distribution of the total interaction energy over the sugar-base and base-base contributions is controlled by the cis-trans isomerism.

Varying Acidity of Aqua Ligands in Dependence on the Microenvironment in Mononucleobase (nb) Complexes of Type cis- and trans-[Pt(NH3)2(nb)(H2O)]n+

Aqua ligands in mixed aqua/nucleobase metal complexes are potential sites of acid-base catalysis and/or, when present as hydroxo ligands, can directly be involved in hydrolysis reactions. pKa values of close to 7 are consequently of particular interest and potential significance. Here we report on the differential acidity of aqua complexes in model nucleobase (nb) complexes of cis- and trans-[Pt(NH3)2 (nb)(H2O)]n+ and discuss reasons as to why the nb in cis complexes influences the pKa (pKa 4.8-7.0), whereas in trans complexes the pKa values are rather constant (pKa approximately 5.2-5.3). The results of DFT calculations of a series of mono(nucleobase) complexes derived from cis-Pt(NH3)2 are critically examined with regard to the role of exocyclic groups of nucleobases in stabilizing aqua/hydroxo ligands through intracomplex hydrogen bond formation. This applies in particular to the exocyclic amino groups of nucleobases, for which gas-phase calculations suggest that they may act as H bond acceptors in certain cases, yet in the condensed phase this appears not to be the case.

Structure of a DNA repair substrate containing an alkyl interstrand cross-link at 1.65 A resolution

Chemotherapeutic alkylating agents, such as bifunctional nitrogen mustards and cisplatins, generate interstrand DNA cross-links that inhibit cell proliferation by arresting DNA transcription and replication. A synthetic N4C-ethyl-N4C interstrand cross-link between opposing cytidines mimics the DNA damage produced by this class of clinically important compounds and can be synthesized in large quantities to study the repair, physical properties, and structures of these DNA adducts. The X-ray structure of a DNA duplex d(CCAAC*GTTGG)2 containing a synthetic N4C-ethyl-N4C interstrand cross-link between the cytosines of the central CpG step (*) has been determined at 1.65 A resolution. This structure reveals that the ethyl cross-link in the CpG major groove does not significantly disrupt the B-form DNA helix. Comparison of the N4C-ethyl-N4C cross-linked structure with the structure of an un-cross-linked oligonucleotide of the same sequence reveals that the cross-link selectively stabilizes a preexisting alternative conformation. The conformation preferred by the cross-linked DNA is constrained by the geometry of the ethyl group bridging the cytosine amines. Characteristics of the cross-linked CpG step include subtle differences in the roll of the base pairs, optimized Watson-Crick hydrogen bonds, and loss of a divalent cation binding site. Given that the N4C-ethyl-N4C cross-link stabilizes a preexisting conformation of the CpG step, this synthetically accessible substrate presents an ideal model system for studying the genomic effects of covalently coupling the DNA strands, independent of gross alterations in DNA structure.

The small planarization barriers for the amino group in the nucleic acid bases

The amino group in the nucleic acid bases frequently interacts with other bases or with other molecular systems. Thus any nonplanarity of the amino group may affect the molecular recognition of nucleic acids. Ab initio Hartree-Fock (HF) and second-order Moller-Plesset perturbation (MP2) levels of theory have been used to obtain the equilibrium geometries of the C(l) and C(s) structures for five common nucleic acid bases. The energy barriers between the C(l) and C(s) structures have also been predicted. A series of correlation consistent basis sets up to cc-pCVQZ and aug-cc-pVQZ has been used to systematically study the dependence of the amino group nonplanarity. The equilibrium geometries of the nucleic acid bases with an amino group, including adenine, guanine, and cytosine, are examined carefully. At the MP2 level of theory, larger basis sets decrease the extent of nonplanarity of the amino group, but the decrease slows down when the QZ basis sets are used, demonstrating the intrinsic property of nonplanarity for guanine. For adenine and cytosine the situation is less clear; as the HF limit is approached, these two structures become planar. Addition of core correlation effects or diffuse functions further decreases the degree of nucleic acid base nonplanarity, in comparison to the original cc-pVXZ (X=D, T, and Q) basis sets. The aug-cc-pVXZ basis shows smaller degrees of nonplanarity than the cc-pCVXZ sets. The aug-cc-pVXZ basis is less size dependent than the cc-pVXZ and cc-pCVXZ sets in the prediction of the amino-group-related bond angles and dihedral angles and energy barriers for adenine, guanine, and cytosine. The cc-pCVQZ and aug-cc-pVQZ MP2 results may be regarded as benchmark predictions for the five common bases. The predicted classical barriers to planarization are 0.02 (adenine), 0.74 (guanine), and 0.03(cytosine) kcal mol(-1).

Netropsin interactions in the minor groove of d(GGCCAATTGG) studied by a combination of resolution enhancement and ab initio calculations

The structure of the complex between the minor groove binder netropsin and d(GGCCAATTGG) was determined via single-crystal X-ray techniques. The structure was refined to completion using refmac5.1.24, resulting in a residual R-factor of 20.0% (including 68 water molecules). Using crystal engineering and cryocooling techniques, the resolution could be enhanced to 1.75 A, resulting in an unambiguous determination of the drug conformation and orientation. As previously noticed, bifurcated hydrogen bonds are formed between the amide nitrogen atoms of the drug and the N3 and O2 atoms of A and T base pairs, respectively, clearly cataloging the structure to class I. As the bulky NH2 group on guanine was believed to prevent binding of the drug in the minor groove, the detailed nature of several of the amidinium and guanidinium end contacts were further investigated by ab initio quantum chemical methods.

Comprehensive Density Functional Theory Study on Serine and Related Ions in Gas Phase: Conformations, Gas Phase Basicities, and Acidities

Density functional theory (DFT) calculations have been performed to investigate the gas-phase conformations of serine and its three related ions (serineH(+), serine(-), and serine(2-)). The full ensemble of possible conformations, 324 conformations for serine, 108 for serineH(+), 162 for serine(-) and 54 for serine(2-), were first surveyed at B3LYP/6-31G level, and then the obtained unique conformations were further refined at B3LYP/6-311+G level. From full optimizations, 74 unique conformations for seine, 14 for serineH(+), 11 for serine(-), and 4 for serine(2-) were located, and their relative energies were also determined at B3LYP/6-311+G level. Atoms in molecules (AIM) analysis was carried out to establish rigorous definition of hydrogen bonds. Six types of intramolecular H-bonds in conformers of serine, six types in serineH(+), three types in serine(-), and two types in serine(2-) were identified within the framework of AIM theory and their relative strengths were determined based on topological properties at bond critical points (BCPs) of H-bonds. The intramolecular H-bonds were demonstrated to play an important role in deciding the relative stability of conformations of amino acids and the related ions. The enthalpies and Gibbs free energies of protonation and deprotonation reactions of serine and its related ions were calculated at B3LYP/6-311+G//B3LYP/6-31G, and B3LYP/6-311+G//B3LYP/6-311+G level. The calculated results are both in excellent agreement with the experimental data. We demonstrate in this study that B3LYP is an efficient and accurate method to predict the thermochemical and structural parameters of amino acids and the related ions.

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg2+ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na+ ions with residency times of individual cations 1-2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.

Complex Formation of Isocytosine Tautomers with PdII and PtII

Isocytosine (ICH) exists in solution as two major tautomers, the keto form with N1 carrying a proton (1a) and the keto form with N3 being protonated (1b). In water, 1a and 1b exist in equilibrium with almost equal amounts of both forms present. Reactions with a series of Pd(II) and Pt(II) am(m)ine species such as (dien)Pd(II), (dien)Pt(II), and trans-(NH(3))(2)Pt(II) reveal, however, a distinct preference of these metals for the N3 site, as determined by (1)H NMR spectroscopy. Individual species have been identified by the pD dependence of the ICH resonances. pK(a) values (calculated for H(2)O) for deprotonation of the individual tautomers complexes are 6.5 and 6.4 for the N3 linkage isomers of dienPd(II) and dienPt(II), respectively, as well as 6.2 and 6.0 for the N1 linkage isomers. The dimetalated species [(dienM)(2)(IC-N1,N3)](3+) (M = Pd(II) or Pt(II)) are insensitive over a wide range of pD. The crystal structure analysis of [(dien)Pd(ICH-N3)](NO(3))(2) is reported. Ab initio calculations have been performed for tautomer compounds of composition [(NH(3))(3)Pt(ICH)](2+), cis- and trans-[(NH(3))(2)PtCl(ICH)](+), as well as trans-[(NH(3))(2)Pt(ICH)(2)](2+). Without exception, N3 linkage isomers are more stable, in agreement with experimental findings. As to the reasons for this binding preference, an NBO (natural bond orbital) analysis for [(NH(3))(3)Pt(ICH-N3)](2+)strongly suggests that intramolecular hydrogen bonding between trans-positioned NH(3) ligands and the two exocyclic groups of the ICH is of prime importance. The calculations furthermore show a marked pyramidalization of the NH(2) group of ICH in the complex once the heterocyclic ligand forms a dihedral angle <90 degrees with the Pt coordination plane.

New parameterization of the Cornell et al. empirical force field covering amino group nonplanarity in nucleic acid bases

The Cornell et al. empirical potential (Cornell et al., J Am Chem Soc 1995, 117, 5197) was modified by the introduction of nonplanarity of the amino group in guanine, adenine, and cytosine. Reparameterization was performed for 12 bond parameters of the amino group (three valence angles (C-N-H(1), C-N-H(2), and H(1)-N-H(2)) and the improper dihedral angle at nitrogen, four dihedral angles (X-C-N-H(1), Y-C-N-H(2), X-X-C-N, Y-Y-C-N), three valence angles (X-C-N, Y-C-N, X-C-Y), and the improper dihedral angle at the adjacent carbon), and was based on correlated ab initio potential energy surfaces. Calculations were performed using the resolution of identity MP2 (RIMP2) method with SVP (3s2p1d/2s1p), TZVP (5s3p1d/3s1p), TZVPP (5s3p2d1f/3s2p1d), and augTZVPP (6s4p3d2f/4s3p2d) basis sets. Results obtained on the latter two levels are practically identical with the literature reference data [MP2/6-311++G(2df,p)] and the RIMP2/augTZVPP data can be considered a new reference set. The potential energy surface used for reparameterization was evaluated at the RIMP2/SVP level because respective geometry as well as energy data were close to the reference ones, and the computational time was very favorable. The modified potential was tested for 19 H-bonded and eight stacked nucleic acid base pairs. Interaction energies as well as geometries were described by the modified potential considerably better than by the original one and significant improvement resulted in the description of the nonplanar H-bonded and stacked complexes.

Theoretical Study of Cisplatin Binding to Purine Bases: Why Does Cisplatin Prefer Guanine over Adenine?

The thermodynamics and kinetics for the monofunctional binding of the antitumor drug cisplatin, cis-diamminedichloroplatinum(II), to a purine base site of DNA were studied computationally using guanine and adenine as model reactants. A dominating preference for initial attack at the N7-position of guanine is established experimentally, which is a crucial first step for the formation of a 1,2-intrastrand cross-link of adjacent guanine bases that leads to bending and unwinding of DNA. These structural distortions are proposed ultimately to be responsible for the anticancer activity of cisplatin. Utilizing density functional theory in combination with a continuum solvation model, we developed a concept for the initial Pt-N7 bond formation to atomic detail. In good agreement with experiments that suggested DeltaG++ = approximately 23 kcal/mol for the monofunctional platination of guanine, our model gives DeltaG++ = 24.6 kcal/mol for guanine, whereas 30.2 kcal/mol is computed when adenine is used. This result predicts that guanine is 3-4 orders of magnitude more reactive toward cisplatin than adenine. A detailed energy decomposition and molecular orbital analysis was conducted to explain the different barrier heights. Two effects are equally important to give the preference for guanine over adenine: First, the transition state is characterized by a strong hydrogen bond between the ammine-hydrogen of cisplatin and the O=C6 moiety of guanine in addition to a stronger electrostatic interaction between the two reacting fragments. When adenine binds, only a weak hydrogen bond forms between the chloride ligand of cisplatin and the H(2)N-C6 group of adenine. Second, a significantly stronger molecular orbital interaction is identified for guanine compared to adenine. A detailed MO analysis is presented to provide an intuitive view into the different electronic features governing the character of the Pt-N7 bond in platinated purine bases.

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.