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.

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.

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.

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.

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.

Manipulating turn residues on de novo designed ß-hairpin peptides for selectivity against drug-resistant bacteria.

Synthetic ß-hairpin antimicrobial peptides (AMPs) offer a useful source for the development of novel antimicrobial agents. ß-hairpin peptides generally consist of two side strands bridged by a reverse turn. In literature, most studies focused on the modifications of the side strands to manipulate the stability and activity of ß-hairpin peptides, and much less is known about the impact of the turn region. By designing a series of de novo ß-hairpin peptides with identical side strands but varied turns, we demonstrated that mutations of only 2 to 4 amino acids at the turn region could impart a wide range of antimicrobial profiles among synthetic ß-hairpin AMPs. BTT2-4 and BTT6 displayed selective potency against Gram-negative bacteria, with minimum inhibitory concentrations (MICs) of 4-8 µM. In contrast, BTT1 exhibited broad-spectrum activity, with MICs of 4-8 µM against both Gram-positive and Gram-negative strains. Additionally, BTT1 was potent against methicillin-resistant Staphylococcus aureus (MRSA) and colistin-resistant Enterobacterales. The antimicrobial potency of BTT1 persisted after 14 days of serial passage. Mechanistic studies revealed that interactions between lipopolysaccharide (LPS) and the peptides were critical to their membranolytic activity against the bacterial inner membrane. Aside from folding stability, we observed that a degree of conformational flexibility was required for disruptive membrane interactions. STATEMENT OF SIGNIFICANCE: By examining the significance of the turn region of ß-hairpin peptides, we present valuable knowledge to the design toolkit of novel antimicrobial peptides as alternative therapeutics to overcome antibiotic resistance. Our de novo designed synthetic peptides displayed selective activity against Gram-negative bacteria and potent activity against clinically relevant antibiotic-resistant strains (e.g. colistin-resistant Enterobacterales and methicillin-resistant Staphylococcus aureus). The bactericidal activity of our peptides was shown to be robust in the presence of proteolytic trypsin and saline, conditions that could suppress peptide activity. Our peptides were also determined to be non-cytotoxic against a human cell line.

The nanotube express: Delivering a stapled peptide to the cell surface.

HYPOTHESIS: Carbon nanotubes (CNTs) represent a novel platform for cellular delivery of therapeutic peptides. Chemically-functionalized CNTs may enhance peptide uptake by improving their membrane targeting properties. EXPERIMENTS: Using coarse-grained (CG) molecular dynamics (MD) simulations, we investigate membrane interactions of a peptide conjugated to pristine and chemically-modified CNTs. As proof of principle, we focus on their interactions with PM2, an amphipathic stapled peptide that inhibits the E3 ubiquitin ligase HDM2 from negatively regulating the p53 tumor suppressor. CNT interaction with both simple planar lipid bilayers as well as spherical lipid vesicles was studied, the latter as a surrogate for curved cellular membranes. FINDINGS: Membrane permeation was rapid and spontaneous for both pristine and oxidized CNTs when unconjugated. This was slowed upon addition of a noncovalently attached peptide surface "sheath", which may be an effective way to slow CNT entry and avert membrane rupture. The CNT conjugates were observed to "desheath" their peptide layer at the bilayer interface upon insertion, leaving their cargo behind in the outer leaflet. This suggests that a synergy may exist to optimize CNT safety whilst enhancing the delivery efficiency of "hitchhiking" therapeutic molecules.

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.

A Benzene-Mapping Approach for Uncovering Cryptic Pockets in Membrane-Bound Proteins.

