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

Salt Effects on Caffeine across Concentration Regimes

Salts affect the solvation thermodynamics of molecules of all sizes; the Hofmeister series is a prime example in which different ions lead to salting-in or salting-out of aqueous proteins. Early work of Tanford led to the discovery that the solvation of molecular surface motifs is proportional to the solvent accessible surface area (SASA), and later studies have shown that the proportionality constant varies with the salt concentration and type. Using multiscale computer simulations combined with vapor-pressure osmometry on caffeine-salt solutions, we reveal that this SASA description captures a rich set of molecular driving forces in tertiary solutions at changing solute and osmolyte concentrations. Central to the theoretical work is a new potential energy function that depends on the instantaneous surface area, salt type, and concentration. Used in, e.g., Monte Carlo simulations, this allows for a highly efficient exploration of many-body interactions and the resulting thermodynamics at elevated solute and salt concentrations.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally Determined Enthalpic and Entropic Contributions to Interactions of Amide sp2O, N, C and sp3C Unified Atoms with Naphthalene sp2C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp2 atoms of amides, aromatics, and other groups occur in protein self-assembly processes including folding, oligomerization, and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp2O with amide sp2N unified atoms (presumably C═O···H-N hydrogen bonds) and amide/aromatic sp2C (lone pair π, n-π*) are particularly favorable. Sp3C-sp3C (hydrophobic), sp3C-sp2C (hydrophobic, CH-π), sp2C-sp2C (hydrophobic, π-π), and sp3C-sp2N interactions are favorable, sp2C-sp2N interactions are neutral, while sp2O-sp2O and sp2N-sp2N self-interactions and sp2O-sp3C interactions are unfavorable. Here, from determinations of favorable effects of 14 amides on naphthalene solubility at 10, 25, and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp2O, sp2N, sp2C, and sp3C unified atoms with aromatic sp2C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp2O-aromatic sp2C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g., lone pair-π), while amide sp3C- and sp2C-aromatic sp2C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp2 atoms in protein processes.

Metal-binding and circular permutation-dependent thermodynamic and kinetic stability of azurin

Native topology is known to determine the folding kinetics and the energy landscape of proteins. Furthermore, the circular permutation (CP) of proteins alters the order of the secondary structure connectivity while retaining the three-dimensional structure, making it an elegant and powerful approach to altering native topology. Previous studies elucidated the influence of CP in proteins with different folds such as Greek key β-barrel, β-sandwich, β-α-β, and all α-Greek key. CP mainly affects the protein stability and unfolding kinetics, while folding kinetics remains mostly unaltered. However, the effect of CP on metalloproteins is yet to be elaborately studied. The active site of metalloproteins poses an additional complexity in studying protein folding. Here, we investigate a CP variant (cpN42) of azurin-in both metal-free and metal-bound (holo) forms. As observed earlier in other proteins, apo-forms of wild-type (WT) and cpN42 fold with similar rates. In contrast, zinc-binding accelerates the folding of WT but decelerates the folding of cpN42. On zinc-binding, the spontaneous folding rate of WT increases by >250 times that of cpN42, which is unprecedented and the highest for any CP to date. On the other hand, zinc-binding reduces the spontaneous unfolding rate of cpN42 by ~100 times, making the WT and CP azurins unfold at similar rates. Our study demonstrates metal binding as a novel way to modulate the unfolding and folding rates of CPs compared to their WT counterparts. We hope our study increases the understanding of the effect of CP on the folding mechanism and energy landscape of metalloproteins.

Clamp Interactions with +3/+6 Duplex and Upstream-to-Downstream Allosteric Effects in Late Steps of Forming a Stable RNA Polymerase-Promoter Open Complex

