Oncogenic and RASopathy-associated K-RAS mutations relieve membrane-dependent occlusion of the effector-binding site

Pioneer in Molecular Biology: Conformational Ensembles in Molecular Recognition, Allostery, and Cell Function

In 1978, for my PhD, I developed the efficient O(n3) dynamic programming algorithm for the-then open problem of RNA secondary structure prediction. This algorithm, now dubbed the "Nussinov algorithm", "Nussinov plots", and "Nussinov diagrams", is still taught across Europe and the U.S. As sequences started coming out in the 1980s, I started seeking genome-encoded functional signals, later becoming a bioinformatics trend. In the early 1990s I transited to proteins, co-developing a powerful computer vision-based docking algorithm. In the late 1990s, I proposed the foundational role of conformational ensembles in molecular recognition and allostery. At the time, conformational ensembles and free energy landscapes were viewed as physical properties of proteins but were not associated with function. The classical view of molecular recognition and binding was based on only two conformations captured by crystallography: open and closed. I proposed that all conformational states preexist. Proteins always have not one folded form-nor two-but many folded forms. Thus, rather than inducing fit, binding can work by shifting the ensembles between states, and this shifting, or redistributing the ensembles to maintain equilibrium, is the origin of the allosteric effect and protein, thus cell, function. This transformative paradigm impacted community views in allosteric drug design, catalysis, and regulation. Dynamic conformational ensemble shifts are now acknowledged as the origin of recognition, allostery, and signaling, underscoring that conformational ensembles-not proteins-are the workhorses of the cell, pioneering the fundamental idea that dynamic ensembles are the driving force behind cellular processes. Nussinov was recognized as pioneer in molecular biology by JMB.

The Effect of Phosphoserine-Containing Membranes on Electrostatic Fields at the Protein-Protein Interface Measured through Vibrational Stark Effect Spectroscopy

In the cell, Ras GTPases function as membrane-bound molecular switches for a variety of cell signaling pathways. Ras isoforms have long been of interest because of the connection between amino acid mutations and tumorigenesis. Much research focused on Ras has used truncated, solubilized constructs, which exclude the membrane-binding domain and therefore ignore the effects of membrane binding on Ras function. Since the membrane is a highly charged surface, it could have a significant impact on the electrostatic environment at or near the protein-protein interface. Here, we use a thiocyanate probe chemically inserted into the Ras-binding domain of RalGDS to investigate the effect of membrane binding at the Ras active site. Changes in the electric field caused by the membrane were measured by the probe as vibrational energy shifts in the infrared (IR) spectrum. For a selection of mutants which caused large shifts at this interface on the soluble H-Ras construct, binding to a 30% phosphatidylserine (PS)/70% phosphatidylcholine (PC) nanodisc caused reduced shifts compared to the solubilized counterparts. Additionally, the vibrational probe bonded to the wildtype (WT) Ras construct demonstrated a shift of 0.7 cm-1 as a PC nanodisc was doped from 0% to 30% PS, but mutations introduced to the Ras active site caused the probe to show no shift across these PS concentrations. These results indicate that the local membrane environment has an effect on the electrostatics at the Ras active site and needs to be considered when investigating the effect of oncogenic mutations on Ras function.

Blocking C-terminal processing of KRAS4b via a direct covalent attack on the CaaX-box cysteine

RAS is the most frequently mutated oncogene in cancer. RAS proteins show high sequence similarities in their G-domains but are significantly different in their C-terminal hypervariable regions (HVR). These regions interact with the cell membrane via lipid anchors that result from posttranslational modifications (PTM) of cysteine residues. KRAS4b is unique as it has only one cysteine that undergoes PTM, C185. Small molecule covalent modification of C185 would block any form of prenylation and subsequently inhibit attachment of KRAS4b to the cell membrane, blocking its biological activity. We translated this concept to the discovery and development of disulfide tethering screen hits into irreversible covalent modifiers of C185. These compounds inhibited proliferation of KRAS4b-driven mouse embryonic fibroblasts, but not cells driven by N-myristoylated KRAS4b that harbor a C185S mutation and are not dependent on C185 prenylation. Top-down proteomics was used to confirm target engagement in cells. These compounds bind in a pocket formed when the HVR folds back between helix 3 and 4 in the G-domain (HVR-α3-α4). This interaction can happen in the absence of small molecules as predicted by molecular dynamics simulations and is stabilized in the presence of C185 binders as confirmed by small-angle X-ray scattering and solution NMR. NOESY-HSQC, an NMR approach that measures internuclear distances of 6 Å or less, and structure analysis identified the critical residues and interactions that define the HVR-α3-α4 pocket. Further development of compounds that bind to this pocket could be the basis of a new approach to targeting KRAS cancers.

The small GTPase MRAS is a broken switch

Intense research on founding members of the RAS superfamily has defined our understanding of these critical signalling proteins, leading to the premise that small GTPases function as molecular switches dependent on differential nucleotide loading. The closest homologs of H/K/NRAS are the three-member RRAS family, and interest in the MRAS GTPase as a regulator of MAPK activity has recently intensified. We show here that MRAS does not function as a classical switch and is unable to exchange GDP-to-GTP in solution or when tethered to a lipid bilayer. The exchange defect is unaffected by inclusion of the GEF SOS1 and is conserved in a distal ortholog from nematodes. Synthetic activating mutations widely used to study the function of MRAS in a presumed GTP-loaded state do not increase exchange, but instead drive effector binding due to sampling of an activated conformation in the GDP-loaded state. This includes nucleation of the SHOC2-PP1Cα holophosphatase complex. Acquisition of NMR spectra from isotopically labeled MRAS in live cells validated the GTPase remains fully GDP-loaded, even a supposed activated mutant. These data show that RAS GTPases, including those most similar to KRAS, have disparate biochemical activities and challenge current dogma on MRAS, suggesting previous data may need reinterpretation.

Motif-guided identification of KRAS-interacting proteins

BACKGROUND

For decades, KRAS has always been a huge challenge to the field of drug discovery for its significance in cancer progression as well as its difficulties in being targeted as an "undruggable" protein. KRAS regulates downstream signaling pathways through protein-protein interactions, whereas many interaction partners of KRAS remain unknown.

RESULTS

We developed a workflow to computationally predict and experimentally validate the potential KRAS-interacting proteins based on the interaction mode of KRAS and its known binding partners. We extracted 17 KRAS-interacting motifs from all experimentally determined KRAS-containing protein complexes as queries to identify proteins containing fragments structurally similar to the queries in the human protein structure database using our in-house protein-protein interaction prediction method, PPI-Miner. Finally, out of the 78 predicted potential interacting proteins of KRAS, 10 were selected for experimental validation, including BRAF, a previously reported interacting protein, which served as the positive control in our validation experiments. Additionally, a known peptide that binds to KRAS, KRpep-2d, was also used as a positive control. The predicted interacting motifs of these 10 proteins were synthesized to perform biolayer interferometry assays, with 4 out of 10 exhibiting binding affinities to KRAS, and the strongest, GRB10, was selected for further validation. Additionally, the interaction between GRB10 (RA-PH domain) and KRAS was confirmed via immunofluorescence and co-immunoprecipitation.

CONCLUSIONS

These results demonstrate the effectiveness of our workflow in predicting potential interacting proteins for KRAS and deepen the understanding of KRAS-driven tumor mechanisms and the development of therapeutic strategies.

