Sequence-specific 1H and 15N resonance assignments and secondary structure of GDP-bound human c-Ha-Ras protein in solution

Conserved allosteric perturbation of the GTPase domains by region 1 of Ras hypervariable regions

Ras proteins are important intracellular signaling hubs that can interact with numerous downstream effectors and upstream regulators through their GTPase domains (G-domains) anchored to plasma membranes by the C-terminal hypervariable regions (HVRs). The biological functions of Ras were proposed to be regulated at multiple levels including the intramolecular G-domain-HVR interactions, of which the exact mechanism and specificity are still controversial. Here, we demonstrate that the HVRs, instead of having direct contacts, can weakly perturb the G-domains via an allosteric interaction that is restricted to a ∼20 Å range and highly conserved in the tested Ras isoforms (HRas and KRas4B) and nucleotide-bound states. The origin of this allosteric perturbation has been localized to a short segment (residues 167-171) coinciding with region 1 of HVRs, which exhibits moderate to weak α-helical propensities. A charge-reversal mutation (E168K) of KRas4B in region 1, previously described in the Catalog of Somatic Mutations in Cancer database, was found to induce similar chemical shift perturbations as truncation of the HVR does. Further membrane paramagnetic relaxation enhancement (mPRE) data show that this region 1 mutation alters the membrane orientations of KRas4B and moderately increases the relative population of the signaling-compatible state.

1H, 15N backbone assignment and comparative analysis of the wild type and G12C, G12D, G12V mutants of K-Ras bound to GDP at physiological pH

K-Ras protein is a membrane-bound small GTPase acting as a molecular switch. It plays a key role in many signal transduction pathways regulating cell proliferation, differentiation, survival, etc. It alternates between its GTP-bound active and the GDP-bound inactive conformers regulated by guanine nucleotide exchange factors and GTPase activating proteins. Its most frequent oncogenic mutants are G12C, G12D, and G12V that have impaired GTPase activity, thus induce malignant tumors. Here we report the resonance assignment of the backbone ¹H and ¹⁵N nuclei of K-Ras wildtype, G12C, G12D and G12V proteins' catalytic G domain (1-169 residues) in GDP-bound state, and ¹³C of backbone and side chains of G12C mutant at physiological pH 7.4. Triple resonance data were used to get secondary structure information and backbone dynamics of G12C, the best-known drug target among K-Ras mutants. Simultaneous investigation of G12C, G12D and G12V mutants, along with the wild type form at the very same conditions allowed us to perform a comprehensive analysis based on the combined chemical shifts to reveal the effect of mutation at G12 position on structure. Intriguingly, the G12C and G12V mutants found to be structurally very similar at the three most important regions of K-Ras (P-loop, Switch-I, Switch-II), while the G12D mutant significantly differs at P-loop and Switch-II from the wildtype as well as G12C and G12V mutants. However, in Switch-I it hardly deviates from the wildtype protein.

NMR 1H,13C, 15N backbone and 13C side chain resonance assignment of the G12C mutant of human K-Ras bound to GDP

K-Ras is a key driver of oncogenesis, accounting for approximately 80% of Ras-driven human cancers. The small GTPase cycles between an inactive, GDP-bound and an active, GTP-bound state, regulated by guanine nucleotide exchange factors and GTPase activating proteins, respectively. Activated K-Ras regulates cell proliferation, differentiation and survival by signaling through several effector pathways, including Raf-MAPK. Oncogenic mutations that impair the GTPase activity of K-Ras result in a hyperactivated state, leading to uncontrolled cellular proliferation and tumorogenesis. A cysteine mutation at glycine 12 is commonly found in K-Ras associated cancers, and has become a recent focus for therapeutic intervention. We report here ¹H^{N, 15}N, and ¹³C resonance assignments for the 19.3 kDa (aa 1-169) human K-Ras protein harboring an oncogenic G12C mutation in the GDP-bound form (K-RAS^G12C-GDP), using heteronuclear, multidimensional NMR spectroscopy. Backbone ¹H-¹⁵N correlations have been assigned for all non-proline residues, except for the first methionine residue.

