CHARMM-GUI Martini Maker for modeling and simulation of complex bacterial membranes with lipopolysaccharides

Chain-Length Dependence of Peptide-Lipid Bilayer Interaction Strength and Binding Kinetics: A Combined Theoretical and Experimental Approach

Physical interactions between polypeptide chains and lipid membranes underlie critical cellular processes. Yet, despite fundamental importance, key mechanistic aspects of these interactions remain elusive. Bulk experiments have revealed a linear relationship between free energy and peptide chain length in a model system, but does this linearity extend to the interaction strength and to the kinetics of lipid binding? To address these questions, we utilized a combination of coarse-grained molecular dynamics (CG MD) simulations, analytical modeling, and atomic force microscopy (AFM)-based single molecule force spectroscopy. Following previous bulk experiments, we focused on interactions between short hydrophobic peptides (WLn, n = 1, ..., 5) with 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayers, a simple system that probes peptide primary structure effects. Potentials of mean force extracted from CG MD recapitulated the linearity of free energy with the chain length. Simulation results were quantitatively connected to bulk biochemical experiments via a single scaling factor of order unity, corroborating the methodology. Additionally, CG MD revealed an increase in the distance to the transition state, a result that weakens the dependence of the dissociation force on the peptide chain length. AFM experiments elucidated rupture force distributions and, through modeling, intrinsic dissociation rates. Taken together, the analysis indicates a rupture force plateau in the WLn-POPC system, suggesting that the final rupture event involves the last 2 or 3 residues. In contrast, the linear dependence on chain length was preserved in the intrinsic dissociation rate. This study advances the understanding of peptide-lipid interactions and provides potentially useful insights for the design of peptides with tailored membrane-interacting properties.

Faster but Not Sweeter: A Model of Escherichia coli Re-level Lipopolysaccharide for Martini 3 and a Martini 2 Version with Accelerated Kinetics

Lipopolysaccharide (LPS) is a complex glycolipid molecule that is the main lipidic component of the outer leaflet of the outer membrane of Gram-negative bacteria. It has very limited lateral motion compared to phospholipids, which are more ubiquitous in biological membranes, including in the inner leaflet of the outer membrane of Gram-negative bacteria. The slow-moving nature of LPS can present a hurdle for molecular dynamics simulations, given that the (pragmatically) accessible timescales to simulations are currently limited to microseconds, during which LPS displays some conformational dynamics but hardly any lateral diffusion. Thus, it is not feasible to observe phenomena such as insertion of molecules, including antibiotics/antimicrobials, directly into the outer membrane from the extracellular side nor to observe LPS dissociating from proteins via molecular dynamics using currently available models at the atomistic and more coarse-grained levels of granularity. Here, we present a model of deep rough LPS compatible with the Martini 2 coarse-grained force field with scaled down nonbonded interactions to enable faster diffusion. We show that the faster-diffusing LPS model is able to reproduce the salient biophysical properties of the standard models, but due to its faster lateral motion, molecules are able to penetrate deeper into membranes containing the faster model. We show that the fast ReLPS model is able to reproduce experimentally determined patterns of interaction with outer membrane proteins while also allowing for LPS to associate and dissociate with proteins within microsecond timescales. We also complete the Martini 3 LPS toolkit for Escherichia coli by presenting a (standard) model of deep rough LPS for this force field.

Lipid Organization by the Caveolin-1 Complex

Caveolins are lipid-binding proteins that can organize membrane remodeling and oligomerize into the 8S-complex. The CAV1 8S-complex comprises a disk-like structure, about 15nm in diameter, with a central beta barrel. Further oligomerization of 8S-complexes remodels the membrane into caveolae vessels, with a dependence on cholesterol concentration. However, the molecular mechanisms behind membrane remodeling and cholesterol filtering are still not understood. Performing atomistic Molecular Dynamics simulations in combination with advanced sampling techniques, we describe how the CAV1-8S complex bends the membrane and accumulates cholesterol. Here, our simulations show an enhancing effect by the palmitoylations of CAV1, and we predict that the CAV1-8S complex can extract cholesterol molecules from the lipid bilayer and accommodate them in its beta barrel. Through backmapping to the all-atom level we also conclude that the Martini v2 coarse-grained forcefield overestimates membrane bending, as the atomistic simulations exhibit only very localized bending.

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Molecular Dynamics Simulations Help Determine the Molecular Mechanisms of Lasioglossin-III and Its Variant Peptides' Membrane Interfacial Interactions

Lasioglossin-III (LL-III) is a potent broad-spectrum antimicrobial peptide used in diverse antimicrobial applications. In this work, coarse-grained and all-atom molecular dynamics simulation strategies were used in tandem to interpret the molecular mechanisms involved in the interfacial dynamics of LL-III and its recombinant variants during interactions with diverse cell membrane systems. Our results indicate that the membrane charges act as the driving force for initiating the membrane-peptide interactions, while the hydrophobic or van der Waals forces help to reinforce the membrane-peptide bindings. The optimized charge-hydrophobicity ratio of the LL-III peptides helps ensure their high specificity toward bacterial membranes compared to mammalian membrane systems, which also helps explain our experimental observations. Overall, we hope that our work gives new insight into the antimicrobial action of LL-III peptides and that the adopted simulation strategy will help other scientists and engineers extract maximal information from complex molecular simulations using minimal computational power.

There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

Manufacturing process of liposomal Formation: A coarse-grained molecular dynamics simulation