From disorder comes function: Regulation of small GTPase function by intrinsically disordered lipidated membrane anchor

The intrinsically disordered, lipid-modified membrane anchor of small GTPases is emerging as a critical modulator of function through its ability to sort lipids in a conformation-dependent manner. We reviewed recent computational and experimental studies that have begun to shed light on the sequence-ensemble-function relationship in this unique class of lipidated intrinsically disordered regions (LIDRs).

KRAS Mutation Subtypes and Their Association with Other Driver Mutations in Oncogenic Pathways

The KRAS mutation stands out as one of the most influential oncogenic mutations, which directly regulates the hallmark features of cancer and interacts with other cancer-causing driver mutations. However, there remains a lack of precise information on their cooccurrence with mutated variants of KRAS and any correlations between KRAS and other driver mutations. To enquire about this issue, we delved into cBioPortal, TCGA, UALCAN, and Uniport studies. We aimed to unravel the complexity of KRAS and its relationships with other driver mutations. We noticed that G12D and G12V are the prevalent mutated variants of KRAS and coexist with the TP53 mutation in PAAD and CRAD, while G12C and G12V coexist with LUAD. We also noticed similar observations in the case of PIK3CA and APC mutations in CRAD. At the transcript level, a positive correlation exists between KRAS and PIK3CA and between APC and KRAS in CRAD. The existence of the co-mutation of KRAS and other driver mutations could influence the signaling pathway in the neoplastic transformation. Moreover, it has immense prognostic and predictive implications, which could help in better therapeutic management to treat cancer.

Reconstitution and characterization of BRAF in complex with 14-3-3 and KRAS4B on nanodiscs

RAF kinases are key components of the RAS-MAPK signaling pathway, which drives cell growth and is frequently overactivated in cancer. Upstream signaling activates the small GTPase RAS, which recruits RAF to the cell membrane, driving a transition of the latter from an auto-inhibited monomeric conformation to an active dimer. Despite recent progress, mechanistic details underlying RAF activation remain unclear, particularly the role of RAS and the membrane in mediating this conformational rearrangement of RAF together with 14-3-3 to permit RAF kinase domain dimerization. Here, we reconstituted an active complex of dimeric BRAF, a 14-3-3 dimer and two KRAS4B on a nanodisc bilayer and verified that its assembly is GTP-dependent. Biolayer interferometry (BLI) was used to compare the binding affinities of monomeric versus dimeric full-length BRAF:14-3-3 complexes for KRAS4B-conjugated nanodiscs (RAS-ND) and to investigate the effects of membrane lipid composition and spatial density of KRAS4B on binding. 1,2-Dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) and higher KRAS4B density enhanced the interaction of BRAF:14-3-3 with RAS-ND to different degrees depending on BRAF oligomeric state. We utilized our reconstituted system to dissect the effects of KRAS4B and the membrane on the kinase activity of monomeric and dimeric BRAF:14-3-3 complexes, finding that KRAS4B or nanodiscs alone were insufficient to stimulate activity, whereas RAS-ND increased activity of both states of BRAF. The reconstituted assembly of full-length BRAF with 14-3-3 and KRAS on a cell-free, defined lipid bilayer offers a more holistic biophysical perspective to probe regulation of this multimeric signaling complex at the membrane surface.

RAS G-domains allosterically contribute to the recognition of lipid headgroups and acyl chains

Mutant RAS are major contributors to cancer and signal primarily from nanoclusters on the plasma membrane (PM). Their C-terminal membrane anchors are main features of membrane association. However, the same RAS isoform bound to different guanine nucleotides spatially segregate. Different RAS nanoclusters all enrich a phospholipid, phosphatidylserine (PS). These findings suggest more complex membrane interactions. Our electron microscopy-spatial analysis shows that wild-types, G12V mutants, and membrane anchors of isoforms HRAS, KRAS4A, and KRAS4B prefer distinct PS species. Mechanistically, reorientation of KRAS4B G-domain exposes distinct residues, such as Arg 135 in orientation state 1 (OS1) and Arg 73/Arg 102 in OS2, to the PM and differentially facilitates the recognition of PS acyl chains. Allele-specific oncogenic mutations of KRAS4B also shift G-domain reorientation equilibrium. Indeed, KRAS4BG12V, KRAS4BG12D, KRAS4BG12C, KRAS4BG13D, and KRAS4BQ61H associate with PM lipids with headgroup and acyl chain specificities. Distribution of these KRAS4B oncogenic mutants favors different nanoscale membrane topography. Thus, RAS G-domains allosterically facilitate membrane lateral distribution.

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.

Membrane-Driven Dimerization of the Peripheral Membrane Protein KRAS: Implications for Downstream Signaling

Transient homo-dimerization of the RAS GTPase at the plasma membrane has been shown to promote the mitogen-activated protein kinase (MAPK) signaling pathway essential for cell proliferation and oncogenesis. To date, numerous crystallographic studies have focused on the well-defined GTPase domains of RAS isoforms, which lack the disordered C-terminal membrane anchor, thus providing limited structural insight into membrane-bound RAS molecules. Recently, lipid-bilayer nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses have revealed several distinct structures of the membrane-anchored homodimers of KRAS, an isoform that is most frequently mutated in human cancers. The KRAS dimerization interface is highly plastic and altered by biologically relevant conditions, including oncogenic mutations, the nucleotide states of the protein, and the lipid composition. Notably, PRE-derived structures of KRAS homodimers on the membrane substantially differ in terms of the relative orientation of the protomers at an "α-α" dimer interface comprising two α4-α5 regions. This interface plasticity along with the altered orientations of KRAS on the membrane impact the accessibility of KRAS to downstream effectors and regulatory proteins. Further, nanodisc platforms used to drive KRAS dimerization can be used to screen potential anticancer drugs that target membrane-bound RAS dimers and probe their structural mechanism of action.

Structural insights into the complex of oncogenic KRas4BG12V and Rgl2, a RalA/B activator

About a quarter of total human cancers carry mutations in Ras isoforms. Accumulating evidence suggests that small GTPases, RalA, and RalB, and their activators, Ral guanine nucleotide exchange factors (RalGEFs), play an essential role in oncogenic Ras-induced signalling. We studied the interaction between human KRas4B and the Ras association (RA) domain of Rgl2 (Rgl2RA), one of the RA-containing RalGEFs. We show that the G12V oncogenic KRas4B mutation changes the interaction kinetics with Rgl2RA The crystal structure of the KRas4BG12V: Rgl2RA complex shows a 2:2 heterotetramer where the switch I and switch II regions of each KRasG12V interact with both Rgl2RA molecules. This structural arrangement is highly similar to the HRasE31K:RALGDSRA crystal structure and is distinct from the well-characterised Ras:Raf complex. Interestingly, the G12V mutation was found at the dimer interface of KRas4BG12V with its partner. Our study reveals a potentially distinct mode of Ras:effector complex formation by RalGEFs and offers a possible mechanistic explanation for how the oncogenic KRas4BG12V hyperactivates the RalA/B pathway.

Characterization of nanodisc-forming peptides for membrane protein studies