NMR Investigation of Structures of G-protein Coupled Receptor Folding Intermediates

Folding of G-protein coupled receptors (GPCRs) according to the two-stage model (Popot, J. L., and Engelman, D. M. (1990) Biochemistry 29, 4031-4037) is postulated to proceed in 2 steps: partitioning of the polypeptide into the membrane followed by diffusion until native contacts are formed. Herein we investigate conformational preferences of fragments of the yeast Ste2p receptor using NMR. Constructs comprising the first, the first two, and the first three transmembrane (TM) segments, as well as a construct comprising TM1-TM2 covalently linked to TM7 were examined. We observed that the isolated TM1 does not form a stable helix nor does it integrate well into the micelle. TM1 is significantly stabilized upon interaction with TM2, forming a helical hairpin reported previously (Neumoin, A., Cohen, L. S., Arshava, B., Tantry, S., Becker, J. M., Zerbe, O., and Naider, F. (2009) Biophys. J. 96, 3187-3196), and in this case the protein integrates into the hydrophobic interior of the micelle. TM123 displays a strong tendency to oligomerize, but hydrogen exchange data reveal that the center of TM3 is solvent exposed. In all GPCRs so-far structurally characterized TM7 forms many contacts with TM1 and TM2. In our study TM127 integrates well into the hydrophobic environment, but TM7 does not stably pack against the remaining helices. Topology mapping in microsomal membranes also indicates that TM1 does not integrate in a membrane-spanning fashion, but that TM12, TM123, and TM127 adopt predominantly native-like topologies. The data from our study would be consistent with the retention of individual helices of incompletely synthesized GPCRs in the vicinity of the translocon until the complete receptor is released into the membrane interior.

The Ras G Domain Lacks the Intrinsic Propensity to Form Dimers

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature.

¹H, ¹³C and ¹⁵N resonance assignment for the human K-Ras at physiological pH

K-Ras, a member of the Ras family of small GTPases, is involved in cell growth, proliferation, differentiation and apoptosis and is frequently mutated in cancer. The activity of Ras is mediated by the inter-conversion between GTP- and GDP- bound states. This conversion is regulated by binding of effector proteins such as guanine nucleotide exchange factors and GTPase activating proteins. Previously, NMR signals from these effector-binding regions of Ras often remained unassigned and largely unobservable due to conformational exchange and polysterism inherent to this protein. In this paper, we report the complete backbone and C(β), as well as partial H(α), H(β) and C(γ), NMR assignment for human K-Ras (residues 1-166) in the GDP-bound form at a physiological pH of 7.4. These data thereby make possible detailed monitoring of the functional cycle of Ras and its interactions with nucleotides and effector proteins through the observation of fingerprint signals from all the functionally important regions of the protein.

NMR characterization of full-length farnesylated and non-farnesylated H-Ras and its implications for Raf activation

The C terminus, also known as the hypervariable region (residues 166-189), of H-, N-, and K-Ras proteins has sequence determinants necessary for full activation of downstream effectors such as Raf kinase and PI-3 kinase as well as for the correct targeting of Ras proteins to lipid rafts and non-raft membranes. There is considerable interest in understanding how residues in the extreme C terminus of the different Ras proteins and farnesylation of the CaaX box cysteine affect Ras membrane localization and allosteric activation of Raf kinase. To provide insights into the structural and dynamic changes that occur in Ras upon farnesylation, we have used NMR spectroscopy to compare the properties of truncated H-Ras (1-166), to non-processed full-length H-Ras (residues 1-185) and full-length (1-189) farnesylated H-Ras. We report that the C-terminal helix alpha-5 extends to residue N172, and the remaining 17 amino acid residues in the C terminus are conformationally averaged in solution. Removal of either 23 or 18 amino acid residues from the C terminus of full length H-Ras generates truncated H-Ras (1-166) and H-Ras (1-171) proteins, respectively, that have been structurally characterized and are biochemical active. Here we report that C-terminal truncation of H-Ras results in minor structural and dynamic perturbations that are propagated throughout the H-Ras protein including increased flexibility of the central beta-sheet and the C-terminal helix alpha-5. Ordering of residues in loop-2, which is involved in Raf CRD binding is also observed. Farnesylation of full-length H-Ras at C186 does not result in detectable conformational changes in H-Ras. Chemical shift mapping studies of farnesylated and non-farnesylated forms of H-Ras with the Raf-CRD show that the farnesyl moiety, the extreme H-Ras C terminus and residues 23-30, contribute to H-Ras:Raf-CRD interactions, thereby increasing the affinity of H-Ras for the Raf-CRD.

