Free energy of WALP23 dimer association in DMPC, DPPC, and DOPC bilayers

Chlorophyll-Induced Lamellar to Nonlamellar Phase Transitions and Dynamical Heterogeneity in Plant Thylakoid Membranes

Chlorophyll a (CLA) pigments and thylakoid membranes are crucial components of plants for photosynthesis. To understand the effect of CLA on the structure and dynamics of thylakoid membranes, coarse-grained molecular dynamics (CG MD) simulations of thylakoid membranes are performed by varying the numbers of CLA at 293 K using MARTINI-2 force fields. The membrane undergoes a lamellar to nonlamellar phase transition above a critical concentration of CLA. The CLAs dynamically form aggregates of different orders and preferentially fetch the least unsaturated nonbilayer-forming lipids around them, resulting in a nonlamellar phase with fused regions. These fused regions cause a structural arrest of CLA and lipids, inducing dynamic heterogeneity manifested by non-Gaussian parameters and van Hove correlation functions. The lamellar to nonlamellar phase transition of the membrane is associated with a drastic reduction in correlation length of the immobile CLA and lipids governed by the fused topology. Such insights into CLA-induced structural transitions in thylakoid membranes are pertinent for understanding nonphotochemical quenching mechanisms and hold promise for designing future artificial photosynthetic materials and applications in photodynamic therapy.

Membrane remodeling by FAM92A1 during brain development regulates neuronal morphology, synaptic function, and cognition

The Bin/Amphiphysin/Rvs (BAR) domain protein FAM92A1 is a multifunctional protein engaged in regulating mitochondrial ultrastructure and ciliogenesis, but its physiological role in the brain remains unclear. Here, we show that FAM92A1 is expressed in neurons starting from embryonic development. FAM92A1 knockout in mice results in altered brain morphology and age-associated cognitive deficits, potentially due to neuronal degeneration and disrupted synaptic plasticity. Specifically, FAM92A1 deficiency impairs diverse neuronal membrane morphology, including the mitochondrial inner membrane, myelin sheath, and synapses, indicating its roles in membrane remodeling and maintenance. By determining the crystal structure of the FAM92A1 BAR domain, combined with atomistic molecular dynamics simulations, we uncover that FAM92A1 interacts with phosphoinositide- and cardiolipin-containing membranes to induce lipid-clustering and membrane curvature. Altogether, these findings reveal the physiological role of FAM92A1 in the brain, highlighting its impact on synaptic plasticity and neural function through the regulation of membrane remodeling and endocytic processes.

Machine Learning Derived Collective Variables for the Study of Protein Homodimerization in Membrane

The accurate calculation of equilibrium constants for protein-protein association is of fundamental importance to quantitative biology and remains an outstanding challenge for computational biophysics. Traditionally, equilibrium constants have been computed from one-dimensional free energy surfaces derived from sampling along a single collective variable. Importantly, recent advances in enhanced sampling methodology have facilitated the characterization of multidimensional free energy landscapes, often exposing multiple thermodynamically important minima missed by more restrictive sampling methods. A key to the effectiveness of this multidimensional sampling approach is the identification of collective variables that effectively define the configurational space of dissociated and associated states. Here we present the application of two machine learning methods for the unbiased determination of collective variables for enhanced sampling and analysis of protein-protein association. Our results both validate prior work, based on intuition derived collective variables, and demonstrate the effectiveness of the machine learning methods for the identification of collective variables for association reactions in complex biomolecular systems.

Energetics of the Transmembrane Peptide Sorting by Hydrophobic Mismatch

Hydrophobic mismatch between a lipid membrane and embedded transmembrane peptides or proteins plays a role in their lateral localization and function. Earlier studies have resolved numerous mechanisms through which the peptides and membrane proteins adapt to mismatch, yet the energetics of lateral sorting due to hydrophobic mismatch have remained elusive due to the lack of suitable computational or experimental protocols. Here, we pioneer a molecular dynamics simulation approach to study the sorting of peptides along a membrane thickness gradient. Peptides of different lengths tilt and diffuse along the membrane to eliminate mismatch with a rate directly proportional to the magnitude of mismatch. We extract the 2-dimensional free energy profiles as a function of local thickness and peptide orientation, revealing the relative contributions of sorting and tilting, and suggesting their thermally accessible regimes. Our approach can readily be applied to study other membrane systems of biological interest where hydrophobic mismatch, or membrane thickness in general, plays a role.

Molecular dynamics simulation of apolipoprotein E3 lipid nanodiscs

Nanodiscs are binary discoidal complexes of a phospholipid bilayer circumscribed by belt-like helical scaffold proteins. Using coarse-grained and all-atom molecular dynamics simulations, we explore the stability, size, and structure of nanodiscs formed between the N-terminal domain of apolipoprotein E3 (apoE3-NT) and variable number of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) molecules. We study both parallel and antiparallel double-belt configurations, consisting of four proteins per nanodisc. Our simulations predict nanodiscs containing between 240 and 420 DMPC molecules to be stable. The antiparallel configurations exhibit an average of 1.6 times more amino acid interactions between protein chains and 2 times more ionic contacts, compared to the parallel configuration. With one exception, DMPC order parameters are consistently larger in the antiparallel configuration than in the parallel one. In most cases, the root mean square deviation of the positions of the protein backbone atoms is smaller in the antiparallel configuration. We further report nanodisc size, thickness, radius of gyration, and solvent accessible surface area. Combining all investigated parameters, we hypothesize the antiparallel protein configuration leading to more stable and more rigid nanodiscs than the parallel one.

Thylakoid Composition Facilitates Chlorophyll a Dimerization through Stronger Interlipid Interactions

Plant thylakoid membrane serves as a crucial matrix for the aggregation of chlororophyll a (CLA) pigments, essential for light harvesting. To understand the role of lipid compositions in the stability of CLA aggregates, dimerization of chlorophyll a molecules (CLA) is studied in the presence of the thylakoid and the bilayers comprising either the least or the highest unsaturated lipids by using coarse-grained molecular dynamics simulations. The thylakoid membrane enhances the stability of the CLA dimer compared with other membranes due to very strong lipid-lipid interactions. The thylakoid exhibits a distinct distribution of lipids around the CLA dimer. Less unsaturated lipids reside in close proximity to the dimer, promoting increased order and efficient packing. Conversely, higher unsaturated lipids are depleted from the dimer, imparting flexibility to the membrane. The combination of tight packing near the dimer and membrane flexibility away from the dimer enhances the stability of the dimer in the thylakoid membrane. Our results suggest that lipid mixing, rather than lipid unsaturation, plays a critical role in facilitating CLA dimerization by modulating the membrane microenvironment through stronger lipid-lipid interactions. These insights will be useful in understanding how lipid compositions affect efficient light absorption and energy transfer during photosynthesis in the future.

