Flexible-body motions of calmodulin and the farnesylated hypervariable region yield a high-affinity interaction enabling K-Ras4B membrane extraction

Two-phase kinetics and cell cortex elastic behavior in Xenopus gastrula cell-cell adhesion

Morphogenetic movements during animal development involve repeated making and breaking of cell-cell contacts. Recent biophysical models of cell-cell adhesion integrate adhesion molecule interactions and cortical cytoskeletal tension modulation, describing equilibrium states for established contacts. We extend this emerging unified concept of adhesion to contact formation kinetics, showing that aggregating Xenopus embryonic cells rapidly achieve Ca2+-independent low-contact states. Subsequent transitions to cadherin-dependent high-contact states show rapid decreases in contact cortical F-actin levels but slow contact area growth. We developed a biophysical model that predicted contact growth quantitatively from known cellular and cytoskeletal parameters, revealing that elastic resistance to deformation and cytoskeletal network turnover are essential determinants of adhesion kinetics. Characteristic time scales of contact growth to low and high states differ by an order of magnitude, being at a few minutes and tens of minutes, respectively, thus providing insight into the timescales of cell-rearrangement-dependent tissue movements.

Cell phenotypes can be predicted from propensities of protein conformations

Proteins exist as dynamic conformational ensembles. Here we suggest that the propensities of the conformations can be predictors of cell function. The conformational states that the molecules preferentially visit can be viewed as phenotypic determinants, and their mutations work by altering the relative propensities, thus the cell phenotype. Our examples include (i) inactive state variants harboring cancer driver mutations that present active state-like conformational features, as in K-Ras4BG12V compared to other K-Ras4BG12X mutations; (ii) mutants of the same protein presenting vastly different phenotypic and clinical profiles: cancer and neurodevelopmental disorders; (iii) alterations in the occupancies of the conformational (sub)states influencing enzyme reactivity. Thus, protein conformational propensities can determine cell fate. They can also suggest the allosteric drugs efficiency.

Higher-order interactions of Bcr-Abl can broaden chronic myeloid leukemia (CML) drug repertoire

Bcr-Abl, a nonreceptor tyrosine kinase, is associated with leukemias, especially chronic myeloid leukemia (CML). Deletion of Abl's N-terminal region, to which myristoyl is linked, renders the Bcr-Abl fusion oncoprotein constitutively active. The substitution of Abl's N-terminal region by Bcr enables Bcr-Abl oligomerization. Oligomerization is critical: it promotes clustering on the membrane, which is essential for potent MAPK signaling and cell proliferation. Here we decipher the Bcr-Abl specific, step-by-step oligomerization process, identify a specific packing surface, determine exactly how the process is structured and identify its key elements. Bcr's coiled coil (CC) domain at the N-terminal controls Bcr-Abl oligomerization. Crystallography validated oligomerization via Bcr-Abl dimerization between two Bcr CC domains, with tetramerization via tight packing between two binary assemblies. However, the structural principles guiding Bcr CC domain oligomerization are unknown, hindering mechanistic understanding and drugs exploiting it. Using molecular dynamics (MD) simulations, we determine that the binary complex of the Bcr CC domain serves as a basic unit in the quaternary complex providing a specific surface for dimer-dimer packing and higher-order oligomerization. We discover that the small α1-helix is the key. In the binary assembly, the helix forms interchain aromatic dimeric packing, and in the quaternary assembly, it contributes to the specific dimer-dimer packing. Our mechanism is supported by the experimental literature. It offers the key elements controlling this process which can expand the drug discovery strategy, including by Bcr CC-derived peptides, and candidate residues for small covalent drugs, toward quenching oligomerization, supplementing competitive and allosteric tyrosine kinase inhibitors.

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

The structural basis of BCR-ABL recruitment of GRB2 in chronic myelogenous leukemia

BCR-ABL drives chronic myeloid leukemia (CML). BCR binding to GRB2 transduces signaling via the Ras/MAPK pathway. Despite considerable data confirming the binding, molecular-level understanding of exactly how the two proteins interact, and especially, what are the determinants of the specificity of the SH2^GRB2 domain-pBCR recognition are still open questions. Yet, this is vastly important for understanding binding selectivity, and for predicting the phosphorylated receptors, or peptides, that are likely to bind. Here, we uncover these determinants and ascertain to what extent they relate to the affinity of the interaction. Toward this end, we modeled the complexes of the phosphorylated BCR (pBCR) and SH2^GRB2 and other pY/Y-peptide-SH2 complexes and compared their specificity and affinity. We observed that pBCR's ₁₇₆FpYVNV₁₈₀ motif is favorable and specific to SH2^GRB2, similar to pEGFR, but not other complexes. SH2^GRB2 contains two binding pockets: pY-binding recognition pocket triggers binding, and the specificity pocket whose interaction is governed by N179 in pBCR and W121 in SH2^GRB2. Our proposed motif with optimal affinity to SH2^GRB2 is E/D-pY-E/V-N-I/L. Collectively, we provide the structural basis of BCR-ABL recruitment of GRB2, outline its specificity hallmarks, and delineate a blueprint for prediction of BCR-binding scaffolds and for therapeutic peptide design.

The mechanism of activation of MEK1 by B-Raf and KSR1

