QM/MM Simulations Reveal the Determinants of Carbapenemase Activity in Class A β-Lactamases

JAK3 Inhibitors: Covalent and Noncovalent Interactions of a Cyanamide Group Investigated by Multiscale Free-Energy Simulations

Janus kinase type 3 (JAK3), an emerging target for treating autoimmune diseases, possesses a front pocket cysteine that is targeted by covalent modifiers, best represented by the marketed drug ritlecitinib (1). Recently, 2,3-dihydro-1H-inden-1-ylcyanamides have been developed as novel JAK3 inhibitors. Among them, the N-(6-(7H-pyrrolo[2,3-d]pyrimidin-4-yl)-2,3-dihydro-1H-inden-1-yl)cyanamide inhibitor (2) and its methylated analogue (3), while being potent inhibitors, displayed different mechanisms of action (covalent vs noncovalent) and binding modes (Casimiro-Garcia et al., J Med Chem 2018). Prompted by this intriguing behavior, we applied a multiscale approach to characterize the reaction mechanism between the JAK3 front-pocket Cys909 and cyanamide-based inhibitors. Quantum mechanics/molecular mechanics simulations showed that 2 can readily form an isothiourea adduct with the Cys909 only when a conserved water molecule assists the reaction as a proton shuttle and that methylation of the 2,3-dihydro-1H-inden-1-ylcyanamide moiety of 2 hampers the isothiourea formation by displacing this water molecule. Metadynamics and thermodynamic integration simulations were applied to investigate the relative abundance of alternative poses accessible to 2,3-dihydro-1H-inden-1-ylcyanamides, explaining the effect of methylation on the relative binding mode preference. This multiscale approach provides new chemical insights into the mechanism of action of cyanamide inhibitors and emerges as an effective protocol to investigate the interaction between drugs and molecular targets.

Dynamical responses predict a distal site that modulates activity in an antibiotic resistance enzyme

β-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-2G89D 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.

Electric Fields Are a Key Determinant of Carbapenemase Activity in Class A β-Lactamases

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.

Cooperative Conformational Transitions Underpin the Activation Heat Capacity in the Temperature Dependence of Enzyme Catalysis

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 (ΔCP⧧) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of ΔCP⧧ (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.

Enhanced activity against a third-generation cephalosporin by destabilization of the active site of a class A beta-lactamase

The β-lactamase BlaC conveys resistance to a broad spectrum of β-lactam antibiotics to its host Mycobacterium tuberculosis but poorly hydrolyzes third-generation cephalosporins, such as ceftazidime. Variants of other β-lactamases have been reported to gain activity against ceftazidime at the cost of the native activity. To understand this trade-off, laboratory evolution was performed, screening for enhanced ceftazidime activity. The variant BlaC Pro167Ser shows faster breakdown of ceftazidime, poor hydrolysis of ampicillin and only moderately reduced activity against nitrocefin. NMR spectroscopy, crystallography and kinetic assays demonstrate that the resting state of BlaC P167S exists in an open and a closed state. The open state is more active in the hydrolysis of ceftazidime. In this state the catalytic residue Glu166, generally believed to be involved in the activation of the water molecule required for deacylation, is rotated away from the active site, suggesting it plays no role in the hydrolysis of ceftazidime. In the closed state, deacylation of the BlaC-ceftazidime adduct is slow, while hydrolysis of nitrocefin, which requires the presence of Glu166 in the active site, is barely affected, providing a structural explanation for the trade-off in activities.

Tautomer-Specific Deacylation and Ω-Loop Flexibility Explain the Carbapenem-Hydrolyzing Broad-Spectrum Activity of the KPC-2 β-Lactamase

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 (kcat), 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-(2R) 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-(2R) 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-(2R) 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.

Mechanistic Modeling of Lys745 Sulfonylation in EGFR C797S Reveals Chemical Determinants for Inhibitor Activity and Discriminates Reversible from Irreversible Agents

Targeted covalent inhibitors hold promise for drug discovery, particularly for kinases. Targeting the catalytic lysine of epidermal growth factor receptor (EGFR) has attracted attention as a new strategy to overcome resistance due to the emergence of C797S mutation. Sulfonyl fluoride derivatives able to inhibit EGFR^{L858R/T790M/C797S} by sulfonylation of Lys745 have been reported. However, atomistic details of this process are still poorly understood. Here, we describe the mechanism of inhibition of an innovative class of compounds that covalently engage the catalytic lysine of EGFR, through a sulfur(VI) fluoride exchange (SuFEx) process, with the help of hybrid quantum mechanics/molecular mechanics (QM/MM) and path collective variables (PCVs) approaches. Our simulations identify the chemical determinants accounting for the irreversible activity of agents targeting Lys745 and provide hints for the further optimization of sulfonyl fluoride agents.

Unveiling the structural features that regulate carbapenem deacylation in KPC-2 through QM/MM and interpretable machine learning

Resistance to carbapenem β-lactams presents major clinical and economical challenges for the treatment of pathogen infections. The fast hydrolysis of carbapenems by carbapenemase-producing bacterial strains enables the effective deactivation of carbapenem antibiotics. In this study, we aim to unravel the structural features that distinguish the notable deacylation activity of carbapenemases. The deacylation reactions between imipenem (IPM) and the KPC-2 class A serine-based β-lactamases (ASβLs) are modeled with combined quantum mechanical/molecular mechanical (QM/MM) minimum energy pathway (MEP) calculations and interpretable machine-learning (ML) methods. We first applied a dual-level computational protocol to achieve fast sampling of QM/MM MEPs. A tree-based ensemble ML model was employed to learn the MEP activation barriers from the conformational features of the KPC-2/IPM active site. The barrier-predicting model was then unboxed using the Shapley additive explanation (SHAP) importance attribution methods to derive mechanistic insights, which were also verified by additional QM/MM wavefunction analysis. Essentially, we show that potential hydrogen bonding interactions of the general base and the tautomerization states of the carbapenem pyrroline ring could concertedly regulate the activation barrier of KPC-2/IPM deacylation. Nonetheless, we demonstrate the efficacy of interpretable ML to assist the analysis of QM/MM simulation data for robust extraction of human-interpretable mechanistic insights.