Simple method for determining nucleobase pK(a) values by indirect labeling and demonstration of a pK(a) of neutrality in dsDNA

Effect of ionic liquids as mobile phase additives on retention behaviors of G-quadruplexes in reversed-phase high performance liquid chromatography

G-quadruplexes (G4s) play an important role in a variety of biological processes and have extensive application prospects. Due to the significance of G4s in physiology and biosensing, studies on G4s have attracted much attention, stimulating the development or improvement of methods for G4 structures and polymorphism analysis. In this work, ionic liquids (ILs) were involved as mobile phase additives in reversed-phase high performance liquid chromatography (RP-HPLC) to analyse G4s with various conformations for the first time. How ILs affected the retention behaviors of G4s was investigated comprehensively. It was found that the addition of ILs markedly enhanced G4 retention, along with obvious amelioration on chromatographic peak shapes and separation. The influence of pH of mobile phase and types of ILs were also included in order to acquire an in-depth understanding. It appeared that the effect of ILs on G4 retention behaviors was the result of a combination of various interactions between G4s with the hydrophobic stationary phase and with the IL-containing mobile phase, where ion pair mechanism and enhanced hydrophobic interaction dominated. The findings of this work revealed that ILs could effectively improve the separation of G4s in RP-HPLC, which was conducive to G4 structural analysis, especially for G4s polymorphism elucidation.

Apurinic/apyrimidinic (AP) endonuclease 1 processing of AP sites with 5' mismatches

Despite the DNA duplex being central to biological functions, many intricacies of this molecule, including the dynamic nature of mismatched base pairing, are still unknown. The unique conformations adopted by DNA mismatches can provide insight into the forces at play between nucleotides. Moreover, DNA-binding proteins apply their own individualized steric and electrochemical influences on the nucleotides that they interact with, further altering base-pairing conformations. Here, seven X-ray crystallographic structures of the human nuclease apurinic/apyrimidinic (AP) endonuclease 1 (APE1) in complex with its substrate target flanked by a 5' mismatch are reported. The structures reveal how APE1 influences the conformations of a variety of different mismatched base pairs. Purine-purine mismatches containing a guanine are stabilized by a rotation of the guanine residue about the N-glycosidic bond to utilize the Hoogsteen edge for hydrogen bonding. Interestingly, no rotation of adenine, the other purine, is observed. Mismatches involving both purine and pyrimidine bases adopt wobble conformations to accommodate the mismatch. Pyrimidine-pyrimidine mismatches also wobble; however, the smaller profile of a pyrimidine base results in a gap between the Watson-Crick faces that is reduced by a C1'-C1' compression. These results advance our understanding of mismatched base pairing and the influence of a bound protein.

Phosphorothioate DNA Stabilized Fluorescent Gold and Silver Nanoclusters

Unmodified single-stranded DNA has recently gained popularity for the templated synthesis of fluorescent noble metal nanoclusters (NCs). Bright, stable, and biocompatible clusters have been developed primarily through optimization of DNA sequence. However, DNA backbone modifications have not yet been investigated. In this work, phosphorothioate (PS) DNAs are evaluated in the synthesis of Au and Ag nanoclusters, and are employed to successfully template a novel emitter using T₁₅ DNA at neutral pH. Mechanistic studies indicate a distinct UV-dependent formation mechanism that does not occur through the previously reported thymine N3. The positions of PS substitution have been optimized. This is the first reported use of a T₁₅ template at physiological pH for AgNCs.

Inverse thio effects in the hepatitis delta virus ribozyme reveal that the reaction pathway is controlled by metal ion charge density

