Aggregation of lipid-anchored full-length H-Ras in lipid bilayers: simulations with the MARTINI force field

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

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.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

RAS GTPases and Interleaflet Coupling in the Plasma Membrane

RAS genes are frequently mutated in cancer. The primary signaling compartment of wild-type and constitutively active oncogenic mutant RAS proteins is the inner leaflet of the plasma membrane (PM). Thus, a better understanding of the unique environment of the PM inner leaflet is important to shed further light on RAS function. Over the past few decades, an integrated approach of superresolution imaging, molecular dynamic simulations, and biophysical assays has yielded new insights into the capacity of RAS proteins to sort lipids with specific headgroups and acyl chains, to assemble signaling nanoclusters on the inner PM. RAS proteins also sense and respond to changes in components of the outer PM leaflet, including glycophosphatidylinositol-anchored proteins, sphingophospholipids, glycosphingolipids, and galectins, as well as cholesterol that translocates between the two leaflets. Such communication between the inner and outer leaflets of the PM, called interleaflet coupling, allows RAS to potentially integrate extracellular mechanical and electrostatic information with intracellular biochemical signaling events, and reciprocally allows mutant RAS-transformed tumor cells to modify tumor microenvironments. Here, we review RAS-lipid interactions and speculate on potential mechanisms that allow communication between the opposing leaflets of the PM.

Measuring pico-Newton Forces with Lipid Anchors as Force Sensors in Molecular Dynamics Simulations

Binding forces between biomolecules are ubiquitous in nature but sometimes as weak as a few pico-Newtons (pN). In many cases, the binding partners are attached to biomembranes with the help of a lipid anchor. One important example are glycolipids that promote membrane adhesion through weak carbohydrate-carbohydrate binding between adjacent membranes. Here, we use molecular dynamics (MD) simulations to quantify the forces generated by bonds involving membrane-anchored molecules. We introduce a method in which the protrusion of the lipid anchors from the membrane acts as the force sensor. Our results with two different glycolipids reveal binding forces of up to 20 pN and corroborate the recent notion that carbohydrate-carbohydrate interactions are generic rather than specific.

Remodeling of the Plasma Membrane by Surface-Bound Protein Monomers and Oligomers: The Critical Role of Intrinsically Disordered Regions

The plasma membrane (PM) of cells is a dynamic structure whose morphology and composition is in constant flux. PM morphologic changes are particularly relevant for the assembly and disassembly of signaling platforms involving surface-bound signaling proteins, as well as for many other mechanochemical processes that occur at the PM surface. Surface-bound membrane proteins (SBMP) require efficient association with the PM for their function, which is often achieved by the coordinated interactions of intrinsically disordered regions (IDRs) and globular domains with membrane lipids. This review focuses on the role of IDR-containing SBMPs in remodeling the composition and curvature of the PM. The ability of IDR-bearing SBMPs to remodel the Gaussian and mean curvature energies of the PM is intimately linked to their ability to sort subsets of phospholipids into nanoclusters. We therefore discuss how IDRs of many SBMPs encode lipid-binding specificity or facilitate cluster formation, both of which increase their membrane remodeling capacity, and how SBMP oligomers alter membrane shape by monolayer surface area expansion and molecular crowding.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

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.

Approaches to inhibiting oncogenic K-Ras

Activating somatic K-Ras mutations are associated with >15% all human tumors and up to 90% of specific tumor types such as pancreatic cancer. Successfully inhibiting abnormal K-Ras signaling would therefore be a game changer in cancer therapy. However, K-Ras has long been considered an undruggable target for various reasons. This view is now changing by the discovery of allosteric inhibitors that directly target K-Ras and inhibit its functions, and by the identification of new mechanisms to dislodge it from the plasma membrane and thereby abrogate its cellular activities. In this review, we will discuss recent progresses and challenges to inhibiting aberrant K-Ras functions by these two approaches. We will also provide a broad overview of other approaches such as inhibition of K-Ras effectors, and offer a brief perspective on the way forward.

Estimating the accuracy of the MARTINI model towards the investigation of peripheral protein–membrane interactions

In this article, we investigate the ability of the MARTINI CG force field, specifically the 3 open-beta version, to reproduce known experimental observations regarding the membrane binding behavior of 12 peripheral membrane proteins and peptides.

