Premelting and the hydrogen-exchange open state in synthetic RNA duplexes

A General Method for Predicting Amino Acid Residues Experiencing Hydrogen Exchange

Information on protein hydrogen exchange can help delineate key regions involved in protein-protein interactions and provides important insight towards determining functional roles of genetic variants and their possible mechanisms in disease processes. Previous studies have shown that the degree of hydrogen exchange is affected by hydrogen bond formations, solvent accessibility, proximity to other residues, and experimental conditions. However, a general predictive method for identifying residues capable of hydrogen exchange transferable to a broad set of proteins is lacking. We have developed a machine learning method based on random forest that can predict whether a residue experiences hydrogen exchange. Using data from the Start2Fold database, which contains information on 13,306 residues (3,790 of which experience hydrogen exchange and 9,516 which do not exchange), our method achieves good performance. Specifically, we achieve an overall out-of-bag (OOB) error, an unbiased estimate of the test set error, of 20.3 percent. Using a randomly selected test data set consisting of 500 residues experiencing hydrogen exchange and 500 which do not, our method achieves an accuracy of 0.79, a recall of 0.74, a precision of 0.82, and an F1 score of 0.78.

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

RNA structure analysis at single nucleotide resolution by selective 2'-hydroxyl acylation and primer extension (SHAPE)

The reactivity of an RNA ribose hydroxyl is shown to be exquisitely sensitive to local nucleotide flexibility because a conformationally constrained adjacent 3'-phosphodiester inhibits formation of the deprotonated, nucleophilic oxyanion form of the 2'-hydroxyl group. Reaction with an appropriate electrophile, N-methylisatoic anhydride, to form a 2'-O-adduct thus can be used to monitor local structure at every nucleotide in an RNA. We develop a quantitative approach involving Selective 2'-Hydroxyl Acylation analyzed by Primer Extension (SHAPE) to map the structure of and to distinguish fine differences in structure for tRNAAsp transcripts at single nucleotide resolution. Modest extensions of the SHAPE approach will allow RNA structure to be monitored comprehensively and at single nucleotide resolution for RNAs of arbitrary sequence and structural complexity and under diverse solution environments.

The opening of a single base without perturbations of neighboring nucleotides: a study on crystal B-DNA duplex d(CGCGAATTCGCG)2

In this work we explore the possibility of the opening of a single base without perturbation of its neighboring nucleotides. Low energy base opening into the grooves can be accomplished by rotation of the relevant backbone and glycosidic bond torsion angles. The pathway has been determined by identifying zeta torsion angle as the reaction coordinate together with the accompanying geometric requirement that guides the displacement of other torsion angles. Our study on Dickerson dodecamer duplex d(CGCGAATTCGCG)2 showed that all bases with normal equilibrium zeta can be rotated by approximately 30 degrees, corresponding to approximately 3.5A base displacement, towards the major groove. Such an opening extent is comparable with estimated amplitudes of local angular motions in DNA bases from NMR experiments, which might facilitate proton exchange. The computed base opening energy barrier is also comparable with measured base pair opening enthalpy. These results indicate possible relevance of the pathway studied in this work with experimentally observed base pair opening process. Our analysis also showed a preference for base opening along the major groove and an abnormal opening behavior for bases with unusual equilibrium zeta torsion angle.

Synergistic effects in the melting of DNA hydration shell: melting of the minor groove hydration spine in poly(dA).poly(dT) and its effect on base pair stability

We propose that water of hydration in contact with the double helix can exist in several states. One state, found in the narrow groove of poly(dA).poly(dT), should be considered as frozen to the helix, i.e., an integral part of the double helix. We find that this enhanced helix greatly effects the stability of that helix against base separation melting. Most water surrounding the helix is, however, melted or disassociated with respect to being an integral part of helix and plays a much less significant role in stabilizing the helix dynamically, although these water molecules play an important role in stabilizing the helix conformation statically. We study the temperature dependence of the melting of the hydration spine and find that narrow groove nonbonded interactions are necessary to stabilize the spine above room temperature and to show the broad transition observed experimentally. This calculation requires that synergistic effects of nonbonded interactions between DNA and its hydration shell affect the state of water-base atom hydrogen bonds. The attraction of waters into narrow groove tends to retain waters in the groove and compress or strain these hydrogen bonds.

The role of a minor groove spine of hydration in stabilizing poly(dA).poly(dT) against fluctuational interbase H-bond disruption in the premelting temperature regime

Experimental estimates of the premelting Adenine-Thymine base pair opening probability for some B-DNA sequences are two orders of magnitude smaller than those of other B-DNA sequences. The AT pairs in the sequence with smaller open probability seem to be those that have a well defined spine of hydration in the minor groove. We show that this spine of hydration can significantly enhance the thermal stability of the base pairs to which they are attached. The effect of this spine of hydration coupled with the possible stabilization effect contributed from neighboring GC pairs can explain the differences in the observed AT pair opening probability for different AT containing B-DNA sequences.

