Refinement of OprH-LPS Interactions by Molecular Simulations

Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

Advances in Development of Novel Therapeutic Strategies against Multi-Drug Resistant Pseudomonas aeruginosa

Pseudomonas aeruginosa (P. aeruginosa) with multi-drug resistance (MDR) is a major cause of serious healthcare-associated infections, leading to high morbidity and mortality. This opportunistic pathogen is responsible for various infectious diseases, such as those seen in cystic fibrosis, ventilator-associated pneumonia, urinary tract infection, otitis externa, and burn and wound injuries. Due to its relatively large genome, P. aeruginosa has great diversity and can use various molecular mechanisms for antimicrobial resistance. For example, outer membrane permeability can contribute to antimicrobial resistance and is determined by lipopolysaccharide (LPS) and porin proteins. Recent findings on the regulatory interaction between peptidoglycan and LPS synthesis provide additional clues against pathogenic P. aeruginosa. This review focuses on recent advances in antimicrobial agents and inhibitors targeting LPS and porin proteins. In addition, we explore current and emerging treatment strategies for MDR P. aeruginosa, including phages, vaccines, nanoparticles, and their combinatorial therapies. Novel strategies and their corresponding therapeutic agents are urgently needed for combating MDR pathogens.

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

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

Measurement of Small Molecule Accumulation into Diderm Bacteria

Some of the most dangerous bacterial pathogens (Gram-negative and mycobacterial) deploy a formidable secondary membrane barrier to reduce the influx of exogenous molecules. For Gram-negative bacteria, this second exterior membrane is known as the outer membrane (OM), while for the Gram-indeterminate Mycobacteria, it is known as the "myco" membrane. Although different in composition, both the OM and mycomembrane are key structures that restrict the passive permeation of small molecules into bacterial cells. Although it is well-appreciated that such structures are principal determinants of small molecule permeation, it has proven to be challenging to assess this feature in a robust and quantitative way or in complex, infection-relevant settings. Herein, we describe the development of the bacterial chloro-alkane penetration assay (BaCAPA), which employs the use of a genetically encoded protein called HaloTag, to measure the uptake and accumulation of molecules into model Gram-negative and mycobacterial species, Escherichia coli and Mycobacterium smegmatis, respectively, and into the human pathogen Mycobacterium tuberculosis. The HaloTag protein can be directed to either the cytoplasm or the periplasm of bacteria. This offers the possibility of compartmental analysis of permeation across individual cell membranes. Significantly, we also showed that BaCAPA can be used to analyze the permeation of molecules into host cell-internalized E. coli and M. tuberculosis, a critical capability for analyzing intracellular pathogens. Together, our results show that BaCAPA affords facile measurement of permeability across four barriers: the host plasma and phagosomal membranes and the diderm bacterial cell envelope.

Conformational dynamics and putative substrate extrusion pathways of the N-glycosylated outer membrane factor CmeC from Campylobacter jejuni

The outer membrane factor CmeC of the efflux machinery CmeABC plays an important role in conferring antibiotic and bile resistance to Campylobacter jejuni. Curiously, the protein is N-glycosylated, with the glycans playing a key role in the effective function of this system. In this work we have employed atomistic equilibrium molecular dynamics simulations of CmeC in a representative model of the C. jejuni outer membrane to characterise the dynamics of the protein and its associated glycans. We show that the glycans are more conformationally labile than had previously been thought. The extracellular loops of CmeC visit the open and closed states freely suggesting the absence of a gating mechanism on this side, while the narrow periplasmic entrance remains tightly closed, regulated via coordination to solvated cations. We identify several cation binding sites on the interior surface of the protein. Additionally, we used steered molecular dynamics simulations to elucidate translocation pathways for a bile acid and a macrolide antibiotic. These, and additional equilibrium simulations suggest that the anionic bile acid utilises multivalent cations to climb the ladder of acidic residues that line the interior surface of the protein.

Pseudomonas aeruginosa: pathogenesis, virulence factors, antibiotic resistance, interaction with host, technology advances and emerging therapeutics

