Search

Glycosaminoglycans (GAGs) such as heparin and heparan sulfate (HS) play crucial roles in inflammation and wound healing, serving as regulators of growth factors and pro-inflammatory mediators. In this study, we investigated the influence of heparin/HS on thrombin proteolysis and its interaction with the generated 11 kDa thrombin-derived C-terminal peptides (TCPs). Employing various biochemical and biophysical methods, we demonstrated that 11 kDa TCPs aggregate in the presence of GAGs, including heparin, heparan sulfate, and chondroitin sulfate-B. Circular dichroism analysis demonstrated that 11 kDa TCPs, in the presence of GAGs, adopt a β-sheet structure, a finding supported by thioflavin T1 (ThT) fluorescence measurements and visualization of 11 kDa TCP-heparin complexes using transmission electron microscopy (TEM). Furthermore, our investigations revealed a stronger binding affinity between 11 kDa TCPs and GAGs with higher sulfate group contents. Congruently, in silico simulations showed that interactions between 11 kDa TCPs and heparin/HS are predominantly electrostatic in nature. Collectively, our study suggests that 11 kDa TCPs have the capacity to aggregate in the presence of GAGs, shedding light on their potential roles in inflammation and wound healing., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Artur Schmidtchen reports a relationship with In2cure AB that includes: board membership. Artur Schmidtchen reports a relationship with Xinnate AB that includes: board membership. Ganna Petruk reports a relationship with Xinnate AB that includes: employment. A.S. is a founder of In2cure AB, a parent company of Xinnate AB, companies that are developing therapies based on thrombin-derived peptides and variants. G.P. is employed part-time by Xinnate AB. The other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors.)

The outermost surface layer of any virus is formed by either a capsid shell or envelope. Such layers have traditionally been thought of as immovable structures, but it is becoming apparent that they cannot be viewed exclusively as static architectures protecting the viral genome. A limited number of proteins on the virion surface must perform a multitude of functions in order to orchestrate the viral life cycle, and allostery can regulate their structures at multiple levels of organization, spanning individual molecules, protomers, large oligomeric assemblies, or entire viral surfaces. Here, we review recent contributions from the molecular simulation field to viral surface allostery, with a particular focus on the trimeric spike glycoprotein emerging from the coronavirus surface, and the icosahedral flaviviral envelope complex. As emerging viral pathogens continue to pose a global threat, an improved understanding of viral dynamics and allosteric regulation will prove crucial in developing novel therapeutic strategies., Competing Interests: Declaration of competing interest The authors declare no competing interests., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Importance: Type IV pili and type II secretion systems are members of the widespread type IV filament (T4F) superfamily of nanomachines that assemble dynamic and versatile surface fibers in archaea and bacteria. The assembly and retraction of T4 filaments with diverse surface properties and functions require the plasma membrane platform proteins of the GspF/PilC superfamily. Generally considered dimeric, platform proteins are thought to function as passive transmitters of the mechanical energy generated by the ATPase motor, to somehow promote insertion of pilin subunits into the nascent pilus fibers. Here, we generate and experimentally validate structural predictions that support the trimeric state of a platform protein PulF from a type II secretion system. The PulF trimers form selective proton or sodium channels which might energize pilus assembly using the membrane potential. The conservation of the channel sequence and structural features implies a common mechanism for all T4F assembly systems. We propose a model of the oligomeric PulF-PulE ATPase complex that provides an essential framework to investigate and understand the pilus assembly mechanism., Competing Interests: The authors declare no conflict of interest.

Surgical site infections (SSI) are a clinical and economic burden. Suture-associated SSI may develop when bacteria colonize the suture surface and form biofilms that are resistant to antibiotics. Thrombin-derived C-terminal peptide (TCP)-25 is a host defense peptide with a unique dual mode of action that can target both bacteria and the excessive inflammation induced by bacterial products. The peptide demonstrates therapeutic potential in preclinical in vivo wound infection models. In this study, the authors set out to explore whether TCP-25 can provide a new bioactive innate immune feature to hydrophilic polyglactin sutures (Vicryl). Using a combination of biochemical, biophysical, antibacterial, biofilm, and anti-inflammatory assays in vitro, in silico molecular modeling studies, along with experimental infection and inflammation models in mice, a proof-of-concept that TCP-25 can provide Vicryl sutures with a previously undisclosed host defense capacity, that enables targeting of bacteria, biofilms, and the accompanying inflammatory response, is shown., (© 2023 The Authors. Advanced Healthcare Materials published by Wiley-VCH GmbH.)