Premelting thermal fluctuational interbase hydrogen-bond disrupted states of a B-DNA guanine-cytosine base pair: significance for amino and imino proton exchange

Modified self-consistent phonon theory when applied to the DNA double helix indicates the existence of fairly long-lived states in which single interbase H bonds are disrupted. One can then postulate a number of situations in which particular disrupted H bonds can enhance particular proton exchange. In this paper we postulate a number of such partially open states for a B-conformation GC base pair and calculate the probability of each of these states for a B-conformation poly(dG).poly(dC). We compare these probabilities to those probabilities needed to explain various observed proton exchange rates. We propose that, for a GC base pair in B conformation, there are two amino proton exchangeable states--a cytosine amino proton exchangeable state and a guanine amino proton exchangeable state; both require the disruption of only the corresponding interbase H bond. The imino proton exchange, however, requires the disruption of all three interbase H bonds and this defines a third open state. Our calculated probabilities for a GC base pair in these three states are in fair agreement with available experimental estimates from measurements of amino and imino proton exchange.

Premelting thermal fluctuational base pair opening probability of poly(dA).poly(dT) as predicted by the modified self-consistent phonon theory

We employ a mean field, modified, self-consistent phonon theory to evaluate the single base-pair opening rate and the probability of a base pair in the amino proton exchangeable state for the homopolymer poly(dA).poly(dT) at temperatures below the helix-coil transition region. Our calculated premelting single base-pair opening probabilities are in general agreement with several available experimental estimates from imino proton exchange and formaldehyde-induced DNA melting measurements. These calculated opening probabilities, however, are in disagreement with the prediction of the helix-coil transition theory. Possible reasons for the differences are discussed, especially the possible different definition of a meaningful open state in the premelting region. The premelting open state of the modified self-consistent phonon approximation theory seems to be appropriate to describe a solvent-accessible open configuration that is sufficient to facilitate important chemical reactions such as imino proton exchange and formaldehyde reaction with the bases. This can be compared with the completely unstacked open state of the helix-coil transition theory originally defined in the helix-coil transition region. We propose that the amino proton exchangeable state is different from the open state associated with melting and only involves the breaking of the amino interbase H bond. The agreement between the calculated and experimentally estimated probability of a base pair in the amino proton exchangeable state seems to support this hypothesis.

Interaction of chloroquine with linear and supercoiled DNAs. Effect on the torsional dynamics, rigidity, and twist energy parameter

The magnitude and uniformity of the torsion elastic constant (alpha) of linear pBR322 DNA and supercoiled pBR322 DNAs with high-twist (sigma = -0.083) and normal-twist (sigma = -0.48) are measured in 0.1 M NaCl as a function of added chloroquine/base-pair ratio (chl/bp) by studying the fluorescence polarization anisotrophy (FPA) of intercalated ethidium dye. The time-resolved FPA is measured by using a picosecond dye laser for excitation and time-correlated single-photon counting detection. A general theory is developed for the binding of ligands that unwind superhelical DNAs, and the simultaneous binding of two different intercalators is treated in detail. The equilibrium constant (K) for binding chloroquine to linear pBR322 DNA and the number (r) of bound chloroquines per base pair are determined from the relative amplitude ratio of the slow (normally intercalated) and fast (free) components in the decay of the (probe) ethidium fluorescence intensity as a function of chl/bp. For chloroquine binding to supercoiled pBR322 DNAs, the intrinsic binding constant is assumed to be the same as for the linear DNA, but the twist energy parameter ET (N times the free energy to change the linking number from 0 to 1 in units of kBT) is regarded as adjustable. Using the best-fit ET, the binding ratios r are calculated for each chl/bp ratio. Twist energy parameters are also determined for ethidium binding to these supercoiled DNAs by competitive dialysis. For chloroquine binding, we obtain ET = 360 and 460 respectively for the normal-twist and high-twist supercoiled DNAs. For ethidium binding the corresponding values are ET = 280 +/- 70 and 347 +/- 50. Like other dye-binding values, these are substantially lower than those obtained by ligation methods. In the absence of chloroquine, the torsion constants of all three DNAs are virtually identical, alpha = (5.0 +/- 0.4) x 10(-12) dyn.cm. For linear pBR322 DNA, the magnitude and uniformity of alpha remain unaltered by intercalated chloroquine up to r = 0.19. This finding argues that the FPA is not significantly relaxed by diffusion of any kinks or solitons. If alpha d denotes the torsion constant between a dye and a base pair and alpha 0 that between two base pairs, then our data imply that alpha d/alpha 0 lies in the range 0.65-1.64, with a most probable value of 1.0.(ABSTRACT TRUNCATED AT 400 WORDS)