Stable 37 °C open complexes (OC) of E. coli RNA polymerase (RNAP) at λP_R and T7A1 promoters form at similar rates but have very different lifetimes. To understand the downstream interactions responsible for OC lifetime, how promoter sequence directs them and when they form, we report lifetimes of stable OC and unstable late (I₂) intermediates for promoters with different combinations of λP_R (L) and T7A1 (T) discriminators, core promoters and UP elements. I₂ lifetimes are similarly short, while stable OC lifetimes differ greatly, determined largely by the discriminator and modulated by core-promoter and UP elements. The free energy change ΔG₃^o for I₂ → stable OC is approximately -4 kcal more favorable for L-discriminator than for T-discriminator promoters. Downstream-truncation at +6 (DT+6) greatly destabilizes OC at L-discriminator but not T-discriminator promoters, making all ΔG₃^o values similar (approximately -4 kcal). Urea reduces OC lifetime greatly by affecting ΔG₃^o. We deduce that urea acts by disfavoring coupled folding of key elements of the β'-clamp, that I₂ is an open-clamp OC, and that clamp-closing in I₂ → stable OC involves coupled folding. Differences in ΔG₃^o between downstream-truncated and full-length promoters yield contributions to ΔG₃^o from interactions with downstream mobile elements (DME) including β-lobe and β'-jaw, more favorable for L-discriminator than for T-discriminator promoters. We deduce how competition between far-downstream DNA and σ⁷⁰ region 1.1 affects ΔG₃^o values. We discuss variant-specific ΔG₃^o contributions in terms of the allosteric network by which differences in discriminator and -10 sequence are sensed and transmitted downstream to affect DME-duplex interactions in I₂ → stable OC.

Experimentally determined strengths of favorable and unfavorable interactions of amide atoms involved in protein self-assembly in water

Folding and other protein self-assembly processes are driven by favorable interactions between O, N, and C unified atoms of the polypeptide backbone and side chains. These processes are perturbed by solutes that interact with these atoms differently than water does. Amide NH···O=C hydrogen bonding and various π-system interactions have been better characterized structurally or by simulations than experimentally in water, and unfavorable interactions are relatively uncharacterized. To address this situation, we previously quantified interactions of alkyl ureas with amide and aromatic compounds, relative to interactions with water. Analysis yielded strengths of interaction of each alkylurea with unit areas of different hybridization states of unified O, N, and C atoms of amide and aromatic compounds. Here, by osmometry, we quantify interactions of 10 pairs of amides selected to complete this dataset. An analysis yields intrinsic strengths of six favorable and four unfavorable atom-atom interactions, expressed per unit area of each atom and relative to interactions with water. The most favorable interactions are sp2O-sp2C (lone pair-π, presumably n-π*), sp2C-sp2C (π-π and/or hydrophobic), sp2O-sp2N (hydrogen bonding) and sp3C-sp2C (CH-π and/or hydrophobic). Interactions of sp3C with itself (hydrophobic) and with sp2N are modestly favorable, while sp2N interactions with sp2N and with amide/aromatic sp2C are modestly unfavorable. Amide sp2O-sp2O interactions and sp2O-sp3C interactions are more unfavorable, indicating the preference of amide sp2O to interact with water. These intrinsic interaction strengths are used to predict interactions of amides with proteins and chemical effects of amides (including urea, N-ethylpyrrolidone [NEP], and polyvinylpyrrolidone [PVP]) on protein stability.

Using Single-Molecule Chemo-Mechanical Unfolding to Simultaneously Probe Multiple Structural Parameters in Protein Folding

While single-molecule force spectroscopy has greatly advanced the study of protein folding, there are limitations to what can be learned from studying the effect of force alone. We developed a novel technique, chemo-mechanical unfolding, that combines multiple perturbants—force and chemical denaturant—to more fully characterize the folding process by simultaneously probing multiple structural parameters—the change in end-to-end distance, and solvent accessible surface area. Here, we describe the theoretical background, experimental design, and data analysis for chemo-mechanical unfolding experiments probing protein folding thermodynamics and kinetics. This technique has been applied to characterize parallel protein folding pathways, the protein denatured state, protein folding on the ribosome, and protein folding intermediates.