The role of structural heterogeneity in the homodimerization of transmembrane proteins

The equilibrium association of transmembrane proteins plays a fundamental role in membrane protein function and cellular signaling. While the study of the equilibrium binding of single pass transmembrane proteins has received significant attention in experiment and simulation, the accurate assessment of equilibrium association constants remains a challenge to experiment and simulation. In experiment, there remain wide variations in association constants derived from experimental studies of the most widely studied transmembrane proteins. In simulation, state-of-the art methods have failed to adequately sample the thermodynamically relevant structures of the dimer state ensembles using coarse-grained models. In addition, all-atom force fields often fail to accurately assess the relative free energies of the dimer and monomer states. Given the importance of this fundamental biophysical process, it is essential to address these shortcomings. In this work, we establish an effective computational protocol for the calculation of equilibrium association constants for transmembrane homodimer formation. A set of transmembrane protein homodimers, used in the parameterization of the MARTINI v3 force field, are simulated using metadynamics, based on three collective variables. The method is found to be accurate and computationally efficient, providing a standard to be used in the future simulation studies using coarse-grained or all-atom models.

The Impact of Antimicrobial Peptides on the Acinetobacter baumannii Inner Membrane Is Modulated by Lipid Polyunsaturation

The Gram-negative pathogen Acinetobacter baumannii is a primary contributor to nosocomial multi-drug-resistant (MDR) infections. To combat the rise of MDR infections, novel features of A. baumannii need to be considered for the development of new treatment options. One such feature is the preferential scavenging of exogenous lipids, including host-derived polyunsaturated fatty acids (PUFAs), for membrane phospholipid synthesis. These alterations in membrane composition impact both the lipid chemistry and the membrane biophysical properties. In this work we examine how antimicrobial peptides (AMPs) interact with the inner membranes of A. baumannii in the presence and absence of polyunsaturated phospholipids. Using coarse-grained molecular dynamics simulations of complex A. baumannii inner membrane models derived from lipidomes of bacteria grown in the presence and absence of PUFAs, we examine the impact of the adsorption of four prototypical AMPs (CAMEL, LL-37, pexiganan, and magainin-2) on the membrane biophysical properties. Our simulations reveal that the impact of AMP adsorption on the membrane biophysical properties was dependent on both the membrane composition and the specific AMP involved. Both lipid headgroup charge and tail unsaturation played important roles in driving the interactions that occurred both within the membrane and between the membrane and AMPs. The changes to the membrane biophysical properties also showed a complex relationship with the AMP's physical properties, such as AMP charge, chain length, and charge-to-mass ratio. Cumulatively, this work highlights the importance of studying AMPs using a complex membrane environment and provides insights into the mechanistic action of AMPs in polyunsaturated lipid-rich bacterial membranes.

Effects of Gly Residue and Cholesterol on the GXXXG-Mediated Parallel Association of Transmembrane Helices: A Single-Pair FRET Study

Small residue-mediated interhelical packing is ubiquitous in helical membrane proteins: however, the lipid dependence of its stability remains unclear. We previously demonstrated that the introduction of a GXXXG sequence in the middle of de novo-designed (AALALAA)₃ helices (AALALAA AGLALGA AALALAA) facilitated their dimerization, which was abolished by cholesterol. Here single-pair FRET measurements revealed that a longer GXXXGXXXG segment (AALALAA A GLALGA AAGALAA) promoted helix dimerization in POPC/cholesterol bilayers, but not without cholesterol. The predicted dimer structures and degrees of helix packing suggested that helix dimers with small (∼10°) and large (∼55°) crossing angles were only stabilized in POPC and POPC/cholesterol membranes, respectively. A steric hindrance in the dimer interface and the large flexibility of helices prevented the formation of stable dimers. Therefore, amino acid sequences and lipid compositions distinctively constrain stable dimer structures in membranes.

Formation of lipid raft nanodomains in homogeneous ternary lipid mixture of POPC/DPSM/cholesterol: Theoretical insights

Lipid rafts, in biological membranes, are cholesterol-rich nanodomains that regulate many protein activities and cellular processes. Understanding the formation of the lipid-raft nanodomains helps us elucidate many complex interactions in the cell. In this study, the formation of lipid-raft nanodomains in a ternary palmitoyl-oleoyl-phosphatidylcholine/stearoyl-sphingomyelin/cholesterol (POPC/DPSM/Chol) lipid mixture, the most realistic surrogate model for biological membranes, has been successfully observed for the first time in-silico using microsecond timescale molecular dynamics simulations. The model reveals the formation of cholesterol-induced nanodomains with raft-like characteristics and their underlying mechanism: the cholesterol molecules segregate themselves into cholesterol nanodomains and then enrich the cholesterol-rich domain with sphingomyelin molecules to form a lipid raft thanks to the weak bonding of cholesterol with sphingomyelin. Besides, it is found that the increase in cholesterol concentration enhances the biophysical properties (e.g., bilayer thickness, area per lipid headgroup, and order parameter) of the lipid raft nanodomains. Such findings suggest that the POPC/DPSM/Chol bilayer is a suitable model to fundamentally extend the nanodomain evolution to investigate their lifetime and kinetics as well as to study protein-lipid interaction, protein-protein interaction, and selection of therapeutic molecules in the presence of lipid rafts.

Polyvinyl Chloride Nanoparticles Affect Cell Membrane Integrity by Disturbing the Properties of the Multicomponent Lipid Bilayer in Arabidopsis thaliana

The ubiquitous presence of nanoplastics (NPs) in natural ecosystems is a serious concern, as NPs are believed to threaten every life form on Earth. Micro- and nanoplastics enter living systems through multiple channels. Cell membranes function as the first barrier of entry to NPs, thus playing an important biological role. However, in-depth studies on the interactions of NPs with cell membranes have not been performed, and effective theoretical models of the underlying molecular details and physicochemical behaviors are lacking. In the present study, we investigated the uptake of polyvinyl chloride (PVC) nanoparticles by Arabidopsis thaliana root cells, which leads to cell membrane leakage and damage to membrane integrity. We performed all-atom molecular dynamics simulations to determine the effects of PVC NPs on the properties of the multicomponent lipid bilayer. These simulations revealed that PVCs easily permeate into model lipid membranes, resulting in significant changes to the membrane, including reduced density and changes in fluidity and membrane thickness. Our exploration of the interaction mechanisms between NPs and the cell membrane provided valuable insights into the effects of NPs on membrane structure and integrity.