A method of producing liposomes has been previously developed using a continuous manufacturing technology that involves a co-axial turbulent jet in co-flow. In this study, coarse-grained molecular dynamics (CG-MD) simulations were used to gain a deeper understanding of how the self-assembly process of liposomes is affected by the material attributes (such as the concentration of ethanol) and the process parameters (such as temperature), while also providing detailed information on a nano-scale molecular level. Specifically, the CG-MD simulations yield a comprehensive internal view of the structure and formation mechanisms of liposomes containing DPPC, DPPG, and cholesterol molecules. The importance of this work is that structural details on the molecular level are proposed, and such detail is not possible to obtain through experimental studies alone. The assessment of structural properties, including the area per lipid, diffusion coefficient, and order parameters, indicated that a thicker bilayer was observed at higher ethanol concentrations, while a thinner bilayer was present at higher temperatures. These conditions led to more water penetrating the interior of the bilayer and an unstable structure, as indicated by a larger contact area between lipids and water, and a higher coefficient of lipid lateral diffusion. However, stable liposomes were found through these evaluations at lower ethanol concentrations and/or lower process temperatures. Furthermore, the CG-MD model was further compared and validated with experimental and computational data including liposomal bilayer thickness and area per lipid measurements.

Inhibition and transport mechanisms of the ABC transporter hMRP5

Human multidrug resistance protein 5 (hMRP5) effluxes anticancer and antivirus drugs, driving multidrug resistance. To uncover the mechanism of hMRP5, we determine six distinct cryo-EM structures, revealing an autoinhibitory N-terminal peptide that must dissociate to permit subsequent substrate recruitment. Guided by these molecular insights, we design an inhibitory peptide that could block substrate entry into the transport pathway. We also identify a regulatory motif, comprising a positively charged cluster and hydrophobic patches, within the first nucleotide-binding domain that modulates hMRP5 localization by engaging with membranes. By integrating our structural, biochemical, computational, and cell biological findings, we propose a model for hMRP5 conformational cycling and localization. Overall, this work provides mechanistic understanding of hMRP5 function, while informing future selective hMRP5 inhibitor development. More broadly, this study advances our understanding of the structural dynamics and inhibition of ABC transporters.

Construction of Multiscale Dissipative Particle Dynamics (DPD) Models from Other Coarse-Grained Models

We present a general scheme for converting coarse-grained models into Dissipative Particle Dynamics (DPD) models. We build the corresponding DPD models by analogy with the de novo DPD coarse-graining scheme suggested by Groot and Warren (J. Chem. Phys., 1997). Electrostatic interactions between charged DPD particles are represented though the addition of a long-range Slater Coulomb potential as suggested by González-Melchor et al. (J. Chem. Phys., 2006). The construction is illustrated by converting MARTINI models for various proteins into a DPD representation, but it not restricted to the usual potential form in the MARTINI model-viz., Lennard-Jones potentials. We further extended the DPD scheme away from the typical use of homogeneous particle sizes, therefore faithfully representing the variations in the particle sizes seen in the underlying MARTINI model. The accuracy of the resulting construction of our generalized DPD models with respect to several structural observables has been benchmarked favorably against all-atom and MARTINI models for a selected set of peptides and proteins, and variations in the scales of the coarse-graining of the water solvent.

Building complex membranes with Martini 3

The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.

Dissecting the mechanism of atlastin-mediated homotypic membrane fusion at the single-molecule level

Homotypic membrane fusion of the endoplasmic reticulum (ER) is mediated by dynamin-like GTPase atlastin (ATL). This fundamental process relies on GTP-dependent domain rearrangements in the N-terminal region of ATL (ATLcyto), including the GTPase domain and three-helix bundle (3HB). However, its conformational dynamics during the GTPase cycle remain elusive. Here, we combine single-molecule FRET imaging and molecular dynamics simulations to address this conundrum. Different from the prevailing model, ATLcyto can form a loose crossover dimer upon GTP binding, which is tightened by GTP hydrolysis for membrane fusion. Furthermore, the α-helical motif between the 3HB and transmembrane domain, which is embedded in the surface of the lipid bilayer and self-associates in the crossover dimer, is required for ATL function. To recycle the proteins, Pi release, which disassembles the dimer, activates frequent relative movements between the GTPase domain and 3HB, and subsequent GDP dissociation alters the conformational preference of the ATLcyto monomer for entering the next reaction cycle. Finally, we found that two disease-causing mutations affect human ATL1 activity by destabilizing GTP binding-induced loose crossover dimer formation and the membrane-embedded helix, respectively. These results provide insights into ATL-mediated homotypic membrane fusion and the pathological mechanisms of related disease.

Cooperative Gating of a K+ Channel by Unmodified Biological Anionic Lipids Viewed by Solid-State NMR Spectroscopy

Lipids adhere to membrane proteins to stimulate or suppress molecular and ionic transport and signal transduction. Yet, the molecular details of lipid-protein interaction and their functional impact are poorly characterized. Here we combine NMR, coarse-grained molecular dynamics (CGMD), and functional assays to reveal classic cooperativity in the binding and subsequent activation of a bacterial inward rectifier potassium (Kir) channel by phosphatidylglycerol (PG), a common component of many membranes. Past studies of lipid activation of Kir channels focused primarily on phosphatidylinositol bisphosphate, a relatively rare signaling lipid that is tightly regulated in space and time. We use solid-state NMR to quantify the binding of unmodified 13C-PG to the K+ channel KirBac1.1 in liposomes. This specific lipid-protein interaction has a dissociation constant (Kd) of ∼7 mol percentage PG (ΧPG) with positive cooperativity (n = 3.8) and approaches saturation near 20% ΧPG. Liposomal flux assays show that K+ flux also increases with PG in a cooperative manner with an EC50 of ∼20% ΧPG, within the physiological range. Further quantitative fitting of these data reveals that PG acts as a partial (80%) agonist with fivefold K+ flux amplification. Comparisons of NMR chemical shift perturbation and CGMD simulations at different ΧPG confirm the direct interaction of PG with key residues, several of which would not be accessible to lipid headgroups in the closed state of the channel. Allosteric regulation by a common lipid is directly relevant to the activation mechanisms of several human ion channels. This study highlights the role of concentration-dependent lipid-protein interactions and tightly controlled protein allostery in the activation and regulation of ion channels.

Bilayer tension-induced clustering of the UPR sensor IRE1

