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

Molecular dynamics simulations of bacterial outer membrane lipid extraction: Adequate sampling?

The outer membrane of Gram-negative bacteria is almost exclusively composed of lipopolysaccharide in its outer leaflet, whereas the inner leaflet contains a mixture of phospholipids. Lipopolysaccharide diffuses at least an order of magnitude slower than phospholipids, which can cause issues for molecular dynamics simulations in terms of adequate sampling. Here, we test a number of simulation protocols for their ability to achieve convergence with reasonable computational effort using the MARTINI coarse-grained force-field. This is tested in the context both of potential of mean force (PMF) calculations for lipid extraction from membranes and of lateral mixing within the membrane phase. We find that decoupling the cations that cross-link the lipopolysaccharide headgroups from the extracted lipid during PMF calculations is the best approach to achieve convergence comparable to that for phospholipid extraction. We also show that lateral lipopolysaccharide mixing/sorting is very slow and not readily addressable even with Hamiltonian replica exchange. We discuss why more sorting may be unrealistic for the short (microseconds) timescales we simulate and provide an outlook for future studies of lipopolysaccharide-containing membranes.

CHARMM-GUI supports the Amber force fields

As part of our ongoing efforts to support diverse force fields and simulation programs in CHARMM-GUI, this work presents the development of FF-Converter to prepare Amber simulation inputs with various Amber force fields within the current CHARMM-GUI workflow. The currently supported Amber force fields are ff14SB/ff19SB (protein), Bsc1 (DNA), OL3 (RNA), GLYCAM06 (carbohydrate), Lipid17 (lipid), GAFF/GAFF2 (small molecule), TIP3P/TIP4P-EW/OPC (water), and 12-6-4 ions, and more will be added if necessary. The robustness and usefulness of this new CHARMM-GUI extension are demonstrated by two exemplary systems: a protein/N-glycan/ligand/membrane system and a protein/DNA/RNA system. Currently, CHARMM-GUI supports the Amber force fields only for the Amber program, but we will expand the FF-Converter functionality to support other simulation programs that support the Amber force fields.

Molecular architecture of the endocytic TPLATE complex

Eukaryotic cells rely on endocytosis to regulate their plasma membrane proteome and lipidome. Most eukaryotic groups, except fungi and animals, have retained the evolutionary ancient TSET complex as an endocytic regulator. Unlike other coatomer complexes, structural insight into TSET is lacking. Here, we reveal the molecular architecture of plant TSET [TPLATE complex (TPC)] using an integrative structural approach. We identify crucial roles for specific TSET subunits in complex assembly and membrane interaction. Our data therefore generate fresh insight into the differences between the hexameric TSET in Dictyostelium and the octameric TPC in plants. Structural elucidation of this ancient adaptor complex represents the missing piece in the coatomer puzzle and vastly advances our functional as well as evolutionary insight into the process of endocytosis.

IFITM3 functions as a PIP3 scaffold to amplify PI3K signalling in B cells

Interferon-induced transmembrane protein 3 (IFITM3) has previously been identified as an endosomal protein that blocks viral infection1-3. Here we studied clinical cohorts of patients with B cell leukaemia and lymphoma, and identified IFITM3 as a strong predictor of poor outcome. In normal resting B cells, IFITM3 was minimally expressed and mainly localized in endosomes. However, engagement of the B cell receptor (BCR) induced both expression of IFITM3 and phosphorylation of this protein at Tyr20, which resulted in the accumulation of IFITM3 at the cell surface. In B cell leukaemia, oncogenic kinases phosphorylate IFITM3 at Tyr20, which causes constitutive localization of this protein at the plasma membrane. In a mouse model, Ifitm3-/- naive B cells developed in normal numbers; however, the formation of germinal centres and the production of antigen-specific antibodies were compromised. Oncogenes that induce the development of leukaemia and lymphoma did not transform Ifitm3-/- B cells. Conversely, the phosphomimetic IFITM3(Y20E) mutant induced oncogenic PI3K signalling and initiated the transformation of premalignant B cells. Mechanistic experiments revealed that IFITM3 functions as a PIP3 scaffold and central amplifier of PI3K signalling. The amplification of PI3K signals depends on IFITM3 using two lysine residues (Lys83 and Lys104) in its conserved intracellular loop as a scaffold for the accumulation of PIP3. In Ifitm3-/- B cells, lipid rafts were depleted of PIP3, which resulted in the defective expression of over 60 lipid-raft-associated surface receptors, and impaired BCR signalling and cellular adhesion. We conclude that the phosphorylation of IFITM3 that occurs after B cells encounter antigen induces a dynamic switch from antiviral effector functions in endosomes to a PI3K amplification loop at the cell surface. IFITM3-dependent amplification of PI3K signalling, which in part acts downstream of the BCR, is critical for the rapid expansion of B cells with high affinity to antigen. In addition, multiple oncogenes depend on IFITM3 to assemble PIP3-dependent signalling complexes and amplify PI3K signalling for malignant transformation.