The hepatitis delta virus (HDV) ribozyme self-cleaves in the presence of a wide range of monovalent and divalent ions. Prior theoretical studies provided evidence that self-cleavage proceeds via a concerted or stepwise pathway, with the outcome dictated by the valency of the metal ion. In the present study, we measure stereospecific thio effects at the nonbridging oxygens of the scissile phosphate under a wide range of experimental conditions, including varying concentrations of diverse monovalent and divalent ions, and combine these with quantum mechanical/molecular mechanical (QM/MM) free energy simulations on the stereospecific thio substrates. The R_P substrate gives large normal thio effects in the presence of all monovalent ions. The S_P substrate also gives normal or no thio effects, but only for smaller monovalent and divalent cations, such as Li⁺, Mg²⁺, Ca²⁺, and Sr²⁺; in contrast, sizable inverse thio effects are found for larger monovalent and divalent cations, including Na⁺, K⁺, NH₄⁺, and Ba²⁺. Proton inventories are found to be unity in the presence of the larger monovalent and divalent ions, but two in the presence of Mg²⁺. Additionally, rate–pH profiles are inverted for the low charge density ions, and only imidazole plus ammonium ions rescue an inactive C75Δ variant in the absence of Mg²⁺. Results from the thio effect experiments, rate–pH profiles, proton inventories, and ammonium/imidazole rescue experiments, combined with QM/MM free energy simulations, support a change in the mechanism of HDV ribozyme self-cleavage from concerted and metal ion-stabilized to stepwise and proton transfer-stabilized as the charge density of the metal ion decreases.

Experimental approaches for measuring pKa's in RNA and DNA

RNA and DNA carry out diverse functions in biology including catalysis, splicing, gene regulation, and storage of genetic information. Interest has grown in understanding how nucleic acids perform such sophisticated functions given their limited molecular repertoire. RNA can fold into diverse shapes that often perturb pKa values and allow it to ionize appreciably under biological conditions, thereby extending its molecular diversity. The goal of this chapter is to enable experimental measurement of pKa's in RNA and DNA. A number of experimental methods for measuring pKa values in RNA and DNA have been developed over the last 10 years, including RNA cleavage kinetics; UV-, fluorescence-, and NMR-detected pH titrations; and Raman crystallography. We begin with general considerations for choosing a pKa assay and then describe experimental conditions, advantages, and disadvantages for these assays. Potential pitfalls in measuring a pKa are provided including the presence of apparent pKa's due to a kinetic pKa or coupled acid- and alkali-promoted RNA unfolding, as well as degradation of RNA, precipitation of metal hydroxides and poor baselines. Use of multiple data fitting procedures and the study of appropriate mutants are described as ways to avoid some of these pitfalls. Application of these experimental methods to RNA and DNA will increase the number of available nucleic acid pKa values in the literature, which should deepen insight into biology and provide benchmarks for pKa calculations. Future directions for measuring pKa's in nucleic acids are discussed.

pKa shifting in double-stranded RNA is highly dependent upon nearest neighbors and bulge positioning

Shifting of pKa's in RNA is important for many biological processes; however, the driving forces responsible for shifting are not well understood. Herein, we determine how structural environments surrounding protonated bases affect pKa shifting in double-stranded RNA (dsRNA). Using (31)P NMR, we determined the pKa of the adenine in an A(+)·C base pair in various sequence and structural environments. We found a significant dependence of pKa on the base pairing strength of nearest neighbors and the location of a nearby bulge. Increasing nearest neighbor base pairing strength shifted the pKa of the adenine in an A(+)·C base pair higher by an additional 1.6 pKa units, from 6.5 to 8.1, which is well above neutrality. The addition of a bulge two base pairs away from a protonated A(+)·C base pair shifted the pKa by only ~0.5 units less than a perfectly base paired hairpin; however, positioning the bulge just one base pair away from the A(+)·C base pair prohibited formation of the protonated base pair as well as several flanking base pairs. Comparison of data collected at 25 °C and 100 mM KCl to biological temperature and Mg(2+) concentration revealed only slight pKa changes, suggesting that similar sequence contexts in biological systems have the potential to be protonated at biological pH. We present a general model to aid in the determination of the roles protonated bases may play in various dsRNA-mediated processes including ADAR editing, miRNA processing, programmed ribosomal frameshifting, and general acid-base catalysis in ribozymes.

