Titratable residues that drive RND efflux: Insights from molecular simulations.

The resistance-nodulation-division efflux machinery confers antimicrobial resistance to Gram-negative bacteria by actively pumping antibiotics out of the cell. The protein complex is powered by proton motive force; however, the proton transfer mechanism itself and indeed even its stoichiometry is still unclear. Here we review computational studies from the last decade that focus on elucidating the number of protons transferred per conformational cycle of the pump. Given the difficulties in studying proton movement using even state-of-the-art structural biology methods, the contributions from computational studies have been invaluable from a mechanistic perspective.

Molecular Crowding Alters the Interactions of Polymyxin Lipopeptides within the Periplasm of E. coli: Insights from Molecular Dynamics.

The cell envelope of Gram-negative bacteria is a crowded tripartite architecture that separates the cell interior from the external environment. Two membranes encapsulate the aqueous periplasm, which contains the cell wall. Little is known about the mechanisms via which antimicrobial peptides move through the periplasm from the outer membrane to their site of action, the inner membrane. We utilize all-atom molecular dynamics to study two antimicrobial peptides, polymyxins B1 and E, within models of the E. coli periplasm crowded to different extents. In a simple chemical environment, both PMB1 and PME bind irreversibly to the cell wall. The presence of specific macromolecules leads to competition with the polymyxins for cell wall interaction sites, resulting in polymyxin dissociation from the cell wall. Chemical complexity also impacts interactions between polymyxins and Braun's lipoprotein; thus, the interaction modes of lipoprotein antibiotics within the periplasm are dependent upon the nature of the other species present.

Activation and substrate specificity of the human P4-ATPase ATP8B1.

Asymmetric distribution of phospholipids in eukaryotic membranes is essential for cell integrity, signaling pathways, and vesicular trafficking. P4-ATPases, also known as flippases, participate in creating and maintaining this asymmetry through active transport of phospholipids from the exoplasmic to the cytosolic leaflet. Here, we present a total of nine cryo-electron microscopy structures of the human flippase ATP8B1-CDC50A complex at 2.4 to 3.1 Å overall resolution, along with functional and computational studies, addressing the autophosphorylation steps from ATP, substrate recognition and occlusion, as well as a phosphoinositide binding site. We find that the P4-ATPase transport site is occupied by water upon phosphorylation from ATP. Additionally, we identify two different autoinhibited states, a closed and an outward-open conformation. Furthermore, we identify and characterize the PI(3,4,5)P3 binding site of ATP8B1 in an electropositive pocket between transmembrane segments 5, 7, 8, and 10. Our study also highlights the structural basis of a broad lipid specificity of ATP8B1 and adds phosphatidylinositol as a transport substrate for ATP8B1. We report a critical role of the sn-2 ester bond of glycerophospholipids in substrate recognition by ATP8B1 through conserved S403. These findings provide fundamental insights into ATP8B1 catalytic cycle and regulation, and substrate recognition in P4-ATPases.

Dual function of OmpM as outer membrane tether and nutrient uptake channel in diderm Firmicutes.

The outer membrane (OM) in diderm, or Gram-negative, bacteria must be tethered to peptidoglycan for mechanical stability and to maintain cell morphology. Most diderm phyla from the Terrabacteria group have recently been shown to lack well-characterised OM attachment systems, but instead have OmpM, which could represent an ancestral tethering system in bacteria. Here, we have determined the structure of the most abundant OmpM protein from Veillonella parvula (diderm Firmicutes) by single particle cryogenic electron microscopy. We also characterised the channel properties of the transmembrane ß-barrel of OmpM and investigated the structure and PG-binding properties of its periplasmic stalk region. Our results show that OM tethering and nutrient acquisition are genetically linked in V. parvula, and probably other diderm Terrabacteria. This dual function of OmpM may have played a role in the loss of the OM in ancestral bacteria and the emergence of monoderm bacterial lineages.

Phospholipid transport to the bacterial outer membrane through an envelope-spanning bridge.

