Mechanisms of Antibiotic Resistance: QM/MM Modeling of the Acylation Reaction of a Class A β-Lactamase with Benzylpenicillin

Decoding Drug Discovery: Exploring A-to-Z In Silico Methods for Beginners

The drug development process is a critical challenge in the pharmaceutical industry due to its time-consuming nature and the need to discover new drug potentials to address various ailments. The initial step in drug development, drug target identification, often consumes considerable time. While valid, traditional methods such as in vivo and in vitro approaches are limited in their ability to analyze vast amounts of data efficiently, leading to wasteful outcomes. To expedite and streamline drug development, an increasing reliance on computer-aided drug design (CADD) approaches has merged. These sophisticated in silico methods offer a promising avenue for efficiently identifying viable drug candidates, thus providing pharmaceutical firms with significant opportunities to uncover new prospective drug targets. The main goal of this work is to review in silico methods used in the drug development process with a focus on identifying therapeutic targets linked to specific diseases at the genetic or protein level. This article thoroughly discusses A-to-Z in silico techniques, which are essential for identifying the targets of bioactive compounds and their potential therapeutic effects. This review intends to improve drug discovery processes by illuminating the state of these cutting-edge approaches, thereby maximizing the effectiveness and duration of clinical trials for novel drug target investigation.

The Evolution of the Acylation Mechanism in β-Lactamase and Rapid Protein Dynamics

β-Lactamases are a class of well-studied enzymes that are known to have existed since billions of years ago, starting as a defense mechanism to stave off competitors and are now enzymes responsible for antibiotic resistance. Using ancestral sequence reconstruction, it is possible to study the crystal structure of a laboratory resurrected 2-3 billion year-old β-lactamase. Comparing the ancestral enzyme to its modern counterpart, a TEM-1 β-lactamase, the structural changes are minor, and it is probable that dynamic effects play an important role in the evolution of function. We used molecular dynamics simulations and employed transition path sampling methods to identify the presence of rate-enhancing dynamics at the femtosecond level in both systems, found that these fast motions are more efficiently coordinated in the modern enzyme, and examined how specific dynamics can pinpoint evolutionary effects that are essential for improving enzymatic catalysis.

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.

Structural Insights into the Catalytic Activity of Cyclobacterium marinum N-Acetylglucosamine Deacetylase

N-Acetylglucosamine deacetylase from Cyclobacterium marinum (CmCBDA) is a highly effective and selective biocatalyst for the production of d-glucosamine (GlcN) from N-acetylglucosamine (GlcNAc). However, the underlying catalytic mechanism remains elusive. Here, we show that CmCBDA is a metalloenzyme with a preference for Ni2+ over Mn2+. Crystal structures of CmCBDA in complex with Ni2+ and Mn2+ revealed slight remodeling of the CmCBDA active site by the metal ions. We also demonstrate that CmCBDA exists as a mixture of homodimers and monomers in solution, and dimerization is indispensable for catalytic activity. A mutagenesis analysis also indicated that the active site residues Asp22, His72, and His143 as well as the residues involved in dimerization, Pro52, Trp53, and Tyr55, are essential for catalytic activity. Furthermore, a mutation on the protein surface, Lys219Glu, resulted in a 2.3-fold improvement in the deacetylation activity toward GlcNAc. Mechanistic insights obtained here may facilitate the development of CmCBDA variants with higher activities.

Mutagenesis and structural analysis reveal the CTX-M β-lactamase active site is optimized for cephalosporin catalysis and drug resistance

CTX-M β-lactamases are a widespread source of resistance to β-lactam antibiotics in Gram-negative bacteria. These enzymes readily hydrolyze penicillins and cephalosporins, including oxyimino-cephalosporins such as cefotaxime. To investigate the preference of CTX-M enzymes for cephalosporins, we examined eleven active-site residues in the CTX-M-14 β-lactamase model system by alanine mutagenesis to assess the contribution of the residues to catalysis and specificity for the hydrolysis of the penicillin, ampicillin, and the cephalosporins cephalothin and cefotaxime. Key active site residues for class A β-lactamases, including Lys73, Ser130, Asn132, Lys234, Thr216, and Thr235, contribute significantly to substrate binding and catalysis of penicillin and cephalosporin substrates in that alanine substitutions decrease both kcat and kcat/KM. A second group of residues, including Asn104, Tyr105, Asn106, Thr215, and Thr216, contribute only to substrate binding, with the substitutions decreasing only kcat/KM. Importantly, calculating the average effect of a substitution across the 11 active-site residues shows that the most significant impact is on cefotaxime hydrolysis while ampicillin hydrolysis is least affected, suggesting the active site is highly optimized for cefotaxime catalysis. Furthermore, we determined X-ray crystal structures for the apo-enzymes of the mutants N106A, S130A, N132A, N170A, T215A, and T235A. Surprisingly, in the structures of some mutants, particularly N106A and T235A, the changes in structure propagate from the site of substitution to other regions of the active site, suggesting that the impact of substitutions is due to more widespread changes in structure and illustrating the interconnected nature of the active site.

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.

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.

Mechanistic Investigation of Lysine-Targeted Covalent Inhibition of PI3Kδ via ONIOM QM:QM Computations