MEK1 interactions with B-Raf and KSR1 are key steps in Ras/Raf/MEK/ERK signaling. Despite this, vital mechanistic details of how these execute signal transduction are still enigmatic. Among these is why, despite B-Raf and KSR1 kinase domains similarity, the B-Raf/MEK1 and KSR1/MEK1 complexes have distinct contributions to MEK1 activation, and broadly, what is KSR1's role. Our molecular dynamics simulations clarify these still unresolved ambiguities. Our results reveal that the proline-rich (P-rich) loop of MEK1 plays a decisive role in MEK1 activation loop (A-loop) phosphorylation. In the inactive B-Raf/MEK1 heterodimer, the collapsed A-loop of B-Raf interacts with the P-rich loop and A-loop of MEK1, minimizing MEK1 A-loop fluctuation and preventing it from phosphorylation. In the active B-Raf/MEK1 heterodimer, the P-rich loop moves in concert with the A-loop of B-Raf as it extends. This reduces the number of residues interacting with MEK1 A-loop, allowing increased A-loop fluctuation, and bringing Ser222 closer to ATP for phosphorylation. B-Raf αG-helix Arg662 promotes MEK1 activation by orienting Ser218 towards ATP. In KSR1/MEK1, the KSR1 αG-helix has Ala826 in place of B-Raf Arg662. This difference results in much fewer interactions between KSR1 αG-helix and MEK1 A-loop, thus a more flexible A-loop. We postulate that if KSR1 were to adopt an active configuration with an extended A-loop as seen in other protein kinases, then the MEK1 P-rich loop would extend in a similar manner, as seen in the active B-Raf/MEK1 heterodimer. This would result in highly flexible MEK1 A-loop, and KSR1 functioning as an active, B-Raf-like, kinase.

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.

The mechanism of activation of monomeric B-Raf V600E

Oncogenic mutations in the serine/threonine kinase B-Raf, particularly the V600E mutation, are frequent in cancer, making it a major drug target. Although much is known about B-Raf's active and inactive states, questions remain about the mechanism by which the protein changes between these two states. Here, we utilize molecular dynamics to investigate both wild-type and V600E B-Raf to gain mechanistic insights into the impact of the Val to Glu mutation. The results show that the wild-type and mutant follow similar activation pathways involving an extension of the activation loop and an inward motion of the αC-helix. The V600E mutation, however, destabilizes the inactive state by disrupting hydrophobic interactions present in the wild-type structure while the active state is stabilized through the formation of a salt bridge between Glu600 and Lys507. Additionally, when the activation loop is extended, the αC-helix is able to move between an inward and outward orientation as long as the DFG motif adopts a specific orientation. In that orientation Phe595 rotates away from the αC-helix, allowing the formation of a salt bridge between Lys483 and Glu501. These mechanistic insights have implications for the development of new Raf inhibitors.

The structural basis of Akt PH domain interaction with calmodulin

Akt plays a key role in the Ras/PI3K/Akt/mTOR signaling pathway. In breast cancer, Akt translocation to the plasma membrane is enabled by the interaction of its pleckstrin homology domain (PHD) with calmodulin (CaM). At the membrane, the conformational change promoted by PIP3 releases CaM and facilitates Thr308 and Ser473 phosphorylation and activation. Here, using modeling and molecular dynamics simulations, we aim to figure out how CaM interacts with Akt's PHD at the atomic level. Our simulations show that CaM-PHD interaction is thermodynamically stable and involves a β-strand rather than an α-helix, in agreement with NMR data, and that electrostatic and hydrophobic interactions are critical. The PHD interacts with CaM lobes; however, multiple modes are possible. IP4, the polar head of PIP3, weakens the CaM-PHD interaction, implicating the release mechanism at the plasma membrane. Recently, we unraveled the mechanism of PI3Kα activation at the atomistic level and the structural basis for Ras role in the activation. Here, our atomistic structural data clarify the mechanism of how CaM interacts, delivers, and releases Akt-the next node in the Ras/PI3K pathway-at the plasma membrane.

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

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.

Mutation-Specific and Common Phosphotyrosine Signatures of KRAS G12D and G13D Alleles

KRAS is one of the most frequently mutated genes across all cancer subtypes. Two of the most frequent oncogenic KRAS mutations observed in patients result in glycine to aspartic acid substitution at either codon 12 (G12D) or 13 (G13D). Although the biochemical differences between these two predominant mutations are not fully understood, distinct clinical features of the resulting tumors suggest involvement of disparate signaling mechanisms. When we compared the global phosphotyrosine proteomic profiles of isogenic colorectal cancer cell lines bearing either G12D or G13D KRAS mutation, we observed both shared as well as unique signaling events induced by the two KRAS mutations. Remarkably, while the G12D mutation led to an increase in membrane proximal and adherens junction signaling, the G13D mutation led to activation of signaling molecules such as nonreceptor tyrosine kinases, MAPK kinases, and regulators of metabolic processes. The importance of one of the cell surface molecules, MPZL1, which was found to be hyperphosphorylated in G12D cells, was confirmed by cellular assays as its knockdown led to a decrease in proliferation of G12D but not G13D expressing cells. Overall, our study reveals important signaling differences across two common KRAS mutations and highlights the utility of our approach to systematically dissect subtle differences between related oncogenic mutants and potentially lead to individualized treatments.

