Parameterization of Palmitoylated Cysteine, Farnesylated Cysteine, Geranylgeranylated Cysteine, and Myristoylated Glycine for the Martini Force Field

Caveolin assemblies displace one bilayer leaflet to organize and bend membranes

Caveolin is a monotopic integral membrane protein, widely expressed in metazoa and responsible for constructing enigmatic membrane invaginations known as caveolae. Recently, the high-resolution structure of a purified human caveolin assembly, the CAV1-8S complex, revealed a unique organization of 11 protomers arranged in a tightly packed, radially symmetric spiral disc. One face and the outer rim of this disc are highly hydrophobic, suggesting that the complex incorporates into membranes by displacing hundreds of lipids from one leaflet. The feasibility of this unique molecular architecture and its biophysical and functional consequences are currently unknown. Using Langmuir film balance measurements, we find that CAV1-8S is highly surface active and intercalates into lipid monolayers. Molecular simulations of biomimetic bilayers support this 'leaflet replacement' model and reveal that while CAV1-8S effectively displaces phospholipids from one bilayer leaflet, it accumulates 40-70 cholesterol molecules into a disordered monolayer between the complex and its distal lipid leaflet. We find that CAV1-8S preferentially associates with positively curved membrane surfaces due to its influence on the conformations of distal leaflet lipids, and that these effects laterally sort lipids of the distal leaflet. Large-scale simulations of multiple caveolin assemblies confirmed their association with large, positively curved membrane morphologies, consistent with the shape of caveolae. Further, association with curved membranes regulates the exposure of caveolin residues implicated in protein-protein interactions. Altogether, the unique structure of CAV1-8S imparts unusual modes of membrane interaction with implications for membrane organization, morphology, and physiology.

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.

Native mass spectrometry of proteoliposomes containing integral and peripheral membrane proteins

Cellular membranes are critical to the function of membrane proteins, whether they are associated (peripheral) or embedded (integral) within the bilayer. While detergents have contributed to our understanding of membrane protein structure and function, there remains challenges in characterizing protein-lipid interactions within the context of an intact membrane. Here, we developed a method to prepare proteoliposomes for native mass spectrometry (MS) studies. We first use native MS to detect the encapsulation of soluble proteins within liposomes. We then find the peripheral Gβ1γ2 complex associated with the membrane can be ejected and analyzed using native MS. Four different integral membrane proteins (AmtB, AqpZ, TRAAK, and TREK2), all of which have previously been characterized in detergent, eject from the proteoliposomes as intact complexes bound to lipids that have been shown to tightly associate in detergent, drawing a correlation between the two approaches. We also show the utility of more complex lipid environments, such as a brain polar lipid extract, and show TRAAK ejects from liposomes of this extract bound to lipids. These findings underscore the capability to eject protein complexes from membranes bound to both lipids and metal ions, and this approach will be instrumental in the identification of key protein-lipid interactions.

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.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Effects of surfactant adsorption on the wettability and friction of biomimetic surfaces

The properties of solid-liquid interfaces can be markedly altered by surfactant adsorption. Here, we use molecular dynamics (MD) simulations to study the adsorption of ionic surfactants at the interface between water and heterogeneous solid surfaces with randomly arranged hydrophilic and hydrophobic regions, which mimic the surface properties of human hair. We use the coarse-grained MARTINI model to describe both the hair surfaces and surfactant solutions. We consider negatively-charged virgin and bleached hair surface models with different grafting densities of neutral octadecyl and anionic sulfonate groups. The adsorption of cationic cetrimonium bromide (CTAB) and anionic sodium dodecyl sulfate (SDS) surfactants from water are studied above the critical micelle concentration. The simulated adsorption isotherms suggest that cationic surfactants adsorb to the surfaces via a two-stage process, initially forming monolayers and then bilayers at high concentrations, which is consistent with previous experiments. Anionic surfactants weakly adsorb via hydrophobic interactions, forming only monolayers on both virgin and medium bleached hair surfaces. We also conduct non-equilibrium molecular dynamics simulations, which show that applying cationic surfactant solutions to bleached hair successfully restores the low friction seen with virgin hair. Friction is controlled by the combined surface coverage of the grafted lipids and the adsorbed CTAB molecules. Treated surfaces containing monolayers and bilayers both show similar friction, since the latter are easily removed by compression and shear. Further wetting MD simulations show that bleached hair treated with CTAB increases the hydrophobicity to similar levels seen for virgin hair. Treated surfaces containing CTAB monolayers with the tailgroups pointing predominantly away from the surface are more hydrophobic than bilayers due to the electrostatic interactions between water molecules and the exposed cationic headgroups.