Guardians of the Cell: State-of-the-Art of Membrane Proteins from a Computational Point-of-View

Membrane proteins (MPs) encompass a large family of proteins with distinct cellular functions, and although representing over 50% of existing pharmaceutical drug targets, their structural and functional information is still very scarce. Over the last years, in silico analysis and algorithm development were essential to characterize MPs and overcome some limitations of experimental approaches. The optimization and improvement of these methods remain an ongoing process, with key advances in MPs' structure, folding, and interface prediction being continuously tackled. Herein, we discuss the latest trends in computational methods toward a deeper understanding of the atomistic and mechanistic details of MPs.

Advances in understanding and in multi-disciplinary methodology used to assess lipid regulation of signalling cascades from the cancer cell plasma membrane

The lipid bilayer is a functional component of cells, forming a stable platform for the initiation of key biological processes, including cell signalling. There are distinct changes in the lipid composition of cell membranes during oncogenic transformation resulting in aberrant activation and inactivation of signalling transduction pathways. Studying the role of the cell membrane in cell signalling is challenging, since techniques are often limited to by timescale, resolution, sensitivity, and averaging. To overcome these limitations, combining 'computational', 'wet-lab' and 'semi-dry' approaches offers the best opportunity to resolving complex biological processes involved in membrane organisation. In this review, we highlight analytical tools that have been applied for the study of cell signalling initiation from the cancer cell membranes through computational microscopy, biological assays, and membrane biophysics. The cancer therapeutic potential of extracellular membrane-modulating agents, such as cholesterol-reducing agents is also discussed, as is the need for future collaborative inter-disciplinary research for studying the role of the cell membrane and its components in cancer therapy.

Probing the Conformational and Energy Landscapes of KRAS Membrane Orientation

Membrane reorientation of oncogenic RAS proteins is emerging as an important modulator of their functions. Previous studies have shown that the most common orientations include those with either the three C-terminal α-helices (OS1) or N-terminal β-strands (OS2) of the catalytic domain facing the membrane. OS1 and OS2 differ by the degree to which the effector-interacting surface is occluded by the membrane. However, the relative stability of these states and the rates of transition between them remained undetermined. How mutations might modulate preferences for specific orientation states is also far from clear. The current work attempted to address these questions through a comprehensive analysis of two 20 μs-long atomistic molecular dynamics simulations. The simulations were conducted on the oncogenic G12D and Q61H KRAS mutants bound to an anionic lipid bilayer. G12D and Q61H are among the most prevalent cancer-causing mutations at the P-loop and switch 2 regions of KRAS, respectively. We found that both mutants fluctuate in a similar manner between OS1 and OS2 via an intermediate orientation OS0, and both favor the signaling competent OS1 and OS0 over the occluded OS2. However, they differ in the details, such as in the extent to which they sample OS1. Analysis of the orientation free-energy landscapes estimated from the simulations indicate that OS1 and OS2 are the most stable states. However, the overall free energy surface is rugged, indicating a large diversity of conformations including at least two substates in each orientation state that differ in stability only by about 0.5-1.0 kcal/mol. Reversible transitions between OS1 and OS2 occur via two well-defined pathways that traverse OS0. In the minimum energy path, helix 4 remains close to the membrane as the angle of the catalytic domain from the membrane plane changes, resulting in a barrier of ∼1 kcal/mol for OS1/OS2 interconversions. Estimation of the rates of the various transitions based on survival probabilities yielded two rate constants in the order of 10⁷ and 10⁶ s^-1, which we attribute to intrinsic protein conformational dynamics and transient protein-lipid interactions, respectively. The faster process dominates every transition, confirming a previous suggestion that RAS membrane reorientation is driven by conformational fluctuations rather than protein-lipid interactions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Interplay between Membrane Curvature and Cholesterol: Role of Palmitoylated Caveolin-1