A small single-domain protein folds through the same pathway on and off the ribosome

In vivo, proteins fold and function in a complex environment subject to many stresses that can modulate a protein's energy landscape. One aspect of the environment pertinent to protein folding is the ribosome, since proteins have the opportunity to fold while still bound to the ribosome during translation. We use a combination of force and chemical denaturant (chemomechanical unfolding), as well as point mutations, to characterize the folding mechanism of the src SH3 domain both as a stalled ribosome nascent chain and free in solution. Our results indicate that src SH3 folds through the same pathway on and off the ribosome. Molecular simulations also indicate that the ribosome does not affect the folding pathway for this small protein. Taken together, we conclude that the ribosome does not alter the folding mechanism of this small protein. These results, if general, suggest the ribosome may exert a bigger influence on the folding of multidomain proteins or protein domains that can partially fold before the entire domain sequence is outside the ribosome exit tunnel.

In Vivo Titration of Folate Pathway Enzymes

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

Salt Effects on Protein Folding Thermodynamics

Salts differ in their ability to stabilize protein conformations, thereby affecting the thermodynamics and kinetics of protein folding. We developed a coarse-grained protein model that can predict salt-induced changes in protein properties by using the transfer free-energy data of various chemical groups from water to salt solutions. Using this model and molecular dynamics simulations, we probed the effect of seven different salts on the folding thermodynamics of the DNA binding domain of lac repressor protein ( lac-DBD) and N-terminal domain of ribosomal protein (NTL9). We show that a salt can act as a protein stabilizing or destabilizing agent depending on the protein sequence and folded state topology. The computed thermodynamic properties, especially the m values for various salts, which reveal the relative ability of a salt to stabilize the protein folded state, are in quantitative agreement with the experimentally measured values. The computations show that the degree of protein compaction in the denatured ensemble strongly depends on the salt identity, and for the same variation in salt concentration, the compaction in the protein dimensions varies from ∼4% to ∼30% depending on the salt. The transition-state ensemble (TSE) of lac-DBD is homogeneous and polarized, while the TSE of NTL9 is heterogeneous and diffusive. Salts induce subtle structural changes in the TSE that are in agreement with Hammond's postulate. The barrier to protein folding tends to disappear in the presence of moderate concentrations (∼3-4 m) of strongly stabilizing salts.

Quantifying Interactions of Nucleobase Atoms with Model Compounds for the Peptide Backbone and Glutamine and Asparagine Side Chains in Water

Alkylureas display hydrocarbon and amide groups, the primary functional groups of proteins. To obtain the thermodynamic information that is needed to analyze interactions of amides and proteins with nucleobases and nucleic acids, we quantify preferential interactions of alkylureas with nucleobases differing in the amount and composition of water-accessible surface area (ASA) by solubility assays. Using an established additive ASA-based analysis, we interpret these thermodynamic results to determine interactions of each alkylurea with five types of nucleobase unified atoms (carbonyl sp²O, amino sp³N, ring sp²N, methyl sp³C, and ring sp²C). All alkylureas interact favorably with nucleobase sp²C and sp³C atoms; these interactions become more favorable with an increasing level of alkylation of urea. Interactions with nucleobase sp²O are most favorable for urea, less favorable for methylurea and ethylurea, and unfavorable for dialkylated ureas. Contributions to overall alkylurea-nucleobase interactions from interactions with each nucleobase atom type are proportional to the ASA of that atom type with proportionality constant (interaction strength) α, as observed previously for urea. Trends in α-values for interactions of alkylureas with nucleobase atom types parallel those for corresponding amide compound atom types, offset because nucleobase α-values are more favorable. Comparisons between ethylated and methylated ureas show interactions of amide compound sp³C with nucleobase sp²C, sp³C, sp²N, and sp³N atoms are favorable while amide sp³C-nucleobase sp²O interactions are unfavorable. Strongly favorable interactions of urea with nucleobase sp²O but weakly favorable interactions with nucleobase sp³N indicate that amide sp²N-nucleobase sp²O and nucleobase sp³N-amide sp²O hydrogen bonding (NH···O═C) interactions are favorable while amide sp²N-nucleobase sp³N interactions are unfavorable. These favorable amide-nucleobase hydrogen bonding interactions are prevalent in specific protein-nucleotide complexes.