Ultrahigh Adhesion Force Between Silica-Binding Peptide SB7 and Glass Substrate Studied by Single-Molecule Force Spectroscopy and Molecular Dynamic Simulation

Many proteins and peptides have been identified to effectively and specifically bind on certain surfaces such as silica, polystyrene and titanium dioxide. It is of great interest, in many areas such as enzyme immobilization, surface functionalization and nanotechnology, to understand how these proteins/peptides bind to solid surfaces. Here we use single-molecule force spectroscopy (SMFS) based on atomic force microscopy to directly measure the adhesion force between a silica-binding peptide SB7 and glass surface at single molecule level. SMFS results show that the adhesion force of a single SB7 detaching from the glass surface distributes in two populations at ~220 pN and 610 pN, which is higher than the unfolding forces of most mechanically stable proteins and the unbinding forces of most stable protein-protein interactions. Molecular dynamics simulation reveals that the electrostatic interactions between positively charged arginine residues and the silica surface dominates the binding of SB7 on silica. Our study provides experimental evidence and molecular mechanism at the single-molecule level for the SB7-based immobilization of proteins on silica-based surface, which is able to withstand high mechanical forces, making it an ideal fusion tag for silica surface immobilization or peptide-base adhesive materials.

Double-Headed Cationic Lipopeptides: An Emerging Class of Antimicrobials

Antimicrobial peptides (AMPs) constitute a promising tool in the development of novel therapeutic agents useful in a wide range of bacterial and fungal infections. Among the modifications improving pharmacokinetic and pharmacodynamic characteristics of natural AMPs, an important role is played by lipidation. This study focuses on the newly designed and synthesized lipopeptides containing multiple Lys residues or their shorter homologues with palmitic acid (C₁₆) attached to the side chain of a residue located in the center of the peptide sequence. The approach resulted in the development of lipopeptides representing a model of surfactants with two polar headgroups. The aim of this study is to explain how variations in the length of the peptide chain or the hydrocarbon side chain of an amino acid residue modified with C₁₆, affect biological functions of lipopeptides, their self-assembling propensity, and their mode of action.

Oxidation degree of a cell membrane model and its response to structural changes, a coarse-grained molecular dynamics approach

Oxidative stress plays an essential role in the regulation of vital processes in living organisms. Reactive oxygen species can react chemically with the constituents of the cells leading to irreversible damage. The first structure of the cell in contact with the environment that surrounds it is the membrane, which protects it and allows the exchange of substances. Some signals manifest when the components of a bilayer are undergoing oxidation, like an increase in the lipid area, decrease in the thickness of the bilayer, and exchange of the oxidized groups toward the bilayer surface. In this investigation, a molecular dynamics simulation was done on a set of Dioleoylphosphatidylcholine membranes with different percentage of oxidized lipids, in order to observe the effect of the oxidation degree on the membrane structure. It was found that, as higher the concentration of oxidized lipids is, the larger the damage of the membrane. This is reflected in the increase in the lipid area and the decrease in the thickness and membrane packing. Also, it was observed that hydrophobicity inside the membrane decreases as the oxidation percentage increases.Communicated by Ramaswamy H. Sarma.

Engineering REST-Specific Synthetic PUF Proteins to Control Neuronal Gene Expression: A Combined Experimental and Computational Study

