A miniaturized technique for assessing protein thermodynamics and function using fast determination of quantitative cysteine reactivity

ThermoBRET: A Ligand-Engagement Nanoscale Thermostability Assay Applied to GPCRs

Measurements of membrane protein thermostability reflect ligand binding. Current thermostability assays often require protein purification or rely on pre-existing radiolabelled or fluorescent ligands, limiting their application to established targets. Alternative methods, such as fluorescence-detection size exclusion chromatography thermal shift, detect protein aggregation but are not amenable to high-throughput screening. Here, we present a ThermoBRET method to quantify the relative thermostability of G protein coupled receptors (GPCRs), using cannabinoid receptors (CB1 and CB2 ) and the β2 -adrenoceptor (β2 AR) as model systems. ThermoBRET reports receptor unfolding, does not need labelled ligands and can be used with non-purified proteins. It uses Bioluminescence Resonance Energy Transfer (BRET) between Nanoluciferase (Nluc) and a thiol-reactive fluorescent dye that binds cysteines exposed by unfolding. We demonstrate that the melting point (Tm ) of Nluc-fused GPCRs can be determined in non-purified detergent solubilised membrane preparations or solubilised whole cells, revealing differences in thermostability for different solubilising conditions and in the presence of stabilising ligands. We extended the range of the assay by developing the thermostable tsNLuc by incorporating mutations from the fragments of split-Nluc (Tm of 87 °C versus 59 °C). ThermoBRET allows the determination of GPCR thermostability, which is useful for protein purification optimisation and drug discovery screening.

Cys-labeling kinetics of membrane protein GlpG: a role for specific SDS binding and micelle changes?

Empirically, α-helical membrane protein folding stability in surfactant micelles can be tuned by varying the mole fraction MF^SDS of anionic (sodium dodecyl sulfate (SDS)) relative to nonionic (e.g., dodecyl maltoside (DDM)) surfactant, but we lack a satisfying physical explanation of this phenomenon. Cysteine labeling (CL) has thus far only been used to study the topology of membrane proteins, not their stability or folding behavior. Here, we use CL to investigate membrane protein folding in mixed DDM-SDS micelles. Labeling kinetics of the intramembrane protease GlpG are consistent with simple two-state unfolding-and-exchange rates for seven single-Cys GlpG variants over most of the explored MF^SDS range, along with exchange from the native state at low MF^SDS (which inconveniently precludes measurement of unfolding kinetics under native conditions). However, for two mutants, labeling rates decline with MF^SDS at 0-0.2 MF^SDS (i.e., native conditions). Thus, an increase in MF^SDS seems to be a protective factor for these two positions, but not for the five others. We propose different scenarios to explain this and find the most plausible ones to involve preferential binding of SDS monomers to the site of CL (based on computational simulations) along with changes in size and shape of the mixed micelle with changing MF^SDS (based on SAXS studies). These nonlinear impacts on protein stability highlights a multifaceted role for SDS in membrane protein denaturation, involving both direct interactions of monomeric SDS and changes in micelle size and shape along with the general effects on protein stability of changes in micelle composition.

GPCR Solubilization and Quality Control

G protein-coupled receptors (GPCRs) are versatile membrane proteins involved in the regulation of many physiological processes and pathological conditions, making them interesting pharmacological targets. In order to study their structure and function, GPCRs are traditionally extracted from membranes using detergents. However, due to their hydrophobic nature, intrinsic instability in aqueous solutions, and their denaturing effects, the isolation of properly folded and functional GPCRs is not trivial. Therefore, it is of crucial importance to solubilize receptors under mild conditions and control the sample quality subsequently. Here we describe widely used methods for small-scale GPCR solubilization, followed by quality control based on fluorescence size-exclusion chromatography, SDS-PAGE, temperature-induced protein unfolding (CPM dye binding) and fluorescent ligand binding assay. These methods can easily be used to assess the thermostability and functionality of a GPCR sample exposed to different conditions, such as the use of various detergents, addition of lipids and ligands, making them valuable for obtaining an optimal sample quality for structural and functional studies.

