Negative allosteric modulation of the glucagon receptor by RAMP2.

Receptor activity-modifying proteins (RAMPs) modulate the activity of many Family B GPCRs. We show that RAMP2 directly interacts with the glucagon receptor (GCGR), a Family B GPCR responsible for blood sugar homeostasis, and broadly inhibits receptor-induced downstream signaling. HDX-MS experiments demonstrate that RAMP2 enhances local flexibility in select locations in and near the receptor extracellular domain (ECD) and in the 6th transmembrane helix, whereas smFRET experiments show that this ECD disorder results in the inhibition of active and intermediate states of the intracellular surface. We determined the cryo-EM structure of the GCGR-Gs complex at 2.9 Å resolution in the presence of RAMP2. RAMP2 apparently does not interact with GCGR in an ordered manner; however, the receptor ECD is indeed largely disordered along with rearrangements of several intracellular hallmarks of activation. Our studies suggest that RAMP2 acts as a negative allosteric modulator of GCGR by enhancing conformational sampling of the ECD.

Stepwise activation of a metabotropic glutamate receptor.

Metabotropic glutamate receptors belong to a family of G protein-coupled receptors that are obligate dimers and possess a large extracellular ligand-binding domain that is linked via a cysteine-rich domain to their 7-transmembrane domain1. Upon activation, these receptors undergo a large conformational change to transmit the ligand binding signal from the extracellular ligand-binding domain to the G protein-coupling 7-transmembrane domain2. In this manuscript, we propose a model for a sequential, multistep activation mechanism of metabotropic glutamate receptor subtype 5. We present a series of structures in lipid nanodiscs, from inactive to fully active, including agonist-bound intermediate states. Further, using bulk and single-molecule fluorescence imaging, we reveal distinct receptor conformations upon allosteric modulator and G protein binding.

Selection pressures on evolution of ribonuclease H explored with rigorous free-energy-based design.

Understanding natural protein evolution and designing novel proteins are motivating interest in development of high-throughput methods to explore large sequence spaces. In this work, we demonstrate the application of multisite λ dynamics (MSλD), a rigorous free energy simulation method, and chemical denaturation experiments to quantify evolutionary selection pressure from sequence-stability relationships and to address questions of design. This study examines a mesophilic phylogenetic clade of ribonuclease H (RNase H), furthering its extensive characterization in earlier studies, focusing on E. coli RNase H (ecRNH) and a more stable consensus sequence (AncCcons) differing at 15 positions. The stabilities of 32,768 chimeras between these two sequences were computed using the MSλD framework. The most stable and least stable chimeras were predicted and tested along with several other sequences, revealing a designed chimera with approximately the same stability increase as AncCcons, but requiring only half the mutations. Comparing the computed stabilities with experiment for 12 sequences reveals a Pearson correlation of 0.86 and root mean squared error of 1.18 kcal/mol, an unprecedented level of accuracy well beyond less rigorous computational design methods. We then quantified selection pressure using a simple evolutionary model in which sequences are selected according to the Boltzmann factor of their stability. Selection temperatures from 110 to 168 K are estimated in three ways by comparing experimental and computational results to evolutionary models. These estimates indicate selection pressure is high, which has implications for evolutionary dynamics and for the accuracy required for design, and suggests accurate high-throughput computational methods like MSλD may enable more effective protein design.

A helical fulcrum in eIF2B coordinates allosteric regulation of stress signaling.

The integrated stress response (ISR) enables cells to survive a variety of acute stresses, but chronic activation of the ISR underlies age-related diseases. ISR signaling downregulates translation and activates expression of stress-responsive factors that promote return to homeostasis and is initiated by inhibition of the decameric guanine nucleotide exchange factor eIF2B. Conformational and assembly transitions regulate eIF2B activity, but the allosteric mechanisms controlling these dynamic transitions and mediating the therapeutic effects of the small-molecule ISR inhibitor ISRIB are unknown. Using hydrogen-deuterium exchange-mass spectrometry and cryo-electron microscopy, we identified a central α-helix whose orientation allosterically coordinates eIF2B conformation and assembly. Biochemical and cellular signaling assays show that this 'switch-helix' controls eIF2B activity and signaling. In sum, the switch-helix acts as a fulcrum of eIF2B conformational regulation and is a highly conserved actuator of ISR signal transduction. This work uncovers a conserved allosteric mechanism and unlocks new therapeutic possibilities for ISR-linked diseases.

Evolution avoids a pathological stabilizing interaction in the immune protein S100A9.