Regulation of gene transcription is an essential mechanism for differentiation and adaptation of organisms. A key actor in this regulation process is the repressor element 1 (RE1)-silencing transcription factor (REST), a transcriptional repressor that controls more than 2000 putative target genes, most of which are neuron-specific. With the purpose of modulating REST expression, we exploited synthetic, ad hoc designed, RNA binding proteins (RBPs) able to specifically target and dock to REST mRNA. Among the various families of RBPs, we focused on the Pumilio and FBF (PUF) proteins, present in all eukaryotic organisms and controlling a variety of cellular functions. Here, a combined experimental and computational approach was used to design and test 8- and 16-repeat PUF proteins specific for REST mRNA. We explored the conformational properties and atomic features of the PUF-RNA recognition code by Molecular Dynamics simulations. Biochemical assays revealed that the 8- and 16-repeat PUF-based variants specifically bind the endogenous REST mRNA without affecting its translational regulation. The data also indicate a key role of stacking residues in determining the binding specificity. The newly characterized REST-specific PUF-based constructs act as excellent RNA-binding modules and represent a versatile and functional platform to specifically target REST mRNA and modulate its endogenous expression.

Assessing Barriers for Antimicrobial Penetration in Complex Asymmetric Bacterial Membranes: A Case Study with Thymol

The bacterial cell envelope is a complex multilayered structure evolved to protect bacteria in hostile environments. An understanding of the molecular basis for the interaction and transport of antibacterial therapeutics with the bacterial cell envelope will enable the development of drug molecules to combat bacterial infections and suppress the emergence of drug-resistant strains. Here we report the successful creation of an in vitro supported lipid bilayer (SLB) platform of the outer membrane (OM) of E. coli, an archetypical Gram-negative bacterium, containing the full smooth lipopolysaccharide (S-LPS) architecture of the membrane. Using this platform, we performed fluorescence correlation spectroscopy (FCS) in combination with molecular dynamics (MD) simulations to measure lipid diffusivities and provide molecular insights into the transport of natural antimicrobial agent thymol. Lipid diffusivities measured on symmetric supported lipid bilayers made up of inner membrane lipids show a distinct increase in the presence of thymol as also corroborated by MD simulations. However, lipid diffusivities in the asymmetric OM consisting of only S-LPS are invariant upon exposure to thymol. Increasing the phospholipid content in the LPS-containing outer leaflet improved the penetration toward thymol as reflected in slightly higher relative diffusivity changes in the inner leaflet when compared with the outer leaflet. Free-energy computations reveal the presence of a barrier (∼6 kT) only in the core-saccharide region of the OM for the translocation of thymol while the external O-antigen part is easily traversed. In contrast, thymol spontaneously inserts into the inner membrane. In addition to providing leaflet-resolved penetration barriers in bacterial membranes, we also assess the ability of small molecules to penetrate various membrane components. With rising bacterial resistance, our study opens up the possibility of screening potential antimicrobial drug candidates using these realistic model platforms for Gram-negative bacteria.

Insights into the mechanism of membrane fusion induced by the plant defense element, plant-specific insert

In plants, many natural defense mechanisms include cellular membrane fusion as a way to resist infection by external pathogens. Several plant proteins mediate membrane fusion, but the detailed mechanism by which they promote fusion is less clear. Understanding this process could provide valuable insights into these proteins' physiological functions and guide bioengineering applications (i.e. the design of antimicrobial proteins). The plant-specific insert (PSI) from Solanum tuberosum can help reduce certain pathogen attack via membrane fusion. To gain new insights into the process of PSI-induced membrane fusion, a combined approach of NMR, FRET, and in silico studies was used. Our results indicate that (i) under acidic conditions, the PSI experiences a monomer-dimer equilibrium, and the dimeric PSI induces membrane fusion below a certain critical pH; (ii) after fusion, the PSI resides in a highly dehydrated environment with limited solvent accessibility, suggesting its capability in reducing repulsive dehydration forces between liposomes to facilitate fusion; and (iii) as shown by molecular dynamics simulations, the PSI dimer can bind stably to membrane surfaces and can bridge liposomes in close proximity, a critical step for the membrane fusion. In summary, this study provides new and unique insights into the mechanisms by which the PSI and similar proteins induce membrane fusion.

Evidence for the Supramolecular Organization of a Bacterial Outer-Membrane Protein from In Vivo Pulse Electron Paramagnetic Resonance Spectroscopy