Dynamics of the Acinetobacter baumannii inner membrane under exogenous polyunsaturated fatty acid stress

Exogenous polyunsaturated fatty acids (PUFAs) are readily incorporated into the synthesis pathways of A. baumannii membrane phospholipids, where they contribute to reduced bacterial fitness and increased antimicrobial susceptibility. Here we examine the impact of PUFA membrane modification on membrane organisation and biophysical properties using coarse grained MARTINI simulations of chemically representative membrane models developed from mass-spectrometry datasets of an untreated, arachidonic acid (AA) treated and docosahexaenoic acid (DHA) treated A. baumannii membranes. Enzymatic integration of AA or DHA into phospholipids of the A. baumannii membrane resulted in modulation of membrane biophysical properties. Membrane thickness decreased slightly following PUFA treatment, concomitant with changes in the lateral area per lipid of each lipid headgroup class. PUFA treatment resulted in a decrease in membrane ordering and an increase in lipid lateral diffusion. Changes in lateral membrane organisation were observed in the PUFA treated membranes, with a concurrent increase in ordered cardiolipin domains and disordered PUFA-containing domains. Notably, separation between ordered and disordered domains was enhanced and was more pronounced for DHA relative to AA, providing a possible mechanism for greater antimicrobial action of DHA relative to AA observed experimentally. Furthermore, the membrane active antimicrobial, pentamidine, preferentially adsorbs to cardiolipin domains of the A. baumannii model membranes. This interaction, and membrane penetration of pentamidine, was enhanced following PUFA treatment. Cumulatively, this work explores the wide-ranging effects of PUFA incorporation on the A. baumannii membrane and provides a molecular basis for bacterial inner membrane disruption by PUFAs.

SPICA Force Field for Proteins and Peptides

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

Elucidating the Structural Features of ABCA1 in its Heterogeneous Membrane Environment

ATP Binding Cassette Transporter A1 (ABCA1) plays an integral part in Reverse Cholesterol Transport (RCT) and is critical for maintaining lipid homeostasis. One theory of lipid efflux by the transporter (alternating access) proposes that ABCA1 harbours two different conformations that provide alternating access for lipid binding and release. This is followed by sequestration via a direct interaction between ABCA1 and its partner, ApoA1. The other theory (lateral access) proposes that ABCA1 obtains lipids laterally from the membrane to form a temporary extracellular “reservoir”. This reservoir contains an isolated lipid monolayer due to the net accumulation of lipids in the exofacial leaflet. Recently, a full-length Cryo-EM structure of this 2,261-residue transmembrane protein showed its discreetly folded domains and have detected the presence of a tunnel enclosed within the extracellular domains (ECDs) but not in the TMDs, giving it an outward-facing conformation. This structure was hypothesized to substantiate the lateral access theory. Utilizing long time-scale multiple replica atomistic molecular dynamics simulations (MDS), we simulated the structure in a large heterogeneous lipid environment and found that the protein undergoes several large conformational changes in its extremities. We observed that the cavity enclosed within ATP unbound form of ABCA1 is narrow at the distal ends of TMD as well as the ECD region substantiating the “lateral access” theory. We have also characterized ABCA1 and the lipid dynamics along with the protein-lipid interactions in the heterogeneous environment, providing novel insights into understanding ABCA1 conformation at an atomistic level.

Partitioning of a Hybrid Lipid in Domains of Saturated and Unsaturated Lipids in a Model Cellular Membrane

The cellular membranes are composed of hundreds of components such as lipids, proteins, and sterols that are chemically and physically distinct from each other. The lipid-lipid and lipid-protein interactions form domains in this membrane, which play vital roles in membrane physiology. The hybrid lipids (HLs) with one saturated and one unsaturated chain can control the shape and size of these domains, ensuring the thermodynamic stability of a membrane. In this study, the thermodynamics of mixing of a HL and its structural effects on the phase separated domains in a model membrane composed of a saturated and an unsaturated lipid have been investigated. The HL is observed to mix into an unsaturated lipid reducing the Gibbs free energy, whereas the mixing is unfavorable in a saturated lipid. The presence of an HL in an unsaturated lipid tends to increase its area fraction, which is reflected in the enhanced correlation length across the bilayers in a multilayered sample. There is a feeble effect on the domain structure of the saturated lipid due to the presence of the HLs at the phase boundary. This study concludes that the HLs preferentially participate in the unsaturated lipid regions compared to that of a saturated lipid.

Genetic, cellular, and structural characterization of the membrane potential-dependent cell-penetrating peptide translocation pore

Cell-penetrating peptides (CPPs) allow intracellular delivery of bioactive cargo molecules. The mechanisms allowing CPPs to enter cells are ill-defined. Using a CRISPR/Cas9-based screening, we discovered that KCNQ5, KCNN4, and KCNK5 potassium channels positively modulate cationic CPP direct translocation into cells by decreasing the transmembrane potential (V_m). These findings provide the first unbiased genetic validation of the role of V_m in CPP translocation in cells. In silico modeling and live cell experiments indicate that CPPs, by bringing positive charges on the outer surface of the plasma membrane, decrease the V_m to very low values (–150 mV or less), a situation we have coined megapolarization that then triggers formation of water pores used by CPPs to enter cells. Megapolarization lowers the free energy barrier associated with CPP membrane translocation. Using dyes of varying dimensions in CPP co-entry experiments, the diameter of the water pores in living cells was estimated to be 2 (–5) nm, in accordance with the structural characteristics of the pores predicted by in silico modeling. Pharmacological manipulation to lower transmembrane potential boosted CPP cellular internalization in zebrafish and mouse models. Besides identifying the first proteins that regulate CPP translocation, this work characterized key mechanistic steps used by CPPs to cross cellular membranes. This opens the ground for strategies aimed at improving the ability of cells to capture CPP-linked cargos in vitro and in vivo.

Before a drug can have its desired effect, it must reach its target tissue or organ, and enter its cells. This is not easy because cells are surrounded by the plasma membrane, a fat-based barrier that separates the cell from its external environment. The plasma membrane contains proteins that act as channels, shuttling specific molecules in and out of the cell, and it also holds charge, with its inside surface being more negatively charged than its outside surface.

