Ab initio QM/MM free-energy studies of arginine deiminase catalysis: the protonation state of the Cys nucleophile

Thermostability and Specific-Activity Enhancement of an Arginine Deiminase from Enterococcus faecalis SK23.001 via Semirational Design for l-Citrulline Production

l-Citrulline is a nonessential amino acid with a variety of physiological functions and can be enzymatically produced by arginine deiminase (ADI, EC 3.5.3.6). The enzymatic-production approach is of immense interest because of its mild conditions, high yield, low cost, and environmental benignity. However, the major hindrances of l-citrulline industrialization are the poor thermostability and enzyme activity of ADI. Hence, in this work, directed evolution and site-directed mutagenesis aided with in silico screening, including the use of b-factor values and HoTMuSiC, were applied to a previously identified ADI from Enterococcus faecalis SK23.001 ( EfADI), and a triple-site variant R15K-F269Y-G292P was obtained. The triple-site variant displays a 2.5-fold higher specific enzyme activity (333 U mg^-1), a lower K_m value of 6.4 mM, and a 6.1-fold longer half-life ( t_1/2,45°C = 86.7 min) than wild-type EfADI. This work provides a protein-engineering strategy to improve enzyme activity and thermostability, which might be transferrable to other ADIs and enzymes.

QM/MM and MM MD Simulations on the Pyrimidine-Specific Nucleoside Hydrolase: A Comprehensive Understanding of Enzymatic Hydrolysis of Uridine

The pyrimidine-specific nucleoside hydrolase Yeik (CU-NH) from Escherichia coli cleaves the N-glycosidic bond of uridine and cytidine with a 10²-10⁴-fold faster rate than that of purine nucleoside substrates, such as inosine. Such a remarkable substrate specificity and the plausible hydrolytic mechanisms of uridine have been explored by using QM/MM and MM MD simulations. The present calculations show that the relatively stronger hydrogen-bond interactions between uridine and the active-site residues Gln227 and Tyr231 in CU-NH play an important role in enhancing the substrate binding and thus promoting the N-glycosidic bond cleavage, in comparison with inosine. The estimated energy barrier of 30 kcal/mol for the hydrolysis of inosine is much higher than 22 kcal/mol for uridine. Extensive MM MD simulations on the transportation of substrates to the active site of CU-NH indicate that the uridine binding is exothermic by ∼23 kcal/mol, more remarkable than inosine (∼12 kcal/mol). All of these arise from the noncovalent interactions between the substrate and the active site featured in CU-NH, which account for the substrate specificity. Quite differing from other nucleoside hydrolases, here the enzymatic N-glycosidic bond cleavage of uridine is less influenced by its protonation.

Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase

QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.

Theoretical insights into the protonation states of active site cysteine and citrullination mechanism of Porphyromonas gingivalis peptidylarginine deiminase

Porphyromonas gingivalis peptidylarginine deiminase (PPAD) catalyzes the citrullination of peptidylarginine, which plays a critical role in the rheumatoid arthritis (RA) and gene regulation. For a better understanding of citrullination mechanism of PPAD, it is required to establish the protonation states of active site cysteine, which is still a controversial issue for the members of guanidino-group-modifying enzyme superfamily. In this work, we first explored the transformation between the two states: State N (both C351 and H236 are neutral) and State I (both residues exist as a thiolate-imidazolium ion pair), and then investigated the citrullination reaction of peptidylarginine, using a combined QM/MM approach. State N is calculated to be more stable than State I by 8.46 kcal/mol, and State N can transform to State I via two steps of substrate-assisted proton transfer. Citrullination of the peptidylarginine contains deamination and hydrolysis. Starting from State N, the deamination reaction corresponds to an energy barrier of 18.82 kcal/mol. The deprotonated C351 initiates the nucleophilic attack to the substrate, which is the key step for deamination reaction. The hydrolysis reaction contains two chemical steps. Both the deprotonated D238 and H236 can act as the bases to activate the hydrolytic water, which correspond to similar energy barriers (∼17 kcal/mol). On the basis of our calculations, C351, D238, and H236 constitute a catalytic triad, and their protonation states are critical for both the deamination and hydrolysis processes. In view of the sequence similarity, these findings may be shared with human PAD1-PAD4 and other guanidino-group-modifying enzymes. Proteins 2017; 85:1518-1528. © 2017 Wiley Periodicals, Inc.

How is acetylcholinesterase phosphonylated by soman? An ab initio QM/MM molecular dynamics study

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nerve system that hydrolyzes acetylcholine (ACh) and terminates synaptic signals by reducing the effective concentration of ACh in the synaptic clefts. Organophosphate compounds irreversibly inhibit AChEs, leading to irreparable damage to nerve cells. By employing Born–Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the covalent inhibition mechanism between AChE and the nerve toxin soman and determined its free energy profile for the first time. Our results indicate that phosphonylation of the catalytic serine by soman employs an addition–elimination mechanism, which is highly associative and stepwise: in the initial addition step, which is also rate-limiting, His440 acts as a general base to facilitate the nucleophilic attack of Ser200 on the soman’s phosphorus atom to form a trigonal bipyrimidal pentacovalent intermediate; in the subsequent elimination step, Try121 of the catalytic gorge stabilizes the leaving fluorine atom prior to its dissociation from the active site. Together with our previous characterization of the aging mechanism of soman inhibited AChE, our simulations have revealed detailed molecular mechanistic insights into the damaging function of the nerve agent soman.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

A QM/MM MD study of the pH-dependent ring-opening catalysis and lid motif flexibility in glucosamine 6-phosphate deaminase

