Modeling and Simulation of Bacterial Outer Membranes with Lipopolysaccharides and Enterobacterial Common Antigen

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

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

Seeing is believing: Illuminating the Gram-negative outer membrane with molecular dynamics simulations

Recent advances in molecular dynamics (MD) simulations have led to rapid improvement in our understanding of the molecular details of the outer membranes (OMs) of Gram-negative bacteria. In this review, we highlight the latest discoveries from MD simulations of OMs, shedding light on the dynamic nature of these bacteria's first line of defense. With the focus on cutting-edge approaches, we explore the OM's sensitivity to structural features, including divalent cations and membrane composition, which have emerged as crucial determinants of antimicrobial passage. Additionally, studies have provided novel insights into outer-membrane proteins (OMPs), revealing their intricate roles in substrate translocation and their distinct interactions with lipopolysaccharides (LPS) in the OM. Finally, we explore the challenging process of β-barrel membrane protein insertion, showcasing recent findings that have enhanced our grasp of this fundamental biological phenomenon.

Molecular Dynamics of Outer Membrane-Embedded Polysaccharide Secretion Porins Reveals Closed Resting-State Surface Gates Targetable by Virtual Fragment Screening for Drug Hotspot Identification

Recent advances in iterative neural network analyses (e.g., AlphaFold2 and RoseTTA fold) have been revolutionary for protein 3D structure prediction, especially for difficult-to-manipulate α-helical/β-barrel integral membrane proteins. These model structures are calculated based on the coevolution of amino acids within the protein of interest and similarities to existing protein structures; the local effects of the membrane on folding and stability of the calculated model structures are not considered. We recently reported the discovery, 3D modeling, and characterization of 18-β-stranded outer-membrane (OM) WzpX, WzpS, and WzpB β-barrel secretion porins for the exopolysaccharide (EPS), major spore coat polysaccharide (MASC), and biosurfactant polysaccharide (BPS) pathways (respectively) in the Gram-negative social predatory bacterium Myxococcus xanthus DZ2. However, information was not obtained regarding the dynamic behavior of surface-gating WzpX/S/B loop domains or on potential treatments to inactivate these porins. Herein, we developed a molecular dynamics (MD) protocol to study the core stability and loop dynamism of neural network-based integral membrane protein structure models embedded in an asymmetric OM bilayer, using the M. xanthus WzpX, WzpS, and WzpB proteins as test candidates. This was accomplished through integration of the CHARMM-graphical user interface (GUI) and Molecular Operating Environment (MOE) workflows to allow for a rapid simulation system setup and facilitate data analysis. In addition to serving as a method of model structure validation, our molecular dynamics simulations revealed a minimal movement of extracellular WzpX/S/B loops in the absence of an external stimulus as well as druggable cavities between the loops. Virtual screening of a commercial fragment library against these cavities revealed putative fragment-binding hotspots on the cell-surface face of each β-barrel, along with key interacting residues, and identified promising hits for the design of potential binders capable of plugging the β-barrels and inhibiting polysaccharide secretion.

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

Martini-3 Coarse-Grained Models for the Bacterial Lipopolysaccharide Outer Membrane of Escherichia coli

The outer lipopolysaccharide (LPS) membrane of Gram-negative bacteria forms the main barrier for transport of antimicrobial molecules into the bacterial cell. In this study we develop coarse-grained models for the outer membrane of Escherichia coli in the Martini-3 framework. The coarse-grained model force field was parametrized and validated using all-atom simulations of symmetric membranes of lipid A and rough LPS as well as a complete asymmetric membrane of LPS with the O-antigen. The bonded parameters were obtained using an iterative refinement procedure with target bonded distributions obtained from all-atom simulations. The membrane thickness, area of the LPS, and density distributions for the different regions as well as the water and ion densities in Martini-3 simulations show excellent agreement with the all-atom data. Additionally the solvent accessible surface area for individual molecules in water was found to be in good agreement. The binding of calcium ions with phosphate and carboxylate moieties of LPS is accurately captured in the Martini-3 model, indicative of the integrity of the highly negatively charged LPS molecules in the outer membranes of Gram-negative bacteria. The melting transition of the coarse-grained lipid A membrane model was found to occur between 300 and 310 K, and the model captured variations in area per LPS, order parameter, and membrane thickness across the melting transition. Our study reveals that the proposed Martini-3 models for LPS are able to capture the physicochemical balance of the complex sugar architecture of the outer membrane of Escherichia coli. The coarse-grained models developed in this study would be useful for determining membrane protein interactions and permeation of potential antimicrobials through bacterial membranes at mesoscopic spatial and temporal scales.

Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy

Glycans in the form of oligosaccharides, polysaccharides, and glycoconjugates are ubiquitous in nature, and their structures range from linear assemblies to highly branched and decorated constructs. Solution state NMR spectroscopy facilitates elucidation of preferred conformations and shapes of the saccharides, motions, and dynamic aspects related to processes over time as well as the study of transient interactions with proteins. Identification of intermolecular networks at the atomic level of detail in recognition events by carbohydrate-binding proteins known as lectins, unraveling interactions with antibodies, and revealing substrate scope and action of glycosyl transferases employed for synthesis of oligo- and polysaccharides may efficiently be analyzed by NMR spectroscopy. By utilizing NMR active nuclei present in glycans and derivatives thereof, including isotopically enriched compounds, highly detailed information can be obtained by the experiments. Subsequent analysis may be aided by quantum chemical calculations of NMR parameters, machine learning-based methodologies and artificial intelligence. Interpretation of the results from NMR experiments can be complemented by extensive molecular dynamics simulations to obtain three-dimensional dynamic models, thereby clarifying molecular recognition processes involving the glycans.

Alternating L4 loop architecture of the bacterial polysaccharide co-polymerase WzzE

Lipopolysaccharides such as the enterobacterial common antigen are important components of the enterobacterial cell envelope that act as a protective barrier against the environment and are often polymerized by the inner membrane bound Wzy-dependent pathway. By employing cryo-electron microscopy we show that WzzE, the co-polymerase component of this pathway that is responsible for the length modulation of the enterobacterial common antigen, is octameric with alternating up-down conformations of its L4 loops. The alternating up-down nature of these essential loops, located at the top of the periplasmic bell, are modulated by clashing helical faces between adjacent protomers that flank the L4 loops around the octameric periplasmic bell. This alternating arrangement and a highly negatively charged binding face create a dynamic environment in which the polysaccharide chain is extended, and suggest a ratchet-type mechanism for polysaccharide elongation.

Gut microbial composition and functionality of school-age Mexican population with metabolic syndrome and type-2 diabetes mellitus using shotgun metagenomic sequencing

Gut metagenome in pediatric subjects with metabolic syndrome (MetS) and type-2 diabetes mellitus (T2DM) has been poorly studied, despite an alarming worldwide increase in the prevalence and incidence of obesity and MetS within this population. The objective of this study was to characterize the gut microbiome taxonomic composition of Mexican pediatric subjects with MetS and T2DM using shotgun metagenomics and analyze the potential relationship with metabolic changes and proinflammatory effects. Paired-end reads of fecal DNA samples were obtained through the Illumina HiSeq X Platform. Statistical analyses and correlational studies were conducted using gut microbiome data and metadata from all individuals. Gut microbial dysbiosis was observed in MetS and T2DM children compared to healthy subjects, which was characterized by an increase in facultative anaerobes (i.e., enteric and lactic acid bacteria) and a decrease in strict anaerobes (i.e., Erysipelatoclostridium, Shaalia, and Actinomyces genera). This may cause a loss of gut hypoxic environment, increased gut microbial nitrogen metabolism, and higher production of pathogen-associated molecular patterns. These metabolic changes may trigger the activation of proinflammatory activity and impair the host's intermediate metabolism, leading to a possible progression of the characteristic risk factors of MetS and T2DM, such as insulin resistance, dyslipidemia, and an increased abdominal circumference. Furthermore, specific viruses (Jiaodavirus genus and Inoviridae family) showed positive correlations with proinflammatory cytokines involved in these metabolic diseases. This study provides novel evidence for the characterization of MetS and T2DM pediatric subjects in which the whole gut microbial composition has been characterized. Additionally, it describes specific gut microorganisms with functional changes that may influence the onset of relevant health risk factors.

Modeling and Simulation of Bacterial Outer Membranes with Lipopolysaccharides and Capsular Polysaccharides

