Catalytic mechanism and metal specificity of bacterial peptide deformylase: a density functional theory QM/MM study

Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine

A recently designed and discovered cocaine hydrolase (CocH), engineered from human butyrylcholinesterase (BChE), has been proven promising as a novel enzyme therapy for treatment of cocaine overdose and addiction because it is highly efficient in catalyzing hydrolysis of naturally occurring (-)-cocaine. It has been known that the CocH also has a high catalytic efficiency against (+)-cocaine, a synthetic enantiomer of cocaine. Reaction pathway and the corresponding free energy profile for the CocH-catalyzed hydrolysis of (+)-cocaine have been determined, in the present study, by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Acordingt to the QM/MM-FE results, the catalytic hydrolysis process is initiated by the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester via hydroxyl oxygen of S198 side chain, and the second reaction step (i.e. dissociation of benzoyl ester) is rate-determining. This finding for CocH-catalyzed hydrolysis of (+)-cocaine is remarkably different from that for the (+)-cocaine hydrolysis catalyzed by bacterial cocaine esterase in which the first reaction step of the deacylation is associated with the highest free energy barrier (~17.9 kcal/mol). The overall free energy barrier (~16.0 kcal/mol) calculated for the acylation stage of CocH-catalyzed hydrolysis of (+)-cocaine is in good agreement with the experimental free energy barrier of ~14.5 kcal/mol derivated from the experimental kinetic data.

Thiol versus hydroxamate as zinc binding group in HDAC inhibition: An ab initio QM/MM molecular dynamics study

Zinc-dependent histone deacetylases (HDACs) play a critical role in transcriptional repression and gene silencing, and are among the most attractive targets for the development of new therapeutics against cancer and various other diseases. Two HDAC inhibitors have been approved by FDA as anti-cancer drugs: one is SAHA whose hydroxamate is directly bound to zinc, the other is FK228 whose active form may use thiol as the zinc binding group. In spite of extensive studies, it remains to be ambiguous regarding how thiol and hydroxamate are bound to the zinc active site of HDACs. In this work, our computational approaches center on Born-Oppenheimer ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics with umbrella sampling, which allow for modeling of the zinc active site with reasonable accuracy while properly including dynamics and effects of protein environment. Meanwhile, an improved short-long effective function (SLEF2) to describe non-bonded interactions between zinc and other atoms has been employed in initial MM equilibrations. Our ab initio QM/MM MD simulations have confirmed that hydroxamate is neutral when it is bound to HDAC8, and found that thiol is deprotonated when directly bound to zinc in the HDAC active site. By comparing thiol and hydroxamate, our results elucidated the differences in their binding environment in the HDAC active sites, and emphasized the importance of the linker design to achieve more specific binding toward class IIa HDACs.

Why does the G117H mutation considerably improve the activity of human butyrylcholinesterase against sarin? Insights from quantum mechanical/molecular mechanical free energy calculations

Human butyrylcholinesterase (BChE) is recognized as the most promising bioscavenger for organophosphorus (OP) warfare nerve agents. The G117H mutant of human BChE has been identified as a potential catalytic bioscavenger with a remarkably improved activity against OP nerve agents such as sarin, but it still does not satisfy the clinical use. For further design of the higher-activity mutants against OP nerve agents, it is essential to understand how the G117H mutation improves the activity. The reaction mechanisms and the free energy profiles for spontaneous reactivation of wild-type BChE and its G117H mutant phosphorylated by sarin have been explored, in this study, by performing first-principles quantum mechanical/molecular mechanical free energy calculations, and the remarkable role of the G117H mutation on the activity has been elucidated. For both the wild-type and G117H mutant enzymes, H438 acts as a general base to initiate the spontaneous reactivation that consists of two reaction steps: the nucleophilic attack at the phosphorus by a water molecule and decomposition of the pentacoordinated phosphorus intermediate. The calculated overall free energy barriers, i.e., 30.2 and 23.9 kcal/mol for the wild type and G117H mutant, respectively, are in good agreement with available experimental kinetic data. On the basis of the calculated results, the mutated residue (H117 in the G117H mutant) cannot initiate the spontaneous reactivation as a general base. Instead, it skews the oxyanion hole and makes the phosphorus more open to the nucleophilic water molecule, resulting in a remarkable change in the rate-determining step and significantly improved catalytic activity of human BChE.

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