Phosphoinositide 3-kinase (PI3K) enzymes are important drug targets, especially in oncology, and several inhibitors are currently under investigation in clinical trials for the treatment of lymphocytic leukemia, follicular lymphoma, breast, thyroid, colorectal, and lung cancer. Targeted covalent inhibitors hold significant promise for drug discovery research especially for kinases. Targeting the lysine residues attracts attention as a new strategy in designing targeted covalent inhibitors, since the lysine residue provides several advantages over the traditional cysteine residue. Recently, new highly selective covalent inhibitors of PI3Kδ with activated ester warheads, targeting the conserved Lys779 residue, were reported. Based on the observed kinetics, a covalent inhibition mechanism was proposed, but the atomistic details of the reaction are still not understood. Therefore, in the present work, we have conducted quantum chemical ONIOM M06-2X/6-31+G(d,p):PM6 calculations on the active site cluster structure of PI3Kδ to elucidate the microscopic details of the mechanism of the aminolysis reaction between Lys779 and the ester inhibitors. Our calculations clearly discriminate the noncovalent methyl ester inhibitor and the covalent inhibitors with activated phenolic esters. For the representative p-NO₂, p-F, p-H, and p-OCH₃ phenolic esters, the Gibbs free energy profiles of alternative mechanistic paths through either Asp782 or Asp911 demonstrate the modulatory role of active site aspartate residues. The most plausible path alters depending on the electron-withdrawing/donating nature of the p-substituted phenolate leaving group. Inhibitors with sufficiently strong electron-withdrawing group prefer direct dissociation of the leaving group from the tetrahedral zwitterion intermediate, while the ones with electron-donating group favor the formation of a neutral tetrahedral intermediate prior to the dissociation. The relative Gibbs free energy barriers of p-NO₂ < p-F < p-H < p-OCH₃ substituted phenyl esters display the same qualitative trend as the experimentally measured k_inact/K_I values. Our results provide in depth insight into the mechanism, which can pave the way for optimizing the inhibitor efficiency.

Breathable and Stretchable Dressings for Accelerating Healing of Infected Wounds

Multidrug-resistant (MDR) bacteria-infected wounds are challenging issues that threaten human health. Herein, L-thioproline (T) and Boc-capped L-thioproline (BT)-decorated gold nanoparticles (TBT-GNPs) with potent antibacterial activity against MDR bacteria are reported. The TBT-GNPs are composited with bacterial cellulose to form wound dressings which show excellent antimicrobial performance both in vitro and in vivo. Moreover, this dressing is both breathable and stretchable which is favorable for gas exchange to accelerate the wound healing. This work is insightful for developing multifunctional dressings to satisfy the clinical requirements.

Mechanistic Insights into Enzyme Catalysis from Explaining Machine-Learned Quantum Mechanical and Molecular Mechanical Minimum Energy Pathways

With the increasing popularity of machine learning (ML) applications, the demand for explainable artificial intelligence techniques to explain ML models developed for computational chemistry has also emerged. In this study, we present the development of the Boltzmann-weighted cumulative integrated gradients (BCIG) approach for effective explanation of mechanistic insights into ML models trained on high-level quantum mechanical and molecular mechanical (QM/MM) minimum energy pathways. Using the acylation reactions of the Toho-1 β-lactamase and two antibiotics (ampicillin and cefalexin) as the model systems, we show that the BCIG approach could quantitatively attribute the energetic contribution in one system and the relative reactivity of individual steps across different systems to specific chemical processes such as the bond making/breaking and proton transfers. The proposed BCIG contribution attribution method quantifies chemistry-interpretable insights in terms of contributions from each elementary chemical process, which is in agreement with the validating QM/MM calculations and our intuitive mechanistic understandings of the model reactions.

Quantum Mechanics/Molecular Mechanics Studies on the Catalytic Mechanism of a Novel Esterase (FmtA) of Staphylococcus aureus

FmtA is a novel esterase that shares the penicillin-binding protein (PBP) core structural folding but found to hydrolyze the removal of d-Ala from teichoic acids. Molecular docking, dynamics, and MM-GBSA of FmtA and its variants S127A, K130A, Y211A, D213A, and K130AY211A, in the presence or absence of wall teichoic acid (WTA), suggest that active site residues S127, K130, Y211, D213, N343, and G344 play a role in substrate binding. Quantum mechanics (QM)/molecular mechanics (MM) calculations reveal that during WTA catalysis, K130 deprotonates S127, and the nucleophilic S127 attacks the carbonyl carbon of d-Ala bound to WTA. The tetrahedral intermediate (TI) complex is stabilized by hydrogen bonding to the oxyanion holes. The TI complex displays a high energy gap and collapses to an energetically favorable acyl-enzyme complex.

Deciphering the Mechanism of Binding Selectivity of Chlorofluoroacetamide-Based Covalent Inhibitors toward L858R/T790M Resistance Mutation

Covalent modification of the oncogenic mutant epidermal growth factor receptor (EGFR) by small molecules is an efficient strategy for achieving an enhanced and sustained pharmacological effect in the treatment of non-small-cell lung cancer. NSP-037 (18), an irreversible inhibitor of the L858R/T790M double-mutant EGFR (EGFR^DM) using α-chlorofluoroacetamide (CFA) as a novel warhead, has seven times the inhibition selectivity for EGFR^DM over the wild type (EGFR^WT), as compared to clinically approved osimertinib (7). Here, we employ multiple computational approaches to elucidate the mechanism underlining this improved selectivity, as well as the effect of CFA on the selectivity enhancement of inhibitor 18 over 7. We find that EGFR^DM undergoes significantly larger conformational changes than EGFR^WT upon binding to 18. The conformational stability of the diamine side chain and the CFA motif of 18 in the orthosteric site of EGFR^DM is identified as key for the disparate binding mechanism and inhibitory prowess of 18 with respect to EGFR^WT and EGFR^DM and 18's higher selectivity than 7. The binding free energy of the 18-bound complexes is -6.38 kcal/mol greater than that of the 7-bound complexes, explaining the difference in selectivity of these inhibitors. Further, free energy decomposition analysis indicates that the electrostatic contribution of key residues plays an important role in the 18-bound complexes. QM/MM calculations show that the most favored mechanism for the Cys797 alkylation reaction is the direct displacement mechanism through a CFA-based inhibitor, producing a reaction with the lowest energy barrier and most stable product.

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5′-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2′-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2′-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.