Characterizing the protonation state of cytosine in transient G·C Hoogsteen base pairs in duplex DNA

G·C Hoogsteen base pairs can form transiently in duplex DNA and play important roles in DNA recognition, replication, and repair. G·C Hoogsteen base pairs are thought to be stabilized by protonation of cytosine N3, which affords a second key hydrogen bond, but experimental evidence for this is sparse because the proton cannot be directly visualized by X-ray crystallography and nuclear magnetic resonance spectroscopy. Here, we combine NMR and constant pH molecular dynamics simulations to directly investigate the pKa of cytosine N3 in a chemically trapped N1-methyl-G·C Hoogsteen base pair within duplex DNA. Analysis of NMR chemical shift perturbations and NOESY data as a function of pH revealed that cytosine deprotonation is coupled to a syn-to-anti transition in N1-methyl-G, which results in a distorted Watson-Crick geometry at pH >9. A four-state analysis of the pH titration profiles yields a lower bound pKa estimate of 7.2 ± 0.1 for the G·C Hoogsteen base pair, which is in good agreement with the pKa value (7.1 ± 0.1) calculated independently using constant pH MD simulations. Based on these results and pH-dependent NMR relaxation dispersion measurements, we estimate that under physiological pH (pH 7-8), G·C Hoogsteen base pairs in naked DNA have a population of 0.02-0.002%, as compared to 0.4% for A·T Hoogsteen base pairs, and likely exist primarily as protonated species.

Towards Accurate Prediction of Protonation Equilibrium of Nucleic Acids

The role of protonated nucleotides in modulating the pH-dependent properties of nucleic acids is one of the emerging frontiers in the field of nucleic acid biology. The recent development of a constant pH molecular dynamics simulation (CPHMD^MSλD) framework for simulating nucleic acids has provided a tool for realistic simulations of pH-dependent dynamics. We enhanced the CPHMD^MSλD framework with pH-based replica exchange (pH-REX), which significantly improves the sampling of both titration and spatial coordinates. The results from our pK_a calculations for the GAAA tetraloop, which was predicted with lower accuracy previously due to sampling challenges, demonstrates that pH-REX reduces the average unsigned error (AUE) to 0.7 pK_a units, and the error of the most poorly predicted residue A17 was drastically reduced from 2.9 to 1.2 pK_a unit. Lastly, we show that pH-REX CPHMD^MSλD simulations can be used to identify the dominant conformation of nucleic acid structures in alternate pH environments. This work suggests that pH-REX CPHMD^MSλD simulations provide a practical tool for predicting nucleic acid protonation equilibrium from first-principles, and offering structural and mechanistic insight into the study of pH-dependent properties of nucleic acids.

Charged nucleobases and their potential for RNA catalysis

Catalysis in living cells is carried out by both proteins and RNA. Protein enzymes have been known for over 200 years, but RNA enzymes, or "ribozymes", were discovered only 30 years ago. Developing insight into RNA enzyme mechanisms is invaluable for better understanding both extant biological catalysis as well as the primitive catalysis envisioned in an early RNA-catalyzed life. Natural ribozymes include large RNAs such as the group I and II introns; small RNAs such as the hepatitis delta virus and the hairpin, hammerhead, VS, and glmS ribozymes; and the RNA portion of the ribosome and spliceosome. RNA enzymes use many of the same catalytic strategies as protein enzymes, but do so with much simpler side chains. Among these strategies are metal ion, general acid-base, and electrostatic catalysis. In this Account, we examine evidence for participation of charged nucleobases in RNA catalysis. Our overall approach is to integrate direct measurements on catalytic RNAs with thermodynamic studies on oligonucleotide model systems. The charged amino acids make critical contributions to the mechanisms of nearly all protein enzymes. Ionized nucleobases should be critical for RNA catalysis as well. Indeed, charged nucleobases have been implicated in RNA catalysis as general acid-bases and oxyanion holes. We provide an overview of ribozyme studies involving nucleobase catalysis and the complications involved in developing these mechanisms. We also consider driving forces for perturbation of the pK(a) values of the bases. Mechanisms for pK(a) values shifting toward neutrality involve electrostatic stabilization and the addition of hydrogen bonding. Both mechanisms couple protonation with RNA folding, which we treat with a thermodynamic formalism and conceptual models. Furthermore, ribozyme reaction mechanisms can be multichannel, which demonstrates the versatility of ribozymes but makes analysis of experimental data challenging. We examine advances in measuring and analyzing perturbed pK(a) values in RNA. Raman crystallography and fluorescence spectroscopy have been especially important for pK(a) measurement. These methods reveal pK(a) values for the nucleobases A or C equal to or greater than neutrality, conferring potential histidine- and lysine/arginine-like behavior on them. Structural support for ionization of the nucleobases also exists: an analysis of RNA structures in the databases conducted herein suggests that charging of the bases is neither especially uncommon nor difficult to achieve under cellular conditions. Our major conclusions are that cationic and anionic charge states of the nucleobases occur in RNA enzymes and that these states make important catalytic contributions to ribozyme activity. We conclude by considering outstanding questions and possible experimental and theoretical approaches for further advances.

