Structure-based Calculations Research Articles

Structure-Based Correlation of Light-Induced Histidine Reactivity in A Model Protein.

Light is known to induce covalently linked aggregates in proteins. These aggregates can be immunogenic and are of concern for drug product development in the biotechnology industry. Histidine (His) is proposed to be a key residue in cross-link generation ( Pattison , D. I. Photochem. Photobiol. Sci. 2012 , 11 , 38 - 53 ). However, the factors that influence the reactivity of His in proteins, especially the intrinsic factors are little known. Here, we used rhDNase, which only forms His-His covalent dimers after light treatment to determine the factors that influence the light-induced reactivity of His. This system allowed us to fully characterize the light-induced covalent dimer and rank the reactivities of the His residues in this protein. The reactivities of these His residues were correlated with solvent accessibility-related parameters both by crystal structure-based calculations of solvent-accessible surface area and by hydrogen-deuterium exchange (HDX) experiments. Through this correlation, we demonstrate that the photoreactivity of His is determined by both solvent accessibility and structural flexibility. This new insight can explain the highly complex chemistry of light-induced aggregation and help predict the aggregation propensity of protein under light treatment.

Elucidating the Structures of the Low- and High-pH Mo(V) Species in Respiratory Nitrate Reductase: A Combined EPR, 14,15N HYSCORE, and DFT Study.

Respiratory nitrate reductases (Nars), members of the prokaryotic Mo/W-bis Pyranopterin Guanosine dinucleotide (Mo/W-bisPGD) enzyme superfamily, are key players in nitrate respiration, a major bioenergetic pathway widely used by microorganisms to cope with the absence of dioxygen. The two-electron reduction of nitrate to nitrite takes place at their active site, where the molybdenum ion cycles between Mo(VI) and Mo(IV) states via a Mo(V) intermediate. The active site shows two distinct pH-dependent Mo(V) electron paramagnetic resonance (EPR) signals whose structure and catalytic relevance have long been debated. In this study, we use EPR and HYSCORE techniques to probe their nuclear environment in Escherichia coli Nar (EcNar). By using samples prepared at different pH and through different enrichment strategies in 98Mo and 15N nuclei, we demonstrate that each of the two Mo(V) species is coupled to a single nitrogen nucleus with similar quadrupole characteristics. Structure-based density functional theory calculations allow us to propose a molecular model of the low-pH Mo(V) species consistent with EPR spectroscopic data. Our results show that the metal ion is coordinated by a monodentate aspartate ligand and permit the assignment of the coupled nitrogen nuclei to the Nδ of Asn52, a residue located ∼3.9 Å to the Mo atom in the crystal structures. This is confirmed by measurements on selectively 15N-Asn labeled EcNar. Further, we propose a Mo-O(H)···HN structure to account for the transfer of spin density onto the interacting nitrogen nucleus deduced from HYSCORE analysis. This work provides a foundation for monitoring the structure of the molybdenum active site in the presence of various substrates or inhibitors in Nars and other molybdenum enzymes.

In Situ Characterization of Protein Adsorption onto Nanoparticles by Fluorescence Correlation Spectroscopy.