Nanoscale friction of biomimetic hair surfaces

We investigate the nanoscale friction between biomimetic hair surfaces using chemical colloidal probe atomic force microscopy experiments and nonequilibrium molecular dynamics simulations. In the experiments, friction is measured between water-lubricated silica surfaces functionalised with monolayers formed from either octadecyl or sulfonate groups, which are representative of the surfaces of virgin and ultimately bleached hair, respectively. In the simulations, friction is monitored between coarse-grained model hair surfaces with different levels of chemical damage, where a specified amount of grafted octadecyl groups are randomly replaced with sulfonate groups. The sliding velocity dependence of friction in the simulations can be described using an extended stress-augmented thermally activation model. As the damage level increases in the simulations, the friction coefficient generally increases, but its sliding velocity-dependence decreases. At low sliding velocities, which are closer to those encountered experimentally and physiologically, we observe a monotonic increase of the friction coefficient with damage ratio, which is consistent with our new experiments using biomimetic surfaces and previous ones using real hair. This observation demonstrates that modified surface chemistry, rather than roughness changes or subsurface damage, control the increase in nanoscale friction of bleached or chemically damaged hair. We expect the methods and biomimetic surfaces proposed here to be useful to screen the tribological performance of hair care formulations both experimentally and computationally.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure-function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck's activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations  using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck's conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

Molecular Mechanism of CD44 Homodimerization Modulated by Palmitoylation and Membrane Environments

The homodimerization of CD44 plays a key role in an intercellular-to-extracellular signal transduction and tumor progression. Acylated modification and specific membrane environments have been reported to mediate translocation and oligomerization of CD44, however, the underlying molecular mechanism remains elusive. In this study, extensive molecular dynamics simulations are performed to characterize the dimerization of palmitoylated CD44 variants in different bilayer environments. CD44 forms homodimer depending on the cysteines on the juxta-membrane domains, and the dimerization efficiency and packing configurations are defected by their palmitoylated modifications. In the phase-segregated (raft included) membrane, homodimerization of the palmitoylated CD44 is hardly observed, whereas PIP2 addition compensates to realize dimerization. However, the structure of CD44 homodimer formed in the phase-segregated bilayer turns susceptive and PIP2 addition allows for an extensive conformation of the cytoplasmic domain, a proposal prerequisite to access the cytoskeleton linker proteins. The results unravel a delicate competitive relationship between PIP2, lipid raft and palmitoylation in mediating protein homodimerization, which helps to clarify the dynamic dimer conformations and involved cellular signaling of the CD44 likewise proteins.

Ultralong recovery time in nanosecond electroporation systems enabled by orientational-disordering processes

The growing importance of applications based on molecular medicine and genetic engineering is driving the need to develop high-performance electroporation technologies. The electroporation phenomenon involves disruption of the cell for increasing membrane permeability. Although there is a multitude of research focused on exploring new electroporation techniques, the engineering of programming schemes suitable for these electroporation methods remains a challenge. Nanosecond stimulations could be promising candidates for these techniques owing to their ability to generate a wide range of biological responses. Here we control the membrane permeabilization of cancer cells using different numbers of electric-field pulses through orientational disordering effects. We then report our exploration of a few-volt nanosecond alternating-current (AC) stimulation method with an increased number of pulses for developing electroporation systems. A recovery time of ∼720 min was achieved, which is above the average of ∼76 min for existing electroporation methods using medium cell populations, as well as a previously unreported increased conductance with an increase in the number of pulses using weak bias amplitudes. All-atom molecular dynamics (MD) simulations reveal the orientation-disordering-facilitated increase in the degree of permeabilization. These findings highlight the potential of few-volt nanosecond AC-stimulation with an increased number of pulse strategies for the development of next-generation low-power electroporation systems.

