Oncogenic K-Ras4B Dimerization Enhances Downstream Mitogen-activated Protein Kinase Signaling

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.

Allostery: Allosteric Cancer Drivers and Innovative Allosteric Drugs

Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.

Allostery, and how to define and measure signal transduction

Here we ask: What is productive signaling? How to define it, how to measure it, and most of all, what are the parameters that determine it? Further, what determines the strength of signaling from an upstream to a downstream node in a specific cell? These questions have either not been considered or not entirely resolved. The requirements for the signal to propagate downstream to activate (repress) transcription have not been considered either. Yet, the questions are pivotal to clarify, especially in diseases such as cancer where determination of signal propagation can point to cell proliferation and to emerging drug resistance, and to neurodevelopmental disorders, such as RASopathy, autism, attention-deficit/hyperactivity disorder (ADHD), and cerebral palsy. Here we propose a framework for signal transduction from an upstream to a downstream node addressing these questions. Defining cellular processes, experimentally measuring them, and devising powerful computational AI-powered algorithms that exploit the measurements, are essential for quantitative science.

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.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Raf promotes dimerization of the Ras G-domain with increased allosteric connections

Ras dimerization is critical for Raf activation. Here we show that the Ras binding domain of Raf (Raf-RBD) induces robust Ras dimerization at low surface densities on supported lipid bilayers and, to a lesser extent, in solution as observed by size exclusion chromatography and confirmed by SAXS. Community network analysis based on molecular dynamics simulations shows robust allosteric connections linking the two Raf-RBD D113 residues located in the Galectin scaffold protein binding site of each Raf-RBD molecule and 85 Å apart on opposite ends of the dimer complex. Our results suggest that Raf-RBD binding and Ras dimerization are concerted events that lead to a high-affinity signaling complex at the membrane that we propose is an essential unit in the macromolecular assembly of higher order Ras/Raf/Galectin complexes important for signaling through the Ras/Raf/MEK/ERK pathway.

Ras isoform-specific expression, chromatin accessibility, and signaling

The anchorage of Ras isoforms in the membrane and their nanocluster formations have been studied extensively, including their detailed interactions, sizes, preferred membrane environments, chemistry, and geometry. However, the staggering challenge of their epigenetics and chromatin accessibility in distinct cell states and types, which we propose is a major factor determining their specific expression, still awaits unraveling. Ras isoforms are distinguished by their C-terminal hypervariable region (HVR) which acts in intracellular transport, regulation, and membrane anchorage. Here, we review some isoform-specific activities at the plasma membrane from a structural dynamic standpoint. Inspired by physics and chemistry, we recognize that understanding functional specificity requires insight into how biomolecules can organize themselves in different cellular environments. Within this framework, we suggest that isoform-specific expression may largely be controlled by the chromatin density and physical compaction, which allow (or curb) access to "chromatinized DNA." Genes are preferentially expressed in tissues: proteins expressed in pancreatic cells may not be equally expressed in lung cells. It is the rule-not an exception, and it can be at least partly understood in terms of chromatin organization and accessibility state. Genes are expressed when they can be sufficiently exposed to the transcription machinery, and they are less so when they are persistently buried in dense chromatin. Notably, chromatin accessibility can similarly determine expression of drug resistance genes.

Mechanistic Differences of Activation of Rac1P29S and Rac1A159V

Rac1 is a small GTPase that plays key roles in actin reorganization, cell motility, and cell survival/growth as well as in various cancer types and neurodegenerative diseases. Similar to other Ras superfamily GTPases, Rac1 switches between active GTP-bound and inactive GDP-bound states. Switch I and II regions open and close during GDP/GTP exchange. P29S and A159V (paralogous to K-RasA146) mutations are the two most common somatic mutations of Rac1. Rac1P29S is a known hotspot for melanoma, whereas Rac1A159V most commonly occurs in head and neck cancer. To investigate how these substitutions induce the Rac1 dynamics, we used atomistic molecular dynamics simulations on the wild-type Rac1 and two mutant systems (P29S and A159V) in the GTP bound state, and on the wild-type Rac1 and P29S mutated system in the GDP bound state. Here, we show that P29S and A159V mutations activate Rac1 with different mechanisms. In Rac1P29S-GTP, the substitution increases the flexibility of Switch I based on RMSF and dihedral angle calculations and leads to an open conformation. We propose that the open Switch I conformation is one of the underlying reasons for rapid GDP/GTP exchange of Rac1P29S. On the other hand, in Rac1A159V-GTP, some of the contacts of the guanosine ring of GTP with Rac1 are temporarily lost, enabling the guanosine ring to move toward Switch I and subsequently close the switch. Rac1A159V-GTP adopts a Ras state 2 like conformation, where both switch regions are in closed conformation and Thr35 forms a hydrogen bond with the nucleotide.

Diffusion-based determination of protein homodimerization on reconstituted membrane surfaces