The endoplasmic reticulum acts as a protein quality control center where a range of chaperones and foldases facilitates protein folding. IRE1 is a sensory transmembrane protein that transduces signals of proteotoxic stress by forming clusters and activating a cellular program called the unfolded protein response (UPR). Recently, membrane thickness variation due to membrane compositional changes have been shown to drive IRE1 cluster formation, activating the UPR even in the absence of proteotoxic stress. Here, we demonstrate a direct relationship between bilayer tension and UPR activation based on IRE1 dimer stability. The stability of the IRE1 dimer in a (50%DOPC-50%POPC) membrane at different applied bilayer tensions was analyzed via molecular dynamics simulations. The potential of mean force for IRE1 dimerization predicts a higher concentration of IRE1 dimers for both tensed and compressed ER membranes. This study shows that IRE1 may be a mechanosensitive membrane protein and establishes a direct biophysical relationship between bilayer tension and UPR activation.

Molecular Modeling and Simulation of the Mycobacterial Cell Envelope: From Individual Components to Cell Envelope Assemblies

Unlike typical Gram-positive bacteria, the cell envelope of mycobacteria is unique and composed of a mycobacterial outer membrane, also known as the mycomembrane, a peptidoglycan layer, and a mycobacterial inner membrane, which is analogous to that of Gram-negative bacteria. Despite its importance, however, our understanding of this complex cell envelope is rudimentary at best. Thus, molecular modeling and simulation of such an envelope can benefit the scientific community by proposing new hypotheses about the biophysical properties of its different layers. In this Perspective, we present recent advances in molecular modeling and simulation of the mycobacterial cell envelope from individual components to cell envelope assemblies. We also show how modeling other types of cell envelopes, such as that of Escherichia coli, may help modeling part of the mycobacterial envelopes. We hope that the studies presented here are just the beginning of the road and more and more new modeling and simulation studies help us to understand crucial questions related to mycobacteria such as antibiotic resistance or bacterial survival in the host.

Photodynamic inactivation of bacteria: Why it is not enough to excite a photosensitizer

BACKGROUND

The development of multidrug resistance (MDR) in infectious agents is one of the most serious global problems facing humanity. Antimicrobial photodynamic therapy (APDT) shows encouraging results in the fight against MDR pathogens, including those in biofilms.

METHODS

Photosensitizers (PS), monocationic methylene blue, polycationic and polyanionic derivatives of phthalocyanines, electroneutral and polycationic derivatives of bacteriochlorin were used to study photodynamic inactivation of Gram-positive and Gram-negative planktonic bacteria and biofilms under LED irradiation. Zeta potential measurements, confocal fluorescence imaging, and coarse-grained modeling were used to evaluate the interactions of PS with bacteria. PS aggregation and photobleaching were studied using absorption and fluorescence spectroscopy.

RESULTS

The main approaches to ensure high efficiency of bacteria photosensitization are analyzed.

CONCLUSIONS

PS must maintain a delicate balance between binding to exocellular and external structures of bacterial cells and penetration through the cell wall so as not to get stuck on the way to photooxidation-sensitive structures of the bacterial cell.

Computational microbiology of bacteria: Advancements in molecular dynamics simulations

Microbiology is traditionally considered within the context of wet laboratory methodologies. Computational techniques have a great potential to contribute to microbiology. Here, we describe our loose definition of "computational microbiology" and provide a short survey focused on molecular dynamics simulations of bacterial systems that fall within this definition. It is our contention that increased compositional complexity and realistic levels of molecular crowding within simulated systems are key for bridging the divide between experimental and computational microbiology.

Probing the Structural Topology and Dynamic Properties of gp28 Using Continuous Wave Electron Paramagnetic Resonance Spectroscopy

Lysis of Gram-negative bacteria by dsDNA phages is accomplished through either the canonical holin-endolysin pathway or the pinholin-SAR endolysin pathway. During lysis, the outer membrane (OM) is disrupted, typically by two-component spanins or unimolecular spanins. However, in the absence of spanins, phages use alternative proteins called Disruptin to disrupt the OM. The Disruptin family includes the cationic antimicrobial peptide gp28, which is found in the virulent podophage φKT. In this study, EPR spectroscopy was used to analyze the dynamics and topology of gp28 incorporated into a lipid bilayer, revealing differences in mobility, depth parameter, and membrane interaction among different segments and residues of the protein. Our results indicate that multiple points of helix 2 and helix 3 interact with the phospholipid membrane, while others are solvent-exposed, suggesting that gp28 is a surface-bound peptide. The CW-EPR power saturation data and helical wheel analysis confirmed the amphipathic-helical structure of gp28. Additionally, course-grain molecular dynamics simulations were further used to develop the structural model of the gp28 peptide associated with the lipid bilayers. Based on the data obtained in this study, we propose a structural topology model for gp28 with respect to the membrane. This work provides important insights into the structural and dynamic properties of gp28 incorporated into a lipid bilayer environment.

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.

Structural basis of dimerization of chemokine receptors CCR5 and CXCR4

G protein-coupled receptors (GPCRs) are prominent drug targets responsible for extracellular-to-intracellular signal transduction. GPCRs can form functional dimers that have been poorly characterized so far. Here, we show the dimerization mechanism of the chemokine receptors CCR5 and CXCR4 by means of an advanced free-energy technique named coarse-grained metadynamics. Our results reproduce binding events between the GPCRs occurring in the minute timescale, revealing a symmetric and an asymmetric dimeric structure for each of the three investigated systems, CCR5/CCR5, CXCR4/CXCR4, and CCR5/CXCR4. The transmembrane helices TM4-TM5 and TM6-TM7 are the preferred binding interfaces for CCR5 and CXCR4, respectively. The identified dimeric states differ in the access to the binding sites of the ligand and G protein, indicating that dimerization may represent a fine allosteric mechanism to regulate receptor activity. Our study offers structural basis for the design of ligands able to modulate the formation of CCR5 and CXCR4 dimers and in turn their activity, with therapeutic potential against HIV, cancer, and immune-inflammatory diseases.