Cell-penetrating peptides are short sequences of amino acids (the building blocks that form proteins) that carry positive charges. These positive charges allow them to cross the membrane easily, but it is not well understood how.

To find out how cell-penetrating peptides cross the membrane, Trofimenko et al. attached them to dyes of different sizes. This revealed that the cell-penetrating peptides enter the cell through temporary holes called water pores, which measure about two nanometres across. The water pores form when the membrane becomes ‘megapolarized’, this is, when the difference in charge between the inside and the outside of the membrane becomes greater than normal. This can happen when the negative charge on the inside surface or the positive charge on the outer surface of the membrane increase. Megapolarization depends on potassium channels, which transport positive potassium ions outside the cell, making the outside of the membrane positive. When cell-penetrating peptides arrive at the outer surface of the cell near potassium channels, they make it even more positive. This increases the charge difference between the inside and the outside of the cell, allowing water pores to form. Once the peptides pass through the pores, the charge difference between the inside and the outside of the cell membrane dissipates, and the pores collapse.

Drug developers are experimenting with attaching cell-penetrating peptides to drugs to help them get inside their target cells. Currently there are several experimental medications of this kind in clinical trials. Understanding how these peptides gain entry, and what size of molecule they could carry with them, provides solid ground for further drug development.

Plant Sterol Clustering Correlates with Membrane Microdomains as Revealed by Optical and Computational Microscopy

Local inhomogeneities in lipid composition play a crucial role in the regulation of signal transduction and membrane traffic. This is particularly the case for plant plasma membrane, which is enriched in specific lipids, such as free and conjugated forms of phytosterols and typical phytosphingolipids. Nevertheless, most evidence for microdomains in cells remains indirect, and the nature of membrane inhomogeneities has been difficult to characterize. We used a new push–pull pyrene probe and fluorescence lifetime imaging microscopy (FLIM) combined with all-atom multiscale molecular dynamics simulations to provide a detailed view on the interaction between phospholipids and phytosterol and the effect of modulating cellular phytosterols on membrane-associated microdomains and phase separation formation. Our understanding of the organization principles of biomembranes is limited mainly by the challenge to measure distributions and interactions of lipids and proteins within the complex environment of living cells. Comparing phospholipids/phytosterol compositions typical of liquid-disordered (Ld) and liquid-ordered (Lo) domains, we furthermore show that phytosterols play crucial roles in membrane homeostasis. The simulation work highlights how state-of-the-art modeling alleviates some of the prior concerns and how unrefuted discoveries can be made through a computational microscope. Altogether, our results support the role of phytosterols in the lateral structuring of the PM of plant cells and suggest that they are key compounds for the formation of plant PM microdomains and the lipid-ordered phase.

Molecular Insights into Pore Formation Mechanism, Membrane Perturbation, and Water Permeation by the Antimicrobial Peptide Pleurocidin: A Combined All-Atom and Coarse-Grained Molecular Dynamics Simulation Study

The antimicrobial peptide (AMP) pleurocidin has a broad antimicrobial activity against Gram-negative and Gram-positive bacteria by perturbation and permeabilizing their membranes; however, understanding the mechanism of action of pleurocidin, a promising AMP for replacing current antibiotic agents, has tremendous importance for future applications. Hence, we applied all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations to provide molecular-level insights into the pore-forming process. The early stages of pore formation were examined by 500 ns AA simulations. The results demonstrated that pleurocidin has the ability to create a pore with two peptides through which water molecules can flow. However, the results of the 25 μs CG simulations indicate that the final pore will be created by accumulation of more than two peptides. The results show that after 2.5 μs of simulations, peptides will aggregate and create a channel-like pore across the membrane. Pleurocidin can construct a more efficient and stable pore in the anionic membranes than in the zwitterionic membranes. Moreover, the structure amphipathicity, polarity, and basic residues play crucial roles in the pore formation and flow of water molecules across the lipid bilayers. In general, the findings revealed that based on the lipid compositions of the membranes, pleurocidin could act by forming either toroidal or disordered toroidal pores with different peptide arrangements.

General Protocol for Constructing Molecular Models of Nanodiscs

Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.

Quantitative Estimation of Cyclotide-Induced Bilayer Membrane Disruption by Lipid Extraction with Mesoscopic Simulation

Cyclotide-induced membrane disruption is studied at the microsecond timescale by dissipative particle dynamics to quantitatively estimate a kinetic rate constant for membrane lipid extraction with a ″sandwich″ interaction model where two bilayer membranes enclose a cyclotide/water compartment. The obtained bioactivity trends for cyclotides Kalata B1, Cycloviolacin O2, and selected mutants with different membrane types are in agreement with experimental findings: For all membranes investigated, Cycloviolacin O2 shows a higher lipid extraction activity than Kalata B1. The presence of cholesterol leads to a decreased cyclotide activity compared to cholesterol-free membranes. Phosphoethanolamine-rich membranes exhibit an increased membrane disruption. A cyclotide's ″hydrophobic patch″ surface area is important for its bioactivity. A replacement of or with charged amino acid residues may lead to super-mutants with above-native activity but without simple charge-activity patterns. Cyclotide mixtures show linearly additive bioactivities without significant sub- or over-additive effects. The proposed method can be applied as a fast and easy-to-use tool for exploring structure-activity relationships of cyclotide/membrane systems: With the open software provided, the rate constant of a single cyclotide/membrane system can be determined in about 1 day by a scientific end-user without programming skills.

The antibiotic darobactin mimics a β-strand to inhibit outer membrane insertase

Antibiotics that target Gram-negative bacteria in new ways are needed to resolve the antimicrobial resistance crisis^1-3. Gram-negative bacteria are protected by an additional outer membrane, rendering proteins on the cell surface attractive drug targets^4,5. The natural compound darobactin targets the bacterial insertase BamA⁶-the central unit of the essential BAM complex, which facilitates the folding and insertion of outer membrane proteins^7-13. BamA lacks a typical catalytic centre, and it is not obvious how a small molecule such as darobactin might inhibit its function. Here we resolve the mode of action of darobactin at the atomic level using a combination of cryo-electron microscopy, X-ray crystallography, native mass spectrometry, in vivo experiments and molecular dynamics simulations. Two cyclizations pre-organize the darobactin peptide in a rigid β-strand conformation. This creates a mimic of the recognition signal of native substrates with a superior ability to bind to the lateral gate of BamA. Upon binding, darobactin replaces a lipid molecule from the lateral gate to use the membrane environment as an extended binding pocket. Because the interaction between darobactin and BamA is largely mediated by backbone contacts, it is particularly robust against potential resistance mutations. Our results identify the lateral gate as a functional hotspot in BamA and will allow the rational design of antibiotics that target this bacterial Achilles heel.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

