Lipid unsaturation promotes BAX and BAK pore activity during apoptosis.

BAX and BAK are proapoptotic members of the BCL2 family that directly mediate mitochondrial outer membrane permeabilition (MOMP), a central step in apoptosis execution. However, the molecular architecture of the mitochondrial apoptotic pore remains a key open question and especially little is known about the contribution of lipids to MOMP. By performing a comparative lipidomics analysis of the proximal membrane environment of BAK isolated in lipid nanodiscs, we find a significant enrichment of unsaturated species nearby BAK and BAX in apoptotic conditions. We then demonstrate that unsaturated lipids promote BAX pore activity in model membranes, isolated mitochondria and cellular systems, which is further supported by molecular dynamics simulations. Accordingly, the fatty acid desaturase FADS2 not only enhances apoptosis sensitivity, but also the activation of the cGAS/STING pathway downstream mtDNA release. The correlation of FADS2 levels with the sensitization to apoptosis of different lung and kidney cancer cell lines by co-treatment with unsaturated fatty acids supports the relevance of our findings. Altogether, our work provides an insight on how local lipid environment affects BAX and BAK function during apoptosis.

A sporulation signature protease is required for assembly of the spore surface layers, germination and host colonization in Clostridioides difficile.

A genomic signature for endosporulation includes a gene coding for a protease, YabG, which in the model organism Bacillus subtilis is involved in assembly of the spore coat. We show that in the human pathogen Clostridioidesm difficile, YabG is critical for the assembly of the coat and exosporium layers of spores. YabG is produced during sporulation under the control of the mother cell-specific regulators σE and σK and associates with the spore surface layers. YabG shows an N-terminal SH3-like domain and a C-terminal domain that resembles single domain response regulators, such as CheY, yet is atypical in that the conserved phosphoryl-acceptor residue is absent. Instead, the CheY-like domain carries residues required for activity, including Cys207 and His161, the homologues of which form a catalytic diad in the B. subtilis protein, and also Asp162. The substitution of any of these residues by Ala, eliminates an auto-proteolytic activity as well as interdomain processing of CspBA, a reaction that releases the CspB protease, required for proper spore germination. An in-frame deletion of yabG or an allele coding for an inactive protein, yabGC207A, both cause misassemby of the coat and exosporium and the formation of spores that are more permeable to lysozyme and impaired in germination and host colonization. Furthermore, we show that YabG is required for the expression of at least two σK-dependent genes, cotA, coding for a coat protein, and cdeM, coding for a key determinant of exosporium assembly. Thus, YabG also impinges upon the genetic program of the mother cell possibly by eliminating a transcriptional repressor. Although this activity has not been described for the B. subtilis protein and most of the YabG substrates vary among sporeformers, the general role of the protease in the assembly of the spore surface is likely to be conserved across evolutionary distance.

Xport-A functions as a chaperone by stabilizing the first five transmembrane domains of rhodopsin-1.

Rhodopsin-1 (Rh1), the main photosensitive protein of Drosophila, is a seven-transmembrane domain protein, which is inserted co-translationally in the endoplasmic reticulum (ER) membrane. Biogenesis of Rh1 occurs in the ER, where various chaperones interact with Rh1 to aid in its folding and subsequent transport from the ER to the rhabdomere, the light-sensing organelle of the photoreceptors. Xport-A has been proposed as a chaperone/transport factor for Rh1, but the exact molecular mechanism for Xport-A activity upon Rh1 is unknown. Here, we propose a model where Xport-A functions as a chaperone during the biogenesis of Rh1 in the ER by stabilizing the first five transmembrane domains (TMDs) of Rh1.

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.

Rotation of the c-Ring Promotes the Curvature Sorting of Monomeric ATP Synthases.

ATP synthases are proteins that catalyse the formation of ATP through the rotatory movement of their membrane-spanning subunit. In mitochondria, ATP synthases are found to arrange as dimers at the high-curved edges of cristae. Here, a direct link is explored between the rotatory movement of ATP synthases and their preference for curved membranes. An active curvature sorting of ATP synthases in lipid nanotubes pulled from giant vesicles is found. Coarse-grained simulations confirm the curvature-seeking behaviour of rotating ATP synthases, promoting reversible and frequent protein-protein contacts. The formation of transient protein dimers relies on the membrane-mediated attractive interaction of the order of 1.5 kB T produced by a hydrophobic mismatch upon protein rotation. Transient dimers are sustained by a conic-like arrangement characterized by a wedge angle of θ ≈ 50°, producing a dynamic coupling between protein shape and membrane curvature. The results suggest a new role of the rotational movement of ATP synthases for their dynamic self-assembly in biological membranes.