Structural basis for triacylglyceride extraction from mycobacterial inner membrane by MFS transporter Rv1410

Mycobacterium tuberculosis is protected from antibiotic therapy by a multi-layered hydrophobic cell envelope. Major facilitator superfamily (MFS) transporter Rv1410 and the periplasmic lipoprotein LprG are involved in transport of triacylglycerides (TAGs) that seal the mycomembrane. Here, we report a 2.7 Å structure of a mycobacterial Rv1410 homologue, which adopts an outward-facing conformation and exhibits unusual transmembrane helix 11 and 12 extensions that protrude ~20 Å into the periplasm. A small, very hydrophobic cavity suitable for lipid transport is constricted by a functionally important ion-lock likely involved in proton coupling. Combining mutational analyses and MD simulations, we propose that TAGs are extracted from the core of the inner membrane into the central cavity via lateral clefts present in the inward-facing conformation. The functional role of the periplasmic helix extensions is to channel the extracted TAG into the lipid binding pocket of LprG.

Computational investigation on lipid bilayer disruption induced by amphiphilic Janus nanoparticles: combined effect of Janus balance and charged lipid concentration

Janus nanoparticles (NPs) with charged/hydrophobic compartments have garnered attention for their potential antimicrobial activity. These NPs have been shown to disrupt lipid bilayers in experimental studies, yet the underlying mechanisms of this disruption at the particle-membrane interface remain unclear. To address this knowledge gap, the present study conducts a computational investigation to systematically examine the disruption of lipid bilayers induced by amphiphilic Janus NPs. The focus of this study is on the combined effects of the hydrophobicity of the Janus NP, referred to as the Janus balance, defined as the ratio of hydrophilic to hydrophobic surface coverage, and the concentration of charged phospholipids on the interactions between Janus NPs and lipid bilayers. Computational simulations were conducted using a coarse-grained molecular dynamics (MD) approach. The results of these MD simulations reveal that while the area change of the bilayer increases monotonically with the Janus balance, the effect of charged lipid concentration in the membrane is not easy to be predicted. Specifically, it was found that the concentration of negatively charged lipids is directly proportional to the intensity of membrane disruption. Conversely, positively charged lipids have a negligible effect on membrane defects. This study provides molecular insights into the significant role of Janus balance in the disruption of lipid bilayers by Janus NPs and supports the selectivity of Janus NPs for negatively charged lipid membranes. Furthermore, the anisotropic properties of Janus NPs were found to play a crucial role in their ability to disrupt the membrane via the combination of hydrophobic and electrostatic interactions. This finding is validated by testing the current Janus NP design on a bacterial membrane-mimicking model. This computational study may serve as a foundation for further studies aimed at optimizing the properties of Janus NPs for specific antimicrobial applications.

Cardiolipin occupancy profiles of YidC paralogs reveal the significance of respective TM2 helix residues in determining paralog-specific phenotypes

YidC belongs to an evolutionarily conserved family of insertases, YidC/Oxa1/Alb3, in bacteria, mitochondria, and chloroplasts, respectively. Unlike Gram-negative bacteria, Gram-positives including Streptococcus mutans harbor two paralogs of YidC. The mechanism for paralog-specific phenotypes of bacterial YidC1 versus YidC2 has been partially attributed to the differences in their cytoplasmic domains. However, we previously identified a W138R gain-of-function mutation in the YidC1 transmembrane helix 2. YidC1W138R mostly phenocopied YidC2, yet the mechanism remained unknown. Primary sequence comparison of streptococcal YidCs led us to identify and mutate the YidC1W138 analog, YidC2S152 to W/A, which resulted in a loss of YidC2- and acquisition of YidC1-like phenotype. The predicted lipid-facing side chains of YidC1W138/YidC2S152 led us to propose a role for membrane phospholipids in specific-residue dependent phenotypes of S. mutans YidC paralogs. Cardiolipin (CL), a prevalent phospholipid in the S. mutans cytoplasmic membrane during acid stress, is encoded by a single gene, cls. We show a concerted mechanism for cardiolipin and YidC2 under acid stress based on similarly increased promoter activities and similar elimination phenotypes. Using coarse grain molecular dynamics simulations with the Martini2.2 Forcefield, YidC1 and YidC2 wild-type and mutant interactions with CL were assessed in silico. We observed substantially increased CL interaction in dimeric versus monomeric proteins, and variable CL occupancy in YidC1 and YidC2 mutant constructs that mimicked characteristics of the other wild-type paralog. Hence, paralog-specific amino acid- CL interactions contribute to YidC1 and YidC2-associated phenotypes that can be exchanged by point mutation at positions 138 or 152, respectively.

Phospholipid transport to the bacterial outer membrane through an envelope-spanning bridge

The outer membrane of Gram-negative bacteria provides a formidable barrier, essential for both pathogenesis and antimicrobial resistance. Biogenesis of the outer membrane requires the transport of phospholipids across the cell envelope. Recently, YhdP was implicated as a major protagonist in the transport of phospholipids from the inner membrane to the outer membrane however the molecular mechanism of YhdP mediated transport remains elusive. Here, utilising AlphaFold, we observe YhdP to form an elongated assembly of 60 β strands that curve to form a continuous hydrophobic groove. This architecture is consistent with our negative stain electron microscopy data which reveals YhdP to be approximately 250 Å in length and thus sufficient to span the bacterial cell envelope. Furthermore, molecular dynamics simulations and in vivo bacterial growth assays indicate essential helical regions at the N- and C-termini of YhdP, that may embed into the inner and outer membranes respectively, reinforcing its envelope spanning nature. Our in vivo crosslinking data reveal phosphate-containing substrates captured along the length of the YhdP groove, providing direct evidence that YhdP transports phospholipids. This finding is congruent with our molecular dynamics simulations which demonstrate the propensity for inner membrane lipids to spontaneously enter the groove of YhdP. Collectively, our results support a model in which YhdP bridges the cell envelope, providing a hydrophobic environment for the transport of phospholipids to the outer membrane.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