Nanotechnology holds great promise for applications in many fields including biology and medicine. Unfortunately, the processes occurring at the interface between nanomaterials and living systems are exceedingly complex and not yet well understood, which has significantly hampered the realization of many nanobiotechnology applications. Whenever nanoparticles (NPs) are incorporated by a living organism, a protein adsorption layer, also known as the "protein corona", forms on the NP surface. Accordingly, living organisms interact with protein-coated rather than bare NPs, and their biological responses depend on the nature of the protein corona. In recent years, a wide variety of biophysical techniques have been employed to elucidate mechanistic aspects of NP-protein interactions. In most studies, NPs are immersed in protein or biofluid (e.g., blood serum) solutions and then separated from the liquid for analysis. Because this approach may modify the composition and structure of the protein corona, our group has pioneered the use of fluorescence correlation spectroscopy (FCS) as an in situ technique, capable of examining NP-protein interactions while the NPs are suspended in biological fluids. FCS allows us to measure, with subnanometer precision and as a function of protein concentration, the increase in hydrodynamic radius of the NPs due to protein adsorption. This Account aims at reviewing recent progress in the exploration of NP-protein interactions by using FCS. In vitro FCS studies of the adsorption of important serum proteins onto water-solubilized luminescent NPs always showed a stepwise increase of the NP radius upon protein binding in the form of a binding isotherm, regardless of the type of NP and its specific surface functionalization. This observation indicates formation of a protein monolayer on the NP. Structure-based calculations of protein surface potentials revealed that positively charged patches on the proteins interact electrostatically with negatively charged NP surfaces, and the observed protein layer thickness always matched the known molecular dimensions of the proteins binding in certain orientations. Temperature and NP surface functionalization have also been identified as important parameters controlling protein corona formation. Notably, while the corona formed from a single type of serum protein was reversible, protein adsorption from complex biological media such as blood serum was entirely irreversible. These quantitative in vitro studies are of great relevance to the bio-nano community and especially to researchers developing engineered nanomaterials for biological and biomedical applications. Future efforts will be directed toward elucidating kinetic aspects of protein corona formation and the detailed structure of the adsorbed proteins at the molecular level. To better appreciate the biological responses triggered by NP exposure, more efforts will be devoted to the exploration of the biomolecular corona as it forms on NPs in contact with living cells, tissues, and even entire model organisms. These studies are challenging when performed in a well-controlled and quantitative fashion and rely on the availability of sophisticated analytical tools, particularly, quantitative optical imaging techniques including FCS and related fluctuation methods.

Structural basis of reversine selectivity in inhibiting Mps1 more potently than aurora B kinase.

Monopolar spindle 1 (Mps1, also known as TTK) is a protein kinase crucial for ensuring that cell division progresses to anaphase only after all chromosomes are connected to spindle microtubules. Incomplete chromosomal attachment leads to abnormal chromosome counts in the daughter cells (aneuploidy), a condition common in many solid cancers. Therefore Mps1 is an established target in cancer therapy. Mps1 kinase inhibitors include reversine (2-(4-morpholinoanilino)-6-cyclohexylaminopurine), a promiscuous compound first recognized as an inhibitor of the Aurora B mitotic kinase. Here, we present the 3.0-Å resolution crystal structure of the Mps1 kinase domain bound to reversine. Structural comparison of reversine bound to Mps1 and Aurora B, indicates a similar binding pose for the purine moiety of reversine making three conserved hydrogen bonds to the protein main chain, explaining the observed promiscuity of this inhibitor. The cyclohexyl and morpholinoaniline moieties of reversine however, have more extensive contacts with the protein in Mps1 than in Aurora B. This is reflected both in structure-based docking energy calculations, and in new experimental data we present here, that both confirm that the affinity of reversine towards Mps1 is about two orders of magnitude higher than towards Aurora B. Thus, our data provides detailed structural understanding of the existing literature that argues reversine inhibits Mps1 more efficiently than Aurora B based on biochemical and in-cell assays. Proteins 2016; 84:1761-1766. © 2016 Wiley Periodicals, Inc.

Spectral Hole Burning in Cyanobacterial Photosystem I with P700 in Oxidized and Neutral States.

