Search

β-Lactamases, which hydrolyse β-lactam antibiotics, are key determinants of antibiotic resistance. Predicting the sites and effects of distal mutations in enzymes is challenging. For β-lactamases, the ability to make such predictions would contribute to understanding activity against, and development of, antibiotics and inhibitors to combat resistance. Here, using dynamical non-equilibrium molecular dynamics (D-NEMD) simulations combined with experiments, we demonstrate that intramolecular communication networks differ in three class A SulpHydryl Variant (SHV)-type β-lactamases. Differences in network architecture and correlated motions link to catalytic efficiency and β-lactam substrate spectrum. Further, the simulations identify a distal residue at position 89 in the clinically important Klebsiella pneumoniae carbapenemase 2 (KPC-2), as a participant in similar networks, suggesting that mutation at this position would modulate enzyme activity. Experimental kinetic, biophysical and structural characterisation of the naturally occurring, but previously biochemically uncharacterised, KPC-2 G89D mutant with several antibiotics and inhibitors reveals significant changes in hydrolytic spectrum, specifically reducing activity towards carbapenems without effecting major structural or stability changes. These results show that D-NEMD simulations can predict distal sites where mutation affects enzyme activity. This approach could have broad application in understanding enzyme evolution, and in engineering of natural and de novo enzymes., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Mycobacterium tuberculosis is the single most important global infectious disease killer and a World Health Organization critical priority pathogen for development of new antimicrobials. M. tuberculosis DNA gyrase is a validated target for anti-TB agents, but those in current use target DNA breakage-reunion, rather than the ATPase activity of the GyrB subunit. Here, virtual screening, subsequently validated by whole-cell and enzyme inhibition assays, was applied to identify candidate compounds that inhibit M. tuberculosis GyrB ATPase activity from the Specs compound library. This approach yielded six compounds: four carbazole derivatives ( 1 , 2 , 3 , and 8 ), the benzoindole derivative 11 , and the indole derivative 14 . Carbazole derivatives can be considered a new scaffold for M. tuberculosis DNA gyrase ATPase inhibitors. IC 50 values of compounds 8 , 11 , and 14 (0.26, 0.56, and 0.08 μM, respectively) for inhibition of M. tuberculosis DNA gyrase ATPase activity are 5-fold, 2-fold, and 16-fold better than the known DNA gyrase ATPase inhibitor novobiocin. MIC values of these compounds against growth of M. tuberculosis H37Ra are 25.0, 3.1, and 6.2 μg/mL, respectively, superior to novobiocin (MIC > 100.0 μg/mL). Molecular dynamics simulations of models of docked GyrB:inhibitor complexes suggest that hydrogen bond interactions with GyrB Asp79 are crucial for high-affinity binding of compounds 8 , 11 , and 14 to M. tuberculosis GyrB for inhibition of ATPase activity. These data demonstrate that virtual screening can identify known and new scaffolds that inhibit both M. tuberculosis DNA gyrase ATPase activity in vitro and growth of M. tuberculosis bacteria.

Methyl carboxylate esters promote the formation of dimethyl ether (DME) from the dehydration of methanol in H-ZSM-5 zeolite. We employ a multilevel quantum method to explore the possible associative and dissociative mechanisms in the presence, and absence, of six methyl ester promoters. This hybrid method combines density functional theory, with dispersion corrections (DFT-D3), for the full periodic system, with second-order Møller-Plesset perturbation theory (MP2) for small clusters representing the reaction site, and coupled cluster with single, double, and perturbative triple substitution (CCSD(T)) for the reacting molecules. The calculated adsorption enthalpy of methanol, and reaction enthalpies of the dehydration of methanol to DME within H-ZSM-5, agree with experiment to within chemical accuracy (∼4 kJ mol -1 ). For the promoters, a reaction pathway via an associative mechanism gives lower overall reaction enthalpies and barriers compared to the reaction with methanol only. Each stage of this mechanism is explored and related to experimental data. We provide evidence that suggests the promoter's adsorption to the Brønsted acid site is the most important factor dictating its efficiency.

DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis , which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A 2 B 2 ) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.

Sire is a Python/C++ library that is used both to prototype new algorithms and as an interoperability engine for exchanging information between molecular simulation programs. It provides a collection of file parsers and information converters that together make it easier to combine and leverage the functionality of many other programs and libraries. This empowers researchers to use sire to write a single script that can, for example, load a molecule from a PDBx/mmCIF file via Gemmi, perform SMARTS searches via RDKit, parameterize molecules using BioSimSpace, run GPU-accelerated molecular dynamics via OpenMM, and then display the resulting dynamics trajectory in a NGLView Jupyter notebook 3D molecular viewer. This functionality is built on by BioSimSpace, which uses sire's molecular information engine to interconvert with programs such as GROMACS, NAMD, Amber, and AmberTools for automated molecular parameterization and the running of molecular dynamics, metadynamics, and alchemical free energy workflows. Sire comes complete with a powerful molecular information search engine, plus trajectory loading and editing, analysis, and energy evaluation engines. This, when combined with an in-built computer algebra system, gives substantial flexibility to researchers to load, search for, edit, and combine molecular information from multiple sources and use that to drive novel algorithms by combining functionality from other programs. Sire is open source (GPL3) and is available via conda and at a free Jupyter notebook server at https://try.openbiosim.org. Sire is supported by the not-for-profit OpenBioSim community interest company., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Resistance to antibiotics is a public health crisis. Although carbapenems are less susceptible to resistance than other β-lactam antibiotics, β-lactamases mediating resistance against these drugs are spreading. Here, we dissect the contributions of electric fields to carbapenemase activity in class A β-lactamases. We perform QM/MM molecular dynamics simulations of meropenem acyl-enzyme hydrolysis that correctly discriminate carbapenemases. Electric field analysis shows that active-site fields in the deacylation transition state and tetrahedral intermediate are important determinants of activity. The active-site fields identify several residues, some distal, that distinguish efficient carbapenemases. Our field analysis script (www.github.com/bunzela/FieldTools) may help in understanding and combating antibiotic resistance., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Introduction: The pathogenic bacterium Helicobacter pylori has evolved glycan-mediated mechanisms to evade host immune defenses. This study tests the hypothesis that genetic disruption of H. pylori glycan biosynthesis alters immune recognition and response by human gastric epithelial cells and monocyte-derived dendritic cells., Methods: To test this hypothesis, human cell lines were challenged with wildtype H. pylori alongside an array of H. pylori glycosylation mutants. The relative levels of immune response were measured via immature dendritic cell maturation and cytokine secretion., Results: Our findings indicate that disruption of lipopolysaccharide biosynthesis diminishes gastric cytokine production, without disrupting dendritic cell recognition and activation. In contrast, variable immune responses were observed in protein glycosylation mutants which prompted us to test the hypothesis that phase variation plays a role in regulating bacterial cell surface glycosylation and subsequent immune recognition. Lewis antigen presentation does not correlate with extent of immune response, while the extent of lipopolysaccharide O-antigen elaboration does., Discussion: The outcomes of this study demonstrate that H. pylori glycans modulate the host immune response. This work provides a foundation to pursue immune-based tailoring of bacterial glycans towards modulating immunogenicity of microbial pathogens., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Barrett, Kassama, Surks, Mulholland, Moulton and Dube.)