Stability constrains evolution. While much is known about constraints on destabilizing mutations, less is known about the constraints on stabilizing mutations. We recently identified a mutation in the innate immune protein S100A9 that provides insight into such constraints. When introduced into human S100A9, M63F simultaneously increases the stability of the protein and disrupts its natural ability to activate Toll-like receptor 4. Using chemical denaturation, we found that M63F stabilizes a calcium-bound conformation of hS100A9. We then used NMR to solve the structure of the mutant protein, revealing that the mutation distorts the hydrophobic binding surface of hS100A9, explaining its deleterious effect on function. Hydrogen-deuterium exchange (HDX) experiments revealed stabilization of the region around M63F in the structure, notably Phe37. In the structure of the M63F mutant, the Phe37 and Phe63 sidechains are in contact, plausibly forming an edge-face π-stack. Mutating Phe37 to Leu abolished the stabilizing effect of M63F as probed by both chemical denaturation and HDX. It also restored the biological activity of S100A9 disrupted by M63F. These findings reveal that Phe63 creates a molecular staple with Phe37 that stabilizes a nonfunctional conformation of the protein, thus disrupting function. Using a bioinformatic analysis, we found that S100A9 proteins from different organisms rarely have Phe at both positions 37 and 63, suggesting that avoiding a pathological stabilizing interaction indeed constrains S100A9 evolution. This work highlights an important evolutionary constraint on stabilizing mutations, namely, that they must avoid inappropriately stabilizing nonfunctional protein conformations.

Mechanistic basis for ubiquitin modulation of a protein energy landscape.

Ubiquitin is a common posttranslational modification canonically associated with targeting proteins to the 26S proteasome for degradation and also plays a role in numerous other nondegradative cellular processes. Ubiquitination at certain sites destabilizes the substrate protein, with consequences for proteasomal processing, while ubiquitination at other sites has little energetic effect. How this site specificity-and, by extension, the myriad effects of ubiquitination on substrate proteins-arises remains unknown. Here, we systematically characterize the atomic-level effects of ubiquitination at various sites on a model protein, barstar, using a combination of NMR, hydrogen-deuterium exchange mass spectrometry, and molecular dynamics simulation. We find that, regardless of the site of modification, ubiquitination does not induce large structural rearrangements in the substrate. Destabilizing modifications, however, increase fluctuations from the native state resulting in exposure of the substrate's C terminus. Both of the sites occur in regions of barstar with relatively high conformational flexibility. Nevertheless, destabilization appears to occur through different thermodynamic mechanisms, involving a reduction in entropy in one case and a loss in enthalpy in another. By contrast, ubiquitination at a nondestabilizing site protects the substrate C terminus through intermittent formation of a structural motif with the last three residues of ubiquitin. Thus, the biophysical effects of ubiquitination at a given site depend greatly on local context. Taken together, our results reveal how a single posttranslational modification can generate a broad array of distinct effects, providing a framework to guide the design of proteins and therapeutics with desired degradation and quality control properties.

Single-molecule diffusometry reveals no catalysis-induced diffusion enhancement of alkaline phosphatase as proposed by FCS experiments.

Theoretical and experimental observations that catalysis enhances the diffusion of enzymes have generated exciting implications about nanoscale energy flow, molecular chemotaxis, and self-powered nanomachines. However, contradictory claims on the origin, magnitude, and consequence of this phenomenon continue to arise. To date, experimental observations of catalysis-enhanced enzyme diffusion have relied almost exclusively on fluorescence correlation spectroscopy (FCS), a technique that provides only indirect, ensemble-averaged measurements of diffusion behavior. Here, using an anti-Brownian electrokinetic (ABEL) trap and in-solution single-particle tracking, we show that catalysis does not increase the diffusion of alkaline phosphatase (ALP) at the single-molecule level, in sharp contrast to the ∼20% enhancement seen in parallel FCS experiments using p-nitrophenyl phosphate (pNPP) as substrate. Combining comprehensive FCS controls, ABEL trap, surface-based single-molecule fluorescence, and Monte Carlo simulations, we establish that pNPP-induced dye blinking at the ∼10-ms timescale is responsible for the apparent diffusion enhancement seen in FCS. Our observations urge a crucial revisit of various experimental findings and theoretical models--including those of our own--in the field, and indicate that in-solution single-particle tracking and ABEL trap are more reliable means to investigate diffusion phenomena at the nanoscale.

Site-specific ubiquitination affects protein energetics and proteasomal degradation.