We explored the rich satellite hole structures emerging as a result of spectral hole burning in cyanobacterial photosystem I (PSI) and demonstrated that hole burning properties of PSI, particularly at high resolution, are strongly affected by the oxidation state of the primary donor P700, as P700+ effectively quenches the excitations of the lowest-energy antenna states responsible for fluorescence. Obtaining better control of this variable will be crucial for high-resolution ensemble experiments on protein energy landscapes in PSI. The separate nonphotochemical spectral hole burning (NPHB) signatures of various red antenna states were obtained, allowing for additional constraints on excitonic structure-based calculations. Preliminary evidence is presented for an additional red state of PSI of T. elongatus peaked at 712.6 nm, distinct from previously reported C708 and C715 states and possibly involving chlorophyll B15. Excitation at wavelengths as long as 800 nm results in charge separation at cryogenic temperatures in PSI also in Synechocystis sp. PCC 6803. Both the "P700+ minus P700" holes and nonphotochemical spectral holes were subjected to thermocycling. The distribution of barriers manifesting in recovery of the "P700+ minus P700" signature contains two components in sample-dependent proportions, likely reflecting the percentages of FA and FB clusters being successfully prereduced before the optical experiment. The barrier distribution for the recovery of the lower-energy nonphotochemical spectral holes resembles those observed for other pigment-protein complexes, suggesting similar structural elements are responsible for NPHB. Higher-energy components exhibit evidence of "domino effects" such as shifts of certain bands persisting past the lower-energy hole recovery. Thus, conformational changes triggered by excitation of one pigment likely can affect multiple pigments in this tightly packed system.

Binding Mechanism of the N-Terminal SH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The N-terminal Src homology 3 (nSH3) domain of a signaling adaptor protein, CT-10 regulator of kinase II (CrkII), recognizes proline-rich motifs (PRMs) of binding partners, such as cAbl kinase. The interaction between CrkII and cAbl kinase is involved in the regulation of cell spreading, microbial pathogenesis, and cancer metastasis. Here, we report the detailed biophysical characterizations of the interactions between the nSH3 domain of CrkII and PRMs in cAbl. We identified that the nSH3 domain of CrkII binds to three PRMs in cAbl with virtually identical affinities. Structural studies, by using x-ray crystallography and NMR spectroscopy, revealed that the binding modes of all three nSH3:PRM complexes are highly similar to each other. Van ’t Hoff analysis revealed that nSH3:PRM interaction is associated with favorable enthalpy and unfavorable entropy change. The combination of experimentally determined thermodynamic parameters, structure-based calculations, and 15N NMR relaxation analysis highlights the energetic contribution of conformational entropy change upon the complex formation, and water molecules structured in the binding interface of the nSH3:PRM complex. Understanding the molecular basis of nSH3:PRM interaction will provide, to our knowledge, new insights for the rational design of small molecules targeting the interaction between CrkII and cAbl.

Crowding Activates Heat Shock Protein 90

Hsp90 is a dimeric ATP-dependent chaperone involved in the folding, maturation, and activation of diverse target proteins. Extensive in vitro structural analysis has led to a working model of Hsp90's ATP-driven conformational cycle. An implicit assumption is that dilute experimental conditions do not significantly perturb Hsp90 structure and function. However, Hsp90 undergoes a dramatic open/closed conformational change, which raises the possibility that this assumption may not be valid for this chaperone. Indeed, here we show that the ATPase activity of Hsp90 is highly sensitive to molecular crowding, whereas the ATPase activities of Hsp60 and Hsp70 chaperones are insensitive to crowding conditions. Polymer crowders activate Hsp90 in a non-saturable manner, with increasing efficacy at increasing concentration. Crowders exhibit a non-linear relationship between their radius of gyration and the extent to which they activate Hsp90. This experimental relationship can be qualitatively recapitulated with simple structure-based volume calculations comparing open/closed configurations of Hsp90. Thermodynamic analysis indicates that crowding activation of Hsp90 is entropically driven, which is consistent with a model in which excluded volume provides a driving force that favors the closed active state of Hsp90. Multiple Hsp90 homologs are activated by crowders, with the endoplasmic reticulum-specific Hsp90, Grp94, exhibiting the highest sensitivity. Finally, we find that crowding activation works by a different mechanism than co-chaperone activation and that these mechanisms are independent. We hypothesize that Hsp90 has a higher intrinsic activity in the cell than in vitro.

