Bartender: Martini 3 Bonded Terms via Quantum Mechanics-Based Molecular Dynamics.

Coarse-grained (CG) molecular dynamics (MD) simulations have grown in applicability over the years. The recently released version of the Martini CG force field (Martini 3) has been successfully applied to simulate many processes, including protein-ligand binding. However, the current ligand parametrization scheme is manual and requires an a priori reference all-atom (AA) simulation for benchmarking. For systems with suboptimal AA parameters, which are often unknown, this translates into a CG model that does not reproduce the true dynamical behavior of the underlying molecule. Here, we present Bartender, a quantum mechanics (QM)/MD-based parametrization tool written in Go. Bartender harnesses the power of QM simulations and produces reasonable bonded terms for Martini 3 CG models of small molecules in an efficient and user-friendly manner. For small, ring-like molecules, Bartender generates models whose properties are indistinguishable from the human-made models. For more complex, drug-like ligands, it is able to fit functional forms beyond simple harmonic dihedrals and thus better captures their dynamical behavior. Bartender has the power to both increase the efficiency and the accuracy of Martini 3-based high-throughput applications by producing numerically stable and physically realistic CG models.

Assessing the Martini 3 protein model: A review of its path and potential.

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

Martini 3 Coarse-Grained Force Field for Cholesterol.

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

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.

Expulsion mechanism of the substrate-translocating subunit in ECF transporters.

Energy-coupling factor (ECF)-type transporters mediate the uptake of micronutrients in many bacteria. They consist of a substrate-translocating subunit (S-component) and an ATP-hydrolysing motor (ECF module) Previous data indicate that the S-component topples within the membrane to alternately expose the binding site to either side of the membrane. In many ECF transporters, the substrate-free S-component can be expelled from the ECF module. Here we study this enigmatic expulsion step by cryogenic electron microscopy and reveal that ATP induces a concave-to-convex shape change of two long helices in the motor, thereby destroying the S-component's docking site and allowing for its dissociation. We show that adaptation of the membrane morphology to the conformational state of the motor may favour expulsion of the substrate-free S-component when ATP is bound and docking of the substrate-loaded S-component after hydrolysis. Our work provides a picture of bilayer-assisted chemo-mechanical coupling in the transport cycle of ECF transporters.

SARS-CoV-2 variants impact RBD conformational dynamics and ACE2 accessibility.

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has killed over 6 million people and is having a devastating social and economic impact around the world. The rise of new variants of concern (VOCs) represents a difficult challenge due to the loss of vaccine and natural immunity, as well as increased transmissibility. All VOCs contain mutations in the spike glycoprotein, which mediates fusion between the viral and host cell membranes. The spike glycoprotein binds to angiotensin-converting enzyme 2 (ACE2) via its receptor binding domain (RBD) initiating the infection process. Attempting to understand the effect of RBD mutations in VOCs, a lot of attention has been given to the RBD-ACE2 interaction. However, this type of analysis ignores more indirect effects, such as the conformational dynamics of the RBD itself. Observing that some mutations occur in residues that are not in direct contact with ACE2, we hypothesized that they could affect the RBD conformational dynamics. To test this, we performed long atomistic (AA) molecular dynamics (MD) simulations to investigate the structural dynamics of wt RBD, and that of four VOCs (Alpha, Beta, Delta, and Omicron). Our results show that the wt RBD presents two distinct conformations: an "open" conformation where it is free to bind ACE2; and a "closed" conformation, where the RBM ridge blocks the binding surface. The Alpha and Beta variants shift the open/closed equilibrium towards the open conformation by roughly 20%, likely increasing ACE2 binding affinity. Simulations of the Delta and Omicron variants showed extreme results, with the closed conformation being rarely observed. The Delta variant also differed substantially from the other variants, alternating between the open conformation and an alternative "reversed" one, with a significantly changed orientation of the RBM ridge. This alternate conformation could provide a fitness advantage due to increased availability for ACE2 binding, and by aiding antibody escape through epitope occlusion. These results support the hypothesis that VOCs, and particularly the Omicron and Delta variants, impact RBD conformational dynamics in a direction that promotes efficient binding to ACE2 and, in the case of Delta, may assist antibody escape.

Impact of Ca2+-Induced PI(4,5)P2 Clusters on PH-YFP Organization and Protein-Protein Interactions.

Despite its low abundance, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) is a key modulator of membrane-associated signaling events in eukaryotic cells. Temporal and spatial regulation of PI(4,5)P2 concentration can achieve localized increases in the levels of this lipid, which are crucial for the activation or recruitment of peripheral proteins to the plasma membrane. The recent observation of the dramatic impact of physiological divalent cation concentrations on PI(4,5)P2 clustering, suggests that protein anchoring to the plasma membrane through PI(4,5)P2 is likely not defined solely by a simple (monomeric PI(4,5)P2)/(protein bound PI(4,5)P2) equilibrium, but instead depends on complex protein interactions with PI(4,5)P2 clusters. The insertion of PI(4,5)P2-binding proteins within these clusters can putatively modulate protein-protein interactions in the membrane, but the relevance of such effects is largely unknown. In this work, we characterized the impact of Ca2+ on the organization and protein-protein interactions of PI(4,5)P2-binding proteins. We show that, in giant unilamellar vesicles presenting PI(4,5)P2, the membrane diffusion properties of pleckstrin homology (PH) domains tagged with a yellow fluorescent protein (YFP) are affected by the presence of Ca2+, suggesting direct interactions between the protein and PI(4,5)P2 clusters. Importantly, PH-YFP is found to dimerize in the membrane in the absence of Ca2+. This oligomerization is inhibited in the presence of physiological concentrations of the divalent cation. These results confirm that cation-dependent PI(4,5)P2 clustering promotes interactions between PI(4,5)P2-binding proteins and has the potential to dramatically influence the organization and downstream interactions of PI(4,5)P2-binding proteins in the plasma membrane.