Molecular dynamics (MD) simulations in combination with small organic probes present in the solvent have previously been used as a method to reveal cryptic pockets that may not have been identified in experimental structures. We report such a method implemented within the CHARMM force field using the GROMACS simulation package to effectively explore cryptic pockets on the surfaces of membrane-embedded proteins using benzene as a probe molecule. This method, for which we have made implementation files freely available, relies on modified nonbonded parameters in addition to repulsive potentials between membrane lipids and benzene molecules. The method was tested on part of the outer shell of the dengue virus (DENV), for which research into a safe and effective neutralizing antibody or drug molecule is still ongoing. In particular, the envelope (E) protein, associated with the membrane (M) protein, is a lipid membrane-embedded complex which forms a dimer in the mature viral envelope. Solvent mapping was performed for the full, membrane-embedded EM protein complex and compared with similar calculations performed for the isolated, soluble E protein ectodomain dimer in the solvent. Ectodomain-only simulations with benzene exhibited unfolding effects not observed in the more physiologically relevant membrane-associated systems. A cryptic pocket which has been experimentally shown to bind n-octyl-ß-d-glucoside detergent was consistently revealed in all benzene-containing simulations. The addition of benzene also enhanced the flexibility and hydrophobic exposure of cryptic pockets at a key, functional interface in the E protein and revealed a novel, potentially druggable pocket that may be targeted to prevent conformational changes associated with viral entry into the cell.

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.

Characterizing the Hydration Properties of Proton Binding Sites in the ATP Synthase c-Rings of Bacillus Species.

The membrane-embedded domain of ATP synthases contains the c-ring, which translocates ions across the membrane, and its resultant rotation is coupled to ATP synthesis in the extramembranous domain. During rotation, the c-ring becomes accessible on both sides of the lipid bilayer to solvent via channels connected to the other membrane-embedded component, the a subunit, and thereby allows the ion to be released into the solvent environment. In recent times, many experimental structures of c-rings from different species have been solved. In some of these, a water molecule with a proposed "structural role" has been identified within the c-ring ion binding site, but in general, the requirement for high resolution to resolve specific water densities complicates their interpretation. In the present study, we use molecular dynamics (MD) simulations and rigorous free energy calculations to characterize the dynamics and energetics of a water molecule within the ion binding site of the c-ring from Bacillus pseudofirmus OF4, in its wild type (WT) and P51A mutant forms, along with the c-ring from thermophilic Bacillus PS3. Our data suggest that a water molecule stably binds to the P51A mutant, as well as helping to identify a bound water molecule in Bacillus PS3 whose presence was previously overlooked due to the limited resolution of the structural data. Sequence analysis further identifies a novel conserved sequence motif that is likely required to harbor a water molecule for stable ion coordination in the binding site of such proteins.

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.

Multiscale modelling and simulation of viruses.

In recent years, advances in structural biology, integrative modelling, and simulation approaches have allowed us to gain unprecedented insights into viral structure and dynamics. In this article we survey recent studies utilizing this wealth of structural information to build computational models of partial or complete viruses and to elucidate mechanisms of viral function. Additionally, the close interplay of viral pathogens with host factors - such as cellular and intracellular membranes, receptors, antibodies, and other host proteins - makes accurate models of viral interactions and dynamics essential. As viruses continue to pose severe challenges in prevention and treatment, enhancing our mechanistic understanding of viral infection is vital to enable the development of novel therapeutic strategies.

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 structural basis for membrane assembly of immunoreceptor signalling complexes.

Immunoreceptors are TM complexes that consist of separate ligand-binding and signal-transducing modules. Mounting evidence suggests that interactions with the local environment may influence the architecture of these TM domains, which assemble via crucial sets of conserved ionisable residues, and also control the peripheral association of immunoreceptor tyrosine-based activation motifs (ITAMs) whose phosphorylation triggers cytoplasmic signalling cascades. We now report a molecular dynamics (MD) simulation study of the archetypal T cell receptor (TCR) and its cluster of differentiation 3 (CD3) signalling partners, along with the analogous DNAX-activation protein of 12 kDa (DAP12)/natural killer group 2C (NKG2C) complex. Based on > 15 µs of explicitly solvated, atomic-resolution sampling, we explore molecular aspects of immunoreceptor complex stability in different functionally relevant states. A novel alchemical approach is used to simulate the cytoplasmic CD3ε tail at different depths within lipid bilayer models, revealing that the conformation and cytoplasmic exposure of ITAMs are highly sensitive to local enrichment by different lipid species and to phosphorylation. Furthermore, simulations of the TCR and DAP12 TM domains in various states of oligomerisation suggest that, during the early stages of assembly, stable membrane insertion is facilitated by the interfacial lipid/solvent environment and/or partial ionisation of charged residues. Collectively, our results indicate that the architecture and mechanisms of signal transduction in immunoreceptor complexes are tightly regulated by interactions with the microenvironment.