Phase behavior of crowded like-charged mixed polyelectrolytes in a cell-sized sphere

We studied the phase behavior of a mixture of two semiflexible negatively charged polyelectrolytes, giant DNA and alginate, under crowded condition in a cell-sized sphere (5-40 μm in diameter), where the persistence length of giant DNA is 50 nm and that of alginate is 5 nm. Through microscopic observation, we found that the polymer mixture exhibits a unique phase behavior, which depends on the size of the sphere, whereas the mixture remains homogeneous and isotropic in bulk solution. When the sphere is small, DNA is completely depleted on the surface. When the sphere is medium sized, a portion of the DNA is depleted on the surface, and the remainder stays within the sphere. By introducing a curvature-dependent term for the interaction between DNA and the surface into the Flory-Huggins model, we interpret the observed characteristics of the phase behavior in terms of the relative importance of the surface-to-volume effect.

Driving forces for nucleic acid pK(a) shifting in an A(+).C wobble: effects of helix position, temperature, and ionic strength

Secondary structure plays critical roles in nucleic acid function. Mismatches in DNA can lead to mutation and disease, and some mismatches involve a protonated base. Among protonated mismatches, A(+).C wobble pairs form near physiological pH and have relatively minor effects on helix geometry, making them especially important in biology. Herein, we investigate effects of helix position, temperature, and ionic strength on pK(a) shifting in A(+).C wobble pairs in DNA. We observe that pK(a) shifting is favored by internal A(+).C wobbles, which have low cooperativities of folding and make large contributions to stability, and disfavored by external A(+).C wobbles, which have high folding cooperativities but make small contributions to stability. An inverse relationship between pK(a) shifting and temperature is also found, which supports a model in which protonation is enthalpically favored overall and entropically correlated with cooperativity of folding. We also observe greater pK(a) shifts as the ionic strength decreases, consistent with anticooperativity between proton binding and counterion-condensed monovalent cation. Under the most favorable temperature and ionic strength conditions tested, a pK(a) of 8.0 is observed for the A(+).C wobble pair, which represents an especially large shift ( approximately 4.5 pK(a) units) from the unperturbed pK(a) value of adenosine. This study suggests that protonated A(+).C wobble pairs exist in DNA under biologically relevant conditions, where they can drive conformational changes and affect replication and transcription fidelity.

Direct Raman measurement of an elevated base pKa in the active site of a small ribozyme in a precatalytic conformation

