Interactions of the osmolyte glycine betaine with molecular surfaces in water: thermodynamics, structural interpretation, and prediction of m-values

Its preferential interactions with biopolymers account for diverse observed effects of trehalose

Biopolymer homeostasis underlies the health of organisms, and protective osmolytes have emerged as one strategy used by Nature to preserve biopolymer homeostasis. However, a great deal remains unknown about the mechanism of action of osmolytes. Trehalose, as a prominent example, stabilizes proteins against denaturation by extreme temperature and denaturants, preserves membrane integrity upon freezing or in dry conditions, inhibits polyQ-mediated protein aggregation, and suppresses the aggregation of denatured proteins. The underlying thermodynamic mechanisms of such diverse effects of trehalose remain unclear or controversial. In this study, we applied the surface-additive method developed in the Record laboratory to attack this issue. We characterized the key features of trehalose-biopolymer preferential interactions and found that trehalose has strong unfavorable interactions with aliphatic carbon and significant favorable interactions with amide/anionic oxygen. This dissection has allowed us to elucidate the diverse effects of trehalose and to identify the crucial functional group(s) responsible for its effects. With (semi)quantitative thermodynamic analysis, we discovered that 1) the unfavorable interaction of trehalose with hydrophobic surfaces is the dominant factor in its effect on protein stability, 2) the favorable interaction of trehalose with polar amides enables it to inhibit polyQ-mediated protein aggregation and the aggregation of denatured protein in general, and 3) the favorable interaction of trehalose with phosphate oxygens, together with its unfavorable interaction with aliphatic carbons, enables trehalose to preserve membrane integrity in aqueous solution. These results provide a basis for a full understanding of the role of trehalose in biopolymer homeostasis and the reason behind its evolutionary selection as an osmolyte, as well as for a better application of trehalose as a chemical chaperone.

Chemical Interactions of Polyethylene Glycols (PEGs) and Glycerol with Protein Functional Groups: Applications to Effects of PEG and Glycerol on Protein Processes

In this work, we obtain the data needed to predict chemical interactions of polyethylene glycols (PEGs) and glycerol with proteins and related organic compounds and thereby interpret or predict chemical effects of PEGs on protein processes. To accomplish this, we determine interactions of glycerol and tetraEG with >30 model compounds displaying the major C, N, and O functional groups of proteins. Analysis of these data yields coefficients (α values) that quantify interactions of glycerol, tetraEG, and PEG end (-CH2OH) and interior (-CH2OCH2-) groups with these groups, relative to interactions with water. TetraEG (strongly) and glycerol (weakly) interact favorably with aromatic C, amide N, and cationic N, but unfavorably with amide O, carboxylate O, and salt ions. Strongly unfavorable O and salt anion interactions help make both small and large PEGs effective protein precipitants. Interactions of tetraEG and PEG interior groups with aliphatic C are quite favorable, while interactions of glycerol and PEG end groups with aliphatic C are not. Hence, tetraEG and PEG300 favor unfolding of the DNA-binding domain of lac repressor (lacDBD), while glycerol and di- and monoethylene glycol are stabilizers. Favorable interactions with aromatic and aliphatic C explain why PEG400 greatly increases the solubility of aromatic hydrocarbons and steroids. PEG400-steroid interactions are unusually favorable, presumably because of simultaneous interactions of multiple PEG interior groups with the fused ring system of the steroid. Using α values reported here, chemical contributions to PEG m-values can be predicted or interpreted in terms of changes in water-accessible surface area (ΔASA) and separated from excluded volume effects.

Thermodynamics and solvent linkage of macromolecule-ligand interactions

Binding involves two steps, desolvation and association. While water is ubiquitous and occurs at high concentration, it is typically ignored. In vitro experiments typically use infinite dilution conditions, while in vivo, the concentration of water is decreased due to the presence of high concentrations of molecules in the cellular milieu. This review discusses isothermal titration calorimetry approaches that address the role of water in binding. For example, use of D2O allows the contribution of solvent reorganization to the enthalpy component to be assessed. Further, the addition of osmolytes will decrease the water activity of a solution and allow effects on Ka to be determined. In most cases, binding becomes tighter in the presence of osmolytes as the desolvation penalty associated with binding is minimized. In other cases, the osmolytes prefer to interact with the ligand or protein, and if their removal is more difficult than shedding water, then binding can be weakened. These complicating layers can be discerned by different slopes in ln(Ka) vs osmolality plots and by differential scanning calorimetry in the presence of the osmolyte.

An osmolyte mitigates the destabilizing effect of protein crowding

Most theories predict that macromolecular crowding stabilizes globular proteins, but recent studies show that weak attractive interactions can result in crowding-induced destabilization. Osmolytes are ubiquitous in biology and help protect cells against stress. Given that dehydration stress adds to the crowded nature of the cytoplasm, we speculated that cells might use osmolytes to overcome the destabilization caused by the increased weak interactions that accompany desiccation. We used NMR-detected amide proton exchange experiments to measure the stability of the test protein chymotrypsin inhibitor 2 under physiologically relevant crowded conditions in the presence and absence of the osmolyte glycine betaine. The osmolyte overcame the destabilizing effect of the cytosol. This result provides a physiologically relevant explanation for the accumulation of osmolytes by dehydration-stressed cells.