The mechanism and high-free-energy transition state of lac repressor-lac operator interaction

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)–lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (k_a) and dissociation (k_d) rate constants. RO-destabilizing solutes (urea, KCl) reduce k_a comparably (urea) or more than (KCl) they increase k_d, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce k_d without affecting k_a significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO.

Partially native intermediates mediate misfolding of SOD1 in single-molecule folding trajectories

Prion-like misfolding of superoxide dismutase 1 (SOD1) is associated with the disease ALS, but the mechanism of misfolding remains unclear, partly because misfolding is difficult to observe directly. Here we study the most misfolding-prone form of SOD1, reduced un-metallated monomers, using optical tweezers to measure unfolding and refolding of single molecules. We find that the folding is more complex than suspected, resolving numerous previously undetected intermediate states consistent with the formation of individual β-strands in the native structure. We identify a stable core of the protein that unfolds last and refolds first, and directly observe several distinct misfolded states that branch off from the native folding pathways at specific points after the formation of the stable core. Partially folded intermediates thus play a crucial role mediating between native and non-native folding. These results suggest an explanation for SOD1's propensity for prion-like misfolding and point to possible targets for therapeutic intervention.

l-Proline and RNA Duplex m-Value Temperature Dependence

The temperature dependence of l-proline interactions with the RNA dodecamer duplex surface exposed after unfolding was quantified using thermal and isothermal titration denaturation monitored by uv-absorbance. The m-value quantifying proline interactions with the RNA duplex surface area exposed after unfolding was measured using RNA duplexes with GC content ranging between 17 and 83%. The m-values from thermal denaturation decreased with increasing GC content signifying increasingly favorable proline interactions with the exposed RNA surface area. However, m-values from isothermal titration denaturation at 25.0 °C were independent of GC content and less negative than those from thermal denaturation. The m-value from isothermal titration denaturation for a 50% GC RNA duplex decreased (became more negative) as the temperature increased and was in nearly exact agreement with the m-value from thermal denaturation. Since RNA duplex transition temperatures increased with GC content, the more favorable proline interactions with the high GC content duplex surface area observed from thermal denaturation resulted from the temperature dependence of proline interactions rather than the RNA surface chemical composition. The enthalpy contribution to the m-value was positive and small (indicating a slight increase in duplex unfolding enthalpy with proline) while the entropic contribution to the m-value was positive and increased with temperature. Our results will facilitate proline's use as a probe of solvent accessible surface area changes during biochemical reactions at different reaction temperatures.

Experimental Atom-by-Atom Dissection of Amide-Amide and Amide-Hydrocarbon Interactions in H2O

