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

Dynamic Coupling and Entropy Changes in KRAS G12D Mutation: Insights into Molecular Flexibility, Allostery and Function

The oncogenic G12D mutation in KRAS is a major driver of cancer progression, yet the complete mechanism by which this mutation alters protein dynamics and function remains incompletely understood. Here, we investigate how the G12D mutation alters KRAS's conformational landscape and residue-residue interactions using molecular dynamics simulations coupled with entropy calculations and mutual information (MI) analysis. We demonstrate that the mutation increases local entropy at key functional residues (D12, Y32, G60, and Q61), and introduces new peaks to the Ramachandran angles, disrupting the precise structural alignment necessary for GTP hydrolysis. Notably, while individual residue entropy increases, joint entropy analysis shows a complex reorganization pattern. MI analysis identifies enhanced dynamic coupling between distant residues, suggesting that the mutation establishes new long-range interactions that stabilize the active state. These findings show how G12D mutation redefines KRAS's dynamic network, leading to persistent activation through enhanced residue coupling rather than mere local disruption. Our results suggest novel therapeutic strategies focused on modulating protein dynamics rather than targeting specific binding sites, potentially offering new approaches to combat KRAS-driven cancers.

Existence of Singularities in NMR Relaxation Dispersion Profiles: Implications for Hidden Dynamics

It is common for NMR relaxation dispersion experiments to suggest an absence of dynamics despite anecdotal indications of conformational dynamics. We explore the potential explanations and approaches to this conundrum. Some inconsistencies have been observed between two relaxation dispersion experiments, Carr-Purcell-Meiboom-Gill (CPMG) and adiabatic relaxation dispersion experiments, in recent dynamic studies of different biomolecules. Theoretical analyses show that such seemingly paradoxical results might come from a complex exchange topology that is concealed by the application of the simple two-site exchange model for interpretation of the relaxation dispersion data. Scenarios are explored and revealed in which the presence of complex millisecond conformational exchange could suppress the amplitude of CPMG relaxation dispersion profiles, even when the exchange rates are within the detectable range of the experiments. With experimental errors, the suppressed relaxation dispersion profiles could lead to the conclusion of "no millisecond conformational exchange". However, such hidden dynamics can potentially be detected by adiabatic relaxation dispersion experiments. Finally, we demonstrate some advantages of adiabatic relaxation dispersion experiments over conventional relaxation dispersion experiments and a simplified computational approach to analyze the adiabatic relaxation dispersion profiles.

Producing recombinant proteins in Vibrio natriegens

The diversity of chemical and structural attributes of proteins makes it inherently difficult to produce a wide range of proteins in a single recombinant protein production system. The nature of the target proteins themselves, along with cost, ease of use, and speed, are typically cited as major factors to consider in production. Despite a wide variety of alternative expression systems, most recombinant proteins for research and therapeutics are produced in a limited number of systems: Escherichia coli, yeast, insect cells, and the mammalian cell lines HEK293 and CHO. Recent interest in Vibrio natriegens as a new bacterial recombinant protein expression host is due in part to its short doubling time of ≤ 10 min but also stems from the promise of compatibility with techniques and genetic systems developed for E. coli. We successfully incorporated V. natriegens as an additional bacterial expression system for recombinant protein production and report improvements to published protocols as well as new protocols that expand the versatility of the system. While not all proteins benefit from production in V. natriegens, we successfully produced several proteins that were difficult or impossible to produce in E. coli. We also show that in some cases, the increased yield is due to higher levels of properly folded protein. Additionally, we were able to adapt our enhanced isotope incorporation methods for use with V. natriegens. Taken together, these observations and improvements allowed production of proteins for structural biology, biochemistry, assay development, and structure-based drug design in V. natriegens that were impossible and/or unaffordable to produce in E. coli.

Molecular Mechanism of Phosphorylation-Mediated Impacts on the Conformation Dynamics of GTP-Bound KRAS Probed by GaMD Trajectory-Based Deep Learning