Effect of cetrimonium carrier micelles on bacterial membranes and extracellular DNA, an in silico study

Microorganisms do not live as dispersed single cells but rather they form aggregates with extracellular polymeric substances at interfaces. Biofilms are considered efficient life forms because they shield bacteria from biocides and collect dilute nutrients. This is a big concern in industry since the microorganisms can colonize a wide range of surfaces, accelerating material deterioration, colonizing medical devices, contaminating ultrapure drinking water, increasing energy costs and creating focus of infection. Conventional biocides that target a specific component of the bacteria are not effective in the presence of biofilms. Efficient biofilm inhibitors are based on a multitarget approach interacting with the bacteria and the biofilm matrix. Their rationale design requires a thorough understanding of inhibitory mechanisms that are still largely lacking today. Herein we uncover via molecular modelling the inhibition mechanism of cetrimonium 4-OH cinnamate (CTA-4OHcinn). Simulations show that CTA-4OH micelles can disrupt symmetric and asymmetric bilayers, representative of inner and outer bacterial membranes, following three stages: adsorption, assimilation, and defect formation. The main driving force for micellar attack is electrostatic interactions. In addition to disrupting the bilayers, the micelles work as carriers facilitating the trapping of 4OH cinnamate anions within the bilayer upper leaflet and overcoming electrostatic repulsion. The micelles also interact with extracellular DNA (e-DNA), which is one of the main components of biofilms. It is observed that CTA-4OHcinn forms spherical micelles on the DNA backbone; which hinders their ability to pack. This is demonstrated by modelling the DNA along the hbb histone-like protein, showing that in the presence of CTA-4OHcinn, DNA does not pack properly around hbb. The abilities of CTA-4OHcinn to cause cell death through membrane disruption and to disperse a mature, multi-species biofilm are also confirmed experimentally.

Nonlocal Effects of Antibiotic-Resistance-Causing Mutations Reveal an Alternative Region for Targeting on FtsW-Penicillin-Binding Protein 3 Complex of Haemophilus influenzae

Currently prescribed antibiotics target the catalytic sites of wild-type bacterial proteins; however, bacteria adopt mutations at this site, eventually leading to the emergence of resistance. Therefore, the identification of alternative drug binding sites is crucial, which requires knowledge of the dynamics of the mutant protein. Here, we set out to investigate the impact of a high-resistance-causing triple mutation (S385T + L389F + N526K) on the dynamics of a prioritized resistant pathogen, Haemophilus influenzae, using computational techniques. We studied penicillin-binding protein 3 (PBP3) and its complex with FtsW, which display resistance toward β-lactam antibiotics. We showed that mutations displayed local and nonlocal effects. In terms of the former, the orientation of the β-sheet, which surrounds the active site of PBP3, was impacted and the catalytic site was exposed to the periplasmic region. In addition, the flexibility of the β3-β4 loop, which modulates the catalysis of the enzyme, increased in the mutant FtsW-PBP3 complex. As for nonlocal effects, the dynamics of the pedestal domain (N-terminal periplasmic modulus (N-t)), i.e., the opening of the fork, was different between the wild-type and mutant enzymes. We showed the closed fork caused a greater number of residues to participate in the hypothesized allosteric communication network connecting N-t to the transpeptidase domain in the mutant enzyme. Finally, we demonstrated that the closed fork results in more favorable binding with β-lactam antibiotics, particularly cefixime, suggesting that small therapeutics that can stabilize the closed fork of mutant PBP3 may lead to the development of more effective molecules to combat resistant bacteria.

HylleraasMD: A Domain Decomposition-Based Hybrid Particle-Field Software for Multiscale Simulations of Soft Matter

We present HylleraasMD (HyMD), a comprehensive implementation of the recently proposed Hamiltonian formulation of hybrid particle-field molecular dynamics. The methodology is based on a tunable, grid-independent length-scale of coarse graining, obtained by filtering particle densities in reciprocal space. This enables systematic convergence of energies and forces by grid refinement, also eliminating nonphysical force aliasing. Separating the time integration of fast modes associated with internal molecular motion from slow modes associated with their density fields, we enable the first time-reversible, energy-conserving hybrid particle-field simulations. HyMD comprises the optional use of explicit electrostatics, which, in this formalism, corresponds to the long-range potential in particle-mesh Ewald. We demonstrate the ability of HyMD to perform simulations in the microcanonical and canonical ensembles with a series of test cases, comprising lipid bilayers and vesicles, surfactant micelles, and polypeptide chains, comparing our results to established literature. An on-the-fly increase of the characteristic coarse-grain length significantly speeds up dynamics, accelerating self-diffusion and leading to expedited aggregation. Exploiting this acceleration, we find that the time scales involved in the self-assembly of polymeric structures can lie in the tens to hundreds of picoseconds instead of the multimicrosecond regime observed with comparable coarse-grained models.

Bilayer lipid membrane formation on surface assemblies with sparsely distributed tethers

A combined computational and experimental study of small unilamellar vesicle (SUV) fusion on mixed self-assembled monolayers (SAMs) terminated with different deuterated tether moieties (-(CD₂)₇CD₃ or -(CD₂)₁₅CD₃) is reported. Tethered bilayer lipid membrane (tBLM) formation of synthetic 1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine was initially probed on SAMs with controlled tether (d-alkyl tail) surface densities and lateral molecular packing using quartz crystal microbalance with dissipation monitoring (QCM-D). Long time-scale coarse-grained molecular dynamics (MD) simulations were then employed to elucidate the mechanisms behind the interaction between the SUVs and the different phases formed by the -(CD₂)₇CD₃ and -(CD₂)₁₅CD₃ tethers. Furthermore, a series of real time kinetics was recorded under different osmotic conditions using QCM-D to determine the accumulated lipid mass and for probing the fusion process. It is shown that the key factors driving the SUV fusion and tBLM formation on this type of surfaces involve tether insertion into the SUVs along with vesicle deformation. It is also evident that surface densities of the tethers as small as a few mol% are sufficient to obtain stable tBLMs with a high reproducibility. The described "sparsely tethered" tBLM system can be advantageous in studying different biophysical phenomena, such as membrane protein insertion, effects of receptor clustering, and raft formation.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Soft Matter under Pressure: Pushing Particle–Field Molecular Dynamics to the Isobaric Ensemble