Improved Parameterization of Phosphatidylinositide Lipid Headgroups for the Martini 3 Coarse-Grain Force Field.

Phosphoinositides are a family of membrane phospholipids that play crucial roles in membrane regulatory events. As such, these lipids are often a key part of molecular dynamics simulation studies of biological membranes, in particular of those employing coarse-grain models because of the potential long times and sizes of the involved membrane processes. Version 3 of the widely used Martini coarse-grain force field has been recently published, greatly refining many aspects of biomolecular interactions. In order to properly use it for lipid membrane simulations with phosphoinositides, we put forth the Martini 3-specific parameterization of inositol, phosphatidylinositol, and seven physiologically relevant phosphorylated derivatives of phosphatidylinositol. Compared to parameterizations for earlier Martini versions, focus was put on a more accurate reproduction of the behavior seen in both atomistic simulations and experimental studies, including the signaling-relevant phosphoinositide interaction with divalent cations. The models that we develop improve upon the conformational dynamics of phosphoinositides in the Martini force field and provide stable topologies at typical Martini time steps. They are able to reproduce experimentally known protein-binding poses as well as phosphoinositide aggregation tendencies. The latter was tested both in the presence and absence of calcium and included correct behavior of PI(4,5)P2 calcium-induced clusters, which can be of relevance for regulation.

Quantitative FRET Microscopy Reveals a Crucial Role of Cytoskeleton in Promoting PI(4,5)P2 Confinement.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) is an essential plasma membrane component involved in several cellular functions, including membrane trafficking and cytoskeleton organization. This function multiplicity is partially achieved through a dynamic spatiotemporal organization of PI(4,5)P2 within the membrane. Here, we use a Förster resonance energy transfer (FRET) approach to quantitatively assess the extent of PI(4,5)P2 confinement within the plasma membrane. This methodology relies on the rigorous evaluation of the dependence of absolute FRET efficiencies between pleckstrin homology domains (PHPLCδ) fused with fluorescent proteins and their average fluorescence intensity at the membrane. PI(4,5)P2 is found to be significantly compartmentalized at the plasma membrane of HeLa cells, and these clusters are not cholesterol-dependent, suggesting that membrane rafts are not involved in the formation of these nanodomains. On the other hand, upon inhibition of actin polymerization, compartmentalization of PI(4,5)P2 is almost entirely eliminated, showing that the cytoskeleton network is the critical component responsible for the formation of nanoscale PI(4,5)P2 domains in HeLa cells.

Acyl-chain saturation regulates the order of phosphatidylinositol 4,5-bisphosphate nanodomains.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) plays a critical role in the regulation of various plasma membrane processes and signaling pathways in eukaryotes. A significant amount of cellular resources are spent on maintaining the dominant 1-stearoyl-2-arachidonyl PI(4,5)P2 acyl-chain composition, while less abundant and more saturated species become more prevalent in response to specific stimuli, stress or aging. Here, we report the impact of acyl-chain structure on the biophysical properties of cation-induced PI(4,5)P2 nanodomains. PI(4,5)P2 species with increasing levels of acyl-chain saturation cluster in progressively more ordered nanodomains, culminating in the formation of gel-like nanodomains for fully saturated species. The formation of these gel-like domains was largely abrogated in the presence of 1-stearoyl-2-arachidonyl PI(4,5)P2. This is, to the best of our knowledge, the first report of the impact of PI(4,5)P2 acyl-chain composition on cation-dependent nanodomain ordering, and provides important clues to the motives behind the enrichment of PI(4,5)P2 with polyunsaturated acyl-chains. We also show how Ca2+-induced PI(4,5)P2 nanodomains are able to generate local negative curvature, a phenomenon likely to play a role in membrane remodeling events.

Structure and Lateral Organization of Phosphatidylinositol 4,5-bisphosphate.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) is a minor but ubiquitous component of the inner leaflet of the plasma membrane of eukaryotic cells. However, due to its particular complex biophysical properties, it stands out from its neighboring lipids as one of the most important regulators of membrane-associated signaling events. Despite its very low steady-state concentration, PI(4,5)P2 is able to engage in a multitude of simultaneous cellular functions that are temporally and spatially regulated through the presence of localized transient pools of PI(4,5)P2 in the membrane. These pools are crucial for the recruitment, activation, and organization of signaling proteins and consequent regulation of downstream signaling. The present review showcases some of the most important PI(4,5)P2 molecular and biophysical properties as well as their impact on its membrane dynamics, lateral organization, and interactions with other biochemical partners.