QM/MM MD and MM MD simulations reveal pH-dependent proton-shuttle ring-opening mechanisms of GlcN6P and dynamical behavior of the lid motif inSmuNagB.

A full picture of enzymatic catalysis by hydroxynitrile lyases from Hevea brasiliensis: protonation dependent reaction steps and residue-gated movement of the substrate and the product

Hydroxynitrile lyases (HNLs) defend plants from herbivores and microbial attack by releasing cyanide from hydroxynitriles.

Benchmark Study for the Cysteine-Histidine Proton Transfer Reaction in a Protein Environment: Gas Phase, COSMO, QM/MM Approaches

Proton transfer reactions are of crucial interest for the investigation of proteins. We have investigated the accuracy of commonly used quantum chemical methods for the description of proton transfer reactions in different environments (gas phase, COSMO, QM/MM) using the proton transfer between the catalytic dyad residues cysteine 145 and histidine 41 of SARS coronavirus main protease as a case study. The test includes thermodynamic, kinetic, and structural properties. The study comprises computationally demanding ab initio approaches (HF, CC2, MP2, SCS-CC2, SCS-MP2, CCSD(T)), popular density functional theories (BLYP, B3LYP, M06-2X), and semiempirical methods (MNDO/d, AM1, RM1, PM3, PM6). The approximated coupled cluster approach LCCSD(T) is taken as a reference method. We find that the robustness of the tested methods with respect to the environment correlates well with the level of theory. As an example HF, CC2, MP2, and their SCS variants show similar errors for gas phase, COSMO, or QM/MM computations. In contrast for semiempirical methods, the errors strongly diversify if one goes from gas phase to COSMO or QM/MM. Particular problems are observed for the recent semiempirical methods PM6 and RM1, which show the best performance for gas phase calculations but possess larger errors in conjunction with COSMO. Finally, a combination of SCS-MP2 and B3LYP or M06-2X allows reliable estimates about remaining errors.

Modeling methods for studying post-translational and transcriptional modifying enzymes

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.

A catalytic mechanism for cysteine N-terminal nucleophile hydrolases, as revealed by free energy simulations

The N-terminal nucleophile (Ntn) hydrolases are a superfamily of enzymes specialized in the hydrolytic cleavage of amide bonds. Even though several members of this family are emerging as innovative drug targets for cancer, inflammation, and pain, the processes through which they catalyze amide hydrolysis remains poorly understood. In particular, the catalytic reactions of cysteine Ntn-hydrolases have never been investigated from a mechanistic point of view. In the present study, we used free energy simulations in the quantum mechanics/molecular mechanics framework to determine the reaction mechanism of amide hydrolysis catalyzed by the prototypical cysteine Ntn-hydrolase, conjugated bile acid hydrolase (CBAH). The computational analyses, which were confirmed in water and using different CBAH mutants, revealed the existence of a chair-like transition state, which might be one of the specific features of the catalytic cycle of Ntn-hydrolases. Our results offer new insights on Ntn-mediated hydrolysis and suggest possible strategies for the creation of therapeutically useful inhibitors.

Insights into the phosphoryl transfer mechanism of cyclin-dependent protein kinases from ab initio QM/MM free-energy studies

Phosphorylation reactions catalyzed by kinases and phosphatases play an indispensible role in cellular signaling, and their malfunctioning is implicated in many diseases. A better understanding of the catalytic mechanism will help design novel and effective mechanism-based inhibitors of these enzymes. In this work, ab initio quantum mechanical/molecular mechanical studies are reported for the phosphoryl transfer reaction catalyzed by a cyclin-dependent kinase, CDK2. Our results suggest that an active-site Asp residue, rather than ATP as previously proposed, serves as the general base to activate the Ser nucleophile. The corresponding transition state features a dissociative, metaphosphate-like structure, stabilized by the Mg(2+) ion and several hydrogen bonds. The calculated free-energy barrier is consistent with experimental values. Implications of our results in this and other protein kinases are discussed.

On the mechanism of dimethylarginine dimethylaminohydrolase inactivation by 4-halopyridines

Small molecules capable of selective covalent protein modification are of significant interest for the development of biological probes and therapeutics. We recently reported that 2-methyl-4-bromopyridine is a quiescent affinity label for the nitric oxide controlling enzyme dimethylarginine dimethylaminohydrolase (DDAH) (Johnson, C. M.; Linsky, T. W.; Yoon, D. W.; Person, M. D.; Fast, W. J. Am. Chem. Soc. 2011, 133, 1553-1562). Discovery of this novel protein modifier raised the possibility that the 4-halopyridine motif may be suitable for wider application. Therefore, the inactivation mechanism of the related compound 2-hydroxymethyl-4-chloropyridine is probed here in more detail. Solution studies support an inactivation mechanism in which the active site Asp66 residue stabilizes the pyridinium form of the inactivator, which has enhanced reactivity toward the active site Cys, resulting in covalent bond formation, loss of the halide, and irreversible inactivation. A 2.18 Å resolution X-ray crystal structure of the inactivated complex elucidates the orientation of the inactivator and its covalent attachment to the active site Cys, but the structural model does not show an interaction between the inactivator and Asp66. Molecular modeling is used to investigate inactivator binding, reaction, and also a final pyridinium deprotonation step that accounts for the apparent differences between the solution-based and structural studies with respect to the role of Asp66. This work integrates multiple approaches to elucidate the inactivation mechanism of a novel 4-halopyridine "warhead," emphasizing the strategy of using pyridinium formation as a "switch" to enhance reactivity when bound to the target protein.