In the outer membrane of Gram-negative bacteria, membrane proteins are thought to be organized into domains or islands that play a role in the segregation, movement, and turnover of membrane components. However, there is presently limited information on the structure of these domains or the molecular interactions that mediate domain formation. In the present work, the Escherichia coli outer membrane vitamin B₁₂ transporter, BtuB, was spin-labeled, and double electron-electron resonance was used to measure the distances between proteins in intact cells. These data together with Monte Carlo simulations provide evidence for the presence of specific intermolecular contacts between BtuB monomers that could drive the formation of string-like oligomers. Moreover, the EPR data provide evidence for the location of the interacting interfaces and indicate that lipopolysaccharide mediates the contacts between BtuB monomers.

Dual Resolution Membrane Simulations Using Virtual Sites

All-atomistic (AA) and coarse-grain (CG) simulations have been successfully applied to investigate a broad range of biomolecular processes. However, the accessible time and length scales of AA simulation are limited and the specific molecular details of CG simulation are simplified. Here, we propose a virtual site (VS) based hybrid scheme that can concurrently couple AA and CG resolutions in a single membrane simulation, mitigating the shortcomings of either representation. With some adjustments to make the AA and CG force fields compatible, we demonstrate that lipid bilayer properties are well kept in our hybrid approach. Our VS hybrid method was also applied to simulate a small lipid vesicle, with the inner leaflet and interior solvent represented in AA, and the outer leaflet together with exterior solvent at the CG level. Our multiscale method opens the way to investigate biomembrane properties at increased computational efficiency, in particular applications involving large solvent filled regions.

Molecular Simulations of Gram-Negative Bacterial Membranes Come of Age

Gram-negative bacteria are protected by a multicompartmental molecular architecture known as the cell envelope that contains two membranes and a thin cell wall. As the cell envelope controls influx and efflux of molecular species, in recent years both experimental and computational studies of such architectures have seen a resurgence due to the implications for antibiotic development. In this article we review recent progress in molecular simulations of bacterial membranes. We show that enormous progress has been made in terms of the lipidic and protein compositions of bacterial systems. The simulations have moved away from the traditional setup of one protein surrounded by a large patch of the same lipid type toward a more bio-logically representative viewpoint. Simulations with multiple cell envelope components are also emerging. We review some of the key method developments that have facilitated recent progress, discuss some current limitations, and offer a perspective on future directions.

Characterizing Membrane Association and Periplasmic Transfer of Bacterial Lipoproteins through Molecular Dynamics Simulations

Escherichia coli lipoprotein precursors at the inner membrane undergo three maturation stages before transport by the Lol system to the outer membrane. Here, we develop a pipeline to simulate the membrane association of bacterial lipoproteins in their four maturation states. This has enabled us to model and simulate 81 of the predicted 114 E. coli lipoproteins and reveal their interactions with the host lipid membrane. As part of this set we characterize the membrane contacts of LolB, the lipoprotein involved in periplasmic translocation. We also consider the means and bioenergetics for lipoprotein localization. Our calculations uncover a preference for LolB over LolA and therefore indicate how a lipoprotein may be favorably transferred from the inner to outer membrane. Finally, we reveal that LolC has a role in membrane destabilization, thereby promoting lipoprotein transfer to LolA.

Nanodiscs as a New Tool to Examine Lipid-Protein Interactions

The interactions between lipids and proteins are one of the most fundamental processes in living organisms, responsible for critical cellular events ranging from replication, cell division, signaling, and movement. Enabling the central coupling responsible for maintaining the functionality of the breadth of proteins, receptors, and enzymes that find their natural home in biological membranes, the fundamental mechanisms of recognition of protein for lipid, and vice versa, have been a focal point of biochemical and biophysical investigations for many decades. Complexes of lipids and proteins, such as the various lipoprotein factions, play central roles in the trafficking of important proteins, small molecules and metabolites and are often implicated in disease states. Recently an engineered lipoprotein particle, termed the nanodisc, a modified form of the human high density lipoprotein fraction, has served as a membrane mimetic for the investigation of membrane proteins and studies of lipid-protein interactions. In this review, we summarize the current knowledge regarding this self-assembling lipid-protein complex and provide examples for its utility in the investigation of a large number of biological systems.

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.

Lipopolysaccharide Simulations Are Sensitive to Phosphate Charge and Ion Parameterization