Caveolin-1 (cav-1) is an important player in cell signaling and endocytosis that has been shown to colocalize with cholesterol-rich membrane domains. Experimental studies with varying cav-1 constructs have suggested that it can induce both cholesterol clustering and membrane curvature. Here, we probe the molecular origin of membrane curvature and cholesterol clustering by cav-1 by using coarse-grain molecular dynamics simulations. We have performed a series of simulations of a functionally important cav-1 construct, comprising the membrane-interacting domains and a C-terminal palmitoyl tail. Our results suggest that cav-1 is able to induce cholesterol clustering in the membrane leaflet to which it is bound as well as the opposing leaflet. A positive membrane curvature is observed upon cav-1 binding in cholesterol-containing bilayers. Interestingly, we observe an interplay between cholesterol clustering and membrane curvature such that cav-1 is able to induce higher membrane curvature in cholesterol-rich membranes. The role of the cav-1 palmitoyl tail is less clear and appears to increase the membrane contacts. Further, we address the importance of the secondary structure of cav-1 domains and show that it could play an important role in membrane curvature and cholesterol clustering. Our work is an important step toward a molecular picture of caveolae and vesicular endocytosis.

Interactions of the EphA2 Kinase Domain with PIPs in Membranes: Implications for Receptor Function

EphA2 is a member of the receptor tyrosine kinase family. Interactions of the cytoplasmic region of EphA2 with the cell membrane are functionally important and yet remain incompletely characterized. Molecular dynamics simulations combined with biochemical studies reveal the interactions of the transmembrane, juxtamembrane (JM), and kinase domains with the membrane. We describe how the kinase domain is oriented relative to the membrane and how the JM region can modulate this interaction. We highlight the role of phosphatidylinositol phosphates (PIPs) in mediating the interaction of the kinase domain with the membrane and, conversely, how positively charged patches at the kinase surface and in the JM region induce the formation of nanoclusters of PIP molecules in the membrane. Integration of these results with those from previous studies enable computational reconstitution of a near complete EphA2 receptor within a membrane, suggesting a role for receptor-lipid interactions in modulation of EphA2.

Self-Assembled Ligands Targeting TLR7: A Molecular Level Investigation

Toll-like receptors (TLRs) are pattern recognition transmembrane proteins that play an important role in innate immunity. In particular, TLR7 plays a role in detecting nucleic acids derived from viruses and bacteria. The huge number of pathologies in which TLR7 is involved has led to an increasing interest in developing new compounds targeting this protein. Several conjugation strategies were proposed for TLR7 agonists to increase the potency while maintaining a low toxicity. In this work, we focus the attention on two promising classes of TLR7 compounds derived from the same pharmacophore conjugated with phospholipid and polyethylene glycol (PEG). A multidisciplinary investigation has been carried out by molecular dynamics (MD), dynamic light scattering (DLS), electron paramagnetic resonance (EPR), and cytotoxicity assessment. DLS and MD indicated how only the phospholipid conjugation provides the compound abilities to self-assemble in an orderly fashion with a maximal pharmacophore exposition to the solvent. Further EPR and cytotoxicity experiments highlighted that phospholipid compounds organize in stable aggregates and well interact with TLR7, whereas PEG conjugation was characterized by poorly stable aggregates at the cells surface. The methodological framework proposed in this study may be used to investigate, at a molecular level, the interactions generally occurring between aggregated ligands, to be used as drugs, and protein receptors.

Spatiotemporal Analysis of K-Ras Plasma Membrane Interactions Reveals Multiple High Order Homo-oligomeric Complexes

Self-assembly of plasma membrane-associated Ras GTPases has major implications to the regulation of cell signaling. However, the structural basis of homo-oligomerization and the fractional distribution of oligomeric states remained undetermined. We have addressed these issues by deciphering the distribution of dimers and higher-order oligomers of K-Ras4B, the most frequently mutated Ras isoform in human cancers. We focused on the constitutively active G12V K-Ras and two of its variants, K101E and K101C/E107C, which respectively destabilize and stabilize oligomers. Using raster image correlation spectroscopy and number and brightness analysis combined with fluorescence recovery after photobleaching, fluorescence correlation spectroscopy and electron microscopy in live cells, we show that G12V K-Ras exists as a mixture of monomers, dimers and larger oligomers, while the K101E mutant is predominantly monomeric and K101C/E107C is dominated by oligomers. This observation demonstrates the ability of K-Ras to exist in multiple oligomeric states whose population can be altered by interfacial mutations. Using molecular modeling and simulations we further show that K-Ras uses two partially overlapping interfaces to form compositionally and topologically diverse oligomers. Our results thus provide the first detailed insight into the multiplicity, structure, and membrane organization of K-Ras homomers.