QM/MM modeling of class A β-lactamases reveals distinct acylation pathways for ampicillin and cefalexin

Efficient mechanism-based design of antibiotics that are not susceptible to β-lactamases is hindered by the lack of comprehensive knowledge on the energetic landscapes for the hydrolysis of various β-lactams. Herein, we adopted efficient quantum mechanics/molecular mechanics simulations to explore the acylation reaction catalyzed by CTX-M-44 (Toho-1) β-lactamase. We show that the catalytic pathways for β-lactam hydrolysis are correlated to substrate scaffolds: using Glu166 as the only general base for acylation is viable for ampicillin but prohibitive for cefalexin. The present computational workflow provides quantitative insights to facilitate the optimization of future β-lactam antibiotics.

Allosteric communication in class A β-lactamases occurs via cooperative coupling of loop dynamics

Understanding allostery in enzymes and tools to identify it offer promising alternative strategies to inhibitor development. Through a combination of equilibrium and nonequilibrium molecular dynamics simulations, we identify allosteric effects and communication pathways in two prototypical class A β-lactamases, TEM-1 and KPC-2, which are important determinants of antibiotic resistance. The nonequilibrium simulations reveal pathways of communication operating over distances of 30 Å or more. Propagation of the signal occurs through cooperative coupling of loop dynamics. Notably, 50% or more of clinically relevant amino acid substitutions map onto the identified signal transduction pathways. This suggests that clinically important variation may affect, or be driven by, differences in allosteric behavior, providing a mechanism by which amino acid substitutions may affect the relationship between spectrum of activity, catalytic turnover, and potential allosteric behavior in this clinically important enzyme family. Simulations of the type presented here will help in identifying and analyzing such differences.

Active site dynamics and catalytic mechanism in arabinan hydrolysis catalyzed by GH43 endo-arabinanase from QM/MM molecular dynamics simulation and potential energy surface

The endo-1,5-α-L-arabinanases, belonging to glycoside hydrolase family 43 (GH43), catalyse the hydrolysis of α-1,5-arabinofuranosidic bonds in arabinose-containing polysaccharides. These enzymes are proposed targets for industrial and medical applications. Here, molecular dynamics (MD), potential energy surface and free energy (potential of mean force) simulations are undertaken using hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to understand the active site dynamics, catalytic mechanism and the electrostatic influence of active site residues of the GH43 endo-arabinanase from G. stearothermophilus. The calculated results give support to the single-displacement mechanism proposed for the inverting GH43 enzymes: first a proton is transferred from the general acid E201 to the substrate, followed by a nucleophilic attack by water, activated by the general base D27, on the anomer carbon. A conformational change (²E ↔E₃ ↔ ⁴E) in the -1 sugar ring is observed involving a transition state featuring an oxocarbenium ion character. Residues D87, K106, H271 are highlighted as potential targets for future mutation experiments in order to increase the efficiency of the reaction. To our knowledge, this is the first QM/MM study providing molecular insights into the glycosidic bond hydrolysis of a furanoside substrate by an inverting GH in solution.Communicated by Ramaswamy H. Sarma.

Antimicrobial Peptides and Copper(II) Ions: Novel Therapeutic Opportunities

The emergence of new pathogens and multidrug resistant bacteria is an important public health issue that requires the development of novel classes of antibiotics. Antimicrobial peptides (AMPs) are a promising platform with great potential for the identification of new lead compounds that can combat the aforementioned pathogens due to their broad-spectrum antimicrobial activity and relatively low rate of resistance emergence. AMPs of multicellular organisms made their debut four decades ago thanks to ingenious researchers who asked simple questions about the resistance to bacterial infections of insects. Questions such as "Do fruit flies ever get sick?", combined with pioneering studies, have led to an understanding of AMPs as universal weapons of the immune system. This review focuses on a subclass of AMPs that feature a metal binding motif known as the amino terminal copper and nickel (ATCUN) motif. One of the metal-based strategies of hosts facing a pathogen, it includes wielding the inherent toxicity of copper and deliberately trafficking this metal ion into sites of infection. The sudden increase in the concentration of copper ions in the presence of ATCUN-containing AMPs (ATCUN-AMPs) likely results in a synergistic interaction. Herein, we examine common structural features in ATCUN-AMPs that exist across species, and we highlight unique features that deserve additional attention. We also present the current state of knowledge about the molecular mechanisms behind their antimicrobial activity and the methods available to study this promising class of AMPs.

Unraveling the energetic significance of chemical events in enzyme catalysis via machine-learning based regression approach

The bacterial enzyme class of β-lactamases are involved in benzylpenicillin acylation reactions, which are currently being revisited using hybrid quantum mechanical molecular mechanical (QM/MM) chain-of-states pathway optimizations. Minimum energy pathways are sampled by reoptimizing pathway geometry under different representative protein environments obtained through constrained molecular dynamics simulations. Predictive potential energy surface models in the reaction space are trained with machine-learning regression techniques. Herein, using TEM-1/benzylpenicillin acylation reaction as the model system, we introduce two model-independent criteria for delineating the energetic contributions and correlations in the predicted reaction space. Both methods are demonstrated to effectively quantify the energetic contribution of each chemical process and identify the rate limiting step of enzymatic reaction with high degrees of freedom. The consistency of the current workflow is tested under seven levels of quantum chemistry theory and three non-linear machine-learning regression models. The proposed approaches are validated to provide qualitative compliance with experimental mutagenesis studies.

Mechanistic Study of Potent Fluorinated EGFR Kinase Inhibitors with a Quinazoline Scaffold against L858R/T790M/C797S Resistance Mutation: Unveiling the Fluorine Substituent Cooperativity Effect on the Inhibitory Activity