Activity modulation of the Escherichia coli F1FO ATP synthase by a designed antimicrobial peptide via cardiolipin sequestering.

Most antimicrobial peptides (AMPs) exert their microbicidal activity through membrane permeabilization. The designed AMP EcDBS1R4 has a cryptic mechanism of action involving the membrane hyperpolarization of Escherichia coli, suggesting that EcDBS1R4 may hinder processes involved in membrane potential dissipation. We show that EcDBS1R4 can sequester cardiolipin, a phospholipid that interacts with several respiratory complexes of E. coli. Among these, F1FO ATP synthase uses membrane potential to fuel ATP synthesis. We found that EcDBS1R4 can modulate the activity of ATP synthase upon partition to membranes containing cardiolipin. Molecular dynamics simulations suggest that EcDBS1R4 alters the membrane environment of the transmembrane FO motor, impairing cardiolipin interactions with the cytoplasmic face of the peripheral stalk that binds the catalytic F1 domain to the FO domain. The proposed mechanism of action, targeting membrane protein function through lipid reorganization may open new venues of research on the mode of action and design of other AMPs.

Studies on the Interaction between Model Proteins and Fluorinated Ionic Liquids.

Proteins are inherently unstable, which limits their use as therapeutic agents. However, the use of biocompatible cosolvents or surfactants can help to circumvent this problem through the stabilization of intramolecular and solvent-mediated interactions. Ionic liquids (ILs) have been known to act as cosolvents or surface-active compounds. In the presence of proteins, ILs can have a beneficial effect on their refolding, shelf life, stability, and enzymatic activities. In the work described herein, we used small-angle X-ray scattering (SAXS) to monitor the aggregation of different concentrations of ILs with protein models, lysozyme (Lys) and bovine serum albumin (BSA), and fluorescence microscopy to assess micelle formation of fluorinated ILs (FILs) with Lys. Furthermore, coarse-grained molecular dynamics (CG-MD) simulations provided a better understanding of Lys-FIL interactions. The results showed that the proteins maintain their globular structures in the presence of FILs, with signs of partial unfolding for Lys and compaction for BSA with increased flexibility at higher FIL concentrations. Lys was encapsulated by FIL, thus reinforcing the potential of ILs to be used in the formulation of protein-based pharmaceuticals.

Temperature-Dependent Re-alignment of the Short Multifunctional Peptide BP100 in Membranes Revealed by Solid-State NMR Spectroscopy and Molecular Dynamics Simulations.

BP100 is a cationic undecamer peptide with antimicrobial and cell-penetrating activities. The orientation of this amphiphilic α-helix in lipid bilayers was examined under numerous conditions using solid-state 19 F, 15 N and 2 H NMR. At high temperatures in saturated phosphatidylcholine lipids, BP100 lies flat on the membrane surface, as expected. Upon lowering the temperature towards the lipid phase transition, the helix is found to flip into an upright transmembrane orientation. In thin bilayers, this inserted state was stable at low peptide concentration, but thicker membranes required higher peptide concentrations. In the presence of lysolipids, the inserted state prevailed even at high temperature. Molecular dynamics simulations suggest that BP100 monomer insertion can be stabilized by snorkeling lysine side chains. These results demonstrate that even a very short helix like BP100 can span (and thereby penetrate through) a cellular membrane under suitable conditions.

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.

Parainfluenza Fusion Peptide Promotes Membrane Fusion by Assembling into Oligomeric Porelike Structures.

Paramyxoviruses are enveloped viruses harboring a negative-sense RNA genome that must enter the host's cells to replicate. In the case of the parainfluenza virus, the cell entry process starts with the recognition and attachment to target receptors, followed by proteolytic cleavage of the fusion glycoprotein (F) protein, exposing the fusion peptide (FP) region. The FP is responsible for binding to the target membrane, and it is believed to play a crucial role in the fusion process, but the mechanism by which the parainfluenza FP (PIFP) promotes membrane fusion is still unclear. To elucidate this matter, we performed biophysical experimentation of the PIFP in membranes, together with coarse grain (CG) and atomistic (AA) molecular dynamics (MD) simulations. The simulation results led to the pinpointing of the most important PIFP amino acid residues for membrane fusion and show that, at high concentrations, the peptide induces the formation of a water-permeable porelike structure. This structure promotes lipid head intrusion and lipid tail protrusion, which facilitates membrane fusion. Biophysical experimental results validate these findings, showing that, depending on the peptide/lipid ratio, the PIFP can promote fusion and/or membrane leakage. Our work furthers the understanding of the PIFP-induced membrane fusion process, which might help foster development in the field of viral entry inhibition.

Coarse-Grain Simulations of Membrane-Adsorbed Helical Peptides.