The high proportion of lipopolysaccharide (LPS) molecules in the outer membrane of Gram-negative bacteria makes it a highly effective barrier to small molecules, antibiotic drugs, and other antimicrobial agents. Given this vital role in protecting bacteria from potentially hostile environments, simulations of LPS bilayers and outer membrane systems represent a critical tool for understanding the mechanisms of bacterial resistance and the development of new antibiotic compounds that circumvent these defenses. The basis of these simulations is parameterizations of LPS, which have been developed for all major molecular dynamics force fields. However, these parameterizations differ in both the protonation state of LPS and how LPS membranes behave in the presence of various ion species. To address these discrepancies and understand the effects of phosphate charge on bilayer properties, simulations were performed for multiple distinct LPS chemotypes with different ion parameterizations in both protonated or deprotonated lipid A states. These simulations show that bilayer properties, such as the area per lipid and inter-lipid hydrogen bonding, are highly influenced by the choice of phosphate group charges, cation type, and ion parameterization, with protonated LPS and monovalent cations with modified nonbonded parameters providing the best match to the experiments. Additionally, alchemical free energy simulations were performed to determine theoretical pK_a values for LPS and subsequently validated by ³¹P solid-state nuclear magnetic resonance experiments. Results from these complementary computational and experimental studies demonstrate that the protonated state dominates at physiological pH, contrary to the deprotonated form modeled by many LPS force fields. Overall, these results highlight the sensitivity of LPS simulations to phosphate charge and ion parameters while offering recommendations for how existing models should be updated for consistency between force fields as well as to best match experiments.

Aggregation of Lipid A Variants: A Hybrid Particle-Field Model

Lipid A is one of the three components of bacterial lipopolysaccharides constituting the outer membrane of Gram-negative bacteria, and is recognized to have an important biological role in the inflammatory response of mammalians. Its biological activity is modulated by the number of acyl-chains that are present in the lipid and by the dielectric medium, i.e., the type of counter-ions, through electrostatic interactions. In this paper, we report on a coarse-grained model of chemical variants of Lipid A based on the hybrid particle-field/molecular dynamics approach (hPF-MD). In particular, we investigate the stability of Lipid A bilayers for two different hexa- and tetra-acylated structures. Comparing particle density profiles along bilayer cross-sections, we find good agreement between the hPF-MD model and reference all-atom simulation for both chemical variants of Lipid A. hPF-MD models of constituted bilayers composed by hexa-acylated Lipid A in water are stable within the simulation time. We further validate our model by verifying that the phase behavior of Lipid A/counterion/water mixtures is correctly reproduced. In particular, hPF-MD simulations predict the correct self-assembly of different lamellar and micellar phases from an initially random distribution of Lipid A molecules with counterions in water. Finally, it is possible to observe the spontaneous formation and stability of Lipid A vesicles by fusion of micellar aggregates.

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.

Vanillin enhances the passive transport rate and absorption of drugs with moderate oral bioavailability in vitro and in vivo by affecting the membrane structure

Vanillin is a popular flavoring agent in the food, tobacco, and perfume industries.

The Contribution of the Ankyrin Repeat Domain of TRPV1 as a Thermal Module

The TRPV1 cation nonselective ion channel plays an essential role in thermosensation and perception of other noxious stimuli. TRPV1 can be activated by low extracellular pH, high temperature, or naturally occurring pungent molecules such as allicin, capsaicin, or resiniferatoxin. Its noxious thermal sensitivity makes it an important participant as a thermal sensor in mammals. However, details of the mechanism of channel activation by increases in temperature remain unclear. Here, we used a combination of approaches to try to understand the role of the ankyrin repeat domain (ARD) in channel behavior. First, a computational modeling approach by coarse-grained molecular dynamics simulation of the whole TRPV1 embedded in a phosphatidylcholine and phosphatidylethanolamine membrane provides insight into the dynamics of this channel domain. Global analysis of the structural ensemble shows that the ARD is a region that sustains high fluctuations during dynamics at different temperatures. We then performed biochemical and thermal stability studies of the purified ARD by the means of circular dichroism and tryptophan fluorescence and demonstrate that this region undergoes structural changes at similar temperatures that lead to TRPV1 activation. Our data suggest that the ARD is a dynamic module and that it may participate in controlling the temperature sensitivity of TRPV1.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

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.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

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.

Binding from Both Sides: TolR and Full-Length OmpA Bind and Maintain the Local Structure of the E. coli Cell Wall