Fluorination has considerable potential with regard to the design of kinase inhibitors for anticarcinoma therapy. It was recently reported that fluorination increases the potency of inhibitors of the epidermal growth factor receptor (EGFR), mutations of which have been linked specifically to nonsmall-cell lung cancer. For the L858R/T790M/C797S triplet mutant (EGFR^TM), a difluorinated inhibitor, 25g, was found to have 4.23 times greater potency against the EGFR^TM than an unfluorinated inhibitor, 25a. This discovery necessitates a rational explanation for the underlying inhibitory mechanisms. Here, we apply multiple computational approaches to explore, validate, and differentiate the binding modes of 25a and 25g in the EGFR^TM and investigate the cooperativity effect of fluorine substituents on the inhibitory activity. Our results showed that the EGFR^TM in the presence of 25g undergoes a series of conformational changes that favor inhibitor binding to both the active and allosteric sites. Further, the cooperativity effect of fluorine substituents is positive: the complex stability is increased by each additional fluorine substituent. Estimated binding free energies show good correlation with the experimental biological activity. Subsequently, the decomposition energy analysis revealed that the van der Waals interaction is the principal force contributing to variations in the binding affinities of 25a and 25g to the EGFR^TM. Per-residue energy-based hierarchical clustering analysis suggests that three hot-spot residues, L718, K745, and D855, are the key in achieving optimal binding modes for 25g with higher affinity in the EGFR^TM compared to 25a. This study provides a rationale for the superior EGFR^TM-inhibitory potency exhibited by 25g over 25a, which is expected to be useful for the future rational structure-based design of novel EGFR^TM inhibitors with improved potency and selectivity.

Acylation and deacylation mechanism and kinetics of penicillin G reaction with Streptomyces R61 DD-peptidase

Two quantum mechanical (QM)-cluster models are built for studying the acylation and deacylation mechanism and kinetics of Streptomyces R61 DD-peptidase with the penicillin G at atomic level detail. DD-peptidases are bacterial enzymes involved in the cross-linking of peptidoglycan to form the cell wall, necessary for bacterial survival. The cross-linking can be inhibited by antibiotic beta-lactam derivatives through acylation, preventing the acyl-enzyme complex from undergoing further deacylation. The deacylation step was predicted to be rate-limiting. Transition state and intermediate structures are found using density functional theory in this study, and thermodynamic and kinetic properties of the proposed mechanism are evaluated. The acyl-enzyme complex is found lying in a deep thermodynamic sink, and deacylation is indeed the severely rate-limiting step, leading to suicide inhibition of the peptidoglycan cross-linking. The usage of QM-cluster models is a promising technique to understand, improve, and design antibiotics to disrupt function of the Streptomyces R61 DD-peptidase.

The hydrolytic water molecule of Class A β-lactamase relies on the acyl-enzyme intermediate ES* for proper coordination and catalysis

Serine-based β-lactamases of Class A, C and D all rely on a key water molecule to hydrolyze and inactivate β-lactam antibiotics. This process involves two conserved catalytic steps. In the first acylation step, the β-lactam antibiotic forms an acyl-enzyme intermediate (ES*) with the catalytic serine residue. In the second deacylation step, an activated water molecule serves as nucleophile (WAT_Nu) to attack ES* and release the inactivated β-lactam. The coordination and activation of WAT_Nu is not fully understood. Using time-resolved x-ray crystallography and QM/MM simulations, we analyzed three intermediate structures of Class A β-lactamase PenP as it slowly hydrolyzed cephaloridine. WAT_Nu is centrally located in the apo structure but becomes slightly displaced away by ES* in the post-acylation structure. In the deacylation structure, WAT_Nu moves back and is positioned along the Bürgi–Dunitz trajectory with favorable energetic profile to attack ES*. Unexpectedly, WAT_Nu is also found to adopt a catalytically incompetent conformation in the deacylation structure forming a hydrogen bond with ES*. Our results reveal that ES* plays a significant role in coordinating and activating WAT_Nu through subtle yet distinct interactions at different stages of the catalytic process. These interactions may serve as potential targets to circumvent β-lactamase-mediated antibiotic resistance.

Computational Simulations to Guide Enzyme-Mediated Prodrug Activation

Prodrugs are designed to improve pharmaceutical/biopharmaceutical characteristics, pharmacokinetic/pharmacodynamic properties, site-specificity, and more. A crucial step in successful prodrug is its activation, which releases the active parent drug, exerting a therapeutic effect. Prodrug activation can be based on oxidation/reduction processes, or through enzyme-mediated hydrolysis, from oxidoreductases (i.e., Cytochrome P450) to hydrolytic enzymes (i.e., carboxylesterase). This study provides an overview of the novel in silico methods for the optimization of enzyme-mediated prodrug activation. Computational methods simulating enzyme-substrate binding can be simpler like molecular docking, or more complex, such as quantum mechanics (QM), molecular mechanics (MM), and free energy perturbation (FEP) methods such as molecular dynamics (MD). Examples for MD simulations used for elucidating the mechanism of prodrug (losartan, paclitaxel derivatives) metabolism via CYP450 enzyme are presented, as well as an MD simulation for optimizing linker length in phospholipid-based prodrugs. Molecular docking investigating quinazolinone prodrugs as substrates for alkaline phosphatase is also presented, as well as QM and MD simulations used for optimal fit of different prodrugs within the human carboxylesterase 1 catalytical site. Overall, high quality computational simulations may show good agreement with experimental results, and should be used early in the prodrug development process.

Assessment of the Cruzain Cysteine Protease Reversible and Irreversible Covalent Inhibition Mechanism