Hamiltonian hybrid particle-field molecular dynamics is a computationally efficient method to study large soft matter systems. In this work, we extend this approach to constant-pressure (NPT) simulations. We reformulate the calculation of internal pressure from the density field by taking into account the intrinsic spread of the particles in space, which naturally leads to a direct anisotropy in the pressure tensor. The anisotropic contribution is crucial for reliably describing the physics of systems under pressure, as demonstrated by a series of tests on analytical and monatomic model systems as well as realistic water/lipid biphasic systems. Using Bayesian optimization, we parametrize the field interactions of phospholipids to reproduce the structural properties of their lamellar phases, including area per lipid, and local density profiles. The resulting model excels in providing pressure profiles in qualitative agreement with all-atom modeling, and surface tension and area compressibility in quantitative agreement with experimental values, indicating the correct description of long-wavelength undulations in large membranes. Finally, we demonstrate that the model is capable of reproducing the formation of lipid droplets inside a lipid bilayer.

UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape

The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.

Insights into membrane association of the SMP domain of extended synaptotagmin

The Synaptotagmin-like Mitochondrial-lipid-binding Protein (SMP) domain is a newly identified lipid transfer module present in proteins that regulate lipid homeostasis at membrane contact sites (MCSs). However, how the SMP domain associates with the membrane to extract and unload lipids is unclear. Here, we performed in vitro DNA brick-assisted lipid transfer assays and in silico molecular dynamics simulations to investigate the molecular basis of the membrane association by the SMP domain of extended synaptotagmin (E-Syt), which tethers the tubular endoplasmic reticulum (ER) to the plasma membrane (PM). We demonstrate that the SMP domain uses its tip region to recognize the extremely curved subdomain of tubular ER and the acidic-lipid-enriched PM for highly efficient lipid transfer. Supporting these findings, disruption of these mechanisms results in a defect in autophagosome biogenesis contributed by E-Syt. Our results suggest a model that provides a coherent picture of the action of the SMP domain at MCSs.

Facilitating CG Simulations with MAD: The MArtini Database Server

