Structural basis of metallo-ß-lactamase resistance to taniborbactam.

The design of inhibitors against metallo-ß-lactamases (MBLs), the largest family of carbapenemases, has been a strategic goal in designing novel antimicrobial therapies. In this regard, the development of bicyclic boronates, such as taniborbactam (TAN) and xeruborbactam, is a major achievement that may help in overcoming the threat of MBL-producing and carbapenem-resistant Gram-negative pathogens. Of concern, a recent report has shown that New Delhi MBL-9 (NDM-9) escapes the inhibitory action of TAN by a single amino acid substitution with respect to New Delhi MBL-1 (NDM-1), the most widely disseminated MBL. Here, we report a docking and computational analysis that identifies that "escape variants" against TAN can arise by disruption of the electrostatic interaction of negative charges in the active site loops of MBLs with the N-(2-aminoethyl)cyclohexylamine side chain of TAN. These changes result in non-productive binding modes of TAN that preclude reaction with the MBLs, a phenomenon that is not restricted to NDM-9. This analysis demonstrates that single amino acid substitutions in non-essential residues in MBL loops can unexpectedly elicit resistance to TAN.

The interaction of the azetidine thiazole side chain with the active site loop (ASL) 3 drives the evolution of IMP metallo-ß-lactamase against tebipenem.

Imipenemase (IMP) metallo-ß-lactamases (MBLs) hydrolyze almost all available ß-lactams including carbapenems and are not inhibited by any commercially available ß-lactamase inhibitor. Tebipenem (TP) pivoxil is the first orally available carbapenem and possesses a unique bicyclic azetidine thiazole moiety located at the R2 position. TP has potent in vitro activity against Enterobacterales producing extended-spectrum and/or AmpC ß-lactamases. Thus far, the activity of TP against IMP-producing strains is understudied. To address this knowledge gap, we explored the structure activity relationships of IMP MBLs by investigating whether IMP-6, IMP-10, IMP-25, and IMP-78 [MBLs with expanded hydrolytic activity against meropenem (MEM)] would demonstrate enhanced activity against TP. Most of the Escherichia coli DH10B strains expressing IMP-1 variants displayed a ≥twofold MIC difference between TP and MEM, while those expressing VIM or NDM variants demonstrated comparable MICs. Catalytic efficiency (kcat/KM) values for the TP hydrolysis by IMP-1, IMP-6, IMP-10, IMP-25, and IMP-78 were significantly lower than those obtained for MEM. Molecular dynamic simulations reveal that V67F and S262G substitutions (found in IMP-78) reposition active site loop 3, ASL-3, to better accommodate the bicyclic azetidine thiazole side chain, allowing microbiological/catalytic activity to approach that of comparison MBLs used in this study. These findings suggest that modifying the R2 side chain of carbapenems can significantly impact hydrolytic stability. Furthermore, changes in conformational dynamics due to single amino acid substitutions should be used to inform drug design of novel carbapenems.

Structural role of K224 in taniborbactam inhibition of NDM-1.

Taniborbactam (TAN; VNRX-5133) is a novel bicyclic boronic acid ß-lactamase inhibitor (BLI) being developed in combination with cefepime (FEP). TAN inhibits both serine and some metallo-ß-lactamases. Previously, the substitution R228L in VIM-24 was shown to increase activity against oxyimino-cephalosporins like FEP and ceftazidime (CAZ). We hypothesized that substitutions at K224, the homologous position in NDM-1, could impact FEP/TAN resistance. To evaluate this, a library of codon-optimized NDM K224X clones for minimum inhibitory concentration (MIC) measurements was constructed; steady-state kinetics and molecular docking simulations were next performed. Surprisingly, our investigation revealed that the addition of TAN restored FEP susceptibility only for NDM-1, as the MICs for the other 19 K224X variants remained comparable to those of FEP alone. Moreover, compared to NDM-1, all K224X variants displayed significantly lower MICs for imipenem, tebipenem, and cefiderocol (32-, 133-, and 33-fold lower, respectively). In contrast, susceptibility to CAZ was mostly unaffected. Kinetic assays with the K224I variant, the only variant with hydrolytic activity to FEP comparable to NDM-1, confirmed that the inhibitory capacity of TAN was modestly compromised (IC50 0.01 µM vs 0.14 µM for NDM-1). Lastly, structural modeling and docking simulations of TAN in NDM-1 and in the K224I variant revealed that the hydrogen bond between TAN's carboxylate with K224 is essential for the productive binding of TAN to the NDM-1 active site. In addition to the report of NDM-9 (E149K) as FEP/TAN resistant, this study demonstrates the fundamental role of single amino acid substitutions in the inhibition of NDM-1 by TAN.