Reversible and irreversible covalent ligands are advanced cysteine protease inhibitors in the drug development pipeline. K777 is an irreversible inhibitor of cruzain, a necessary enzyme for the survival of the Trypanosoma cruzi (T. cruzi) parasite, the causative agent of Chagas disease. Despite their importance, irreversible covalent inhibitors are still often avoided due to the risk of adverse effects. Herein, we replaced the K777 vinyl sulfone group with a nitrile moiety to obtain a reversible covalent inhibitor (Neq0682) of cysteine protease. Then, we used advanced experimental and computational techniques to explore details of the inhibition mechanism of cruzain by reversible and irreversible inhibitors. The isothermal titration calorimetry (ITC) analysis shows that inhibition of cruzain by an irreversible inhibitor is thermodynamically more favorable than by a reversible one. The hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) simulations were used to explore the mechanism of the reaction inhibition of cruzain by K777 and Neq0682. The calculated free energy profiles show that the Cys25 nucleophilic attack and His162 proton transfer occur in a single step for a reversible inhibitor and two steps for an irreversible covalent inhibitor. The hybrid QM/MM calculated free energies for the inhibition reaction correspond to -26.7 and -5.9 kcal mol^-1 for K777 and Neq0682 at the MP2/MM level, respectively. These results indicate that the ΔG of the reaction is very negative for the process involving K777, consequently, the covalent adduct cannot revert to a noncovalent protein-ligand complex, and its binding tends to be irreversible. Overall, the present study provides insights into a covalent inhibition mechanism of cysteine proteases.

Dynamical Behavior of β-Lactamases and Penicillin- Binding Proteins in Different Functional States and Its Potential Role in Evolution

β-Lactamases are enzymes produced by bacteria to hydrolyze β-lactam-based antibiotics, and pose serious threat to public health through related antibiotic resistance. Class A β-lactamases are structurally and functionally related to penicillin-binding proteins (PBPs). Despite the extensive studies of the structures, catalytic mechanisms and dynamics of both β-lactamases and PBPs, the potentially different dynamical behaviors of these proteins in different functional states still remain elusive in general. In this study, four evolutionarily related proteins, including TEM-1 and TOHO-1 as class A β-lactamases, PBP-A and DD-transpeptidase as two PBPs, are subjected to molecular dynamics simulations and various analyses to characterize their dynamical behaviors in different functional states. Penicillin G and its ring opening product serve as common ligands for these four proteins of interest. The dynamic analyses of overall structures, the active sites with penicillin G, and three catalytically important residues commonly shared by all four proteins reveal unexpected cross similarities between Class A β-lactamases and PBPs. These findings shed light on both the hidden relations among dynamical behaviors of these proteins and the functional and evolutionary relations among class A β-lactamases and PBPs.

A Mechanistic Study of a Potent and Selective Epidermal Growth Factor Receptor Inhibitor against the L858R/T790M Resistance Mutation

Covalent targeting is a promising strategy for increasing the potency and selectivity of potential drug candidates. This therapeutic approach was recently reported for the epidermal growth factor receptor (EGFR), wherein a covalent binder, 20g [N-(3-{7-[2-methoxy-4-(4-methylpiperazin-1-yl)phenylamino]-3,4-dihydro-3-isopropyl-2,4-dioxopyrimido[4,5-d]pyrimidin-1(2H)-yl}phenyl)acrylamide], demonstrated significant selectivity and inhibitory activity toward the EGFR L858R/T790M double mutant (EGFR^DM) relative to the EGFR wild-type form (EGFR^WT). The enhanced therapeutic potency of 20g against EGFR^DM is 263 times greater than that against EGFR^WT, which necessitates a rational explanation for the underlying selective and inhibitory mechanisms. In this work, we investigate the differential binding modes of 20g with EGFR^WT and EGFR^DM using molecular dynamics simulations coupled with free energy calculations and further identify key residues involved in the selective targeting, binding, and inhibitory mechanisms mediated by 20g. We find that systematic orientational and conformational changes in the α-loop, p-loop, active loop, and αC-helix are responsible for the disparate binding mechanisms and inhibitory prowess of 20g with respect to EGFR^WT and EGFR^DM. The calculated binding free energies show good correlation with the experimental biological activity. The total binding free energy difference between EGFR^WT-20g and EGFR^DM-20g is -11.47 kcal/mol, implying that 20g binds more strongly to EGFR^DM. This enhanced binding affinity of 20g for EGFR^DM is a result of a large increase in the van der Waals and electrostatic interactions with three critical residues (Met790, Gln791, and Met793) that are chiefly responsible for the high-affinity interactions mediated by 20g with EGFR^DM relative to EGFR^WT.

β-Lactamases and β-Lactamase Inhibitors in the 21st Century

The β-lactams retain a central place in the antibacterial armamentarium. In Gram-negative bacteria, β-lactamase enzymes that hydrolyze the amide bond of the four-membered β-lactam ring are the primary resistance mechanism, with multiple enzymes disseminating on mobile genetic elements across opportunistic pathogens such as Enterobacteriaceae (e.g., Escherichia coli) and non-fermenting organisms (e.g., Pseudomonas aeruginosa). β-Lactamases divide into four classes; the active-site serine β-lactamases (classes A, C and D) and the zinc-dependent or metallo-β-lactamases (MBLs; class B). Here we review recent advances in mechanistic understanding of each class, focusing upon how growing numbers of crystal structures, in particular for β-lactam complexes, and methods such as neutron diffraction and molecular simulations, have improved understanding of the biochemistry of β-lactam breakdown. A second focus is β-lactamase interactions with carbapenems, as carbapenem-resistant bacteria are of grave clinical concern and carbapenem-hydrolyzing enzymes such as KPC (class A) NDM (class B) and OXA-48 (class D) are proliferating worldwide. An overview is provided of the changing landscape of β-lactamase inhibitors, exemplified by the introduction to the clinic of combinations of β-lactams with diazabicyclooctanone and cyclic boronate serine β-lactamase inhibitors, and of progress and strategies toward clinically useful MBL inhibitors. Despite the long history of β-lactamase research, we contend that issues including continuing unresolved questions around mechanism; opportunities afforded by new technologies such as serial femtosecond crystallography; the need for new inhibitors, particularly for MBLs; the likely impact of new β-lactam:inhibitor combinations and the continuing clinical importance of β-lactams mean that this remains a rewarding research area.