Structure and dynamics of the SARS-CoV-2 envelope protein monomer

Coronaviruses, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), present an ongoing threat to human wellbeing. Consequently, elucidation of molecular determinants of their function and interaction with the host is an important task. Whereas some of the coronaviral proteins are extensively characterized, others remain understudied. Here, we use molecular dynamics simulations to analyze the structure and dynamics of the SARS-CoV-2 envelope (E) protein (a viroporin) in the monomeric form. The protein consists of the hydrophobic α-helical transmembrane domain (TMD) and amphiphilic α-helices H2 and H3, connected by flexible linkers. We show that TMD has a preferable orientation in the membrane, while H2 and H3 reside at the membrane surface. Orientation of H2 is strongly influenced by palmitoylation of cysteines Cys40, Cys43, and Cys44. Glycosylation of Asn66 affects the orientation of H3. We also observe that the monomeric E protein both generates and senses the membrane curvature, preferably localizing with the C-terminus at the convex regions of the membrane; the protein in the pentameric form displays these properties as well. Localization to curved regions may be favorable for assembly of the E protein oligomers, whereas induction of curvature may facilitate the budding of the viral particles. The presented results may be helpful for a better understanding of the function of the coronaviral E protein and viroporins in general, and for overcoming the ongoing SARS-CoV-2 pandemic.

Coarse-grained molecular models of the surface of hair

We present a coarse-grained molecular model of the surface of human hair, which consists of a supported lipid monolayer, in the MARTINI framework. Using coarse-grained molecular dynamics (MD) simulations, we identify a lipid grafting distance that yields a monolayer thickness consistent with both atomistic MD simulations and experimental measurements of the hair surface. Coarse-grained models for fully-functionalised, partially damaged, and fully damaged hair surfaces are created by randomly replacing neutral thioesters with anionic sulfonate groups. This mimics the progressive removal of fatty acids from the hair surface by bleaching and leads to chemically heterogeneous surfaces. Using molecular dynamics (MD) simulations, we study the island structures formed by the lipid monolayers at different degrees of damage in vacuum and in the presence of polar (water) and non-polar (n-hexadecane) solvents. We also use MD simulations to compare the wetting behaviour of water and n-hexadecane droplets on the model surfaces through contact angle measurements, which are compared to experiments using virgin and bleached hair. The model surfaces capture the experimentally-observed transition of the hair surface from hydrophobic (and oleophilic) to hydrophilic (and oleophobic) as the level of bleaching damage increases. By selecting surfaces with specific damage ratios, we obtain contact angles from the MD simulations that are in good agreement with experiments for both solvents on virgin and bleached human hairs. To negate the possible effects of microscale curvature and roughness of real hairs on wetting, we also conduct additional experiments using biomimetic surfaces that are co-functionalised with fatty acids and sulfonate groups. In both the MD simulations and experiments, the cosine of the water contact angle increases linearly with the sulfonate group surface coverage with a similar slope. We expect that the proposed systems will be useful for future molecular dynamics simulations of the adsorption and tribological behaviour of hair, as well as other chemically heterogeneous surfaces.

Structures of the TMC-1 complex illuminate mechanosensory transduction

The initial step in the sensory transduction pathway underpinning hearing and balance in mammals involves the conversion of force into the gating of a mechanosensory transduction channel¹. Despite the profound socioeconomic impacts of hearing disorders and the fundamental biological significance of understanding mechanosensory transduction, the composition, structure and mechanism of the mechanosensory transduction complex have remained poorly characterized. Here we report the single-particle cryo-electron microscopy structure of the native transmembrane channel-like protein 1 (TMC-1) mechanosensory transduction complex isolated from Caenorhabditis elegans. The two-fold symmetric complex is composed of two copies each of the pore-forming TMC-1 subunit, the calcium-binding protein CALM-1 and the transmembrane inner ear protein TMIE. CALM-1 makes extensive contacts with the cytoplasmic face of the TMC-1 subunits, whereas the single-pass TMIE subunits reside on the periphery of the complex, poised like the handles of an accordion. A subset of complexes additionally includes a single arrestin-like protein, arrestin domain protein (ARRD-6), bound to a CALM-1 subunit. Single-particle reconstructions and molecular dynamics simulations show how the mechanosensory transduction complex deforms the membrane bilayer and suggest crucial roles for lipid-protein interactions in the mechanism by which mechanical force is transduced to ion channel gating.