Capsule is one of the common virulence factors in Gram-negative bacteria protecting pathogens from host defenses and consists of long-chain capsular polysaccharides (CPS) anchored in the outer membrane (OM). Elucidating structural properties of CPS is important to understand its biological functions as well as the OM properties. However, the outer leaflet of the OM in current simulation studies is represented exclusively by LPS due to the complexity and diversity of CPS. In this work, representative Escherichia coli CPS, K_LPS (a lipid A-linked form) and K_PG (a phosphatidylglycerol-linked form), are modeled and incorporated into various symmetric bilayers with co-existing LPS in different ratios. All-atom molecular dynamics simulations of these systems have been conducted to characterize various bilayer properties. Incorporation of K_LPS makes the acyl chains of LPS more rigid and ordered, while incorporation of K_PG makes them less ordered and flexible. These results are consistent with the calculated area per lipid (APL) of LPS, in which the APL of LPS becomes smaller when K_LPS is incorporated, whereas it gets larger when K_PG is included. Torsional analysis reveals that the influence of the CPS presence on the conformational distributions of the glycosidic linkages of LPS is small, and minor differences are also detected for the inner and outer regions of the CPS. Combined with previously modeled enterobacterial common antigens (ECAs) in the form of mixed bilayers, this work provides more realistic OM models as well as the basis for characterization of interactions between the OM and OM proteins.

Phosphatidylglycerol Is the Lipid Donor for Synthesis of Phospholipid-Linked Enterobacterial Common Antigen

The Gram-negative outer membrane (OM) is an asymmetric bilayer with phospholipids in its inner leaflet and mainly lipopolysaccharide (LPS) in its outer leaflet and is largely impermeable to many antibiotics. In Enterobacterales (e.g., Escherichia, Salmonella, Klebsiella, and Yersinia), the outer leaflet of the OM also contains phosphoglyceride-linked enterobacterial common antigen (ECA_PG). This molecule consists of the conserved ECA carbohydrate linked to diacylglycerol-phosphate (DAG-P) through a phosphodiester bond. ECA_PG contributes to the OM permeability barrier and modeling suggests that it may alter the packing of LPS molecules in the OM. Here, we investigate, in Escherichia coli K-12, the reaction synthesizing ECA_PG from ECA precursor linked to an isoprenoid carrier to identify the lipid donor that provides the DAG-P moiety to ECA_PG. Through overexpression of phospholipid biosynthesis genes, we observed alterations expected to increase levels of phosphatidylglycerol (PG) increased the synthesis of ECA_PG, whereas alterations expected to decrease levels of PG decreased the synthesis of ECA_PG. We discovered depletion of PG levels in strains that could synthesize ECA_PG, but not other forms of ECA, causes additional growth defects, likely due to the buildup of ECA precursor on the isoprenoid carrier inhibiting peptidoglycan biosynthesis. Our results demonstrate ECA_PG can be synthesized in the absence of the other major phospholipids (phosphatidylethanolamine and cardiolipin). Overall, these results conclusively demonstrate PG is the lipid donor for the synthesis of ECA_PG and provide a key insight into the reaction producing ECA_PG. In addition, these results provide an interesting parallel to lipoprotein acylation, which also uses PG as its DAG donor. IMPORTANCE The Gram-negative outer membrane is a permeability barrier preventing cellular entry of antibiotics. However, outer membrane biogenesis pathways are targets for small molecule development. Here, we investigate the synthesis of a form of enterobacterial common antigen (ECA), ECA_PG, found in the outer membrane of Enterobacterales (e.g., Escherichia, Salmonella, and Klebsiella). ECA_PG consists of the conserved ECA carbohydrate unit linked to diacylglycerol-phosphate-ECA is a phospholipid headgroup. The details of the reaction forming this molecule from polymerized ECA precursor are unknown. We determined the lipid donor providing the phospholipid moiety is phosphatidylglycerol. Understanding the synthesis of outer membrane constituents such as ECA_PG provides the opportunity for development of molecules to increase outer membrane permeability, expanding the antibiotics available to treat Gram-negative infections.

Computer-Aided Drug Design: An Update

Computer-aided drug design (CADD) approaches are playing an increasingly important role in understanding the fundamentals of ligand-receptor interactions and helping medicinal chemists design therapeutics. About 5 years ago, we presented a chapter devoted to an overview of CADD methods and covered typical CADD protocols including structure-based drug design (SBDD) and ligand-based drug design (LBDD) approaches that were frequently used in the antibiotic drug design process. Advances in computational hardware and algorithms and emerging CADD methods are enhancing the accuracy and ability of CADD in drug design and development. In this chapter, an update to our previous chapter is provided with a focus on new CADD approaches from our laboratory and other peers that can be employed to facilitate the development of antibiotic therapeutics.

The Role of O-antigen in P1 Transduction of Shigella flexneri and Escherichia coli with its Alternative S' Tail Fibre