Pseudomonas aeruginosa (P. aeruginosa) is a Gram-negative opportunistic pathogen that infects patients with cystic fibrosis, burn wounds, immunodeficiency, chronic obstructive pulmonary disorder (COPD), cancer, and severe infection requiring ventilation, such as COVID-19. P. aeruginosa is also a widely-used model bacterium for all biological areas. In addition to continued, intense efforts in understanding bacterial pathogenesis of P. aeruginosa including virulence factors (LPS, quorum sensing, two-component systems, 6 type secretion systems, outer membrane vesicles (OMVs), CRISPR-Cas and their regulation), rapid progress has been made in further studying host-pathogen interaction, particularly host immune networks involving autophagy, inflammasome, non-coding RNAs, cGAS, etc. Furthermore, numerous technologic advances, such as bioinformatics, metabolomics, scRNA-seq, nanoparticles, drug screening, and phage therapy, have been used to improve our understanding of P. aeruginosa pathogenesis and host defense. Nevertheless, much remains to be uncovered about interactions between P. aeruginosa and host immune responses, including mechanisms of drug resistance by known or unannotated bacterial virulence factors as well as mammalian cell signaling pathways. The widespread use of antibiotics and the slow development of effective antimicrobials present daunting challenges and necessitate new theoretical and practical platforms to screen and develop mechanism-tested novel drugs to treat intractable infections, especially those caused by multi-drug resistance strains. Benefited from has advancing in research tools and technology, dissecting this pathogen's feature has entered into molecular and mechanistic details as well as dynamic and holistic views. Herein, we comprehensively review the progress and discuss the current status of P. aeruginosa biophysical traits, behaviors, virulence factors, invasive regulators, and host defense patterns against its infection, which point out new directions for future investigation and add to the design of novel and/or alternative therapeutics to combat this clinically significant pathogen.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

Model architectures for bacterial membranes

The complex composition of bacterial membranes has a significant impact on the understanding of pathogen function and their development towards antibiotic resistance. In addition to the inherent complexity and biosafety risks of studying biological pathogen membranes, the continual rise of antibiotic resistance and its significant economical and clinical consequences has motivated the development of numerous in vitro model membrane systems with tuneable compositions, geometries, and sizes. Approaches discussed in this review include liposomes, solid-supported bilayers, and computational simulations which have been used to explore various processes including drug-membrane interactions, lipid-protein interactions, host–pathogen interactions, and structure-induced bacterial pathogenesis. The advantages, limitations, and applicable analytical tools of all architectures are summarised with a perspective for future research efforts in architectural improvement and elucidation of resistance development strategies and membrane-targeting antibiotic mechanisms.

Lipopolysaccharides at Solid and Liquid Interfaces: Models for Biophysical Studies of the Gram-negative Bacterial Outer Membrane

Lipopolysaccharides (LPSs) are a constitutive element of the cell envelope of Gram-negative bacteria, representing the main lipid in the external leaflet of their outer membrane (OM) lipid bilayer. These unique surface-exposed glycolipids play a central role in the interactions of Gram-negative organisms with their surrounding environment and represent a key element for protection against antimicrobials and the development of antibiotic resistance. The biophysical investigation of a wide range of different types of in vitro model membranes containing reconstituted LPS has revealed functional and structural properties of these peculiar membrane lipids, providing molecular-level details of their interaction with antimicrobial compounds. LPS assemblies reconstituted at interfaces represent a versatile tool to study the properties of the Gram-negative OM by exploiting several surface-sensitive techniques, in particular X-ray and neutron scattering, which can probe the structure of thin films with sub-nanometer resolution. This review provides an overview of different approaches employed to investigate structural and biophysical properties of LPS, focusing on studies on Langmuir monolayers of LPS at the air/liquid interface and a range of supported LPS-containing model membranes reconstituted at solid/liquid interfaces.

Membrane Potentials Trigger Molecular-Scale Rearrangements in the Outer Membrane of Gram-Negative Bacteria