Lipid-bilayer nanodiscs provide a stable, native-like membrane environment for the functional and structural studies of membrane proteins and other membrane-binding molecules. Peptide-based nanodiscs having unique properties are developed for membrane protein studies and other biological applications. While the self-assembly process rendering the formation of peptide-nanodiscs is attractive, it is important to understand the stability and suitability of these nanodisc systems for membrane protein studies. In this study, we investigated the nanodiscs formation by the anti-inflammatory and tumor-suppressing peptide AEM28. AEM28 is a chimeric peptide containing a cationic-rich heparan sulfate proteoglycan- (HSPG)-binding domain from human apolipoprotein E (hapoE) (141-150) followed by the 18A peptide's amino acid sequence. AEM28-based nanodiscs made with different types of lipids were characterized using various biophysical techniques and compared with the nanodiscs formed using 2F or 4F peptides. Variable temperature dynamic light-scattering and 31P NMR experiments indicated the fusion and size heterogeneity of nanodiscs at high temperatures. The suitability of AEM28 and Ac-18A-NH2- (2F-) based nanodiscs for studying membrane proteins is demonstrated by reconstituting and characterizing a drug-metabolizing enzyme, cytochrome-P450 (CYP450), or the redox complex CYP450-CYP450 reductase. AEM28 and 2F were also tested for their efficacies in solubilizing E. coli membranes to understand the possibility of using them for detergent-free membrane protein isolation. Our experimental results suggest that AEM28 nanodiscs are suitable for studying membrane proteins with a net positive charge, whereas 2F-based nanodiscs are compatible with any membrane proteins and their complexes irrespective of their charge. Furthermore, both peptides solubilized E. coli cell membranes, indicating their use in membrane protein isolation and other applications related to membrane solubilization.

In silico design of a lipid-like compound targeting KRAS4B-G12D through non-covalent bonds

One of the most common drivers in human cancer is the peripheral membrane protein KRAS4B, able to promote oncogenic signalling. To signal, oncogenic KRAS4B not only requires a sufficient nucleotide exchange, but also needs to recruit effectors by exposing its effector-binding sites while anchoring to the phospholipid bilayer where KRAS4B-mediated signalling events occur. The enzyme phosphodiesterase-δ plays an important role in sequestering KRAS4B from the cytoplasm and targeting it to cellular membranes of different cell species. In this work, we present an in silico design of a lipid-like compound that has the remarkable feature of being able to target both an oncogenic KRAS4B-G12D mutant and the phosphodiesterase-δ enzyme. This double action is accomplished by adding a lipid tail (analogous to the farnesyl group of the KRAS4B protein) to an previously known active compound (2H-1,2,4-benzothiadiazine, 3,4-dihydro-,1,1-dioxide). The proposed lipid-like molecule was found to lock KRAS4B-G12D in its GDP-bound state by adjusting the effector-binding domain to be blocked by the interface of the lipid bilayer. Meanwhile, it can tune GTP-bound KRAS4B-G12D to shift from the active orientation state to the inactive state. The proposed compound is also observed to stably accommodate itself in the prenyl-binding pocket of phosphodiesterase-δ, which impairs KRAS4B enrichment at the lipid bilayer, potentially reducing the proliferation of KRAS4B inside the cytoplasm and its anchoring at the bilayer. In conclusion, we report a potential inhibitor of KRAS4B-G12D with a lipid tail attached to a specific warhead, a compound which has not yet been considered for drugs targeting RAS mutants. Our work provides new ways to target KRAS4B-G12D and can also foster drug discovery efforts for the targeting of oncogenes of the RAS family and beyond.

Biochemical and biophysical characterization of the RAS family small GTPase protein DiRAS3

DiRAS3, also called ARHI, is a RAS (sub)family small GTPase protein that shares 50-60% sequence identity with H-, K-, and N-RAS, with substitutions in key conserved G-box motifs and a unique 34 amino acid extension at its N-terminus. Unlike the RAS proto-oncogenes, DiRAS3 exhibits tumor suppressor properties. DiRAS3 function has been studied through genetics and cell biology, but there has been a lack of understanding of the biochemical and biophysical properties of the protein, likely due to its instability and poor solubility. To overcome this solubility issue, we engineered a DiRAS3 variant (C75S/C80S), which significantly improved soluble protein expression in E. coli. Recombinant DiRAS3 was purified by Ni-NTA and size exclusion chromatography (SEC). Concentration dependence of the SEC chromatogram indicated that DiRAS3 exists in monomer-dimer equilibrium. We then produced truncations of the N-terminal (ΔN) and both (ΔNC) extensions to the GTPase domain. Unlike full-length DiRAS3, the SEC profiles showed that ΔNC is monomeric while ΔN was monomeric with aggregation, suggesting that the N and/or C-terminal tail(s) contribute to dimerization and aggregation. The 1H-15N HSQC NMR spectrum of ΔNC construct displayed well-dispersed peaks similar to spectra of other GTPase domains, which enabled us to demonstrate that DiRAS3 has a GTPase domain that can bind GDP and GTP. Taken together, we conclude that, despite the substitutions in the G-box motifs, DiRAS3 can switch between nucleotide-bound states and that the N- and C-terminal extensions interact transiently with the GTPase domain in intra- and inter-molecular fashions, mediating weak multimerization of this unique small GTPase.

Myr-Arf1 conformational flexibility at the membrane surface sheds light on the interactions with ArfGAP ASAP1

ADP-ribosylation factor 1 (Arf1) interacts with multiple cellular partners and membranes to regulate intracellular traffic, organelle structure and actin dynamics. Defining the dynamic conformational landscape of Arf1 in its active form, when bound to the membrane, is of high functional relevance and key to understanding how Arf1 can alter diverse cellular processes. Through concerted application of nuclear magnetic resonance (NMR), neutron reflectometry (NR) and molecular dynamics (MD) simulations, we show that, while Arf1 is anchored to the membrane through its N-terminal myristoylated amphipathic helix, the G domain explores a large conformational space, existing in a dynamic equilibrium between membrane-associated and membrane-distal conformations. These configurational dynamics expose different interfaces for interaction with effectors. Interaction with the Pleckstrin homology domain of ASAP1, an Arf-GTPase activating protein (ArfGAP), restricts motions of the G domain to lock it in what seems to be a conformation exposing functionally relevant regions.

Real-time monitoring of the reaction of KRAS G12C mutant specific covalent inhibitor by in vitro and in-cell NMR spectroscopy

KRAS mutations are major drivers of various cancers. Recently, allele-specific inhibitors of the KRAS G12C mutant were developed that covalently modify the thiol of Cys12, thereby trapping KRAS in an inactive GDP-bound state. To study the mechanism of action of the covalent inhibitors in both in vitro and intracellular environments, we used real-time NMR to simultaneously observe GTP hydrolysis and inhibitor binding. In vitro NMR experiments showed that the rate constant of ARS-853 modification is identical to that of GTP hydrolysis, indicating that GTP hydrolysis is the rate-limiting step for ARS-853 modification. In-cell NMR analysis revealed that the ARS-853 reaction proceeds significantly faster than that in vitro, reflecting acceleration of GTP hydrolysis by endogenous GTPase proteins. This study demonstrated that the KRAS covalent inhibitor is as effective in the cell as in vitro and that in-cell NMR is a valuable validation tool for assessing the pharmacological properties of the drug in the intracellular context.

Conformational ensemble-dependent lipid recognition and segregation by prenylated intrinsically disordered regions in small GTPases