A Non-Canonical Calmodulin Target Motif Comprising a Polybasic Region and Lipidated Terminal Residue Regulates Localization

Calmodulin (CaM) is a Ca²⁺-sensor that regulates a wide variety of target proteins, many of which interact through short basic helical motifs bearing two hydrophobic ‘anchor’ residues. CaM comprises two globular lobes, each containing a pair of EF-hand Ca²⁺-binding motifs that form a Ca²⁺-induced hydrophobic pocket that binds an anchor residue. A central flexible linker allows CaM to accommodate diverse targets. Several reported CaM interactors lack these anchors but contain Lys/Arg-rich polybasic sequences adjacent to a lipidated N- or C-terminus. Ca²⁺-CaM binds the myristoylated N-terminus of CAP23/NAP22 with intimate interactions between the lipid and a surface comprised of the hydrophobic pockets of both lobes, while the basic residues make electrostatic interactions with the negatively charged surface of CaM. Ca²⁺-CaM binds farnesylcysteine, derived from the farnesylated polybasic C-terminus of KRAS4b, with the lipid inserted into the C-terminal lobe hydrophobic pocket. CaM sequestration of the KRAS4b farnesyl moiety disrupts KRAS4b membrane association and downstream signaling. Phosphorylation of basic regions of N-/C-terminal lipidated CaM targets can reduce affinity for both CaM and the membrane. Since both N-terminal myristoylated and C-terminal prenylated proteins use a Singly Lipidated Polybasic Terminus (SLIPT) for CaM binding, we propose these polybasic lipopeptide elements comprise a non-canonical CaM-binding motif.

Calmodulin disrupts plasma membrane localization of farnesylated KRAS4b by sequestering its lipid moiety

KRAS4b is a small guanosine triphosphatase (GTPase) protein that regulates several signal transduction pathways that underlie cell proliferation, differentiation, and survival. KRAS4b function requires prenylation of its C terminus and recruitment to the plasma membrane, where KRAS4b activates effector proteins including the RAF family of kinases. The Ca²⁺-sensing protein calmodulin (CaM) has been suggested to regulate the localization of KRAS4b through direct, Ca²⁺-dependent interaction, but how CaM and KRAS4b functionally interact is controversial. Here, we determined a crystal structure, which was supported by solution nuclear magnetic resonance (NMR), that revealed the sequestration of the prenyl moiety of KRAS4b in the hydrophobic pocket of the C-terminal lobe of Ca²⁺-bound CaM. Our engineered fluorescence resonance energy transfer (FRET)-based biosensor probes (CaMeRAS) showed that, upon stimulation of Ca²⁺ influx by extracellular ligands, KRAS4b reversibly translocated in a Ca²⁺-CaM-dependent manner from the plasma membrane to the cytoplasm in live HeLa and HEK293 cells. These results reveal a mechanism underlying the inhibition of KRAS4b activity by Ca²⁺ signaling pathways.

The quaternary assembly of KRas4B with Raf-1 at the membrane

Proximally located in the membrane, oncogenic Ras dimers (or nanoclusters) can recruit and promote Raf dimerization and MAPK (Raf/MEK/ERK) signaling. Among Ras isoforms, KRas4B is the most frequently mutated. Recent data on the binary KRas4B–Raf-1 complex suggested that Raf-1 CRD not only executes membrane anchorage, but also supports the high-affinity interaction of Raf-1 RBD with KRas4B catalytic domain. For a detailed mechanistic picture of Raf activation at the membrane, we employ explicit MD simulations of the quaternary KRas4B–Raf-1 complex. The complex contains two active GTP-bound KRas4B proteins forming a dimer through the allosteric lobe interface and two tandem RBD-CRD segments of Raf-1 interacting with the effector lobes at both ends of the KRas4B dimer. We show that Raf-1 RBD-CRD supports stable KRas4B dimer at preferred interface and orientation at the membrane, thereby cooperatively enhancing the affinity of the KRas4B–Raf-1 interaction. We propose that a Ras dimer at the membrane can increase the population of proximal Raf kinase domains, promoting kinase domain dimerization in the cytoplasm. Collectively, the dynamic Ras–Raf assembly promotes Raf activation not by allostery; instead, Ras activates Raf by shifting its ensemble toward kinase domain-accessible states through enhanced affinity at the membrane.

Ca2+-Dependent Switch of Calmodulin Interaction Mode with Tandem IQ Motifs in the Scaffolding Protein IQGAP1

IQ domain GTPase-activating scaffolding protein 1 (IQGAP1) mediates cytoskeleton, cell migration, proliferation, and apoptosis events. Calmodulin (CaM) modulates IQGAP1 functions by binding to its four tandem IQ motifs. Exactly how CaM binds the IQ motifs and which functions of IQGAP1 CaM regulates and how are fundamental mechanistic questions. We combine experimental pull-down assays, mutational data, and molecular dynamics simulations to understand the IQ-CaM complexes with and without Ca2+ at the atomic level. Apo-CaM favors the IQ3 and IQ4 motifs but not the IQ1 and IQ2 motifs that lack two hydrophobic residues for interactions with apo-CaM's hydrophobic pocket. Ca2+-CaM binds all four IQ motifs, with both N- and C-lobes tightly wrapped around each motif. Ca2+ promotes IQ-CaM interactions and increases the amount of IQGAP1-loaded CaM for IQGAP1-mediated signaling. Collectively, we describe IQ-CaM binding in atomistic detail and feature the emergence of Ca2+ as a key modulator of the CaM-IQGAP1 interactions.