The structural complexity of the cell envelope of Gram-negative bacteria limits the fabrication of realistic models of bacterial cell membranes. A vertical Langmuir-Blodgett withdrawing was used to deposit a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) monolayer on the Au(111) surface. The second leaflet composed of di[3-deoxy-D-manno-octulosonyl]-lipid A (KLA) was deposited using Langmuir-Schaefer transfer. The use of an electrode material as a support for the POPE-KLA bilayer allowed electrochemical control of the membrane's stability, compactness, and structure. Capacitance-potential curves showed a typical pattern for the supported lipid bilayers electrochemical characteristic. The minimum membrane capacitance was ∼4 μF cm^-2 and did not change in the following desorption-adsorption cycles, indicating the presence of a stable bilayer structure with an asymmetric composition of both leaflets. However, at a molecular scale, as elucidated in spectroelectrochemical experiments, large differences in the response of both leaflets to electric potentials were observed. The acyl chains in POPE and KLA existed in a liquid state. The quantitative analysis of the CH stretching modes indicated potential-driven reorientations in the hydrophobic fragment of the bilayer, already in the adsorbed state. To assign observed rearrangements to POPE and KLA lipids in both leaflets, per-deuterated d₃₁-POPE was transferred into the inner leaflet. Since no potential-dependent changes of the CD₂ stretching modes in the d₃₁-POPE-KLA bilayer were observed, reorientations in the acyl chain region were assigned to the KLA molecules. Mg²⁺ ions were bound to the polar head groups of KLA. The strength of electrostatic interactions in the polar head group region of KLA was dependent on the direction of the electric field. At negative electric potentials, the binding of divalent cations weakened, which gave the KLA molecules increased orientational flexibility. This behavior in electric fields is peculiar for the outer membrane and indicates that the microbial cell membranes have different electrochemical properties than phospholipid bilayers.

In silico investigation and surmounting of Lipopolysaccharide barrier in Gram-Negative Bacteria: How far has molecular dynamics Come?

Lipopolysaccharide (LPS), a main component of the outer membrane of Gram-negative bacteria, has crucial implications on both antibiotic resistance and the overstimulation of the host innate immune system. Fighting against these global concerns calls for the molecular understanding of the barrier function and immunostimulatory ability of LPS. Molecular dynamics (MD) simulations have become an invaluable tool for uncovering important findings in LPS research. While the reach of MD simulations for investigating the immunostimulatory ability of LPS has been already outlined, little attention has been paid to the role of MD simulations for exploring its barrier function and synthesis. Herein, we give an overview about the impact of MD simulations on gaining insight into the shield role and synthesis pathway of LPS, which have attracted considerable attention to discover molecules able to surmount antibiotic resistance, either circumventing LPS defenses or disrupting its synthesis. We specifically focus on the enhanced sampling and free energy calculation methods that have been combined with MD simulations to address such research. We also highlight the use of special-purpose MD supercomputers, the importance of appropriate LPS and ions parameterization to obtain reliable results, and the complementary views that MD and wet-lab experiments provide. Thereby, this work, which covers the last five years of research, apart from outlining the phenomena and strategies that are being explored, evidences the valuable insights that are gained by MD, which may be useful to advance antibiotic design, and what the prospects of this in silico method could be in LPS research.

Fighting Against Bacterial Lipopolysaccharide-Caused Infections through Molecular Dynamics Simulations: A Review

Lipopolysaccharide (LPS) is the primary component of the outer leaflet of Gram-negative bacterial outer membranes. LPS elicits an overwhelming immune response during infection, which can lead to life-threatening sepsis or septic shock for which no suitable treatment is available so far. As a result of the worldwide expanding multidrug-resistant bacteria, the occurrence and frequency of sepsis are expected to increase; thus, there is an urge to develop novel strategies for treating bacterial infections. In this regard, gaining an in-depth understanding about the ability of LPS to both stimulate the host immune system and interact with several molecules is crucial for fighting against LPS-caused infections and allowing for the rational design of novel antisepsis drugs, vaccines and LPS sequestration and detection methods. Molecular dynamics (MD) simulations, which are understood as being a computational microscope, have proven to be of significant value to understand LPS-related phenomena, driving and optimizing experimental research studies. In this work, a comprehensive review on the methods that can be combined with MD simulations, recently applied in LPS research, is provided. We focus especially on both enhanced sampling methods, which enable the exploration of more complex systems and access to larger time scales, and free energy calculation approaches. Thereby, apart from outlining several strategies for surmounting LPS-caused infections, this work reports the current state-of-the-art of the methods applied with MD simulations for moving a step forward in the development of such strategies.