We studied diverse prenylated intrinsically disordered regions (PIDRs) of Ras and Rho family small GTPases using long timescale atomistic molecular dynamics simulations in an asymmetric model membrane of phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. Here we show that conformational plasticity is a key determinant of lipid sorting by polybasic PIDRs and provide evidence for lipid sorting based on both headgroup and acyl chain structures. We further show that conformational ensemble-based lipid recognition is generalizable to all polybasic PIDRs, and that the sequence outside the polybasic domain (PBD) modulates the conformational plasticity, bilayer adsorption, and interactions of PIDRs with membrane lipids. Specifically, we find that palmitoylation, the ratio of basic to acidic residues, and the hydrophobic content of the sequence outside the PBD significantly impact the diversity of conformational substates and hence the extent of conformation-dependent lipid interactions. We thus propose that the PBD is required but not sufficient for the full realization of lipid sorting by prenylated PBD-containing membrane anchors, and that the membrane anchor is not only responsible for high affinity membrane binding but also directs the protein to the right target membrane where it participates in lipid sorting.

Revealing KRas4b topology on the membrane surface

KRas4b is a membrane-bound regulatory protein belonging to the family of small GTPases that function as a molecular switch, facilitating signal transduction from activated membrane receptors to intracellular pathways controlling cell growth and proliferation. Oncogenic mutations locking KRas4b in the active GTP state are responsible for nearly 85% of all Ras-driven cancers. Understanding the membrane-bound state of KRas4b is crucial for designing new therapeutic approaches targeting oncogenic KRas-driven signaling pathways. Extensive research demonstrates the significant involvement of the membrane bilayer in Ras-effector interactions, with anionic lipids playing a critical role in determining protein conformations The preferred topology of KRas4b for interacting with signaling partners has been a long-time question. Computational studies suggest a membrane-proximal conformation, while other biophysical methods like neutron reflectivity propose a membrane-distal conformation. To address these gaps, we employed FRET measurements to investigate the conformation of KRas4b. Using fully post-translationally modified KRas4b, we designed a Nanodisc based FRET assay to study KRas4b-membrane interactions. We suggest an extended conformation of KRas4b relative to the membrane surface. Measurement of FRET donor - acceptor distances reveal that a negatively charged membrane surface weakly favors closer association with the membrane surface. Our findings provide insights into the role of anionic lipids in determining the dynamic conformations of KRas4b and shed light on the predominant conformation of its topology on lipid headgroups.

Lipid nanodiscs as a template for high-resolution cryo-EM structures of peripheral membrane proteins

Peripheral membrane proteins are ubiquitous throughout cell biology and are required for a variety of cellular processes such as signal transduction, membrane trafficking, and autophagy. Transient binding to the membrane has a profound impact on protein function, serving to induce conformational changes and alter biochemical and biophysical parameters by increasing the local concentration of factors and restricting diffusion to two dimensions. Despite the centrality of the membrane in serving as a template for cell biology, there are few reported high-resolution structures of peripheral membrane proteins bound to the membrane. We analyzed the utility of lipid nanodiscs to serve as a template for cryo-EM analysis of peripheral membrane proteins. We tested a variety of nanodiscs and we report a 3.3 Å structure of the AP2 clathrin adaptor complex bound to a 17-nm nanodisc, with sufficient resolution to visualize a bound lipid head group. Our data demonstrate that lipid nanodiscs are amenable to high-resolution structure determination of peripheral membrane proteins and provide a framework for extending this analysis to other systems.

Small-Molecule Inhibition of KRAS through Conformational Selection