Changes in the cellular environment modulate protein energy landscapes to drive important biology, with consequences for signaling, allostery and other vital processes. The effects of ubiquitination are particularly important because of their potential influence on degradation by the 26S proteasome. Moreover, proteasomal engagement requires unstructured initiation regions that many known proteasome substrates lack. To assess the energetic effects of ubiquitination and how these manifest at the proteasome, we developed a generalizable strategy to produce isopeptide-linked ubiquitin within structured regions of a protein. The effects on the energy landscape vary from negligible to dramatic, depending on the protein and site of ubiquitination. Ubiquitination at sensitive sites destabilizes the native structure and increases the rate of proteasomal degradation. In well-folded proteins, ubiquitination can even induce the requisite unstructured regions needed for proteasomal engagement. Our results indicate a biophysical role of site-specific ubiquitination as a potential regulatory mechanism for energy-dependent substrate degradation.

Exploring the Evolutionary History of Kinetic Stability in the α-Lytic Protease Family.

In addition to encoding the tertiary fold and stability, the primary sequence of a protein encodes the folding trajectory and kinetic barriers that determine the speed of folding. How these kinetic barriers are encoded is not well understood. Here, we use evolutionary sequence variation in the α-lytic protease (αLP) protein family to probe the relationship between sequence and energy landscape. αLP has an unusual energy landscape: the native state of αLP is not the most thermodynamically favored conformation and, instead, remains folded due to a large kinetic barrier preventing unfolding. To fold, αLP utilizes an N-terminal pro region similar in size to the protease itself that functions as a folding catalyst. Once folded, the pro region is removed, and the native state does not unfold on a biologically relevant time scale. Without the pro region, αLP folds on the order of millennia. A phylogenetic search uncovers αLP homologs with a wide range of pro region sizes, including some with no pro region at all. In the resulting phylogenetic tree, these homologs cluster by pro region size. By studying homologs naturally lacking a pro region, we demonstrate they can be thermodynamically stable, fold much faster than αLP, yet retain the same fold as αLP. Key amino acids thought to contribute to αLP's extreme kinetic stability are lost in these homologs, supporting their role in kinetic stability. This study highlights how the entire energy landscape plays an important role in determining the evolutionary pressures on the protein sequence.

The folding and unfolding behavior of ribonuclease H on the ribosome.

The health of a cell depends on accurate translation and proper protein folding, whereas misfolding can lead to aggregation and disease. The first opportunity for a protein to fold occurs during translation, when the ribosome and surrounding environment can affect the nascent chain energy landscape. However, quantifying these environmental effects is challenging because ribosomal proteins and rRNA preclude most spectroscopic measurements of protein energetics. Here, we have applied two gel-based approaches, pulse proteolysis and force-profile analysis, to probe the folding and unfolding pathways of RNase H (RNH) nascent chains stalled on the prokaryotic ribosome in vitro We found that ribosome-stalled RNH has an increased unfolding rate compared with free RNH. Because protein stability is related to the ratio of the unfolding and folding rates, this increase completely accounts for the observed change in protein stability and indicates that the folding rate is unchanged. Using arrest peptide-based force-profile analysis, we assayed the force generated during the folding of RNH on the ribosome. Surprisingly, we found that population of the RNH folding intermediate is required to generate sufficient force to release a stall induced by the SecM stalling sequence and that readthrough of SecM directly correlates with the stability of the RNH folding intermediate. Together, these results imply that the folding pathway of RNH is unchanged on the ribosome. Furthermore, our findings indicate that the ribosome promotes RNH unfolding while the nascent chain is proximal to the ribosome, which may limit the deleterious effects of RNH misfolding and assist in folding fidelity.

The heat released during catalytic turnover enhances the diffusion of an enzyme.

Recent studies have shown that the diffusivity of enzymes increases in a substrate-dependent manner during catalysis. Although this observation has been reported and characterized for several different systems, the precise origin of this phenomenon is unknown. Calorimetric methods are often used to determine enthalpies from enzyme-catalysed reactions and can therefore provide important insight into their reaction mechanisms. The ensemble averages involved in traditional bulk calorimetry cannot probe the transient effects that the energy exchanged in a reaction may have on the catalyst. Here we obtain single-molecule fluorescence correlation spectroscopy data and analyse them within the framework of a stochastic theory to demonstrate a mechanistic link between the enhanced diffusion of a single enzyme molecule and the heat released in the reaction. We propose that the heat released during catalysis generates an asymmetric pressure wave that results in a differential stress at the protein-solvent interface that transiently displaces the centre-of-mass of the enzyme (chemoacoustic effect). This novel perspective on how enzymes respond to the energy released during catalysis suggests a possible effect of the heat of reaction on the structural integrity and internal degrees of freedom of the enzyme.

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.

Evolutionary trend toward kinetic stability in the folding trajectory of RNases H.