The changing landscape of SARS-CoV-2 Spike protein is linked to the emergence of variants, immune-escape and reduced efficacy of the existing repertoire of anti-viral antibodies. The functional activity of neutralizing antibodies is linked to their quaternary changes occurring as a result of antibody-Spike trimer interactions. Here, we reveal the conformational dynamics and allosteric perturbations linked to binding of novel human antibodies and the viral Spike protein. We identified epitope hotspots, and associated changes in Spike dynamics that distinguish weak, moderate and strong neutralizing antibodies. We show the impact of mutations in Wuhan-Hu-1, Delta, and Omicron variants on differences in the antibody-induced conformational changes in Spike and illustrate how these render certain antibodies ineffective. Antibodies with similar binding affinities may induce destabilizing or stabilizing allosteric effects on Spike, with implications for neutralization efficacy. Our results provide mechanistic insights into the functional modes and synergistic behavior of human antibodies against COVID-19 and may assist in designing effective antiviral strategies., (© 2023. The Author(s).)

There is a clinical need for conceptually new treatments that target the excessive activation of inflammatory pathways during systemic infection. Thrombin-derived C-terminal peptides (TCPs) are endogenous anti-infective immunomodulators interfering with CD14-mediated TLR-dependent immune responses. Here we describe the development of a peptide-based compound for systemic use, sHVF18, expressing the evolutionarily conserved innate structural fold of natural TCPs. Using a combination of structure- and in silico-based design, nuclear magnetic resonance spectroscopy, biophysics, mass spectrometry, cellular, and in vivo studies, we here elucidate the structure, CD14 interactions, protease stability, transcriptome profiling, and therapeutic efficacy of sHVF18. The designed peptide displays a conformationally stabilized, protease resistant active innate fold and targets the LPS-binding groove of CD14. In vivo, it shows therapeutic efficacy in experimental models of endotoxin shock in mice and pigs and increases survival in mouse models of systemic polymicrobial infection. The results provide a drug class based on Nature´s own anti-infective principles., (© 2023. Springer Nature Limited.)

Bats are the reservoir for numerous human pathogens, including coronaviruses. Despite many coronaviruses having descended from bat ancestors, little is known about virus-host interactions and broader evolutionary history involving bats. Studies have largely focused on the zoonotic potential of coronaviruses with few infection experiments conducted in bat cells. To determine genetic changes derived from replication in bat cells and possibly identify potential novel evolutionary pathways for zoonotic virus emergence, we serially passaged six human 229E isolates in a newly established Rhinolophus lepidus (horseshoe bat) kidney cell line. Here, we observed extensive deletions within the spike and open reading frame 4 (ORF4) genes of five 229E viruses after passaging in bat cells. As a result, spike protein expression and infectivity of human cells was lost in 5 of 6 viruses, but the capability to infect bat cells was maintained. Only viruses that expressed the spike protein could be neutralized by 229E spike-specific antibodies in human cells, whereas there was no neutralizing effect on viruses that did not express the spike protein inoculated on bat cells. However, one isolate acquired an early stop codon, abrogating spike expression but maintaining infection in bat cells. After passaging this isolate in human cells, spike expression was restored due to acquisition of nucleotide insertions among virus subpopulations. Spike-independent infection of human coronavirus 229E may provide an alternative mechanism for viral maintenance in bats that does not rely on the compatibility of viral surface proteins and known cellular entry receptors. IMPORTANCE Many viruses, including coronaviruses, originated from bats. Yet, we know little about how these viruses switch between hosts and enter human populations. Coronaviruses have succeeded in establishing in humans at least five times, including endemic coronaviruses and the recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In an approach to identify requirements for host switches, we established a bat cell line and adapted human coronavirus 229E viruses by serial passage. The resulting viruses lost their spike protein but maintained the ability to infect bat cells, but not human cells. Maintenance of 229E viruses in bat cells appears to be independent of a canonical spike receptor match, which in turn might facilitate cross-species transmission in bats., Competing Interests: The authors declare no conflict of interest.