Quantitative information about amide interactions in water is needed to understand their contributions to protein folding and amide effects on aqueous processes and to compare with computer simulations. Here we quantify interactions of urea, alkylated ureas, and other amides by osmometry and amide-aromatic hydrocarbon interactions by solubility. Analysis of these data yields strengths of interaction of ureas and naphthalene with amide sp2O, amide sp2N, aliphatic sp3C, and amide and aromatic sp2C unified atoms in water. Interactions of amide sp2O with urea and naphthalene are favorable, while amide sp2O-alkylurea interactions are unfavorable, becoming more unfavorable with increasing alkylation. Hence, amide sp2O-amide sp2N interactions (proposed n-σ* hydrogen bond) and amide sp2O-aromatic sp2C (proposed n-π*) interactions are favorable in water, while amide sp2O-sp3C interactions are unfavorable. Interactions of all ureas with sp3C and amide sp2N are favorable and increase in strength with increasing alkylation, indicating favorable sp3C-amide sp2N and sp3C-sp3C interactions. Naphthalene results show that aromatic sp2C-amide sp2N interactions in water are unfavorable while sp2C-sp3C interactions are favorable. These results allow interactions of amide and hydrocarbon moieties and effects of urea and alkylureas on aqueous processes to be predicted or interpreted in terms of structural information. We predict strengths of favorable urea-benzene and N-methylacetamide interactions from experimental information to compare with simulations and indicate how amounts of hydrocarbon and amide surfaces buried in protein folding and other biopolymer processes and transition states can be determined from analysis of urea and diethylurea effects on equilibrium and rate constants.

Basis of Protein Stabilization by K Glutamate: Unfavorable Interactions with Carbon, Oxygen Groups

Potassium glutamate (KGlu) is the primary Escherichia coli cytoplasmic salt. After sudden osmotic upshift, cytoplasmic KGlu concentration increases, initially because of water efflux and subsequently by K⁺ transport and Glu^- synthesis, allowing water uptake and resumption of growth at high osmolality. In vitro, KGlu ranks with Hofmeister salts KF and K₂SO₄ in driving protein folding and assembly. Replacement of KCl by KGlu stabilizes protein-nucleic acid complexes. To interpret and predict KGlu effects on protein processes, preferential interactions of KGlu with 15 model compounds displaying six protein functional groups-sp³ (aliphatic) C; sp² (aromatic, amide, carboxylate) C; amide and anionic (carboxylate) O; and amide and cationic N-were determined by osmometry or solubility assays. Analysis of these data yields interaction potentials (α-values) quantifying non-Coulombic chemical interactions of KGlu with unit area of these six groups. Interactions of KGlu with the 15 model compounds predicted from these six α-values agree well with experimental data. KGlu interactions with all carbon groups and with anionic (carboxylate) and amide oxygen are unfavorable, while KGlu interactions with cationic and amide nitrogen are favorable. These α-values, together with surface area information, provide quantitative predictions of why KGlu is an effective E. coli cytoplasmic osmolyte (because of the dominant effect of unfavorable interactions of KGlu with anionic and amide oxygens and hydrocarbon groups on the water-accessible surface of cytoplasmic biopolymers) and why KGlu is a strong stabilizer of folded proteins (because of the dominant effect of unfavorable interactions of KGlu with hydrocarbon groups and amide oxygens exposed in unfolding).

TMAO-Protein Preferential Interaction Profile Determines TMAO's Conditional In Vivo Compatibility

Trimethylamine N-oxide (TMAO) exemplifies how Nature uses the solute effect as a simple chemical strategy to cope with hydrodynamic pressure or urea stress to maintain proteostasis. It is a gut-microbe-generated metabolite that strongly promotes the development of atherosclerosis. It remains unclear how TMAO exerts its effects. In this study, we experimentally characterized the profile of the preferential interaction potential of TMAO with proteins, a thermodynamic key to understanding the effects of TMAO on protein processes and the distinction of TMAO among osmolytes. TMAO is thus found to be highly preferentially excluded from most types of protein surface, which explains why TMAO is a strong globular protein stabilizer and identifies the dominant stabilizing factor as the unfavorable interaction of TMAO with the hydrophobic surface exposed upon unfolding. We dissected the mechanism of the counteracting effects of TMAO and urea: the contrary feature of the interaction profiles of the two solutes maximizes the possibility for them to offset each other's perturbing effect on protein processes. The interaction profile also predicts that TMAO promotes aggregation of amyloidogenic intrinsically disordered peptide, as demonstrated here in Aβ42, and that TMAO has a strong potential to impact protein processes in the absence of stressors. Our data suggest that although TMAO is an evolutionally selected chemical chaperone for some organisms or organs, its compatibility in vivo is conditional and determined by its interaction profile with biopolymers and the nature of the essential biopolymer processes. Our thermodynamic framework plus the TMAO-protein interaction profile provides a basis for exploring the broad biological significance of TMAO, including its pathological impact in the absence of stressors. We argue for the general importance of controlling in vivo background solutes and the pathological significance of a control failure.