High-Throughput Kinetic Analysis for Target-Directed Covalent Ligand Discovery

Cysteine‐reactive small molecules are used as chemical probes of biological systems and as medicines. Identifying high‐quality covalent ligands requires comprehensive kinetic analysis to distinguish selective binders from pan‐reactive compounds. Quantitative irreversible tethering (qIT), a general method for screening cysteine‐reactive small molecules based upon the maximization of kinetic selectivity, is described. This method was applied prospectively to discover covalent fragments that target the clinically important cell cycle regulator Cdk2. Crystal structures of the inhibitor complexes validate the approach and guide further optimization. The power of this technique is highlighted by the identification of a Cdk2‐selective allosteric (type IV) kinase inhibitor whose novel mode‐of‐action could be exploited therapeutically.

Tyrosine phosphorylation switching of a G protein

Heterotrimeric G protein complexes are molecular switches relaying extracellular signals sensed by G protein-coupled receptors (GPCRs) to downstream targets in the cytoplasm, which effect cellular responses. In the plant heterotrimeric GTPase cycle, GTP hydrolysis, rather than nucleotide exchange, is the rate-limiting reaction and is accelerated by a receptor-like regulator of G signaling (RGS) protein. We hypothesized that posttranslational modification of the Gα subunit in the G protein complex regulates the RGS-dependent GTPase cycle. Our structural analyses identified an invariant phosphorylated tyrosine residue (Tyr¹⁶⁶ in the Arabidopsis Gα subunit AtGPA1) located in the intramolecular domain interface where nucleotide binding and hydrolysis occur. We also identified a receptor-like kinase that phosphorylates AtGPA1 in a Tyr¹⁶⁶-dependent manner. Discrete molecular dynamics simulations predicted that phosphorylated Tyr¹⁶⁶ forms a salt bridge in this interface and potentially affects the RGS protein-accelerated GTPase cycle. Using a Tyr¹⁶⁶ phosphomimetic substitution, we found that the cognate RGS protein binds more tightly to the GDP-bound Gα substrate, consequently reducing its ability to accelerate GTPase activity. In conclusion, we propose that phosphorylation of Tyr¹⁶⁶ in AtGPA1 changes the binding pattern with AtRGS1 and thereby attenuates the steady-state rate of the GTPase cycle. We coin this newly identified mechanism "substrate phosphoswitching."

Coordinated regulation of intracellular pH by two glucose-sensing pathways in yeast

The yeast Saccharomyces cerevisiae employs multiple pathways to coordinate sugar availability and metabolism. Glucose and other sugars are detected by a G protein-coupled receptor, Gpr1, as well as a pair of transporter-like proteins, Rgt2 and Snf3. When glucose is limiting, however, an ATP-driven proton pump (Pma1) is inactivated, leading to a marked decrease in cytoplasmic pH. Here we determine the relative contribution of the two sugar-sensing pathways to pH regulation. Whereas cytoplasmic pH is strongly dependent on glucose abundance and is regulated by both glucose-sensing pathways, ATP is largely unaffected and therefore cannot account for the changes in Pma1 activity. These data suggest that the pH is a second messenger of the glucose-sensing pathways. We show further that different sugars differ in their ability to control cellular acidification, in the manner of inverse agonists. We conclude that the sugar-sensing pathways act via Pma1 to invoke coordinated changes in cellular pH and metabolism. More broadly, our findings support the emerging view that cellular systems have evolved the use of pH signals as a means of adapting to environmental stresses such as those caused by hypoxia, ischemia, and diabetes.

Regulation of Ras Paralog Thermostability by Networks of Buried Ionizable Groups