Enterobacteria phage P1 expresses two types of tail fibre, S and S'. Despite the wide usage of phage P1 for transduction, the host range and the receptor for its alternative S' tail fibre was never determined. Here, a ΔS-cin Δpac E. coli P1 lysogenic strain was generated to allow packaging of phagemid DNA into P1 phage having either S or S' tail fibre. P1(S') could transduce phagemid DNA into Shigella flexneri 2a 2457O, Shigella flexneri 5a M90T and Escherichia coli O3 efficiently. Mutational analysis of the O-antigen assembly genes and LPS inhibition assays indicated that P1(S') transduction requires at least one O-antigen unit. E. coli O111:B4 LPS produced a high neutralising effect against P1(S') transduction, indicating that this E. coli strain could be susceptible to P1(S')-mediated transduction. Mutations in the O-antigen modification genes of S. flexneri 2a 2457O and S. flexneri 5a M90T did not cause significant changes to P1(S') transduction efficiency. A higher transduction efficiency of P1(S') improved the delivery of a cas9 antimicrobial phagemid into both S. flexneri 2457O and M90T. These findings provide novel insights into P1 tropism-switching, by identifying the bacterial strains which are susceptible to P1(S')-mediated transduction, as well as demonstrating its potential for delivering a DNA sequence-specific Cas9 antimicrobial into clinically relevant S. flexneri.

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

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

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Physical properties of the bacterial outer membrane

It has long been appreciated that the Gram-negative outer membrane acts as a permeability barrier, but recent studies have uncovered a more expansive and versatile role for the outer membrane in cellular physiology and viability. Owing to recent developments in microfluidics and microscopy, the structural, rheological and mechanical properties of the outer membrane are becoming apparent across multiple scales. In this Review, we discuss experimental and computational studies that have revealed key molecular factors and interactions that give rise to the spatial organization, limited diffusivity and stress-bearing capacity of the outer membrane. These physical properties suggest broad connections between cellular structure and physiology, and we explore future prospects for further elucidation of the implications of outer membrane construction for cellular fitness and survival.

Commensal and Pathogenic Bacterial-Derived Extracellular Vesicles in Host-Bacterial and Interbacterial Dialogues: Two Sides of the Same Coin

Extracellular vesicles (EVs) cause effective changes in various domains of life. These bioactive structures are essential to the bidirectional organ communication. Recently, increasing research attention has been paid to EVs derived from commensal and pathogenic bacteria in their potential role to affect human disease risk for cancers and a variety of metabolic, gastrointestinal, psychiatric, and mental disorders. The present review presents an overview of both the protective and harmful roles of commensal and pathogenic bacteria-derived EVs in host-bacterial and interbacterial interactions. Bacterial EVs could impact upon human health by regulating microbiota–host crosstalk intestinal homeostasis, even in distal organs. The importance of vesicles derived from bacteria has been also evaluated regarding epigenetic modifications and applications. Generally, the evaluation of bacterial EVs is important towards finding efficient strategies for the prevention and treatment of various human diseases and maintaining metabolic homeostasis.

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.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-021-00913-7.

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.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

A Journey from Structure to Function of Bacterial Lipopolysaccharides

Lipopolysaccharide (LPS) is a crucial constituent of the outer membrane of most Gram-negative bacteria, playing a fundamental role in the protection of bacteria from environmental stress factors, in drug resistance, in pathogenesis, and in symbiosis. During the last decades, LPS has been thoroughly dissected, and massive information on this fascinating biomolecule is now available. In this Review, we will give the reader a third millennium update of the current knowledge of LPS with key information on the inherent peculiar carbohydrate chemistry due to often puzzling sugar residues that are uniquely found on it. Then, we will drive the reader through the complex and multifarious immunological outcomes that any given LPS can raise, which is strictly dependent on its chemical structure. Further, we will argue about issues that still remain unresolved and that would represent the immediate future of LPS research. It is critical to address these points to complete our notions on LPS chemistry, functions, and roles, in turn leading to innovative ways to manipulate the processes involving such a still controversial and intriguing biomolecule.

Solving the structural puzzle of bacterial glycome

The analysis of the bacterial glycome (glycomics) is among the complex 'omics' analysis owing to the inherent difficulties in structural and functional characterization of glycans. The complexity and variability of bacterial glycans, spanning from simple carbohydrates to complex glycolipids, glycopeptides and glycoproteins, make their study a challenging research area. The last two decades have witnessed tremendous advances and development of highly sophisticated methods, in combination with optimized protocols and hyphenate techniques for the understanding of structure, conformations, dynamics and organization of microbial glycans. We here present an overview of the novel approaches that have massively improved our understanding of the carbohydrate-based world of bacteria.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.