S-acylation controls SARS-CoV-2 membrane lipid organization and enhances infectivity

SARS-CoV-2 virions are surrounded by a lipid bilayer that contains membrane proteins such as spike, responsible for target-cell binding and virus fusion. We found that during SARS-CoV-2 infection, spike becomes lipid modified, through the sequential action of the S-acyltransferases ZDHHC20 and 9. Particularly striking is the rapid acylation of spike on 10 cytosolic cysteines within the ER and Golgi. Using a combination of computational, lipidomics, and biochemical approaches, we show that this massive lipidation controls spike biogenesis and degradation, and drives the formation of localized ordered cholesterol and sphingolipid-rich lipid nanodomains in the early Golgi, where viral budding occurs. Finally, S-acylation of spike allows the formation of viruses with enhanced fusion capacity. Our study points toward S-acylating enzymes and lipid biosynthesis enzymes as novel therapeutic anti-viral targets.

Molecular Modeling of the HR2 and Transmembrane Domains of the SARS-CoV-2 S Protein in the Prefusion State

SARS-CoV-2, the causative agent of COVID-19, remains the focus of research worldwide. SARS-CoV-2 entry into the cell starts with its S protein binding to the angiotensin-converting enzyme-2 (ACE2) expressed on the cell surface. The knowledge of the S protein’s spatial structure is indispensable for understanding the molecular principles of its work. The S protein structure has been almost fully described using experimental approaches with the only exception for the protein’s endodomain, the transmembrane domain, and the ectodomain parts adjacent to the latter. The paper reports molecular modelling of the S protein fragment corresponding to its coiled coil HR2 domain and fully palmitoylated transmembrane domain. Model stability in lipid bilayer was confirmed by all-atom and coarse-grained molecular dynamics simulations. It has been demonstrated that palmitoylation leads to a significant decrease in transmembrane domain mobility and local bilayer thickening, which may be relevant for protein trimerization.

Coarse-Grained Molecular Dynamics Simulations of Membrane Proteins: A Practical Guide

Current computer architectures, coupled with state-of-the-art molecular dynamics simulation software, facilitate the in-depth study of large biomolecular systems at high levels of detail. However, biological phenomena take place at various time and length scales and as a result a multiscale approach must be adopted. One such approach is coarse-graining, where biochemical accuracy is sacrificed for computational efficiency. Here, we present a practical guide to setting up and carrying out coarse-grained molecular dynamics simulations.

Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG

Understanding the molecular mechanisms by which amyloidogenic proteins interact with membranes is a challenging task. Amyloid accumulates from many human diseases have been observed to contain membrane lipids. In this work, coarse-grained molecular dynamics simulations have been used to inspect hen egg white lysozyme (HEWL) aggregation and membrane association in the presence of a pure POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayer and a POPC and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol) mixed bilayer. It was observed that, in both cases, two HEWLs formed aggregates. In the presence of a mixed bilayer, after aggregation, the aggregated system started to interact with the membrane. It has been found that one of the lysozymes which came closer to the mixed bilayer unfolded more. The process of the initial insertion of an aggregated system in the mixed bilayer has been analyzed. The structural rearrangements of the protein and lipids were analyzed as well along the course of the simulation. Although with a pure POPC bilayer, aggregation was observed, the aggregated system moved away from the membrane. We believe that our study will provide considerable insights into lysozyme aggregation in the presence of a membrane environment.

Effect of palmitoylation on the dimer formation of the human dopamine transporter

