Determinants of Membrane Orientation Dynamics in Lipid-Modified Small GTPases

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.

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.

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 these IDRs are key determinants of membrane localization, and interactions between the IDR and the membrane have been shown to affect signaling in RAS proteins through modulation of dynamic membrane organization. Here we utilized atomistic molecular dynamics simulations to study the membrane interactions, 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 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 generate semi-stable conformational sub-states, and a lack of these residues leads to greater conformational diversity. Finally, we show that small GTPase IDRs with a polybasic domain alter local lipid composition by segregating anionic membrane lipids, and, in some cases, excluding other lipids from their immediate proximity in favor of anionic lipids.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

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.

Direct Observation of Membrane-Associated H-Ras in the Native Cellular Environment by In-Cell 19F-NMR Spectroscopy

Ras acts as a molecular switch to control intracellular signaling on the plasma membrane (PM). Elucidating how Ras associates with PM in the native cellular environment is crucial for understanding its control mechanism. Here, we used in-cell nuclear magnetic resonance (NMR) spectroscopy combined with site-specific ¹⁹F-labeling to explore the membrane-associated states of H-Ras in living cells. The site-specific incorporation of p-trifluoromethoxyphenylalanine (OCF₃Phe) at three different sites of H-Ras, i.e., Tyr32 in switch I, Tyr96 interacting with switch II, and Tyr157 on helix α5, allowed the characterization of their conformational states depending on the nucleotide-bound states and an oncogenic mutational state. Exogenously delivered ¹⁹F-labeled H-Ras protein containing a C-terminal hypervariable region was assimilated via endogenous membrane-trafficking, enabling proper association with the cell membrane compartments. Despite poor sensitivity of the in-cell NMR spectra of membrane-associated H-Ras, the Bayesian spectral deconvolution identified distinct signal components on three ¹⁹F-labeled sites, thus offering the conformational multiplicity of H-Ras on the PM. Our study may be helpful in elucidating the atomic-scale picture of membrane-associated proteins in living cells.

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 guidance and adhesion protein FLRT2 dimerizes in cis via dual small-X3-small transmembrane motifs

Fibronectin Leucine-rich Repeat Transmembrane (FLRT 1-3) proteins are a family of broadly expressed single-spanning transmembrane receptors that play key roles in development. Their extracellular domains mediate homotypic cell-cell adhesion and heterotypic protein interactions with other receptors to regulate cell adhesion and guidance. These in trans FLRT interactions determine the formation of signaling complexes of varying complexity and function. Whether FLRTs also interact at the surface of the same cell, in cis, remains unknown. Here, molecular dynamics simulations reveal two dimerization motifs in the FLRT2 transmembrane helix. Single particle tracking experiments show that these Small-X₃-Small motifs synergize with a third dimerization motif encoded in the extracellular domain to permit the cis association and co-diffusion patterns of FLRT2 receptors on cells. These results may point to a competitive switching mechanism between in cis and in trans interactions, which suggests that homotypic FLRT interaction mirrors the functionalities of classic adhesion molecules.