Investigating the lipid fingerprint of SLC6 neurotransmitter transporters: a comparison of dDAT, hDAT, hSERT, and GlyT2

The local lipid annulus, or "fingerprint", of four SLC6 transporters (dDAT, hDAT, hSERT, and GlyT2) embedded in a complex neuronal membrane were compared and characterised using molecular dynamics. Our analysis included the development of new tools to improve membrane leaflet detection and the analysis of leaflet-dependent properties. Overall, the lipid fingerprints of the four transporters are comprised of similar lipids when grouped by headgroup or tail saturation. The enrichment and depletion of specific lipids, including sites of cholesterol contacts, varies between transporters. The subtle differences in lipid fingerprints results in varying membrane biophysical properties near the transporter. Our results highlight that the lipid-fingerprint of SLC6 transporters in complex membranes is highly dependent on membrane composition. Our results further characterize how the presence and identity of membrane proteins affects the complex interplay of lipid-protein interactions, influencing the local lipid environment and membrane biophysical properties.

Thermodynamic and kinetic stabilities of transmembrane helix bundles as revealed by single-pair FRET analysis: Effects of the number of membrane-spanning segments and cholesterol

The tertiary structures and conformational dynamics of transmembrane (TM) helical proteins are maintained by the interhelical interaction network in membranes, although it is complicated to analyze the underlying driving forces because the amino acid sequences can involve multiple and various types of interactions. To obtain insights into basal and common effects of the number of membrane-spanning segments and membrane cholesterol, we measured stabilities of helix bundles composed of simple TM helices (AALALAA)₃ (1TM) and (AALALAA)₃-G₅-(AALALAA)₃ (2TM). Association-dissociation dynamics for 1TM-1TM, 1TM-2TM, and 2TM-2TM pairs were monitored to compare stabilities of 2-, 3-, and 4-helical bundles, respectively, with single-pair fluorescence resonance energy transfer (sp-FRET) in liposome membranes. Both thermodynamic and kinetic stabilities of the helix bundles increased with a greater number of membrane-spanning segments in POPC. The presence of 30 mol% cholesterol strongly enhanced the formation of 1TM-1TM and 1TM-2TM bundles (~ - 9 kJ mol^-1), whereas it only weakly stabilized the 2TM-2TM bundle (~ - 3 kJ mol^-1). Fourier transform infrared-polarized attenuated total reflection (ATR-FTIR) spectroscopy revealed an ~30° tilt of the helix axis relative to bilayer normal for the 1TM-2TM pair in the presence of cholesterol, suggesting the formation of a tilted helix bundle to release high lateral pressure at the center of cholesterol-containing membranes. These results demonstrate that the number of membrane-spanning segments affects the stability and structure of the helix bundle, and their cholesterol-dependences. Such information is useful to understand the basics of folding and assembly of multispanning TM proteins.

Assessing the Perturbing Effects of Drugs on Lipid Bilayers Using Gramicidin Channel-Based In Silico and In Vitro Assays

Partitioning of bioactive molecules, including drugs, into cell membranes may produce indiscriminate changes in membrane protein function. As a guide to safe drug development, it therefore becomes important to be able to predict the bilayer-perturbing potency of hydrophobic/amphiphilic drugs candidates. Toward this end, we exploited gramicidin channels as molecular force probes and developed in silico and in vitro assays to measure drugs' bilayer-modifying potency. We examined eight drug-like molecules that were found to enhance or suppress gramicidin channel function in a thick 1,2-dierucoyl-sn-glycero-3-phosphocholine (DC22:1PC) but not in thin 1,2-dioleoyl-sn-glycero-3-phosphocholine (DC18:1PC) lipid bilayer. The mechanism underlying this difference was attributable to the changes in gramicidin dimerization free energy by drug-induced perturbations of lipid bilayer physical properties and bilayer-gramicidin interactions. The combined in silico and in vitro approaches, which allow for predicting the perturbing effects of drug candidates on membrane protein function, have implications for preclinical drug safety assessment.

Piezo1 Forms Specific, Functionally Important Interactions with Phosphoinositides and Cholesterol

Touch, hearing, and blood pressure regulation require mechanically gated ion channels that convert mechanical stimuli into electrical currents. One such channel is Piezo1, which plays a key role in the transduction of mechanical stimuli in humans and is implicated in diseases, such as xerocytosis and lymphatic dysplasia. There is building evidence that suggests Piezo1 can be regulated by the membrane environment, with the activity of the channel determined by the local concentration of lipids, such as cholesterol and phosphoinositides. To better understand the interaction of Piezo1 with its environment, we conduct simulations of the protein in a complex mammalian bilayer containing more than 60 different lipid types together with electrophysiology and mutagenesis experiments. We find that the protein alters its local membrane composition, enriching specific lipids and forming essential binding sites for phosphoinositides and cholesterol that are functionally relevant and often related to Piezo1-mediated pathologies. We also identify a number of key structural connections between the propeller and pore domains located close to lipid-binding sites.

Triacylglycerols sequester monotopic membrane proteins to lipid droplets

Triacylglycerols (TG) are synthesized at the endoplasmic reticulum (ER) bilayer and packaged into organelles called lipid droplets (LDs). LDs are covered by a single phospholipid monolayer contiguous with the ER bilayer. This connection is used by several monotopic integral membrane proteins, with hydrophobic membrane association domains (HDs), to diffuse between the organelles. However, how proteins partition between ER and LDs is not understood. Here, we employed synthetic model systems and found that HD-containing proteins strongly prefer monolayers and returning to the bilayer is unfavorable. This preference for monolayers is due to a higher affinity of HDs for TG over membrane phospholipids. Protein distribution is regulated by PC/PE ratio via alterations in monolayer packing and HD-TG interaction. Thus, HD-containing proteins appear to non-specifically accumulate to the LD surface. In cells, protein editing mechanisms at the ER membrane would be necessary to prevent unspecific relocation of HD-containing proteins to LDs.

Triacylglycerols (TG) are synthesized at the endoplasmic reticulum (ER) bilayer and packaged into monolayer lipid droplets (LDs), but how proteins partition between ER and LDs is poorly understood. Here authors use synthetic model systems and find that proteins containing hydrophobic membrane association domains strongly prefer monolayers and that returning to the bilayer is unfavorable.

