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

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

Antibody affinity versus dengue morphology influences neutralization.

Different strains within a dengue serotype (DENV1-4) can have smooth, or "bumpy" surface morphologies with different antigenic characteristics at average body temperature (37°C). We determined the neutralizing properties of a serotype cross-reactive human monoclonal antibody (HMAb) 1C19 for strains with differing morphologies within the DENV1 and DENV2 serotypes. We mapped the 1C19 epitope to E protein domain II by hydrogen deuterium exchange mass spectrometry, cryoEM and molecular dynamics simulations, revealing that this epitope is likely partially hidden on the virus surface. We showed the antibody has high affinity for binding to recombinant DENV1 E proteins compared to those of DENV2, consistent with its strong neutralizing activities for all DENV1 strains tested regardless of their morphologies. This finding suggests that the antibody could out-compete E-to-E interaction for binding to its epitope. In contrast, for DENV2, HMAb 1C19 can only neutralize when the epitope becomes exposed on the bumpy-surfaced particle. Although HMAb 1C19 is not a suitable therapeutic candidate, this study with HMAb 1C19 shows the importance of choosing a high-affinity antibody that could neutralize diverse dengue virus morphologies for therapeutic purposes.

Computational modelling of flavivirus dynamics: The ins and outs.

Enveloped viruses such as the flaviviruses represent a significant burden to human health around the world, with hundreds of millions of people each year affected by dengue alone. In an effort to improve our understanding of the molecular basis for the infective mechanisms of these viruses, extensive computational modelling approaches have been applied to elucidate their conformational dynamics. Multiscale protocols have been developed to simulate flavivirus envelopes in close accordance with biophysical data, in particular derived from cryo-electron microscopy, enabling high-resolution refinement of their structures and elucidation of the conformational changes associated with adaptation both to host environments and to immunological factors such as antibodies. Likewise, integrative modelling efforts combining data from biophysical experiments and from genome sequencing with chemical modification are providing unparalleled insights into the architecture of the previously unresolved nucleocapsid complex. Collectively, this work provides the basis for the future rational design of new antiviral therapeutics and vaccine development strategies targeting enveloped viruses.

Molecular basis of dengue virus serotype 2 morphological switch from 29°C to 37°C.

The ability of DENV2 to display different morphologies (hence different antigenic properties) complicates vaccine and therapeutics development. Previous studies showed most strains of laboratory adapted DENV2 particles changed from smooth to "bumpy" surfaced morphology when the temperature is switched from 29°C at 37°C. Here we identified five envelope (E) protein residues different between two alternative passage history DENV2 NGC strains exhibiting smooth or bumpy surface morphologies. Several mutations performed on the smooth DENV2 infectious clone destabilized the surface, as observed by cryoEM. Molecular dynamics simulations demonstrated how chemically subtle substitution at various positions destabilized dimeric interactions between E proteins. In contrast, three out of four DENV2 clinical isolates showed a smooth surface morphology at 37°C, and only at high fever temperature (40°C) did they become "bumpy". These results imply vaccines should contain particles representing both morphologies. For prophylactic and therapeutic treatments, this study also informs on which types of antibodies should be used at different stages of an infection, i.e., those that bind to monomeric E proteins on the bumpy surface or across multiple E proteins on the smooth surfaced virus.

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

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

Extending the Martini Coarse-Grained Force Field to N-Glycans.

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.

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

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

The architecture of the OmpC-MlaA complex sheds light on the maintenance of outer membrane lipid asymmetry in Escherichia coli.

A distinctive feature of the Gram-negative bacterial cell envelope is the asymmetric outer membrane (OM), where lipopolysaccharides and phospholipids (PLs) reside in the outer and inner leaflets, respectively. This unique lipid asymmetry renders the OM impermeable to external insults, including antibiotics and bile salts. In Escherichia coli, the complex comprising osmoporin OmpC and the OM lipoprotein MlaA is believed to maintain lipid asymmetry by removing mislocalized PLs from the outer leaflet of the OM. How this complex performs this function is unknown. Here, we defined the molecular architecture of the OmpC-MlaA complex to gain insights into its role in PL transport. Using in vivo photo-cross-linking and molecular dynamics simulations, we established that MlaA interacts extensively with OmpC and is located entirely within the lipid bilayer. In addition, MlaA forms a hydrophilic channel, likely enabling PL translocation across the OM. We further showed that flexibility in a hairpin loop adjacent to the channel is critical in modulating MlaA activity. Finally, we demonstrated that OmpC plays a functional role in maintaining OM lipid asymmetry together with MlaA. Our work offers glimpses into how the OmpC-MlaA complex transports PLs across the OM and has important implications for future antibacterial drug development.