Structure and function of crocodilian hemoglobins and allosteric regulation by chloride, ATP, and CO2.

Hemoglobins (Hbs) of crocodilians are reportedly characterized by unique mechanisms of allosteric regulatory control, but there are conflicting reports regarding the importance of different effectors, such as chloride ions, organic phosphates, and CO2. Progress in understanding the unusual properties of crocodilian Hbs has also been hindered by a dearth of structural information. Here, we present the first comparative analysis of blood properties and Hb structure and function in a phylogenetically diverse set of crocodilian species. We examine mechanisms of allosteric regulation in the Hbs of 13 crocodilian species belonging to the families Crocodylidae and Alligatoridae. We also report new amino acid sequences for the α- and ß-globins of these taxa, which, in combination with structural analyses, provide insights into molecular mechanisms of allosteric regulation. All crocodilian Hbs exhibited a remarkably strong sensitivity to CO2, which would permit effective O2 unloading to tissues in response to an increase in metabolism during intense activity and diving. Although the Hbs of all crocodilians exhibit similar intrinsic O2-affinities, there is considerable variation in sensitivity to Cl- ions and ATP, which appears to be at least partly attributable to variation in the extent of NH2-terminal acetylation. Whereas chloride appears to be a potent allosteric effector of all crocodile Hbs, ATP has a strong, chloride-independent effect on Hb-O2 affinity only in caimans. Modeling suggests that allosteric ATP binding has a somewhat different structural basis in crocodilian and mammalian Hbs.

Study of the role of Mg2+ in dsRNA processing mechanism by bacterial RNase III through QM/MM simulations.

The ribonuclease III (RNase III) cleaves dsRNA in specific positions generating mature RNAs. RNase III enzymes play important roles in RNA processing, post-transcriptional gene expression, and defense against viral infection. The enzyme's active site contains Mg2+ ions bound by a network of acidic residues and water molecules, but there is a lack of information about their specific roles. In this work, multiple steered molecular dynamics simulations at QM/MM level were performed to explore the hydrolysis reaction carried out by the enzyme. Free energy profiles modifying the features of the active site are obtained and the role of Mg2+ ions, the solvent molecules and the residues of the active site are discussed in detail. Our results show that Mg2+ ions carry out different roles in the hydrolysis process positioning the substrate for the attack from a coordinated nucleophile and activating it to perform hydrolysis reaction, cleaving the dsRNA backbone in a SN2 substitution. In addition, water molecules present in the active site lower the energy barrier of the process. RNase III hydrolyzes dsRNA to generate mature RNAs. For this purpose, its active site contains Mg2+ which has an important role during the reaction. Results show that the Mg2+ activates the solvent molecule that produces the nucleophilic attack and the surrounding waters contribute significantly to the hydrolysis process.

Conformational and Reaction Dynamic Coupling in Histidine Kinases: Insights from Hybrid QM/MM Simulations.

Histidine kinases (HK) of bacterial two-component systems represent a hallmark of allosterism in proteins, being able to detect a signal through the sensor domain and transmit this information through the protein matrix to the kinase domain which, once active, autophosphorylates a specific histidine residue. Inactive-to-active transition results in a large conformational change that moves the kinase on top of the histidine. In the present work, we use several molecular simulation techniques (Molecular Dynamics, Hybrid QM/MM, and constant pH molecular dynamics) to study the activation and autophosphorylation reactions in L. plantarum WalK, a cis-acting HK. In agreement with previous results, we show that the chemical step requires tight coupling with the conformational step in order to maintain the histidine phosphoacceptor in the correct tautomeric state, with a reactive δ-nitrogen. During the conformational transition, the kinase domain is never released and walks along the HK helix axis, breaking and forming several conserved residue-based contacts. The phosphate transfer reaction is concerted in the transition state region and is catalyzed through the stabilization of the negative developing charge of transferring phosphate along the reaction.

Crystal Structure of the Metallo-ß-Lactamase GOB in the Periplasmic Dizinc Form Reveals an Unusual Metal Site.

