Probing the Conformational and Energy Landscapes of KRAS Membrane Orientation

KRas4b-Calmodulin Interaction with Membrane Surfaces: Role of Headgroup, Acyl Chain, and Electrostatics

KRas4b is a small plasma membrane-bound G-protein that regulates signal transduction pathways. The interaction of KRas4b with the plasma membrane is governed by both its basic C-terminus, which is farnesylated and methylated, and the lipid composition of the membrane itself. The signaling activity of KRas4b is intricately related to its interaction with various binding partners at the plasma membrane, underlining the critical role played by the lipid environment. The calcium-binding protein calmodulin binds farnesylated KRas4b and plays an important role in the dynamic spatial cycle of KRas4b trafficking in the cell. We utilize Biolayer Interferometry to assay the role of lipid headgroup, chain length, and electrostatics in the dissociation kinetics of fully post-translationally modified KRas4b from Nanodisc bilayers with defined lipid compositions. Our results suggest that calmodulin promotes the dissociation of KRas4b from an anionic membrane, with a comparatively slower displacement of KRas4b from PIP2 relative to PS containing bilayers. In addition to this headgroup dependence, KRas4b dissociation appears to be slower from Nanodiscs wherein the lipid composition contains mismatched, unsaturated acyl chains as compared to lipids with a matched acyl chain length. These findings contribute to understanding the role of the lipid composition in the binding of KRas4b and release from lipid bilayers, showing that the overall charge of the bilayer, the identity of the headgroups present, and the length and saturation of the acyl chains play key roles in KRas4b release from the membrane, potentially providing insights in targeting Ras-membrane interactions for therapeutic interventions.

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).

Intrinsically Disordered Membrane Anchors of Rheb, RhoA, and DiRas3 Small GTPases: Molecular Dynamics, Membrane Organization, and Interactions

Protein structure has been well established to play a key role in determining function; however, intrinsically disordered proteins and regions (IDPs and IDRs) defy this paradigm. IDPs and IDRs exist as an ensemble of structures rather than a stable 3D structure yet play essential roles in many cell-signaling processes. Nearly all Ras superfamily GTPases are tethered to membranes by a lipid tail at the end of a flexible IDR. The sequence of the IDR is a key determinant of membrane localization, and interaction between the IDR and the membrane has been shown to affect signaling in RAS proteins through the modulation of dynamic membrane organization. Here, we utilized atomistic molecular dynamics simulations to study the membrane interaction, conformational dynamics, and lipid sorting of three IDRs from small GTPases Rheb, RhoA, and DiRas3 in model membranes representing their physiological target membranes. We found that complementarity between the lipidated IDR sequence and target membrane lipid composition is a determinant of conformational plasticity. We also show that electrostatic interactions between anionic lipids and basic residues on IDRs are correlated with sampling of semistable conformational substates, and lack of these interactions is associated with greater conformational diversity. Finally, we show that small GTPase IDRs with a polybasic domain alter local lipid composition by segregating anionic lipids and, in some cases, excluding other lipids from their immediate vicinity in favor of anionic lipids.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

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.

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.

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.

Degradation of Rat Sarcoma Proteins Targeting the Post-Translational Prenyl Modifications via Cascade Azidation/Fluorination and Click Reaction

Protein degradation is emerging as a powerful strategy to modulate protein functions and alter cellular signaling pathways. Proteolysis-targeting chimeras (PROTACs) have been used to degrade a range of "undruggable" proteins in cells. Here, we present a type of chemically catalyzed PROTAC to induce rat sarcoma (RAS) degradation based on the chemistry of post-translational prenyl modification. Trimethylsilyl azide and Selectfluor were used to chemically tag the prenyl modification on Caax motif of RAS protein, and a sequential click reaction was applied using the propargyl pomalidomide probe to degrade the prenylated RAS in several cells. Thus, this approach was successfully applied to degrade RAS in multiple cancer cell lines including HeLa, HEK 293T, A549, MCF-7, and HT-29. This novel approach targeting RAS's post-translational prenyl modification to induce RAS degradation by employing the sequential azidation/fluorination and click reaction has been demonstrated efficiently and highly selectively, expanding PROTAC toolsets in the study of disease-relevant protein targets.

