Abstract:The discovery of unique substrates is important for developing potential applications of enzymes. However, the experimental procedures for substrate identification are laborious, time-consuming, and expensive. Although in silico structure based approaches show great promise, recent extensive studies have shown that these approaches remain a formidable challenge for current biocomputational methodologies. Here we present an open-source, extensible, and flexible software platform for predicting enzyme substrates called THEMIS, which performs in silico virtual screening for potential catalytic targets of an enzyme on the basis of the enzyme's catalysis mechanism. On the basis of a generalized transition state theory of enzyme catalysis, we introduce a modified docking procedure called "mechanism-based restricted docking" (MBRD) for novel substrate recognition from molecular docking. Comprising a series of utilities written in C/Python, THEMIS automatically executes parallel-computing MBRD tasks and evaluates the results with various molecular mechanics (MM) criteria such as energy, distance, angle, and dihedral angle to help identify desired substrates. Exhaustive sampling and statistical measures were used to improve the robustness and reproducibility of the method. We used Candida antarctica lipase B (CALB) as a test system to demonstrate the effectiveness of our computational prediction of (non)substrates. A novel MM score function for CALB substrate identification derived from the near-attack conformation was used to evaluate the possibility of chemical transformation. A highly positive rate of 93.4% was achieved from a CALB substrate library with 61 known substrates and 35 nonsubstrates, and the screening rate has reached 103 compounds/day (96 CPU cores, 100 samples/compound). The performance shows that the present method is perhaps the first reported scheme to meet the requirement for practical applicability to enzyme studies. An additional study was performed to validate the universality of our method. In this verification we employed two distinct enzymes, nitrilase Nit6803 and SDR Gox2181, where the correct rates of both enzymes exceeded 90%. The source code used will be released under the GNU General Public License (GPLv3) and will be free to download. We believe that the present method will provide new insights into enzyme research and accelerate the development of novel enzyme applications.

Predicting enzyme-substrate specificity with QM/MM methods: a case study of the stereospecificity of (D)-glucarate dehydratase

Quantitative prediction of enantioselectivity of Candida antarctica lipase B by combining docking simulations and quantitative structure-activity relationship (QSAR) analysis

Prediction of the Enantioselectivity of Lipases and Esterases by Molecular Docking Method with Modified Force Field Parameters.

Rational Design of Dehydrogenase/Reductases Based on Comparative Structural Analysis of Prereaction-State and Free-State Simulations for Efficient Asymmetric Reduction of Bulky Aryl Ketones

Quantum Mechanics/Molecular Mechanics Modeling of Regioselectivity of Drug Metabolism in Cytochrome P450 2C9

A Semiautomated Structure-Based Method to Predict Substrates of Enzymes Via Molecular Docking: A Case Study with Candida Antarctica Lipase B

Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450

Molecular Recognition of Metabotropic Glutamate Receptor Type 1 (mglur1): Synergistic Understanding with Free Energy Perturbation and Linear Response Modeling

Toward On-The-Fly Quantum Mechanical/Molecular Mechanical (QM/MM) Docking: Development and Benchmark of a Scoring Function

Understanding the Effectiveness of Enzyme Pre-Reaction State by a Quantum-Based Machine Learning Model

Systematic Quantum Mechanical Region Determination in QM/MM Simulation

Quantum chemical modeling of enzymatic reactions: The case of decarboxylation

Identification of specificity determining residues in enzymes using environment specific substitution tables

How large should the QM region be in QM/MM calculations? The case of catechol O-methyltransferase

Substrate Selectivity for UDP-glucuronosyltransferase1A8 Using the Pharmacophore Approach

A Computational Approach to Enzyme Design: Predicting ω-Aminotransferase Catalytic Activity Using Docking and MM-GBSA Scoring

A Case Study of the Glycoside Hydrolase Enzyme Mechanism Using an Automated QM-Cluster Model Building Toolkit

Understanding Substrate Binding and Reactivity of Stearoyl-CoA Desaturase (SCD1) through Classical and Multiscale Molecular Dynamics Simulations

Harder, better, faster, stronger: large-scale QM and QM/MM for predictive modeling in enzymes and proteins

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations

Two Sides of Quantum-Based Modeling of Enzyme-Catalyzed Reactions: Mechanistic and Electronic Structure Aspects of the Hydrolysis by Glutamate Carboxypeptidase