The amphipathic α-helix is a common motif for peptide adsorption to membranes. Many physiologically relevant events involving membrane-adsorbed peptides occur over time and size scales readily accessible to coarse-grain molecular dynamics simulations. This methodological suitability, however, comes with a number of pitfalls. Here, I exemplify a multi-step adsorption equilibration procedure on the antimicrobial peptide Magainin 2. It involves careful control of peptide freedom to promote optimal membrane adsorption before other interactions are allowed. This shortens preparation times prior to production simulations while avoiding divergence into unrealistic or artifactual configurations.

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.

Overlapping Properties of the Short Membrane-Active Peptide BP100 With (i) Polycationic TAT and (ii) α-helical Magainin Family Peptides.

BP100 is a short, designer-made membrane-active peptide with multiple functionalities: antimicrobial, cell-penetrating, and fusogenic. Consisting of five lysines and 6 hydrophobic residues, BP100 was shown to bind to lipid bilayers as an amphipathic α-helix, but its mechanism of action remains unclear. With these features, BP100 embodies the characteristics of two distinctly different classes of membrane-active peptides, which have been studied in detail and where the mechanism of action is better understood. On the one hand, its amphiphilic helical structure is similar to the pore forming magainin family of antimicrobial peptides, though BP100 is much too short to span the membrane. On the other hand, its length and high charge density are reminiscent of the HIV-TAT family of cell penetrating peptides, for which inverted micelles have been postulated as translocation intermediates, amongst other mechanisms. Assays were performed to test the antimicrobial and hemolytic activity, the induced leakage and fusion of lipid vesicles, and cell uptake. From these results the functional profiles of BP100, HIV-TAT, and the magainin-like peptides magainin 2, PGLa, MSI-103, and MAP were determined and compared. It is observed that the activity of BP100 resembles most closely the much longer amphipathic α-helical magainin-like peptides, with high antimicrobial activity along with considerable fusogenic and hemolytic effects. In contrast, HIV-TAT shows almost no antimicrobial, fusogenic, or hemolytic effects. We conclude that the amphipathic helix of BP100 has a similar membrane-based activity as magainin-like peptides and may have a similar mechanism of action.

Coarse-Grained Parameterization of Nucleotide Cofactors and Metabolites: Protonation Constants, Partition Coefficients, and Model Topologies.

Nucleotides are structural units relevant not only in nucleic acids but also as substrates or cofactors in key biochemical reactions. The size- and timescales of such nucleotide-protein interactions fall well within the scope of coarse-grained molecular dynamics, which holds promise of important mechanistic insight. However, the lack of specific parameters has prevented accurate coarse-grained simulations of protein interactions with most nucleotide compounds. In this work, we comprehensively develop coarse-grained parameters for key metabolites/cofactors (FAD, FMN, riboflavin, NAD, NADP, ATP, ADP, AMP, and thiamine pyrophosphate) in different oxidation and protonation states as well as for smaller molecules derived from them (among others, nicotinamide, adenosine, adenine, ribose, thiamine, and lumiflavin), summing up a total of 79 different molecules. In line with the Martini parameterization methodology, parameters were tuned to reproduce octanol-water partition coefficients. Given the lack of existing data, we set out to experimentally determine these partition coefficients, developing two methodological approaches, based on 31P-NMR and fluorescence spectroscopy, specifically tailored to the strong hydrophilicity of most of the parameterized compounds. To distinguish the partition of each relevant protonation species, we further potentiometrically characterized the protonation constants of key molecules. This work successfully builds a comprehensive and relevant set of computational models that will boost the biochemical application of coarse-grained simulations. It does so based on the measurement of partition and acid-base physicochemical data that, in turn, covers important gaps in nucleotide characterization.

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.

Localization Preference of Antimicrobial Peptides on Liquid-Disordered Membrane Domains.

We performed coarse-grained simulations of the antimicrobial peptides Magainin-2, BP100, MSI-103, and MSI-78 on a phase-separated membrane to study their preference for the different domains. All the peptides displayed a clear preference for the liquid-disordered (Ld) phase over the liquid-ordered (Lo) one. For BP100, MSI-103, and MSI-78 there was a further preference of the peptides for the domain interface. The peptides' preference toward the disordered phase was shown to reflect a penalization of lipid-lipid interaction enthalpy in the Lo phase, when in the vicinity of peptides. Similar results were observed at the two studied concentrations, although Ld phase saturation at the higher concentration drove some of the peptide excess to the Lo phase. Magainin-2 and MSI-103 were found to dimerize, in agreement with available experimental data. Interestingly, at high concentrations of Magainin-2 toroidal pores spontaneously formed in the Ld phase. We performed additional simulations to characterize this phenomenon, which is likely related to Magainin-2's membranolytic action.