Glucosylation of Ras by Clostridium sordellii lethal toxin: consequences for effector loop conformations observed by NMR spectroscopy

The lethal toxin (LT) from Clostridium sordellii, which belongs to the family of large clostridial cytotoxins, acts as a monoglucosyltransferase for the Rho subfamily GTPase Rac and also modifies Ras. In the present study we investigated structural changes of H-Ras in its di- and triphosphate form that occur upon glucosylation of the effector domain amino acid threonine-35 by LT. (31)P NMR experiments recorded during the enzymatic glucosylation process, using UDP-glucose as a cosubstrate, show that the modification of the threonine side chain influences the chemical shifts of the phosphate groups of the bound nucleotides. In the diphosphate-bound form (Ras.GDP) glucosylation of Thr35 induces only small changes in the chemical environment of the active center. In the triphosphate form with the GTP analogue GppNHp bound (Ras.GppNHp) Ras shows at least two different conformations in the active center that exchange on a medium-range time scale (10 to 0.1 ms). Glucosylation selectively stabilizes one distinct conformation of the effector loop (state 1) with tyrosine-32 probably apart from the nucleotide and threonine-35 not involved in magnesium ion coordination. This conformation is known to have a low affinity to effector proteins such as Raf-1, AF-6, or Byr2 and thus prevents the transduction of the activation signal in the Ras-mediated pathway. NMR correlation spectra of Ras(T35glc).GDP and denaturation experiments with urea indicate that the glucose is bound in the alpha-anomeric form to the hydroxyl group of the threonine-35 side chain. Inhibition of the glucosylation reaction by 1,5-gluconolactone suggests a stereospecific reaction mechanism with a glucosyl oxonium ion transition state for the enzymatic activity of LT.

Nuclear magnetic resonance and molecular dynamics studies on the interactions of the Ras-binding domain of Raf-1 with wild-type and mutant Ras proteins

The Ras protein and its homolog, Rap1A, have an identical "effector region" (residues 32-40) preceded by Asp30-Glu31 and Glu30-Lys31, respectively. In the complex of the "Ras-like" E30D/K31E mutant Rap1A with the Ras-binding domain (RBD), residues 51-131 of Raf-1, Glu31 in Rap1A forms a tight salt bridge with Lys84 in Raf-1. However, we have recently found that Raf-1 RBD binding of Ras is indeed reduced by the E31K mutation, but is not affected by the E31A mutation. Here, the "Rap1A-like" D30E/E31K mutant of Ras was prepared and shown to bind the Raf-1 RBD less strongly than wild-type Ras, but slightly more tightly than the E31K mutant. The backbone 1H, 13C, and 15N magnetic resonances of the Raf-1 RBD were assigned in complexes with the wild-type and D30E/E31K mutant Ras proteins in the guanosine 5'-O-(beta,gamma-imidotriphosphate)-bound form. The Lys84 residue in the Raf-1 RBD exhibited a large change in chemical shift upon binding wild-type Ras, suggesting that Lys84 interacts with wild-type Ras. The D30E/E31K mutant of Ras caused nearly the same perturbations in Raf-1 chemical shifts, including that of Lys84. We hypothesized that Glu31 in Ras may not be the major salt bridge partner of Lys84 in Raf-1. A molecular dynamics simulation of a model structure of the Raf-1 RBD.Ras.GTP complex suggested that Lys84 in Raf-1 might instead form a tight salt bridge with Asp33 in Ras. Consistent with this, the D33A mutation in Ras greatly reduced its Raf-I RBD binding activity. We conclude that the major salt bridge partner of Lys84 in Raf-1 may be Asp33 in Ras.