Investigation of osmolyte effects on FolM: comparison with other dihydrofolate reductases

A weak association between osmolytes and dihydrofolate (DHF) decreases the affinity of the substrate for the Escherichia coli chromosomal and R67 plasmid dihydrofolate reductase (DHFR) enzymes. To test whether the osmolyte-DHF association also interferes with binding of DHF to FolM, an E. coli enzyme that possesses weak DHFR activity, ligand binding was monitored in the presence of osmolytes. The affinity of FolM for DHF, measured by kcat/Km(DHF), was decreased by the addition of an osmolyte. Additionally, binding of the antifolate drug, methotrexate, to FolM was weakened by the addition of an osmolyte. The changes in ligand binding with water activity were unique for each osmolyte, indicating preferential interaction between the osmolyte and folate and its derivatives; however, additional evidence provided support for further interactions between FolM and osmolytes. Binding of the reduced nicotinamide adenine dinucleotide phosphate (NADPH) cofactor to FolM was monitored by isothermal titration calorimetry as a control for protein-osmolyte association. In the presence of betaine (proposed to be the osmolyte most excluded from protein surfaces), the NADPH Kd decreased, consistent with dehydration effects. However, other osmolytes did not tighten binding to the cofactor. Rather, dimethyl sulfoxide (DMSO) had no effect on the NADPH Kd, while ethylene glycol and polyethylene glycol 400 weakened cofactor binding. Differential scanning calorimetry of FolM in the presence of osmolytes showed that both DMSO and ethylene glycol decreased the stability of FolM, while betaine increased the stability of the protein. These results suggest that some osmolytes can destabilize FolM by preferentially interacting with the protein. Further, these weak attractions can impede ligand binding. These various contributions have to be considered when interpreting osmotic pressure results.

Quantifying the temperature dependence of glycine-betaine RNA duplex destabilization

Glycine-betaine (GB) stabilizes folded protein structure because of its unfavorable thermodynamic interactions with amide oxygen and aliphatic carbon surface area exposed during protein unfolding. However, GB can attenuate nucleic acid secondary structure stability, although its mechanism of destabilization is not currently understood. Here we quantify GB interactions with the surface area exposed during thermal denaturation of nine RNA dodecamer duplexes with guanine-cytosine (GC) contents of 17-100%. Hyperchromicity values indicate increasing GB molality attenuates stacking. GB destabilizes higher-GC-content RNA duplexes to a greater extent than it does low-GC-content duplexes due to greater accumulation at the surface area exposed during unfolding. The accumulation is very sensitive to temperature and displays characteristic entropy-enthalpy compensation. Since the entropic contribution to the m-value (used to quantify GB interaction with the RNA solvent-accessible surface area exposed during denaturation) is more dependent on temperature than is the enthalpic contribution, higher-GC-content duplexes with their larger transition temperatures are destabilized to a greater extent than low-GC-content duplexes. The concentration of GB at the RNA surface area exposed during unfolding relative to bulk was quantified using the solute-partitioning model. Temperature correction predicts a GB concentration at 25 °C to be nearly independent of GC content, indicating that GB destabilizes all sequences equally at this temperature.

Liberate and grab it, ingest and digest it: the GbdR regulon of the pathogen Pseudomonas aeruginosa

The compatible solute glycine betaine is a powerful osmostress protectant, but many microorganisms can also use it as a nutrient. K. J. Hampel et al. (J. Bacteriol. 196:7-15, 2014) defined a regulon in the notorious pathogen Pseudomonas aeruginosa that comprises modules for the harvest and import of the glycine betaine biosynthetic precursor choline and its subsequent catabolism to pyruvate. The reported data link the GbdR activator with the metabolism of host-derived compounds (e.g., phosphocholine) and virulence traits of P. aeruginosa.

Probing the protein-folding mechanism using denaturant and temperature effects on rate constants

Protein folding has been extensively studied, but many questions remain regarding the mechanism. Characterizing early unstable intermediates and the high-free-energy transition state (TS) will help answer some of these. Here, we use effects of denaturants (urea, guanidinium chloride) and temperature on folding and unfolding rate constants and the overall equilibrium constant as probes of surface area changes in protein folding. We interpret denaturant kinetic m-values and activation heat capacity changes for 13 proteins to determine amounts of hydrocarbon and amide surface buried in folding to and from TS, and for complete folding. Predicted accessible surface area changes for complete folding agree in most cases with structurally determined values. We find that TS is advanced (50-90% of overall surface burial) and that the surface buried is disproportionately amide, demonstrating extensive formation of secondary structure in early intermediates. Models of possible pre-TS intermediates with all elements of the native secondary structure, created for several of these proteins, bury less amide and hydrocarbon surface than predicted for TS. Therefore, we propose that TS generally has both the native secondary structure and sufficient organization of other regions of the backbone to nucleate subsequent (post-TS) formation of tertiary interactions. The approach developed here provides proof of concept for the use of denaturants and other solutes as probes of amount and composition of the surface buried in coupled folding and other large conformational changes in TS and intermediates in protein processes.