We present a molecular modeling and simulation study of the E. coli cell envelope, with a particular focus on the role of TolR, a native protein of the E. coli inner membrane, in interactions with the cell wall. TolR has been proposed to bind to peptidoglycan, but the only structure of this protein thus far is in a conformation in which the putative peptidoglycan binding domain is not accessible. We show that a model of the extended conformation of the protein in which this domain is exposed binds peptidoglycan largely through electrostatic interactions. Non-covalent interactions of TolR and OmpA with the cell wall, from the inner membrane and outer membrane sides, respectively, maintain the position of the cell wall even in the absence of Braun's lipoprotein. The charged residues that mediate the cell-wall interactions of TolR in our simulations are conserved across a number of species of gram-negative bacteria.

Atomistic and Coarse Grain Simulations of the Cell Envelope of Gram-Negative Bacteria: What Have We Learned?

Bacterial membranes, and those of Gram-negative bacteria in particular, are some of the most biochemically diverse membranes known. They incorporate a wide range of lipid types and proteins of varying sizes, architectures, and functions. While simpler biological membranes have been the focus of myriad simulation studies over the years that have yielded invaluable details to complement, and often to direct, ongoing experimental studies, simulations of complex bacterial membranes have been slower to emerge. However, the past few years have seen tremendous activity in this area, leading to advances such as the development of atomistic and coarse-grain models of the lipopolysaccharide (LPS) component of the outer membrane that are compatible with widely used simulation codes. In this Account, we review our contributions to the field of molecular simulations of the bacterial cell envelope, including the development of models of both membranes and the cell wall of Gram-negative bacteria, with a predominant focus on E. coli. At the atomistic level, simulations of chemically accurate models of both membranes have revealed the tightly cross-linked nature of the LPS headgroups and have shown that penetration of solutes through these regions is not as straightforward as the route through phospholipids. The energetic differences between the two routes have been calculated. Simulations of native outer membrane proteins in LPS-containing membranes have shown that the conformational dynamics of the proteins is not only slower in LPS but also different compared to in simpler models of phospholipid bilayers. These chemically more complex and consequently biologically more relevant models are leading to details of conformational dynamics that were previously inaccessible from simulations. Coarse-grain models have enabled simulations of multiprotein systems on time scales of microseconds, leading to insights not only into the rates of protein and lipid diffusion but also into the trends in their respective directions of flow. We find that the motions of LPS molecules are highly correlated with each other but also with outer membrane proteins embedded within the membrane. We have shown that the two leaflets of the outer membrane exhibit communication, whereby regions of low disorder in one leaflet correspond to regions of high disorder in the other. The cell wall remains a comparatively neglected component, although models of the E. coli peptidoglycan are now emerging, particularly at the atomistic level. Our simulations of Braun's lipoprotein have shown that bending and tilting of this protein afford a degree of variability in the gap between the cell wall and the OM. The noncovalent interactions with the cell wall of proteins such as OmpA can further influence the width of this gap by extension or contraction of their linker domains. Overall we have shown that the dynamics of proteins, lipids, and other molecular species within the outer membrane cannot be approximated using simpler phospholipid bilayers, if one is addressing questions regarding the in vivo behavior of Gram-negative bacteria. These membranes have their own unique chemical characteristics that cannot be decoupled from their biological functions.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

A Practical View of the Martini Force Field

Martini is a coarse-grained (CG) force field suitable for molecular dynamics (MD) simulations of (bio)molecular systems. It is based on mapping of two to four heavy atoms to one CG particle. The effective interactions between the CG particles are parametrized to reproduce partitioning free energies of small chemical compounds between polar and apolar phases. In this chapter, a summary of the key elements of this CG force field is presented, followed by an example of practical application: a lipoplex-membrane fusion experiment. Formulated as hands-on practice, this chapter contains guidelines to build CG models of important biological systems, such as asymmetric bilayers and double-stranded DNA. Finally, a series of notes containing useful information, limitations, and tips are described in the last section.

Protein assemblies ejected directly from native membranes yield complexes for mass spectrometry