The Self-Association of the KRAS4b Protein is Altered by Lipid-Bilayer Composition and Electrostatics

KRAS is a peripheral membrane protein that regulates multiple signaling pathways, and is mutated in ≈30 % of cancers. Transient self-association of KRAS is essential for activation of the downstream effector RAF and oncogenicity. The presence of anionic phosphatidylserine (PS) lipids in the membrane was shown to promote KRAS self-assembly, however, the structural mechanisms remain elusive. Here, we employed nanodisc bilayers of defined lipid compositions, and probed the impact of PS concentration on KRAS self-association. Paramagnetic NMR experiments demonstrated the existence of two transient dimer conformations involving alternate electrostatic contacts between R135 and either D153 or E168 on the "α4/5-α4/5" interface, and revealed that lipid composition and salt modulate their dynamic equilibrium. These dimer interfaces were validated by charge-reversal mutants. This plasticity demonstrates how the dynamic KRAS dimerization interface responds to the environment, and likely extends to the assembly of other signaling complexes on the membrane.

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.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

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.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

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.

Intrinsically disordered proteins and membranes: a marriage of convenience for cell signalling?

The structure-function paradigm has guided investigations into the molecules involved in cellular signalling for decades. The peripheries of this paradigm, however, start to unravel when considering the co-operation between proteins and the membrane in signalling processes. Intrinsically disordered regions hold distinct advantages over folded domains in terms of their binding promiscuity, sensitivity to their particular environment and their ease of modulation through post-translational modifications. Low sequence complexity and bias towards charged residues are also favourable for the multivalent electrostatic interactions that occur at the surfaces of lipid bilayers. This review looks at the principles behind the successful marriage between protein disorder and membranes in addition to the role of this partnership in modifying and regulating signalling in cellular processes. The HVR (hypervariable region) of small GTPases is highlighted as a well-studied example of the nuanced role a short intrinsically disordered region can play in the fine-tuning of signalling pathways.

The dynamic nature of the K-Ras/calmodulin complex can be altered by oncogenic mutations

Oncogenic mutant K-Ras promotes cancer cell proliferation, migration, invasion, and survival by assembling signaling complexes. To date, the functional and structural roles of K-Ras mutations within these complexes are incompletely understood despite their mechanistic and therapeutic significance. Here, we review recent advances in understanding specific binding between K-Ras and the calcium sensor calmodulin. This interaction positively and negatively regulates diverse functions of K-Ras in cancer, suggesting flexibility in K-Ras/calmodulin complex formation. Also, structural data suggest that oncogenic K-Ras likely samples several conformational states, influencing its distinct assemblies with calmodulin and with other proteins. Understanding how K-Ras interacts with calmodulin and with other partners is essential to discovering novel inhibitors of K-Ras in cancer.

RAS Nanoclusters Selectively Sort Distinct Lipid Headgroups and Acyl Chains

RAS proteins are lipid-anchored small GTPases that switch between the GTP-bound active and GDP-bound inactive states. RAS isoforms, including HRAS, NRAS and splice variants KRAS4A and KRAS4B, are some of the most frequently mutated proteins in cancer. In particular, constitutively active mutants of KRAS comprise ∼80% of all RAS oncogenic mutations and are found in 98% of pancreatic, 45% of colorectal and 31% of lung tumors. Plasma membrane (PM) is the primary location of RAS signaling in biology and pathology. Thus, a better understanding of how RAS proteins localize to and distribute on the PM is critical to better comprehend RAS biology and to develop new strategies to treat RAS pathology. In this review, we discuss recent findings on how RAS proteins sort lipids as they undergo macromolecular assembly on the PM. We also discuss how RAS/lipid nanoclusters serve as signaling platforms for the efficient recruitment of effectors and signal transduction, and how perturbing the PM biophysical properties affect the spatial distribution of RAS isoforms and their functions.