The outer membrane of Gram-negative bacteria provides a formidable barrier, essential for both pathogenesis and antimicrobial resistance. Biogenesis of the outer membrane requires the transport of phospholipids across the cell envelope. Recently, YhdP was implicated as a major protagonist in the transport of phospholipids from the inner membrane to the outer membrane however the molecular mechanism of YhdP mediated transport remains elusive. Here, utilising AlphaFold, we observe YhdP to form an elongated assembly of 60 ß strands that curve to form a continuous hydrophobic groove. This architecture is consistent with our negative stain electron microscopy data which reveals YhdP to be approximately 250 Å in length and thus sufficient to span the bacterial cell envelope. Furthermore, molecular dynamics simulations and in vivo bacterial growth assays indicate essential helical regions at the N- and C-termini of YhdP, that may embed into the inner and outer membranes respectively, reinforcing its envelope spanning nature. Our in vivo crosslinking data reveal phosphate-containing substrates captured along the length of the YhdP groove, providing direct evidence that YhdP transports phospholipids. This finding is congruent with our molecular dynamics simulations which demonstrate the propensity for inner membrane lipids to spontaneously enter the groove of YhdP. Collectively, our results support a model in which YhdP bridges the cell envelope, providing a hydrophobic environment for the transport of phospholipids to the outer membrane.

Computational microbiology of bacteria: Advancements in molecular dynamics simulations.

Microbiology is traditionally considered within the context of wet laboratory methodologies. Computational techniques have a great potential to contribute to microbiology. Here, we describe our loose definition of "computational microbiology" and provide a short survey focused on molecular dynamics simulations of bacterial systems that fall within this definition. It is our contention that increased compositional complexity and realistic levels of molecular crowding within simulated systems are key for bridging the divide between experimental and computational microbiology.

The solute carrier SPNS2 recruits PI(4,5)P2 to synergistically regulate transport of sphingosine-1-phosphate.

Solute carrier spinster homolog 2 (SPNS2), one of only four known major facilitator superfamily (MFS) lysolipid transporters in humans, exports sphingosine-1-phosphate (S1P) across cell membranes. Here, we explore the synergistic effects of lipid binding and conformational dynamics on SPNS2's transport mechanism. Using mass spectrometry, we discovered that SPNS2 interacts preferentially with PI(4,5)P2. Together with functional studies and molecular dynamics (MD) simulations, we identified potential PI(4,5)P2 binding sites. Mutagenesis of proposed lipid binding sites and inhibition of PI(4,5)P2 synthesis reduce S1P transport, whereas the absence of the N terminus renders the transporter essentially inactive. Probing the conformational dynamics of SPNS2, we show how synergistic binding of PI(4,5)P2 and S1P facilitates transport, increases dynamics of the extracellular gate, and stabilizes the intracellular gate. Given that SPNS2 transports a key signaling lipid, our results have implications for therapeutic targeting and also illustrate a regulatory mechanism for MFS transporters.

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.

Structural basis for bacterial energy extraction from atmospheric hydrogen.

Diverse aerobic bacteria use atmospheric H2 as an energy source for growth and survival1. This globally significant process regulates the composition of the atmosphere, enhances soil biodiversity and drives primary production in extreme environments2,3. Atmospheric H2 oxidation is attributed to uncharacterized members of the [NiFe] hydrogenase superfamily4,5. However, it remains unresolved how these enzymes overcome the extraordinary catalytic challenge of oxidizing picomolar levels of H2 amid ambient levels of the catalytic poison O2 and how the derived electrons are transferred to the respiratory chain1. Here we determined the cryo-electron microscopy structure of the Mycobacterium smegmatis hydrogenase Huc and investigated its mechanism. Huc is a highly efficient oxygen-insensitive enzyme that couples oxidation of atmospheric H2 to the hydrogenation of the respiratory electron carrier menaquinone. Huc uses narrow hydrophobic gas channels to selectively bind atmospheric H2 at the expense of O2, and 3 [3Fe-4S] clusters modulate the properties of the enzyme so that atmospheric H2 oxidation is energetically feasible. The Huc catalytic subunits form an octameric 833 kDa complex around a membrane-associated stalk, which transports and reduces menaquinone 94 Å from the membrane. These findings provide a mechanistic basis for the biogeochemically and ecologically important process of atmospheric H2 oxidation, uncover a mode of energy coupling dependent on long-range quinone transport, and pave the way for the development of catalysts that oxidize H2 in ambient air.

#COVIDisAirborne: AI-enabled multiscale computational microscopy of delta SARS-CoV-2 in a respiratory aerosol.