Cell-free production and stable-isotope labeling of milligram quantities of proteins

We have improved the productivity of an Escherichia coli cell-free protein synthesis system. First, creatine phosphate and creatine kinase were used as the energy source regeneration system, and the other components of the reaction mixture were optimized. Second, the E. coli S30 cell extract was condensed by dialysis against a polyethylene glycol solution to increase the rate of synthesis. Third, during the protein synthesis, the reaction mixture was dialyzed against a low-molecular-weight substrate solution to prolong the reaction. Thus, the yield of chloramphenicol acetyltransferase was raised to 6 mg/ml of reaction mixture. Stable-isotope labeling of a protein with 13C/15N-labeled amino acids for NMR spectroscopy was achieved by this method.

Solution structure of the Ras-binding domain of RGL

The RGL protein, a homolog of the Ral GDP dissociation stimulator (RalGDS), has been identified as a downstream effector of Ras. In the present study, the solution structure of the Ras-binding domain of RGL (RGL-RBD) was determined by NMR spectroscopy. The overall fold of RGL-RBD consists of a five-stranded beta-sheet and two alpha-helices, which is the same topology as that of RalGDS-RBD. The backbone chemical shift perturbation of RGL-RBD upon interaction with the GTP analog-bound Ras was also examined. The solution structure of RGL-RBD, especially around some of the Ras-interacting residues, is appreciably different from that of RalGDS-RBD.

Regional polysterism in the GTP-bound form of the human c-Ha-Ras protein

The backbone 1H, 13C, and 15N resonances of the c-Ha-Ras protein [a truncated version consisting of residues 1-171, Ras(1-171)] bound with GMPPNP (a slowly hydrolyzable analogue of GTP) were assigned and compared with those of the GDP-bound Ras(1-171). The backbone amide resonances of amino acid residues 10-13, 21, 31-39, 57-64, and 71 of Ras(1-171).GMPPNP, but not those of Ras(1-171).GDP, were extremely broadened, whereas other residues of Ras(1-171).GMPPNP exhibited amide resonances nearly as sharp as those of Ras(1-171). GDP. The residues exhibiting the extreme broadening, except for residues 21 and 71, are localized in three functional loop regions [loops L1, L2 (switch I), and L4 (switch II)], which are involved in hydrolysis of GTP and interactions with other proteins. From the temperature and magnetic field strength dependencies of the backbone amide resonance intensities, the extreme broadening was ascribed to the exchange at an intermediate rate on the NMR time scale. It was shown that the Ras(1-171) protein bound with GTP or GTPgammaS (another slowly hydrolyzable analogue of GTP) exhibits the same type of broadening. Therefore, it is a characteristic feature of the GTP-bound form of Ras that the L1, L2, and L4 loop regions, but not other regions, are in a rather slow interconversion between two or more stable conformers. This phenomenon, termed a "regional polysterism", of these loop regions may be related with their multifunctionality: the GTP-dependent interactions with several downstream target groups such as the Raf and RalGDS families and also with the GTPase activating protein (GAP) family. In fact, the binding of Ras(1-171).GMPPNP with the Ras-binding domain (residues 51-131) of c-Raf-1 was shown to eliminate the regional polysterism nearly completely. It was indicated, therefore, that each target/regulator selects its appropriate conformer among those presented by the "polysteric" binding interface of Ras. As the downstream target groups exhibit no apparent sequence homology to each other, it is possible that one target group prefers a conformer different from that preferred by another group. The involvement of loop L1 in the regional polysterism might suggest that the negative regulators, GAPs, bind to the polysteric binding interface (loops L2 and L4) of Ras and cooperatively select a conformer suitable for transition of the GTPase catalytic center, involving loops L1 and L4, into the highly active state.