Monoubiquitination of KRAS at Lysine104 and Lysine147 Modulates Its Dynamics and Interaction with Partner Proteins

KRAS, a 21 kDa guanine nucleotide-binding protein that functions as a molecular switch, plays a key role in regulating cellular growth. Dysregulation of this key signaling node leads to uncontrolled cell growth, a hallmark of cancer cells. KRAS undergoes post-translational modification by monoubiquitination at various locations, including at lysine104 (K104) and lysine147 (K147). Previous studies have suggested that K104 stabilizes helix-2/helix-3 interactions and K147 is involved in nucleotide binding. However, the impact of monoubiquitination at these residues on the overall structure, dynamics, or function of KRAS is not fully understood. In this study, we examined KRAS monoubiquitination at these sites using data from extensive (12 μs aggregate time) molecular dynamics simulations complemented by nuclear magnetic resonance spectroscopy data. We found that ubiquitin forms dynamic nonspecific interactions with various regions of KRAS and that ubiquitination at both sites modulates conformational fluctuations. In both cases, ubiquitin samples a broad range of conformational space and does not form long-lasting noncovalent contacts with KRAS but it adopts several preferred orientations relative to KRAS. To examine the functional impact of these preferred orientations, we performed a systematic comparison of the dominant configurations of the ubiquitin/KRAS simulated complex with experimental structures of KRAS bound to regulatory and effector proteins as well as a model membrane. Results from these analyses suggest that conformational selection and population shift may minimize the deleterious effects of KRAS ubiquitination at K104 and K147 on binding to some but not all interaction partners. Our findings thus provide new insights into the steric effects of ubiquitin and suggest a potential avenue for therapeutic targeting.

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms.

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.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Influence of Cholesterol on the Orientation of the Farnesylated GTP-Bound KRas-4B Binding with Anionic Model Membranes

The Ras family of proteins is tethered to the inner leaflet of the cell membranes which plays an essential role in signal transduction pathways that promote cellular proliferation, survival, growth, and differentiation. KRas-4B, the most mutated Ras isoform in different cancers, has been under extensive study for more than two decades. Here we have focused our interest on the influence of cholesterol on the orientations that KRas-4B adopts with respect to the plane of the anionic model membranes. How cholesterol in the bilayer might modulate preferences for specific orientation states is far from clear. Herein, after analyzing data from in total 4000 ns-long molecular dynamics (MD) simulations for four KRas-4B systems, properties such as the area per lipid and thickness of the membrane as well as selected radial distribution functions, penetration of different moieties of KRas-4B, and internal conformational fluctuations of flexible moieties in KRas-4B have been calculated. It has been shown that high cholesterol content in the plasma membrane (PM) favors one orientation state (OS1), exposing the effector-binding loop for signal transduction in the cell from the atomic level. We confirm that high cholesterol in the PM helps KRas-4B mutant stay in its constitutively active state, which suggests that high cholesterol intake can increase mortality and may promote cancer progression for cancer patients. We propose that during the treatment of KRas-4B-related cancers, reducing the cholesterol level in the PM and sustaining cancer progression by controlling the plasma cholesterol intake might be taken into account in anti-cancer therapies.

Mechanisms of Ras Membrane Organization and Signaling: Ras Rocks Again