Protein folding is governed by a variety of molecular forces including hydrophobic and ionic interactions. Less is known about the molecular determinants of protein stability. Here we used a recently developed computer algorithm (pHinder) to investigate the relationship between buried charge and thermostability. Our analysis revealed that charge networks in the protein core are generally smaller in thermophilic organisms as compared to mesophilic organisms. To experimentally test whether core network size influences protein thermostability, we purified 18 paralogous Ras superfamily GTPases from yeast and determined their melting temperatures (Tm, or temperature at which 50% of the protein is unfolded). This analysis revealed a wide range of Tm values (35-63 °C) that correlated significantly (R = 0.87) with core network size. These results suggest that thermostability depends in part on the arrangement of ionizable side chains within a protein core. An improved capacity to predict protein thermostability may be useful for selecting the best candidates for protein crystallography, the development of protein-based therapeutics, as well as for industrial enzyme applications.

The experimental power of FR900359 to study Gq-regulated biological processes

Despite the discovery of heterotrimeric αβγ G proteins ∼25 years ago, their selective perturbation by cell-permeable inhibitors remains a fundamental challenge. Here we report that the plant-derived depsipeptide FR900359 (FR) is ideally suited to this task. Using a multifaceted approach we systematically characterize FR as a selective inhibitor of Gq/11/14 over all other mammalian Gα isoforms and elaborate its molecular mechanism of action. We also use FR to investigate whether inhibition of Gq proteins is an effective post-receptor strategy to target oncogenic signalling, using melanoma as a model system. FR suppresses many of the hallmark features that are central to the malignancy of melanoma cells, thereby providing new opportunities for therapeutic intervention. Just as pertussis toxin is used extensively to probe and inhibit the signalling of Gi/o proteins, we anticipate that FR will at least be its equivalent for investigating the biological relevance of Gq.

Redox regulation of Rac1 by thiol oxidation

The Rac1 GTPase is an essential and ubiquitous protein that signals through numerous pathways to control critical cellular processes, including cell growth, morphology, and motility. Rac1 deletion is embryonic lethal, and its dysregulation or mutation can promote cancer, arthritis, cardiovascular disease, and neurological disorders. Rac1 activity is highly regulated by modulatory proteins and posttranslational modifications. Whereas much attention has been devoted to guanine nucleotide exchange factors that act on Rac1 to promote GTP loading and Rac1 activation, cellular oxidants may also regulate Rac1 activation by promoting guanine nucleotide exchange. Herein, we show that Rac1 contains a redox-sensitive cysteine (Cys(18)) that can be selectively oxidized at physiological pH because of its lowered pKa. Consistent with these observations, we show that Rac1 is glutathiolated in primary chondrocytes. Oxidation of Cys(18) by glutathione greatly perturbs Rac1 guanine nucleotide binding and promotes nucleotide exchange. As aspartate substitutions have been previously used to mimic cysteine oxidation, we characterized the biochemical properties of Rac1(C18D). We also evaluated Rac1(C18S) as a redox-insensitive variant and found that it retains structural and biochemical properties similar to those of Rac1(WT) but is resistant to thiol oxidation. In addition, Rac1(C18D), but not Rac1(C18S), shows greatly enhanced nucleotide exchange, similar to that observed for Rac1 oxidation by glutathione. We employed Rac1(C18D) in cell-based studies to assess whether this fast-cycling variant, which mimics Rac1 oxidation by glutathione, affects Rac1 activity and function. Expression of Rac1(C18D) in Swiss 3T3 cells showed greatly enhanced GTP-bound Rac1 relative to Rac1(WT) and the redox-insensitive Rac1(C18S) variant. Moreover, expression of Rac1(C18D) in HEK-293T cells greatly promoted lamellipodia formation. Our results suggest that Rac1 oxidation at Cys(18) is a novel posttranslational modification that upregulates Rac1 activity.

Phosphorylation at Y1065 in vinculin mediates actin bundling, cell spreading, and mechanical responses to force