Geodesic curvature driven surface microdomain formation

Lipid bilayer membranes are not uniform and clusters of lipids in a more ordered state exist within the generally disorder lipid milieu of the membrane. These clusters of ordered lipids microdomains are now referred to as lipid rafts. Recent reports attribute the formation of these microdomains to the geometrical and molecular mechanical mismatch of lipids of different species on the boundary. Here we introduce the geodesic curvature to characterize the geometry of the domain boundary, and develop a geodesic curvature energy model to describe the formation of these microdomains as a result of energy minimization. Our model accepts the intrinsic geodesic curvature of any binary lipid mixture as an input, and will produce microdomains of the given geodesic curvature as demonstrated by three sets of numerical simulations. Our results are in contrast to the surface phase separation predicted by the classical surface Cahn-Hilliard equation, which tends to generate large domains as a result of the minimizing line tension. Our model provides a direct and quantified description of the structure inhomogeneity of lipid bilayer membrane, and can be coupled to the investigations of biological processes on membranes for which such inhomogeneity plays essential roles.

Lipoprotein insertion into membranes of various complexity: lipid sorting, interfacial adsorption and protein clustering

In a combined chemical-biological and biophysical approach we explored the membrane partitioning of the lipidated signaling proteins N-Ras and K-Ras4B into membrane systems of different complexity, ranging from one-component lipid bilayers and anionic binary and ternary heterogeneous membrane systems even up to partitioning studies on protein-free and protein-containing giant plasma membrane vesicles (GPMVs). To yield a pictorial view of the localization process, imaging using confocal laser scanning and atomic force microscopy was performed. The results reveal pronounced isoform-specific differences regarding the lateral distribution and formation of protein-rich membrane domains. Line tension is one of the key parameters controlling not only the size and dynamic properties of segregated lipid domains but also the partitioning process of N-Ras that acts as a lineactant. The formation of N-Ras protein clusters is even recorded for almost vanishing hydrophobic mismatch. Conversely, for K-Ras4B, selective localization and clustering are electrostatically mediated by its polybasic farnesylated C-terminus. The formation of K-Ras4B clusters is also observed for the multi-component GPMV membrane, i.e., it seems to be a general phenomenon, largely independent of the details of the membrane composition, including the anionic charge density of lipid headgroups. Our data indicate that unspecific and entropy-driven membrane-mediated interactions play a major role in the partitioning behavior, thus relaxing the need for a multitude of fine-tuned interactions. Such a scenario seems also to be reasonable recalling the high dynamic nature of cellular membranes. Finally, we note that even relatively simple models of heterogeneous membranes are able to reproduce many of the properties of much more complex biological membranes.

Domain Stability in Biomimetic Membranes Driven by Lipid Polyunsaturation

Biological membranes contain a broad variety of lipid species whose individual physicochemical properties and collective interactions ultimately determine membrane organization. A key aspect of the organization of cellular membranes is their lateral subdivision into domains of distinct structure and composition. The most widely studied membrane domains are lipid rafts, which are the biological manifestations of liquid-ordered phases that form in sterol-containing membranes. Detailed studies of biomimetic membrane mixtures have yielded wide-ranging insights into the physical principles behind lipid rafts; however, these simplified models do not fully capture the diversity and complexity of the mammalian lipidome, most notably in their exclusion of polyunsaturated lipids. Here, we assess the role of lipid acyl chain unsaturation as a driving force for phase separation using coarse-grained molecular dynamics (CGMD) simulations validated by model membrane experiments. The clear trends in our observations and good qualitative agreements between simulations and experiments support the conclusions that highly unsaturated lipids promote liquid-liquid domain stability by enhancing the differences in cholesterol content and lipid chain order between the coexisting domains. These observations reveal the important role of noncanonical biological lipids in the physical properties of membranes, showing that lipid polyunsaturation is a driving force for liquid-liquid phase separation.

Ras Dimer Formation as a New Signaling Mechanism and Potential Cancer Therapeutic Target