Perplexing cooperative folding and stability of a low-sequence complexity, polyproline 2 protein lacking a hydrophobic core

The burial of hydrophobic side chains in a protein core generally is thought to be the major ingredient for stable, cooperative folding. Here, we show that, for the snow flea antifreeze protein (sfAFP), stability and cooperativity can occur without a hydrophobic core, and without α-helices or β-sheets. sfAFP has low sequence complexity with 46% glycine and an interior filled only with backbone H-bonds between six polyproline 2 (PP2) helices. However, the protein folds in a kinetically two-state manner and is moderately stable at room temperature. We believe that a major part of the stability arises from the unusual match between residue-level PP2 dihedral angle bias in the unfolded state and PP2 helical structure in the native state. Additional stabilizing factors that compensate for the dearth of hydrophobic burial include shorter and stronger H-bonds, and increased entropy in the folded state. These results extend our understanding of the origins of cooperativity and stability in protein folding, including the balance between solvent and polypeptide chain entropies.

Subcellular modulation of protein VlsE stability and folding kinetics

The interior of a cell interacts differently with proteins than a dilute buffer because of a wide variety of macromolecules, chaperones, and osmolytes that crowd and interact with polypeptide chains. We compare folding of fluorescent constructs of protein VlsE among three environments inside cells. The nucleus increases the stability of VlsE relative to the cytoplasm, but slows down folding kinetics. VlsE is also more stable in the endoplasmic reticulum, but unlike PGK, tends to aggregate there. Although fluorescent-tagged VlsE and PGK show opposite stability trends from in vitro to the cytoplasm, their trends from cytoplasm to nucleus are similar.

Positioning the Intracellular Salt Potassium Glutamate in the Hofmeister Series by Chemical Unfolding Studies of NTL9

In vitro, replacing KCl with potassium glutamate (KGlu), the Escherichia coli cytoplasmic salt and osmolyte, stabilizes folded proteins and protein-nucleic acid complexes. To understand the chemical basis for these effects and rank Glu- in the Hofmeister anion series for protein unfolding, we quantify and interpret the strong stabilizing effect of KGlu on the ribosomal protein domain NTL9, relative to the effects of other stabilizers (KCl, KF, and K2SO4) and destabilizers (GuHCl and GuHSCN). GuHSCN titrations at 20 ° C, performed as a function of the concentration of KGlu or another salt and monitored by NTL9 fluorescence, are analyzed to obtain R-values quantifying the Hofmeister salt concentration (m3) dependence of the unfolding equilibrium constant K(obs) [r-value = −d ln K(obs)/dm3 = (1/RT) dΔG(obs) ° /dm3 = m-value/RT]. r-Values for both stabilizing K+ salts and destabilizing GuH+ salts are compared with predictions from model compound data. For two-salt mixtures, we find that contributions of stabilizing and destabilizing salts to observed r-values are additive and independent. At 20 ° C, we determine a KGlu r-value of 3.22 m(−1) and K2SO4, KF, KCl, GuHCl, and GuHSCN r-values of 5.38, 1.05, 0.64, −1.38, and −3.00 m(−1), respectively. The KGlu r-value represents a 25-fold (1.9 kcal) stabilization per molal KGlu added. KGlu is much more stabilizing than KF, and the stabilizing effect of KGlu is larger in magnitude than the destabilizing effect of GuHSCN. Interpretation of the data reveals good agreement between predicted and observed relative r-values and indicates the presence of significant residual structure in GuHSCN-unfolded NTL9 at 20 ° C.