Vinculin is an essential structural adaptor protein that localizes to sites of adhesion and is involved in a number of cell processes including adhesion, spreading, motility, force transduction, and cell survival. The C-terminal vinculin tail domain (Vt) contains the necessary structural components to bind and cross-link actin filaments. Actin binding to Vt induces a conformational change that promotes dimerization through the C-terminal hairpin of Vt and enables actin filament cross-linking. Here we show that Src phosphorylation of Y1065 within the C-terminal hairpin regulates Vt-mediated actin bundling and provide a detailed characterization of Y1065 mutations. Furthermore, we show that phosphorylation at Y1065 plays a role in cell spreading and the response to the application of mechanical force.

Energetics-based methods for protein folding and stability measurements

Over the past 15 years, a series of energetics-based techniques have been developed for the thermodynamic analysis of protein folding and stability. These techniques include Stability of Unpurified Proteins from Rates of amide H/D Exchange (SUPREX), pulse proteolysis, Stability of Proteins from Rates of Oxidation (SPROX), slow histidine H/D exchange, lysine amidination, and quantitative cysteine reactivity (QCR). The above techniques, which are the subject of this review, all utilize chemical or enzymatic modification reactions to probe the chemical denaturant- or temperature-induced equilibrium unfolding properties of proteins and protein-ligand complexes. They employ various mass spectrometry-, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE)-, and optical spectroscopy-based readouts that are particularly advantageous for high-throughput and in some cases multiplexed analyses. This has created the opportunity to use protein folding and stability measurements in new applications such as in high-throughput screening projects to identify novel protein ligands and in mode-of-action studies to identify protein targets of a particular ligand.

Guanine nucleotide-binding protein (Gα) endocytosis by a cascade of ubiquitin binding domain proteins is required for sustained morphogenesis and proper mating in yeast

Heterotrimeric G proteins are well known to transmit signals from cell surface receptors to intracellular effector proteins. There is growing appreciation that G proteins are also present at endomembrane compartments, where they can potentially interact with a distinct set of signaling proteins. Here, we examine the cellular trafficking function of the G protein α subunit in yeast, Gpa1. Gpa1 contains a unique 109-amino acid insert within the α-helical domain that undergoes a variety of posttranslational modifications. Among these is monoubiquitination, catalyzed by the NEDD4 family ubiquitin ligase Rsp5. Using a newly optimized method for G protein purification together with biophysical measures of structure and function, we show that the ubiquitination domain does not influence enzyme activity. By screening a panel of 39 gene deletion mutants, each lacking a different ubiquitin binding domain protein, we identify seven that are necessary to deliver Gpa1 to the vacuole compartment including four proteins (Ede1, Bul1, Ddi1, and Rup1) previously not known to be involved in this process. Finally, we show that proper endocytosis of the G protein is needed for sustained cellular morphogenesis and mating in response to pheromone stimulation. We conclude that a cascade of ubiquitin-binding proteins serves to deliver the G protein to its final destination within the cell. In this instance and in contrast to the previously characterized visual system, endocytosis from the plasma membrane is needed for proper signal transduction rather than for signal desensitization.

Biophysical and proteomic characterization strategies for cysteine modifications in Ras GTPases

Cysteine is one of the most reactive amino acids and is modified by a number of oxidants. The reactivity of cysteines is dependent on the thiol pK a; however, measuring cysteine pK a values is nontrivial. Ras family GTPases have been shown to contain a free cysteine that is sensitive to oxidation, and free radical-mediated oxidation of this cysteine has been shown to be activating. Here, we present a new technique that allows for measuring cysteine pK a values using a fluorescent detection system with the molecule 4-fluoro-7-aminosulfonylbenzofurazan (ABD-F). In addition, we also describe how to generate several oxidants. Lastly, we describe several mass spectrometry-based experiments and the necessary adjustments to the experiments to detect cysteine oxidation.

Differences in the regulation of K-Ras and H-Ras isoforms by monoubiquitination