Accumulating evidence indicates a potential role for bacterial lipopolysaccharide (LPS) in the overactivation of the immune response during SARS-CoV-2 infection. LPS is recognized by Toll-like receptor 4, mediating proinflammatory effects. We previously reported that LPS directly interacts with SARS-CoV-2 spike (S) protein and enhances proinflammatory activities. Using native gel electrophoresis and hydrogen-deuterium exchange mass spectrometry, we showed that LPS binds to multiple hydrophobic pockets spanning both the S1 and S2 subunits of the S protein. Molecular simulations validated by a microscale thermophoresis binding assay revealed that LPS binds to the S2 pocket with a lower affinity compared to S1, suggesting a role as an intermediate in LPS transfer. Congruently, nuclear factor-kappa B (NF-κB) activation in monocytic THP-1 cells is strongly boosted by S2. Using NF-κB reporter mice followed by bioimaging, a boosting effect was observed for both S1 and S2, with the former potentially facilitated by proteolysis. The Omicron S variant binds to LPS, but with reduced affinity and LPS boosting in vitro and in vivo. Taken together, the data provide a molecular mechanism by which S protein augments LPS-mediated hyperinflammation., (© The Author(s) (2022). Published by Oxford University Press on behalf of Journal of Molecular Cell Biology, CEMCS, CAS.)