Introduction of a polar core into the de novo designed protein Top7

Design of polar interactions is a current challenge for protein design. The de novo designed protein Top7, like almost all designed proteins, has an entirely nonpolar core. Here we describe the replacing of a sizable fraction (5 residues) of this core with a designed polar hydrogen bond network. The polar core design is expressed at high levels in E. coli, has a folding free energy of 10 kcal/mol, and retains the multiphasic folding kinetics of the original Top7. The NMR structure of the design shows that conformations of three of the five residues, and the designed hydrogen bonds between them, are very close to those in the design model. The remaining two residues, which are more solvent exposed, sample a wide range of conformations in the NMR ensemble. These results show that hydrogen bond networks can be designed in protein cores, but also highlight challenges that need to be overcome when there is competition with solvent.

Even with nonnative interactions, the updated folding transition states of the homologs Proteins G & L are extensive and similar

Experimental and computational folding studies of Proteins L & G and NuG2 typically find that sequence differences determine which of the two hairpins is formed in the transition state ensemble (TSE). However, our recent work on Protein L finds that its TSE contains both hairpins, compelling a reassessment of the influence of sequence on the folding behavior of the other two homologs. We characterize the TSEs for Protein G and NuG2b, a triple mutant of NuG2, using ψ analysis, a method for identifying contacts in the TSE. All three homologs are found to share a common and near-native TSE topology with interactions between all four strands. However, the helical content varies in the TSE, being largely absent in Proteins G & L but partially present in NuG2b. The variability likely arises from competing propensities for the formation of nonnative β turns in the naturally occurring proteins, as observed in our TerItFix folding algorithm. All-atom folding simulations of NuG2b recapitulate the observed TSEs with four strands for 5 of 27 transition paths [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517-520]. Our data support the view that homologous proteins have similar folding mechanisms, even when nonnative interactions are present in the transition state. These findings emphasize the ongoing challenge of accurately characterizing and predicting TSEs, even for relatively simple proteins.

Single-molecule chemo-mechanical unfolding reveals multiple transition state barriers in a small single-domain protein

A fundamental question in protein folding is whether proteins fold through one or multiple trajectories. While most experiments indicate a single pathway, simulations suggest proteins can fold through many parallel pathways. Here, we use a combination of chemical denaturant, mechanical force and site-directed mutations to demonstrate the presence of multiple unfolding pathways in a simple, two-state folding protein. We show that these multiple pathways have structurally different transition states, and that seemingly small changes in protein sequence and environment can strongly modulate the flux between the pathways. These results suggest that in vivo, the crowded cellular environment could strongly influence the mechanisms of protein folding and unfolding. Our study resolves the apparent dichotomy between experimental and theoretical studies, and highlights the advantage of using a multipronged approach to reveal the complexities of a protein's free-energy landscape.

Although most protein folding experiments can be explained by a single pathway, theoretical evidence suggests the presence of multiple pathways. Here, the authors resolve this using a combination of force, chemical denaturation and mutagenesis to modulate the flux between parallel pathways.

Dispersion interactions between urea and nucleobases contribute to the destabilization of RNA by urea in aqueous solution

Urea has long been used to investigate protein folding and, more recently, RNA folding. Studies have proposed that urea denatures RNA by participating in stacking interactions and hydrogen bonds with nucleic acid bases. In this study, the ability of urea to form unconventional stacking interactions with RNA bases is investigated using ab initio calculations (RI-MP2 and CCSD(T) methods with the aug-cc-pVDZ basis set). A total of 29 stable nucleobase-urea stacked complexes are identified in which the intermolecular interaction energies (up to -14 kcal/mol) are dominated by dispersion effects. Natural bond orbital (NBO) and atoms in molecules (AIM) calculations further confirm strong interactions between urea and nucleobases. Calculations on model systems with multiple urea and water molecules interacting with a guanine base lead to a hypothesis that urea molecules along with water are able to form cage-like structures capable of trapping nucleic acid bases in extrahelical states by forming both hydrogen-bonded and dispersion interactions, thereby contributing to the unfolding of RNA in the presence of urea in aqueous solution.