Ras is the most frequently mutated oncogene and recent drug development efforts have spurred significant new research interest. Here we review progress toward understanding how Ras functions in nanoscale, proteo-lipid signaling complexes on the plasma membrane, called nanoclusters. We discuss how G-domain reorientation is plausibly linked to Ras-nanoclustering and -dimerization. We then look at how these mechanistic features could cooperate in the engagement and activation of RAF by Ras. Moreover, we show how this structural information can be integrated with microscopy data that provide nanoscale resolution in cell biological experiments. Synthesizing the available data, we propose to distinguish between two types of Ras nanoclusters, an active, immobile RAF-dependent type and an inactive/neutral membrane anchor-dependent. We conclude that it is possible that Ras reorientation enables dynamic Ras dimerization while the whole Ras/RAF complex transits into an active state. These transient di/oligomer interfaces of Ras may be amenable to pharmacological intervention. We close by highlighting a number of open questions including whether all effectors form active nanoclusters and whether there is an isoform specific composition of Ras nanocluster.

Membrane Tethering Potency of Rab-Family Small GTPases Is Defined by the C-Terminal Hypervariable Regions

Membrane tethering is a crucial step to determine the spatiotemporal specificity of secretory and endocytic trafficking pathways in all eukaryotic endomembrane systems. Recent biochemical studies by a chemically-defined reconstitution approach reveal that, in addition to the structurally-diverse classic tethering factors such as coiled-coil tethering proteins and multisubunit tethering complexes, Rab-family small GTPases also retain the inherent membrane tethering functions to directly and physically bridge two distinct lipid bilayers by themselves. Although Rab-mediated membrane tethering reactions are fairly efficient and specific in the physiological context, its mechanistic basis is yet to be understood. Here, to explore whether and how the intrinsic tethering potency of Rab GTPases is controlled by their C-terminal hypervariable region (HVR) domains that link the conserved small GTPase domains (G-domains) to membrane anchors at the C-terminus, we quantitatively compared tethering activities of two representative Rab isoforms in humans (Rab5a, Rab4a) and their HVR-deleted mutant forms. Strikingly, deletion of the HVR linker domains enabled both Rab5a and Rab4a isoforms to enhance their intrinsic tethering potency, exhibiting 5- to 50-fold higher initial velocities of tethering for the HVR-deleted mutants than those for the full-length, wild-type Rabs. Furthermore, we revealed that the tethering activity of full-length Rab5a was significantly reduced by the omission of anionic lipids and cholesterol from membrane lipids and, however, membrane tethering driven by HVR-deleted Rab5a mutant was completely insensitive to the headgroup composition of lipids. Reconstituted membrane tethering assays with the C-terminally-truncated mutants of Rab4a further uncovered that the N-terminal residues in the HVR linker, located adjacent to the G-domain, are critical for regulating the intrinsic tethering activity. In conclusion, our current findings establish that the non-conserved, flexible C-terminal HVR linker domains define membrane tethering potency of Rab-family small GTPases through controlling the close attachment of the globular G-domains to membrane surfaces, which confers the active tethering-competent state of the G-domains on lipid bilayers.

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as "membrane-distal" and estimate to comprise ∼90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.

How to make an undruggable enzyme druggable: lessons from ras proteins

Significant advances have been made toward discovering allosteric inhibitors for challenging drug targets such as the Ras family of membrane-associated signaling proteins. Malfunction of Ras proteins due to somatic mutations is associated with up to a quarter of all human cancers. Computational techniques have played critical roles in identifying and characterizing allosteric ligand-binding sites on these proteins, and to screen ligand libraries against those sites. These efforts, combined with a wide range of biophysical, structural, biochemical and cell biological experiments, are beginning to yield promising inhibitors to treat malignancies associated with mutated Ras proteins. In this chapter, we discuss some of these developments and how the lessons learned from Ras might be applied to similar other challenging drug targets.

A regulatory role of membrane by direct modulation of the catalytic kinase domain

Cell membrane modulates the function and activity of specific proteins and acts more than just a non-specific scaffolding machinery. In this review, I focus on studies that highlight a direct membrane-mediated modulation of the catalytic kinase domain of a variety of kinases thereby regulating the kinase activity. It emerges that membrane provides a second level of regulation once kinase domain is relieved of its inactive auto-inhibitory state. For the first time a generalized regulatory role of membrane is proposed that governs the kinase activity by modulating the catalytic kinase domain. Striking similarities among a variety of multi-domain kinases as well as single-domain lipidated enzymes such as RAS proteins are presented.