Multiscale modeling of innate immune receptors: Endotoxin recognition and regulation by host defense peptides.

The innate immune system provides a first line of defense against foreign microorganisms, and is typified by the Toll-like receptor (TLR) family. TLR4 is of particular interest, since over-stimulation of its pathway by excess lipopolysaccharide (LPS) molecules from the outer membranes of Gram-negative bacteria can result in sepsis, which causes millions of deaths each year. In this review, we outline our use of molecular simulation approaches to gain a better understanding of the determinants of LPS recognition, towards the search for novel immunotherapeutics. We first describe how atomic-resolution simulations have enabled us to elucidate the regulatory conformational changes in TLR4 associated with different LPS analogues, and hence a means to rationalize experimental structure-activity data. Furthermore, multiscale modelling strategies have provided a detailed description of the thermodynamics and intermediate structures associated with the entire TLR4 relay - which consists of a number of transient receptor/coreceptor complexes - allowing us trace the pathway of LPS transfer from bacterial membranes to the terminal receptor complex at the plasma membrane surface. Finally, we describe our efforts to leverage these computational models, in order to elucidate previously undisclosed anti-inflammatory mechanisms of endogenous host-defense peptides found in wounds. Collectively, this work represents a promising avenue for the development of novel anti-septic treatments, inspired by nature's innate defense strategies.

A polar SxxS motif drives assembly of the transmembrane domains of Toll-like receptor 4.

Like all members of the Toll-like receptor (TLR) family, TLR4 comprises of a large ectodomain (ECD) involved in ligand recognition at the cell-surface, and a cytosolic Toll/interleukin-1 receptor (TIR) signalling domain, linked by a lipid membrane-anchored transmembrane (TM) domain (TMD). Binding of immunostimulatory pathogen-associated molecular patterns (PAMPs) such as bacterial lipopolysaccharide (LPS) to myeloid differentiation factor 2 (MD-2) coreceptor-complexed TLR4 leads to its dimerization, resulting in productive juxtaposition of TIR domains and the initiation of pro-inflammatory innate immune responses. Whilst the process of PAMP recognition is relatively well understood, thanks to numerous high-resolution crystallographic structures of ECDs, the mechanism by which such recognition is translated into TMD dimerization and activating conformational changes is less clear. Based on available biophysical and biochemical experimental data, ab initio modelling, and multiscale molecular dynamics (MD) simulations entailing a total of >13µs and >200µs of atomistic and coarse-grained sampling, respectively, we investigate the structural basis for TLR4 TMD dimerization within a biologically relevant lipid membrane environment. A key polar-xx-polar (637SxxS640) motif is shown to drive association of the TLR4 TMDs, and to maintain a flexible interface, which may be disrupted by selected point mutations. Furthermore, MD simulations of various TMD+ECD constructs have been used to investigate the coupling between domains, revealing that flexible linkers abrogate dimerization via aggregation of ECDs at the membrane surface, explaining previous biochemical observations. These results improve our understanding of the assembly and signalling mechanisms of TLR4, and pave the way for rational structure-based development of membrane-associated immunomodulatory molecules.

South-east Asian Zika virus strain linked to cluster of cases in Singapore, August 2016.

Zika virus (ZIKV) is an ongoing global public health emergency with 70 countries and territories reporting evidence of ZIKV transmission since 2015. On 27 August 2016, Singapore reported its first case of local ZIKV transmission and identified an ongoing cluster. Here, we report the genome sequences of ZIKV strains from two cases and find through phylogenetic analysis that these strains form an earlier branch distinct from the recent large outbreak in the Americas.

Free energy predictions of ligand binding to an α-helix using steered molecular dynamics and umbrella sampling simulations.