A New Tool for NMR Crystallography: Complete (13)C/(15)N Assignment of Organic Molecules at Natural Isotopic Abundance Using DNP-Enhanced Solid-State NMR.

NMR crystallography of organic molecules at natural isotopic abundance (NA) strongly relies on the comparison of assigned experimental and computed NMR chemical shifts. However, a broad applicability of this approach is often hampered by the still limited (1)H resolution and/or difficulties in assigning (13)C and (15)N resonances without the use of structure-based chemical shift calculations. As shown here, such difficulties can be overcome by (13)C-(13)C and for the first time (15)N-(13)C correlation experiments, recorded with the help of dynamic nuclear polarization. We present the complete de novo (13)C and (15)N resonance assignment at NA of a self-assembled 2'-deoxyguanosine derivative presenting two different molecules in the asymmetric crystallographic unit cell. This de novo assignment method is exclusively based on aforementioned correlation spectra and is an important addition to the NMR crystallography approach, rendering firstly (1)H assignment straightforward, and being secondly a prerequisite for distance measurements with solid-state NMR.

Conformational Reorganization Coupled to the Ionization of Internal Lys Residues in Proteins.

Ionizable groups buried in the hydrophobic interior of proteins are essential for energy transduction and catalysis. Because the protein interior is usually neither as polar nor as polarizable as water, these groups tend to have anomalous pKa values, and their ionization tends to be coupled to conformational reorganization. To elucidate mechanisms of energy transduction in proteins, it is necessary to understand the structural determinants of the pKa values of these buried groups, including the range and character of the conformational reorganization that the ionization of these buried groups can elicit. The L25K and L125K variants of staphylococcal nuclease (SNase) were used to characterize the diverse types of structural reorganization that can be promoted by the ionization of buried groups. NMR relaxation dispersion and ZZ-exchange experiments were used to identify the locations and measure the time scales and extent of pH-dependent conformational exchange in these two proteins. The buried Lys-25 and Lys-125 residues titrate with pKa of 6.3 and 6.2, respectively. The L25K protein fluctuates between the native state and an ensemble of locally unfolded states on the 400 μs to 7 ms time scale. On the 100 to 500 ms time scale the native state exchanges with a subglobally unfolded state in which the β-barrel is partially reorganized. The equilibrium between the native state and this alternative state is highly pH dependent; at pH values below the pKa of Lys-25 the state with the partially reorganized β-barrel is the dominant state. In contrast, the L125K protein only exhibited pH-independent fluctuation in the microsecond to millisecond time scale in the region near Lys-125. The study illustrates how diverse and how localized the coupling between conformational reorganization and ionization of buried groups can be. The pH-sensitive exchange between the fully native and subglobally or locally unfolded states in time scales well into hundreds of milliseconds will challenge all computational methods for structure-based calculations of pKa values.

Arginine: Its pKa value revisited.

Using complementary approaches of potentiometry and NMR spectroscopy, we have determined that the equilibrium acid dissociation constant (pKa value) of the arginine guanidinium group is 13.8 ± 0.1. This is substantially higher than that of ∼ 12 often used in structure-based electrostatics calculations and cited in biochemistry textbooks. The revised intrinsic pKa value helps explains why arginine side chains in proteins are always predominantly charged, even at pH values as great as 10. The high pKa value also reinforces the observation that arginine side chains are invariably protonated under physiological conditions of near neutral pH. This occurs even when the guanidinium moiety is buried in a hydrophobic micro-environment, such as that inside a protein or a lipid membrane, thought to be incompatible with the presence of a charged group.

PKA Values of Buried Groups in Proteins are Sensitive to the Global Thermodynamic Stability