Reaction mechanism for cocaine esterase-catalyzed hydrolyses of (+)- and (-)-cocaine: unexpected common rate-determining step

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to examine the catalytic mechanism for cocaine esterase (CocE)-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. It has been shown that the acylation of (+)-cocaine consists of nucleophilic attack of the hydroxyl group of Ser117 on the carbonyl carbon of (+)-cocaine benzoyl ester and the dissociation of (+)-cocaine benzoyl ester. The first reaction step of deacylation of (+)-cocaine, which is identical to that of (-)-cocaine, is rate-determining, indicating that CocE-catalyzed hydrolyses of (+)- and (-)-cocaine have a common rate-determining step. The computational results predict that the catalytic rate constant of CocE against (+)-cocaine should be the same as that of CocE against (-)-cocaine, in contrast with the remarkable difference between human butyrylcholinesterase-catalyzed hydrolyses of (+)- and (-)-cocaine. The prediction has been confirmed by experimental kinetic analysis on CocE-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. The determined common rate-determining step indicates that rational design of a high-activity mutant of CocE should be focused on the first reaction step of the deacylation. Furthermore, the obtained mechanistic insights into the detailed differences in the acylation between the (+)- and (-)-cocaine hydrolyses provide indirect clues for rational design of amino acid mutations that could more favorably stabilize the rate-determining transition state in the deacylation and, thus, improve the catalytic activity of CocE. This study provides a valuable mechanistic base for rational design of an improved esterase for therapeutic treatment of cocaine abuse.

Zinc chelation with hydroxamate in histone deacetylases modulated by water access to the linker binding channel

It is of significant biological interest and medical importance to develop class- and isoform-selective histone deacetylase (HDAC) modulators. The impact of the linker component on HDAC inhibition specificity has been revealed but is not understood. Using Born-Oppenheimer ab initio QM/MM MD simulations, a state-of-the-art approach to simulating metallo-enzymes, we have found that the hydroxamic acid remains to be protonated upon its binding to HDAC8, and thus disproved the mechanistic hypothesis that the distinct zinc-hydroxamate chelation modes between two HDAC subclasses come from different protonation states of the hydroxamic acid. Instead, our simulations suggest a novel mechanism in which the chelation mode of hydroxamate with the zinc ion in HDACs is modulated by water access to the linker binding channel. This new insight into the interplay between the linker binding and the zinc chelation emphasizes its importance and gives guidance regarding linker design for the development of new class-IIa-specific HDAC inhibitors.

A Transferable Non-bonded Pairwise Force Field to Model Zinc Interactions in Metalloproteins

Herein we introduce a novel practical strategy to overcome the well-known challenge of modeling the divalent zinc cation in metalloproteins. The main idea is to design short-long effective functions (SLEF) to describe charge interactions between the zinc ion and all other atoms. This SLEF approach has the following desired features: (1). It is pairwise, additive and compatible with widely used atomic pair-wise force fields for modeling biomolecules; (2). It only changes interactions between the zinc ion and other atoms, and does not affect force field parameters that model other interactions in the system; (3). It is a non-bonded model that is inherently capable to describe different zinc ligands and coordination modes. By optimizing two SLEF parameters as well as zinc vdW parameters through force matching based on Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, we have successfully developed the first SLEF force field (SLEF1) to describe zinc interactions. Extensive molecular dynamics simulations of seven zinc enzyme systems with different coordination ligands and distinct chelation modes (4-,5-,6-fold), including the binuclear zinc active site, yielded zinc coordination numbers and binding distances in good agreement with the corresponding crystal structures as well as ab initio QM/MM MD results. This not only demonstrates the transferability and adequacy of the new SLEF1 force field in describing a variety of zinc proteins, but also indicates that this novel SLEF approach is a promising direction to explore for improving force field description of metal ion interactions.

A proton-shuttle reaction mechanism for histone deacetylase 8 and the catalytic role of metal ions

Zinc-dependent histone deacetylase 8 (HDAC8) catalyzes the removal of acetyl moieties from histone tails, and is critically involved in regulating chromatin structure and gene expression. The detailed knowledge of its catalytic process is of high importance since it has been established as a most promising target for the development of new antitumor drugs. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations and umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have provided further evidence against the originally proposed general acid-base catalytic pair mechanism for Zinc-dependent histone deacetylases. Instead, our results indicated that HDAC8 employs a proton-shuttle catalytic mechanism, in which a neutral His143 first serves as the general base to accept a proton from the zinc-bound water molecule in the initial rate-determining nucleophilic attack step, and then shuttles it to the amide nitrogen atom to facilitate the cleavage of the amide bond. During the deacetylation process, the Zn(2+) ion changes its coordination mode and plays multiple catalytic roles. For the K(+) ion, which is located about 7 A from the catalytic Zn(2+) ion and conserved in class I and II HDACs, our simulations indicated that its removal would lead to the different transition state structure and a higher free energy reaction barrier for the rate-determining step. It is found that the existence of this conserved K(+) ion would enhance the substrate binding, increase the basicity of His143, strengthen the catalytic role of zinc ion, and improve the transition state stabilization by the enzyme environment.

Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase

The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.

Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study

The different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly due to the well known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio QM/MM molecular dynamics simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in Thermolysin (TLN) and Histone Deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's time-scale. For TLN, the coordination of the carboxylate group of Glu166 to Zinc is found to continuously change between monodentate and bidentate manner dynamically; while for HDAC8, the flexibility mainly comes from the coordination to a non-amino-acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of Zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.

Fundamental reaction mechanism and free energy profile for (-)-cocaine hydrolysis catalyzed by cocaine esterase