Luciferin Regeneration in Firefly Bioluminescence via Proton-Transfer-Facilitated Hydrolysis, Condensation and Chiral Inversion

Firefly bioluminescence is produced via luciferin enzymatic reactions in luciferase. Luciferin has to be unceasingly replenished to maintain bioluminescence. How is the luciferin reproduced after it has been exhausted? In the early 1970s, Okada proposed the hypothesis that the oxyluciferin produced by the previous bioluminescent reaction could be converted into new luciferin for the next bioluminescent reaction. To some extent, this hypothesis was evidenced by several detected intermediates. However, the detailed process and mechanism of luciferin regeneration remained largely unknown. For the first time, we investigated the entire process of luciferin regeneration in firefly bioluminescence by density functional theory calculations. This theoretical study suggests that luciferin regeneration consists of three sequential steps: the oxyluciferin produced from the last bioluminescent reaction generates 2-cyano-6-hydroxybenzothiazole (CHBT) in the luciferin regenerating enzyme (LRE) via a hydrolysis reaction; CHBT combines with L-cysteine in vivo to form L-luciferin via a condensation reaction; and L-luciferin inverts into D-luciferin in luciferase and thioesterase. The presently proposed mechanism not only supports the sporadic evidence from previous experiments but also clearly describes the complete process of luciferin regeneration. This work is of great significance for understanding the long-term flashing of fireflies without an in vitro energy supply.

Influence of substrates and inhibitors on the structure of Klebsiella pneumoniae carbapenemase-2

The hydrolysis of last resort carbapenem antibiotics by Klebsiella pneumoniae carbapenemase-2 (KPC-2) presents a significant danger to global health. Combined with horizontal gene transfer, the emergence KPC-2 threatens to quickly expand carbapenemase activity to ever increasing numbers of pathogens. Our understanding of KPC-2 has greatly increased over the past decade thanks, in great part, to 20 crystal structures solved by groups around the world. These include apo KPC-2 structures, along with structures featuring a library of 10 different inhibitors representing diverse structural and functional classes. Herein we focus on cataloging the available KPC-2 structures and presenting a discussion of key aspects of each structure and important relationships between structures. Although the available structures do not provide information on dynamic motions with KPC-2, and the family of structures indicates small conformational changes across a wide array of bound inhibitors, substrates, and products, the structures provide a strong foundation for additional studies in the coming years to discover new KPC-2 inhibitors.

Impact statement

The work herein is important to the field as it provides a clear and succinct accounting of available KPC-2 structures. The work advances the field by collecting and analyzing differences and similarities across the available structures. This work features new analyses and interpretations of the existing structures which will impact the field in a positive way by making structural insights more widely available among the beta-lactamase community.

Semi-Empirical Born-Oppenheimer Molecular Dynamics (SEBOMD) within the Amber Biomolecular Package

Semi-empirical quantum methods from the neglect of differential diatomic overlap (NDDO) family such as MNDO, AM1, or PM3 are fast albeit approximate quantum methods. By combining them with linear scaling methods like the divide & conquer (D&C) method, it is possible to quickly evaluate the energy of systems containing hundreds to thousands of atoms. We here present our implementation in the Amber biomolecular package of a SEBOMD module that provides a way to run semi-empirical Born-Oppenheimer molecular dynamics. At each step of a SEBOMD, a fully converged self-consistent field (SCF) calculation is performed to obtain the semiempirical quantum potential energy of a molecular system encaged or not in periodic boundary conditions. We describe the implementation and the features of our SEBOMD implementation. We show the requirements to conserve the total energy in NVE simulations, and how to accelerate SCF convergence through density matrix extrapolation. Specific ways of handling periodic boundary conditions using mechanical embedding or electrostatic embedding through a tailored quantum Ewald summation is developed. The parallel performance of SEBOMD simulations using the D&C scheme are presented for liquid water systems of various sizes, and a comparison between the traditional full diagonalization scheme and the D&C approach for the reproduction of the structure of liquid water illustrates the potentiality of SEBOMD to simulate molecular systems containing several hundreds of atoms for hundreds of picoseconds with a quantum mechanical potential in a reasonable amount of CPU time.

Defining the architecture of KPC-2 Carbapenemase: identifying allosteric networks to fight antibiotics resistance

The rise of multi-drug resistance in bacterial pathogens is one of the grand challenges facing medical science. A major concern is the speed of development of β-lactamase-mediated resistance in Gram-negative species, thus putting at risk the efficacy of the most recently approved antibiotics and inhibitors, including carbapenems and avibactam, respectively. New strategies to overcome resistance are urgently required, which will ultimately be facilitated by a deeper understanding of the mechanisms that regulate the function of β-lactamases such as the Klebsiella Pneumoniae carbapenemases (KPCs). Using enhanced sampling computational methods together with site-directed mutagenesis, we report the identification of two “hydrophobic networks” in the KPC-2 enzyme, the integrity of which has been found to be essential for protein stability and corresponding resistance. Present throughout the structure, these networks are responsible for the structural integrity and allosteric signaling. Disruption of the networks leads to a loss of the KPC-2 mediated resistance phenotype, resulting in restored susceptibility to different classes of β-lactam antibiotics including carbapenems and cephalosporins. The ”hydrophobic networks” were found to be highly conserved among class-A β-lactamases, which implies their suitability for exploitation as a potential target for therapeutic intervention.