Ionizable groups buried in hydrophobic environments in proteins are essential for all forms of biological energy transduction. The molecular determinants of the pKa values of these internal groups are poorly understood. It is increasingly apparent that conformational reorganization coupled to the ionization of the buried group is a major determinant of these pKa values. Specifically, the creation of charge in hydrophobic environments can trigger a shift from the fully folded state to local or partially unfolded states in which the charge can gain access to water or to an environment where the charge can be solvated. These alternative conformational states are not normally populated owing to the large free energy difference between the alternative and fully-folded native states; however, the partially unfolded states can become the new ground state under pH conditions where the internal group is charged. If the ionization of an internal group promotes the transition to a new conformational state then its pKa should be sensitive to the global thermodynamic stability (ΔG°) of the protein because this determines the energy gap between the ground and the alternative states. This was tested by measuring the pKa of two internal Lys residues in variants of staphylococcal nuclease with thermodynamic stabilities ranging from 8.4 to 13.8 kcal/mol. The magnitude of the shift in the pKa of the internal Lys residues was found to be sensitive to the ΔG° of the protein confirming that the pKa values of these Lys residues are determined by the probability of structural reorganization more than by local dielectric properties of their microenvironments. These observations imply that structure-based pKa calculations for buried groups and other electrostatic processes in hydrophobic environments require accurate treatment of conformational reorganization, which remains an extremely challenging proposition.

Light-Induced Switching of HAMP Domain Conformation and Dynamics Revealed by Time-Resolved EPR Spectroscopy

In microbial photo- and chemotaxis a two-component signaling cascade mediates a regulated response of the flagellar motor to environmental conditions. Upon activation, photo- and chemoreceptors transfer a signal across the plasma membrane to activate the histidine kinase CheA. Successive regulation of the CheY-phosphorylation level controls the flagellar motor. In Natronomonas pharaonis a sensory rhodopsin II - transducer complex (SRII/HtrII) mediates negative phototaxis. As the initial signal, a light-induced outward movement of receptor helix F leads to a conformational change of transducer helix TM2, which in turn propagates the signal to the adjacent HAMP domain. For the HAMP domain, a widely abundant signaling module, several mechanisms were suggested, all comprising two distinct conformational states which we previously observed by two-component cw-EPR spectra at ambient temperatures. Here, we trace the conformational signal and it's propagation throughout the elongated transducer.(1) We applied cw- and pulse-EPR spectroscopy in conjunction with nitroxide spin labeling. We follow transient changes by time-resolved cw-EPR spectroscopy and compare the resulting spectral changes to simulated EPR difference spectra revealing a shift in the thermodynamic equilibrium between the two states. Structure-based calculations of the expected spectral differences shows agreement with a shift towards a more compact state of the HAMP domain. To extend the current signaling models to the whole complex, a trimer of NpSRII/NpHtrII dimers, we carried out molecular dynamics simulations and observed differences between the deactivated and activated complex, ultimately leading to a signaling model that can now be tested experimentally. [1] Klose, D. et al., FEBS Lett. (2014) http://dx.doi.org/10.1016/j.febslet.2014.09.012 (in press)

Electric fingerprint of voltage sensor domains.

A dynamic transmembrane voltage field has been suggested as an intrinsic element in voltage sensor (VS) domains. Here, the dynamic field contribution to the VS energetics was analyzed via electrostatic calculations applied to a number of atomistic structures made available recently. We find that the field is largely static along with the molecular motions of the domain, and more importantly, it is minimally modified across VS variants. This finding implies that sensor domains transfer approximately the same amount of gating charges when moving the electrically charged S4 helix between fixed microscopic configurations. Remarkably, the result means that the observed operational diversity of the domain, including the extension, rate, and voltage dependence of the S4 motion, as dictated by the free energy landscape theory, must be rationalized in terms of dominant variations of its chemical free energy.

Light-induced switching of HAMP domain conformation and dynamics revealed by time-resolved EPR spectroscopy

HAMP domains are widely abundant signaling modules. The putative mechanism of their function comprises switching between two distinct states. To unravel these conformational transitions, we apply site-directed spin labeling and time-resolved EPR spectroscopy to the phototactic receptor/transducer complex NpSRII/NpHtrII. We characterize the kinetic coupling of NpHtrII to NpSRII along with the activation period of the transducer and follow the transient conformational signal. The observed transient shift towards a more compact state of the HAMP domain upon light-activation agrees with structure-based calculations. It thereby validates the two modeled signaling states and integrates the domain’s dynamics into the current model.