Charge-dependent interactions of monomeric and filamentous actin with lipid bilayers.

The cytoskeletal protein actin polymerizes into filaments that are essential for the mechanical stability of mammalian cells. In vitro experiments showed that direct interactions between actin filaments and lipid bilayers are possible and that the net charge of the bilayer as well as the presence of divalent ions in the buffer play an important role. In vivo, colocalization of actin filaments and divalent ions are suppressed, and cells rely on linker proteins to connect the plasma membrane to the actin network. Little is known, however, about why this is the case and what microscopic interactions are important. A deeper understanding is highly beneficial, first, to obtain understanding in the biological design of cells and, second, as a possible basis for the building of artificial cortices for the stabilization of synthetic cells. Here, we report the results of coarse-grained molecular dynamics simulations of monomeric and filamentous actin in the vicinity of differently charged lipid bilayers. We observe that charges on the lipid head groups strongly determine the ability of actin to adsorb to the bilayer. The inclusion of divalent ions leads to a reversal of the binding affinity. Our in silico results are validated experimentally by reconstitution assays with actin on lipid bilayer membranes and provide a molecular-level understanding of the actin-membrane interaction.

Pitfalls of the Martini Model.

The computational and conceptual simplifications realized by coarse-grain (CG) models make them a ubiquitous tool in the current computational modeling landscape. Building block based CG models, such as the Martini model, possess the key advantage of allowing for a broad range of applications without the need to reparametrize the force field each time. However, there are certain inherent limitations to this approach, which we investigate in detail in this work. We first study the consequences of the absence of specific cross Lennard-Jones parameters between different particle sizes. We show that this lack may lead to artificially high free energy barriers in dimerization profiles. We then look at the effect of deviating too far from the standard bonded parameters, both in terms of solute partitioning behavior and solvent properties. Moreover, we show that too weak bonded force constants entail the risk of artificially inducing clustering, which has to be taken into account when designing elastic network models for proteins. These results have implications for the current use of the Martini CG model and provide clear directions for the reparametrization of the Martini model. Moreover, our findings are generally relevant for the parametrization of any other building block based force field.

Ceramides bind VDAC2 to trigger mitochondrial apoptosis.

Ceramides draw wide attention as tumor suppressor lipids that act directly on mitochondria to trigger apoptotic cell death. However, molecular details of the underlying mechanism are largely unknown. Using a photoactivatable ceramide probe, we here identify the voltage-dependent anion channels VDAC1 and VDAC2 as mitochondrial ceramide binding proteins. Coarse-grain molecular dynamics simulations reveal that both channels harbor a ceramide binding site on one side of the barrel wall. This site includes a membrane-buried glutamate that mediates direct contact with the ceramide head group. Substitution or chemical modification of this residue abolishes photolabeling of both channels with the ceramide probe. Unlike VDAC1 removal, loss of VDAC2 or replacing its membrane-facing glutamate with glutamine renders human colon cancer cells largely resistant to ceramide-induced apoptosis. Collectively, our data support a role of VDAC2 as direct effector of ceramide-mediated cell death, providing a molecular framework for how ceramides exert their anti-neoplastic activity.

Structure-Stability-Function Mechanistic Links in the Anti-Measles Virus Action of Tocopherol-Derivatized Peptide Nanoparticles.

Measles remains one of the leading causes of child mortality worldwide and is re-emerging in some countries due to poor vaccine coverage, concomitant with importation of measles virus (MV) from endemic areas. The lack of specific chemotherapy contributes to negative outcomes, especially in infants or immunodeficient individuals. Fusion inhibitor peptides derived from the MV Fusion protein C-terminal Heptad Repeat (HRC) targeting MV envelope fusion glycoproteins block infection at the stage of entry into host cells, thus preventing viral multiplication. To improve efficacy of such entry inhibitors, we have modified a HRC peptide inhibitor by introducing properties of self-assembly into nanoparticles (NP) and higher affinity for both viral and cell membranes. Modification of the peptide consisted of covalent grafting with tocopherol to increase amphipathicity and lipophilicity (HRC5). One additional peptide inhibitor consisting of a peptide dimer grafted to tocopherol was also used (HRC6). Spectroscopic, imaging, and simulation techniques were used to characterize the NP and explore the molecular basis for their antiviral efficacy. HRC5 forms micellar stable NP while HRC6 aggregates into amorphous, loose, unstable NP. Interpeptide cluster bridging governs NP assembly into dynamic metastable states. The results are consistent with the conclusion that the improved efficacy of HRC6 relative to HRC5 can be attributed to NP instability, which leads to more extensive partition to target membranes and binding to viral target proteins.