We seek to completely revise current models of airborne transmission of respiratory viruses by providing never-before-seen atomic-level views of the SARS-CoV-2 virus within a respiratory aerosol. Our work dramatically extends the capabilities of multiscale computational microscopy to address the significant gaps that exist in current experimental methods, which are limited in their ability to interrogate aerosols at the atomic/molecular level and thus obscure our understanding of airborne transmission. We demonstrate how our integrated data-driven platform provides a new way of exploring the composition, structure, and dynamics of aerosols and aerosolized viruses, while driving simulation method development along several important axes. We present a series of initial scientific discoveries for the SARS-CoV-2 Delta variant, noting that the full scientific impact of this work has yet to be realized.

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

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

A novel fold for acyltransferase-3 (AT3) proteins provides a framework for transmembrane acyl-group transfer.

Acylation of diverse carbohydrates occurs across all domains of life and can be catalysed by proteins with a membrane bound acyltransferase-3 (AT3) domain (PF01757). In bacteria, these proteins are essential in processes including symbiosis, resistance to viruses and antimicrobials, and biosynthesis of antibiotics, yet their structure and mechanism are largely unknown. In this study, evolutionary co-variance analysis was used to build a computational model of the structure of a bacterial O-antigen modifying acetyltransferase, OafB. The resulting structure exhibited a novel fold for the AT3 domain, which molecular dynamics simulations demonstrated is stable in the membrane. The AT3 domain contains 10 transmembrane helices arranged to form a large cytoplasmic cavity lined by residues known to be essential for function. Further molecular dynamics simulations support a model where the acyl-coA donor spans the membrane through accessing a pore created by movement of an important loop capping the inner cavity, enabling OafB to present the acetyl group close to the likely catalytic resides on the extracytoplasmic surface. Limited but important interactions with the fused SGNH domain in OafB are identified, and modelling suggests this domain is mobile and can both accept acyl-groups from the AT3 and then reach beyond the membrane to reach acceptor substrates. Together this new general model of AT3 function provides a framework for the development of inhibitors that could abrogate critical functions of bacterial pathogens.

The fatty membrane that surrounds cells is an essential feature of all living things. It is a selective barrier, only allowing certain substances to enter and exit the cell, and it contains the proteins and carbohydrates that the cell uses to interact with its environment. In bacteria, the carbohydrates on the outer side of the membrane can become 'tagged' or modified with small chemical entities which often prove useful for the cell. Acyl groups, for example, allow disease-causing bacteria to evade the immune system and contribute to infections persisting in the body. As a rule, activated acyl groups are only found inside the cell, so they need to move across the membrane before they can be attached onto the carbohydrates at the surface. This transfer is performed by a group of proteins that sit within the membrane called the acyltransferase-3 (AT3) family. The structure of these proteins and the mechanism by which they facilitate membrane crossing have remained unclear. Newman, Tindall et al. combined computational and structural modelling techniques with existing experimental data to establish how this family of proteins moves acyl groups across the membrane. They focused on OafB, an AT3 protein from the foodborne bacterial pathogen Salmonella typhimurium. The experimental data used by the team included information about which parts of OafB are necessary for this protein to acylate carbohydrates molecules. In their experiments, Newman, Tindall et al. studied how different parts of OafB move, how they interact with the molecules that carry an acyl group to the membrane, and how the acyl group is then transferred to the carbohydrate acceptor. Their results suggest that AT3 family proteins have a central pore or hole, plugged by a loop. This loop moves and therefore 'unplug' the pore, resulting in the emergence of a channel across the membrane. This channel can accommodate the acyl-donating molecule, presenting the acyl group to the outer surface of the membrane where it can be transferred to the acceptor carbohydrate. The AT3 family of proteins participates in many cellular processes involving the membrane, and a range of bacterial pathogens rely on these proteins to successfully infect human hosts. The results of Newman Tindall et al. could therefore be used across the biological sciences to provide more detailed understanding of the membrane, and to inform the design of drugs to fight bacterial diseases.

Donor-strand exchange drives assembly of the TasA scaffold in Bacillus subtilis biofilms.

