M. Rachel
Mehaffey
a,
Christopher L.
Schardon
b,
Elisa T.
Novelli
a,
Michael B.
Cammarata
a,
Lauren J.
Webb
a,
Walter
Fast
b and
Jennifer S.
Brodbelt
*a
aDepartment of Chemistry, University of Texas at Austin, Austin, TX 78712-0165, USA. E-mail: jbrodbelt@cm.utexas.edu; Tel: +1-512-471-0028
bDivision of Chemical Biology and Medicinal Chemistry, College of Pharmacy, University of Texas at Austin, Austin, TX 78712, USA
First published on 15th July 2019
Mutations in the GTPase enzyme K-Ras, specifically at codon G12, remain the most common genetic alterations in human cancers. The mechanisms governing activation of downstream signaling pathways and how they relate back to the identity of the mutation have yet to be completely defined. Here we use native mass spectrometry (MS) combined with ultraviolet photodissociation (UVPD) to investigate the impact of three G12X mutations (G12C, G12V, G12S) on the homodimerization of K-Ras as well as heterodimerization with a downstream effector protein, Raf. Electrospray ionization (ESI) was used to transfer complexes of WT or G12X K-Ras bound to guanosine 5′-diphosphate (GDP) or GppNHp (non-hydrolyzable analogue of GTP) into the gas phase. Relative abundances of homo- or hetero-dimer complexes were estimated from ESI-MS spectra. K-Ras + Raf heterocomplexes were activated with UVPD to probe structural changes responsible for observed differences in the amount of heterocomplex formed for each variant. Holo (ligand-bound) fragment ions resulting from photodissociation suggest the G12X mutants bind Raf along the expected effector binding region (β-interface) but may interact with Raf via an alternative α-interface as well. Variations in backbone cleavage efficiencies during UV photoactivation of each variant were used to relate mutation identity to structural changes that might impact downstream signaling. Specifically, oncogenic upregulation for hydrogen-bonding amino acid substitutions (G12C, G12S) is achieved by stabilizing β-interface interactions with Raf, while a bulkier, hydrophobic G12V substitution leads to destabilization of this interface and instead increases the proximity of residues along the α-helical bundles. This study deciphers new pieces of the complex puzzle of how different K-Ras mutations exert influence in downstream signaling.
In recent years, mass spectrometry (MS) has proven to be a useful experimental technique to probe theoretical models and address structural biology questions in a sensitive, rapid manner.37,38 The use of covalent chemical probes, such as hydrogen/deuterium exchange (HDX) or cross-linking reagents, in conjunction with bottom-up tandem mass spectrometry (MS/MS) is now a standard approach for evaluating solvent accessibilities and mapping protein interfaces.39–41 Top-down techniques have more recently gained traction for characterization of protein structure with the advent of native MS and advances in MS/MS methods.37,38 The framework for native MS entails the transfer of intact protein complexes to the gas phase by electrospray ionization of solutions of high ionic strength which largely preserves tertiary and quaternary structures of macromolecules.42 Consequently, native MS is now routinely used for identification of oligomeric states of proteins and elucidation of binding stoichiometry with other proteins or ligands.38 To further discern three-dimensional shapes of protein complexes in the gas phase, ion mobility (IM) mass spectrometry has proved to be an indispensable technique for providing information complementary to static high resolution structures resulting from X-ray crystallography, NMR spectroscopy, or cryo-EM.43
MS/MS methods sensitive to protein structure have been developed to study the architecture of protein–protein and protein–ligand complexes. Traditional collisional activation of native-like complexes mainly yields ejection of ligands or subunits,44,45 which is often uninformative and has spurred the application of alternative activation methods. Electron-based methods, including electron transfer dissociation (ETD) and electron capture dissociation (ECD), and surface induced dissociation (SID) represent the most widely used alternative methods to characterize protein structure.46–52 ETD and ECD yield mostly diagnostic sequence ions with abundances that can be correlated with crystallographic B-factors and used to probe the topology of protein assemblies for flexible regions, thus revealing insight into the higher order structure of the complexes.46–51 Conversely, SID causes disassembly of protein–protein complexes into intact subunits in a way that conserves the symmetry of charge distribution among the separating proteins or subunits.52 Coupled with IM, this MS/MS technique allows for the rapid determination of stoichiometry and topology for computationally designed protein complexes.53 Ultraviolet photodissociation (UVPD) represents a third alternative MS/MS method that utilizes fast, high-energy excitation via 193 nm photons and yields the widest array of diagnostic sequence ions for protein complexes.54–58 The abundances of holo (ligand-bound) and apo (ligand-free) product ions resulting from activation of protein-ligand complexes by 193 nm UVPD have been demonstrated to reflect secondary and tertiary structural features.59–62 The suppression or enhancement of UVPD at each position along the backbone is thought to be modulated by non-covalent interactions that stabilize structural features and prevent separation of fragments after bond cleavage. Recent studies have demonstrated this approach is sensitive to loop movements upon ligand binding in dihydrofolate reductase (DHFR),60 dynamic plasticity throughout a catalytic cycle of the active enzyme adenylate kinase,61 and even long-range conformational changes resulting from single residue mutations in K-Ras.62
Previously variations in UVPD backbone cleavage efficiencies were monitored for monomeric protein-ligand complexes of K-Ras (wild-type (WT) and three G12X mutants: G12C, G12V, G12S) upon exchange of GDP for a non-hydrolyzable GTP analogue, guanosine 5′-[β,γ-imido]triphosphate (GppNHp).62 Switching from the inactive diphosphate bound state to the active triphosphate bound protein yielded unique structural changes for each mutant. Based on the examination of monomeric K-Ras and a series of G12X variants, the way each mutation modulated homo- or hetero-dimer formation with downstream effectors was inferred. Given that the G12 position lies on the surface of the protein, longer-range conformational changes in areas of the protein related to dimer formation are likely involved. Specifically, observations included stabilization of the α-helical bundle for G12C K-Ras, stabilization of the β-sheet region for G12V K-Ras, and increased flexibility of the β-sheet region for G12S K-Ras.62 Although the mechanism is not yet well defined, this supports that the identity of the G12 substitution impacts downstream signaling and patient survival differently.63–69
Here we report the use of native MS and UVPD to directly interrogate homo- and hetero-dimers of WT and three clinically relevant G12X mutants of K-Ras (G12C, G12S, G12V). We use native MS to define the relative abundance of homodimer formation for each K-Ras G12X mutant bound to either GDP or GppNHp (the non-hydrolyzable analogue of GTP), and compare the amount of heterodimer formed between each of the K-Ras mutants and the Ras binding domain (RBD) of a downstream effector kinase, rapidly accelerated fibrosarcoma (Raf). Comparing UVPD cleavage efficiencies of K-Ras·GppNHp + RafRBD heterodimers for each of the mutants to the WT suggests hydrogen-bonding amino acid substitutions (G12C, G12S) rely on directly stabilizing the β-interface interactions with Raf for oncogenic upregulation, whereas a bulkier, hydrophobic G12V substitution causes destabilization of this interface and instead results in tightened α-helical bundles. Direct structural interrogation of intact dimers by UVPD-MS represents an advance in unraveling the mutation-dependent interplay of structure and binding interactions for K-Ras signaling phenotypes.
Owing to the fact that UVPD-MS is a new approach for characterization of structural variations in protein complexes and providing insight into potential changes in conformation and binding motifs, there is need for critical evaluation of the method via appropriate controls and statistical tests. Some of the statistical tests and controls are described above and shown in Fig. S2, S9 and S16.† Moreover, the K-Ras protein used in the present study was also the subject of a prior study (ref. 62). Additional controls and statistical tests were performed previously, as provided in Fig. S2–S4† in ref. 62, including the assessment of the statistical significance of backbone cleavage efficiency for K-Ras and the G12C, G12V, and G12S mutants, both as native-like and denatured forms. One relevant control entailed the examination of UVPD of the denatured proteins, confirming that the relative fragment ion abundances were identical for 166 out of 169 comparisons of specific backbone positions. In contrast, the UVPD fragment ion abundances of the native-like proteins showed numerous statistically significant differences in fragment ion abundances. These types of control experiments support the premise that UVPD is sensitive to protein structure.
The distribution of homodimers versus monomers for each variant are estimated by comparing normalized abundances of the ions in the native ESI mass spectra shown in Fig. S4.†Fig. 1 illustrates the relative abundances of each observed K-Ras homodimer normalized to the total ion current for WT, G12C, G12V, and G12S K-Ras bound to GDP or GppNHp for protein concentrations of 20 or 80 μM. A larger bar indicates a greater portion of the dimeric form of the protein relative to the monomeric protein. GppNHp-bound WT K-Ras was observed to form significantly more abundant homodimers compared to its GDP-bound counterpart. This result is in line with the significant role of dimerization in the activation of K-Ras. The same general trend indicating enhanced formation of dimers containing GppNHp over GDP was observed for both G12C and G12S K-Ras, and for each of these two variants the overall abundances of homodimers containing GDP or GppNHp were substantially higher than the abundances of WT homodimers. Since the G12C and G12S mutations result in constitutively active variants of K-Ras, higher abundances of the homodimers support the idea that oligomerization is important for activation of the enzyme. Interestingly, the abundances of G12V K-Ras homodimers were low in general and were only observed at the higher (80 μM) protein concentration, suggesting dimerization is not a significant factor in maintaining an active state when the bulky Val is substituted for Gly.
Fig. 1 Graph displaying the abundance of each K-Ras homodimer normalized to the total ion current (TIC) for the native ESI mass spectra (Fig. S4†) for WT (blue), G12C (red), G12V (green), and G12S (yellow) bound to two GDP or two GppNHp (1:1 ratio) at a protein concentration of 20 μM (solid bars) or 80 μM (hashed bars). GppNHp is a non-hydrolyzable analogue of GTP and represents the active state ligand of K-Ras. The theoretical and measured monoisotopic masses for each observed homodimer are reported in Table S1.† |
In general, the relatively low percentage of dimers (<1% to 11% in Fig. 1) compared to the corresponding monomers is not surprising given that the protein is in solution with no membrane localization to aid in self-assembly.23,28 Previous dynamic light scattering experiments observed higher order oligomers only for K-Ras bound to another GTP-analogue, guanosine 5′-O-[γ-thio]triphosphate (GTP-γ-S), not GDP nor GppNHP.26 Although GTP-γ-S is expected to be a better mimic of GTP, for the experimental conditions tested herein this compound hydrolyzed too quickly to allow confident measurements of protein complexation (Fig. S5†).
There is accumulating evidence that WT K-Ras can act as a tumor suppressor and counteract the activating effects of G12X mutations.72 To mimic cells heterozygous for G12X mutations, each of the three mutants bound to GDP or GppNHp was mixed in a 1:1 ratio with WT K-Ras at 20 or 80 μM total protein concentration prior to ESI. No dimers were observed in the ESI-MS spectra, thus suggesting that homodimers composed of two protein variants (WT and G12X K-Ras) are less stable or membrane localization is required for formation (Fig. S6†).
Fig. 2 Native ESI mass spectra of 1:1 K-Ras:RafRBD solution of WT K-Ras bound to (A) GDP or (B) GppNHp (1:1 ratio). Observed species are labelled as colored circles (6+, 7+, 8+, 9+ RafRBD monomers, 9+, 10+, 11+ K-Ras monomers, and 10+, 11+, 12+, 13+, 14+ heterodimers). Corresponding spectra for G12C, G12V, and G12S K-Ras are shown in Fig. S7.† K-Ras:RafRBD heterodimers were only observed for solutions containing GppNHp (non-hydrolyzable analogue of GTP), not GDP. The relative abundances of each species are estimated from SEC-MS and reported in Fig. 3. The 12+ heterodimer was selectively isolated and activated to yield the UVPD spectra shown in Fig. S10.† |
Given that the observed monomeric and dimeric protein complexes span a broad m/z range, on-line size-exclusion chromatography was used to separate the various complexes in order to allow optimization of MS tuning parameters and improve the accurate measurement of abundances. Fig. S8† illustrates extracted ion chromatograms (EICs) derived for the m/z values corresponding to each key species (K-Ras·GppNHp, RafRBD, and K-Ras·GppNHp + RafRBD). Peak areas from the EICs were used to estimate the distributions of the monomeric proteins (K-Ras·GppNHp and RafRBD) and heterodimers (K-Ras·GppNHp + RafRBD) for WT K-Ras and each of the three variants (Fig. 3). The portions of K-Ras·GppNHp + RafRBD heterodimers (green segments) are the focal point of Fig. 3. G12V yields a significantly lower percentage of heterodimers than WT K-Ras, and conversely G12C and G12S display a larger percentage of heterodimers than WT K-Ras. This general trend in the formation of K-Ras·GppNHp + RafRBD heterodimers mirrors the results reported above for K-Ras homodimer formation. These results suggest that amino acid substitutions (G12C, G12S) that have enhanced hydrogen-bonding capabilities contribute to stabilizing the β-interface interactions with Raf. In contrast, the bulkier, hydrophobic G12V substitution causes destabilization of this interface.
Fig. 3 Distribution of the normalized relative abundances of K-Ras·GppNHp (red), RafRBD (blue), and K-Ras·GppNHp + RafRBD heterodimer (green) from on-line SEC of solutions containing WT, G12C, G12V, or G12S K-Ras·GppNHp with an equimolar amount of RafRBD. Peak areas of the EIC traces of m/z values corresponding to each species were used to calculate the reported percentages (Fig. S8†). |
The fragment ions generated upon UVPD can be further categorized as apo (RafRBD-free) or holo (containing RafRBD) product ions. While traditional HCD typically causes ejection of the intact ligand prior to backbone fragmentation owing to the preferential cleavage of non-covalent interactions during collisional activation, UV photoactivation of protein-ligand complexes yields diagnostic sequence ions that retain non-covalent interactions with the bound ligands, termed holo fragment ions.55,60–62,73 HCD has been reported to yield holo ions consisting of a portion of a non-covalently bound nucleotide74 but UVPD consistently allows retention of intact nucleotides like GTP or GppNHp.62 The high energy imparted to a protein via absorption of 193 nm photons results in activation to excited electronic states and affords preferential cleavage of backbone bonds rather than disruption of electrostatic interactions with bound ligands. Consequently, mapping the location of observed holo ions along the protein backbone affords insight into ligand-binding interfaces. Fig. 4A displays the sequence of K-Ras with indicators above specific residues to show the sites of backbone cleavages that lead to production of holo ions. The indicators are color-coded to reflect whether the backbone cleavage sites correspond to formation of N-terminal ions (a,b,c), C-terminal ions (x,y,z), or bi-directional ions (i.e. backbone cleavages resulting in complementary N- and C-terminal holo ions). Although UVPD of the K-Ras·RafRBD complexes caused backbone fragmentation of both the K-Ras and RafRBD subunits, for streamlined data interpretation RafRBD was treated as one large ligand and only holo ions containing the intact RafRBD protein were considered. The crystal structure of K-Ras bound to RafRBD is shown in Fig. 4B with residues corresponding to the observed sites of backbone cleavages color coded to aid in visualization for WT K-Ras (and corresponding structures for the three G12X mutants given in Fig. S13†).
Fig. 4 (A) The backbone cleavage sites upon UVPD of the K-Ras·GppNHp + RafRBD heterodimer (12+) are mapped above the sequence of K-Ras (X13 = G, C, V, or S) as lines or symbols for WT (solid), G12C (dashes), G12V (dots), and G12S (pattern). The backbone cleavages that result in N-terminal ions (a,b,c) are colored blue; those that result in C-terminal ions (x,y,z) are colored red; and those that yield complementary N-terminal and C-terminal ions are categorized as bi-directional and are colored green. The β-interface region is shaded in purple. (B) Crystal structure of K-Ras bound to GppNHp and complexed with RafRBD (PDB ID: 4G0N) with residues colored corresponding to N-terminal (a,b,c-type ions; blue), C-terminal (x,y,z-type ions; red), and bi-directional (green) holo fragment ions observed during UVPD of the 12+ WT heterodimer. Corresponding structures for the G12X mutants are given in Fig. S13.† The original UVPD spectra are shown in Fig. S10.† The GppNHp ligand (a non-hydrolyzable analogue of GTP) is shown as pink sticks and labelled in (B). |
Of the two known interfaces along K-Ras (α- and β-interfaces), effector binding is expected to occur at the stronger β-interface containing Switch-I (residues 27–36, 41–45, shaded in Fig. 4A).26 The observed backbone cleavage sites that lead to RafRBD-bound holo ions upon UVPD of the K-Ras·RafRBD complexes support the β-interface as the primary binding region for WT K-Ras as well as for the three G12X mutants. In particular, backbone cleavages that result in C-terminal (x,y,z) holo ions only occur in the region of K-Ras that is N-terminal to the β-interface, and backbone cleavages that lead to N-terminal (a,b,c) ions occur exclusively in the region C-terminal to the β-interface (Fig. 4A). Moreover, the bi-directional holo ions that result in complementary C-terminal and N-terminal product ions originate from backbone cleavages along or near this expected interface. Further evidence is provided by the apo sequence coverage maps (Fig. S12†) in which coverage of the β-interfacial region (shaded in purple) is relatively sparse, suggesting two possibilities. One possibility is that backbone cleavages in the β-interfacial region might instead preferentially produce the mass-shifted holo fragment ions showcased in Fig. 4A (thus depleting the abundances of apo sequence ions). Alternatively, fragmentation along the β-interfacial region might be suppressed owing to its engagement in interactions with RafRBD (thus stabilizing the β-interfacial region and preventing separation of fragment ions).
Interestingly, there are several C-terminal and bi-directional holo fragment ions originating from backbone cleavages adjacent to residues N86-E92 (C-terminal to the β-interface) which are observed for all three G12X mutants but not WT K-Ras (Fig. 4A). This region is part of the weaker α-interface (residues 86–105, 126–138) where K-Ras homodimerization is thought to occur.16,26 Perhaps G12X mutants also use the α-interface to bind effector proteins and maintain an active state. However, the fact that fewer holo ions originate from this region suggests either it is a lower population effector binding region compared to the β-interface or the non-covalent interactions stabilizing this alternative binding mode are too weak to survive photodissociation (thus preventing survival of holo ions). The observations about the involvement of the α-interface merit future investigation using other traditional biophysical methods.
Fig. 5 Heat maps of the enhancement (red) and suppression (blue) of UVPD for RafRBD-bound heterodimers shown as a linear map across the sequence (A) or displayed along the crystal structure of K-Ras bound to GppNHp (non-hydrolyzable analogue of GTP) and complexed with RafRBD for (B) G12C, (C) G12V, and (D) G12S K-Ras variants relative to WT K-Ras. These heat maps are difference plots: (i.e. G12X – WT). The colored regions represent statistically significant changes in cleavage efficiency of the protein backbone during UVPD based on the difference plots in Fig. S14.† The GppNHp ligand (pink) is shown as sticks and labelled in (B). The β-interface region is highlighted in purple in (B). Relevant loops, α-helices, and interfaces are labelled underneath the x-axis of (A). |
To gain insight into how the identity of the G12X mutations impacts the effects on signaling, we focus on conformational differences observed in three key regions: β-interface (switch I: residues 27–36, 41–45), switch II (residues 58–64), and α-interface (α-C and α-D helices: residues 86–105 and 126–138). The three key regions of the difference plots shown in Fig. S15† are expanded in Fig. S16,† and the statistically meaningful variations at the 99% confidence level (p ≤ 0.01) are demarcated to help highlight the trends in UVPD backbone cleavage efficiencies for each of the three variants compared to WT K-Ras. Fig. S17† summarizes the observed enhancement or suppression in UVPD backbone cleavage efficiency of those three regions for RafRBD-bound heterodimers of G12C, G12V, and G12S compared to WT K-Ras. G12C and G12S result in similar changes in these regions, mainly stabilization of the β-interface (e.g., suppression of UVPD). Conversely the G12V substitution favors increased proximity of the helices in the α-interfacial region (observed as suppression of UVPD) and destabilization of the β-interface (resulting in enhancement of UVPD).
These observations bring new light to a previous study that used UVPD-MS to evaluate conformational changes that occurred for K-Ras and its three G12X mutants during the GTP loading step of K-Ras activation.62 In the previous study, a more stable α-helical region was inferred for G12C K-Ras in comparison to WT K-Ras based on suppression of UVPD in that region, whereas G12S K-Ras exhibited more flexibility in the β-interface region based on enhancement of UVPD.62 Both of these findings are consistent with the adoption of a more stable β-interface upon RafRBD binding in the present study (i.e. given the β-interface was found to be more flexible, more intramolecular contacts could be formed upon Raf binding). For G12V K-Ras compared to WT K-Ras, the converse was observed previously: the β-interface appeared to be more rigid (stabilized by more interactions; lower UVPD fragmentation efficiency) after substitution of GppNHp for GDP.62 Again this finding correlates with the more stabilized α-interface of G12V after RafRBD binding as evidenced by the decreased UVPD fragmentation efficiency in the present study. This outcome implies engagement of fewer interactions between RafRBD and the β-interface, since G12V K-Ras was more rigidly pre-organized prior to effector binding. Collectively these results suggest that mutations of K-Ras which introduce hydrogen-bonding groups (G12C, G12S) result in a stabilized β-interface, although it is also acknowledged that thiols form weaker hydrogen bonds than alcohols.75 In contrast, a bulkier, more hydrophobic substitution (G12V) relies preferentially on contacts along the α-interface for oncogenic activation of K-Ras.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sc01032g |
This journal is © The Royal Society of Chemistry 2019 |