Structure, Assembly, and Function of Tripartite Efflux and Type 1 Secretion Systems in Gram-Negative Bacteria

Tripartite efflux pumps and the related type 1 secretion systems (T1SSs) in Gram-negative organisms are diverse in function, energization, and structural organization. They form continuous conduits spanning both the inner and the outer membrane and are composed of three principal components-the energized inner membrane transporters (belonging to ABC, RND, and MFS families), the outer membrane factor channel-like proteins, and linking the two, the periplasmic adaptor proteins (PAPs), also known as the membrane fusion proteins (MFPs). In this review we summarize the recent advances in understanding of structural biology, function, and regulation of these systems, highlighting the previously undescribed role of PAPs in providing a common architectural scaffold across diverse families of transporters. Despite being built from a limited number of basic structural domains, these complexes present a staggering variety of architectures. While key insights have been derived from the RND transporter systems, a closer inspection of the operation and structural organization of different tripartite systems reveals unexpected analogies between them, including those formed around MFS- and ATP-driven transporters, suggesting that they operate around basic common principles. Based on that we are proposing a new integrated model of PAP-mediated communication within the conformational cycling of tripartite systems, which could be expanded to other types of assemblies.

Pseudomonas aeruginosa: An Audacious Pathogen with an Adaptable Arsenal of Virulence Factors

Pseudomonas aeruginosa is a dominant pathogen in people with cystic fibrosis (CF) contributing to morbidity and mortality. Its tremendous ability to adapt greatly facilitates its capacity to cause chronic infections. The adaptability and flexibility of the pathogen are afforded by the extensive number of virulence factors it has at its disposal, providing P. aeruginosa with the facility to tailor its response against the different stressors in the environment. A deep understanding of these virulence mechanisms is crucial for the design of therapeutic strategies and vaccines against this multi-resistant pathogen. Therefore, this review describes the main virulence factors of P. aeruginosa and the adaptations it undergoes to persist in hostile environments such as the CF respiratory tract. The very large P. aeruginosa genome (5 to 7 MB) contributes considerably to its adaptive capacity; consequently, genomic studies have provided significant insights into elucidating P. aeruginosa evolution and its interactions with the host throughout the course of infection.

Structural insights into the membrane receptor ShuA in DDM micelles and in a model of gram-negative bacteria outer membrane as seen by SAXS and MD simulations

Successful crystallization of membrane proteins in detergent micelles depends on key factors such as conformational stability of the protein in micellar assemblies, the protein-detergent complex (PDC) monodispersity and favorable protein crystal contacts by suitable shielding of the protein hydrophobic surface by the detergent belt. With the aim of studying the influence of amphiphilic environment on membrane protein structure, stability and crystallizability, we combine molecular dynamics (MD) simulations with SEC-MALLS and SEC-SAXS (Size Exclusion Chromatography in line with Multi Angle Laser Light Scattering or Small Angle X-ray Scattering) experiments to describe the protein-detergent interactions that could help to rationalize PDC crystallization. In this context, we compare the protein-detergent interactions of ShuA from Shigella dysenteriae in n-Dodecyl-β-D-Maltopyranoside (DDM) with ShuA inserted in a realistic model of gram-negative bacteria outer membrane (OM) containing a mixture of bacterial lipopolysaccharide and phospholipids. To evaluate the quality of the PDC models, we compute the corresponding SAXS curves from the MD trajectories and compare with the experimental ones. We show that computed SAXS curves obtained from the MD trajectories reproduce better the SAXS obtained from the SEC-SAXS experiments for ShuA surrounded by 268 DDM molecules. The MD results show that the DDM molecules form around ShuA a closed belt whose the hydrophobic thickness appears slightly smaller (~22 Å) than the hydrophobic transmembrane domain of the protein (24.6 Å) suggested by Orientations of Proteins in Membranes (OPM) database. The simulations also show that ShuA transmembrane domain is remarkably stable in all the systems except for the extracellular and periplasmic loops that exhibit larger movements due to specific molecular interactions with lipopolysaccharides (LPS). We finally point out that this detergent behavior may lead to the occlusion of the periplasmic hydrophilic surface and poor crystal contacts leading to difficulties in crystallization of ShuA in DDM.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Mutually constructive roles of Ail and LPS in Yersinia pestis serum survival

