Direct Modulators of K-Ras-Membrane Interactions

Molecular dynamics simulations for the structure-based drug design: targeting small-GTPases proteins

INTRODUCTION

Molecular Dynamics (MD) simulations can support mechanism-based drug design. Indeed, MD simulations by capturing biomolecule motions at finite temperatures can reveal hidden binding sites, accurately predict drug-binding poses, and estimate the thermodynamics and kinetics, crucial information for drug discovery campaigns. Small-Guanosine Triphosphate Phosphohydrolases (GTPases) regulate a cascade of signaling events, that affect most cellular processes. Their deregulation is linked to several diseases, making them appealing drug targets. The broad roles of small-GTPases in cellular processes and the recent approval of a covalent KRas inhibitor as an anticancer agent renewed the interest in targeting small-GTPase with small molecules.

AREA COVERED

This review emphasizes the role of MD simulations in elucidating small-GTPase mechanisms, assessing the impact of cancer-related variants, and discovering novel inhibitors.

EXPERT OPINION

The application of MD simulations to small-GTPases exemplifies the role of MD simulations in the structure-based drug design process for challenging biomolecular targets. Furthermore, AI and machine learning-enhanced MD simulations, coupled with the upcoming power of quantum computing, are promising instruments to target elusive small-GTPases mutations and splice variants. This powerful synergy will aid in developing innovative therapeutic strategies associated to small-GTPases deregulation, which could potentially be used for personalized therapies and in a tissue-agnostic manner to treat tumors with mutations in small-GTPases.

Lysophosphatidylcholine acyltransferase 1 suppresses nanoclustering and function of KRAS

KRAS is frequently mutated in cancer, contributing to 20% of all human cancer especially pancreatic, colorectal and lung cancer. Signaling of the constitutively active KRAS oncogenic mutants is mostly compartmentalized to proteolipid nanoclusters on the plasma membrane (PM). Signaling nanoclusters of many KRAS mutants selectively enrich phosphatidylserine (PS) lipids with unsaturated sn-2 acyl chains, but not the fully saturated PS species. Thus, remodeling PS acyl chains may suppress KRAS oncogenesis. Lysophosphatidylcholine acyltransferases (LPCATs) remodel sn-2 acyl chains of phospholipids, with LPCAT1 preferentially generating the fully saturated lipids. Here, we show that stable expression of LPCAT1 depletes major PS species with unsaturated sn-2 chains while decreasing minor phosphatidylcholine (PC) species with the corresponding acyl chains. LPCAT1 expression more effectively disrupts the nanoclustering of oncogenic GFP-KRASG12V, which is restored by acute addback of exogenous unsaturated PS. LPCAT1 expression compromises signaling and oncogenic activities of the KRAS-dependent pancreatic tumor lines. LPCAT1 expression sensitizes human pancreatic tumor MiaPaCa-2 cells to KRASG12C specific inhibitor, Sotorasib. Statistical analyses of patient data further reveal that pancreatic cancer patients with KRAS mutations express less LPCAT1. Higher LPCAT1 expression also improves survival probability of pancreatic and lung adenocarcinoma patients with KRAS mutations. Thus, PS acyl chain remodeling selectively suppresses KRAS oncogenesis.

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

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

A Top-Down Proteomic Assay to Evaluate KRAS4B-Compound Engagement

Development of new targeted inhibitors for oncogenic KRAS mutants may benefit from insight into how a given mutation influences the accessibility of protein residues and how compounds interact with mutant or wild-type KRAS proteins. Targeted proteomic analysis, a key validation step in the KRAS inhibitor development process, typically involves both intact mass- and peptide-based methods to confirm compound localization or quantify binding. However, these methods may not always provide a clear picture of the compound binding affinity for KRAS, how specific the compound is to the target KRAS residue, and how experimental conditions may impact these factors. To address this, we have developed a novel top-down proteomic assay to evaluate in vitro KRAS4B-compound engagement while assessing relative quantitation in parallel. We present two applications to demonstrate the capabilities of our assay: maleimide-biotin labeling of a KRAS4BG12D cysteine mutant panel and treatment of three KRAS4B proteins (WT, G12C, and G13C) with small molecule compounds. Our results show the time- or concentration-dependence of KRAS4B-compound engagement in context of the intact protein molecule while directly mapping the compound binding site.

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

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

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.