Proper folding of proteins is critical to producing the biological machinery essential for cellular function. The rates and energetics of a protein's folding process, which is described by its energy landscape, are encoded in the amino acid sequence. Over the course of evolution, this landscape must be maintained such that the protein folds and remains folded over a biologically relevant time scale. How exactly a protein's energy landscape is maintained or altered throughout evolution is unclear. To study how a protein's energy landscape changed over time, we characterized the folding trajectories of ancestral proteins of the ribonuclease H (RNase H) family using ancestral sequence reconstruction to access the evolutionary history between RNases H from mesophilic and thermophilic bacteria. We found that despite large sequence divergence, the overall folding pathway is conserved over billions of years of evolution. There are robust trends in the rates of protein folding and unfolding; both modern RNases H evolved to be more kinetically stable than their most recent common ancestor. Finally, our study demonstrates how a partially folded intermediate provides a readily adaptable folding landscape by allowing the independent tuning of kinetics and thermodynamics.

Force-dependent switch in protein unfolding pathways and transition-state movements.

Although it is known that single-domain proteins fold and unfold by parallel pathways, demonstration of this expectation has been difficult to establish in experiments. Unfolding rate, [Formula: see text], as a function of force f, obtained in single-molecule pulling experiments on src SH3 domain, exhibits upward curvature on a [Formula: see text] plot. Similar observations were reported for other proteins for the unfolding rate [Formula: see text]. These findings imply unfolding in these single-domain proteins involves a switch in the pathway as f or [Formula: see text] is increased from a low to a high value. We provide a unified theory demonstrating that if [Formula: see text] as a function of a perturbation (f or [Formula: see text]) exhibits upward curvature then the underlying energy landscape must be strongly multidimensional. Using molecular simulations we provide a structural basis for the switch in the pathways and dramatic shifts in the transition-state ensemble (TSE) in src SH3 domain as f is increased. We show that a single-point mutation shifts the upward curvature in [Formula: see text] to a lower force, thus establishing the malleability of the underlying folding landscape. Our theory, applicable to any perturbation that affects the free energy of the protein linearly, readily explains movement in the TSE in a ß-sandwich (I27) protein and single-chain monellin as the denaturant concentration is varied. We predict that in the force range accessible in laser optical tweezer experiments there should be a switch in the unfolding pathways in I27 or its mutants.

Quantitative determination of ribosome nascent chain stability.

Accurate protein folding is essential for proper cellular and organismal function. In the cell, protein folding is carefully regulated; changes in folding homeostasis (proteostasis) can disrupt many cellular processes and have been implicated in various neurodegenerative diseases and other pathologies. For many proteins, the initial folding process begins during translation while the protein is still tethered to the ribosome; however, most biophysical studies of a protein's energy landscape are carried out in isolation under idealized, dilute conditions and may not accurately report on the energy landscape in vivo. Thus, the energy landscape of ribosome nascent chains and the effect of the tethered ribosome on nascent chain folding remain unclear. Here we have developed a general assay for quantitatively measuring the folding stability of ribosome nascent chains, and find that the ribosome exerts a destabilizing effect on the polypeptide chain. This destabilization decreases as a function of the distance away from the peptidyl transferase center. Thus, the ribosome may add an additional layer of robustness to the protein-folding process by avoiding the formation of stable partially folded states before the protein has completely emerged from the ribosome.

An in vitro tag-and-modify protein sample generation method for single-molecule fluorescence resonance energy transfer.

Biomolecular systems exhibit many dynamic and biologically relevant properties, such as conformational fluctuations, multistep catalysis, transient interactions, folding, and allosteric structural transitions. These properties are challenging to detect and engineer using standard ensemble-based techniques. To address this drawback, single-molecule methods offer a way to access conformational distributions, transient states, and asynchronous dynamics inaccessible to these standard techniques. Fluorescence-based single-molecule approaches are parallelizable and compatible with multiplexed detection; to date, however, they have remained limited to serial screens of small protein libraries. This stems from the current absence of methods for generating either individual dual-labeled protein samples at high throughputs or protein libraries compatible with multiplexed screening platforms. Here, we demonstrate that by combining purified and reconstituted in vitro translation, quantitative unnatural amino acid incorporation via AUG codon reassignment, and copper-catalyzed azide-alkyne cycloaddition, we can overcome these challenges for target proteins that are, or can be, methionine-depleted. We present an in vitro parallelizable approach that does not require laborious target-specific purification to generate dual-labeled proteins and ribosome-nascent chain libraries suitable for single-molecule FRET-based conformational phenotyping. We demonstrate the power of this approach by tracking the effects of mutations, C-terminal extensions, and ribosomal tethering on the structure and stability of three protein model systems: barnase, spectrin, and T4 lysozyme. Importantly, dual-labeled ribosome-nascent chain libraries enable single-molecule co-localization of genotypes with phenotypes, are well suited for multiplexed single-molecule screening of protein libraries, and should enable the in vitro directed evolution of proteins with designer single-molecule conformational phenotypes of interest.