Quantifying additive interactions of the osmolyte proline with individual functional groups of proteins: comparisons with urea and glycine betaine, interpretation of m-values

To quantify interactions of the osmolyte l-proline with protein functional groups and predict their effects on protein processes, we use vapor pressure osmometry to determine chemical potential derivatives dμ2/dm3 = μ23, quantifying the preferential interactions of proline (component 3) with 21 solutes (component 2) selected to display different combinations of aliphatic or aromatic C, amide, carboxylate, phosphate or hydroxyl O, and amide or cationic N surface. Solubility data yield μ23 values for four less-soluble solutes. Values of μ23 are dissected using an ASA-based analysis to test the hypothesis of additivity and obtain α-values (proline interaction potentials) for these eight surface types and three inorganic ions. Values of μ23 predicted from these α-values agree with the experiment, demonstrating additivity. Molecular interpretation of α-values using the solute partitioning model yields partition coefficients (Kp) quantifying the local accumulation or exclusion of proline in the hydration water of each functional group. Interactions of proline with native protein surfaces and effects of proline on protein unfolding are predicted from α-values and ASA information and compared with experimental data, with results for glycine betaine and urea, and with predictions from transfer free energy analysis. We conclude that proline stabilizes proteins because of its unfavorable interactions with (exclusion from) amide oxygens and aliphatic hydrocarbon surfaces exposed in unfolding and that proline is an effective in vivo osmolyte because of the osmolality increase resulting from its unfavorable interactions with anionic (carboxylate and phosphate) and amide oxygens and aliphatic hydrocarbon groups on the surface of cytoplasmic proteins and nucleic acids.

Diversity in the mechanisms of cosolute action on biomolecular processes

Numerous cellular cosolutes significantly impact the way that proteins and other biomacromolecules act and interact. We have followed the thermodynamic effect of several cosolute classes, including polymers, cellular osmolytes, and inorganic salts, on the stability of biomolecular folding and complexation. By comparing changes in free energy, enthalpy, and entropy upon cosolutes addition for these processes, we identify several thermodynamically distinct mechanisms. Surprisingly, even while many cosolutes display similar scaling of the change in stabilizing free energy with their concentration, a breakdown of this free energy into enthalpic and entropic contributions distinguishes different families of cosolutes. We discuss how these "thermodynamic fingerprints" can direct towards possible underlying mechanisms that govern the cosolute effect.

Quantifying functional group interactions that determine urea effects on nucleic acid helix formation