PH dependence of conformational fluctuations of the protein backbone.

Proton binding equilibria (pK(a) values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK(a) values measured for individual ionizable residues with NMR spectroscopy are actually population-weighted averages of the pK(a) in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK(a) values of surface Glu and Asp residues are affected by pH-sensitive fluctuations of the backbone between folded and locally unfolded conformations. (15)N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps-ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub-global unfolding near carboxylic groups. This study has implications for structure-based pKa calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH-sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co-operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK(a) values of ionizable residues in proteins.

Predicting proton titration in cationic micelle and bilayer environments.

Knowledge of the protonation behavior of pH-sensitive molecules in micelles and bilayers has significant implications in consumer product development and biomedical applications. However, the calculation of pKa's in such environments proves challenging using traditional structure-based calculations. Here we apply all-atom constant pH molecular dynamics with explicit ions and titratable water to calculate the pKa of a fatty acid molecule in a micelle of dodecyl trimethylammonium chloride and liquid as well as gel-phase bilayers of diethyl ester dimethylammonium chloride. Interestingly, the pKa of the fatty acid in the gel bilayer is 5.4, 0.4 units lower than that in the analogous liquid bilayer or micelle, despite the fact that the protonated carboxylic group is significantly more desolvated in the gel bilayer. This work illustrates the capability of all-atom constant pH molecular dynamics in capturing the delicate balance in the free energies of desolvation and Coulombic interactions. It also shows the importance of the explicit treatment of ions in sampling the protonation states. The ability to model dynamics of pH-responsive substrates in a bilayer environment is useful for improving fabric care products as well as our understanding of the side effects of anti-inflammatory drugs.

Protein Dielectric Constants Determined from NMR Chemical Shift Perturbations

Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.

Bayesian model aggregation for ensemble‐based estimates of protein pKavalues

This article investigates an ensemble-based technique called Bayesian Model Averaging (BMA) to improve the performance of protein amino acid pKa predictions. Structure-based pKa calculations play an important role in the mechanistic interpretation of protein structure and are also used to determine a wide range of protein properties. A diverse set of methods currently exist for pKa prediction, ranging from empirical statistical models to ab initio quantum mechanical approaches. However, each of these methods are based on a set of conceptual assumptions that can effect a model's accuracy and generalizability for pKa prediction in complicated biomolecular systems. We use BMA to combine eleven diverse prediction methods that each estimate pKa values of amino acids in staphylococcal nuclease. These methods are based on work conducted for the pKa Cooperative and the pKa measurements are based on experimental work conducted by the García-Moreno lab. Our cross-validation study demonstrates that the aggregated estimate obtained from BMA outperforms all individual prediction methods with improvements ranging from 45 to 73% over other method classes. This study also compares BMA's predictive performance to other ensemble-based techniques and demonstrates that BMA can outperform these approaches with improvements ranging from 27 to 60%. This work illustrates a new possible mechanism for improving the accuracy of pKa prediction and lays the foundation for future work on aggregate models that balance computational cost with prediction accuracy.

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.

Coupled effect of salt and pH on proteins probed with NMR spectroscopy

The coupled effect of ionic strength (50–400mM) and pH (2–8) on ionization and conformation equilibria of lysozyme was studied using NMR spectroscopy. Observed changes in pKa values of the ionizable groups were found to originate from perturbations in the geometry of hydrogen bonds rather than screening of electric fields. Moreover, at the ionic strengths used here, salt-induced local conformational changes had a dominant effect on chemical shifts measured on 1HN and 15N amide nuclei. Accurate modeling of these localized perturbations in structure-based energy calculations is a necessary prerequisite on the way to complete understanding of any salt-induced processes in proteins.