DNA deformability at the base pair level

A complete set of harmonic force constants describing the DNA deformation energetics at the base pair level was obtained using unrestrained atomic-resolution molecular dynamics simulations of selected duplex oligonucleotides and subsequent analysis of structural fluctuations from the simulated trajectories. The deformation was described by the six base pair conformational parameters (buckle, propeller, opening, shear, stretch, stagger). The results for 13 AT pairs and 11 GC pairs in different sequence contexts suggest that buckle and propeller are very flexible (more than roll in TA dinucleotide steps), while stretch is exceptionally stiff. Only stretch and opening stiffness were found to depend unambiguously on the base pair identity (AT vs GC). The relationship of the results to a simple plates-and-springs model of base-base interactions is discussed.     

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.     

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.     

Structure, stability, and dynamics of canonical and noncanonical base pairs: quantum chemical studies

The importance of non-Watson-Crick base pairs in the three-dimensional structure of RNA is now well established. The structure and stability of these noncanonical base pairs are, however, poorly understood. We have attempted to understand structural features of 33 frequently occurring base pairs using density functional theory. These are of three types, namely (i) those stabilized by two or more polar hydrogen bonds between the bases, (ii) those having one polar and another C-H...O/N type interactions, and (iii) those having one H-bond between the bases and another involving one of the sugars linked to the bases. We found that the base pairs having two polar H-bonds are very stable as compared to those having one C-H...O/N interaction. Our quantitatively analysis of structures of these optimized base pairs indicates that they possess a different amount of nonplanarity with large propeller or buckle values as also observed in the crystal structures. We further found that geometry optimization does not modify the hydrogen-bonding pattern, as values of shear and open angle of the base pairs remain conserved. The structures of initial crystal geometry and final optimized geometry of some base pairs having only one polar H-bond and a C-H...O/N interaction, however, are significantly different, indicating the weak nature of the nonpolar interaction. The base pair flexibility, as measured from normal-mode analysis, in terms of the intrinsic standard deviations of the base pair structural parameters are in conformity with those calculated from RNA crystal structures. We also noticed that deformation of a base pair along the stretch direction is impossible for all of the base pairs, and movements of the base pairs along shear and open are also quite restricted. The base pair opening mode through alteration of propeller or buckle is considerably less restricted for most of the 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.     

Molecular Mechanics Investigation of an Adenine-Adenine Non-Canonical Pair Conformational Change

Conformational changes are important in RNA for binding and catalysis and understanding these changes is important for understanding how RNA functions. Computational techniques using all-atom molecular models can be used to characterize conformational changes in RNA. These techniques are applied to an RNA conformational change involving a single base pair within a nine base pair RNA duplex. The Adenine-Adenine (AA) non-canonical pair in the sequence 5'GGUGAAGGCU3' paired with 3'PCCGAAGCCG5', where P is Purine, undergoes conformational exchange between two conformations on the timescale of tens of microseconds, as demonstrated in a previous NMR solution structure [Chen, G., et al., Biochemistry, 2006. 45: 6889-903]. The more populated, major, conformation was estimated to be 0.5 to 1.3 kcal/mol more stable at 30 °C than the less populated, minor, conformation. Both conformations are trans-Hoogsteen/sugar edge pairs, where the interacting edges on the adenines change with the conformational change. Targeted Molecular Dynamics (TMD) and Nudged Elastic Band (NEB) were used to model the pathway between the major and minor conformations using the AMBER software package. The adenines were predicted to change conformation via intermediates in which they are stacked as opposed to hydrogen-bonded. The predicted pathways can be described by an improper dihedral angle reaction coordinate. Umbrella sampling along the reaction coordinate was performed to model the free energy profile for the conformational change using a total of 1800 ns of sampling. Although the barrier height between the major and minor conformations was reasonable, the free energy difference between the major and minor conformations was the opposite of that expected based on the NMR experiments. Variations in the force field applied did not improve the misrepresentation of the free energies of the major and minor conformations. As an alternative, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approximation was applied to predict free energy differences between the two conformations using a total of 800 ns of sampling. MM-PBSA also incorrectly predicted the major conformation to be higher in free energy than the minor conformation.     

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

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.     

Addition of lithiated 5-hydroxymethyl-1,3-dithiane to benzaldehyde: HMPA-controlled trans stereoselectivity

Addition of 5-substituted dithianyl anions to carbonyl compounds normally produces trans adducts. The presence of a nucleophilic hydroxymethyl group in position 5 dramatically decreases the trans stereoselectivity of the reaction in THF. The trans/cis ratio shows a bell curve dependence on HMPA, fitted to a quantitative model involving a series of equilibrated ion pairs, of which an intermediate contact ion pair possessing three (effective) HMPA molecules yields the trans adduct with much higher stereoselectivity. [structure: see text]     

First-principle computational study on the full conformational space of L-threonine diamide, the energetic stability of cis and trans isomers