The phosphorylation of different sites produces a significant effect on the conformational dynamics of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations were combined with deep learning (DL) to explore the molecular mechanism of the phosphorylation-mediated effect on conformational dynamics of the GTP-bound KRAS. The DL finds that the switch domains are involved in obvious differences in conformation contacts and suggests that the switch domains play a key role in the function of KRAS. The analyses of free energy landscapes (FELs) reveal that the phosphorylation of pY32, pY64, and pY137 leads to more disordered states of the switch domains than the wild-type (WT) KRAS and induces conformational transformations between the closed and open states. The results from principal component analysis (PCA) indicate that principal motions PC1 and PC2 are responsible for the closed and open states of the phosphorylated KRAS. Interaction networks were analyzed and the results verify that the phosphorylation alters interactions of GTP and magnesium ion Mg2+ with the switch domains. It is concluded that the phosphorylation pY32, pY64, and pY137 tune the activity of KRAS through changing conformational dynamics and interactions of the switch domains. We anticipated that this work could provide theoretical aids for deeply understanding the function of KRAS.

Comparative analysis of KRAS4a and KRAS4b splice variants reveals distinctive structural and functional properties

KRAS, the most frequently mutated oncogene in human cancer, produces two isoforms, KRAS4a and KRAS4b, through alternative splicing. These isoforms differ in exon 4, which encodes the final 15 residues of the G-domain and hypervariable regions (HVRs), vital for trafficking and membrane localization. While KRAS4b has been extensively studied, KRAS4a has been largely overlooked. Our multidisciplinary study compared the structural and functional characteristics of KRAS4a and KRAS4b, revealing distinct structural properties and thermal stability. Position 151 influences KRAS4a's thermal stability, while position 153 affects binding to RAF1 CRD protein. Nuclear magnetic resonance analysis identified localized structural differences near sequence variations and provided a solution-state conformational ensemble. Notably, KRAS4a exhibits substantial transcript abundance in bile ducts, liver, and stomach, with transcript levels approaching KRAS4b in the colon and rectum. Functional disparities were observed in full-length KRAS variants, highlighting the impact of HVR variations on interaction with trafficking proteins and downstream effectors like RAF and PI3K within cells.

Perspectives on Nuclear Magnetic Resonance Spectroscopy in Drug Discovery Research

The drug discovery landscape has undergone a significant transformation over the past decade, owing to research endeavors in a wide range of areas leading to strategies for pursuing new drug targets and the emergence of novel drug modalities. NMR spectroscopy has been a technology of fundamental importance to these research pursuits and has seen its use expanded both within and outside of traditional medicinal chemistry applications. In this perspective, we will present advancement of NMR-derived methods that have facilitated the characterization of small molecules and novel drug modalities including macrocyclic peptides, cyclic dinucleotides, and ligands for protein degradation. We will discuss innovations in NMR spectroscopy at the chemistry and biology interface that have broadened NMR's utility from hit identification through lead optimization activities. We will also discuss the promise of emerging NMR approaches in bridging our understanding and addressing challenges in the pursuit of the therapeutic agents of the future.

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.

Production of Isotopically Labeled KRAS4b

Isotopically labelled proteins are important reagents in structural biology as well as in targeted drug development. The field continues to advance with complex multi-isotope labeling. We have combined our experience in high-level soluble KRAS4b expression with protocols for isotope incorporation, to achieve reliable and efficient approaches for several labeling strategies. Typical experiments achieve nearly 100% 15N incorporation, with yields in the range of 1.3-24.6 mg/L (median = 6.4 mg/L, n = 53). Improvements in the growth parameters in the presence of deuterium reduce the standard time of fermentation from 5 days to 3 days by modifying the medium used during the weaning process. The methods described are compatible with multi-isotope labeling and site-specific labeling.

1D and 2D NMR for KRAS:Ligand Binding

Fragment-based screening by ligand-observed 1D NMR and binding interface mapping by protein-observed 2D NMR are popular methods used in drug discovery. These methods allow researchers to detect compound binding over a wide range of affinities and offer a simultaneous assessment of solubility, purity, and chemical formula accuracy of the target compounds and the 15N-labeled protein when examined by 1D and 2D NMR, respectively. These methods can be applied for screening fragment binding to the active (GMPPNP-bound) and inactive (GDP-bound) states of oncogenic KRAS mutants.

Targeting the RAS/RAF/MAPK pathway for cancer therapy: from mechanism to clinical studies