Free energy prediction of ligand binding to macromolecules using explicit solvent molecular dynamics (MD) simulations is computationally very expensive. Recently, we reported a linear correlation between the binding free energy obtained via umbrella sampling (US) versus the rupture force from steered molecular dynamics (SMD) simulations for epigallocatechin-3-gallate (EGCG) binding to α-helical-rich keratin. This linear correlation suggests a potential route for fast free energy predictions using SMD alone. In this work, the generality of the linear correlation is further tested for several ligands interacting with the α-helical motif of keratin. These molecules have significantly varying properties, i.e., octanol/water partition coefficient (log P), and/or overall charges (oleic acid, catechin, Fe(2+), citric acid, hydrogen citrate, dihydrogen citrate, and citrate). Using the constant loading rate of our previous study of the keratin-EGCG system, we observe that the linear correlation for keratin-EGCG can be extended to other uncharged molecules where interactions are governed by hydrogen bonds and/or a combination of hydrogen bonds and hydrophobic forces. For molecules where interactions with the keratin helix are governed primarily by electrostatics between charged molecules, a second, alternative linear correlation model is derived. While further investigations are needed to expand the molecular space and build a fully predictive model, the current approach represents a promising methodology for fast free energy predictions based on short SMD simulations (requiring picoseconds to nanoseconds of sampling) for defined biomolecular systems.

Mechanisms of allostery at the viral surface through the eyes of molecular simulation.

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.

Coarse-Grained Model of Glycosaminoglycans for Biomolecular Simulations.

Proteoglycans contain glycosaminoglycans (GAGs) which are negatively charged linear polymers made of repeating disaccharide units of uronic acid and hexosamine units. They play vital roles in numerous physiological and pathological processes, particularly in governing cellular communication and attachment. Depending on their sulfonation state, acetylation, and glycosidic linkages, GAGs belong to different families. The high molecular weight, heterogeneity, and flexibility of GAGs hamper their characterization at atomic resolution, but this may be circumvented via coarse-grained (CG) approaches. In this work, we report a CG model for a library of common GAG types in their isolated or proteoglycan-linked states compatible with version 2.2 (v2.2) of the widely popular CG Martini force field. The model reproduces conformational and thermodynamic properties for a wide variety of GAGs, as well as matching structural and binding data for selected proteoglycan test systems. The parameters developed here may thus be employed to study a range of GAG-containing biomolecular systems, thereby benefiting from the efficiency and broad applicability of the Martini framework.

Impact on S. aureus and E. coli Membranes of Treatment with Chlorhexidine and Alcohol Solutions: Insights from Molecular Simulations and Nuclear Magnetic Resonance.

Membranes form the first line of defence of bacteria against potentially harmful molecules in the surrounding environment. Understanding the protective properties of these membranes represents an important step towards development of targeted anti-bacterial agents such as sanitizers. Use of propanol, isopropanol and chlorhexidine can significantly decrease the threat imposed by bacteria in the face of growing anti-bacterial resistance via mechanisms that include membrane disruption. Here we have employed molecular dynamics simulations and nuclear magnetic resonance to explore the impact of chlorhexidine and alcohol on the S. aureus cell membrane, as well as the E. coli inner and outer membranes. We identify how sanitizer components partition into these bacterial membranes, and show that chlorhexidine is instrumental in this process.

A pH-dependent cluster of charges in a conserved cryptic pocket on flaviviral envelopes.

Flaviviruses are enveloped viruses which include human pathogens that are predominantly transmitted by mosquitoes and ticks. Some, such as dengue virus, exhibit the phenomenon of antibody-dependent enhancement (ADE) of disease, making vaccine-based routes of fighting infections problematic. The pH-dependent conformational change of the envelope (E) protein required for fusion between the viral and endosomal membranes is an attractive point of inhibition by antivirals as it has the potential to diminish the effects of ADE. We examined six flaviviruses by employing large-scale molecular dynamics (MD) simulations of raft systems that represent a substantial portion of the flaviviral envelope. We utilised a benzene-mapping approach that led to a discovery of shared hotspots and conserved cryptic sites. A cryptic pocket previously shown to bind a detergent molecule exhibited strain-specific characteristics. An alternative conserved cryptic site at the E protein domain interfaces showed a consistent dynamic behaviour across flaviviruses and contained a conserved cluster of ionisable residues. Constant-pH simulations revealed cluster and domain-interface disruption under low pH conditions. Based on this, we propose a cluster-dependent mechanism that addresses inconsistencies in the histidine-switch hypothesis and highlights the role of cluster protonation in orchestrating the domain dissociation pivotal for the formation of the fusogenic trimer.