Mutations in KRAS account for about 20% of human cancers. Despite the major progress in recent years toward the development of KRAS inhibitors, including the discovery of covalent inhibitors of the G12C KRAS variant for the treatment of non-small-cell lung cancer, much work remains to be done to discover broad-acting inhibitors to treat many other KRAS-driven cancers. In a previous report, we showed that a 308.4 Da small-molecule ligand [(2R)-2-(N'-(1H-indole-3-carbonyl)hydrazino)-2-phenyl-acetamide] binds to KRAS with low micro-molar affinity [Chem. Biol. Drug Des.2019; 94(2):1441-1456]. Binding of this ligand, which we call ACA22, to the p1 pocket of KRAS and its interactions with residues at beta-strand 1 and the switch loops have been supported by data from nuclear magnetic resonance spectroscopy and microscale thermophoresis experiments. However, the inhibitory potential of the compound was not demonstrated. Here, we show that ACA22 inhibits KRAS-mediated signal transduction in cells expressing wild type (WT) and G12D mutant KRAS and reduces levels of guanosine triphosphate-loaded WT KRAS more effectively than G12D KRAS. We ruled out the direct effect on nucleotide exchange or effector binding as possible mechanisms of inhibition using a variety of biophysical assays. Combining these observations with binding data that show comparable affinities of the compound for the active and inactive forms of the mutant but not the WT, we propose conformational selection as a possible mechanism of action of ACA22.

Conformational ensemble dependent lipid recognition and segregation by prenylated intrinsically disordered regions in small GTPases

We studied diverse prenylated intrinsically disordered regions (PIDRs) of Ras and Rho family small GTPases using long timescale atomistic molecular dynamics simulations in an asymmetric model membrane of phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. We show that conformational plasticity is a key determinant of lipid sorting by polybasic PIDRs and provide evidence for lipid sorting based on both headgroup and acyl chain structures. We further show that conformational ensemble-based lipid recognition is generalizable to all polybasic PIDRs, and that the sequence outside the polybasic domain (PBD) modulates the conformational plasticity, bilayer adsorption, and interactions of PIDRs with membrane lipids. Specifically, we found that palmitoylation, the ratio of basic to acidic residues, and the hydrophobic content of the sequence outside the PBD significantly impact the diversity of conformational substates and hence the extent of conformation-dependent lipid interactions. We thus propose that the PBD is required but not sufficient for the full realization of lipid sorting by prenylated PBD-containing membrane anchors, and that the membrane anchor is not only responsible for high affinity membrane binding but also directs the protein to the right target membrane where it participates in lipid sorting.

Direct Observation of Membrane-Associated H-Ras in the Native Cellular Environment by In-Cell 19F-NMR Spectroscopy

Ras acts as a molecular switch to control intracellular signaling on the plasma membrane (PM). Elucidating how Ras associates with PM in the native cellular environment is crucial for understanding its control mechanism. Here, we used in-cell nuclear magnetic resonance (NMR) spectroscopy combined with site-specific ¹⁹F-labeling to explore the membrane-associated states of H-Ras in living cells. The site-specific incorporation of p-trifluoromethoxyphenylalanine (OCF₃Phe) at three different sites of H-Ras, i.e., Tyr32 in switch I, Tyr96 interacting with switch II, and Tyr157 on helix α5, allowed the characterization of their conformational states depending on the nucleotide-bound states and an oncogenic mutational state. Exogenously delivered ¹⁹F-labeled H-Ras protein containing a C-terminal hypervariable region was assimilated via endogenous membrane-trafficking, enabling proper association with the cell membrane compartments. Despite poor sensitivity of the in-cell NMR spectra of membrane-associated H-Ras, the Bayesian spectral deconvolution identified distinct signal components on three ¹⁹F-labeled sites, thus offering the conformational multiplicity of H-Ras on the PM. Our study may be helpful in elucidating the atomic-scale picture of membrane-associated proteins in living cells.

RAS-targeted cancer therapy: Advances in drugging specific mutations

Rat sarcoma (RAS), as a frequently mutated oncogene, has been studied as an attractive target for treating RAS-driven cancers for over four decades. However, it is until the recent success of kirsten-RAS (KRAS)G12C inhibitor that RAS gets rid of the title "undruggable". It is worth noting that the therapeutic effect of KRASG12C inhibitors on different RAS allelic mutations or even different cancers with KRASG12C varies significantly. Thus, deep understanding of the characteristics of each allelic RAS mutation will be a prerequisite for developing new RAS inhibitors. In this review, the structural and biochemical features of different RAS mutations are summarized and compared. Besides, the pathological characteristics and treatment responses of different cancers carrying RAS mutations are listed based on clinical reports. In addition, the development of RAS inhibitors, either direct or indirect, that target the downstream components in RAS pathway is summarized as well. Hopefully, this review will broaden our knowledge on RAS-targeting strategies and trigger more intensive studies on exploiting new RAS allele-specific inhibitors.

Quantitative Paramagnetic NMR-Based Analysis of Protein Orientational Dynamics on Membranes: Dissecting the KRas4B-Membrane Interactions

Peripheral membrane proteins can adopt distinct orientations on the surfaces of lipid bilayers that are often short-lived and challenging to characterize by conventional experimental methods. Here we describe a robust approach for mapping protein orientational landscapes through quantitative interpretation of paramagnetic relaxation enhancement (PRE) data arising from membrane mimetics with spin-labeled lipids. Theoretical analysis, followed by experimental verification, reveals insights into the distinct properties of the PRE observables that are generally distorted in the case of stably membrane-anchored proteins. To suppress the artifacts, we demonstrate that undistorted Γ₂ values can be obtained via transient membrane anchoring, based on which a computational framework is established for deriving accurate orientational ensembles obeying Boltzmann statistics. Application of the approach to KRas4B, a classical peripheral membrane protein whose orientations are critical for its functions and drug design, reveals four distinct orientational states that are close but not identical to those reported previously. Similar orientations are also found for a truncated KRas4B without the hypervariable region (HVR) that can sample a broader range of orientations, suggesting a confinement role of the HVR geometrically prohibiting severe tilting. Comparison of the KRas4B Γ₂ rates measured using nanodiscs containing different types of anionic lipids reveals identical Γ₂ patterns for the G-domain but different ones for the HVR, indicating only the latter is able to selectively interact with anionic lipids.

Eliminating oncogenic RAS: back to the future at the drawing board

RAS drug development has made enormous strides in the past ten years, with the first direct KRAS inhibitor being approved in 2021. However, despite the clinical success of covalent KRAS-G12C inhibitors, we are immediately confronted with resistances as commonly found with targeted drugs. Previously believed to be undruggable due to its lack of obvious druggable pockets, a couple of new approaches to hit this much feared oncogene have now been carved out. We here concisely review these approaches to directly target four druggable sites of RAS from various angles. Our analysis focuses on the lessons learnt during the development of allele-specific covalent and non-covalent RAS inhibitors, the potential of macromolecular binders to facilitate the discovery and validation of targetable sites on RAS and finally an outlook on a future that may engage more small molecule binders to become drugs. We foresee that the latter could happen mainly in two ways: First, non-covalent small molecule inhibitors may be derived from the development of covalent binders. Second, reversible small molecule binders could be utilized for novel targeting modalities, such as degraders of RAS. Provided that degraders eliminate RAS by recruiting differentially expressed E3-ligases, this approach could enable unprecedented tissue- or developmental stage-specific destruction of RAS with potential advantages for on-target toxicity. We conclude that novel creative ideas continue to be important to exterminate RAS in cancer and other RAS pathway-driven diseases, such as RASopathies.

Cellular Labeling of Phosphatidylserine Using Clickable Serine Probes

Phosphatidylserine (PS) is a key lipid that plays important roles in disease-related biological processes, and therefore, the means to track PS in live cells are invaluable. Herein, we describe the metabolic labeling of PS in Saccharomyces cerevisiae cells using analogues of serine, a PS precursor, derivatized with azide moieties at either the amino (N-l-SerN3) or carbonyl (C-l-SerN3) groups. The conservative click tag modification enabled these compounds to infiltrate normal lipid biosynthetic pathways, thereby producing tagged PS molecules as supported by mass spectrometry studies, thin-layer chromatography (TLC) analysis, and further derivatization with fluorescent reporters via click chemistry to enable imaging in yeast cells. This approach shows strong prospects for elucidating the complex biosynthetic and trafficking pathways involving PS.

Designer Nanodiscs to Probe and Reprogram Membrane Biology in Synapses

Signal transduction at the synapse is mediated by a variety of protein-lipid interactions, which are vital for the spatial and temporal regulation of synaptic vesicle biogenesis, neurotransmitter release, and postsynaptic receptor activation. Therefore, our understanding of synaptic transmission cannot be completed until the elucidation of these critical protein-lipid interactions. On this front, recent advances in nanodiscs have vastly expanded our ability to probe and reprogram membrane biology in synapses. Here, we summarize the progress of the nanodisc toolbox and discuss future directions in this exciting field.

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α2 and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α4-α5 KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).

Exploring CRD mobility during RAS/RAF engagement at the membrane

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Molecular and clinical profile of patients referred as Noonan or Noonan-like syndrome in Greece: a cohort of 86 patients

Noonan syndrome (NS) is an autosomal dominant disorder characterized by clinical and genetic heterogeneity. It belongs to a wider group of pathologies, known as Rasopathies, due to the implication of genes encoding components of the Ras/MAPK signalling pathway. Recording the genetic alterations across populations helps assessing specific features to specific genes which is essential for better disease's recognition, prognosis and monitoring. Herein, we report the clinical and molecular data of a Greek cohort comprising of 86 NS or NS-like patients admitted at a single tertiary Centre in Athens, Greece. The analysis was performed using Sanger and next-generation sequencing, comprising 14 different genes. The mutational rates of the confirmed NS-associated genes in the Greek NS population are as follows: PTPN11 32.5%; RIT1 5.8%; SOS1 4.7%; BRAF 1.2%; CBL 1.2%; KRAS 1.2%; MAP2K1 1.2%; RAF1 1.2%; SHOC2 1.2%, corresponding to 50% of positivity in total NS population. The genotype-phenotype analysis showed statistically significant differences in craniofacial dysmorphisms (p = 0.005) and pulmonary valve stenosis (PS) (p < 0.001) frequencies between patients harbouring a pathogenic variant and patients without pathogenic variant in any of the tested genes. Patients with at least a pathogenic variant had 6.71 times greater odds to develop PS compared to pathogenic variant-negative patients (OR = 6.71, 95%; CI = (2.61, 17.27)). PTPN11 positive patients showed higher frequency of epicanthal folds (p = 0.004), ptosis (p = 0.001) and coarseness (p = 0.001) and lower frequency of neurological findings (p = 0.006), compared to patients carrying pathogenic variants in other genes.

CONCLUSION

Craniofacial dysmorphism and PS prevail among pathogenic variant positive compared to pathogenic variant negative NS and NS-like patients while neurological defects are less common in PTPN11-affected NS patients compared to patients harbouring pathogenic variants in other genes. The significant prevalence of the Ras/MAPK pathogenic variants (17.4%), other than PTPN11, in Greek NS patients, highlights the necessity of a wider spectrum of molecular diagnosis.

WHAT IS KNOWN

• Noonan syndrome (NS) has been associated with pathogenic variants in molecules-components of the Ras/MAPK pathway. • Clinical and genetic description of NS patients worldwide helps establishing personalized monitoring.

WHAT IS NEW

• NS and NS-like mutational rate in Greece reaches 50% with pathogenic variants identified mostly in PTPN11 (32.5%), RIT1 (6%) and SOS1 (4.7%) genes. • The risk for pulmonary stenosis increases 6.71-fold in NS patients with a pathogenic variant compared to patients without genetic alterations.

The current state of the art and future trends in RAS-targeted cancer therapies

Despite being the most frequently altered oncogenic protein in solid tumours, KRAS has historically been considered ‘undruggable’ owing to a lack of pharmacologically targetable pockets within the mutant isoforms. However, improvements in drug design have culminated in the development of inhibitors that are selective for mutant KRAS in its active or inactive state. Some of these inhibitors have proven efficacy in patients with KRAS^G12C-mutant cancers and have become practice changing. The excitement associated with these advances has been tempered by drug resistance, which limits the depth and/or duration of responses to these agents. Improvements in our understanding of RAS signalling in cancer cells and in the tumour microenvironment suggest the potential for several novel combination therapies, which are now being explored in clinical trials. Herein, we provide an overview of the RAS pathway and review the development and current status of therapeutic strategies for targeting oncogenic RAS, as well as their potential to improve outcomes in patients with RAS-mutant malignancies. We then discuss challenges presented by resistance mechanisms and strategies by which they could potentially be overcome.

The RAS oncogenes are among the most common drivers of tumour development and progression but have historically been considered undruggable. The development of direct KRAS inhibitors has changed this paradigm, although currently clinical use of these novel therapeutics is limited to a select subset of patients, and intrinsic or acquired resistance presents an inevitable challenge to cure. Herein, the authors provide an overview of the RAS pathway in cancer and review the ongoing efforts to develop effective therapeutic strategies for RAS-mutant cancers. They also discuss the current understanding of mechanisms of resistance to direct KRAS inhibitors and strategies by which they might be overcome.

Owing to intrinsic and extrinsic factors, KRAS and other RAS isoforms have until recently been impervious to targeting with small-molecule inhibitors.

Inhibitors of the KRAS^G12C variant constitute a potential breakthrough in the treatment of many cancer types, particularly non-small-cell lung cancer, for which such an agent has been approved by the FDA.

Several forms of resistance to KRAS inhibitors have been defined, including primary, adaptive and acquired resistance; these resistance mechanisms are being targeted in studies that combine KRAS inhibitors with inhibitors of horizontal or vertical signalling pathways.

Mutant KRAS has important effects on the tumour microenvironment, including the immunological milieu; these effects must be considered to fully understand resistance to KRAS inhibitors and when designing novel treatment strategies.

Conformational Dynamics Allows Sampling of an "Active-like" State by Oncogenic K-Ras-GDP

Mutations in K-Ras GTPase replacing Gly12 with either Asp or Val are common in cancer. These mutations decelerate intrinsic and catalyzed GTP hydrolysis, leading to accumulation of K-Ras-GTP in cells. Signaling cascades initiated by K-Ras-GTP promote cell proliferation, survival, and invasion. Despite functional differences between the most frequent G12D mutation and the most aggressive and chemotherapy resistant G12V mutation, their long-suspected distinct structural features remain elusive. Using NMR, X-ray structures, and computational methods, we found that oncogenic mutants of K-Ras4B, the predominant splice variant of K-Ras, exhibit distinct conformational dynamics when GDP-bound, visiting the "active-like" conformational state similar to the one observed in GTP-bound K-Ras. This behavior distinguishes G12V from wild type and G12D K-Ras4B-GDP. The likely reason is hydrophobic interactions between the aliphatic sidechain of V12 and the Switch II region of K-Ras4B^G12V-GDP, which are distinct in K-Ras4B^G12D-GDP. In the X-ray structures, crystal contacts reduce the dynamics of the sidechain at position 12 by stabilizing the Switch I region of the protein. This explains why structural differences between G12V and G12D K-Ras have yet not been reported. Together, our results suggest a previously unknown mechanism of K-Ras activation. This mechanism relies on conformational dynamics caused by specific oncogenic mutations in the GDP-bound state. Our findings also imply that the therapeutic strategies decreasing the level of K-Ras-GTP by interfering with nucleotide exchange or by expediting GTP hydrolysis may work differently in different oncogenic mutants.

1 H, 15 N and 13 C resonance assignments of the Q61H mutant of human KRAS bound to GDP

KRAS proteins are small GTPases binding to the cell membrane and playing important roles in signal transduction. KRAS proteins form complexes with GTP and GDP to result in active and inactive conformations favouring interactions with different proteins. Mutations in KRAS have impact on the GTPase activity and some mutants are related to certain types of cancers. In addition to mutation at position 12, the Q61H mutant is also identified as an oncogenic mutant. Here, we describe resonance assignment for Q61H mutant of human KRAS-4B. A construct containing 1-169 residues of KRAS with a point mutation at position 61 (Q to H) was made for solution NMR studies. The backbone and some side chain resonance assignments were obtained using conventional multi-dimensional experiments. The secondary structures were analysed based on the assigned residues. As NMR is a powerful tool in probing target and ligand interactions, the assignment will be useful for later compound binding studies.

Classification of KRAS-Activating Mutations and the Implications for Therapeutic Intervention

Members of the family of RAS proto-oncogenes, discovered just over 40 years ago, were among the first cancer-initiating genes to be discovered. Of the three RAS family members, KRAS is the most frequently mutated in human cancers. Despite intensive biological and biochemical study of RAS proteins over the past four decades, we are only now starting to devise therapeutic strategies to target their oncogenic properties. Here, we highlight the distinct biochemical properties of common and rare KRAS alleles, enabling their classification into functional subtypes. We also discuss the implications of this functional classification for potential therapeutic avenues targeting mutant subtypes.

SIGNIFICANCE

Efforts in the recent past to inhibit KRAS oncogenicity have focused on kinases that function in downstream signal transduction cascades, although preclinical successes have not translated to patients with KRAS-mutant cancer. Recently, clinically effective covalent inhibitors of KRASG12C have been developed, establishing two principles that form a foundation for future efforts. First, KRAS is druggable. Second, each mutant form of KRAS is likely to have properties that make it uniquely druggable.

A Structure is Worth a Thousand Words: New Insights for RAS and RAF Regulation

The RAS GTPases are frequently mutated in human cancer, with KRAS being the predominant tumor driver. For many years, it has been known that the structure and function of RAS are integrally linked, as structural changes induced by GTP binding or mutational events determine the ability of RAS to interact with regulators and effectors. Recently, a wealth of information has emerged from structures of specific KRAS mutants and from structures of multiprotein complexes containing RAS and/or RAF, an essential effector of RAS. These structures provide key insights regarding RAS and RAF regulation as well as promising new strategies for therapeutic intervention.

SIGNIFICANCE

The RAS GTPases are major drivers of tumorigenesis, and for RAS proteins to exert their full oncogenic potential, they must interact with the RAF kinases to initiate ERK cascade signaling. Although binding to RAS is typically a prerequisite for RAF to become an activated kinase, determining the molecular mechanisms by which this interaction results in RAF activation has been a challenging task. A major advance in understanding this process and RAF regulation has come from recent structural studies of various RAS and RAF multiprotein signaling complexes, revealing new avenues for drug discovery.

Identification of functional substates of KRas during GTP hydrolysis with enhanced sampling simulations

As the hub of major signaling pathways, Ras proteins are implicated in 19% of tumor-caused cancers due to perturbations in their conformational and/or catalytic properties. Despite numerous studies, the functions of the conformational substates for the most important isoform, KRas, remain elusive. In this work, we perform an extensive simulation analysis on the conformational landscape of KRas in its various chemical states during the GTP hydrolysis cycle: the reactant state KRasGTP·Mg²⁺, the intermediate state KRasGDP·Pi·Mg²⁺ and the product state KRasGDP·Mg²⁺. The results from enhanced sampling simulations reveal that State 1 of KRasGTP·Mg²⁺ has multiple stable substates in solution, one of which might account for interacting with GEFs. State 2 of KRasGTP·Mg²⁺ features two substates "Tyr32_in" and "Tyr32_out", which are poised to interact with effectors and GAPs, respectively. For the intermediate state KRasGDP·Pi·Mg²⁺, Gln61 and Pi are found to assume a broad set of conformations, which might account for the weak oncogenic effect of Gln61 mutations in KRas in contrast to the situation in HRas and NRas. Finally, the product state KRasGDP·Mg²⁺ has more than two stable substates in solution, pointing to a conformation-selection mechanism for complexation with GEFs. Based on these results, some specific inhibition strategies for targeting the binding sites of the high-energy substates of KRas during GTP hydrolysis are discussed.

A saturation-mutagenesis analysis of the interplay between stability and activation in Ras

Cancer mutations in Ras occur predominantly at three hotspots: Gly 12, Gly 13, and Gln 61. Previously, we reported that deep mutagenesis of H-Ras using a bacterial assay identified many other activating mutations (Bandaru et al., 2017). We now show that the results of saturation mutagenesis of H-Ras in mammalian Ba/F3 cells correlate well with the results of bacterial experiments in which H-Ras or K-Ras are co-expressed with a GTPase-activating protein (GAP). The prominent cancer hotspots are not dominant in the Ba/F3 data. We used the bacterial system to mutagenize Ras constructs of different stabilities and discovered a feature that distinguishes the cancer hotspots. While mutations at the cancer hotspots activate Ras regardless of construct stability, mutations at lower-frequency sites (e.g. at Val 14 or Asp 119) can be activating or deleterious, depending on the stability of the Ras construct. We characterized the dynamics of three non-hotspot activating Ras mutants by using NMR to monitor hydrogen-deuterium exchange (HDX). These mutations result in global increases in HDX rates, consistent with destabilization of Ras. An explanation for these observations is that mutations that destabilize Ras increase nucleotide dissociation rates, enabling activation by spontaneous nucleotide exchange. A further stability decrease can lead to insufficient levels of folded Ras - and subsequent loss of function. In contrast, the cancer hotspot mutations are mechanism-based activators of Ras that interfere directly with the action of GAPs. Our results demonstrate the importance of GAP surveillance and protein stability in determining the sensitivity of Ras to mutational activation.

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-μs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

Predicting the conformational variability of oncogenic GTP-bound G12D mutated KRas-4B proteins at zwitterionic model cell membranes

KRas proteins are the largest family of mutated Ras isoforms, participating in a wide variety of cancers. Due to their importance, large effort is being carried out on drug development by small-molecule inhibitors. However, understanding protein conformational variability remains a challenge in drug discovery. In the case of the Ras family, their multiple conformational states can affect the binding of potential drug inhibitors. To overcome this challenge, we propose a computational framework based on combined all-atom Molecular Dynamics and Metadynamics simulations in order to accurately access conformational variants of the target protein. We tested the methodology using a G12D mutated GTP bound oncogenic KRas-4B protein located at the interface of a DOPC/DOPS/cholesterol model anionic cell membrane. Two main orientations of KRas-4B at the anionic membrane have been determined. The corresponding torsional angles are taken as reliable reaction coordinates so that free-energy landscapes are obtained by well-tempered metadynamics simulations, revealing local and global minima of the free-energy hypersurface and unveiling reactive paths of the system between the two preferential orientations. We have observed that GTP-binding to KRas-4B has huge influence on the stabilisation of the protein and it can potentially help to open Switch I/II druggable pockets, lowering energy barriers between stable states and resulting in cumulative conformers of KRas-4B. This may highlight new opportunities for targeting the unique meta-stable states through the design of new efficient drugs.

Membrane Composition and Raf[CRD]-Membrane Attachment Are Driving Forces for K-Ras4B Dimer Stability

Ras proteins are membrane-anchored GTPases that regulate key cellular signaling networks. It has been recently shown that different anionic lipid types can affect the properties of Ras in terms of dimerization/clustering on the cell membrane. To understand the effects of anionic lipids on key spatiotemporal properties of dimeric K-Ras4B, we perform all-atom molecular dynamics simulations of the dimer K-Ras4B in the presence and absence of Raf[RBD/CRD] effectors on two model anionic lipid membranes: one containing 78% mol DOPC, 20% mol DOPS, and 2% mol PIP2 and another one with enhanced concentration of anionic lipids containing 50% mol DOPC, 40% mol DOPS, and 10% mol PIP2. Analysis of our results unveils the orientational space of dimeric K-Ras4B and shows that the stability of the dimer is enhanced on the membrane containing a high concentration of anionic lipids in the absence of Raf effectors. This enhanced stability is also observed in the presence of Raf[RBD/CRD] effectors although it is not influenced by the concentration of anionic lipids in the membrane, but rather on the ability of Raf[CRD] to anchor to the membrane. We generate dominant K-Ras4B conformations by Markov state modeling and yield the population of states according to the K-Ras4B orientation on the membrane. For the membrane containing anionic lipids, we observe correlations between the diffusion of K-Ras4B and PIP2 and anchoring of anionic lipids to the Raf[CRD] domain. We conclude that the presence of effectors with the Raf[CRD] domain anchoring on the membrane as well as the membrane composition both influence the conformational stability of the K-Ras4B dimer, enabling the preservation of crucial interface interactions.

Ras Multimers on the Membrane: Many Ways for a Heart-to-Heart Conversation

Formation of Ras multimers, including dimers and nanoclusters, has emerged as an exciting, new front of research in the 'old' field of Ras biomedicine. With significant advances made in the past few years, we are beginning to understand the structure of Ras multimers and, albeit preliminary, mechanisms that regulate their formation in vitro and in cells. Here we aim to synthesize the knowledge accrued thus far on Ras multimers, particularly the presence of multiple globular (G-) domain interfaces, and discuss how membrane nanodomain composition and structure would influence Ras multimer formation. We end with some general thoughts on the potential implications of Ras multimers in basic and translational biology.

Determinants of Membrane Orientation Dynamics in Lipid-Modified Small GTPases

The transient membrane engagement and reorientation of the soluble catalytic domain of Ras proteins has emerged as an important modulator of their functions. However, there has been limited information on whether this phenomenon is applicable to other members of the Ras superfamily. To address this issue, we conducted long-time-scale atomistic molecular dynamics simulations (55 μs aggregate simulation time) on representatives of the Ras, Rho, and Arf family proteins that differ in sequence, lipid modification, and the rigidity of the linker between the lipid anchor and the catalytic G-domain. The results show that the concept of membrane reorientation is generalizable to most but not all members of the Ras superfamily. Specifically, C-terminally prenylated small GTPases that are anchored to membranes via a single flexible linker adopt multiple orientations, whereas those that are N-terminally myristoylated and harbor a rigid linker experience limited orientational dynamics. Combined with published reports on Ras proteins, these observations provide insights into the common principles and determinants of the orientational dynamics of lipidated small GTPases on membrane surfaces and offer new ways of thinking about the regulation and druggability of the Ras superfamily proteins.

Early-stage structure-based drug discovery for small GTPases by NMR spectroscopy

Small GTPase or Ras superfamily, including Ras, Rho, Rab, Ran and Arf, are fundamental in regulating a wide range of cellular processes such as growth, differentiation, migration and apoptosis. They share structural and functional similarities for binding guanine nucleotides and hydrolyzing GTP. Dysregulations of Ras proteins are involved in the pathophysiology of multiple human diseases, however there is still a stringent need for effective treatments targeting these proteins. For decades, small GTPases were recognized as 'undruggable' targets due to their complex regulatory mechanisms and lack of deep pockets for ligand binding. NMR has been critical in deciphering the structural and dynamic properties of the switch regions that are underpinning molecular switch functions of small GTPases, which pave the way for developing new effective inhibitors. The recent progress of drug or lead molecule development made for small GTPases profoundly delineated how modern NMR techniques reshape the field of drug discovery. In this review, we will summarize the progress of structural and dynamic studies of small GTPases, the NMR techniques developed for structure-based drug screening and their applications in early-stage drug discovery for small GTPases.

Machine learning-driven multiscale modeling reveals lipid-dependent dynamics of RAS signaling proteins

RAS is a signaling protein associated with the cell membrane that is mutated in up to 30% of human cancers. RAS signaling has been proposed to be regulated by dynamic heterogeneity of the cell membrane. Investigating such a mechanism requires near-atomistic detail at macroscopic temporal and spatial scales, which is not possible with conventional computational or experimental techniques. We demonstrate here a multiscale simulation infrastructure that uses machine learning to create a scale-bridging ensemble of over 100,000 simulations of active wild-type KRAS on a complex, asymmetric membrane. Initialized and validated with experimental data (including a new structure of active wild-type KRAS), these simulations represent a substantial advance in the ability to characterize RAS-membrane biology. We report distinctive patterns of local lipid composition that correlate with interfacially promiscuous RAS multimerization. These lipid fingerprints are coupled to RAS dynamics, predicted to influence effector binding, and therefore may be a mechanism for regulating cell signaling cascades.

Lipid Nanodiscs for High-Resolution NMR Studies of Membrane Proteins

Membrane proteins (MPs) play essential roles in numerous cellular processes. Because around 70% of the currently marketed drugs target MPs, a detailed understanding of their structure, binding properties, and functional dynamics in a physiologically relevant environment is crucial for a more detailed understanding of this important protein class. We here summarize the benefits of using lipid nanodiscs for NMR structural investigations and provide a detailed overview of the currently used lipid nanodisc systems as well as their applications in solution-state NMR. Despite the increasing use of other structural methods for the structure determination of MPs in lipid nanodiscs, solution NMR turns out to be a versatile tool to probe a wide range of MP features, ranging from the structure determination of small to medium-sized MPs to probing ligand and partner protein binding as well as functionally relevant dynamical signatures in a lipid nanodisc setting. We will expand on these topics by discussing recent NMR studies with lipid nanodiscs and work out a key workflow for optimizing the nanodisc incorporation of an MP for subsequent NMR investigations. With this, we hope to provide a comprehensive background to enable an informed assessment of the applicability of lipid nanodiscs for NMR studies of a particular MP of interest.

Oncogenic KRAS G12D mutation promotes dimerization through a second, phosphatidylserine-dependent interface: a model for KRAS oligomerization

KRAS forms transient dimers and higher-order multimers (nanoclusters) on the plasma membrane, which drive MAPK signaling and cell proliferation. KRAS is a frequently mutated oncogene, and while it is well known that the most prevalent mutation, G12D, impairs GTP hydrolysis, thereby increasing KRAS activation, G12D has also been shown to enhance nanoclustering. Elucidating structures of dynamic KRAS assemblies on a membrane has been challenging, thus we have refined our NMR approach that uses nanodiscs to study KRAS associated with membranes. We incorporated paramagnetic relaxation enhancement (PRE) titrations and interface mutagenesis, which revealed that, in addition to the symmetric ‘α–α’ dimerization interface shared with wild-type KRAS, the G12D mutant also self-associates through an asymmetric ‘α–β’ interface. The ‘α–β’ association is dependent on the presence of phosphatidylserine lipids, consistent with previous reports that this lipid promotes KRAS self-assembly on the plasma membrane in cells. Experiments using engineered mutants to spoil each interface, together with PRE probes attached to the membrane or free in solvent, suggest that dimerization through the primary ‘α–α’ interface releases β interfaces from the membrane promoting formation of the secondary ‘α–β’ interaction, potentially initiating nanoclustering. In addition, the small molecule BI-2852 binds at a β–β interface, stabilizing a new dimer configuration that outcompetes native dimerization and blocks the effector-binding site. Our data indicate that KRAS self-association involves a delicately balanced conformational equilibrium between transient states, which is sensitive to disease-associated mutation and small molecule inhibitors. The methods developed here are applicable to biologically important transient interactions involving other membrane-associated proteins.

Studies of membrane-dependent dimerization of KRAS on nanodiscs using paramagnetic NMR titrations and mutagenesis revealed a novel asymmetric ‘α–β’ interface that provides a potential mechanism for the enhanced assembly of KRAS–G12D nanoclusters.

A structural model of a Ras-Raf signalosome

The protein K-Ras functions as a molecular switch in signaling pathways regulating cell growth. In the human mitogen-activated protein kinase (MAPK) pathway, which is implicated in many cancers, multiple K-Ras proteins are thought to assemble at the cell membrane with Ras effector proteins from the Raf family. Here we propose an atomistic structural model for such an assembly. Our starting point was an asymmetric guanosine triphosphate-mediated K-Ras dimer model, which we generated using unbiased molecular dynamics simulations and verified with mutagenesis experiments. Adding further K-Ras monomers in a head-to-tail fashion led to a compact helical assembly, a model we validated using electron microscopy and cell-based experiments. This assembly stabilizes K-Ras in its active state and presents composite interfaces to facilitate Raf binding. Guided by existing experimental data, we then positioned C-Raf, the downstream kinase MEK1 and accessory proteins (Galectin-3 and 14-3-3σ) on and around the helical assembly. The resulting Ras-Raf signalosome model offers an explanation for a large body of data on MAPK signaling.

Membrane-bound KRAS approximates an entropic ensemble of configurations

KRAS4B is a membrane-anchored signaling protein and a primary target in cancer research. Predictions from molecular dynamics simulations that have previously shaped our mechanistic understanding of KRAS signaling disagree with recent experimental results from neutron reflectometry, NMR, and thermodynamic binding studies. To gain insight into these discrepancies, we compare this body of biophysical data to back-calculated experimental results from a series of molecular simulations that implement different subsets of molecular interactions. Our results show that KRAS4B approximates an entropic ensemble of configurations at model membranes containing 30% phosphatidylserine lipids, which is not significantly shaped by interactions between the globular G-domain of KRAS4B and the lipid membrane. These findings revise our understanding of KRAS signaling and promote a model in which the protein samples the accessible conformational space in a near-uniform manner while being available to bind to effector proteins.