The outer membrane is a key virulence determinant of gram-negative bacteria. In Yersinia pestis, the deadly agent that causes plague, the protein Ail and lipopolysaccharide (LPS)⁶ enhance lethality by promoting resistance to human innate immunity and antibiotics, enabling bacteria to proliferate in the human host. Their functions are highly coordinated. Here we describe how they cooperate to promote pathogenesis. Using a multidisciplinary approach, we identify mutually constructive interactions between Ail and LPS that produce an extended conformation of Ail at the membrane surface, cause thickening and rigidification of the LPS membrane, and collectively promote Y. pestis survival in human serum, antibiotic resistance, and cell envelope integrity. The results highlight the importance of the Ail-LPS assembly as an organized whole, rather than its individual components, and provide a handle for targeting Y. pestis pathogenesis.

Broadening Activity of Polymyxin by Quaternary Ammonium Grafting

Bacterial pathogens continue to impose a tremendous health burden across the globe. Here, we describe a novel series of polymyxin-based agents grafted with membrane-active quaternary ammonium warheads to combine two important classes of Gram-negative antimicrobial scaffolds. The goal was to deliver a targeted quaternary ammonium warhead onto the surface of bacterial pathogens using the outer membrane homing properties of polymyxin. The most potent agents resulted in new scaffolds that retained the ability to target Gram-negative bacteria and had limited toxicity toward mammalian cells. We showed, using a molecular dynamics approach, that the new agents retained their ability to engage in specific interactions with lipopolysaccharide molecules. Significantly, the combination of quaternary ammonium and polymyxin widens the activity to the pathogen Staphylococcus aureus. Our results serve as an example of how two membrane-active agents can be combined to produce a class of novel scaffolds with potent biological activity.

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.

Antibody Binding to the O-Specific Antigen of Pseudomonas aeruginosa O6 Inhibits Cell Growth

Pseudomonas aeruginosa is an opportunistic pathogen that is inherently resistant to many antibiotics and represents an increasing threat due to the emergence of drug-resistant strains. There is a pressing need to develop innovative antimicrobials against this pathogen. In this study, we identified the O-specific antigen (OSA) of P. aeruginosa serotype O6 as a novel target for therapeutic intervention. Binding of monoclonal antibodies and antigen-binding fragments therefrom to O6 OSA leads to rapid outer membrane destabilization and inhibition of cell growth. The antimicrobial effect correlated directly with antibody affinity. Antibody binding to the O antigen of a second lipopolysaccharide (LPS) type present in P. aeruginosa or to the LPS core did not affect cell viability. Atomic force microscopy showed that antibody binding to OSA resulted in early flagellum loss, formation of membrane blebs, and eventually complete outer membrane loss. We hypothesize that antibody binding to OSA disrupts a key interaction in the P. aeruginosa outer membrane.

The current structural glycome landscape and emerging technologies

Carbohydrates represent one of the building blocks of life, along with nucleic acids, proteins and lipids. Although glycans are involved in a wide range of processes from embryogenesis to protein trafficking and pathogen infection, we are still a long way from deciphering the glycocode. In this review, we aim to present a few of the challenges that researchers working in the area of glycobiology can encounter and what strategies can be utilised to overcome them. Our goal is to paint a comprehensive picture of the current saccharide landscape available in the Protein Data Bank (PDB). We also review recently updated repositories relevant to the topic proposed, the impact of software development on strategies to structurally solve carbohydrate moieties, and state-of-the-art molecular and cellular biology methods that can shed some light on the function and structure of glycans.

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.

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.