Using solutes and kinetics to probe large conformational changes in the steps of transcription initiation

Small solutes are useful probes of large conformational changes in RNA polymerase-promoter interactions and other biopolymer processes. In general, a large effect of a solute on an equilibrium constant (or rate constant) indicates a large change in water-accessible biopolymer surface area in the corresponding step (or transition state), resulting from conformational changes, interface formation, or both. Here, we describe nitrocellulose filter binding assays from series used to determine the urea dependence of open complex formation and dissociation with Escherichia coli RNA polymerase and phage λPR promoter DNA. Then, we describe the subsequent data analysis and interpretation of these solute effects.

Benchmarking all-atom simulations using hydrogen exchange

Long-time molecular dynamics (MD) simulations are now able to fold small proteins reversibly to their native structures [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517-520]. These results indicate that modern force fields can reproduce the energy surface near the native structure. To test how well the force fields recapitulate the other regions of the energy surface, MD trajectories for a variant of protein G are compared with data from site-resolved hydrogen exchange (HX) and other biophysical measurements. Because HX monitors the breaking of individual H-bonds, this experimental technique identifies the stability and H-bond content of excited states, thus enabling quantitative comparison with the simulations. Contrary to experimental findings of a cooperative, all-or-none unfolding process, the simulated denatured state ensemble, on average, is highly collapsed with some transient or persistent native 2° structure. The MD trajectories of this protein G variant and other small proteins exhibit excessive intramolecular H-bonding even for the most expanded conformations, suggesting that the force fields require improvements in describing H-bonding and backbone hydration. Moreover, these comparisons provide a general protocol for validating the ability of simulations to accurately capture rare structural fluctuations.

Energetically significant networks of coupled interactions within an unfolded protein

Unfolded and partially unfolded proteins participate in a wide range of biological processes from pathological aggregation to the regulation of normal cellular activity. Unfolded states can be populated under strongly denaturing conditions, but the ensemble which is relevant for folding, stability, and aggregation is that populated under physiological conditions. Characterization of nonnative states is critical for the understanding of these processes, yet comparatively little is known about their energetics and their structural propensities under native conditions. The standard view is that energetically significant coupled interactions involving multiple residues are generally not present in the denatured state ensemble (DSE) or in intrinsically disordered proteins. Using the N-terminal domain of the ribosomal protein L9, a small α-β protein, as an experimental model system, we demonstrate that networks of energetically significant, coupled interactions can form in the DSE of globular proteins, and can involve residues that are distant in sequence and spatially well separated in the native structure. X-ray crystallography, NMR, dynamics studies, native state pKa measurements, and thermodynamic analysis of more than 25 mutants demonstrate that residues are energetically coupled in the DSE. Altering these interactions by mutation affects the stability of the domain. Mutations that alter the energetics of the DSE can impact the analysis of cooperativity and folding, and may play a role in determining the propensity to aggregate.

General mechanism of two-state protein folding kinetics

We describe here a general model of the kinetic mechanism of protein folding. In the Foldon Funnel Model, proteins fold in units of secondary structures, which form sequentially along the folding pathway, stabilized by tertiary interactions. The model predicts that the free energy landscape has a volcano shape, rather than a simple funnel, that folding is two-state (single-exponential) when secondary structures are intrinsically unstable, and that each structure along the folding path is a transition state for the previous structure. It shows how sequential pathways are consistent with multiple stochastic routes on funnel landscapes, and it gives good agreement with the 9 order of magnitude dependence of folding rates on protein size for a set of 93 proteins, at the same time it is consistent with the near independence of folding equilibrium constant on size. This model gives estimates of folding rates of proteomes, leading to a median folding time in Escherichia coli of about 5 s.

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.