KRAS: A Promising Therapeutic Target for Cancer Treatment

Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is the most commonly mutated oncogene in human cancer. The developments of many cancers depend on sustained expression and signaling of KRAS, which makes KRAS a high-priority therapeutic target. Scientists have not successfully developed drugs that target KRAS, although efforts have been made last three decades. In this review, we highlight the emerging experimental strategies of impairing KRAS membrane localization and the direct targeting of KRAS. We also conclude the combinatorial therapies and RNA interference technology for the treatment of KRAS mutant cancers. Moreover, the virtual screening approach to discover novel KRAS inhibitors and synthetic lethality interactors of KRAS are discussed in detail.

Does Ras Activate Raf and PI3K Allosterically?

The mechanism through which oncogenic Ras activates its effectors is vastly important to resolve. If allostery is at play, then targeting allosteric pathways could help in quelling activation of MAPK (Raf/MEK/ERK) and PI3K (PI3K/Akt/mTOR) cell proliferation pathways. On the face of it, allosteric activation is reasonable: Ras binding perturbs the conformational ensembles of its effectors. Here, however, we suggest that at least for Raf, PI3K, and NORE1A (RASSF5), that is unlikely. Raf's long disordered linker dampens effective allosteric activation. Instead, we suggest that the high-affinity Ras–Raf binding relieves Raf's autoinhibition, shifting Raf's ensemble from the inactive to the nanocluster-mediated dimerized active state, as Ras also does for NORE1A. PI3K is recruited and allosterically activated by RTK (e.g., EGFR) at the membrane. Ras restrains PI3K's distribution and active site orientation. It stabilizes and facilitates PIP₂ binding at the active site and increases the PI3K residence time at the membrane. Thus, RTKs allosterically activate PI3Kα; however, merging their action with Ras accomplishes full activation. Here we review their activation mechanisms in this light and draw attention to implications for their pharmacology.

The Structural Basis of the Farnesylated and Methylated KRas4B Interaction with Calmodulin

Ca²⁺-calmodulin (CaM) extracts KRas4B from the plasma membrane, suggesting that KRas4B/CaM interaction plays a role in regulating Ras signaling. To gain mechanistic insight, we provide a computational model, supported by experimental structural data, of farnesylated/methylated KRas4B_1-185 interacting with CaM in solution and at anionic membranes including signaling lipids. Due to multiple interaction modes, we observe diverse conformational ensembles of the KRas4B-CaM complex. A highly populated conformation reveals the catalytic domain interacting with the N-lobe and the hypervariable region (HVR) wrapping around the linker with the farnesyl docking to the extended CaM's C-lobe pocket. Alternatively, KRas4B can interact with collapsed CaM with the farnesyl penetrating CaM's center. At anionic membranes, CaM interacts with the catalytic domain with large fluctuations, drawing the HVR. Signaling lipids establishing strong salt bridges with CaM prevent membrane departure. Membrane-interacting KRas4B-CaM complex can productively recruit phosphatidylinositol 3-kinase α (PI3Kα) to the plasma membrane, serving as a coagent in activating PI3Kα/Akt signaling.

The structural basis for Ras activation of PI3Kα lipid kinase

PI3Kα is a principal Ras effector that phosphorylates PIP2 to PIP3 in the PI3K/Akt/mTOR pathway. How Ras activates PI3K has been unclear: is Ras' role confined to PI3K recruitment to the membrane or does Ras activation also involve allostery? Recently, we determined the mechanism of PI3Kα activation at the atomic level. We showed the vital role and significance of conformational change in PI3Kα activation. Here, by a 'best-match for hydrogen-bonding pair' (BMHP) computational protocol and molecular dynamics (MD) simulations, we model the atomic structure of KRas4B in complex with the Ras binding domain (RBD) of PI3Kα, striving to understand the mechanism of PI3Kα activation by Ras. Point mutations T208D, K210E, and K227E disrupt the KRas4B-RBD interface in the models, in line with the experiments. We identify allosteric signaling pathways connecting Ras to RBD in the p110α subunit. However, the observed weak allosteric signals coupled with the detailed mechanism of PI3Kα activation make us conclude that the dominant mechanistic role of Ras is likely to be recruitment and restriction of the PI3Kα population at the membrane. Thus, RTK recruits the PI3Kα to the membrane and activates it by relieving its autoinhibition exerted by the nSH2 domain, leading to exposure of the kinase domain, which permits PIP2 binding. Ras recruitment can shift the PI3Kα ensemble toward a population where the kinase domain surface and the active site position and orientation favor PIP2 insertion. This work helps elucidate Ras-mediated PI3K activation and explores the structural basis for Ras-PI3Kα drug discovery.

KRAS Prenylation Is Required for Bivalent Binding with Calmodulin in a Nucleotide-Independent Manner