First-principle computations were carried out on the conformational space of trans and cis peptide bond isomers of HCO-Thr-NH2. Using the concept of multidimensional conformational analysis (MDCA), geometry optimizations were performed at the B3LYP/6-31G(d) level of theory, and single-point energies as well as thermodynamic functions were calculated at the G3MP2B3 level of theory for the corresponding optimized structures. Two backbone Ramachandran-type potential energy surfaces (PESs) were computed, one each for the cis and trans isomers, keeping the side chain at the fully extended orientation (chi1=chi2=anti). Similarly, two side chain PESs for the cis and trans isomers were generated for the (phi=psi=anti) orientation corresponding to approximately the betaL backbone conformation. Besides correlating the relative Gibbs free energy of the various stable conformations with the number of stabilizing hydrogen bonds, the process of trans-->cis isomerization is discussed in terms of intrinsic stabilities as measured by the computed thermodynamic functions.     

Analysis of stacking overlap in nucleic acid structures: algorithm and application

RNA contains different secondary structural motifs like pseudo-helices, hairpin loops, internal loops, etc. in addition to anti-parallel double helices and random coils. The secondary structures are mainly stabilized by base-pairing and stacking interactions between the planar aromatic bases. The hydrogen bonding strength and geometries of base pairs are characterized by six intra-base pair parameters. Similarly, stacking can be represented by six local doublet parameters. These dinucleotide step parameters can describe the quality of stacking between Watson–Crick base pairs very effectively. However, it is quite difficult to understand the stacking pattern for dinucleotides consisting of non canonical base pairs from these parameters. Stacking interaction is a manifestation of the interaction between two aromatic bases or base pairs and thus can be estimated best by the overlap area between the planar aromatic moieties. We have calculated base pair overlap between two consecutive base pairs as the buried van der Waals surface between them. In general, overlap values show normal distribution for the Watson–Crick base pairs in most double helices within a range from 45 to 50 Å² irrespective of base sequence. The dinucleotide steps with non-canonical base pairs also are seen to have high overlap value, although their twist and few other parameters are rather unusual. We have analyzed hairpin loops of different length, bulges within double helical structures and pseudo-continuous helices using our algorithm. The overlap area analyses indicate good stacking between few looped out bases especially in GNRA tetraloop, which was difficult to quantitatively characterise from analysis of the base pair or dinucleotide step parameters. This parameter is also seen to be capable to distinguish pseudo-continuous helices from kinked helix junctions.     

Photoinduced conformational switch of enantiopure azobenzenes controlled by a sulfoxide

Two series of enantiopure azobenzenes with a p-tolylsulfoxide at the ortho or meta position with respect to the azo group, have been regioselectively synthesized. Both can act as enantiopure molecular switches showing different structural features owing to the presence of the stereogenic sulfur. The photoisomerization process, studied by UV-vis, circular dichroism (CD), NMR, and chiral HPLC evidenced a double role of the sulfoxide. A transfer of chirality from the sulfoxide to the azo system was observed by CD in both cis and trans-isomers of the meta sulfinyl derivatives 3, whereas this perturbation was evident for the ortho sulfinyl series 7 only in the cis isomer. The NMR study evidenced that the s-cis rigid conformation of the bisaromatic sulfoxide was fixing a different orientation of the overall system in each series both in the trans and cis isomers, by forcing a final U-shaped structure in cis-3 and an S-shaped structure in cis-7. Very different values of specific optical rotations were measured in both trans and cis isomers, also reflecting the existence of distinct chiral entities in the photostationary states. The easy and reversible changes occurring between different conformational states could find applications in the photocontrol of several molecular switches.     

Conformational study of tyramine and its water clusters by laser spectroscopy

Tyramine and its monohydrated clusters have been investigated by several laser spectroscopic methods in a pulsed molecular beam. The conformational structures and their effects on hydration have been revealed by resonant two-photon ionization (R2PI), UV-UV ion-dip, and ab initio calculations. UV rotational band contour spectra of the S1 <-- S0 origin bands enabled determination of ethylamine side chain conformations for all seven stable conformers of tyramine. When coexpanding tyramine with a mixture of Ar and water vapor, we have found two kinds of conformational effects on hydration. One is sensitive to conformation of the ethylamine chain and the other to the orientation of the OH group, particularly in the most stable pair of conformers. UV-UV ion-dip spectra detected seven stable conformers of the monohydrated clusters, of which hydrogen-bonding structures, spectral shifts, and origin band intensity distributions are well explained by considering tyramine as a hybrid of phenylethylamine (PEA) and phenol. Monohydration of the most stable gauche conformer pair (cis and trans) of tyramine leads to more detailed conformational assignments regarding the orientation of the phenolic OH group. Cyclic hydrogen-bonding linkage formed in the monohydrated cluster pair is found to be sensitive to the orientation of the phenolic OH group. One of the cluster pair, in which tyramine has the gauche-cis conformation, is more stabilized by the cyclic hydrogen bonding and its origin band intensity becomes stronger than that of the other.     