Membrane proteins reside in lipid bilayers and are typically extracted from this environment for study, which often compromises their integrity. In this work, we ejected intact assemblies from membranes, without chemical disruption, and used mass spectrometry to define their composition. From Escherichia coli outer membranes, we identified a chaperone-porin association and lipid interactions in the β-barrel assembly machinery. We observed efflux pumps bridging inner and outer membranes, and from inner membranes we identified a pentameric pore of TonB, as well as the protein-conducting channel SecYEG in association with F1FO adenosine triphosphate (ATP) synthase. Intact mitochondrial membranes from Bos taurus yielded respiratory complexes and fatty acid-bound dimers of the ADP (adenosine diphosphate)/ATP translocase (ANT-1). These results highlight the importance of native membrane environments for retaining small-molecule binding, subunit interactions, and associated chaperones of the membrane proteome.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

How nanoscale protein interactions determine the mesoscale dynamic organisation of bacterial outer membrane proteins

The spatiotemporal organisation of membranes is often characterised by the formation of large protein clusters. In Escherichia coli, outer membrane protein (OMP) clustering leads to OMP islands, the formation of which underpins OMP turnover and drives organisation across the cell envelope. Modelling how OMP islands form in order to understand their origin and outer membrane behaviour has been confounded by the inherent difficulties of simulating large numbers of OMPs over meaningful timescales. Here, we overcome these problems by training a mesoscale model incorporating thousands of OMPs on coarse-grained molecular dynamics simulations. We achieve simulations over timescales that allow direct comparison to experimental data of OMP behaviour. We show that specific interaction surfaces between OMPs are key to the formation of OMP clusters, that OMP clusters present a mesh of moving barriers that confine newly inserted proteins within islands, and that mesoscale simulations recapitulate the restricted diffusion characteristics of OMPs.

In Escherichia coli, outer membrane protein (OMP) cluster and form islands, but the origin and behaviour of those clusters remains poorly understood. Here authors use coarse grained molecular dynamics simulation and show that their mesoscale simulations recapitulate the restricted diffusion characteristics of OMPs.

Dynamics of OmpF Trimer Formation in the Bacterial Outer Membrane of Escherichia coli

The self-assembly of outer membrane protein F (OmpF) in the outer membrane of Escherichia coli Gram-negative bacteria was studied using multiscale molecular dynamics simulations. To accommodate the long time scale required for protein assembly, coarse-grained parametrization of E. coli outer membrane lipids was first developed. The OmpF monomers formed stable dimers at specific protein-protein interactions sites irrespective of the lipid membrane environment. The dimer intermediate was asymmetric but provided a template to form a symmetric trimer. Superposition analysis of the self-assembled trimer with the X-ray crystal structure of the trimer available in the protein data bank showed excellent agreement with global root-mean-square deviation of less than 2.2 Å. The free energy change associated with dimer formation was -26 ± 1 kcal mol^-1, and for a dimer to bind to a monomer and to form a trimer yielded -56 ± 4 kcal mol^-1. Based on thermodynamic data, an alternate path to trimer formation via interaction of two dimers is also presented.

In Silico Investigations of Calcium Phosphate Mineralization in Extracellular Vesicles

Calcification in bone, cartilage, and cardiovascular tissues involves the release of specialized extracellular vesicles (EVs) that promote mineral nucleation. The small size of the EVs, however, makes molecular level studies difficult, and consequently uncertainty exists on the role and function of these structures in directing mineralization. The lack of mechanistic understanding associated with the initiators of ectopic mineral deposition has severely hindered the development of potential therapeutic options. Here, we used multiscale molecular dynamics simulations to investigate the calcification within the EVs. Results show that Ca²⁺-HPO₄^2- and phosphatidylserine complexes facilitate the early nucleation. Use of coarse-grained simulations allows investigations of Ca²⁺-PO₄^3- nucleation and crystallization in the EVs. Systematic variation in the ion-to-water ratio shows that the crystallization and growth strongly depend on the enrichment of the ions and dehydration inside the EVs. Our investigations provide insights into the role of EVs on calcium phosphate mineral nucleation and growth in both physiological and pathological mineralization.

Communication between the leaflets of asymmetric membranes revealed from coarse-grain molecular dynamics simulations

We use coarse-grain molecular simulations to investigate the structural and dynamics differences between an asymmetric and a symmetrical membrane, both containing beta barrel transmembrane proteins. We find in where the dynamics of the two leaflets differ greatly, the slowest leaflet dominates the structural effects and importance of protein-lipid interactions.