Ras GTPases are signaling switches that control critical cellular processes including gene expression, differentiation, and apoptosis. The major Ras isoforms (K, H, and N) contain a conserved core GTPase domain, but have distinct biological functions. Among the three Ras isoforms there are clear differences in post-translational regulation, which contribute to differences in localization and signaling output. Modification by ubiquitination was recently reported to activate Ras signaling in cells, but the mechanisms of activation are not well understood. Here, we show that H-Ras is activated by monoubiquitination and that ubiquitination at Lys-117 accelerates intrinsic nucleotide exchange, thereby promoting GTP loading. This mechanism of Ras activation is distinct from K-Ras monoubiquitination at Lys-147, which leads to impaired regulator-mediated GTP hydrolysis. These findings reveal that different Ras isoforms are monoubiquitinated at distinct sites, with distinct mechanisms of action, but with a common ability to chronically activate the protein in the absence of a receptor signal or oncogenic mutation.

Protons as second messenger regulators of G protein signaling

In response to environmental stress, cells often generate pH signals that serve to protect vital cellular components and reprogram gene expression for survival. A major barrier to our understanding of this process has been the identification of signaling proteins that detect changes in intracellular pH. To identify candidate pH sensors, we developed a computer algorithm that searches proteins for networks of proton-binding sidechains. This analysis indicates that Gα subunits, the principal transducers of G protein-coupled receptor (GPCR) signals, are pH sensors. Our structure-based calculations and biophysical investigations reveal that Gα subunits contain networks of pH-sensing sidechains buried between their Ras and helical domains. Further, we show that proton binding induces changes in conformation that promote Gα phosphorylation and suppress receptor-initiated signaling. Together, our computational, biophysical, and cellular analyses reveal an unexpected function for G proteins as mediators of stress-response signaling.

Expanding the number of 'druggable' targets: non-enzymes and protein-protein interactions

Following sequencing and assembly of the human genome, the preferred methods for identification of new drug targets have changed dramatically. Modern tactics such as genome-wide association studies (GWAS) and deep sequencing are fundamentally different from the pharmacology-guided approaches used previously, in which knowledge of small molecule ligands acting at their cellular targets was the primary discovery engine. A consequence of the 'target-first, pharmacology-second' strategy is that many predicted drug targets are non-enzymes, such as scaffolding, regulatory or structural proteins, and their activities are often dependent on protein-protein interactions (PPIs). These types of targets create unique challenges to drug discovery efforts because enzymatic turnover cannot be used as a convenient surrogate for compound potency. Moreover, it is often challenging to predict how ligand binding to non-enzymes might affect changes in protein function and/or pathobiology. Thus, in the postgenomic era, targets might be strongly implicated by molecular biology-based methods, yet they often later earn the designation of 'undruggable'. Can the scope of available targets be widened to include these promising, but challenging, non-enzymes? In this review, we discuss advances in high-throughput screening (HTS) technology and chemical library design that are emerging to deal with these challenges.

Evolution-based design of proteins

Statistical analysis of protein sequences indicates an architecture for natural proteins in which amino acids are engaged in a sparse, hierarchical pattern of interactions in the tertiary structure. This architecture might be a key and distinguishing feature of evolved proteins-a design principle providing not only for foldability and high-performance function but also for robustness to perturbation and the capacity for rapid adaptation to new selection pressures. Here, we describe an approach for systematically testing this design principle for natural-like proteins by (1) computational design of synthetic sequences that gradually add or remove constraints along the hierarchy of interacting residues and (2) experimental testing of the designed sequences for folding and biochemical function. By this process, we hope to understand how the constraints on fold, function, and other aspects of fitness are organized within natural proteins, a first step in understanding the process of "design" by evolution.

Site-specific monoubiquitination activates Ras by impeding GTPase-activating protein function