Many bacteria in nature exist in multicellular communities termed biofilms, where cells are embedded in an extracellular matrix that provides rigidity to the biofilm and protects cells from chemical and mechanical stresses. In the Gram-positive model bacterium Bacillus subtilis, TasA is the major protein component of the biofilm matrix, where it has been reported to form functional amyloid fibres contributing to biofilm structure and stability. Here, we present electron cryomicroscopy structures of TasA fibres, which show that, rather than forming amyloid fibrils, TasA monomers assemble into fibres through donor-strand exchange, with each subunit donating a ß-strand to complete the fold of the next subunit along the fibre. Combining electron cryotomography, atomic force microscopy, and mutational studies, we show how TasA fibres congregate in three dimensions to form abundant fibre bundles that are essential for B. subtilis biofilm formation. Our study explains the previously observed biochemical properties of TasA and shows how a bacterial extracellular globular protein can assemble from monomers into ß-sheet-rich fibres, and how such fibres assemble into bundles in biofilms.

Lipids mediate supramolecular outer membrane protein assembly in bacteria.

ß Barrel outer membrane proteins (OMPs) cluster into supramolecular assemblies that give function to the outer membrane (OM) of Gram-negative bacteria. How such assemblies form is unknown. Here, through photoactivatable cross-linking into the Escherichia coli OM, coupled with simulations, and biochemical and biophysical analysis, we uncover the basis for OMP clustering in vivo. OMPs are typically surrounded by an annular shell of asymmetric lipids that mediate higher-order complexes with neighboring OMPs. OMP assemblies center on the abundant porins OmpF and OmpC, against which low-abundance monomeric ß barrels, such as TonB-dependent transporters, are packed. Our study reveals OMP-lipid-OMP complexes to be the basic unit of supramolecular OMP assembly that, by extending across the entire cell surface, couples the requisite multifunctionality of the OM to its stability and impermeability.

Structural basis for host recognition and superinfection exclusion by bacteriophage T5.

A key but poorly understood stage of the bacteriophage life cycle is the binding of phage receptor-binding proteins (RBPs) to receptors on the host cell surface, leading to injection of the phage genome and, for lytic phages, host cell lysis. To prevent secondary infection by the same or a closely related phage and nonproductive phage adsorption to lysed cell fragments, superinfection exclusion (SE) proteins can prevent the binding of RBPs via modulation of the host receptor structure in ways that are also unclear. Here, we present the cryogenic electron microscopy (cryo-EM) structure of the phage T5 outer membrane (OM) receptor FhuA in complex with the T5 RBP pb5, and the crystal structure of FhuA complexed to the OM SE lipoprotein Llp. Pb5 inserts four loops deeply into the extracellular lumen of FhuA and contacts the plug but does not cause any conformational changes in the receptor, supporting the view that DNA translocation does not occur through the lumen of OM channels. The FhuA-Llp structure reveals that Llp is periplasmic and binds to a nonnative conformation of the plug of FhuA, causing the inward folding of two extracellular loops via "reverse" allostery. The inward-folded loops of FhuA overlap with the pb5 binding site, explaining how Llp binding to FhuA abolishes further infection of Escherichia coli by phage T5 and suggesting a mechanism for SE via the jamming of TonB-dependent transporters by small phage lipoproteins.

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.

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

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

Corrigendum to "Atomistic level characterisation of ssDNA translocation through the E. coli proteins CsgG and CsgF for nanopore sequencing" [Comput. Struct. Biotechnol. J. 19 (2021) 6417-6430].

[This corrects the article DOI: 10.1016/j.csbj.2021.11.014.].

Lipid bilayers as potential ice nucleating agents.

Cellular damage is a key issue in the context of cryopreservation. Much of this damage is believed to be caused by extracellular ice formation at temperatures well above the homogeneous freezing point of pure water. Hence the question: what initiates ice nucleation during cryopreservation? In this paper, we assess whether cellular membranes could be responsible for facilitating the ice nucleation process, and what characteristics would make them good or bad ice nucleating agents. By means of molecular dynamics simulations, we investigate a number of phospholipids and lipopolysaccharide bilayers at the interface with supercooled liquid water. While these systems certainly appear to act as ice nucleating agents, it is likely that other impurities might also play a role in initiating extracellular ice nucleation. Furthermore, we elucidate the factors which affect a bilayer's ability to act as an ice nucleating agent; these are complex, with specific reference to both chemical and structural factors. These findings represent a first attempt to pinpoint the origin of extracellular ice nucleation, with important implications for the cryopreservation process.