The MArtini Database (MAD - https://mad.ibcp.fr) is a web server designed for the sharing of structures and topologies of molecules parametrized with the Martini coarse-grained (CG) force field. MAD can also convert atomistic structures into CG structures and prepare complex systems (including proteins, lipids, etc.) for molecular dynamics (MD) simulations at the CG level. It is dedicated to the generation of input files for Martini 3, the most recent version of this popular CG force field. Specifically, the MAD server currently includes tools to submit or retrieve CG models of a wide range of molecules (lipids, carbohydrates, nanoparticles, etc.), transform atomistic protein structures into CG structures and topologies, with fine control on the process and assemble biomolecules into large systems, and deliver all files necessary to start simulations in the GROMACS MD engine.

Capturing coacervate formation and protein partition by molecular dynamics simulation

Biomolecules localize and function in microenvironments where their local concentration, spatial organization, and biochemical reactivity are regulated. To compartmentalize and control the local properties of the native microenvironment, cellular mimics and artificial bioreactors have been developed to study the properties of membraneless organelles or mimic the bio-environment for life origin. Here, we carried out molecular dynamics simulation with the Martini 3.0 model to reproduce the experimental salt concentration and pH dependency of different complex coacervates. We showed that coacervates inside vesicles are able to change their shape. In addition, we used these coacervate systems to explore the partitioning of the ubiquitous cytoskeletal protein actin and found that actin spontaneously partitions to all the coacervate peripheries. Therefore, we believe that our study can provide a better understanding of the versatile coacervate platform, where biomolecules partition and gather to fulfill their biological duties.

Structural Determinants and Functional Significance of Dimerization for Osmosensing Transporter ProP in Escherichia coli

Osmosensing transporter ProP forestalls cellular dehydration by detecting environments with high osmotic pressure and mediating the accumulation of organic osmolytes by bacterial cells. It is composed of 12 transmembrane helices with cytoplasmic N- and C-termini. In Escherichia coli, dimers form when the C-terminal domains of ProP molecules form homodimeric, antiparallel, α-helical coiled coils. No dominant negative effect was detected when inactive and active ProP molecules formed heterodimers in vivo. Purification of ProP in detergent dodecylmaltoside yielded monomers, which were functional after reconstitution in proteoliposomes. With other evidence, this suggests that ProP monomers function independently whether in the monomeric or dimeric state. Amino acid replacements that disrupted or reversed the coiled coil did not prevent in vivo dimerization of ProP detected with a bacterial two-hybrid system. Maleimide labeling detected no osmolality-dependent variation in the reactivities of cysteine residues introduced to transmembrane helix (TM) XII. In contrast, coarse-grained molecular dynamic simulations detected deformation of the lipid around TMs III and VI, on the lipid-exposed protein surface opposite to TM XII. This suggests that the dimer interface of ProP includes the surfaces of TMs III and VI, not of TM XII as previously suggested by crosslinking data. Homology modeling suggested that coiled-coil formation and dimerization via such an interface are not mutually exclusive. In previous work, alterations to the C-terminal coiled coil blocked co-localization of ProP with phospholipid cardiolipin at E. coli cell poles. Thus, dimerization may contribute to ProP targeting, adjust its lipid environment, and hence indirectly modify its osmotic stress response.

A brief history of visualizing membrane systems in molecular dynamics simulations

Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.

Martini 3 Coarse-Grained Force Field for Carbohydrates

The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme which decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals were developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the polyglucoses dextran (water-soluble) and cellulose (water insoluble but soluble in ionic liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids. We show they are able to reproduce membrane properties and induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.

Mutational Insight into Allosteric Regulation of Kir Channel Activity

Potassium (K⁺) channels are regulated in part by allosteric communication between the helical bundle crossing, or inner gate, and the selectivity filter, or outer gate. This network is triggered by gating stimuli. In concert, there is an allosteric network which is a conjugated set of interactions which correlate long-range structural rearrangements necessary for channel function. Inward-rectifier K⁺ (Kir) channels favor inward K⁺ conductance, are ligand-gated, and help establish resting membrane potentials. KirBac1.1 is a bacterial Kir (KirBac) channel homologous to human Kir (hKir) channels. Additionally, KirBac1.1 is gated by the anionic phospholipid ligand phosphatidylglycerol (PG). In this study, we use site-directed mutagenesis to investigate residues involved in the KirBac1.1 gating mechanism and allosteric network we previously proposed using detailed solid-state NMR (SSNMR) measurements. Using fluorescence-based K⁺ and sodium (Na⁺) flux assays, we identified channel mutants with impaired function that do not alter selectivity of the channel. In tandem, we performed coarse grain molecular dynamics simulations, observing changes in PG-KirBac1.1 interactions correlated with mutant channel activity and contacts between the two transmembrane helices and pore helix tied to this behavior. Lipid affinity is closely tied to the proximity of two tryptophan residues on neighboring subunits which lure anionic lipids to a cationic pocket formed by a cluster of arginine residues. Thus, these simulations establish a structural and functional basis for the role of each mutated site in the proposed allosteric network. The experimental and simulated data provide insight into key functional residues involved in gating and lipid allostery of K⁺ channels. Our findings also have direct implications on the physiology of hKir channels due to conservation of many of the residues identified in this work from KirBac1.1.

Donor-strand exchange drives assembly of the TasA scaffold in Bacillus subtilis biofilms

Many bacteria in nature exist in multicellular communities termed biofilms, where cells are embedded in an extracellular matrix that provides rigidity to the biofilm and protects cells from chemical and mechanical stresses. In the Gram-positive model bacterium Bacillus subtilis, TasA is the major protein component of the biofilm matrix, where it has been reported to form functional amyloid fibres contributing to biofilm structure and stability. Here, we present electron cryomicroscopy structures of TasA fibres, which show that, rather than forming amyloid fibrils, TasA monomers assemble into fibres through donor-strand exchange, with each subunit donating a β-strand to complete the fold of the next subunit along the fibre. Combining electron cryotomography, atomic force microscopy, and mutational studies, we show how TasA fibres congregate in three dimensions to form abundant fibre bundles that are essential for B. subtilis biofilm formation. Our study explains the previously observed biochemical properties of TasA and shows how a bacterial extracellular globular protein can assemble from monomers into β-sheet-rich fibres, and how such fibres assemble into bundles in biofilms.

Lipids mediate supramolecular outer membrane protein assembly in bacteria

β Barrel outer membrane proteins (OMPs) cluster into supramolecular assemblies that give function to the outer membrane (OM) of Gram-negative bacteria. How such assemblies form is unknown. Here, through photoactivatable cross-linking into the Escherichia coli OM, coupled with simulations, and biochemical and biophysical analysis, we uncover the basis for OMP clustering in vivo. OMPs are typically surrounded by an annular shell of asymmetric lipids that mediate higher-order complexes with neighboring OMPs. OMP assemblies center on the abundant porins OmpF and OmpC, against which low-abundance monomeric β barrels, such as TonB-dependent transporters, are packed. Our study reveals OMP-lipid-OMP complexes to be the basic unit of supramolecular OMP assembly that, by extending across the entire cell surface, couples the requisite multifunctionality of the OM to its stability and impermeability.

Understanding the Role of Self-Assembly and Interaction with Biological Membranes of Short Cationic Lipopeptides in the Effective Design of New Antibiotics

This study investigates short cationic antimicrobial lipopeptides composed of 2-4 amino acid residues and C12-C18 fatty acids attached to the N-terminal part of the peptides. The findings were discussed in the context of the relationship among biological activity, self-assembly, stability, and membrane interactions. All the lipopeptides showed the ability to self-assemble in PBS solution. In most cases, the critical aggregation concentration (CAC) much surpassed the minimal inhibitory concentration (MIC) values, suggesting that monomers are the main active form of lipopeptides. The introduction of β-alanine into the peptide sequence resulted in a compound with a high propensity to fibrillate, which increased the peptide stability and activity against S. epidermidis and C. albicans and reduced the cytotoxicity against human keratinocytes. The results of our study indicated that the target of action of lipopeptides is the bacterial membrane. Interestingly, the type of peptide counterion may affect the degree of penetration of the lipid bilayer. In addition, the binding of the lipopeptide to the membrane of Gram-negative bacteria may lead to the release of calcium ions necessary for stabilization of the lipopolysaccharide layer.

Differential interactions of resting, activated, and desensitized states of the α7 nicotinic acetylcholine receptor with lipidic modulators

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that modulates neuronal excitability, largely by allowing Ca2+ permeation. Agonist binding promotes transition from a resting state to an activated state, and then rapidly to a desensitized state. Recently, cryogenic electron microscopy (cryo-EM) structures of the human α7 receptor in nanodiscs were reported in multiple conformations. These were selectively stabilized by inhibitory, activating, or potentiating compounds. However, the functional annotation of these structures and their differential interactions with unresolved lipids and ligands remain incomplete. Here, we characterized their ion permeation, membrane interactions, and ligand binding using computational electrophysiology, free-energy calculations, and coarse-grained molecular dynamics. In contrast to nonconductive structures in apparent resting and desensitized states, the structure determined in the presence of the potentiator PNU-120596 was consistent with an activated state permeable to Ca2+. Transition to this state was associated with compression and rearrangement of the membrane, particularly in the vicinity of the peripheral MX helix. An intersubunit transmembrane site was implicated in selective binding of either PNU-120596 in the activated state or cholesterol in the desensitized state. This substantiates functional assignment of all three lipid-embedded α7-receptor structures with ion-permeation simulations. It also proposes testable models of their state-dependent interactions with lipophilic ligands, including a mechanism for allosteric modulation at the transmembrane subunit interface.

Cecropin a Improves the Antibacterial Activity of Hen Egg White Lysozyme against Challenging Salmonella enterica Serovars

The prevalence of multidrug-resistant Salmonella enterica among animal- and plant-derived food products threatens global healthcare and economic sectors. Hen egg white lysozyme is widely exploited as a food preservative against Gram-positive pathogens. Nevertheless, its limited penetration of the outer membrane renders it ineffective against Gram-negative bacteria. Herein, we present a safe and effective approach to facilitate HEWL access to peptidoglycan layers using cecropin A. In silico analysis of cecropin A peptide revealed an amphipathic α-helical peptide with potential outer membrane permeabilizing activity through its interaction with both hydrophobic and ionic stabilizing forces. Evaluation of HEWL/cecropin A combination showed a cecropin A dose-dependent bacterial count reduction up to 4.16 and 3.18 ± 0.26 log units against Salmonella enterica ATCC 35664 at the logarithmic and stationary growth phases, respectively. Moreover, the combination displayed antibacterial activity of 2.1 ± 0.31 and ~1 log-unit reductions against Salmonella enterica serovars Kentucky, Typhimurium, and Enteritidis, respectively, whereas Hato and Shangani were found irresponsive. The cytotoxicity assay revealed compatibility of cecropin A with oral epithelial cells. These observations suggest HEWL/cecropin A combination as an effective and safe alternative to lysozyme against Salmonella enterica.

Simulation Study of the Effect of Antimicrobial Peptide Associations on the Mechanism of Action with Bacterial and Eukaryotic Membranes

Antimicrobial peptides (AMPs) can be directed to specific membranes based on differences in lipid composition. In this study, we performed atomistic and coarse-grained simulations of different numbers of the designed AMP adepantin-1 with a eukaryotic membrane, cytoplasmic Gram-positive and Gram-negative membranes, and an outer Gram-negative membrane. At the core of adepantin-1's behavior is its amphipathic α-helical structure, which was implemented in its design. The amphipathic structure promotes rapid self-association of peptide in water or upon binding to bacterial membranes. Aggregates initially make contact with the membrane via positively charged residues, but with insertion, the hydrophobic residues are exposed to the membrane's hydrophobic core. This adaptation alters the aggregate's stability, causing the peptides to diffuse in the polar region of the membrane, mostly remaining as a single peptide or pairing up to form an antiparallel dimer. Thus, the aggregate's proposed role is to aid in positioning the peptide into a favorable conformation for insertion. Simulations revealed the molecular basics of adepantin-1 binding to various membranes, and highlighted peptide aggregation as an important factor. These findings contribute to the development of novel anti-infective agents to combat the rapidly growing problem of bacterial resistance to antibiotics.

Design and Evaluation of Short Bovine Lactoferrin-Derived Antimicrobial Peptides against Multidrug-Resistant Enterococcus faecium

Enterococcus faecium has become an important drug-resistant nosocomial pathogen because of widespread antibiotic abuse. We developed short and chemically simple antimicrobial peptides (AMPs) with a selective amino acid composition, fixed charge, and hydrophobicity ratio based on the core antimicrobial motif of bovine lactoferrin (LfcinB6). Among these peptides, 5L and 6L (both 12 residues long) demonstrated a narrow spectrum and high antibacterial activity against drug-resistant E. faecium isolates with a minimal inhibitory concentration (MIC) that ranged from 4–16 µg/mL. At 32 µg/mL, peptides 5L and 6L inhibited E. faecium strain C68 biofilm formation by 90% and disrupted established biofilms by 75%. At 40 µg/mL, 5L reduced 1 × 107E. faecium persister cells by 3 logs within 120 min of exposure, whereas 6L eliminated all persister cells within 60 min. At 0.5× MIC, 5L and 6L significantly downregulated the expression of a crucial biofilm gene ace by 8 folds (p = 0.02) and 4 folds (p = 0.01), respectively. At 32 µg/mL, peptides 5L and 6L both depolarized the E. faecium membrane, increased fluidity, and eventually ruptured the membrane. Physiologically, 5L (at 8 µg/mL) altered the tricarboxylic acid cycle, glutathione, and purine metabolism. Interestingly, in an ex vivo model of porcine skin infection, compared to no treatment, 5L (at 10× MIC) effectively eliminated all 1 × 106 exponential (p = 0.0045) and persister E. faecium cells (p = 0.0002). In conclusion, the study outlines a roadmap for developing narrow-spectrum selective AMPs and presents peptide 5L as a potential therapeutic candidate to be explored against E. faecium.

Interaction between Antimicrobial Peptide Magainin 2 and Nonlipid Components in the Bacterial Outer Envelope

Antimicrobial peptides (AMPs) offer advantages over conventional antibiotics; for example, bacteria develop more resistance to small-molecule antibiotics than to AMPs. The interaction of the AMPs with the lipopolysaccharide (LPS) layer of the Gram-negative bacteria cell envelope is not well understood. A MARTINI model was constructed of a Gram-negative bacterial outer membrane interacting with the AMP Magainin 2. In a 20 μs molecular dynamics (MD) simulation, the AMP diffused to the LPS layer of the cell envelope and remained there, suggesting interactions between the Magainin 2 and the LPS layer, causing the AMP to concentrate at that position. The free energy profile for the insertion of the Magainin 2 into the membrane was also calculated using umbrella sampling, which showed that the AMP positioned such that the cationic side chains of the AMP coordinated to the negatively charged phosphate groups of the LPS layer. These simulations indicate that the AMP Magainin 2 partition into the LPS layer of a bacterial membrane.