Metallo-beta-lactamases (MBLs) are broad-spectrum, Zn(II)-dependent lactamases able to confer resistance to virtually every ß-lactam antibiotic currently available. The large diversity of active-site structures and metal content among MBLs from different sources has limited the design of a pan-MBL inhibitor. GOB-18 is a divergent MBL from subclass B3 that is expressed by the opportunistic Gram-negative pathogen Elizabethkingia meningoseptica This MBL is atypical, since several residues conserved in B3 enzymes (such as a metal ligand His) are substituted in GOB enzymes. Here, we report the crystal structure of the periplasmic di-Zn(II) form of GOB-18. This enzyme displays a unique active-site structure, with residue Gln116 coordinating the Zn1 ion through its terminal amide moiety, replacing a ubiquitous His residue. This situation contrasts with that of B2 MBLs, where an equivalent His116Asn substitution leads to a di-Zn(II) inactive species. Instead, both the mono- and di-Zn(II) forms of GOB-18 are active against penicillins, cephalosporins, and carbapenems. In silico docking and molecular dynamics simulations indicate that residue Met221 is not involved in substrate binding, in contrast to Ser221, which otherwise is conserved in most B3 enzymes. These distinctive features are conserved in recently reported GOB orthologues in environmental bacteria. These findings provide valuable information for inhibitor design and also posit that GOB enzymes have alternative functions.

dsRNA-protein interactions studied by molecular dynamics techniques. Unravelling dsRNA recognition by DCL1.

Double stranded RNA (dsRNA) participates in several biological processes, where RNA molecules acquire secondary structure inside the cell through base complementarity. The double stranded RNA binding domain (dsRBD) is one of the main protein folds that is able to recognize and bind to dsRNA regions. The N-terminal dsRBD of DCL1 in Arabidopsis thaliana (DCL1-1), in contrast to other studied dsRBDs, lacks a stable structure, behaving as an intrinsically disordered protein. DCL1-1 does however recognize dsRNA by acquiring a canonical fold in the presence of its substrate. Here we present a detailed modeling and molecular dynamics study of dsRNA recognition by DCL1-1. We found that DCL1-1 forms stable complexes with different RNAs and we characterized the residues involved in binding. Although the domain shows a binding loop substantially shorter than other homologs, it can still interact with the dsRNA and results in bending of the dsRNA A-type helix. Furthermore, we found that R8, a non-conserved residue located in the first dsRNA binding region, recognizes preferentially mismatched base pairs. We discuss our findings in the context of the function of DCL1-1 within the microRNA processing complex.

Structural modeling of a novel membrane-bound globin-coupled sensor in Geobacter sulfurreducens.

Globin-coupled sensors (GCS) usually consist of three domains: a sensor/globin, a linker, and a transmitter domain. The globin domain (GD), activated by ligand binding and/or redox change, induces an intramolecular signal transduction resulting in a response of the transmitter domain. Depending on the nature of the transmitter domain, GCSs can have different activities and functions, including adenylate and di-guanylate cyclase, histidine kinase activity, aerotaxis and/or oxygen sensing function. The gram-negative delta-proteobacterium Geobacter sulfurreducens expresses a protein with a GD covalently linked to a four transmembrane domain, classified, by sequence similarity, as GCS (GsGCS). While its GD is fully characterized, not so its transmembrane domain, which is rarely found in the globin superfamily. In the present work, GsGCS was characterized spectroscopically and by native ion mobility-mass spectrometry in combination with cryo-electron microscopy. Although lacking high resolution, the oligomeric state and the electron density map were valuable for further rational modeling of the full-length GsGCS structure. This model demonstrates that GsGCS forms a transmembrane domain-driven tetramer with minimal contact between the GDs and with the heme groups oriented outward. This organization makes an intramolecular signal transduction less likely. Our results, including the auto-oxidation rate and redox potential, suggest a potential role for GsGCS as redox sensor or in a membrane-bound e-/H+ transfer. As such, GsGCS might act as a player in connecting energy production to the oxidation of organic compounds and metal reduction. Database searches indicate that GDs linked to a four or seven helices transmembrane domain occur more frequently than expected.

Identification of key sequence features required for microRNA biogenesis in plants.

MicroRNAs (miRNAs) are endogenous small RNAs of ∼21 nt that regulate multiple biological pathways in multicellular organisms. They derive from longer transcripts that harbor an imperfect stem-loop structure. In plants, the ribonuclease type III DICER-LIKE1 assisted by accessory proteins cleaves the precursor to release the mature miRNA. Numerous studies highlight the role of the precursor secondary structure during plant miRNA biogenesis; however, little is known about the relevance of the precursor sequence. Here, we analyzed the sequence composition of plant miRNA primary transcripts and found specifically located sequence biases. We show that changes in the identity of specific nucleotides can increase or abolish miRNA biogenesis. Most conspicuously, our analysis revealed that the identity of the nucleotides at unpaired positions of the precursor plays a crucial role during miRNA biogenesis in Arabidopsis.