Characterizing the Hydration Properties of Proton Binding Sites in the ATP Synthase c-Rings of Bacillus Species

The membrane-embedded domain of ATP synthases contains the c-ring, which translocates ions across the membrane, and its resultant rotation is coupled to ATP synthesis in the extramembranous domain. During rotation, the c-ring becomes accessible on both sides of the lipid bilayer to solvent via channels connected to the other membrane-embedded component, the a subunit, and thereby allows the ion to be released into the solvent environment. In recent times, many experimental structures of c-rings from different species have been solved. In some of these, a water molecule with a proposed "structural role" has been identified within the c-ring ion binding site, but in general, the requirement for high resolution to resolve specific water densities complicates their interpretation. In the present study, we use molecular dynamics (MD) simulations and rigorous free energy calculations to characterize the dynamics and energetics of a water molecule within the ion binding site of the c-ring from Bacillus pseudofirmus OF4, in its wild type (WT) and P51A mutant forms, along with the c-ring from thermophilic Bacillus PS3. Our data suggest that a water molecule stably binds to the P51A mutant, as well as helping to identify a bound water molecule in Bacillus PS3 whose presence was previously overlooked due to the limited resolution of the structural data. Sequence analysis further identifies a novel conserved sequence motif that is likely required to harbor a water molecule for stable ion coordination in the binding site of such proteins.

Lipid composition modulates ATP hydrolysis and calcium phosphate mineral propagation by TNAP-harboring proteoliposomes

Bone biomineralization is mediated by a special class of extracellular vesicles, named matrix vesicles (MVs), released by osteogenic cells. The MV membrane is enriched in sphingomyelin (SM), cholesterol (Chol) and tissue non-specific alkaline phosphatase (TNAP) compared with the parent cells' plasma membrane. TNAP is an ATP phosphohydrolase bound to cell and MV membranes via a glycosylphosphatidylinositol (GPI) anchor. Previous studies have shown that the lipid microenvironment influences the catalytic activity of enzymes incorporated into lipid bilayers. However, there is a lack of information about how the lipid microenvironment controls the ability of MV membrane-bound enzymes to induce mineral precipitation. Herein, we used TNAP-harboring proteoliposomes made of either pure dimyristoylphosphatidylcholine (DMPC) or DMPC mixed with either Chol, SM or both of them as MV biomimetic systems to evaluate how the composition modulates the lipid microenvironment and, in turn, TNAP incorporation into the lipid bilayer by means of calorimetry. These results were correlated with the proteoliposomes' catalytic activity and ability to induce the precipitation of amorphous calcium phosphate (ACP) in vitro. DMPC:SM proteoliposomes displayed the highest efficiency of mineral propagation, apparent affinity for ATP and substrate hydrolysis efficiency, which correlated with their highest degree of membrane organization (highest ΔH), among the tested proteoliposomes. Results obtained from turbidimetry and Fourier transformed infrared (FTIR) spectroscopy showed that the tested proteoliposomes induced ACP precipitation with the order DMPC:SM>DMPC:Chol:SM≈DMPC:Chol>DMPC which correlated with the lipid organization and the presence of SM in the proteoliposome membrane. Our study arises important insights regarding the physical properties and role of lipid organization in MV-mediated mineralization.

Interactions of a Bacterial RND Transporter with a Transmembrane Small Protein in a Lipid Environment

The small protein AcrZ in Escherichia coli interacts with the transmembrane portion of the multidrug efflux pump AcrB and increases resistance of the bacterium to a subset of the antibiotic substrates of that transporter. It is not clear how the physical association of the two proteins selectively changes activity of the pump for defined substrates. Here, we report cryo-EM structures of AcrB and the AcrBZ complex in lipid environments, and comparisons suggest that conformational changes occur in the drug-binding pocket as a result of AcrZ binding. Simulations indicate that cardiolipin preferentially interacts with the AcrBZ complex, due to increased contact surface, and we observe that chloramphenicol sensitivity of bacteria lacking AcrZ is exacerbated when combined with cardiolipin deficiency. Taken together, the data suggest that AcrZ and lipid cooperate to allosterically modulate AcrB activity. This mode of regulation by a small protein and lipid may occur for other membrane proteins.

Size-dependent aggregation of hydrophobic nanoparticles in lipid membranes

The aggregation of nanoparticles affects their reactivity, transport across biological membranes, uptake into cells, toxicity, and fate in the environment. In the case of membrane-embedded, hydrophobic nanoparticles the relationship between size and aggregation pattern is not well understood. Here, we explore this relationship for the case of spherically symmetrical nanoparticles using the MARTINI coarse-grained force field. We find that the free energy of dimerization depends strongly on nanoparticle size: the smallest molecules (mimicking C60 fullerene) aggregate only weakly, the largest ones form large three-dimensional aggregates causing major deformations in the host membrane, and the intermediate-sized particles show a tendency to form linear aggregates. Suppressing membrane undulations reduces very significantly aggregation, and substantially abolishes linear aggregation, suggesting a relationship between membrane curvature and aggregation geometry. At low concentration, when placed on membranes of variable curvature, the intermediate size nanoparticles move rapidly to high curvature regions - suggesting that they can sense membrane curvature. At high concentration, the same nanoparticles induce massive membrane deformations, without affecting the mechanical stability of the membrane - suggesting that they can generate membrane curvature.

Lipid-dependent conformational landscape of the ErbB2 growth factor receptor dimers

Altered lipid metabolism has been linked to cancer development and progression. Several roles have been attributed to the increased saturation and length of lipid acyl tails observed in tumors, but its effect on signaling receptors is still emerging. In this work, we have analyzed the lipid dependence of the ErbB2 growth factor receptor dimerization that plays an important role in the pathogenesis of breast cancer. We have performed coarse-grain ensemble molecular dynamics simulations to comprehensively sample the ErbB2 monomer-dimer association. Our results indicate a dynamic dimer state with a complex conformational landscape that is modulated with increasing lipid tail length. We resolve the native N-terminal "active" and C-terminal "inactive" conformations in all membrane compositions. However, the relative population of the N-terminal and C-terminal conformers is dependent on length of the saturated lipid tails. In short-tail membranes, additional non-specific dimers are observed which are reduced or absent in long-tailed bilayers. Our results indicate that the relative population as well as the structure of the dimer state is modulated by membrane composition. We have correlated these differences to local perturbations of the membrane around the receptor. Our work is an important step in characterizing ErbB dimers in healthy and diseased states and emphasize the importance of sampling lipid dynamics in understanding receptor association.