Parallel DNA double helices incorporating isoG or m5isoC bases studied by FTIR, CD and molecular modeling

FTIR spectroscopy has been used to follow the formation of parallel stranded DNA duplexes incorporating isoG or m5isoC bases and determine their base pairing scheme. The results are discussed in comparison with data concerning anti-parallel duplexes with comparable base composition and sequence. In duplexes containing A-T and isoG-C or m5isoC-G base pairs shifts of the thymine C2=O2 and C4=O4 carbonyl stretching vibrations (to lower and higher wavenumbers, respectively, when compared to their positions in classical cis Watson-Crick (WC) base pairs) reflect the formation of trans Watson-Crick A-T base pairs. All carbonyl groups of cytosines, m5isocytosines, guanines and isoguanines are found to be involved in hydrogen bonds, indicative of the formation of isoG-C and m5isoC-G base pairs with three hydrogen bonds. Molecular modeling shows that both structures form regular right handed helices with C2'endo sugar puckers. The role of the water content on the helical conformation of the parallel duplexes has been studied by FTIR and CD. It is found that a conformational transition similar to the B --> A transition observed for anti-parallel duplexes induced by a decrease of the water content of the samples can occur for these parallel duplexes. Their helical flexibility has been evidenced by FTIR studies on hydrated films by the emergence of absorption bands characteristic of A type geometry, in particular by an S-type --> N-type repuckering of the deoxyribose. All sugars in the parallel duplex with alternating d(isoG-A)/d(C-T) sequence can adopt an N-type geometry in low water content conditions. The conformational transition of the parallel hairpin duplex with alternating d(isoG-A)/d(C-T) sequence was followed by circular dichroism in water/trifluoroethanol solutions and its free energy at 0 degrees C was estimated to be 6.6 +/- 0.3 kcal mol(-1).     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

Formation of an amino-acid-binding pocket through adaptive zippering-up of a large DNA hairpin loop

Background: In vitro selection has identified DNA aptamers that target cofactors, amino acids, peptides and proteins. Structure determination of such ligand-DNA aptamer complexes should elucidate the details of adaptive DNA structural transitions, binding-pocket architectures and ligand recognition. We have determined the solution structure of the complex of a DNA aptamer containing a guanine-rich 18-residue hairpin loop that binds l-argininamide with ∼ 100μM affinity.Results: The DNA aptamer generates its l-argininamide-binding pocket by adaptive zippering up the 18-residue loop through formation of Watson-Crick pairs, mismatch pairs and base triples, while maximizing stacking interactions. Three of the four base triples involve minor-groove recognition through sheared G·A mismatch formation. The unique fold is also achieved through positioning of an adenine residue deep within the minor groove and through nestling of a smaller loop within the larger loop on complex formation. The accessibility to the unique l-argininamide-binding pocket is restricted by a base pair that bridges across one side of the major-groove-binding site. The guanidinium group of the bound l-argininamide aligns through intermolecular hydrogen-bond formation with the base edges of nonadjacent guanine and cytosine residues while being sandwiched between the planes of nonadjacent guanine residues.Conclusions: The available structures of l-arginine/l-argininamide bound to their DNA and RNA targets define the common principles and patterns associated with molecular recognition, as well as the diversity of intermolecular hydrogen-bonding alignments associated with the distinct binding pockets.     

Correlation effects on cis—trans energy differences in open-shell systems: The hoco radical as an example

Ab initio SCF and Cl calculations have shown that the most important factor in determining the relative energies of conformational isomers in free radicals is the movement of electron density into the singly occupied orbital. This is most favorable when an electron pair is situated trans to the radical site (trans correlation effect). In the HOCO radical we find that the trans isomer is more stable than the cis by 3.3 kcal/mol while SCF calculations yield virtually no energy difference between the two isomers.     

Solvent-assisted conformational isomerization and the conformationally-pure REMPI spectrum of 3-aminophenol

3-Aminophenol (3AP) has two conformers, cis and trans, depending on the orientation of the OH group relative to the NH(2) group. While both conformers are found in the jet-cooled spectra of 3AP, only the trans isomer was found in the REMPI spectrum of the 3AP(NH(3))(1) cluster. It was suggested that the cis conformer of the cluster isomerizes to the more stable trans conformer in the ground state during supersonic expansion. Solvent-assisted conformational isomerization (SACI) is believed to drive the population into the more stable trans isomer. SACI also occurs for the 3AP monomer, reducing 50% of the cis/trans ratio when the ammonia concentration in the expansion is higher than 0.1%. Depending on the expansion condition, the cis conformer can be completely depleted. When other solvents were introduced in the expansion, SACI occurred with only certain solvents whose binding energy is higher than the isomerization barrier. SACI can be used as a means to prepare the most stable conformer of gas phase biomolecules.