Physical Properties of Bacterial Outer Membrane Models: Neutron Reflectometry & Molecular Simulation

The outer membrane (OM) of Gram-negative bacteria is an asymmetric bilayer having phospholipids in the inner leaflet and lipopolysaccharides in the outer leaflet. This unique asymmetry and the complex carbohydrates in lipopolysaccharides make it a daunting task to study the asymmetrical OM structure and dynamics, its interactions with OM proteins, and its roles in translocation of substrates, including antibiotics. In this study, we combine neutron reflectometry and molecular simulation to explore the physical properties of OM mimetics. There is excellent agreement between experiment and simulation, allowing experimental testing of the conclusions from simulations studies and also atomistic interpretation of the behavior of experimental model systems, such as the degree of lipid asymmetry, the lipid component (tail, head, and sugar) profiles along the bilayer normal, and lateral packing (i.e., average surface area per lipid). Therefore, the combination of both approaches provides a powerful new means to explore the biological and biophysical behavior of the bacterial OM.

Fosfomycin Permeation through the Outer Membrane Porin OmpF

Fosfomycin is a frequently prescribed drug in the treatment of acute urinary tract infections. It enters the bacterial cytoplasm and inhibits the biosynthesis of peptidoglycans by targeting the MurA enzyme. Despite extensive pharmacological studies and clinical use, the permeability of fosfomycin across the bacterial outer membrane is largely unexplored. Here, we investigate the fosfomycin permeability across the outer membrane of Gram-negative bacteria by electrophysiology experiments as well as by all-atom molecular dynamics simulations including free-energy and applied-field techniques. Notably, in an electrophysiological zero-current assay as well as in the molecular simulations, we found that fosfomycin can rapidly permeate the abundant Escherichia coli porin OmpF. Furthermore, two triple mutants in the constriction region of the porin have been investigated. The permeation rates through these mutants are slightly lower than that of the wild type but fosfomycin can still permeate. Altogether, this work unravels molecular details of fosfomycin permeation through the outer membrane porin OmpF of E. coli and moreover provides hints for understanding the translocation of phosphonic acid antibiotics through other outer membrane pores.

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.

CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and Lipoglycans