The human dopamine transporter (hDAT) is one in three members of the monoamine transporter family (MAT). hDAT is essential for regulating the dopamine concentration in the synaptic cleft through dopamine reuptake into the presynaptic neuron; thereby controlling hDAT dopamine signaling. Dysfunction of the transporter is linked to several psychiatric disorders. hDAT and the other MATs have been shown to form oligomers in the plasma membrane, but only limited data exists on which dimeric and higher order oligomeric states are accessible and energetically favorable. In this work, we present several probable dimer conformations using computational coarse-grained self-assembly simulations and assess the relative stability of the different dimer conformations using umbrella sampling replica exchange molecular dynamics. Overall, the dimer conformations primarily involve TM9 and/or TM11 and/or TM12 at the interface. Furthermore, we show that a palmitoyl group (palm) attached to hDAT on TM12 modifies the free energy of separation for interfaces involving TM12, suggesting that S-palmitoylation may change the relative abundance of dimers involving TM12 in a biological context. Finally, a comparison of the identified interfaces of hDAT and palmitoylated hDAT to the human serotonin transporter interfaces and the leucine transporter interface, suggests similar dimer conformations across these protein family.

Molecular Bases of the Membrane Association Mechanism Potentiating Antibiotic Resistance by New Delhi Metallo-β-lactamase 1

Resistance to last-resort carbapenem antibiotics is an increasing threat to human health, as it critically limits therapeutic options. Metallo-β-lactamases (MBLs) are the largest family of carbapenemases, enzymes that inactivate these drugs. Among MBLs, New Delhi metallo-β-lactamase 1 (NDM-1) has experienced the fastest and largest worldwide dissemination. This success has been attributed to the fact that NDM-1 is a lipidated protein anchored to the outer membrane of bacteria, while all other MBLs are soluble periplasmic enzymes. By means of a combined experimental and computational approach, we show that NDM-1 interacts with the surface of bacterial membranes in a stable, defined conformation, in which the active site is not occluded by the bilayer. Although the lipidation is required for a long-lasting interaction, the globular domain of NDM-1 is tuned to interact specifically with the outer bacterial membrane. In contrast, this affinity is not observed for VIM-2, a natively soluble MBL. Finally, we identify key residues involved in the membrane interaction with NDM-1, which constitute potential targets for developing therapeutic strategies able to combat resistance granted by this enzyme.

Capturing Choline-Aromatics Cation-π Interactions in the MARTINI Force Field

Cation−π interactions play an important role in biomolecular recognition, including interactions between membrane phosphatidylcholine lipids and aromatic amino acids of peripheral proteins. While molecular mechanics coarse grain (CG) force fields are particularly well suited to simulate membrane proteins in general, they are not parameterized to explicitly reproduce cation−π interactions. We here propose a modification of the polarizable MARTINI coarse grain (CG) model enabling it to model membrane binding events of peripheral proteins whose aromatic amino acid interactions with choline headgroups are crucial for their membrane binding. For this purpose, we first collected and curated a dataset of eight peripheral proteins from different families. We find that the MARTINI CG model expectedly underestimates aromatics–choline interactions and is unable to reproduce membrane binding of the peripheral proteins in our dataset. Adjustments of the relevant interactions in the polarizable MARTINI force field yield significant improvements in the observed binding events. The orientation of each membrane-bound protein is comparable to reference data from all-atom simulations and experimental binding data. We also use negative controls to ensure that choline–aromatics interactions are not overestimated. We finally check that membrane properties, transmembrane proteins, and membrane translocation potential of mean force (PMF) of aromatic amino acid side-chain analogues are not affected by the new parameter set. This new version “MARTINI 2.3P” is a significant improvement over its predecessors and is suitable for modeling membrane proteins including peripheral membrane binding of peptides and proteins.

Post-Translational Modifications at the Coarse-Grained Level with the SIRAH Force Field

Post-translational modifications (PTMs) on proteins significantly enlarge the physicochemical diversity present in biological macromolecules, altering function, localization, and interactions. Despite their critical role in regulating cellular processes, theoretical methods are not yet fully capable of coping with this diversity. These limitations are particularly more marked for coarse-grained (CG) models, in which comprehensive and self-consistent parametrizations are less frequent. Here we present a set of topologies and interaction parameters for the most common PTMs, fully compatible with the SIRAH force field. The PTMs introduced here reach the same level of structural description of the existing SIRAH force field, expanding the chemical spectrum with promising applications in dynamical protein-protein interactions in large and complex cellular environments.