The burst-phase folding intermediate of ribonuclease H changes conformation over evolutionary history.

The amino acid sequence encodes the energy landscape of a protein. Therefore, we expect evolutionary mutations to change features of the protein energy landscape, including the conformations adopted by a polypeptide as it folds to its native state. Ribonucleases H (RNase H) from Escherichia coli and Thermus thermophilus both fold via a partially folded intermediate in which the core region of the protein (helices A-D and strands 4-5) is structured. Strand 1, however, uniquely contributes to the T. thermophilus RNase H folding intermediate (Icore+1 ), but not the E. coli RNase H intermediate (Icore ) (Rosen & Marqusee, PLoS One 2015). We explore the origin of this difference by characterizing the folding intermediate of seven ancestral RNases H spanning the evolutionary history of these two homologs. Using fragment models with or without strand 1 and FRET probes to characterize the folding intermediate of each ancestor, we find a distinct evolutionary trend across the family-the involvement of strand 1 in the folding intermediate is an ancestral feature that is maintained in the thermophilic lineage and is gradually lost in the mesophilic lineage. Evolutionary sequence changes indeed modulate the conformations present on the folding landscape and altered the folding trajectory of RNase H.

Discovery of multiple hidden allosteric sites by combining Markov state models and experiments.

The discovery of drug-like molecules that bind pockets in proteins that are not present in crystallographic structures yet exert allosteric control over activity has generated great interest in designing pharmaceuticals that exploit allosteric effects. However, there have only been a small number of successes, so the therapeutic potential of these pockets--called hidden allosteric sites--remains unclear. One challenge for assessing their utility is that rational drug design approaches require foreknowledge of the target site, but most hidden allosteric sites are only discovered when a small molecule is found to stabilize them. We present a means of decoupling the identification of hidden allosteric sites from the discovery of drugs that bind them by drawing on new developments in Markov state modeling that provide unprecedented access to microsecond- to millisecond-timescale fluctuations of a protein's structure. Visualizing these fluctuations allows us to identify potential hidden allosteric sites, which we then test via thiol labeling experiments. Application of these methods reveals multiple hidden allosteric sites in an important antibiotic target--TEM-1 ß-lactamase. This result supports the hypothesis that there are many as yet undiscovered hidden allosteric sites and suggests our methodology can identify such sites, providing a starting point for future drug design efforts. More generally, our results demonstrate the power of using Markov state models to guide experiments.

Thermodynamic system drift in protein evolution.

Proteins from thermophiles are generally more thermostable than their mesophilic homologs, but little is known about the evolutionary process driving these differences. Here we attempt to understand how the diverse thermostabilities of bacterial ribonuclease H1 (RNH) proteins evolved. RNH proteins from Thermus thermophilus (ttRNH) and Escherichia coli (ecRNH) share similar structures but differ in melting temperature (T(m)) by 20 °C. ttRNH's greater stability is caused in part by the presence of residual structure in the unfolded state, which results in a low heat capacity of unfolding (ΔCp) relative to ecRNH. We first characterized RNH proteins from a variety of extant bacteria and found that Tm correlates with the species' growth temperatures, consistent with environmental selection for stability. We then used ancestral sequence reconstruction to statistically infer evolutionary intermediates along lineages leading to ecRNH and ttRNH from their common ancestor, which existed approximately 3 billion years ago. Finally, we synthesized and experimentally characterized these intermediates. The shared ancestor has a melting temperature between those of ttRNH and ecRNH; the T(m)s of intermediate ancestors along the ttRNH lineage increased gradually over time, while the ecRNH lineage exhibited an abrupt drop in Tm followed by relatively little change. To determine whether the underlying mechanisms for thermostability correlate with the changes in T(m), we measured the thermodynamic basis for stabilization--ΔCp and other thermodynamic parameters--for each of the ancestors. We observed that, while the T(m) changes smoothly, the mechanistic basis for stability fluctuates over evolutionary time. Thus, even while overall stability appears to be strongly driven by selection, the proteins explored a wide variety of mechanisms of stabilization, a phenomenon we call "thermodynamic system drift." This suggests that even on lineages with strong selection to increase stability, proteins have wide latitude to explore sequence space, generating biophysical diversity and potentially opening new evolutionary pathways.