Catalytic RNA molecules can achieve rate acceleration by shifting base pK(a) values toward neutrality. Prior evidence has suggested that base A38 of the hairpin ribozyme plays an important role in phosphoryl transfer, possibly functioning as a general acid, or by orienting a specific water molecule for proton transfer. To address the role of A38, we used Raman spectroscopy to measure directly the pK(a) of the N1-imino moiety in the context of hairpin ribozyme crystals representative of a "precatalytic" conformation. The results revealed that the pK(a) of A38 is shifted to 5.46 +/- 0.05 relative to 3.68 +/- 0.06 derived from a reference solution of the nucleotide AMP. The elevated pK(a) correlates well with the first titration point of the macroscopic pH-rate profile of the hairpin ribozyme in solution and strongly supports A38 as a general acid catalyst in bond scission. The results confirm that A38 is protonated before the transition state, which would promote phosphorane development. Overall, the results establish a cogent structure-function paradigm that expands our understanding of how RNA structure can enhance nucleobase reactivity to catalyze biological reactions.

Base pairing configuration and stability of an oligonucleotide duplex containing a 5-chlorouracil-adenine base pair

Inflammation-mediated reactive molecules can damage DNA by oxidation and chlorination. The biological consequences of this damage are as yet incompletely understood. In this paper, we have constructed oligonucleotides containing 5-chlorouracil (ClU), one of the known inflammation damage products. The thermodynamic stability, base pairing configuration, and duplex conformation of oligonucleotides containing ClU paired opposite adenine have been examined. NMR spectra reveal that the ClU-A base pair adopts a geometry similar to that of the T-A base pair, and the ClU-A base pair-containing duplex adopts a normal B-form conformation. The line width of the imino proton of the ClU residue is substantially greater than that of the corresponding T imino proton; however, this difference is not attributed to a reduced thermal or thermodynamic stability or to an increased level of proton exchange with solvent. While the NMR studies reveal an increased level of chemical exchange for the ClU imino proton of the ClU-A base pair, the ClU residue is not a target for removal by the Escherichia coli mispaired uracil glycosylase, which senses damage-related helix instability. The results of this study are consistent with previous reports indicating that the DNA of replicating cells can tolerate substantial substitution with ClU. The fraudulent, pseudo-Watson-Crick ClU-A base pair is sufficiently stable to avoid glycosylase removal and, therefore, might constitute a persistent form of cellular DNA damage.

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNA^Phe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules.

Fluorine substituted adenosines as probes of nucleobase protonation in functional RNAs

Ionized nucleobases are required for folding, conformational switching, or catalysis in a number of functional RNAs. A common strategy to study these sites employs nucleoside analogues with perturbed pKa, but the interpretation of these studies is often complicated by the chemical modification introduced, in particular modifications that add, remove, or translocate hydrogen bonding groups in addition to perturbing pKa values. In the present study we present a series of fluorine substituted adenosine analogues that produce large changes in N1 pKa values with minimal structural perturbation. These analogues include fluorine for hydrogen substitutions in the adenine ring of adenosine and 7-deaza-adenosine with resulting N1 pKa values spanning more than 4 pKa units. To demonstrate the utility of these analogues we have conducted a nucleotide analogue interference mapping (NAIM) study on a self-ligating construct of the Varkud Satellite (VS) ribozyme. We find that each of the analogues is readily incorporated by T7 RNA polymerase and produces fully active transcripts when substituted at the majority of sites. Strong interferences are observed for three sites known to be critical for VS ribozyme function, most notably A756. Substitutions at A756 lead to slight enhancements in activity for elevated pKa analogues and dramatic interferences in activity for reduced pKa analogues, supporting the proposed catalytic role for this base. The structural similarity of these analogues, combined with their even incorporation and selective interference, provides an improved method for identifying sites of adenosine protonation in a variety of systems.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Active Site Labeling of G8 in the Hairpin Ribozyme: Implications for Structure and Mechanism