The fundamental reaction mechanism of cocaine esterase (CocE)-catalyzed hydrolysis of (-)-cocaine and the corresponding free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, the entire hydrolysis reaction consists of four reaction steps, including the nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by the hydroxyl group of Ser117, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally dissociation between the (-)-cocaine benzoyl group and Ser117 of CocE. The third reaction step involving the nucleophilic attack of a water molecule was found to be rate-determining, which is remarkably different from (-)-cocaine hydrolysis catalyzed by wild-type butyrylcholinesterase (BChE; where the formation of the prereactive BChE-(-)-cocaine complex is rate-determining) or its mutants containing Tyr332Gly or Tyr332Ala mutation (where the first chemical reaction step is rate-determining). Besides, the role of Asp259 in the catalytic triad of CocE does not follow the general concept of the "charge-relay system" for all serine esterases. The free energy barrier calculated for the rate-determining step of CocE-catalyzed hydrolysis of (-)-cocaine is 17.9 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 16.2 kcal/mol. In the present study, where many sodium ions are present, the effects of counterions are found to be significant in determining the free energy barrier. The finding of the significant effects of counterions on the free energy barrier may also be valuable in guiding future mechanistic studies on other charged enzymes.

Peptide deformylase--a promising therapeutic target for tuberculosis and antibacterial drug discovery

BACKGROUND

Tuberculosis (TB) remains the most important infectious disease causing morbidity and death, due to the human pathogen Mycobacterium tuberculosis. The emergence of multi-drug-resistant and extensively-drug-resistant forms of TB have resulted in an increase in the number of TB cases. Thus, there is an urgent need to identify new drugs with novel targets to ensure future therapeutic success. Studies have indicated that peptide deformylase is an interesting potential candidate for discovering antimicrobial agents.

OBJECTIVE

To explore the role of peptide deformylase, a highly conserved metalloprotease and an essential enzyme in bacterial life cycle, as a target for antibacterial as well as antimycobacterial drug development.

METHODS

This review is based on recent published literature and online resources related to peptide deformylase inhibitors and their antibacterial potential.

RESULTS/CONCLUSION

Peptide deformylase is an emerging therapeutic target for the treatment of tuberculosis and peptide deformylase inhibitors can act as potential future antibacterial agents.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Highly dissociative and concerted mechanism for the nicotinamide cleavage reaction in Sir2Tm enzyme suggested by ab initio QM/MM molecular dynamics simulations

Sir2 enzymes catalyze the NAD+-dependent protein deacetylation and play critical roles in epigenetics, cell death, and lifespan regulation. In spite of a current flurry of experimental studies, the catalytic mechanism for this unique and important class of enzymes remains elusive. Employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized the initial step of the Sir2Tm-catalyzed reaction, which is also the most controversial portion of its mechanism. Our results indicate that the nicotinamide cleavage reaction employs a highly dissociative and concerted displacement mechanism: the cleavage of the glycosidic bond is facilitated by the nucleophilic participation of the acetyl-lysine, and the dissociative transition state has a significant oxocarbenium ion character. During this step of the reaction, the Sir2Tm enzyme strongly stabilizes the covalent O-alkylamidate intermediate whereas its effect on the transition state is quite minimal. In addition, functional roles of key residues and motifs have been elucidated. This work further demonstrates the feasibility and applicability of the state-of-the-art ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods

Combined quantum mechanics/molecular mechanics (QM/MM) methods provide an accurate and efficient energetic description of complex chemical and biological systems, leading to significant advances in the understanding of chemical reactions in solution and in enzymes. Here we review progress in QM/MM methodology and applications, focusing on ab initio QM-based approaches. Ab initio QM/MM methods capitalize on the accuracy and reliability of the associated quantum-mechanical approaches, however, at a much higher computational cost compared with semiempirical quantum-mechanical approaches. Thus reaction-path and activation free-energy calculations based on ab initio QM/MM methods encounter unique challenges in simulation timescales and phase-space sampling. This review features recent developments overcoming these challenges and enabling accurate free-energy determination for reaction processes in solution and in enzymes, along with applications.

Design-atom approach for the quantum mechanical/molecular mechanical covalent boundary: a design-carbon atom with five valence electrons

A critical issue underlying the accuracy and applicability of the combined quantum mechanical/molecular mechanical (QM/MM) methods is how to describe the QM/MM boundary across covalent bonds. Inspired by the ab initio pseudopotential theory, here we introduce a novel design atom approach for a more fundamental and transparent treatment of this QM/MM covalent boundary problem. The main idea is to replace the boundary atom of the active part with a design atom, which has a different number of valence electrons but very similar atomic properties. By modifying the Troullier-Martins scheme, which has been widely employed to construct norm-conserving pseudopotentials for density functional calculations, we have successfully developed a design-carbon atom with five valence electrons. Tests on a series of molecules yield very good structural and energetic results and indicate its transferability in describing a variety of chemical bonds, including double and triple bonds.