Deregulation of KRAS4b signaling pathway has been implicated in 30% of all cancers. Membrane localization of KRAS4b is an essential step for the initiation of the downstream signaling cascades that guide various cellular mechanisms. KRAS4b plasma membrane (PM) binding is mediated by the insertion of a prenylated moiety that is attached to the terminal carboxy-methylated cysteine, in addition to electrostatic interactions of its positively charged hypervariable region with anionic lipids. Calmodulin (CaM) has been suggested to selectively bind KRAS4b to act as a negative regulator of the RAS/mitogen-activated protein kinase (MAPK) signaling pathway by displacing KRAS4b from the membrane. However, the mechanism by which CaM can recognize and displace KRAS4b from the membrane is not well understood. In this study, we employed biophysical and structural techniques to characterize this mechanism in detail. We show that KRAS4b prenylation is required for binding to CaM and that the hydrophobic pockets of CaM can accommodate the prenylated region of KRAS4b, which might represent a novel CaM-binding motif. Remarkably, prenylated KRAS4b forms a 2:1 stoichiometric complex with CaM in a nucleotide-independent manner. The interaction between prenylated KRAS4b and CaM is enthalpically driven, and electrostatic interactions also contribute to the formation of the complex. The prenylated KRAS4b terminal KSKTKC-farnesylation and carboxy-methylation is sufficient for binding and defines the minimal CaM-binding motif. This is the same region implicated in membrane and phosphodiesterase6-δ binding. Finally, we provide a structure-based docking model by which CaM binds to prenylated KRAS4b. Our data provide new insights into the KRAS4b-CaM interaction and suggest a possible mechanism whereby CaM can regulate KRAS4b membrane localization.

Calmodulin (CaM) Activates PI3Kα by Targeting the "Soft" CaM-Binding Motifs in Both the nSH2 and cSH2 Domains of p85α

PI3Kα is a key lipid kinase in the PI3K/Akt pathway. Its frequent oncogenic mutations make it a primary drug target. Calmodulin (CaM) activates PI3Kα independently of extracellular signals, indicating a significant role in oncogenic PI3Kα activation. Here, we reveal the atomic-scale structures of CaM in complexes with the nSH2 and cSH2 domains of the regulatory p85α subunit of PI3Kα, and illustrate how CaM activates PI3Kα by targeting the "soft 1-5-10" CaM-binding motifs in both nSH2 and cSH2 domains. Experiment observed CaM binding cSH2 first, followed by nSH2 binding hours later. CaM typically prefers binding helical peptides. Here we observe that, unlike in cSH2, the CaM-binding motif in nSH2 populates a mixed β-sheet/α-helix/random coil structure. The population shift from a β-sheet toward CaM's favored α-helical conformation explains why the nSH2 domain needs a longer time for CaM binding in the experiments. The "soft" CaM-binding motifs in both nSH2 and cSH2 domains establish strong CaM-PI3Kα interactions, collectively facilitating PI3Kα activation. This work uncovers the structural basis for CaM-driven PI3Kα activation.

Mutations in LZTR1 drive human disease by dysregulating RAS ubiquitination

The leucine zipper-like transcriptional regulator 1 (LZTR1) protein, an adaptor for cullin 3 (CUL3) ubiquitin ligase complex, is implicated in human disease, yet its mechanism of action remains unknown. We found that Lztr1 haploinsufficiency in mice recapitulates Noonan syndrome phenotypes, whereas LZTR1 loss in Schwann cells drives dedifferentiation and proliferation. By trapping LZTR1 complexes from intact mammalian cells, we identified the guanosine triphosphatase RAS as a substrate for the LZTR1-CUL3 complex. Ubiquitome analysis showed that loss of Lztr1 abrogated Ras ubiquitination at lysine-170. LZTR1-mediated ubiquitination inhibited RAS signaling by attenuating its association with the membrane. Disease-associated LZTR1 mutations disrupted either LZTR1-CUL3 complex formation or its interaction with RAS proteins. RAS regulation by LZTR1-mediated ubiquitination provides an explanation for the role of LZTR1 in human disease.

Modeling and subtleties of K-Ras and Calmodulin interaction

K-Ras, one of the most common small GTPases of the cell, still presents many riddles, despite the intense efforts to unveil its mysteries. Such is the case of its interaction with Calmodulin, a small acidic protein known for its role as a calcium ion sensor. Although the interaction between these two proteins and its biological implications have been widely studied, a model of their interaction has not been performed. In the present work we analyse this intriguing interaction by computational means. To do so, both conventional molecular dynamics and scaled molecular dynamics have been used. Our simulations suggest a model in which Calmodulin would interact with both the hypervariable region and the globular domain of K-Ras, using a lobe to interact with each of them. According to the presented model, the interface of helixes α4 and α5 of the globular domain of K-Ras would be relevant for the interaction with a lobe of Calmodulin. These results were also obtained when bringing the proteins together in a step wise manner with the umbrella sampling methodology. The computational results have been validated using SPR to determine the relevance of certain residues. Our results demonstrate that, when mutating residues of the α4-α5 interface described to be relevant for the interaction with Calmodulin, the interaction of the globular domain of K-Ras with Calmodulin diminishes. However, it is to be considered that our simulations indicate that the bulk of the interaction would fall on the hypervariable region of K-Ras, as many more interactions are identified in said region. All in all our simulations present a suitable model in which K-Ras could interact with Calmodulin at membrane level using both its globular domain and its hypervariable region to stablish an interaction that leads to an altered signalling.