The aggregation of striped nanoparticles in mixed phospholipid bilayers

The unique and adjustable properties of nanoparticles present enormous opportunities for their use as targeted drug delivery vectors. For example, nanoparticles functionalized with key surface ligands have been shown to pass through phospholipid bilayers without causing localised disruption. However, the further effects nanoparticles have on multi-component phospholipid bilayers remain unclear. We use coarse-grained computational models to investigate the structural properties of mixed phospholipid bilayers in the presence of ligand-functionalized nanoparticles. Model bilayers are composed of DPPC, DUPC, DFPC and cholesterol, and the nanoparticles are striped with a hydrophobic-ligand band and charged-ligand spherical caps. Our results show that nanoparticles aggregate near unsaturated phospholipid regions, phospholipid bilayer phase-separation is promoted in the presence of nanoparticles, and the heterogeneous components of a phospholipid bilayer play a significant role in the lateral organization of nanoparticles. This study highlights the need for considering the complexity of realistic phospholipid bilayers when optimising ligand functionalized nanoparticles for efficient drug delivery vectors.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Atomistic and coarse-grained simulations of membrane proteins: A practical guide

Membrane proteins are amphipathic macromolecules whose exposed hydrophobic surfaces promote interactions with lipid membranes. Membrane proteins are remarkably diverse in terms of chemical composition and correspondingly, their biological functions and general biophysical behavior. Conventional experimental techniques provide an approach to study specific properties of membrane proteins e.g. their surface features, the nature and abundance of stabilizing intramolecular forces, preferred bilayer orientation, and the characteristics of their annular lipid shells. Molecular modeling software-and in particular, the suite of molecular dynamics algorithms-enables a more comprehensive exploration of dynamic membrane protein behavior. Molecular dynamics methods enable users to produce stepwise trajectories of proteins on arbitrary spatiotemporal scales that enable the easy identification of dynamic interactions that are beyond the scope of conventional analytical techniques. This article explains the molecular dynamics theoretical framework and popular step-by-step approaches for simulating membrane proteins in planar, and to a lesser extent, nonplanar lipid geometries. We detail popular procedures and computational tools that produce well-packed configurations of lipids and proteins and additionally, the efficient molecular dynamics simulation algorithms that reproduce their dynamic interactions.

Computational Modeling of Claudin Structure and Function

Tight junctions form a barrier to control passive transport of ions and small molecules across epithelia and endothelia. In addition to forming a barrier, some of claudins control transport properties of tight junctions by forming charge- and size-selective ion channels. It has been suggested claudin monomers can form or incorporate into tight junction strands to form channels. Resolving the crystallographic structure of several claudins in recent years has provided an opportunity to examine structural basis of claudins in tight junctions. Computational and theoretical modeling relying on atomic description of the pore have contributed significantly to our understanding of claudin pores and paracellular transport. In this paper, we review recent computational and mathematical modeling of claudin barrier function. We focus on dynamic modeling of global epithelial barrier function as a function of claudin pores and molecular dynamics studies of claudins leading to a functional model of claudin channels.

Computational Nanoscopy of Tight Junctions at the Blood-Brain Barrier Interface

The selectivity of the blood-brain barrier (BBB) is primarily maintained by tight junctions (TJs), which act as gatekeepers of the paracellular space by blocking blood-borne toxins, drugs, and pathogens from entering the brain. The BBB presents a significant challenge in designing neurotherapeutics, so a comprehensive understanding of the TJ architecture can aid in the design of novel therapeutics. Unraveling the intricacies of TJs with conventional experimental techniques alone is challenging, but recently developed computational tools can provide a valuable molecular-level understanding of TJ architecture. We employed the computational methods toolkit to investigate claudin-5, a highly expressed TJ protein at the BBB interface. Our approach started with the prediction of claudin-5 structure, evaluation of stable dimer conformations and nanoscale assemblies, followed by the impact of lipid environments, and posttranslational modifications on these claudin-5 assemblies. These led to the study of TJ pores and barriers and finally understanding of ion and small molecule transport through the TJs. Some of these in silico, molecular-level findings, will need to be corroborated by future experiments. The resulting understanding can be advantageous towards the eventual goal of drug delivery across the BBB. This review provides key insights gleaned from a series of state-of-the-art nanoscale simulations (or computational nanoscopy studies) performed on the TJ architecture.

Discovery of Novel Multi-target Inhibitor of angiotensin type 1 receptor and neprilysin inhibitors from Traditional Chinese Medicine

Angiotensin II type-1 receptor–neprilysin inhibitor (ARNi) is consisted of Angiotensin II type-1 receptor (AT1) antagonist and neprilysin (NEP) inhibitor, which could simultaneously increase the vasodilators of the natriuretic peptides and antagonize vasoconstrictors of Ang II. ARNi has been proved a superior effect and lower risks of death on chronic heart failure (CHF) and hypertension. In this paper, ARNi from Traditional Chinese Medicines (TCM) was discovered based on target combination of AT1 and NEP by virtual screening, biological assay and molecular dynamics (MD) simulations. Two customized strategies of combinatorial virtual screening were implemented to discover AT1 antagonist and NEP inhibitor based on pharmacophore modeling and docking computation respectively. Gyrophoric acid (PubChem CID: 135728) from Parmelia saxatilis was selected as AT1 antagonist and assayed with IC₅₀ of 29.76 μM by calcium influx assay. And 3,5,3′-triiodothyronine (PubChem CID: 861) from Bos taurus domesticus was screened as NEP inhibitor and has a dose dependent inhibitory activity by biochemistry fluorescence assay. Combined with MD simulations, these compounds can generate interaction with the target, key interactive residues of ARG167, TRP84, and VAL108 in AT1, and HIS711 in NEP were also identified respectively. This study designs the combinatorial strategy to discover novel frames of ARNi from TCM, and gyrophoric acid and 3,5,3′-triiodothyronine could provide the clues and revelations of drug design and therapeutic method of CHF and hypertension for TCM clinical applications.

Molecular dynamics study of membrane permeabilization by wild-type and mutant lytic peptides from the non-enveloped Flock House virus