The transient interactions between cellular components, particularly on membrane surfaces, are critical in the proper function of many biochemical reactions. For example, many signaling pathways involve dimerization, oligomerization, or other types of clustering of signaling proteins as a key step in the signaling cascade. However, it is often experimentally challenging to directly observe and characterize the molecular mechanisms such interactions—the greatest difficulty lies in the fact that living cells have an unknown number of background processes that may or may not participate in the molecular process of interest, and as a consequence, it is usually impossible to definitively correlate an observation to a well-defined cellular mechanism. One of the experimental methods that can quantitatively capture these interactions is through membrane reconstitution, whereby a lipid bilayer is fabricated to mimic the membrane environment, and the biological components of interest are systematically introduced, without unknown background processes. This configuration allows the extensive use of fluorescence techniques, particularly fluorescence fluctuation spectroscopy and single-molecule fluorescence microscopy. In this review, we describe how the equilibrium diffusion of two proteins, K-Ras4B and the PH domain of Bruton’s tyrosine kinase (Btk), on fluid lipid membranes can be used to determine the kinetics of homodimerization reactions.

RAS Nanoclusters: Dynamic Signaling Platforms Amenable to Therapeutic Intervention

RAS proteins are mutated in approximately 20% of all cancers and are generally associated with poor clinical outcomes. RAS proteins are localized to the plasma membrane and function as molecular switches, turned on by partners that receive extracellular mitogenic signals. In the on-state, they activate intracellular signal transduction cascades. Membrane-bound RAS molecules segregate into multimers, known as nanoclusters. These nanoclusters, held together through weak protein-protein and protein-lipid associations, are highly dynamic and respond to cellular input signals and fluctuations in the local lipid environment. Disruption of RAS nanoclusters results in downregulation of RAS-mediated mitogenic signaling. In this review, we discuss the propensity of RAS proteins to display clustering behavior and the interfaces that are associated with these assemblies. Strategies to therapeutically disrupt nanocluster formation or the stabilization of signaling incompetent RAS complexes are discussed.

Inhibition of Nonfunctional Ras

Intuitively, functional states should be targeted; not nonfunctional ones. So why could drugging the inactive K-Ras4BG12Cwork-but drugging the inactive kinase will likely not? The reason is the distinct oncogenic mechanisms. Kinase driver mutations work by stabilizing the active state and/or destabilizing the inactive state. Either way, oncogenic kinases are mostly in the active state. Ras driver mutations work by quelling its deactivation mechanisms, GTP hydrolysis, and nucleotide exchange. Covalent inhibitors that bind to the inactive GDP-bound K-Ras4BG12C conformation can thus work. By contrast, in kinases, allosteric inhibitors work by altering the active-site conformation to favor orthosteric drugs. From the translational standpoint this distinction is vital: it expedites effective pharmaceutical development and extends the drug classification based on the mechanism of action. Collectively, here we postulate that drug action relates to blocking the mechanism of activation, not to whether the protein is in the active or inactive state.

Millisecond Allosteric Dynamics of Activated Ras Reproduced with a Slowly Hydrolyzable GTP Analogue

The millisecond timescale dynamics of activated Ras transiently sample a low-populated conformational state that has distinct surface property from the major state and represents a promising target for binding of small-molecule compounds. To avoid the complications of hydrolysis, dynamics and other properties of active Ras have so far been routinely investigated by using non-hydrolyzable GTP analogues, which, however, were previously reported to alter both the kinetics and distribution of the conformational exchange. In this study, we quantitatively measured and validated the internal dynamics of Ras complexed with a slowly hydrolyzable GTP analogue, GTPγS, which increases the lifetime of active Ras by 23 times relative to that of native GTP. It was found that GTPγS, in addition to its better mimicking of the exchange kinetics than the commonly used non-hydrolyzable analogues GppNHp and GppCH₂ p, can rigorously reproduce the natural dynamics network in active Ras, thus indicating its fitness for use in the development of allosteric inhibitors.

KRAS(G12C)-AMG 510 interaction dynamics revealed by all-atom molecular dynamics simulations

The first KRAS(G12C) targeting inhibitor in clinical development, AMG 510, has shown promising antitumor activity in clinical trials. On the molecular level, however, the interaction dynamics of this covalently bound drug–protein complex has been undetermined. Here, we disclose the interaction dynamics of the KRAS(G12C)–AMG 510 complex by long timescale all-atom molecular dynamics (MD) simulations (total of 75 μs). Moreover, we investigated the influence of the recently reported post-translational modification (PTM) of KRAS’ N-terminus, removal of initiator methionine (iMet1) with acetylation of Thr2, to this complex. Our results demonstrate that AMG 510 does not entrap KRAS into a single conformation, as one would expect based on the crystal structure, but rather into an ensemble of conformations. AMG 510 binding is extremely stable regardless of highly dynamic interface of KRAS’ switches. Overall, KRAS(G12C)–AMG 510 complex partially mimic the native dynamics of GDP bound KRAS; however, AMG 510 stabilizes the α3-helix region. N-terminally modified KRAS displays similar interaction dynamics with AMG 510 as when Met1 is present, but this PTM appears to stabilize β2–β3-loop. These results provide novel conformational insights on the molecular level to KRAS(G12C)–AMG 510 interactions and dynamics, providing new perspectives to RAS related drug discovery.