Urea destabilizes helical and folded conformations of nucleic acids and proteins, as well as protein-nucleic acid complexes. To understand these effects, extend previous characterizations of interactions of urea with protein functional groups, and thereby develop urea as a probe of conformational changes in protein and nucleic acid processes, we obtain chemical potential derivatives (μ23 = dμ2/dm3) quantifying interactions of urea (component 3) with nucleic acid bases, base analogues, nucleosides, and nucleotide monophosphates (component 2) using osmometry and hexanol-water distribution assays. Dissection of these μ23 values yields interaction potentials quantifying interactions of urea with unit surface areas of nucleic acid functional groups (heterocyclic aromatic ring, ring methyl, carbonyl and phosphate O, amino N, sugar (C and O); urea interacts favorably with all these groups, relative to interactions with water. Interactions of urea with heterocyclic aromatic rings and attached methyl groups (as on thymine) are particularly favorable, as previously observed for urea-homocyclic aromatic ring interactions. Urea m-values determined for double helix formation by DNA dodecamers near 25 °C are in the range of 0.72-0.85 kcal mol(-1)m(-1) and exhibit little systematic dependence on nucleobase composition (17-42% GC). Interpretation of these results using the urea interaction potentials indicates that extensive (60-90%) stacking of nucleobases in the separated strands in the transition region is required to explain the m-value. Results for RNA and DNA dodecamers obtained at higher temperatures, and literature data, are consistent with this conclusion. This demonstrates the utility of urea as a quantitative probe of changes in surface area (ΔASA) in nucleic acid processes.

Estimating contribution of anthocyanin pigments to osmotic adjustment during winter leaf reddening

The association between plant water stress and synthesis of red, anthocyanin pigments in leaves has led some plant biologists to propose an osmotic function of leaf reddening. According to this hypothesis, anthocyanins function as a solute in osmotic adjustment (OA), contributing to depression of osmotic potential (Ψ(π)) and maintenance of turgor pressure during drought-stressed conditions. Here we calculate the percent contribution of anthocyanin to leaf Ψ(π) during OA in two angiosperm evergreen species, Galax urceolata and Gaultheria procumbens. Both species exhibit dramatic leaf reddening under high light during winter, concomitant with declines in leaf water potential and accumulation of solutes. Data previously published by the authors on osmotic potential at full turgor (Ψ(π,100)) of G. urceolata and G. procumbens leaves before and after leaf reddening were used to estimate OA. In vivo molar concentrations of anthocyanin, glucose, fructose, and sucrose measured from the same individuals were converted to pressure equivalents using the Ideal Gas Law, and percent contribution to OA was estimated. Estimated mean OA during winter was -0.7MPa for G. urceolata and -0.8MPa for G. procumbens. In vivo concentrations of anthocyanin (3-10mM) were estimated to account for ∼2% of OA during winter, and comprised <0.7% of Ψ(π,100) in both species. Glucose, fructose, and sucrose combined accounted for roughly 50 and 80% of OA for G. urceolata and G. procumbens, respectively, and comprised ∼20% of Ψ(π,100). We observed that a co-occurring, acyanic species (Vinca minor) achieved similar OA without synthesizing anthocyanin. We conclude that anthocyanins represent a measurable, albeit meager, component of OA in red-leafed evergreen species during winter. However, due to their low concentrations, metabolic costliness relative to other osmolytes, and striking red color (unnecessary for an osmotic function), it is unlikely that they are synthesized solely for an osmoprotectant role.

Introductory lecture: interpreting and predicting Hofmeister salt ion and solute effects on biopolymer and model processes using the solute partitioning model

Understanding how Hofmeister salt ions and other solutes interact with proteins, nucleic acids, other biopolymers and water and thereby affect protein and nucleic acid processes as well as model processes (e.g. solubility of model compounds) in aqueous solution is a longstanding goal of biophysical research. Empirical Hofmeister salt and solute "m-values" (derivatives of the observed standard free energy change for a model or biopolymer process with respect to solute or salt concentration m3) are equal to differences in chemical potential derivatives: m-value = delta(dmu2/dm3) = delta mu23, which quantify the preferential interactions of the solute or salt with the surface of the biopolymer or model system (component 2) exposed or buried in the process. Using the solute partitioning model (SPM), we dissect mu23 values for interactions of a solute or Hofmeister salt with a set of model compounds displaying the key functional groups of biopolymers to obtain interaction potentials (called alpha-values) that quantify the interaction of the solute or salt per unit area of each functional group or type of surface. Interpreted using the SPM, these alpha-values provide quantitative information about both the hydration of functional groups and the competitive interaction of water and the solute or salt with functional groups. The analysis corroborates and quantifies previous proposals that the Hofmeister anion and cation series for biopolymer processes are determined by ion-specific, mostly unfavorable interactions with hydrocarbon surfaces; the balance between these unfavorable nonpolar interactions and often-favorable interactions of ions with polar functional groups determine the series null points. The placement of urea and glycine betaine (GB) at opposite ends of the corresponding series of nonelectrolytes results from the favorable interactions of urea, and unfavorable interactions of GB, with many (but not all) biopolymer functional groups. Interaction potentials and local-bulk partition coefficients quantifying the distribution of solutes (e.g. urea, glycine betaine) and Hofmeister salt ions in the vicinity of each functional group make good chemical sense when interpreted in terms of competitive noncovalent interactions. These interaction potentials allow solute and Hofmeister (noncoulombic) salt effects on protein and nucleic acid processes to be interpreted or predicted, and allow the use of solutes and salts as probes of

Denaturation of RNA secondary and tertiary structure by urea: simple unfolded state models and free energy parameters account for measured m-values

To investigate the mechanism by which urea destabilizes RNA structure, urea-induced unfolding of four different RNA secondary and tertiary structures was quantified in terms of an m-value, the rate at which the free energy of unfolding changes with urea molality. From literature data and our osmometric study of a backbone analogue, we derived average interaction potentials (per square angstrom of solvent accessible surface) between urea and three kinds of RNA surfaces: phosphate, ribose, and base. Estimates of the increases in solvent accessible surface areas upon RNA denaturation were based on a simple model of unfolded RNA as a combination of helical and single-strand segments. These estimates, combined with the three interaction potentials and a term to account for interactions of urea with released ions, yield calculated m-values that are in good agreement with experimental values (200 mm monovalent salt). Agreement was obtained only if single-stranded RNAs were modeled in a highly stacked, A-form conformation. The primary driving force for urea-induced denaturation is the strong interaction of urea with the large surface areas of bases that become exposed upon denaturation of either RNA secondary or tertiary structure, though interactions of urea with backbone and released ions may account for up to a third of the m-value. Urea m-values for all four RNAs are salt-dependent, which we attribute to an increased extension (or decreased charge density) of unfolded RNAs with an increased urea concentration. The sensitivity of the urea m-value to base surface exposure makes it a potentially useful probe of the conformations of RNA unfolded states.

Molecular mechanism for the preferential exclusion of TMAO from protein surfaces

Trimethylamine N-oxide (TMAO) is a naturally occurring protecting osmolyte that stabilizes the folded state of proteins and also counteracts the destabilizing effect of urea on protein stability. Experimentally, it has been inferred that TMAO is preferentially excluded from the vicinity of protein surfaces. Here, we combine computer modeling and experimental measurements to gain an understanding of the mechanism of the protecting effect of TMAO on proteins. We have developed an all-atom molecular model for TMAO that captures the exclusion of TMAO from model compounds and protein surfaces, as a consequence of incorporating realistic TMAO-water interactions through osmotic pressure measurements. Osmotic pressure measurements also suggest no significant attraction between urea and TMAO molecules in solution. To obtain an accurate potential for molecular simulations of protein stability in TMAO solutions, we have explored different ways of parametrizing the protein/osmolyte and osmolyte/osmolyte interactions by scaling charges and the strength of Lennard-Jones interactions and carried out equilibrium folding experiments of Trp-cage miniprotein in the presence of TMAO to guide the parametrization. Our calculations suggest a general principle for preferential interaction behavior of cosolvents with protein surfaces--preferentially excluded osmolytes have repulsive self-interaction given by osmotic coefficient φ > 1, while denaturants, in addition to having attractive interactions with the proteins, have favorable self-interaction given by osmotic coefficient φ < 1, to enable preferential accumulation in the vicinity of proteins.

Computational investigation on microsolvation of the osmolyte glycine betaine [GB (H(2)O)(1-7)]

The preferential interactions of glycine betaine (GB) with solvent components and the effect of solvent on its stability have been examined. In particular, the microsolvation of organic osmolyte and widely important osmoprotectant in nature as glycine betaine has been reported by using M06 method. A number of configurations (b(X) (a-z)) of the clusters for one to seven water molecules (× = 1-7) have been considered for the microsolvation. Structures of stable conformers are obtained and denoted as b1a, b2a, b3a, b4a, b5a, b6a and b7a. It is observed from the interaction energy difference (∆E) that only seven water molecules can be accommodated in the first solvation shell to stabilize GB. It is also observed that the calculated relative energy using M06 is in close agreement with calculations at the MP2 level of theory.

Weak interactions between folate and osmolytes in solution

Previous osmotic stress studies on the role of solvent in two structurally unrelated dihydrofolate reductases (DHFRs) found weaker binding of dihydrofolate (DHF) to either enzyme in the presence of osmolytes. To explain these unusual results, weak interactions between DHF and osmolytes were proposed, with a competition between osmolyte and DHFR for DHF. High osmolyte concentrations will inhibit binding of the cognate pair. To evaluate this hypothesis, we devised a small molecule approach. Dimerization of folate, monitored by nuclear magnetic resonance, was weakened 2-3-fold upon addition of betaine or dimethyl sulfoxide (DMSO), supporting preferential interaction of either osmolyte with the monomer (as it possesses a larger surface area). Nuclear Overhauser effect (NOE) spectroscopy experiments found a positive NOE for the interaction of the C3'/C5' benzoyl ring protons with the C9 proton in buffer; however, a negative NOE was observed upon addition of betaine or DMSO. This change indicated a decreased tumbling rate, consistent with osmolyte interaction. Osmotic stress experiments also showed that betaine, DMSO, and sucrose preferentially interact with folate. Further, studies with the folate fragments, p-aminobenzoic acid and pterin 6-carboxylate, revealed interactions for both model compounds with betaine and sucrose. In contrast, DMSO was strongly excluded from the pterin ring but preferentially interacted with the p-aminobenzoyl moiety. These interactions are likely to be important in vivo because of the crowded conditions of the cell where weak contacts can more readily compete with specific binding interactions.

Bicarbonate stabilizes isolated D1/D2/cytochrome b559 complex of photosystem 2 against thermoinactivation

It has been shown that thermoinactivation of the isolated D1/D2/cytochrome b(559) complex (RC) of photosystem 2 (PS-2) from pea under anaerobic conditions at 35°C in 20 mM Tris-HCl buffer (pH 7.2) depleted of HCO(3)(-), with 35 mM NaCl and 0.05% n-dodecyl-β-maltoside, results in a decrease in photochemical activity measured by photoreduction of the PS-2 primary electron acceptor, pheophytin (by 50% after 3 min of heating), which is accompanied by aggregation of the D1 and D2 proteins. Bicarbonate, formate, and acetate anions added to the sample under these conditions differently influence the maintenance of photochemical activity: a 50% loss of photochemical activity occurs in 11.5 min of heating in the presence of bicarbonate and in 4 and 4.6 min in the presence of formate and acetate, respectively. The addition of bicarbonate completely prevents aggregation of the D1 and D2 proteins as opposed to formate and acetate (their presence has no effect on the aggregation during thermoinactivation). Since the isolated RCs have neither inorganic Mn/Ca-containing core of the water-oxidizing complex nor nonheme Fe(2+), it is supposed that bicarbonate specifically interacts with the hydrophilic domains of the D1 and D2 proteins, which prevents their structural modification that is a signal for aggregation of these proteins and the loss of photochemical activity.

Quantifying why urea is a protein denaturant, whereas glycine betaine is a protein stabilizer

To explain the large, opposite effects of urea and glycine betaine (GB) on stability of folded proteins and protein complexes, we quantify and interpret preferential interactions of urea with 45 model compounds displaying protein functional groups and compare with a previous analysis of GB. This information is needed to use urea as a probe of coupled folding in protein processes and to tune molecular dynamics force fields. Preferential interactions between urea and model compounds relative to their interactions with water are determined by osmometry or solubility and dissected using a unique coarse-grained analysis to obtain interaction potentials quantifying the interaction of urea with each significant type of protein surface (aliphatic, aromatic hydrocarbon (C); polar and charged N and O). Microscopic local-bulk partition coefficients K(p) for the accumulation or exclusion of urea in the water of hydration of these surfaces relative to bulk water are obtained. K(p) values reveal that urea accumulates moderately at amide O and weakly at aliphatic C, whereas GB is excluded from both. These results provide both thermodynamic and molecular explanations for the opposite effects of urea and glycine betaine on protein stability, as well as deductions about strengths of amide NH--amide O and amide NH--amide N hydrogen bonds relative to hydrogen bonds to water. Interestingly, urea, like GB, is moderately accumulated at aromatic C surface. Urea m-values for protein folding and other protein processes are quantitatively interpreted and predicted using these urea interaction potentials or K(p) values.

The stability of intramolecular DNA G-quadruplexes compared with other macromolecules

DNA quadruplexes are often conceived as very stable structures. However, most of the free energy of stabilization derives from specific ion binding via inner sphere coordination of the GO6 of the guanine residues comprising the basic quartet. When compared with other nucleic acid structures such as DNA or RNA duplexes and hairpins, or proteins of the same number of atoms, metal-coordinated intramolecular quadruplexes are found to be of comparable or lower thermodynamic stability under similar solution conditions. Furthermore, intramolecular quadruplexes are actually less stable kinetically, than DNA duplexes or hairpins of the same size. Although the literature is incomplete, it is clear that polyelectrolyte ion effects, the influence of solvation and steric crowding on stability are qualitatively different between intramolecular quadruplexes and DNA duplexes. For example, decreasing water activity destabilizes DNA duplexes, whereas quadruplexes are stabilized. The variety of folded conformations accessible to a single sequence further implies strong sensitivity of the conformational ensemble to the solution conditions, compared with DNA duplexes or small single domain proteins. These considerations may have relevance to the conditions prevailing inside cell nuclei and therefore the structures that potentially might form in vivo.

Separation of preferential interaction and excluded volume effects on DNA duplex and hairpin stability

Small solutes affect protein and nucleic acid processes because of favorable or unfavorable chemical interactions of the solute with the biopolymer surface exposed or buried in the process. Large solutes also exclude volume and affect processes where biopolymer molecularity and/or shape changes. Here, we develop an analysis to separate and interpret or predict excluded volume and chemical effects of a flexible coil polymer on a process. We report a study of the concentration-dependent effects of the full series from monomeric to polymeric PEG on intramolecular hairpin and intermolecular duplex formation by 12-nucleotide DNA strands. We find that chemical effects of PEG on these processes increase in proportion to the product of the amount of DNA surface exposed on melting and the amount of PEG surface that is accessible to this DNA, and these effects are completely described by two interaction terms that quantify the interactions between this DNA surface and PEG end and interior groups. We find that excluded volume effects, once separated from these chemical effects, are quantitatively described by the analytical theory of Hermans, which predicts the excluded volume between a flexible polymer and a rigid molecule. From this analysis, we show that at constant concentration of PEG monomer, increasing PEG size increases the excluded volume effect but decreases the chemical interaction effect, because in a large PEG coil a smaller fraction of the monomers are accessible to the DNA. Volume exclusion by PEG has a much larger effect on intermolecular duplex formation than on intramolecular hairpin formation.

Interactions of the DNA polymerase X from African Swine Fever Virus with the ssDNA. Properties of the total DNA-binding site and the strong DNA-binding subsite

Interactions of the polymerase X from the African Swine Fever Virus with the ssDNA have been studied, using quantitative fluorescence titration and fluorescence resonance energy transfer techniques. The primary DNA-binding subsite of the enzyme, independent of the DNA conformation, is located on the C-terminal domain. Association of the bound DNA with the catalytic N-terminal domain finalizes the engagement of the total DNA-binding site of the enzyme and induces a large topological change in the structure of the bound ssDNA. The free energy of binding includes a conformational transition of the protein. Large positive enthalpy changes accompanying the ASFV pol X-ssDNA association indicate that conformational changes of the complex are induced by the engagement of the N-terminal domain. The enthalpy changes are offset by large entropy changes accompanying the DNA binding to the C-terminal domain and the total DNA-binding site, predominantly resulting from the release of water molecules.

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ∼5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.

The primary DNA-binding subsite of the rat pol β. Energetics of interactions of the 8-kDa domain of the enzyme with the ssDNA

Interactions of the 8-kDa domain of the rat pol β and the intact enzyme with the ssDNA have been studied, using the quantitative fluorescence titration technique. The 8-kDa domain induces large topological changes in the bound DNA structure and engages much larger fragments of the DNA than when embedded in the intact enzyme. The DNA affinity of the domain is predominantly driven by entropy changes, dominated by the water release from the protein. The thermodynamic characteristics dramatically change when the domain is embedded in the intact polymerase, indicating the presence of significant communication between the 8-kDa domain and the catalytic 31-kDa domain. The diminished water release from the 31-kDa domain strongly contributes to its dramatically lower DNA affinity, as compared to the 8-kDa domain. Unlike the 8-kDa domain, the DNA binding of the intact pol β is driven by entropy changes, originating from the structural changes of the formed complexes.

The BCCT family of carriers: from physiology to crystal structure

Increases in the environmental osmolarity are key determinants for the growth of microorganisms. To ensure a physiologically acceptable level of cellular hydration and turgor at high osmolarity, many bacteria accumulate compatible solutes. Osmotically controlled uptake systems allow the scavenging of these compounds from scarce environmental sources as effective osmoprotectants. A number of these systems belong to the BCCT family (betaine-choline-carnitine-transporter), sodium- or proton-coupled transporters (e.g. BetP and BetT respectively) that are ubiquitous in microorganisms. The BCCT family also contains CaiT, an L-carnitine/γ-butyrobetaine antiporter that is not involved in osmotic stress responses. The glycine betaine transporter BetP from Corynebacterium glutamicum is a representative for osmoregulated symporters of the BCCT family and functions both as an osmosensor and osmoregulator. The crystal structure of BetP in an occluded conformation in complex with its substrate glycine betaine and two crystal structures of CaiT in an inward-facing open conformation in complex with L-carnitine and γ-butyrobetaine were reported recently. These structures and the wealth of biochemical data on the activity control of BetP in response to osmotic stress enable a correlation between the sensing of osmotic stress by a transporter protein with the ensuing regulation of transport activity. Molecular determinants governing the high-affinity binding of the compatible solutes by BetP and CaiT, the coupling in symporters and antiporters, and the osmoregulatory properties are discussed in detail for BetP and various BCCT carriers.

Crowding alone cannot account for cosolute effect on amyloid aggregation

Amyloid fiber formation is a specific form of protein aggregation, often resulting from the misfolding of native proteins. Aimed at modeling the crowded environment of the cell, recent experiments showed a reduction in fibrillation halftimes for amyloid-forming peptides in the presence of cosolutes that are preferentially excluded from proteins and peptides. The effect of excluded cosolutes has previously been attributed to the large volume excluded by such inert cellular solutes, sometimes termed "macromolecular crowding". Here, we studied a model peptide that can fold to a stable monomeric β-hairpin conformation, but under certain solution conditions aggregates in the form of amyloid fibrils. Using Circular Dichroism spectroscopy (CD), we found that, in the presence of polyols and polyethylene glycols acting as excluded cosolutes, the monomeric β-hairpin conformation was stabilized with respect to the unfolded state. Stabilization free energy was linear with cosolute concentration, and grew with molecular volume, as would also be predicted by crowding models. After initiating the aggregation process with a pH jump, fibrillation in the presence and absence of cosolutes was followed by ThT fluorescence, transmission electron microscopy, and CD spectroscopy. Polyols (glycerol and sorbitol) increased the lag time for fibril formation and elevated the amount of aggregated peptide at equilibrium, in a cosolute size and concentration dependent manner. However, fibrillation rates remained almost unaffected by a wide range of molecular weights of soluble polyethylene glycols. Our results highlight the importance of other forces beyond the excluded volume interactions responsible for crowding that may contribute to the cosolute effects acting on amyloid formation.

The osmolyte TMAO stabilizes native RNA tertiary structures in the absence of Mg2+: evidence for a large barrier to folding from phosphate dehydration

The stabilization of RNA tertiary structures by ions is well known, but the neutral osmolyte trimethylamine oxide (TMAO) can also effectively stabilize RNA tertiary structure. To begin to understand the physical basis for the effects of TMAO on RNA, we have quantitated the TMAO-induced stabilization of five RNAs with known structures. So-called m values, the increment in unfolding free energy per molal of osmolyte at constant KCl activity, are ∼0 for a hairpin secondary structure and between 0.70 and 1.85 kcal mol(-1)m(-1) for four RNA tertiary structures (30-86 nt). Further analysis of two RNAs by small-angle X-ray scattering and hydroxyl radical probing shows that TMAO reduces the radius of gyration of the unfolded ensemble to the same endpoint as seen in titration with Mg(2+) and that the structures stabilized by TMAO and Mg(2+) are indistinguishable. Remarkably, TMAO induces the native conformation of a Mg(2+) ion chelation site formed in part by a buried phosphate, even though Mg(2+) is absent. TMAO interacts weakly, if at all, with KCl, ruling out the possibility that TMAO stabilizes RNA indirectly by increasing salt activity. TMAO is, however, strongly excluded from the vicinity of dimethylphosphate (unfavorable interaction free energy, +211 cal mol(-1)m(-1) for the potassium salt), an ion that mimics the RNA backbone phosphate. We suggest that formation of RNA tertiary structure is accompanied by substantial phosphate dehydration (loss of 66-173 water molecules in the RNA structures studied) and that TMAO works principally by reducing the energetic penalty associated with this dehydration. The strong parallels we find between the effects of TMAO and Mg(2+) suggest that RNA sequence is more important than specific ion interactions in specifying the native structure.

Zwitterion modulation of O(2)-evolving activity of cyanobacterial photosystem II

Photosystem II (PSII) is the only enzyme in nature that can catalyze the challenging catalytic photooxidation of H(2)O into four protons, four electrons, and O(2). Slowing down turnover of the O(2)-evolving complex (OEC) is a plausible approach to gain mechanistic information on the reaction. However, modulating the kinetics of the reaction without perturbing the active site is a challenge. In this study, it is shown that the steady-state activity of cyanobacterial PSII is inhibited by small zwitterions, such as glycine betaine and β-alanine. We show that the binding of zwitterions is nondenaturing, is highly reversible, and results in the decrease of the rate of catalytic turnover by ∼50% in the presence of excess zwitterion. Control measurements of photoinduced electron transfer in O(2)-inactive PSII show that the inhibition by zwitterions is the result of a specific decrease in the rate of catalytic turnover of the OEC. Recovery of activity upon addition of an exogenous proton carrier (HCO(3)(-)) provides evidence that proton-transfer pathways, thought to be essential for the relay of protons from the OEC to the lumen, are affected. Interestingly, no inhibition is observed for spinach PSII, suggesting that zwitterions act specifically by binding to the extrinsic proteins on the lumenal side of PSII, which differ significantly between plants and cyanobacteria, to slow proton transfer on the electron donor side of PSII.

Why Hofmeister effects of many salts favor protein folding but not DNA helix formation

The majority ( approximately 70%) of surface buried in protein folding is hydrocarbon, whereas in DNA helix formation, the majority ( approximately 65%) of surface buried is relatively polar nitrogen and oxygen. Our previous quantification of salt exclusion from hydrocarbon (C) accessible surface area (ASA) and accumulation at amide nitrogen (N) and oxygen (O) ASA leads to a prediction of very different Hofmeister effects on processes that bury mostly polar (N, O) surface compared to the range of effects commonly observed for processes that bury mainly nonpolar (C) surface, e.g., micelle formation and protein folding. Here we quantify the effects of salts on folding of the monomeric DNA binding domain (DBD) of lac repressor (lac DBD) and on formation of an oligomeric DNA duplex. In accord with this prediction, no salt investigated has a stabilizing Hofmeister effect on DNA helix formation. Our ASA-based analyses of model compound data and estimates of the surface area buried in protein folding and DNA helix formation allow us to predict Hofmeister effects on these processes. We observe semiquantitative to quantitative agreement between these predictions and the experimental values, obtained from a novel separation of coulombic and Hofmeister effects. Possible explanations of deviations, including salt-dependent unfolded ensembles and interactions with other types of surface, are discussed.

Preferential interactions between small solutes and the protein backbone: a computational analysis

To improve our understanding of the effects of small solutes on protein stability, we conducted atomistic simulations to quantitatively characterize the interactions between two broadly used small solutes, urea and glycine betaine (GB), and a triglycine peptide, which is a good model for a protein backbone. Multiple solute concentrations were analyzed, and each solute-peptide-water ternary system was studied with approximately 200-300 ns of molecular dynamics simulations with the CHARMM force field. The comparison between calculated preferential interaction coefficients (Gamma(23)) and experimentally measured values suggests that semiquantitative agreement with experiments can be obtained if care is exercised to balance interactions among the solute, protein, and water. On the other hand, qualitatively incorrect (i.e., wrong sign in Gamma(23)) results can be obtained if a solute model is constructed by directly taking parameters for chemically similar groups from an existing force field. Such sensitivity suggests that small solute thermodynamic data can be valuable in the development of accurate force field models of biomolecules. Further decomposition of Gamma(23) into group contributions leads to additional insights regarding the effects of small solutes on protein stability. For example, use of the CHARMM force field predicts that urea preferentially interacts with not only amide groups in the peptide backbone but also aliphatic groups, suggesting a role for these interactions in urea-induced protein denaturation; quantitatively, however, it is likely that the CHARMM force field overestimates the interaction between urea and aliphatic groups. The results with GB support a simple thermodynamic model that assumes additivity of preferential interaction between GB and various biomolecular surfaces.

Contributions of the histidine side chain and the N-terminal alpha-amino group to the binding thermodynamics of oligopeptides to nucleic acids as a function of pH

Interactions of histidine with nucleic acid phosphates and histidine pK(a) shifts make important contributions to many protein-nucleic acid binding processes. To characterize these phenomena in simplified systems, we quantified binding of a histidine-containing model peptide HWKK ((+)NH(3)-His-Trp-Lys-Lys-NH(2)) and its lysine analogue KWKK ((+)NH(3)-Lys-Trp-Lys-Lys-NH(2)) to a single-stranded RNA model, polyuridylate (polyU), by changes in tryptophan fluorescence as a function of salt concentration and pH. For both HWKK and KWKK, equilibrium binding constants, K(obs), and magnitudes of log-log salt derivatives, SK(obs) identical with (partial differential logK(obs)/partial differential log[Na(+)]), decreased with increasing pH in the manner expected for a titration curve model in which deprotonation of the histidine and alpha-amino groups weakens binding and reduces its salt-dependence. Fully protonated HWKK and KWKK exhibit the same K(obs) and SK(obs) within uncertainty, and these SK(obs) values are consistent with limiting-law polyelectrolyte theory for +4 cationic oligopeptides binding to single-stranded nucleic acids. The pH-dependence of HWKK binding to polyU provides no evidence for pK(a) shifts nor any requirement for histidine protonation, in stark contrast to the thermodynamics of coupled protonation often seen for these cationic residues in the context of native protein structure where histidine protonation satisfies specific interactions (e.g., salt-bridge formation) within highly complementary binding interfaces. The absence of pK(a) shifts in our studies indicates that additional Coulombic interactions across the nonspecific-binding interface between RNA and protonated histidine or the alpha-amino group are not sufficient to promote proton uptake for these oligopeptides. We present our findings in the context of hydration models for specific vs nonspecific nucleic acid binding.