Metastatic dissemination of solid tumors, a leading cause of cancer-related mortality, underscores the urgent need for enhanced insights into the molecular and cellular mechanisms underlying metastasis, chemoresistance, and the mechanistic backgrounds of individuals whose cancers are prone to migration. The most prevalent signaling cascade governed by multi-kinase inhibitors is the mitogen-activated protein kinase (MAPK) pathway, encompassing the RAS-RAF-MAPK kinase (MEK)-extracellular signal-related kinase (ERK) pathway. RAF kinase is a primary mediator of the MAPK pathway, responsible for the sequential activation of downstream targets, such as MEK and the transcription factor ERK, which control numerous cellular and physiological processes, including organism development, cell cycle control, cell proliferation and differentiation, cell survival, and death. Defects in this signaling cascade are associated with diseases such as cancer. RAF inhibitors (RAFi) combined with MEK blockers represent an FDA-approved therapeutic strategy for numerous RAF-mutant cancers, including melanoma, non-small cell lung carcinoma, and thyroid cancer. However, the development of therapy resistance by cancer cells remains an important barrier. Autophagy, an intracellular lysosome-dependent catabolic recycling process, plays a critical role in the development of RAFi resistance in cancer. Thus, targeting RAF and autophagy could be novel treatment strategies for RAF-mutant cancers. In this review, we delve deeper into the mechanistic insights surrounding RAF kinase signaling in tumorigenesis and RAFi-resistance. Furthermore, we explore and discuss the ongoing development of next-generation RAF inhibitors with enhanced therapeutic profiles. Additionally, this review sheds light on the functional interplay between RAF-targeted therapies and autophagy in cancer.

Discovery of Hit Compounds Targeting the P4 Allosteric Site of K-RAS, Identified through Ensemble-Based Virtual Screening

Mutants of Ras are oncogenic drivers of a large number of human tumors. Despite being recognized as an attractive target for the treatment of cancer, the high affinity for its substrate tagged the protein as undruggable for a few years. The identification of cryptic pockets on the protein surface gave the opportunity to identify molecules capable of acting as allosteric modulators. Several molecules were disclosed in recent years, with sotorasib and adagrasib already approved for clinical use. The present study makes use of computational methods to characterize eight prospective allosteric pockets (P1-P8) in K-Ras, four of which (P1-P4) were previously characterized in the literature. The present study also describes the results of a virtual screening study focused on the discovery of hit compounds, binders of the P4 site that can be considered as peptidomimetics of a fragment of the SOS αI helix, a guanine exchange factor of Ras. After a detailed description of the computational procedure followed, we disclose five hit compounds, prospective binders of the P4 allosteric site that exhibit an inhibitory capability higher than 30% in a cell proliferation assay at 50 μM.

Allosteric site variants affect GTP hydrolysis on Ras

RAS GTPases are proto-oncoproteins that regulate cell growth, proliferation, and differentiation in response to extracellular signals. The signaling functions of RAS, and other small GTPases, are dependent on their ability to cycle between GDP-bound and GTP-bound states. Structural analyses suggest that GTP hydrolysis catalyzed by HRAS can be regulated by an allosteric site located between helices 3, 4, and loop 7. Here we explore the relationship between intrinsic GTP hydrolysis on HRAS and the position of helix 3 and loop 7 through manipulation of the allosteric site, showing that the two sites are functionally connected. We generated several hydrophobic mutations in the allosteric site of HRAS to promote shifts in helix 3 relative to helix 4. By combining crystallography and enzymology to study these mutants, we show that closure of the allosteric site correlates with increased hydrolysis of GTP on HRAS in solution. Interestingly, binding to the RAS binding domain of RAF kinase (RAF-RBD) inhibits GTP hydrolysis in the mutants. This behavior may be representative of a cluster of mutations found in human tumors, which potentially cooperate with RAF complex formation to stabilize the GTP-bound state of RAS.

Direct Modulators of K-Ras-Membrane Interactions