Conformational and dynamic differences between N-ras P21 bound to GTPgammaS and to GMPPNP as studied by NMR

Heteronuclear-edited proton-detected NMR methods are used to study the nucleotide-dependent conformational changes between the GMPPNP form of human N-ras P21 as compared to GDP and GTPgammaS forms. Full-length N-ras P21 was also compared with protein truncated beyond residue 167, to search for interaction points between the more invariant part of the protein and the variable C-terminal section. In both cases, the reporter was the 15N-H 2D spectrum of aspartate amide groups labeled with 15N. Small truncation-induced changes were seen in the spectrum at the resonances of Asp-54, -108, and -109 which are not far from the C-terminal and, surprisingly, at Asp-57 which is more remote. The spectrum obtained for the GMPPNP-ligated form is similar to that of the GTPgammaS form, except that peaks of several residues are weak at low temperature, and strongly temperature-dependent in their intensity, and a new resonance appears at 15 degrees C and above. The observations are discussed in terms of a two-state model for the GMPPNP-ligated protein, previously proposed by Geyer et al. [(1996) Biochemistry 35, 10308-10320].

Conformational transitions in p21ras and in its complexes with the effector protein Raf-RBD and the GTPase activating protein GAP

31P NMR revealed that the complex of p21ras with the GTP analog GppNHp.Mg2+ exists in two conformational states, states 1 and 2. In wild-type p21ras the equilibrium constant K1(12) between the two states is 1.09. The population of these states is different for various mutants but independent of temperature. The activation enthalpy delta H ++ and activation entropy delta S ++ for the conformational transitions were determined by full-exchange matrix analysis for wild-type p21ras and p21ras(S65P). For the wild-type protein one obtains delta H ++ = 89 +/- 2 kJ mol-1 and delta S ++ = 102 +/- 20 J mol-1 K-1 and for the mutant protein delta H ++ = 93 +/- 7 kJ mol-1 and delta S ++ = 138 +/- 30 J mol-1 K-1. The study of various p21ras mutants suggests that the two states correspond to different conformations of loop L2, with Tyr-32 in two different positions relative to the bound nucleotide. High-field EPR at 95 GHz suggest that the observed conformational transition does not directly influence the coordination sphere of the protein-bound metal ion. The influence of this transition on loop L4 was studied by 1H NMR with mutants E62H and E63H. There was no indication that L4 takes part in the transition described in L2, although a reversible conformational change could be induced by decreasing the pH value. The exchange between the two states is slow on the NMR time scale (< 10 s-1): at approximately pH 5 the population of the two states is equal. The interaction of p21ras-triphosphate complexes with the Ras-binding domain (RBD) of the effector protein c-Raf-1, Raf-RBD, and with the GTPase activating protein GAP was studied by 31P NMR spectroscopy. In complex with Raf-RBD the second conformation of p21ras (state 2) is stabilized. In this conformation Tyr-32 is located in close proximity to the phosphate groups of the nucleotide, and the beta-phosphate resonance is shifted upfield by 0.7 ppm. Spectra obtained in the presence of GAP suggest that in the ground state GAP does not interact directly with the nucleotide bound to p21ras and does not induce larger conformational changes in the neighborhood of the nucleotide. The experimental data are consistent with a picture where GAP accelerates the exchange process between the two states and simultaneously increases the population of state 1 at higher temperature.

Cell-free synthesis and amino acid-selective stable isotope labeling of proteins for NMR analysis