The K-, N-, and HRas small GTPases are key regulators of cell physiology and are frequently mutated in human cancers. Despite intensive research, previous efforts to target hyperactive Ras based on known mechanisms of Ras signaling have been met with little success. Several studies have provided compelling evidence for the existence and biological relevance of Ras dimers, establishing a new mechanism for regulating Ras activity in cells additionally to GTP-loading and membrane localization. Existing data also start to reveal how Ras proteins dimerize on the membrane. We propose a dimer model to describe Ras-mediated effector activation, which contrasts existing models of Ras signaling as a monomer or as a 5-8 membered multimer. We also discuss potential implications of this model in both basic and translational Ras biology.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.

Specific cancer-associated mutations in the switch III region of Ras increase tumorigenicity by nanocluster augmentation

Hotspot mutations of Ras drive cell transformation and tumorigenesis. Less frequent mutations in Ras are poorly characterized for their oncogenic potential. Yet insight into their mechanism of action may point to novel opportunities to target Ras. Here, we show that several cancer-associated mutations in the switch III region moderately increase Ras activity in all isoforms. Mutants are biochemically inconspicuous, while their clustering into nanoscale signaling complexes on the plasma membrane, termed nanocluster, is augmented. Nanoclustering dictates downstream effector recruitment, MAPK-activity, and tumorigenic cell proliferation. Our results describe an unprecedented mechanism of signaling protein activation in cancer.

DOI:http://dx.doi.org/10.7554/eLife.08905.001

Cancer is a disease that develops when cells within the body acquire genetic mutations that allow them to grow and divide rapidly. Many human cancers have mutations in a gene that encodes a protein called Ras, which promotes cell growth and division by controlling the activities of other proteins.

Ras congregates at the membrane that surrounds the cell and can assemble into clusters (called nanoclusters) that each contain six to eight Ras proteins. The tight packing of the proteins in these nanoclusters increases the amount of Ras in the membrane locally, which allows Ras to interact with other proteins more efficiently to promote growth and cell division.

In normal cells, other proteins control when Ras is active. However, in many cancer cells, Ras is active all the time due to mutations that occur in three ‘hotspots’ within its gene. Other mutations in the gene that encodes Ras are also found in cancer cells, but these are less common and it is not clear how they alter the activity of the protein.

Here, Solman et al. used microscopy and biochemical techniques to study the effects of some of the less common mutations on Ras activity in human cells. The experiments show that several mutations that alter a region of Ras called the ‘switch III region’ moderately increase the activity of Ras. The mutations probably alter the way that Ras sits in the membrane, which in turn changes the way it interacts with other proteins and the membrane so that more Ras nanoclusters form.

Solman et al.'s findings reveal a new way that Ras can be activated in cancer cells. The next challenge is to develop drugs that block the formation of Ras nanoclusters and to find out if they have the potential to be used to treat cancer.

DOI:http://dx.doi.org/10.7554/eLife.08905.002

Modulation of a small two-domain lipid vesicle by linactants

Linactants, molecules that preferentially localize at the boundary of lipid membrane domains, are attracting considerable attention in recent years due to the recognition that they might regulate lipid-phase separation and thereby modulate membrane morphology. Recent studies have also shown that clustering of some line active agents enhances their ability to modulate membrane curvature. However, the molecular origin of this phenomenon, and the degree to which it impacts biological membranes, remains poorly understood. In this work, we have investigated how linactants induce shape change in multidomain small unilamallar vesicles (SUVs) using extensive dissipative particle dynamics simulations. The linactant was modeled as a two-tailed hybrid lipid with the two tails differing in preference for different lipid domains. We found that addition of a small amount of linactants (∼1%) to a two-domain vesicle leads to substantial reduction in the line tension and neck curvature at the domain boundary. Using cross-linking as a surrogate for clustering, we further show that linactant clusters substantially enhance the boundary preference and therefore the reduction in neck curvature. Moreover, on the basis of analyses of the corresponding changes in the membrane energetics, we highlight how linactants might stabilize nanoscale domains. These results have important implications for the potential existence and physical explanations of nanosized domains in biological membranes.

Simulations of Protein Aggregation: Insights from Atomistic and Coarse-Grained Models