Many enzymes display non-Arrhenius behavior with curved Arrhenius plots in the absence of denaturation. There has been significant debate about the origin of this behavior and recently the role of the activation heat capacity (Δ C P ⧧ ) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of Δ C P ⧧ (arising from narrowing of the conformational space along the reaction coordinate), then curved Arrhenius plots are a consequence. To investigate these phenomena in detail, we have collected high precision temperature-rate data over a wide temperature interval for a model glycosidase enzyme MalL, and a series of mutants that change the temperature-dependence of the enzyme-catalyzed rate. We use these data to test a range of models including macromolecular rate theory (MMRT) and an equilibrium model. In addition, we have performed extensive molecular dynamics (MD) simulations to characterize the conformational landscape traversed by MalL in the enzyme-substrate complex and an enzyme-transition state complex. We have crystallized the enzyme in a transition state-like conformation in the absence of a ligand and determined an X-ray crystal structure at very high resolution (1.10 Å). We show (using simulation) that this enzyme-transition state conformation has a more restricted conformational landscape than the wildtype enzyme. We coin the term "transition state-like conformation (TLC)" to apply to this state of the enzyme. Together, these results imply a cooperative conformational transition between an enzyme-substrate conformation (ES) and a transition-state-like conformation (TLC) that precedes the chemical step. We present a two-state model as an extension of MMRT (MMRT-2S) that describes the data along with a convenient approximation with linear temperature dependence of the activation heat capacity (MMRT-1L) that can be used where fewer data points are available. Our model rationalizes disparate behavior seen for MalL and previous results for a thermophilic alcohol dehydrogenase and is consistent with a raft of data for other enzymes. Our model can be used to characterize the conformational changes required for enzyme catalysis and provides insights into the role of cooperative conformational changes in transition state stabilization that are accompanied by changes in heat capacity for the system along the reaction coordinate. TLCs are likely to be of wide importance in understanding the temperature dependence of enzyme activity and other aspects of enzyme catalysis., Competing Interests: The authors declare no competing financial interest., (© 2024 American Chemical Society.)

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise ( e.g. , mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g. , detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

The tunable design of protein redox potentials promises to open a range of applications in biotechnology and catalysis. Here, we introduce a method to calculate redox potential changes by combining fluctuation relations with molecular dynamics simulations. It involves the simulation of reduced and oxidized states, followed by the instantaneous conversion between them. Energy differences introduced by the perturbations are obtained using the Kubo-Onsager approach. Using a detailed fluctuation relation coupled with Bayesian inference, these are postprocessed into estimates for the redox potentials in an efficient manner. This new method, denoted MD + CB, is tested on a de novo four-helix bundle heme protein (the m4D2 "maquette") and five designed mutants, including some mutants characterized experimentally in this work. The MD + CB approach is found to perform reliably, giving redox potential shifts with reasonably good correlation (0.85) to the experimental values for the mutants. The MD + CB approach also compares well with redox potential shift predictions using a continuum electrostatic method. The estimation method employed within the MD + CB approach is straightforwardly transferable to standard equilibrium MD simulations and holds promise for redox protein engineering and design applications.

Epistasis, the non-additive effect of mutations, can provide combinatorial improvements to enzyme activity that substantially exceed the gains from individual mutations. Yet the molecular mechanisms of epistasis remain elusive, undermining our ability to predict pathogen evolution and engineer biocatalysts. Here we reveal how directed evolution of a β-lactamase yielded highly epistatic activity enhancements. Evolution selected four mutations that increase antibiotic resistance 40-fold, despite their marginal individual effects (≤2-fold). Synergistic improvements coincided with the introduction of super-stochiometric burst kinetics, indicating that epistasis is rooted in the enzyme's conformational dynamics. Our analysis reveals that epistasis stemmed from distinct effects of each mutation on the catalytic cycle. The initial mutation increased protein flexibility and accelerated substrate binding, which is rate-limiting in the wild-type enzyme. Subsequent mutations predominantly boosted the chemical steps by fine-tuning substrate interactions. Our work identifies an overlooked cause for epistasis: changing the rate-limiting step can result in substantial synergy that boosts enzyme activity., Competing Interests: Competing interestsThe authors declare no competing interests., (© The Author(s) 2024, corrected publication 2024.)

Interactive molecular dynamics simulation in virtual reality (iMD-VR) is emerging as a promising technique in molecular science. Here, we demonstrate its use in a range of fifteen applications in materials science and heterogeneous catalysis. In this work, the iMD-VR package Narupa is used with the MD package, DL_POLY [1]. We show how iMD-VR can be used to: (i) investigate the mechanism of lithium fast ion conduction by directing the formation of defects showing that vacancy transport is favoured over interstitialcy mechanisms, and (ii) guide a molecule through a zeolite pore to explore diffusion within zeolites, examining in detail the motion of methyl n-hexanoate in H-ZSM-5 zeolite and identifying bottlenecks restricting diffusion. iMD-VR allows users to manipulate these systems intuitively, to drive changes in them and observe the resulting changes in structure and dynamics. We make these simulations available, as a resource for both teaching and research. All simulation files, with videos, can be found online (https://doi.org/10.5281/zenodo.8252314) and are provided as open-source material., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:Joe Crossley-Lewis reports financial support was provided by BP Plc. Josh Dunn reports financial support was provided by BP Plc., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Due to COVID-19, pandemic preparedness emerges as a key imperative, necessitating new approaches to accelerate development of reagents against infectious pathogens., Methods: Here, we developed an integrated approach combining synthetic, computational and structural methods with in vitro antibody selection and in vivo immunization to design, produce and validate nature-inspired nanoparticle-based reagents against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)., Results: Our approach resulted in two innovations: (i) a thermostable nasal vaccine called ADDoCoV, displaying multiple copies of a SARS-CoV-2 receptor binding motif derived epitope and (ii) a multivalent nanoparticle superbinder, called Gigabody, against SARS-CoV-2 including immune-evasive variants of concern (VOCs). In vitro generated neutralizing nanobodies and electron cryo-microscopy established authenticity and accessibility of epitopes displayed by ADDoCoV. Gigabody comprising multimerized nanobodies prevented SARS-CoV-2 virion attachment with picomolar EC 50 . Vaccinating mice resulted in antibodies cross-reacting with VOCs including Delta and Omicron., Conclusion: Our study elucidates Adenovirus-derived dodecamer (ADDomer)-based nanoparticles for use in active and passive immunization and provides a blueprint for crafting reagents to combat respiratory viral infections., Competing Interests: C.S., K.G. and I.B. report shareholding in Halo Therapeutics Ltd unrelated to this Correspondence. I.B. reports shareholding in Geneva Biotech SARL, unrelated to this correspondence. F.G., J.H. and I.B. report shareholding in Imophoron Ltd, related to this Correspondence. Patents and patent applications have been filed related to ADDomer vaccines and therapeutics (WO2017167988A, EP22191583.8). The other authors do not declare competing interests. ADDomer is a registered trademark of Imophoron Ltd., (© The Author(s) 2023. Published by Oxford University Press on behalf of Antibody Therapeutics. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight., (© 2023. Springer Nature Limited.)

Owing to the emergence of antibiotic resistance, the polymyxin colistin has been recently revived to treat acute, multidrug-resistant Gram-negative bacterial infections. Positively charged colistin binds to negatively charged lipids and damages the outer membrane of Gram-negative bacteria. However, the MCR-1 protein, encoded by the mobile colistin resistance (mcr) gene, is involved in bacterial colistin resistance by catalysing phosphoethanolamine (PEA) transfer onto lipid A, neutralising its negative charge, and thereby reducing its interaction with colistin. Our preliminary results showed that treatment with a reference pyrazolone compound significantly reduced colistin minimal inhibitory concentrations in Escherichia coli expressing mcr-1 mediated colistin resistance (Hanpaibool et al. in ACS Omega, 2023). A docking-MD combination was used in an ensemble-based docking approach to identify further pyrazolone compounds as candidate MCR-1 inhibitors. Docking simulations revealed that 13/28 of the pyrazolone compounds tested are predicted to have lower binding free energies than the reference compound. Four of these were chosen for in vitro testing, with the results demonstrating that all the compounds tested could lower colistin MICs in an E. coli strain carrying the mcr-1 gene. Docking of pyrazolones into the MCR-1 active site reveals residues that are implicated in ligand-protein interactions, particularly E246, T285, H395, H466, and H478, which are located in the MCR-1 active site and which participate in interactions with MCR-1 in ≥ 8/10 of the lowest energy complexes. This study establishes pyrazolone-induced colistin susceptibility in E. coli carrying the mcr-1 gene, providing a method for the development of novel treatments against colistin-resistant bacteria., (© 2023. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted., (© The Author(s) (2023). Published by Oxford University Press on behalf of Journal of Molecular Cell Biology, CEMCS, CAS.)

The target for humoral immunity, SARS-CoV-2 spike glycoprotein, has become the focus of vaccine research and development. Previous work demonstrated that the N-terminal domain (NTD) of SARS-CoV-2 spike binds biliverdin-a product of heme catabolism-causing a strong allosteric effect on the activity of a subset of neutralizing antibodies. Herein, we show that the spike glycoprotein is also able to bind heme (K D  = 0.5 ± 0.2 μM). Molecular modeling indicated that the heme group fits well within the same pocket on the SARS-CoV-2 spike NTD. Lined by aromatic and hydrophobic residues (W104, V126, I129, F192, F194, I203, and L226), the pocket provides a suitable environment to stabilize the hydrophobic heme. Mutagenesis of N121 has a substantive effect on heme binding (K D  = 3000 ± 220 μM), confirming the pocket as a major heme binding location of the viral glycoprotein. Coupled oxidation experiments in the presence of ascorbate indicated that the SARS-CoV-2 glycoprotein can catalyze the slow conversion of heme to biliverdin. The heme trapping and oxidation activities of the spike may allow the virus to reduce levels of free heme during infection to facilitate evasion of the adaptive and innate immunity., Competing Interests: Conflicts of interest The authors declare that they have no conflicts of interest with the contents of this article., (Crown Copyright © 2023. Published by Elsevier Inc. All rights reserved.)

The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.

The SARS-CoV-2 main protease (M pro ) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. M pro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how M pro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of M pro to the presence and absence of a substrate. The results highlight communication between the M pro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of M pro . More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

Is scientific partnership between the UK and EU on the Horizon?, (© 2023 The Authors.)

Mutations in DNA gyrase confer resistance to fluoroquinolones, second-line antibiotics for Mycobacterium tuberculosis infections. Identification of new agents that inhibit M. tuberculosis DNA gyrase ATPase activity is one strategy to overcome this. Here, bioisosteric designs using known inhibitors as templates were employed to define novel inhibitors of M. tuberculosis DNA gyrase ATPase activity. This yielded the modified compound R3-13 with improved drug-likeness compared to the template inhibitor that acted as a promising ATPase inhibitor against M. tuberculosis DNA gyrase. Utilization of compound R3-13 as a virtual screening template, supported by subsequent biological assays, identified seven further M. tuberculosis DNA gyrase ATPase inhibitors with IC 50 values in the range of 0.42-3.59 μM. The most active compound 1 showed an IC 50 value of 0.42 μM, 3-fold better than the comparator ATPase inhibitor novobiocin (1.27 μM). Compound 1 showed noncytotoxicity to Caco-2 cells at concentrations up to 76-fold higher than its IC 50 value. Molecular dynamics simulations followed by decomposition energy calculations identified that compound 1 occupies the binding pocket utilized by the adenosine group of the ATP analogue AMPPNP in the M. tuberculosis DNA gyrase GyrB subunit. The most prominent contribution to the binding of compound 1 to M. tuberculosis GyrB subunit is made by residue Asp79, which forms two hydrogen bonds with the OH group of this compound and also participates in the binding of AMPPNP. Compound 1 represents a potential new scaffold for further exploration and optimization as a M. tuberculosis DNA gyrase ATPase inhibitor and candidate anti-tuberculosis agent.

KPC-2 ( Klebsiella pneumoniae carbapenemase-2) is a globally disseminated serine-β-lactamase (SBL) responsible for extensive β-lactam antibiotic resistance in Gram-negative pathogens. SBLs inactivate β-lactams via a mechanism involving a hydrolytically labile covalent acyl-enzyme intermediate. Carbapenems, the most potent β-lactams, evade the activity of many SBLs by forming long-lived inhibitory acyl-enzymes; however, carbapenemases such as KPC-2 efficiently deacylate carbapenem acyl-enzymes. We present high-resolution (1.25-1.4 Å) crystal structures of KPC-2 acyl-enzymes with representative penicillins (ampicillin), cephalosporins (cefalothin), and carbapenems (imipenem, meropenem, and ertapenem) obtained utilizing an isosteric deacylation-deficient mutant (E166Q). The mobility of the Ω-loop (residues 165-170) negatively correlates with antibiotic turnover rates ( k cat ), highlighting the role of this region in positioning catalytic residues for efficient hydrolysis of different β-lactams. Carbapenem-derived acyl-enzyme structures reveal the predominance of the Δ1-(2 R ) imine rather than the Δ2 enamine tautomer. Quantum mechanics/molecular mechanics molecular dynamics simulations of KPC-2:meropenem acyl-enzyme deacylation used an adaptive string method to differentiate the reactivity of the two isomers. These identify the Δ1-(2 R ) isomer as having a significantly (7 kcal/mol) higher barrier than the Δ2 tautomer for the (rate-determining) formation of the tetrahedral deacylation intermediate. Deacylation is therefore likely to proceed predominantly from the Δ2, rather than the Δ1-(2 R ) acyl-enzyme, facilitated by tautomer-specific differences in hydrogen-bonding networks involving the carbapenem C-3 carboxylate and the deacylating water and stabilization by protonated N-4, accumulating a negative charge on the Δ2 enamine-derived oxyanion. Taken together, our data show how the flexible Ω-loop helps confer broad-spectrum activity upon KPC-2, while carbapenemase activity stems from efficient deacylation of the Δ2-enamine acyl-enzyme tautomer.

The polymyxin colistin is a last line antibiotic for extensively resistant Gram-negative bacteria. Colistin binding to lipid A disrupts the Gram-negative outer membrane, but mobile colistin resistance ( mcr ) gene family members confer resistance by catalyzing phosphoethanolamine (PEA) transfer onto lipid A, neutralizing its negative charge to reduce colistin interactions. Multiple mcr isoforms have been identified in clinical and environmental isolates, with mcr-1 being the most widespread and mcr-3 being common in South and East Asia. Preliminary screening revealed that treatment with pyrazolones significantly reduced mcr-1 , but not mcr-3 , mediated colistin resistance. Molecular dynamics (MD) simulations of the catalytic domains of MCR-1 and a homology model of MCR-3, in different protonation states of active site residues H395/H380 and H478/H463, indicate that the MCR-1 active site has greater water accessibility than MCR-3, but that this is less influenced by changes in protonation. MD-optimized structures of MCR-1 and MCR-3 were used in virtual screening of 20 pyrazolone derivatives. Docking of these into the MCR-1/MCR-3 active sites identifies common residues likely to be involved in protein-ligand interactions, specifically the catalytic threonine (MCR-1 T285, MCR-3 T277) site of PEA addition, as well as differential interactions with adjacent amino acids. Minimal inhibitory concentration assays showed that the pyrazolone with the lowest predicted binding energy (ST3f) restores colistin susceptibility of mcr-1 , but not mcr-3 , expressing Escherichia coli . Thus, simulations indicate differences in the active site structure between MCR-1 and MCR-3 that may give rise to differences in pyrazolone binding and so relate to differential effects upon producer E. coli . This work identifies pyrazolones as able to restore colistin susceptibility of mcr-1 -producing bacteria, laying the foundation for further investigations of their activity as phosphoethanolamine transferase inhibitors as well as of their differential activity toward mcr isoforms., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

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., Competing Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2022.)

Mycobacterium tuberculosis protein kinase B (PknB) is essential to mycobacterial growth and has received considerable attention as an attractive target for novel anti-tuberculosis drug development. Here, virtual screening, validated by biological assays, was applied to select candidate inhibitors of M. tuberculosis PknB from the Specs compound library (www.specs.net). Fifteen compounds were identified as hits and selected for in vitro biological assays, of which three indoles ( 2 , AE-848/42799159; 4 , AH-262/34335013; 10 , AP-124/40904362) inhibited growth of M. tuberculosis H37Rv with minimal inhibitory concentrations of 6.2, 12.5, and 6.2 μg/mL, respectively. Two compounds, 2 and 10 , inhibited M. tuberculosis PknB activity in vitro , with IC 50 values of 14.4 and 12.1 μM, respectively, suggesting this to be the likely basis of their anti-tubercular activity. In contrast, compound 4 displayed anti-tuberculosis activity against M. tuberculosis H37Rv but showed no inhibition of PknB activity (IC 50 > 128 μM). We hypothesize that hydrolysis of its ethyl ester to a carboxylate moiety generates an active species that inhibits other M. tuberculosis enzymes. Molecular dynamics simulations of modeled complexes of compounds 2 , 4 , and 10 bound to M. tuberculosis PknB indicated that compound 4 has a lower affinity for M. tuberculosis PknB than compounds 2 and 10 , as evidenced by higher calculated binding free energies, consistent with experiment. Compounds 2 and 10 therefore represent candidate inhibitors of M. tuberculosis PknB that provide attractive starting templates for optimization as anti-tubercular agents.

The α7 nicotinic acetylcholine receptor (nAChR) is present in neuronal and non-neuronal cells and has anti-inflammatory actions. Molecular dynamics simulations suggested that α7 nAChR interacts with a region of the SARS-CoV-2 spike protein (S), and a potential contribution of nAChRs to COVID-19 pathophysiology has been proposed. We applied whole-cell and single-channel recordings to determine whether a peptide corresponding to the Y674-R685 region of the S protein can directly affect α7 nAChR function. The S fragment exerts a dual effect on α7. It activates α7 nAChRs in the presence of positive allosteric modulators, in line with our previous molecular dynamics simulations showing favourable binding of this accessible region of the S protein to the nAChR agonist binding site. The S fragment also exerts a negative modulation of α7, which is evidenced by a profound concentration-dependent decrease in the durations of openings and activation episodes of potentiated channels and in the amplitude of macroscopic responses elicited by ACh. Our study identifies a potential functional interaction between α7 nAChR and a region of the S protein, thus providing molecular foundations for further exploring the involvement of nAChRs in COVID-19 pathophysiology., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Serine/threonine protein kinase B (PknB) is essential to Mycobacterium tuberculosis (M. tuberculosis) cell division and metabolism and a potential anti-tuberculosis drug target. Here we apply Hologram Quantitative Structure Activity Relationship (HQSAR) and three-dimensional QSAR (Comparative Molecular Similarity Indices Analysis (CoMSIA)) methods to investigate structural requirements for PknB inhibition by a series of previously described quinazoline derivatives. PknB binding of quinazolines was evaluated by molecular dynamics (MD) simulations of the catalytic domain and binding energies calculated by Molecular Mechanics/Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM-GBSA) methods. Evaluation of a training set against experimental data showed both HQSAR and CoMSIA models to reliably predict quinazoline binding to PknB, and identified the quinazoline core and overall hydrophobicity as the major contributors to affinity. Calculated binding energies also agreed with experiment, and MD simulations identified hydrogen bonds to Glu93 and Val95, and hydrophobic interactions with Gly18, Phe19, Gly20, Val25, Thr99 and Met155, as crucial to PknB binding. Based on these results, additional quinazolines were designed and evaluated in silico, with HQSAR and CoMSIA models identifying sixteen compounds, with predicted PknB binding superior to the template, whose activity spectra and physicochemical, pharmacokinetic, and anti-M. tuberculosis properties were assessed. Compound, D060, bearing additional ortho- and meta-methyl groups on its R 2 substituent, was superior to template regarding PknB inhibition and % caseum fraction unbound, and equivalent in other aspects, although predictions identified hepatotoxicity as a likely issue with the quinazoline series. These data provide a structural basis for rational design of quinazoline derivatives with more potent PknB inhibitory activity as candidate anti-tuberculosis agents., (Copyright © 2022 Elsevier Inc. All rights reserved.)

β-lactam antibiotic resistance in Gram-negative bacteria, primarily caused by β-lactamase enzymes that hydrolyze the β-lactam ring, has become a serious clinical problem. Carbapenems were formerly considered "last resort" antibiotics because they escaped breakdown by most β-lactamases, due to slow deacylation of the acyl-enzyme intermediate. However, an increasing number of Gram-negative bacteria now produce β-lactamases with carbapenemase activity: these efficiently hydrolyze the carbapenem β-lactam ring, severely limiting the treatment of some bacterial infections. Here, we use quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reactions of acyl-enzyme complexes of eight β-lactamases of class A (the most widely distributed β-lactamase group) with the carbapenem meropenem to investigate differences between those inhibited by carbapenems (TEM-1, SHV-1, BlaC, and CTX-M-16) and those that hydrolyze them (SFC-1, KPC-2, NMC-A, and SME-1). QM/MM molecular dynamics simulations confirm the two enzyme groups to differ in the preferred acyl-enzyme orientation: carbapenem-inhibited enzymes favor hydrogen bonding of the carbapenem hydroxyethyl group to deacylating water (DW). QM/MM simulations of deacylation give activation free energies in good agreement with experimental hydrolysis rates, correctly distinguishing carbapenemases. For the carbapenem-inhibited enzymes, free energies for deacylation are significantly higher than for the carbapenemases, even when the hydroxyethyl group was restrained to prevent interaction with the DW. Analysis of these simulations, and additional simulations of mutant enzymes, shows how factors including the hydroxyethyl orientation, the active site volume, and architecture (conformations of Asn170 and Asn132; organization of the oxyanion hole; and the Cys69-Cys238 disulfide bond) collectively determine catalytic efficiency toward carbapenems.

We test a range of standard generalized Born (GB) models and protein force fields for a set of five experimentally characterized, designed peptides comprising alternating blocks of glutamate and lysine, which have been shown to differ significantly in α-helical content. Sixty-five combinations of force fields and GB models are evaluated in >800 μs of molecular dynamics simulations. GB models generally do not reproduce the experimentally observed α-helical content, and none perform well for all five peptides. These results illustrate that these models are not usefully predictive in this context. These peptides provide a useful test set for simulation methods.

Introduction: The potential of virtual reality (VR) to contribute to drug design and development has been recognized for many years. A recent advance is to use VR not only to visualize and interact with molecules, but also to interact with molecular dynamics simulations 'on the fly' (interactive molecular dynamics in VR, IMD-VR), which is useful for flexible docking and examining binding processes and conformational changes., Areas Covered: The authors use the term 'interactive VR' to refer to software where interactivity is an inherent part of the user VR experience e.g . in making structural modifications or interacting with a physically rigorous molecular dynamics (MD) simulation, as opposed to using VR controllers to rotate and translate the molecule for enhanced visualization. Here, they describe these methods and their application to problems relevant to drug discovery, highlighting the possibilities that they offer in this arena., Expert Opinion: The ease of viewing and manipulating molecular structures and dynamics, using accessible VR hardware, and the ability to modify structures on the fly ( e.g . adding or deleting atoms) - and for groups of researchers to work together in the same virtual environment - makes modern interactive VR a valuable tool to add to the armory of drug design and development methods.

Aim: In silico screening approaches were performed to discover novel InhA inhibitors. Methods: Candidate InhA inhibitors were obtained from the combination of virtual screening and pharmacokinetic prediction. In addition, molecular mechanics Poisson-Boltzmann surface area, molecular mechanics Generalized Born surface area and WaterSwap methods were performed to investigate the binding interactions and binding energy of candidate compounds. Results: Four candidate compounds with suitable physicochemical, pharmacokinetic and antibacterial properties are proposed. The crucial interactions of the candidate compounds were H-bond, pi-pi and sigma-pi interactions observed in the InhA binding site. The binding affinity of these compounds was improved by hydrophobic interactions with hydrophobic side chains in the InhA pocket. Conclusion: The four newly identified InhA inhibitors reported in this study could serve as promising hit compounds against Mycobacterium tuberculosis and may be considered for further experimental studies.

OXA-48 β-lactamases are frequently encountered in bacterial infections caused by carbapenem-resistant Gram-negative bacteria. Due to the importance of carbapenems in the treatment of healthcare-associated infections and the increasingly wide dissemination of OXA-48-like enzymes on plasmids, these β-lactamases are of high clinical significance. Notably, OXA-48 hydrolyzes imipenem more efficiently than other commonly used carbapenems, such as meropenem. Here, we use extensive multiscale simulations of imipenem and meropenem hydrolysis by OXA-48 to dissect the dynamics and to explore differences in the reactivity of the possible conformational substates of the respective acylenzymes. Quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reaction for both substrates demonstrate that deacylation is favored when the 6α-hydroxyethyl group is able to hydrogen bond to the water molecule responsible for deacylation but disfavored by the increasing hydration of either oxygen of the carboxylated Lys73 general base. Differences in free energy barriers calculated from the QM/MM simulations correlate well with the experimentally observed differences in hydrolytic efficiency between meropenem and imipenem. We conclude that the impaired breakdown of meropenem, compared to imipenem, which arises from a subtle change in the hydrogen bonding pattern between the deacylating water molecule and the antibiotic, is most likely induced by the meropenem 1β-methyl group. In addition to increased insights into carbapenem breakdown by OXA β-lactamases, which may aid in future efforts to design antibiotics or inhibitors, our approach exemplifies the combined use of atomistic simulations in determining the possible different enzyme-substrate substates and their influence on enzyme reaction kinetics., Competing Interests: The authors declare no competing financial interest., (© 2022 American Chemical Society.)

Mycobacterium tuberculosis DNA gyrase manipulates the DNA topology using controlled breakage and religation of DNA driven by ATP hydrolysis. DNA gyrase has been validated as the enzyme target of fluoroquinolones (FQs), second-line antibiotics used for the treatment of multidrug-resistant tuberculosis. Mutations around the DNA gyrase DNA-binding site result in the emergence of FQ resistance in M. tuberculosis ; inhibition of DNA gyrase ATPase activity is one strategy to overcome this. Here, virtual screening, subsequently validated by biological assays, was applied to select candidate inhibitors of the M. tuberculosis DNA gyrase ATPase activity from the Specs compound library (www.specs.net). Thirty compounds were identified and selected as hits for in vitro biological assays, of which two compounds, G24 and G26 , inhibited the growth of M. tuberculosis H37Rv with a minimal inhibitory concentration of 12.5 μg/mL. The two compounds inhibited DNA gyrase ATPase activity with IC 50 values of 2.69 and 2.46 μM, respectively, suggesting this to be the likely basis of their antitubercular activity. Models of complexes of compounds G24 and G26 bound to the M. tuberculosis DNA gyrase ATP-binding site, generated by molecular dynamics simulations followed by pharmacophore mapping analysis, showed hydrophobic interactions of inhibitor hydrophobic headgroups and electrostatic and hydrogen bond interactions of the polar tails, which are likely to be important for their inhibition. Decreasing compound lipophilicity by increasing the polarity of these tails then presents a likely route to improving the solubility and activity. Thus, compounds G24 and G26 provide attractive starting templates for the optimization of antitubercular agents that act by targeting DNA gyrase.

Antibiotic resistance is a major threat to global public health. β-lactamases, which catalyze breakdown of β-lactam antibiotics, are a principal cause. Metallo β-lactamases (MBLs) represent a particular challenge because they hydrolyze almost all β-lactams and to date no MBL inhibitor has been approved for clinical use. Molecular simulations can aid drug discovery, for example, predicting inhibitor complexes, but empirical molecular mechanics (MM) methods often perform poorly for metalloproteins. Here we present a multiscale approach to model thiol inhibitor binding to IMP-1, a clinically important MBL containing two catalytic zinc ions, and predict the binding mode of a 2-mercaptomethyl thiazolidine (MMTZ) inhibitor. Inhibitors were first docked into the IMP-1 active site, testing different docking programs and scoring functions on multiple crystal structures. Complexes were then subjected to molecular dynamics (MD) simulations and subsequently refined through QM/MM optimization with a density functional theory (DFT) method, B3LYP/6-31G(d), increasing the accuracy of the method with successive steps. This workflow was tested on two IMP-1:MMTZ complexes, for which it reproduced crystallographically observed binding, and applied to predict the binding mode of a third MMTZ inhibitor for which a complex structure was crystallographically intractable. We also tested a 12-6-4 nonbonded interaction model in MD simulations and optimization with a SCC-DFTB QM/MM approach. The results show the limitations of empirical models for treating these systems and indicate the need for higher level calculations, for example, DFT/MM, for reliable structural predictions. This study demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems., (© 2021 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Bacterial luciferase (Lux) catalyzes oxidation of reduced flavin mononucleotide (FMN) and aldehyde to form oxidized FMN and carboxylic acid via molecular oxygen with concomitant light generation. The enzyme is useful for various detection applications in biomedical experiments. Upon reacting with oxygen, the reduced FMN generates C4a-peroxy-FMN (FMNH-C4a-OO - ) as a reactive intermediate, which is required for light generation. However, the mechanism and control of FMNH-C4a-OO - formation are not clear. This work investigated the reaction of FMNH-C4a-OO - formation in Lux using QM/MM methods. The B3LYP/6-31G*/CHARMM27 calculations indicate that Lux controls the formation of FMNH-C4a-OO - via the conserved His44 residue. The steps in intermediate formation are found to be as follows: (i) H + reacts with O 2 to generate + OOH. (ii) + OOH attacks C4a of FMNH - to generate FMNH-C4a-OOH. (iii) H + is transferred from FMNH-C4a-OOH to His44 to generate FMNH-C4a-OO - while His44 stabilizes FMNH-C4a-OO - by forming a hydrogen bond to an oxygen atom. This controlling key mechanism for driving the change from FMNH-C4a-OOH to the FMNH-C4a-OO - adduct is confirmed because FMNH-C4a-OO - is more stable than FMNH-C4a-OOH in the luciferase active site.

The cationic porphyrin TMPyP4 is a well-established DNA G-quadruplex (G4) binding ligand that can stabilize different topologies via multiple binding modes. However, TMPyP4 can have both a stabilizing and destabilizing effect on RNA G4 structures. The structural mechanisms that mediate RNA G4 unfolding remain unknown. Here, we report on the TMPyP4-induced RNA G4 unfolding mechanism studied by well-tempered metadynamics (WT-MetaD) with supporting biophysical experiments. The simulations predict a two-state mechanism of TMPyP4 interaction via a groove-bound and a top-face-bound conformation. The dynamics of TMPyP4 stacking on the top tetrad disrupts Hoogsteen H-bonds between guanine bases, resulting in the consecutive TMPyP4 intercalation from top-to-bottom G-tetrads. The results reveal a striking correlation between computational and experimental approaches and validate WT-MetaD simulations as a powerful tool for studying RNA G4-ligand interactions.

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants with increased transmissibility and pathology. In addition to this entrenched diversity, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriSΔ variant, originally identified as a viral subpopulation from SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike an eight amino-acid deletion encompassing a furin recognition motif and S1/S2 cleavage site. We elucidate the structure, function and molecular dynamics of this spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Our results reveal long-range allosteric communication between functional domains that differ in the wild-type and the deletion variant and support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host., (© 2022. The Author(s).)

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, shows the need for effective antiviral treatments. Here, we present a simulation study of the inhibition of the SARS-CoV-2 main protease (M pro ), a cysteine hydrolase essential for the life cycle of the virus. The free energy landscape for the mechanism of the inhibition process is explored by QM/MM umbrella sampling and free energy perturbation simulations at the M06-2X/MM level of theory for two proposed peptidyl covalent inhibitors that share the same recognition motif but feature distinct cysteine-targeting warheads. Regardless of the intrinsic reactivity of the modeled inhibitors, namely a Michael acceptor and a hydroxymethyl ketone activated carbonyl, our results confirm that the inhibitory process takes place by means of a two-step mechanism, in which the formation of an ion pair C145/H41 dyad precedes the protein-inhibitor covalent bond formation. The nature of this second step is strongly dependent on the functional groups in the warhead: while the nucleophilic attack of the C145 sulfur atom on the Cα of the double bond of the Michael acceptor takes place concertedly with the proton transfer from H41 to C β , in the compound with an activated carbonyl, the sulfur attacks the carbonyl carbon concomitant with a proton transfer from H41 to the carbonyl oxygen via the hydroxyl group. An analysis of the free energy profiles, structures along the reaction path, and interactions between the inhibitors and the different pockets of the active site on the protein shows a measurable effect of the warhead on the kinetics and thermodynamics of the process. These results and QM/MM methods can be used as a guide to select warheads to design efficient irreversible and reversible inhibitors of SARS-CoV-2 M pro ., Competing Interests: The authors declare no competing financial interest., (© 2021 American Chemical Society.)

Widespread bacterial resistance to carbapenem antibiotics is an increasing global health concern. Resistance has emerged due to carbapenem-hydrolyzing enzymes, including metallo-β-lactamases (MβLs), but despite their prevalence and clinical importance, MβL mechanisms are still not fully understood. Carbapenem hydrolysis by MβLs can yield alternative product tautomers with the potential to access different binding modes. Here, we show that a combined approach employing crystallography and quantum mechanics/molecular mechanics (QM/MM) simulations allow tautomer assignment in MβL:hydrolyzed antibiotic complexes. Molecular simulations also examine (meta)stable species of alternative protonation and tautomeric states, providing mechanistic insights into β-lactam hydrolysis. We report the crystal structure of the hydrolyzed carbapenem ertapenem bound to the L1 MβL from Stenotrophomonas maltophilia and model alternative tautomeric and protonation states of both hydrolyzed ertapenem and faropenem (a related penem antibiotic), which display different binding modes with L1. We show how the structures of both complexed β-lactams are best described as the (2 S )-imine tautomer with the carboxylate formed after β-lactam ring cleavage deprotonated. Simulations show that enamine tautomer complexes are significantly less stable (e.g., showing partial loss of interactions with the L1 binuclear zinc center) and not consistent with experimental data. Strong interactions of Tyr32 and one zinc ion (Zn1) with ertapenem prevent a C6 group rotation, explaining the different binding modes of the two β-lactams. Our findings establish the relative stability of different hydrolyzed (carba)penem forms in the L1 active site and identify interactions important to stable complex formation, information that should assist inhibitor design for this important antibiotic resistance determinant.

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication., 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., (© 2021 The Authors. Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.)

Zinc metalloproteins are ubiquitous, with protein zinc centers of structural and functional importance, involved in interactions with ligands and substrates and often of pharmacological interest. Biomolecular simulations are increasingly prominent in investigations of protein structure, dynamics, ligand interactions, and catalysis, but zinc poses a particular challenge, in part because of its versatile, flexible coordination. A computational workflow generating reliable models of ligand complexes of biological zinc centers would find broad application. Here, we evaluate the ability of alternative treatments, using (nonbonded) molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) at semiempirical (DFTB3) and density functional theory (DFT) levels of theory, to describe the zinc centers of ligand complexes of six metalloenzyme systems differing in coordination geometries, zinc stoichiometries (mono- and dinuclear), and the nature of interacting groups (specifically the presence of zinc-sulfur interactions). MM molecular dynamics (MD) simulations can overfavor octahedral geometries, introducing additional water molecules to the zinc coordination shell, but this can be rectified by subsequent semiempirical (DFTB3) QM/MM MD simulations. B3LYP/MM geometry optimization further improved the accuracy of the description of coordination distances, with the overall effectiveness of the approach depending upon factors, including the presence of zinc-sulfur interactions that are less well described by semiempirical methods. We describe a workflow comprising QM/MM MD using DFTB3 followed by QM/MM geometry optimization using DFT (e.g., B3LYP) that well describes our set of zinc metalloenzyme complexes and is likely to be suitable for creating accurate models of zinc protein complexes when structural information is more limited.

The emerging fields of citizen science and gamification reformulate scientific problems as games or puzzles to be solved. Through engaging the wider non-scientific community, significant breakthroughs may be made by analyzing citizen-gathered data. In parallel, recent advances in virtual reality (VR) technology are increasingly being used within a scientific context and the burgeoning field of interactive molecular dynamics in VR (iMD-VR) allows users to interact with dynamical chemistry simulations in real time. Here, we demonstrate the utility of iMD-VR as a medium for gamification of chemistry research tasks. An iMD-VR "game" was designed to encourage users to explore the reactivity of a particular chemical system, and a cohort of 18 participants was recruited to playtest this game as part of a user study. The reaction game encouraged users to experiment with making chemical reactions between a propyne molecule and an OH radical, and "molecular snapshots" from each game session were then compiled and used to map out reaction pathways. The reaction network generated by users was compared to existing literature networks demonstrating that users in VR capture almost all the important reaction pathways. Further comparisons between humans and an algorithmic method for guiding molecular dynamics show that through using citizen science to explore these kinds of chemical problems, new approaches and strategies start to emerge.

Activation heat capacity is emerging as a crucial factor in enzyme thermoadaptation, as shown by the non-Arrhenius behaviour of many natural enzymes. However, its physical origin and relationship to the evolution of catalytic activity remain uncertain. Here we show that directed evolution of a computationally designed Kemp eliminase reshapes protein dynamics, which gives rise to an activation heat capacity absent in the original design. These changes buttress transition-state stabilization. Extensive molecular dynamics simulations show that evolution results in the closure of solvent-exposed loops and a better packing of the active site. Remarkably, this gives rise to a correlated dynamical network that involves the transition state and large parts of the protein. This network tightens the transition-state ensemble, which induces a negative activation heat capacity and non-linearity in the activity-temperature dependence. Our results have implications for understanding enzyme evolution and suggest that selectively targeting the conformational dynamics of the transition-state ensemble by design and evolution will expedite the creation of novel enzymes., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Aurora kinase B plays an important role in the cell cycle to orchestrate the mitotic process. The amplification and overexpression of this kinase have been implicated in several human malignancies. Therefore, Aurora kinase B is a potential drug target for anticancer therapies. Here, we combine atom-based 3D-QSAR analysis and pharmacophore model generation to identify the principal structural features of acylureidoindolin derivatives that could potentially be responsible for the inhibition of Aurora kinase B. The selected CoMFA and CoMSIA model showed significant results with cross-validation values (q 2 ) of 0.68, 0.641 and linear regression values (r 2 ) of 0.971, 0.933 respectively. These values support the statistical reliability of our model. A pharmacophore model was also generated, incorporating features of reported crystal complex structures of Aurora kinase B. The pharmacophore model was used to screen commercial databases to retrieve potential lead candidates. The resulting hits were analyzed at each stage for diversity based on the pharmacophore model, followed by molecular docking and filtering based on their interaction with active site residues and 3D-QSAR predictions. Subsequently, MD simulations and binding free energy calculations were performed to test the predictions and to characterize interactions at the molecular level. The results suggested that the identified compounds retained the interactions with binding residues. Binding energy decomposition identified residues Glu155, Trp156 and Ala157 of site B and Leu83 and Leu207 of site C as major contributors to binding affinity, complementary to 3D-QSAR results. To best of our knowledge, this is the first comparison of WaterSwap field and 3D-QSAR maps. Overall, this integrated strategy provides a basis for the development of new and potential AK-B inhibitors and is applicable to other protein targets., (© 2021. The Author(s).)

Conformational sampling profoundly impacts the overall activity and temperature dependence of enzymes. Peroxidases have emerged as versatile platforms for high-value biocatalysis owing to their broad palette of potential biotransformations. Here, we explore the role of conformational sampling in mediating activity in the de novo peroxidase C45. We demonstrate that 2,2,2-triflouoroethanol (TFE) affects the equilibrium of enzyme conformational states, tending toward a more globally rigid structure. This is correlated with increases in both stability and activity. Notably, these effects are concomitant with the emergence of curvature in the temperature-activity profile, trading off activity gains at ambient temperature with losses at high temperatures. We apply macromolecular rate theory (MMRT) to understand enzyme temperature dependence data. These data point to an increase in protein rigidity associated with a difference in the distribution of protein dynamics between the ground and transition states. We compare the thermodynamics of the de novo enzyme activity to those of a natural peroxidase, horseradish peroxidase. We find that the native enzyme resembles the rigidified de novo enzyme in terms of the thermodynamics of enzyme catalysis and the putative distribution of protein dynamics between the ground and transition states. The addition of TFE apparently causes C45 to behave more like the natural enzyme. Our data suggest robust, generic strategies for improving biocatalytic activity by manipulating protein rigidity; for functional de novo protein catalysts in particular, this can provide more enzyme-like catalysts without further rational engineering, computational redesign, or directed evolution., Competing Interests: The authors declare no competing financial interest., (© 2021 The Authors. Published by American Chemical Society.)