Glycolipids (such as glycoglycerolipids, glycosphingolipids, and glycosylphosphatidylinositol) and lipoglycans (such as lipopolysaccharides (LPS), lipooligosaccharides (LOS), mycobacterial lipoarabinomannan, and mycoplasma lipoglycans) are typically found on the surface of cell membranes and play crucial roles in various cellular functions. Characterizing their structure and dynamics at the molecular level is essential to understand their biological roles, but systematic generation of glycolipid and lipoglycan structures is challenging because of great variations in lipid structures and glycan sequences (i.e., carbohydrate types and their linkages). To facilitate the generation of all-atom glycolipid/LPS/LOS structures, we have developed Glycolipid Modeler and LPS Modeler in CHARMM-GUI ( http://www.charmm-gui.org ), a web-based interface that simplifies building of complex biological simulation systems. In addition, we have incorporated these modules into Membrane Builder so that users can readily build a complex symmetric or asymmetric biological membrane system with various glycolipids and LPS/LOS. These tools are expected to be useful in innovative and novel glycolipid/LPS/LOS modeling and simulation research by easing tedious and intricate steps in modeling complex biological systems and shall provide insight into structures, dynamics, and underlying mechanisms of complex glycolipid-/LPS-/LOS-containing biological membrane systems.

Synthetic Immunotherapeutics against Gram-negative Pathogens

While traditional drug discovery continues to be an important platform for the search of new antibiotics, alternative approaches should also be pursued to complement these efforts. We herein designed a class of molecules that decorate bacterial cell surfaces with the goal of re-engaging components of the immune system toward Escherichia coli and Pseudomonas aeruginosa. More specifically, conjugates were assembled using polymyxin B (an antibiotic that inherently attaches to the surface of Gram-negative pathogens) and antigenic epitopes that recruit antibodies found in human serum. We established that the spacer length played a significant role in hapten display within the bacterial cell surface, a result that was confirmed both experimentally and via molecular dynamics simulations. Most importantly, we demonstrated the specific killing of bacteria by our agent in the presence of human serum. By enlisting the immune system, these agents have the potential to pave the way for a potent antimicrobial modality.

Structural Insights into the Yersinia pestis Outer Membrane Protein Ail in Lipid Bilayers

Yersinia pestis the causative agent of plague, is highly pathogenic and poses very high risk to public health. The outer membrane protein Ail (Adhesion invasion locus) is one of the most highly expressed proteins on the cell surface of Y. pestis, and a major target for the development of medical countermeasures. Ail is essential for microbial virulence and is critical for promoting the survival of Y. pestis in serum. Structures of Ail have been determined by X-ray diffraction and solution NMR spectroscopy, but the protein's activity is influenced by the detergents in these samples, underscoring the importance of the surrounding environment for structure-activity studies. Here we describe the backbone structure of Ail, determined in lipid bilayer nanodiscs, using solution NMR spectroscopy. We also present solid-state NMR data obtained for Ail in membranes containing lipopolysaccharide (LPS), a major component of the bacterial outer membranes. The protein in lipid bilayers, adopts the same eight-stranded β-barrel fold observed in the crystalline and micellar states. The membrane composition, however, appears to have a marked effect on protein dynamics, with LPS enhancing conformational order and slowing down the ¹⁵N transverse relaxation rate. The results provide information about the way in which an outer membrane protein inserts and functions in the bacterial membrane.

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

A complex cell envelope, composed of a mixture of lipid types including lipopolysaccharides, protects bacteria from the external environment. Clearly, the proteins embedded within the various components of the cell envelope have an intricate relationship with their local environment. Therefore, to obtain meaningful results, molecular simulations need to mimic as far as possible this chemically heterogeneous system. However, setting up such systems for computational studies is far from trivial, and consequently the vast majority of simulations of outer membrane proteins still rely on oversimplified phospholipid membrane models. This work presents an update of CHARMM-GUI Martini Maker for coarse-grained modeling and simulation of complex bacterial membranes with lipopolysaccharides. The qualities of the outer membrane systems generated by Martini Maker are validated by simulating them in bilayer, vesicle, nanodisc, and micelle environments (with and without outer membrane proteins) using the Martini force field. We expect this new feature in Martini Maker to be a useful tool for modeling large, complicated bacterial outer membrane systems in a user-friendly manner. © 2017 Wiley Periodicals, Inc.

High resolution solid-state NMR spectroscopy of the Yersinia pestis outer membrane protein Ail in lipid membranes

The outer membrane protein Ail (Adhesion invasion locus) is one of the most abundant proteins on the cell surface of Yersinia pestis during human infection. Its functions are expressed through interactions with a variety of human host proteins, and are essential for microbial virulence. Structures of Ail have been determined by X-ray diffraction and solution NMR spectroscopy, but those samples contained detergents that interfere with functionality, thus, precluding analysis of the structural basis for Ail's biological activity. Here, we demonstrate that high-resolution solid-state NMR spectra can be obtained from samples of Ail in detergent-free phospholipid liposomes, prepared with a lipid to protein molar ratio of 100. The spectra, obtained with 13C or 1H detection, have very narrow line widths (0.40-0.60 ppm for 13C, 0.11-0.15 ppm for 1H, and 0.46-0.64 ppm for 15N) that are consistent with a high level of sample homogeneity. The spectra enable resonance assignments to be obtained for N, CO, CA and CB atomic sites from 75 out of 156 residues in the sequence of Ail, including 80% of the transmembrane region. The 1H-detected solid-state NMR 1H/15N correlation spectra obtained for Ail in liposomes compare very favorably with the solution NMR 1H/15N TROSY spectra obtained for Ail in nanodiscs prepared with a similar lipid to protein molar ratio. These results set the stage for studies of the molecular basis of the functional interactions of Ail with its protein partners from human host cells, as well as the development of drugs targeting Ail.