SARS-CoV-2 spike (S) protein is crucial for virus invasion in COVID-19. Here, we showed that lipopolysaccharide (LPS) can trigger S protein aggregation at high doses of LPS and S protein. We demonstrated the formation of S protein aggregates by microscopy analyses, aggregation and gel shift assays. LPS at high levels boosts the formation of S protein aggregates as detected by amytracker and thioflavin T dyes that specifically bind to aggregating proteins. We validated the role of LPS by blocking the formation of aggregates by the endotoxin-scavenging thrombin-derived peptide TCP-25. Aggregation-prone sequences in S protein are predicted to be nearby LPS binding sites, while molecular simulations showed stable formation of S protein-LPS higher-order oligomers. Collectively, our results provide evidence of LPS-induced S protein aggregation., (© 2022 The Authors. FEBS Letters published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

The COVID-19 pandemic has prompted a rapid response in vaccine and drug development. Herein, we modeled a complete membrane-embedded SARS-CoV-2 spike glycoprotein and used molecular dynamics simulations with benzene probes designed to enhance discovery of cryptic pockets. This approach recapitulated lipid and host metabolite binding sites previously characterized by cryo-electron microscopy, revealing likely ligand entry routes, and uncovered a novel cryptic pocket with promising druggable properties located underneath the 617-628 loop. A full representation of glycan moieties was essential to accurately describe pocket dynamics. A multi-conformational behavior of the 617-628 loop in simulations was validated using hydrogen-deuterium exchange mass spectrometry experiments, supportive of opening and closing dynamics. The pocket is the site of multiple mutations associated with increased transmissibility found in SARS-CoV-2 variants of concern including Omicron. Collectively, this work highlights the utility of the benzene mapping approach in uncovering potential druggable sites on the surface of SARS-CoV-2 targets., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Recent developments in the SARS-CoV-2 pandemic point to its inevitable transformation into an endemic disease, urging both refinement of diagnostics for emerging variants of concern (VOCs) and design of variant-specific drugs in addition to vaccine adjustments. Exploring the structure and dynamics of the SARS-CoV-2 Spike protein, we argue that the high-mutability characteristic of RNA viruses coupled with the remarkable flexibility and dynamics of viral proteins result in a substantial involvement of allosteric mechanisms. While allosteric effects of mutations should be considered in predictions and diagnostics of new VOCs, allosteric drugs advantageously avoid escape mutations via non-competitive inhibition originating from alternative distal locations. The exhaustive allosteric signaling and probing maps presented herein provide a comprehensive picture of allostery in the spike protein, making it possible to locate potential mutations that could work as new VOC "drivers" and to determine binding patches that may be targeted by newly developed allosteric drugs., Competing Interests: Declaration of interests The authors declare no competing financial interests., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Extensive glycosylation of viral glycoproteins is a key feature of the antigenic surface of viruses and yet glycan processing can also be influenced by the manner of their recombinant production. The low yields of the soluble form of the trimeric spike (S) glycoprotein from SARS-CoV-2 has prompted advances in protein engineering that have greatly enhanced the stability and yields of the glycoprotein. The latest expression-enhanced version of the spike incorporates six proline substitutions to stabilize the prefusion conformation (termed SARS-CoV-2 S HexaPro). Although the substitutions greatly enhanced expression whilst not compromising protein structure, the influence of these substitutions on glycan processing has not been explored. Here, we show that the site-specific N-linked glycosylation of the expression-enhanced HexaPro resembles that of an earlier version containing two proline substitutions (2P), and that both capture features of native viral glycosylation. However, there are site-specific differences in glycosylation of HexaPro when compared to 2P. Despite these discrepancies, analysis of the serological reactivity of clinical samples from infected individuals confirmed that both HexaPro and 2P protein are equally able to detect IgG, IgA, and IgM responses in all sera analysed. Moreover, we extend this observation to include an analysis of glycan engineered S protein, whereby all N-linked glycans were converted to oligomannose-type and conclude that serological activity is not impacted by large scale changes in glycosylation. These observations suggest that variations in glycan processing will not impact the serological assessments currently being performed across the globe., Competing Interests: Declaration of interests J.S.M. is an inventor on the following U.S. patent applications: no. 62/412,703 (“Prefusion Coronavirus Spike Proteins and Their Use”); no. 62/972,886 (“2019-nCoV Vaccine”); no. 63/032,502 (“Engineered Coronavirus Spike (S) Protein and Methods of Use Thereof”)., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

A central tenet in the design of vaccines is the display of native-like antigens in the elicitation of protective immunity. The abundance of N-linked glycans across the SARS-CoV-2 spike protein is a potential source of heterogeneity among the many different vaccine candidates under investigation. Here, we investigate the glycosylation of recombinant SARS-CoV-2 spike proteins from five different laboratories and compare them against S protein from infectious virus, cultured in Vero cells. We find patterns that are conserved across all samples, and this can be associated with site-specific stalling of glycan maturation that acts as a highly sensitive reporter of protein structure. Molecular dynamics simulations of a fully glycosylated spike support a model of steric restrictions that shape enzymatic processing of the glycans. These results suggest that recombinant spike-based SARS-CoV-2 immunogen glycosylation reproducibly recapitulates signatures of viral glycosylation.

The spike (S) protein is the main handle for SARS-CoV-2 to enter host cells via surface angiotensin-converting enzyme 2 (ACE2) receptors. How ACE2 binding activates proteolysis of S protein is unknown. Here, using amide hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations, we have mapped the S:ACE2 interaction interface and uncovered long-range allosteric propagation of ACE2 binding to sites necessary for host-mediated proteolysis of S protein, critical for viral host entry. Unexpectedly, ACE2 binding enhances dynamics at a distal S1/S2 cleavage site and flanking protease docking site ~27 Å away while dampening dynamics of the stalk hinge (central helix and heptad repeat [HR]) regions ~130 Å away. This highlights that the stalk and proteolysis sites of the S protein are dynamic hotspots in the prefusion state. Our findings provide a dynamics map of the S:ACE2 interface in solution and also offer mechanistic insights into how ACE2 binding is allosterically coupled to distal proteolytic processing sites and viral-host membrane fusion. Thus, protease docking sites flanking the S1/S2 cleavage site represent alternate allosteric hotspot targets for potential therapeutic development., Competing Interests: PR, NT, FS, XQ, KP, GY, MK, WH, JL, PB, PM, GA No competing interests declared, (© 2021, Raghuvamsi et al.)

The high mutation rate in retroviruses is one of the leading causes of drug resistance. In human immunodeficiency virus type-1 (HIV-1), synergistic mutations in its protease and the protease substrate - the Group-specific antigen (Gag) polyprotein - work together to confer drug resistance against protease inhibitors and compensate the mutations affecting viral fitness. Some Gag mutations can restore Gag-protease binding, yet most Gag-protease correlated mutations occur outside of the Gag cleavage site. To investigate the molecular basis for this, we now report multiscale modelling approaches to investigate various sequentially cleaved Gag products in the context of clinically relevant mutations that occur outside of the cleavage sites, including simulations of the largest Gag proteolytic product in its viral membrane-bound state. We found that some mutations, such as G123E and H219Q, involve direct interaction with cleavage site residues to influence their local environment, while certain mutations in the matrix domain lead to the enrichment of lipids important for Gag targeting and assembly. Collectively, our results reveal why non-cleavage site mutations have far-reaching implications outside of Gag proteolysis, with important consequences for drugging Gag maturation intermediates and tackling protease inhibitor resistance., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

Peptide oligomerization dynamics affects peptide structure, activity, and pharmacodynamic properties. The thrombin C-terminal peptide, TCP-25 (GKYGFYTHVFRLKKWIQKVIDQFGE), is currently in preclinical development for improved wound healing and infection prevention. It exhibits turbidity when formulated at pH 7.4, particularly at concentrations of 0.3 mM or more. We used biochemical and biophysical approaches to explore whether the peptide self-associates and forms oligomers. The peptide showed a dose-dependent increase in turbidity as well as α-helical structure at pH 7.4, a phenomenon not observed at pH 5.0. By analyzing the intrinsic tryptophan fluorescence, we demonstrate that TCP-25 is more stable at high concentrations (0.3 mM) when exposed to high temperatures or a high concentration of denaturant agents, which is compatible with oligomer formation. The denaturation process was reversible above 100 µM of peptide. Dynamic light scattering demonstrated that TCP-25 oligomerization is sensitive to changes in pH, time, and temperature. Computational modeling with an active 18-mer region of TCP-25 showed that the peptide can form pH-dependent higher-order end-to-end oligomers and micelle-like structures, which is in agreement with the experimental data. Thus, TCP-25 exhibits pH- and temperature-dependent dynamic changes involving helical induction and reversible oligomerization, which explains the observed turbidity of the pharmacologically developed formulation.

There is a link between high lipopolysaccharide (LPS) levels in the blood and the metabolic syndrome, and metabolic syndrome predisposes patients to severe COVID-19. Here, we define an interaction between SARS-CoV-2 spike (S) protein and LPS, leading to aggravated inflammation in vitro and in vivo. Native gel electrophoresis demonstrated that SARS-CoV-2 S protein binds to LPS. Microscale thermophoresis yielded a KD of ∼47 nM for the interaction. Computational modeling and all-atom molecular dynamics simulations further substantiated the experimental results, identifying a main LPS-binding site in SARS-CoV-2 S protein. S protein, when combined with low levels of LPS, boosted nuclear factor-kappa B (NF-κB) activation in monocytic THP-1 cells and cytokine responses in human blood and peripheral blood mononuclear cells, respectively. The in vitro inflammatory response was further validated by employing NF-κB reporter mice and in vivo bioimaging. Dynamic light scattering, transmission electron microscopy, and LPS-FITC analyses demonstrated that S protein modulated the aggregation state of LPS, providing a molecular explanation for the observed boosting effect. Taken together, our results provide an interesting molecular link between excessive inflammation during infection with SARS-CoV-2 and comorbidities involving increased levels of bacterial endotoxins., (© The Author(s) (2020). Published by Oxford University Press on behalf of Journal of Molecular Cell Biology, IBCB, SIBS, CAS.)

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

Coordination of outer membrane constriction with septation is critical to faithful division in Gram-negative bacteria and vital to the barrier function of the membrane. This coordination requires the recruitment of the peptidoglycan-binding outer-membrane lipoprotein Pal at division sites by the Tol system. Here, we show that Pal accumulation at Escherichia coli division sites is a consequence of three key functions of the Tol system. First, Tol mobilises Pal molecules in dividing cells, which otherwise diffuse very slowly due to their binding of the cell wall. Second, Tol actively captures mobilised Pal molecules and deposits them at the division septum. Third, the active capture mechanism is analogous to that used by the inner membrane protein TonB to dislodge the plug domains of outer membrane TonB-dependent nutrient transporters. We conclude that outer membrane constriction is coordinated with cell division by active mobilisation-and-capture of Pal at division septa by the Tol system.

The therapeutic potential of immunoglobulin M (IgM) is of considerable interest in immunotherapy due to its complement-activating and cell-agglutinating abilities. Pertuzumab and Trastuzumab are monoclonal antibodies used to treat human epidermal growth factor receptor 2 (HER2)-positive breast cancer but exhibit significantly different binding affinities as IgM when compared to its IgG isotype. Using integrative multiscale modelling and simulations of complete antibody assemblies, we show that Pertuzumab IgM is able to utilize all of its V-regions to bind multiple HER2 receptors simultaneously, while similar binding in Trastuzumab IgM is prohibited by steric clashes caused by the large globular domain of HER2. This is subsequently validated by confirming that Pertuzumab IgM inhibits proliferation in HER2 over-expressing live cells more effectively than its IgG counterpart and Trastuzumab IgM. Our study highlights the importance of understanding the molecular details of antibody-antigen interactions for the design and isotype selection of therapeutic antibodies., (This journal is © The Royal Society of Chemistry 2020.)

The outer membrane channel OprD from Pseudomonas aeruginosa transports basic amino acids and clinically relevant carbapenem antibiotics. Understanding the molecular basis of substrate permeation across this channel will therefore lead to better therapeutic designs to treat infections. Using umbrella sampling simulations, we calculated the potential of mean force (PMF) for the arginine permeation pathway through OprD. The PMF reveals a deep free energy well of ∼6 kT around the putative substrate binding site followed by a shallower well of ∼4 kT close to the most constricted region of the pore. Despite becoming partially dehydrated during translocation, some water molecules are retained to shield the guanidinium side chain of arginine from the ladder of basic residues in the protein. Sugars of the lipopolysaccharide headgroups form contacts with arginine and could potentially play an important role in transferring substrate from the external medium to OprD. The PMF through bulk membrane shows a large energetic barrier of ∼45 kT within the hydrophobic core of the membrane, suggesting that spontaneous translocation without OprD is highly unlikely. This significant energetic penalty is likely caused by the extensive distortion of the lower leaflet of the outer membrane as phospholipid headgroups sink inwards to interact with charged groups of arginine. Our results provide quantitative insights into solute permeation across the bacterial outer membrane.

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., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

The cell envelope of Gram-negative bacteria is synthesized and maintained via mechanisms that are targets for development of novel antibiotics. Here we focus on the process of moving Braun's lipoprotein (BLP) from the periplasmic space to the outer membrane of E. coli, via the LolA protein. In contrast to current thinking, we show that binding of multiple inhibitor molecules inside the hydrophobic cavity of LolA does not prevent subsequent binding of BLP inside the same cavity. Rather, based on our atomistic simulations we propose the theory that once inhibitors and BLP are bound inside the cavity of LolA, driven by hydrophobic interactions, they become entangled with each other. Our umbrella sampling calculations show that on the basis of energetics, it is more difficult to dislodge BLP from the cavity of LolA when it is uncomplexed compared to complexed with inhibitor. Thus the inhibitor reduces the affinity of BLP for the LolA cavity.

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.

Lipid membranes play a crucial role in living systems by compartmentalizing biological processes and forming a barrier between these processes and the environment. Naturally, a large apparatus of biomolecules is responsible for construction, maintenance, transport, and degradation of these lipid barriers. Additional classes of biomolecules are tasked with transport of specific substances or transduction of signals from the environment across lipid membranes. In this article, we intend to describe a set of techniques that enable one to build accurate models of lipid systems and their associated proteins, and to simulate their dynamics over a variety of time and length scales. We discuss the methods and challenges that allow us to derive structural, mechanistic, and thermodynamic information from these models. We also show how these models have recently been applied in research to study some of the most complex lipid-protein systems to date, including bacterial and viral envelopes, neuronal membranes, and mammalian signaling systems.

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 F 1 F O 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., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Protein bacteriocins are potent narrow spectrum antibiotics that exploit outer membrane porins to kill bacteria by poorly understood mechanisms. Here, we determine how colicins, bacteriocins specific for Escherichia coli, engage the trimeric porin OmpF to initiate toxin entry. The N-terminal ∼80 residues of the nuclease colicin ColE9 are intrinsically unstructured and house two OmpF binding sites (OBS1 and OBS2) that reside within the pores of OmpF and which flank an epitope that binds periplasmic TolB. Using a combination of molecular dynamics simulations, chemical trimerization, isothermal titration calorimetry, fluorescence microscopy, and single channel recording planar lipid bilayer measurements, we show that this arrangement is achieved by OBS2 binding from the extracellular face of OmpF, while the interaction of OBS1 occurs from the periplasmic face of OmpF. Our study shows how the narrow pores of oligomeric porins are exploited by colicin disordered regions for direction-specific binding, which ensures the constrained presentation of an activating signal within the bacterial periplasm.

Bacterial peptidoglycan glycosyltransferases (PGT) catalyse the essential polymerization of lipid II into linear glycan chains required for peptidoglycan biosynthesis. The PGT domain is composed of a large head subdomain and a smaller jaw subdomain and can be potently inhibited by the antibiotic moenomycin A (MoeA). We present an X-ray structure of the MoeA-bound Staphylococcus aureus monofunctional PGT enzyme, revealing electron density for a second MoeA bound to the jaw subdomain as well as the PGT donor site. Isothermal titration calorimetry confirms two drug-binding sites with markedly different affinities and positive cooperativity. Hydrophobic cluster analysis suggests that the membrane-interacting surface of the jaw subdomain has structural and physicochemical properties similar to amphipathic cationic α -helical antimicrobial peptides for lipid II recognition and binding. Furthermore, molecular dynamics simulations of the drug-free and -bound forms of the enzyme demonstrate the importance of the jaw subdomain movement for lipid II selection and polymerization process and provide molecular-level insights into the mechanism of peptidoglycan biosynthesis by PGTs.

The cell envelope of Gram-negative bacteria is composed of two membranes separated by a soluble region. Here, we report microsecond time scale coarse-grained molecular dynamics simulations of models of the Escherichia coli cell envelope that incorporate both membranes and various native membrane proteins. Our results predict that both the inner and outer membranes curve in a manner dependent on the size of the embedded proteins. The tightly cross-linked lipopolysaccharide molecules (LPS) of the outer membrane cause a strong coupling between the movement of proteins and lipids. While the flow of phospholipids is more random, their diffusion is nevertheless influenced by nearby proteins. Our results reveal protein-induced lipid sorting, whereby cardiolipin is significantly enriched within the vicinity of the water channel AqpZ and the multidrug efflux pump AcrBZ. In summary, our results provide unprecedented details of the intricate relationship between both membranes of E. coli and the proteins embedded within them.

Gram-negative bacteria such as Escherichia coli are protected by a complex cell envelope. The development of novel therapeutics against these bacteria necessitates a molecular level understanding of the structure-dynamics-function relationships of the various components of the cell envelope. We use atomistic MD simulations to reveal the details of covalent and noncovalent protein interactions that link the outer membrane to the aqueous periplasmic region. We show that the Braun's lipoprotein tilts and bends, and thereby lifts the cell wall closer to the outer membrane. Both monomers and dimers of the outer membrane porin OmpA can interact with peptidoglycan in the presence of Braun's lipoprotein, but in the absence of the latter, only dimers of OmpA show a propensity to form contacts with peptidoglycan. Our study provides a glimpse of how the molecular components of the bacterial cell envelope interact with each other to mediate cell wall attachment in E. coli., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Bacteria are protected by complex molecular architectures known as the cell envelope. The cell envelope is composed of regions with distinct chemical compositions and physical properties, namely, membranes and a cell wall. To develop novel antibiotics to combat pathogenic bacteria, molecular level knowledge of the structure, dynamics, and interplay between the chemical components of the cell envelope that surrounds bacterial cells is imperative. In addition, conserved molecular patterns associated with the bacterial envelope are recognized by receptors as part of the mammalian defensive response to infection, and an improved understanding of bacteria-host interactions would facilitate the search for novel immunotherapeutics. This Perspective introduces an emerging area of computational biology: multiscale molecular dynamics simulations of chemically complex models of bacterial lipids and membranes. We discuss progress to date, and identify areas for future development that will enable the study of aspects of the membrane components that are as yet unexplored by computational methods.

The envelope of Gram-negative bacteria is highly complex, containing separate outer and inner membranes and an intervening periplasmic space encompassing a peptidoglycan (PGN) cell wall. The PGN scaffold is anchored non-covalently to the outer membrane via globular OmpA-like domains of various proteins. We report atomically detailed simulations of PGN bound to OmpA in three different states, including the isolated C-terminal domain (CTD), the full-length monomer, or the complete full-length dimeric form. Comparative analysis of dynamics of OmpA CTD from different bacteria helped to identify a conserved PGN-binding mode. The dynamics of full-length OmpA, embedded within a realistic representation of the outer membrane containing full-rough (Ra) lipopolysaccharide, phospholipids, and cardiolipin, suggested how the protein may provide flexible mechanical support to the cell wall. An accurate model of the heterogeneous bacterial cell envelope should facilitate future efforts to develop antibacterial agents., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

OmpA is a multidomain protein found in the outer membranes of most Gram-negative bacteria. Despite a wealth of reported structural and biophysical studies, the structure-function relationships of this protein remain unclear. For example, it is still debated whether it functions as a pore, and the precise molecular role it plays in attachment to the peptidoglycan of the periplasm is unknown. The absence of a consensus view is partly due to the lack of a complete structure of the full-length protein. To address this issue, we performed molecular-dynamics simulations of the full-length model of the OmpA dimer proposed by Robinson and co-workers. The N-terminal domains were embedded in an asymmetric model of the outer membrane, with lipopolysaccharide molecules in the outer leaflet and phospholipids in the inner leaflet. Our results reveal a large dimerization interface within the membrane environment, ensuring that the dimer is stable over the course of the simulations. The linker is flexible, expanding and contracting to pull the globular C-terminal domain up toward the membrane or push it down toward the periplasm, suggesting a possible mechanism for providing mechanical stability to the cell. The external loops were more stabilized than was observed in previous studies due to the extensive dimerization interface and presence of lipopolysaccharide molecules in our outer-membrane model, which may have functional consequences in terms of OmpA adhesion to host cells. In addition, the pore-gating behavior of the protein was modulated compared with previous observations, suggesting a possible role for dimerization in channel regulation., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Membrane transporters are critical modulators of drug pharmacokinetics, efficacy, and safety. One example is the proton-dependent oligopeptide transporter PepT1, also known as SLC15A1, which is responsible for the uptake of the ?-lactam antibiotics and various peptide-based prodrugs. In this study, we modeled the binding of various peptides to a bacterial homolog, PepT(St), and evaluated a range of computational methods for predicting the free energy of binding. Our results show that a hybrid approach (endpoint methods to classify peptides into good and poor binders and a theoretically exact method for refinement) is able to accurately predict affinities, which we validated using proteoliposome transport assays. Applying the method to a homology model of PepT1 suggests that the approach requires a high-quality structure to be accurate. Our study provides a blueprint for extending these computational methodologies to other pharmaceutically important transporter families., (Copyright © 2016 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Mammals obtain nitrogen via the uptake of di- and tri-peptides in the gastrointestinal tract through the action of PepT1 and PepT2, which are members of the POT family of proton-coupled oligopeptide transporters. PepT1 and PepT2 also play an important role in drug transport in the human body. Recent crystal structures of bacterial homologs revealed a conserved peptide-binding site and mechanism of transport. However, a key structural difference exists between bacterial and mammalian homologs with only the latter containing a large extracellular domain, the function of which is currently unknown. Here, we present the crystal structure of the extracellular domain from both PepT1 and PepT2 that reveal two immunoglobulin-like folds connected in tandem, providing structural insight into mammalian peptide transport. Functional and biophysical studies demonstrate that these domains interact with the intestinal protease trypsin, suggesting a role in clustering proteolytic activity to the site of peptide transport in eukaryotic cells., (Copyright © 2015 The Authors. Published by Elsevier Ltd.. All rights reserved.)