K-Ras is one of the most mutated oncogenes in human cancer. Although several studies validate K-Ras protein as good candidate for direct therapeutic targeting, pharmacologic targeting has not been successful. During the last years increasing evidences demonstrate that oncogenic K-Ras activity can be modulated in vivo by dimerization, nanoclustering at the plasma membrane or interaction with non-effector proteins, consequently opening new therapeutic strategies. We have previously demonstrated that Calmodulin, an ubiquitous Ca2+-binding protein, is one of this K-Ras interacting proteins and that it negatively modulates K-Ras signaling. Although experimental data were available showing the relevant regions for this interaction, a model of K-Ras and Calmodulin interaction was missing. In the present work by using different computational modeling techniques we obtained a model for this interaction that agrees with the experimental data. We believe the present model will help to better understand K-Ras regulation, and to design new inhibitors. For instance, base on our model, we can predict that the interaction can take place at the plasma membrane, and that since the surface of K-Ras that interact with Calmodulin is the same that it uses for dimerization, that Calmodulin could be inhibiting K-Ras dimerization.

Autoinhibition in Ras effectors Raf, PI3Kα, and RASSF5: a comprehensive review underscoring the challenges in pharmacological intervention

Autoinhibition is an effective mechanism that guards proteins against spurious activation. Despite its ubiquity, the distinct organizations of the autoinhibited states and their release mechanisms differ. Signaling is most responsive to the cell environment only if a small shift in the equilibrium is required to switch the system from an inactive (occluded) to an active (exposed) state. Ras signaling follows this paradigm. This underscores the challenge in pharmacological intervention to exploit and enhance autoinhibited states. Here, we review autoinhibition and release mechanisms at the membrane focusing on three representative Ras effectors, Raf protein kinase, PI3Kα lipid kinase, and NORE1A (RASSF5) tumor suppressor, and point to the ramifications to drug discovery. We further touch on Ras upstream and downstream signaling, Ras activation, and the Ras superfamily in this light, altogether providing a broad outlook of the principles and complexities of autoinhibition.

Allosteric KRas4B Can Modulate SOS1 Fast and Slow Ras Activation Cycles

Membrane-anchored Ras family proteins are activated by guanine nucleotide exchange factors such as SOS1. The CDC25 domain of SOS1 catalyzes GDP-to-GTP exchange, thereby activating Ras. Here, we aim to decipher the activation mechanism of KRas4B, a significantly mutated oncogene. We perform large-scale molecular dynamics simulations on 12 SOS1 systems, scrutinizing each step in two possible KRas4B activation cycles, fast and slow. To activate KRas4B at the CDC25 catalytic site, the allosteric site in the Ras exchanger motif (REM) domain of SOS1 needs to recruit a (nucleotide-bound) KRas4B molecule. Our simulations indicate that KRas4B-GTP interacts with the REM allosteric site more strongly than with the CDC25 catalytic site, consistent with its allosteric role in the GDP-to-GTP exchange. In the fast cycle, the allosteric KRas4B-GTP induces conformational change at the catalytic site. The conformational change facilitates loading KRas4B-GDP at the catalytic site and opening the KRas4B nucleotide-binding site for GDP release and GTP binding. GTP binding reduces the affinity of KRas4B-GTP to the CDC25 catalytic site, resulting in its release. By contrast, in the slow cycle, KRas4B-GDP binds at the allosteric REM site. The limited, altered conformational change that it induces prevents the exact alignments of switch I and II of KRas4B. The increasing binding strength at both binding sites due to interactions of regions other than switch I and II retards GDP release from the catalytic KRas4B, thus KRas4B activation. The accelerated activation cycle supports a positive feedback loop with allosteric signals communicating between the two Ras molecules and is the predominant, native function of SOS. SOS1 activation details may help drug discovery to inhibit Ras activation.

Interaction of Calmodulin with the cSH2 Domain of the p85 Regulatory Subunit

Calmodulin (CaM) is a calcium sensor protein that directly interacts with the dual-specificity (lipid and protein) kinase PI3Kα through the SH2 domains of the p85 regulatory subunit. In adenocarcinomas, the CaM interaction removes the autoinhibition of the p110 catalytic subunit of PI3Kα, leading to activation of PI3Kα and promoting cell proliferation, survival, and migration. Here we demonstrate that the cSH2 domain of p85α engages its two CaM-binding motifs in the interaction with the N- and C-lobes of CaM as well as the flexible central linker, and our nuclear magnetic resonance experiments provide structural details. We show that in response to binding CaM, cSH2 exposes its tryptophan residue at the N-terminal region to the solvent. Because of the flexible nature of both CaM and cSH2, multiple binding modes of the interactions are possible. Binding of CaM to the cSH2 domain can help release the inhibition imposed on the p110 subunit, similar to the binding of the phosphorylated motif of RTK, or phosphorylated CaM (pCaM), to the SH2 domains. Amino acid sequence analysis shows that CaM-binding motifs are common in SH2 domains of non-RTKs. We speculate that CaM can also activate these kinases through similar mechanisms.

Raf-1 Cysteine-Rich Domain Increases the Affinity of K-Ras/Raf at the Membrane, Promoting MAPK Signaling