Cell growth and differentiation are controlled by growth factor receptors coupled to the GTPase Ras. Oncogenic mutations disrupt GTPase activity leading to persistent Ras signaling and cancer progression. Recent evidence indicates that monoubiquitination of Ras leads to Ras activation. Mutation of the primary site of monoubiquitination impairs the ability of activated K–Ras to promote tumor growth. To determine the mechanism of human Ras activation we chemically ubiquitinated the protein and analyzed its function by NMR, computational modeling, and biochemical activity measurements. We established that monoubiquitination has little effect on Ras GTP binding, GTP hydrolysis, or exchange factor activation, but severely abrogates the response to GTPase activating proteins in a site–specific manner. These findings reveal a new mechanism by which Ras can trigger persistent signaling in the absence of receptor activation or an oncogenic mutation.

Differences in intradomain and interdomain motion confer distinct activation properties to structurally similar Gα proteins

Proteins with similar crystal structures can have dissimilar rates of substrate binding and catalysis. Here we used molecular dynamics simulations and biochemical analysis to determine the role of intradomain and interdomain motions in conferring distinct activation rates to two Gα proteins, Gα(i1) and GPA1. Despite high structural similarity, GPA1 can activate itself without a receptor, whereas Gα(i1) cannot. We found that motions in these proteins vary greatly in type and frequency. Whereas motion is greatest in the Ras domain of Gα(i1), it is greatest in helices αA and αB from the helical domain of GPA1. Using protein chimeras, we show that helix αA from GPA1 is sufficient to confer rapid activation to Gα(i1). Gα(i1) has less intradomain motion than GPA1 and instead displays interdomain displacement resembling that observed in a receptor-heterotrimer crystal complex. Thus, structurally similar proteins can have distinct atomic motions that confer distinct activation mechanisms.

Quantitation of protein-protein interactions by thermal stability shift analysis

Thermal stability shift analysis is a powerful method for examining binding interactions in proteins. We demonstrate that under certain circumstances, protein-protein interactions can be quantitated by monitoring shifts in thermal stability using thermodynamic models and data analysis methods presented in this work. This method relies on the determination of protein stabilities from thermal unfolding experiments using fluorescent dyes such as SYPRO Orange that report on protein denaturation. Data collection is rapid and straightforward using readily available real-time polymerase chain reaction instrumentation. We present an approach for the analysis of the unfolding transitions corresponding to each partner to extract the affinity of the interaction between the proteins. This method does not require the construction of a titration series that brackets the dissociation constant. In thermal shift experiments, protein stability data are obtained at different temperatures according to the affinity- and concentration-dependent shifts in unfolding transition midpoints. Treatment of the temperature dependence of affinity is, therefore, intrinsic to this method and is developed in this study. We used the interaction between maltose-binding protein (MBP) and a thermostable synthetic ankyrin repeat protein (Off7) as an experimental test case because their unfolding transitions overlap minimally. We found that MBP is significantly stabilized by Off7. High experimental throughput is enabled by sample parallelization, and the ability to extract quantitative binding information at a single partner concentration. In a single experiment, we were able to quantify the affinities of a series of alanine mutants, covering a wide range of affinities (∼ 100 nM to ∼ 100 μM).

Protein stability by number: high-throughput and statistical approaches to one of protein science's most difficult problems

Most proteins are only barely stable, which impedes research, complicates therapeutic applications, and makes proteins susceptible to pathologically destabilizing mutations. Our ability to predict the thermodynamic consequences of even single point mutations is still surprisingly limited, and established methods of measuring stability are slow. Recent advances are bringing protein stability studies into the high-throughput realm. Some methods are based on inferential read-outs such as activity, proteolytic resistance or split-protein fragment reassembly. Other methods use miniaturization of direct measurements, such as intrinsic fluorescence, H/D exchange, cysteine reactivity, aggregation and hydrophobic dye binding (DSF). Protein engineering based on statistical analysis (consensus and correlated occurrences of amino acids) is promising, but much work remains to understand and implement these methods.