Protein-membrane interactions (PMIs) are ubiquitous in cellular signaling. Initial steps of signal transduction cascades often rely on transient and dynamic interactions with the inner plasma membrane leaflet to populate and regulate signaling hotspots. Methods to target and modulate these interactions could yield attractive tool compounds and drug candidates. Here, we demonstrate that the conjugation of a medium-chain lipid tail to the covalent K-Ras(G12C) binder MRTX849 at a solvent-exposed site enables such direct modulation of PMIs. The conjugated lipid tail interacts with the tethered membrane and changes the relative membrane orientation and conformation of K-Ras(G12C), as shown by molecular dynamics (MD) simulation-supported NMR studies. In cells, this PMI modulation restricts the lateral mobility of K-Ras(G12C) and disrupts nanoclusters. The described strategy could be broadly applicable to selectively modulate transient PMIs.

Reduced dynamic complexity allows structure elucidation of an excited state of KRASG13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS^G13D, which suggests an exchange between two conformational states in solution. Methyl and ³¹P NMR spectra of active KRAS^G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major P^γ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS^G13D and KRAS^G13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS^G13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

Targeting small GTPases: emerging grasps on previously untamable targets, pioneered by KRAS

Small GTPases including Ras, Rho, Rab, Arf, and Ran are omnipresent molecular switches in regulating key cellular functions. Their dysregulation is a therapeutic target for tumors, neurodegeneration, cardiomyopathies, and infection. However, small GTPases have been historically recognized as "undruggable". Targeting KRAS, one of the most frequently mutated oncogenes, has only come into reality in the last decade due to the development of breakthrough strategies such as fragment-based screening, covalent ligands, macromolecule inhibitors, and PROTACs. Two KRAS^G12C covalent inhibitors have obtained accelerated approval for treating KRAS^G12C mutant lung cancer, and allele-specific hotspot mutations on G12D/S/R have been demonstrated as viable targets. New methods of targeting KRAS are quickly evolving, including transcription, immunogenic neoepitopes, and combinatory targeting with immunotherapy. Nevertheless, the vast majority of small GTPases and hotspot mutations remain elusive, and clinical resistance to G12C inhibitors poses new challenges. In this article, we summarize diversified biological functions, shared structural properties, and complex regulatory mechanisms of small GTPases and their relationships with human diseases. Furthermore, we review the status of drug discovery for targeting small GTPases and the most recent strategic progress focused on targeting KRAS. The discovery of new regulatory mechanisms and development of targeting approaches will together promote drug discovery for small GTPases.

Dynamic regulation of RAS and RAS signaling

RAS proteins regulate most aspects of cellular physiology. They are mutated in 30% of human cancers and 4% of developmental disorders termed Rasopathies. They cycle between active GTP-bound and inactive GDP-bound states. When active, they can interact with a wide range of effectors that control fundamental biochemical and biological processes. Emerging evidence suggests that RAS proteins are not simple on/off switches but sophisticated information processing devices that compute cell fate decisions by integrating external and internal cues. A critical component of this compute function is the dynamic regulation of RAS activation and downstream signaling that allows RAS to produce a rich and nuanced spectrum of biological outputs. We discuss recent findings how the dynamics of RAS and its downstream signaling is regulated. Starting from the structural and biochemical properties of wild-type and mutant RAS proteins and their activation cycle, we examine higher molecular assemblies, effector interactions and downstream signaling outputs, all under the aspect of dynamic regulation. We also consider how computational and mathematical modeling approaches contribute to analyze and understand the pleiotropic functions of RAS in health and disease.

A simple protocol for the production of highly deuterated proteins for biophysical studies

Highly deuterated protein samples expand the biophysics and biological tool kit by providing, among other qualities, contrast matching in neutron diffraction experiments and reduction of dipolar spin interactions from normally protonated proteins in magnetic resonance studies, impacting both electron paramagnetic resonance and NMR spectroscopy. In NMR applications, deuteration is often combined with other isotopic labeling patterns to expand the range of conventional NMR spectroscopy research in both solution and solid-state conditions. However, preparation of deuterated proteins is challenging. We present here a simple, effective, and user-friendly protocol to produce highly deuterated proteins in Escherichia coli cells. The protocol utilizes the common shaker flask growth method and the well-known pET system (which provides expression control via the T7 promotor) for large-scale recombinant protein expression. One liter expression typically yields 5 to 50 mg of highly deuterated protein. Our data demonstrate that the optimized procedure produces a comparable quantity of protein in deuterium (²H₂O) oxide M9 medium compared with that in ¹H₂O M9 medium. The protocol will enable a broader utilization of deuterated proteins in a number of biophysical techniques.