There is mounting evidence that suggests that general acid/base catalysis is operative in the hairpin ribozyme, with analogy to the protein enzyme RNaseA. Nevertheless, the extent of general base catalysis as well as the identity of the specific chemical groups responsible remains the subject of some controversy. An affinity label has previously been used to alkylate histidine 12 (His12), the active general base in RNaseA. To date, no such experiment has been applied to a ribozyme. We have synthesized the analogous affinity label for the hairpin ribozyme with an electrophilic 2'-bromoacetamide group in lieu of the 2'-hydroxyl (2'OH) at the substrate cleavage site and show that guanosine 8 (G8) of the hairpin ribozyme is specifically alkylated, most likely at the N1 position. This evidence strongly implicates N1 of G8 in active site chemistry. By direct analogy to RNase A, these findings could be consistent with the hypothesis that deprotonated G8 residue functions as a general base in the hairpin ribozyme. Other mechanistic possibilities for N1 of G8 such as indirect general base catalysis mediated by a water molecule or transition state stabilization could also be consistent with our findings.

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.

Computational and empirical trans-hydrogen bond deuterium isotope shifts suggest that N1-N3 A:U hydrogen bonds of RNA are shorter than those of A:T hydrogen bonds of DNA

Density functional theory calculations of isolated Watson-Crick A:U and A:T base pairs predict that adenine 13C2 trans-hydrogen bond deuterium isotope shifts due to isotopic substitution at the pyrimidine H3, (2h)Delta13C2, are sensitive to the hydrogen-bond distance between the N1 of adenine and the N3 of uracil or thymine, which supports the notion that (2h)Delta13C2 is sensitive to hydrogen-bond strength. Calculated (2h)Delta13C2 values at a given N1-N3 distance are the same for isolated A:U and A:T base pairs. Replacing uridine residues in RNA with 5-methyl uridine and substituting deoxythymidines in DNA with deoxyuridines do not statistically shift empirical (2h)Delta13C2 values. Thus, we show experimentally and computationally that the C7 methyl group of thymine has no measurable affect on (2h)Delta13C2 values. Furthermore, (2h)Delta13C2 values of modified and unmodified RNA are more negative than those of modified and unmodified DNA, which supports our hypothesis that RNA hydrogen bonds are stronger than those of DNA. It is also shown here that (2h)Delta13C2 is context dependent and that this dependence is similar for RNA and DNA.

Linkage between proton binding and folding in RNA: implications for RNA catalysis

Small ribozymes use their nucleobases to catalyse phosphodiester bond cleavage. The hepatitis delta virus ribozyme employs C75 as a general acid to protonate the 5′-bridging oxygen leaving group, and to accomplish this task efficiently, it shifts its pKa towards neutrality. Simulations and thermodynamic experiments implicate linkage between folding and protonation in nucleobase pKa shifting. Even small oligonucleotides are shown to fold in a highly co-operative manner, although they do so in a context-specific fashion. Linkage between protonation and co-operativity of folding may drive pKa shifting and provide for enhanced function in RNA.

Linkage between proton binding and folding in RNA: a thermodynamic framework and its experimental application for investigating pKa shifting

Perturbation of pKa values can change the favored protonation states of the nucleobases at biological pH and thereby modulate the function of RNA and DNA molecules. In an effort to understand the driving forces for pKa shifting specific to nucleic acids, we developed a thermodynamic framework that relates proton binding to the nucleobases and the helix-coil transition. Key features that emerge from the treatment are a comprehensive description of all the actions of proton binding on RNA folding: acid and alkaline denaturation of the helix and pKa shifting in the folded state. Practical experimental approaches for measuring pKas from thermal denaturation experiments are developed. Microscopic pka values (where ka is the acid dissociation constant) for the unfolded state were determined directly by experiments on unstructured oligonucleotides, which led to a macroscopic pKa for the ensemble of unfolded states shifted toward neutrality. The formalism was then applied to pH-dependent UV melting data for model DNA oligonucleotides. Folded-state pka) values were in good agreement with the outcome of pH titrations, and the acid and alkaline denaturation regions were well described. The formalism developed here is similar to that of Draper and coworkers for Mg2+ binding to RNA, except that the unfolded state is described explicitly owing to the presence of specific proton-binding sites on the bases. A principal conclusion is that it should be possible to attain large pKa shifts by designing RNA molecules that fold cooperatively.