K-Ras4B preferentially activates Raf-1. The high-affinity interaction of Ras-binding domain (RBD) of Raf with Ras was solved, but the relative position of Raf's cysteine-rich domain (CRD) in the Ras/Raf complex at the membrane and key question of exactly how it affects Raf signaling are daunting. We show that CRD stably binds anionic membranes inserting a positively charged loop into the amphipathic interface. Importantly, when in complex with Ras/RBD, covalently connected CRD presents the same membrane interaction mechanism, with CRD locating at the space between the RBD and membrane. To date, CRD's role was viewed in terms of stabilizing Raf-membrane interaction. Our observations argue for a key role in reducing Ras/RBD fluctuations at the membrane, thereby increasing Ras/RBD affinity. Even without K-Ras, via CRD, Raf-1 can recruit to the membrane; however, by reducing the Ras/RBD fluctuations and enhancing Ras/RBD affinity at the membrane, CRD promotes Raf's activation and MAPK signaling over other pathways.

Oncogenic Ras Isoforms Signaling Specificity at the Membrane

How do Ras isoforms attain oncogenic specificity at the membrane? Oncogenic KRas, HRas, and NRas (K-Ras, H-Ras, and N-Ras) differentially populate distinct cancers. How they selectively activate effectors and why is KRas4B the most prevalent are highly significant questions. Here, we consider determinants that may bias isoform-specific effector activation and signaling at the membrane. We merge functional data with a conformational view to provide mechanistic insight. Cell-specific expression levels, pathway cross-talk, and distinct interactions are the key, but conformational trends can modulate selectivity. There are two major pathways in oncogenic Ras-driven proliferation: MAPK (Raf/MEK/ERK) and PI3Kα/Akt/mTOR. All membrane-anchored, proximally located, oncogenic Ras isoforms can promote Raf dimerization and fully activate MAPK signaling. So why the differential statistics of oncogenic isoforms in distinct cancers and what makes KRas so highly oncogenic? Many cell-specific factors may be at play, including higher KRAS mRNA levels. As a key factor, we suggest that because only KRas4B binds calmodulin, only KRas can fully activate PI3Kα/Akt signaling. We propose that full activation of both MAPK and PI3Kα/Akt proliferative pathways by oncogenic KRas4B-but not by HRas or NRas-may help explain why the KRas4B isoform is especially highly populated in certain cancers. We further discuss pharmacologic implications. Cancer Res; 78(3); 593-602. ©2017 AACR.

Oncogenic Ras Isoforms Signaling Specificity at the Membrane

How do Ras isoforms attain oncogenic specificity at the membrane? Oncogenic KRas, HRas, and NRas (K-Ras, H-Ras, and N-Ras) differentially populate distinct cancers. How they selectively activate effectors and why is KRas4B the most prevalent are highly significant questions. Here, we consider determinants that may bias isoform-specific effector activation and signaling at the membrane. We merge functional data with a conformational view to provide mechanistic insight. Cell-specific expression levels, pathway cross-talk, and distinct interactions are the key, but conformational trends can modulate selectivity. There are two major pathways in oncogenic Ras-driven proliferation: MAPK (Raf/MEK/ERK) and PI3Kα/Akt/mTOR. All membrane-anchored, proximally located, oncogenic Ras isoforms can promote Raf dimerization and fully activate MAPK signaling. So why the differential statistics of oncogenic isoforms in distinct cancers and what makes KRas so highly oncogenic? Many cell-specific factors may be at play, including higher KRAS mRNA levels. As a key factor, we suggest that because only KRas4B binds calmodulin, only KRas can fully activate PI3Kα/Akt signaling. We propose that full activation of both MAPK and PI3Kα/Akt proliferative pathways by oncogenic KRas4B-but not by HRas or NRas-may help explain why the KRas4B isoform is especially highly populated in certain cancers. We further discuss pharmacologic implications. Cancer Res; 78(3); 593-602. ©2017 AACR.

Is Nanoclustering essential for all oncogenic KRas pathways? Can it explain why wild-type KRas can inhibit its oncogenic variant?

Membrane-anchored oncogenic KRas can dimerize, form nanoclusters, and signal through the MAPK (Raf/MEK/ERK) and PI3Kα/Akt/mTOR. Both pathways are needed in KRAS-driven proliferation. Here we ask: Is oncogenic KRas nanoclustering (or dimerization) essential for all KRas signaling pathways? Raf kinase domain dimerization, thus MAPK activation, requires KRas nanoclusters. By contrast, the PI3Kα heterodimer acts as a monomeric unit; thus, does PI3Kα activation and PI3Kα/Akt/mTOR signaling require nanoclustering? Further, calmodulin binds only to oncogenic KRas4B. Here we ask: Does calmodulin downregulate KRas4B cancer development as suggested early on, or promote it? We also ask: Why is oncogenic KRas4B the most abundant isoform? Does wild-type Ras indeed inhibit its oncogenic variants as data appeared to suggest? And related to the last question, why is wild-type KRas a more potent inhibitor of its oncogenic form than wild-type NRas of its oncogenic form? Resolving these cardinal questions, and others, such as how exactly does RASSF5 (NORE1A) act as tumor suppressor, and why Ras isoforms tend to occur in distinct cancer types are crucial for effective pharmacology. In this review, we take a nanoclustering/dimerization-centric outlook and show that many questions can be explained by simply considering Ras nanoclustering.