Flock House virus (FHV) serves as a model system for understanding infection mechanisms utilized by non-enveloped viruses to transport across cellular membranes. During the infection cycle of FHV, a fundamental stage involves disruption of the endosomal membrane by membrane active peptides, following externalization of the peptides from the capsid interior. The FHV lytic agents are the 44 C-terminal amino acids residues of the capsid protein, which are auto-catalytically cleaved during the capsid maturation process. The cleaved peptides are termed γ peptides. In this study, we perform multi-scale molecular dynamics simulations including 40 μs all-atom molecular dynamics simulations to study the behavior of pre-inserted transmembrane lytic peptides at a high concentration in a neutral membrane. We study the dynamical organization among peptides to form oligomeric bundles in four systems including the wild-type γ peptide and three mutant forms; namely, a truncation mutant in which the 23 C-terminal residues are deleted (γ₁), a construct where the 8 C-terminal residues of γ are fused to γ₁ (Δ385-399 γ) and a single-point mutant (F402A γ), all of which have been experimentally shown to drastically affect infectivity and lytic activity compared to the wild-type γ. Our results shed light on the actions of varied forms of the FHV lytic peptide including membrane insertion, trans-membrane stability, peptide oligomerization, water permeation activity and dynamic pore formation. Findings from this study provide detailed structural information and rationale for the differences in lytic activity among variants of FHV γ.

Membrane Remodeling by the Lytic Fragment of SticholysinII: Implications for the Toroidal Pore Model

Sticholysins are pore-forming toxins of biomedical interest and represent a prototype of proteins acting through the formation of protein-lipid or toroidal pores. Peptides spanning the N-terminus of sticholysins can mimic their permeabilizing activity and, together with the full-length toxins, have been used as a tool to understand the mechanism of pore formation in membranes. However, the lytic mechanism of these peptides and the lipid shape modulating their activity are not completely clear. In this article, we combine molecular dynamics simulations and experimental biophysical tools to dissect different aspects of the pore-forming mechanism of StII_1-30, a peptide derived from the N-terminus of sticholysin II (StII). With this combined approach, membrane curvature induction and flip-flop movement of the lipids were identified as two important membrane remodeling steps mediated by StII_1-30. Pore formation by this peptide was enhanced by the presence of the negatively curved lipid phosphatidylethanolamine in membranes. This lipid emerged not only as a facilitator of membrane interactions but also as a structural element of the StII_1-30 pore that is recruited to the ring upon its assembly. Collectively, these, to our knowledge, new findings support a toroidal model for the architecture of the pore formed by StII_1-30 and provide new molecular insight into the role of phosphatidylethanolamine as a membrane component that can easily integrate into the ring of toroidal pores, thus probably aiding in their stabilization. This study contributes to a better understanding of the molecular mechanism underlying the permeabilizing activity of StII_1-30 and peptides or proteins acting via a toroidal pore mechanism and offers an informative framework for the optimization of the biomedical application of this and similar molecules.

Cholic Acid-Derived Amphiphile which Combats Gram-Positive Bacteria-Mediated Infections via Disintegration of Lipid Clusters

Inappropriate and uncontrolled use of antibiotics results in the emergence of antibiotic resistance, thereby threatening the present clinical regimens to treat infectious diseases. Therefore, new antimicrobial agents that can prevent bacteria from developing drug resistance are urgently needed. Selective disruption of bacterial membranes is the most effective strategy for combating microbial infections as accumulation of genetic mutations will not allow for the emergence of drug resistance against these antimicrobials. In this work, we tested cholic acid (CA) derived amphiphiles tethered with different alkyl chains for their ability to combat Gram-positive bacterial infections. In-depth biophysical and biomolecular simulation studies suggested that the amphiphile with a hexyl chain (6) executes more effective interactions with Gram-positive bacterial membranes as compared to other hydrophobic counterparts. Amphiphile 6 is effective against multidrug resistant Gram-positive bacterial strains as well and does not allow the adherence of S. aureus on amphiphile 6 coated catheters implanted in mice. Further, treatment of wound infections with amphiphile 6 clears the bacterial infections. Therefore, the current study presents strategic guidelines in design and development of CA-derived membrane-targeting antimicrobials for Gram-positive bacterial infections.

Computational Modeling of Realistic Cell Membranes

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

Outer Membrane Proteins OmpA, FhuA, OmpF, EstA, BtuB, and OmpX Have Unique Lipopolysaccharide Fingerprints

The outer membrane of Gram-negative bacteria has a highly complex asymmetrical architecture, containing a mixture of phospholipids in the inner leaflet and almost exclusively lipopolysaccharide (LPS) molecules in the outer leaflet. In E. coli, the outer membrane contains a wide range of proteins with a β barrel architecture, that vary in size from the smallest having eight strands to larger barrels composed of 22 strands. Here we report coarse-grained molecular dynamics simulations of six proteins from the E. coli outer membrane OmpA, OmpX, BtuB, FhuA, OmpF, and EstA in a range of membrane environments, which are representative of the in vivo conditions for different strains of E. coli. We show that each protein has a unique pattern of interaction with the surrounding membrane, which is influenced by the composition of the protein, the level of LPS in the outer leaflet, and the differing mobilities of the lipids in the two leaflets of the membrane. Overall we present analyses from over 200 μs of simulation for each protein.

The disordered plant dehydrin Lti30 protects the membrane during water-related stress by cross-linking lipids

Dehydrins are intrinsically disordered proteins, generally expressed in plants as a response to embryogenesis and water-related stress. Their suggested functions are in membrane stabilization and cell protection. All dehydrins contain at least one copy of the highly conserved K-segment, proposed to be a membrane-binding motif. The dehydrin Lti30 (Arabidopsis thaliana) is up-regulated during cold and drought stress conditions and comprises six K-segments, each with two adjacent histidines. Lti30 interacts with the membrane electrostatically via pH-dependent protonation of the histidines. In this work, we seek a molecular understanding of the membrane interaction mechanism of Lti30 by determining the diffusion and molecular organization of Lti30 on model membrane systems by imaging total internal reflection- fluorescence correlation spectroscopy (ITIR-FCS) and molecular dynamics (MD) simulations. The dependence of the diffusion coefficient explored by ITIR-FCS together with MD simulations yields insights into Lti30 binding, domain partitioning, and aggregation. The effect of Lti30 on membrane lipid diffusion was studied on fluorescently labeled supported lipid bilayers of different lipid compositions at mechanistically important pH conditions. In parallel, we compared the mode of diffusion for short individual K-segment peptides. The results indicate that Lti30 binds the lipid bilayer via electrostatics, which restricts the mobility of lipids and bound protein molecules. At low pH, Lti30 binding induced lipid microdomain formation as well as protein aggregation, which could be correlated with one another. Moreover, at physiological pH, Lti30 forms nanoscale aggregates when proximal to the membrane suggesting that Lti30 may protect the cell by "cross-linking" the membrane lipids.