Multiscale Simulations of Clavulanate Inhibition Identify the Reactive Complex in Class A β-Lactamases and Predict the Efficiency of Inhibition

Clavulanate is used as an effective drug in combination with β-lactam antibiotics to treat infections of some antibiotic resistant bacteria. Here, we perform combined quantum mechanics/molecular mechanics simulations of several covalent complexes of clavulanate with class A β-lactamases KPC-2 and TEM-1. Simulations of the deacylation reactions identify the decarboxylated trans-enamine complex as being responsible for inhibition. Further, the obtained free energy barriers discriminate clinically relevant inhibition (TEM-1) from less effective inhibition (KPC-2).

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

QM/MM simulations identify the determinants of catalytic activity differences between type II dehydroquinase enzymes

Type II dehydroquinase enzymes (DHQ2), recognized targets for antibiotic drug discovery, show significantly different activities dependent on the species: DHQ2 from Mycobacterium tuberculosis (MtDHQ2) and Helicobacter pylori (HpDHQ2) show a 50-fold difference in catalytic efficiency. Revealing the determinants of this activity difference is important for our understanding of biological catalysis and further offers the potential to contribute to tailoring specificity in drug design. Molecular dynamics simulations using a quantum mechanics/molecular mechanics potential, with correlated ab initio single point corrections, identify and quantify the subtle determinants of the experimentally observed difference in efficiency. The rate-determining step involves the formation of an enolate intermediate: more efficient stabilization of the enolate and transition state of the key step in MtDHQ2, mainly by the essential residues Tyr24 and Arg19, makes it more efficient than HpDHQ2. Further, a water molecule, which is absent in MtDHQ2 but involved in generation of the catalytic Tyr22 tyrosinate in HpDHQ2, was found to destabilize both the transition state and the enolate intermediate. The quantification of the contribution of key residues and water molecules in the rate-determining step of the mechanism also leads to improved understanding of higher potencies and specificity of known inhibitors, which should aid ongoing inhibitor design.

Mechanisms of proton relay and product release by Class A β-lactamase at ultrahigh resolution

The β-lactam antibiotics inhibit penicillin-binding proteins (PBPs) by forming a stable, covalent, acyl-enzyme complex. During the evolution from PBPs to Class A β-lactamases, the β-lactamases acquired Glu166 to activate a catalytic water and cleave the acyl-enzyme bond. Here we present three product complex crystal structures of CTX-M-14 Class A β-lactamase with a ruthenocene-conjugated penicillin-a 0.85 Å resolution structure of E166A mutant complexed with the penilloate product, a 1.30 Å resolution complex structure of the same mutant with the penicilloate product, and a 1.18 Å resolution complex structure of S70G mutant with a penicilloate product epimer-shedding light on the catalytic mechanisms and product inhibition of PBPs and Class A β-lactamases. The E166A-penilloate complex captured the hydrogen bonding network following the protonation of the leaving group and, for the first time, unambiguously show that the ring nitrogen donates a proton to Ser130, which in turn donates a proton to Lys73. These observations indicate that in the absence of Glu166, the equivalent lysine would be neutral in PBPs and therefore capable of serving as the general base to activate the catalytic serine. Together with previous results, this structure suggests a common proton relay network shared by Class A β-lactamases and PBPs, from the catalytic serine to the lysine, and ultimately to the ring nitrogen. Additionally, the E166A-penicilloate complex reveals previously unseen conformational changes of key catalytic residues during the release of the product, and is the first structure to capture the hydrolyzed product in the presence of an unmutated catalytic serine.

DATABASE

Structural data are available in the PDB database under the accession numbers 5TOP, 5TOY, and 5VLE.

Predicting allosteric mutants that increase activity of a major antibiotic resistance enzyme

Allosteric mutations increasing k_cat in a beta lactamase act by changing conformational ensembles of active-site residues identified by machine learning.

The CTX-M family of beta lactamases mediate broad-spectrum antibiotic resistance and are present in the majority of drug-resistant Gram-negative bacterial infections worldwide. Allosteric mutations that increase catalytic rates of these drug resistance enzymes have been identified in clinical isolates but are challenging to predict prospectively. We have used molecular dynamics simulations to predict allosteric mutants increasing CTX-M9 drug resistance, experimentally testing top mutants using multiple antibiotics. Purified enzymes show an increase in catalytic rate and efficiency, while mutant crystal structures show no detectable changes from wild-type CTX-M9. We hypothesize that increased drug resistance results from changes in the conformational ensemble of an acyl intermediate in hydrolysis. Machine-learning analyses on the three top mutants identify changes to the binding-pocket conformational ensemble by which these allosteric mutations transmit their effect. These findings show how molecular simulation can predict how allosteric mutations alter active-site conformational equilibria to increase catalytic rates and thus resistance against common clinically used antibiotics.

Quantum mechanics/molecular mechanics modeling of covalent addition between EGFR-cysteine 797 and N-(4-anilinoquinazolin-6-yl) acrylamide

Irreversible epidermal growth factor receptor (EGFR) inhibitors can circumvent resistance to first-generation ATP-competitive inhibitors in the treatment of nonsmall-cell lung cancer. They covalently bind a noncatalytic cysteine (Cys797) at the surface of EGFR active site by an acrylamide warhead. Herein, we used a hybrid quantum mechanics/molecular mechanics (QM/MM) potential in combination with umbrella sampling in the path-collective variable space to investigate the mechanism of alkylation of Cys797 by the prototypical covalent inhibitor N-(4-anilinoquinazolin-6-yl) acrylamide. Calculations show that Cys797 reacts with the acrylamide group of the inhibitor through a direct addition mechanism, with Asp800 acting as a general base/general acid in distinct steps of the reaction. The obtained reaction free energy is negative (ΔA = -12 kcal/mol) consistent with the spontaneous and irreversible alkylation of Cys797 by N-(4-anilinoquinazolin-6-yl) acrylamide. Our calculations identify desolvation of Cys797 thiolate anion as a key step of the alkylation process, indicating that changes in the intrinsic reactivity of the acrylamide would have only a minor impact on the inhibitor potency.