Computational Insights into the Interactions between Calmodulin and the c/nSH2 Domains of p85α Regulatory Subunit of PI3Kα: Implication for PI3Kα Activation by Calmodulin

Calmodulin (CaM) and phosphatidylinositide-3 kinase (PI3Kα) are well known for their multiple roles in a series of intracellular signaling pathways and in the progression of several human cancers. Crosstalk between CaM and PI3Kα has been an area of intensive research. Recent experiments have shown that in adenocarcinoma, K-Ras4B is involved in the CaM-PI3Kα crosstalk. Based on experimental results, we have recently put forward a hypothesis that the coordination of CaM and PI3Kα with K-Ras4B forms a CaM-PI3Kα-K-Ras4B ternary complex, which leads to the formation of pancreatic ductal adenocarcinoma. However, the mechanism for the CaM-PI3Kα crosstalk is unresolved. Based on molecular modeling and molecular dynamics simulations, here we explored the potential interactions between CaM and the c/nSH2 domains of p85α subunit of PI3Kα. We demonstrated that CaM can interact with the c/nSH2 domains and the interaction details were unraveled. Moreover, the possible modes for the CaM-cSH2 and CaM-nSH2 interactions were uncovered and we used them to construct a complete CaM-PI3Kα complex model. The structural model of CaM-PI3Kα interaction not only offers a support for our previous ternary complex hypothesis, but also is useful for drug design targeted at CaM-PI3Kα protein-protein interactions.

Calmodulin fishing with a structurally disordered bait triggers CyaA catalysis

Once translocated into the cytosol of target cells, the catalytic domain (AC) of the adenylate cyclase toxin (CyaA), a major virulence factor of Bordetella pertussis, is potently activated by binding calmodulin (CaM) to produce supraphysiological levels of cAMP, inducing cell death. Using a combination of small-angle X-ray scattering (SAXS), hydrogen/deuterium exchange mass spectrometry (HDX-MS), and synchrotron radiation circular dichroism (SR-CD), we show that, in the absence of CaM, AC exhibits significant structural disorder, and a 75-residue-long stretch within AC undergoes a disorder-to-order transition upon CaM binding. Beyond this local folding, CaM binding induces long-range allosteric effects that stabilize the distant catalytic site, whilst preserving catalytic loop flexibility. We propose that the high enzymatic activity of AC is due to a tight balance between the CaM-induced decrease of structural flexibility around the catalytic site and the preservation of catalytic loop flexibility, allowing for fast substrate binding and product release. The CaM-induced dampening of AC conformational disorder is likely relevant to other CaM-activated enzymes.

Calmodulin is a widespread and highly conserved protein that interacts with a wide variety of eukaryotic proteins and enzymes, controlling their activities in response to calcium. The adenylate cyclase toxin (CyaA) of Bordetella pertussis, the causative agent of whooping cough, is one such calmodulin target. Once transported across the plasma membrane of eukaryotic cells, the catalytic domain (AC) of CyaA is activated by calmodulin, producing high levels of cAMP, which can induce cell death. We use an integrative structural biology approach combining several biophysical techniques to characterize the structural rearrangements in AC upon calmodulin binding and to elucidate their relationship to CyaA activation. We show that a disordered stretch of 75 amino acid residues in AC serves as a bait for calmodulin capture. Binding induces significant folding within this region, a prerequisite for CyaA activation. Calmodulin binding promotes the stabilization of the distant catalytic site, whilst maintaining its catalytic loop in a flexible and exposed state. Both phenomena contribute to the high enzymatic activity of AC, allowing for fast substrate binding and cAMP release. The calmodulin-induced reduction of AC conformational disorder is likely relevant to other calmodulin-activated enzymes.

Phosphorylated Calmodulin Promotes PI3K Activation by Binding to the SH2 Domains

How calmodulin (CaM) acts in KRAS-driven cancers is a vastly important question. CaM binds to and stimulates PI3Kα/Akt signaling, promoting cell growth and proliferation. Phosphorylation of CaM at Tyr99 (pY99) enhances PI3Kα activation. PI3Kα is a lipid kinase. It phosphorylates PIP2 to produce PIP3, to which Akt binds. PI3Kα has two subunits: the regulatory p85 and the catalytic p110. Here, exploiting explicit-solvent MD simulations we unveil key interactions between phosphorylated CaM (pCaM) and the two SH2 domains in the p85 subunit, confirm experimental observations, and uncover PI3Kα's mechanism of activation. pCaMs form strong and stable interactions with both nSH2 and cSH2 domains, with pY99 being the dominant contributor. Despite the high structural similarity between the two SH2 domains, we observe that nSH2 prefers an extended CaM conformation, whereas cSH2 prefers a collapsed conformation. Notably, collapsed CaM is observed after binding of an extended CaM to K-Ras4B. Thus, the more populated extended pCaM conformation targets nSH2 to release its autoinhibition of p110 catalytic sites. This executes the key activation step of PI3Kα. Independently, K-Ras4B allosterically activates p110. These events are at the cell membrane, which contributes to tighten the PI3Kα Ras binding domain/K-Ras4B interaction, to accomplish K-Ras4B allosteric activation, with a minor contribution from cSH2.