How Anionic Lipids Affect Spatiotemporal Properties of KRAS4B on Model Membranes

RAS proteins are small membrane-anchored GTPases that regulate key cellular signaling networks. It has been recently shown that different anionic lipid types can affect the spatiotemporal properties of RAS through dimerization/clustering and signaling fidelity. To understand the effects of anionic lipids on key spatiotemporal properties of RAS, we dissected 1 ms of data from all-atom molecular dynamics simulations for KRAS4B on two model anionic lipid membranes that have 30% of POPS mixed with neutral POPC and 8% of PIP2 mixed with POPC. We unveiled the orientation space of KRAS4B, whose kinetics were slower and more distinguishable on the membrane containing PIP2 than the membrane containing POPS. Particularly, the PIP2-mixed membrane can differentiate a third kinetic orientation state from the other two known orientation states. We observed that each orientation state may yield different binding modes with an RAF kinase, which is required for activating the MAPK/ERK signaling pathway. However, an overall occluded probability, for which RAF kinases cannot bind KRAS4B, remains unchanged on the two different membranes. We identified rare fast diffusion modes of KRAS4B that appear coupled with orientations exposed to cytosolic RAF. Particularly, on the membrane having PIP2, we found nonlinear correlations between the orientation states and the conformations of the cationic farnesylated hypervariable region, which acts as an anchor in the membrane. Using diffusion coefficients estimated from the all-atom simulations, we quantified the effect of PIP2 and POPS on the KRAS4B dimerization via Green's function reaction dynamics simulations, in which the averaged dimerization rate is 12.5% slower on PIP2-mixed membranes.

The Plasma Membrane as a Competitive Inhibitor and Positive Allosteric Modulator of KRas4B Signaling

Mutant Ras proteins are important drivers of human cancers, yet no approved drugs act directly on this difficult target. Over the last decade, the idea has emerged that oncogenic signaling can be diminished by molecules that drive Ras into orientations in which effector-binding interfaces are occluded by the cell membrane. To support this approach to drug discovery, we characterize the orientational preferences of membrane-bound K-Ras4B in 1.45-ms aggregate time of atomistic molecular dynamics simulations. Individual simulations probe active or inactive states of Ras on membranes with or without anionic lipids. We find that the membrane orientation of Ras is relatively insensitive to its bound guanine nucleotide and activation state but depends strongly on interactions with anionic phosphatidylserine lipids. These lipids slow Ras' translational and orientational diffusion and promote a discrete population in which small changes in orientation control Ras' competence to bind multiple regulator and effector proteins. Our results suggest that compound-directed conversion of constitutively active mutant Ras into functionally inactive forms may be accessible via subtle perturbations of Ras' orientational preferences at the membrane surface.

The current understanding of KRAS protein structure and dynamics

One of the most common drivers in human cancer is the mutant KRAS protein. Not so long ago KRAS was considered as an undruggable oncoprotein. After a long struggle, however, we finally see some light at the end of the tunnel as promising KRAS targeted therapies are in or approaching clinical trials. In recent years, together with the promising progress in RAS drug discovery, our understanding of KRAS has increased tremendously. This progress has been accompanied with a resurgence of publicly available KRAS structures, which were limited to nine structures less than ten years ago. Furthermore, the ever-increasing computational capacity has made biologically relevant timescales accessible, enabling molecular dynamics (MD) simulations to study the dynamics of KRAS protein in more detail at the atomistic level. In this minireview, my aim is to provide the reader an overview of the publicly available KRAS structural data, insights to conformational dynamics revealed by experiments and what we have learned from MD simulations. Also, I will discuss limitations of the current data and provide suggestions for future research related to KRAS, which would fill out the existing gaps in our knowledge and provide guidance in deciphering this enigmatic oncoprotein.