This Perspective highlights recent computational approaches to protein aggregation, from coarse-grained models to atomistic simulations, using the islet amyloid polypeptide (IAPP) as a case study. We review salient open questions where simulations can make an impact, discuss the successes and challenges met by simulations, and explore new directions.

Membrane Remodeling by Surface-Bound Protein Aggregates: Insights from Coarse-Grained Molecular Dynamics Simulation

The mechanism of curvature generation in membranes has been studied for decades due to its important role in many cellular functions. However, it is not clear if, or how, aggregates of lipid-anchored proteins might affect the geometry and elastic property of membranes. As an initial step toward addressing this issue, we performed structural, geometrical, and stress field analyses of coarse-grained molecular dynamics trajectories of a domain-forming bilayer in which an aggregate of lipidated proteins was asymmetrically bound. The results suggest a general mechanism whereby asymmetric incorporation of lipid-modified protein aggregates curve multidomain membranes primarily by expanding the surface area of the monolayer in which the lipid anchor is inserted.

Overview of simulation studies on the enzymatic activity and conformational dynamics of the GTPase Ras

Over the last 40 years, we have learnt a great deal about the Ras onco-proteins. These intracellular molecular switches are essential for the function of a variety of physiological processes, including signal transduction cascades responsible for cell growth and proliferation. Molecular simulations and free energy calculations have played an essential role in elucidating the conformational dynamics and energetics underlying the GTP hydrolysis reaction catalysed by Ras. Here we present an overview of the main lessons from molecular simulations on the GTPase reaction and conformational dynamics of this important anti-cancer drug target. In the first part, we summarise insights from quantum mechanical and combined quantum mechanical/molecular mechanical simulations as well as other free energy methods and highlight consensus viewpoints as well as remaining controversies. The second part provides a very brief overview of new insights emerging from large-scale molecular dynamics simulations. We conclude with a perspective regarding future studies of Ras where computational approaches will likely play an active role.

The efficacy of Raf kinase recruitment to the GTPase H-ras depends on H-ras membrane conformer-specific nanoclustering

Solution structures and biochemical data have provided a wealth of mechanistic insight into Ras GTPases. However, information on how much the membrane organization of these lipid-modified proteins impacts on their signaling is still scarce. Ras proteins are organized into membrane nanoclusters, which are necessary for Ras-MAPK signaling. Using quantitative conventional and super-resolution fluorescence methods, as well as mathematical modeling, we investigated nanoclustering of H-ras helix α4 and hypervariable region mutants that have different bona fide conformations on the membrane. By following the emergence of conformer-specific nanoclusters in the plasma membrane of mammalian cells, we found that conformers impart distinct nanoclustering responses depending on the cytoplasmic levels of the nanocluster scaffold galectin-1. Computational modeling revealed that complexes containing H-ras conformers and galectin-1 affect both the number and lifetime of nanoclusters and thus determine the specific Raf effector recruitment. Our results show that mutations in Ras can affect its nanoclustering response and thus allosterically effector recruitment and downstream signaling. We postulate that cancer- and developmental disease-linked mutations that are associated with the Ras membrane conformation may exhibit so far unrecognized Ras nanoclustering and therefore signaling alterations.

NRas slows the rate at which a model lipid bilayer phase separates

The Ras family of small membrane-associated GTP-ases are important components in many different cell signalling cascades. They are thought to cluster on the cell membrane through association with cholesterol-rich nanodomains. This process remains poorly understood. Here we test the effect of adding multiple copies of NRas, one of the canonical Ras proteins, to a three-component lipid bilayer that rapidly undergoes spinodal decomposition (i.e. unmixing), thereby creating ordered and disordered phases. Coarse-grained molecular dynamics simulations of a large bilayer containing 6000 lipids, with and without protein, are compared. NRas preferentially localises to the interface between the domains and slows the rate at which the domains grow. We infer that this doubly-lipidated cell signalling protein is reducing the line tension between the ordered and disordered regions. This analysis is facilitated by our use of techniques borrowed from image-processing. The conclusions above are contingent upon several assumptions, including the use of a model lipid with doubly unsaturated tails and the limited structural data available for the C-terminus of NRas, which is where the lipid anchors are found.