Evidence for substrate binding-induced zwitterion formation in the catalytic Cys-His dyad of the SARS-CoV main protease

The coronavirus main protease (M(pro)) represents an attractive drug target for antiviral therapy of coronavirus (CoV) infections, including severe acute respiratory syndrome (SARS). The SARS-CoV M(pro) and related CoV proteases have several distinct features, such as an uncharged Cys-His catalytic dyad embedded in a chymotrypsin-like protease fold, that clearly separate these enzymes from archetypical cysteine proteases. To further characterize the catalytic system of CoV main proteases and to obtain information about improved inhibitors, we performed comprehensive simulations of the proton-transfer reactions in the SARS-CoV M(pro) active site that lead to the Cys(-)/His(+) zwitterionic state required for efficient proteolytic activity. Our simulations, comprising the free enzyme as well as substrate-enzyme and inhibitor-enzyme complexes, lead us to predict that zwitterion formation is fostered by substrate binding but not inhibitor binding. This indicates that M(pro) employs a substrate-induced catalytic mechanism that further enhances its substrate specificity. Our computational data are in line with available experimental results, such as X-ray geometries, measured pKa values, mutagenesis experiments, and the measured differences between the kinetic parameters of substrates and inhibitors. The data also provide an atomistic picture of the formerly postulated electrostatic trigger involved in SARS-CoV M(pro) activity. Finally, they provide information on how a specific microenvironment may finely tune the activity of M(pro) toward specific viral protein substrates, which is known to be required for efficient viral replication. Our simulations also indicate that the low inhibition potencies of known covalently interacting inhibitors may, at least in part, be attributed to insufficient fostering of the proton-transfer reaction. These findings suggest ways to achieve improved inhibitors.

Can molecular dynamics and QM/MM solve the penicillin binding protein protonation puzzle?

Benzylpenicillin, a member of the β-lactam antibiotic class, has been widely used to combat bacterial infections since 1947. The general mechanism is well-known: a serine protease enzyme (i.e., DD-peptidase) forms a long lasting intermediate with the lactam ring of the antibiotic known as acylation, effectively preventing biosynthesis of the bacterial cell wall. Despite this overall mechanistic understanding, many details of binding and catalysis are unclear. Specifically, there is ongoing debate about active site protonation states and the role of general acids/bases in the reaction. Herein, a unique combination of MD simulations, QM/MM minimizations, and QM/MM orbital analyses is combined with systematic variation of active site residue protonation states. Critical interactions that maximize the stability of the bound inhibitor are examined and used as metrics. This approach was validated by examining cefoxitin interactions in the CTX-M β-lactamase from E. coli and compared to an ultra high-resolution (0.88 Å) crystal structure. Upon confirming the approach used, an investigation of the preacylated Streptomyces R61 active site with bound benzylpenicillin was performed, varying the protonation states of His298 and Lys65. We concluded that protonated His298 and deprotonated Lys65 are most likely to exist in the R61 active site.

Nylon-Oligomer Hydrolase Promoting Cleavage Reactions in Unnatural Amide Compounds

The active site of 6-aminohexanoate-dimer hydrolase, a nylon-6 byproduct-degrading enzyme with a β-lactamase fold, possesses a Ser112/Lys115/Tyr215 catalytic triad similar to the one of penicillin-recognizing family of serine-reactive hydrolases but includes a unique Tyr170 residue. By using a reactive quantum mechanics/molecular mechanics (QM/MM) approach, we work out its catalytic mechanism and related functional/structural specificities. At variance with other peptidases, we show that the involvement of Tyr170 in the enzyme-substrate interactions is responsible for a structural variation in the substrate-binding state. The acylation via a tetrahedral intermediate is the rate-limiting step, with a free-energy barrier of ∼21 kcal/mol, driven by the catalytic triad Ser112, Lys115, and Tyr215, acting as a nucleophile, general base, and general acid, respectively. The functional interaction of Tyr170 with this triad leads to an efficient disruption of the tetrahedral intermediate, promoting a conformational change of the substrate favorable for proton donation from the general acid.

Catalytic mechanism of histone acetyltransferase p300: from the proton transfer to acetylation reaction

The transcriptional coactivator and histone acetyltransferase (HAT) p300 acetylates the four core histones and other transcription factors to regulate a plethora of fundamental biological processes including cell growth, development, oncogenesis and apoptosis. Recent structural and biochemical studies on the p300 HAT domain revealed a Theorell-Chance, or "hit-and-run", catalytic mechanism. Nonetheless, the chemical mechanism of the entire reaction process including the proton transfer (PT) scheme and consequent acetylation reaction route remains unclear. In this study, a combined computational strategy consisting of molecular modeling, molecular dynamic (MD) simulation, and quantum mechanics/molecular mechanics (QM/MM) simulation was applied to elucidate these important issues. An initial p300/H3/Ac-CoA complex structure was modeled and optimized using a 100 ns MD simulation. Residues that play important roles in substrate binding and the acetylation reaction were comprehensively investigated. For the first time, these studies reveal a plausible PT scheme consisting of Y1394, D1507, and a conserved crystallographic water molecule, with all components of the scheme being stable during the MD simulation and the energy barrier low for PT to occur. The two-dimensional potential energy surface for the nucleophilic attack process was also calculated. The comparison of potential energies for two possible elimination half-reaction mechanisms revealed that Y1467 reprotonates the coenzyme-A leaving group to form product. This study provides new insights into the detailed catalytic mechanism of p300 and has important implications for the discovery of novel small molecule regulators for p300.