For the application of multidimensional NMR spectroscopy to larger proteins, it would be useful to perform selective labeling of one of the 20 amino acids. For some amino acids, however, amino acid metabolism drastically reduces the efficiency and selectivity of labeling in in vivo expression systems. In the present study, a cell-free protein synthesis system was optimized, so that highly efficient and selective stable isotope labeling of proteins can be achieved in the absence of amino acid metabolism. The productivity of the E. coli cell-free coupled transcription-translation system was first improved, by about fivefold, by using the T7 RNA polymerase for transcription and also by improving the translation conditions. Thus, about 0.1 mg protein per 1 ml reaction mixture was synthesized. Then, this improved cell-free system was used for Asp- or Ser-selective 15N-labeling of the human c-Ha-Ras protein. With a 15 ml cell-free reaction, using less than 1 mg of 15N-labeled amino acid, 1 mg of the Ras protein was obtained. 1H-15N HSQC experiments confirmed that the Ras protein was efficiently labeled with high selectivity. These results indicate that this cell-free protein synthesis system is useful for NMR studies.

Site-directed mutagenesis, fluorescence, and two-dimensional NMR studies on microenvironments of effector region aromatic residues of human c-Ha-Ras protein

The Tyr residues in positions 32 and 40 of human c-Ha-Ras protein were replaced by site-directed mutagenesis (Y32F, Y32W, Y40K, and Y40W) to examine their roles in the signal-transducing activity and the sensitivity to the GTPase activating protein (GAP). The signal-transducing activity of the oncogenic Ras protein in PC12 cells was lost upon mutations Y32F and Y40K, but retained upon mutations Y32W and Y40W. These results suggest that residues 32 and 40 are both required to have aromatic groups and residue 32 is further required to have a hydrogen donor. On the other hand, three mutations (Y32F, Y32W, and Y40W) caused no appreciable reduction in either GAP-binding affinity or GAP sensitivity. By the Y40K mutation, GAP-binding affinity was slightly lowered, while GAP sensitivity was drastically impaired. Therefore, for residues 32 and 40 of Ras, interactions with GAP appear to be different from those with the target of signal transduction in the PC12 cell. As for the Y32W-Ras protein bound with an unhydrolyzable GTP analogue (GMPPNP), the Trp32 fluorescence is appreciably red-shifted, weaker, and more susceptible to KI quenching as compared to that of the GDP-bound form. Two-dimensional NMR spectroscopy with selectively deuterated Ras proteins revealed fewer and weaker nuclear Overhauser effects on the aromatic protons of Trp32 in the GMPPNP-bound form than in the GDP-bound form. This indicates that the side chain of Trp32 is more exposed to the solvent in the GMPPNP-bound form than in the GDP-bound form.(ABSTRACT TRUNCATED AT 250 WORDS)

Mapping the nucleotide-dependent conformational change of human N-ras p21 in solution by heteronuclear-edited proton-observed NMR methods

Heteronuclear-edited proton-detected NMR methods are used to study the nucleotide-dependent conformational change between GDP- and GTP gamma S-bound forms of human N-ras p21. Amide groups of Asp are used as sensitive probes. When GTP gamma S is substituted for GDP in cellular N-ras p21, the chemical shifts of resonances Asp-47, -126, -154, and Asn-172, as well as Gly-77 and -151, are not sensitive to nucleotide exchange, whereas Asp-30, -33, -38, -54, -57, -69, -92, -105, and -119 are affected. Distinct chemical shift changes of Asp-33, -38, and -69 indicate that substantial structure changes occur in the effector-binding region and the switch II region. Crystallographic studies of H-ras p21 have indicated that the conformational differences are confined to residues 32-38 and 60-76. Our observations indicate that the nucleotide-dependent structural transitions of the protein in solution may not be identical to those in the crystal. They suggest that the peptide beyond Glu-76 participates in a conformational switch, and possibly is involved in effector function. We propose that the region roughly from Asp-92 to -105, and the region of guanine base-binding motif(s), e.g., 116NKXD, are candidate sites recognized by either a GDP/GTP release factor or a GTPase-affected protein.