Antibacterial and Anti-Inflammatory Effects of Apolipoprotein E.

Apolipoprotein E (APOE) is a lipid-transport protein that functions as a key mediator of lipid transport and cholesterol metabolism. Recent studies have shown that peptides derived from human APOE display anti-inflammatory and antimicrobial effects. Here, we applied in vitro assays and fluorescent microscopy to investigate the anti-bacterial effects of full-length APOE. The interaction of APOE with endotoxins from Escherichia coli was explored using surface plasmon resonance, binding assays, transmission electron microscopy and all-atom molecular dynamics (MD) simulations. We also studied the immunomodulatory activity of APOE using in vitro cell assays and an in vivo mouse model in combination with advanced imaging techniques. We observed that APOE exhibits anti-bacterial activity against several Gram-negative bacterial strains of Pseudomonas aeruginosa and Escherichia coli. In addition, we showed that APOE exhibits a significant binding affinity for lipopolysaccharide (LPS) and lipid A as well as heparin. MD simulations identified the low-density lipoprotein receptor (LDLR) binding region in helix 4 of APOE as a primary binding site for these molecules via electrostatic interactions. Together, our data suggest that APOE may have an important role in controlling inflammation during Gram-negative bacterial infection.

Motional clustering in supra-τc conformational exchange influences NOE cross-relaxation rate.

Biomolecular spin relaxation processes, such as the NOE, are commonly modeled by rotational τc-tumbling combined with fast motions on the sub-τc timescale. Motions on the supra-τc timescale, in contrast, are considered to be completely decorrelated to the molecular tumbling and therefore invisible. Here, we show how supra-τc dynamics can nonetheless influence the NOE build-up between methyl groups. This effect arises because supra-τc motions can cluster the fast-motion ensembles into discrete states, affecting distance averaging as well as the fast-motion order parameter and hence the cross-relaxation rate. We present a computational approach to estimate methyl-methyl cross-relaxation rates from extensive (>100×τc) all-atom molecular dynamics (MD) trajectories on the example of the 723-residue protein Malate Synthase G. The approach uses Markov state models (MSMs) to resolve transitions between metastable states and thus to discriminate between sub-τc and supra-τc conformational exchange. We find that supra-τc exchange typically increases NOESY cross-peak intensities. The methods described in this work extend the theory of modeling sub-µs dynamics in spin relaxation and thus contribute to a quantitative estimation of NOE cross-relaxation rates from MD simulations, eventually leading to increased precision in structural and functional studies of large proteins.

Bridging the N-terminal and middle domains in FliG of the flagellar rotor.

Flagella are necessary for bacterial movement and contribute to various aspects of virulence. They are complex cylindrical structures built of multiple molecular rings with self-assembly properties. The flagellar rotor is composed of the MS-ring and the C-ring. The FliG protein of the C-ring is central to flagellar assembly and function due to its roles in linking the C-ring with the MS-ring and in torque transmission from stator to rotor. No high-resolution structure of an assembled C-ring has been resolved to date, and the conformation adopted by FliG within the ring is unclear due to variations in available crystallographic data. Here, we use molecular dynamics (MD) simulations to study the conformation and dynamics of FliG in different states of assembly, including both in physiologically relevant and crystallographic lattice environments. We conclude that the linker between the FliG N-terminal and middle domain likely adopts an extended helical conformation in vivo, in contrast with the contracted conformation observed in some previous X-ray studies. We further support our findings with integrative model building of full-length FliG and a FliG ring model that is compatible with cryo-electron tomography (cryo-ET) and electron microscopy (EM) densities of the C-ring. Collectively, our study contributes to a better mechanistic understanding of the flagellar rotor assembly and its function.

Uncovering cryptic pockets in the SARS-CoV-2 spike glycoprotein.

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.