Abstract:Enzyme catalysis generally includes multiple physical and chemical processes, such as the substrate transportion and binding to the active site, the chemical reaction in the protein confined environment, and the product release, where any chemical or nonchemical step is probably to dominate the activity of enzyme. In order to understand the enzymatic efficiency, the classical molecular dynamics (MM MD) and combined quantum and classical molecular dynamics (QM/MM MD) simulations, in combination with the umbrella sampling, have been used to explore mechanisms for the enzymatic catalysis and the substrate transportion in several biological systems, including the detoxification of the G-type nerve agent sarin by phosphotriesterase (PTE), the hydrolysis of uridine by the pyrimidine-specific nucleoside hydrolase Yeik (CU-NH), and the transportation of naphthalimide-polyamine compounds into the drug sites by the bovine serum albumin (BSA).Extensive QM/MM MD and MM MD simulations show that the wild-type PTE enzyme is capable of hydrolyzing sarin but with a low catalytic efficiency. Dissociation of the degraded product from the binuclear zinc center is predicted to be the rate-determining step with the free-energy barrier of 21.0 kcal/mol, and the multiple chemical steps are involved in the dissociation process. The P-F cleavage of the substrate sarin follows a two-step mechanism, in which the nucleophilic attack of the bridging hydroxide leads to the P-O-mu formation with the free-energy barrier of 9.8 kcal/mol, and the subsequent P-F cleavage and the proton transfer from the mu-OH to the residue Asp301 yields the degraded product with the free-energy span of 12.3 kcal/mol. The relatively strong hydrogen-bond interactions between the substrate uridine and key residues in the active site of CU-NH play an important role in the substrate specificity. The enzymatic N-glycosidic bond cleavage of uridine by CU-NH is the rate-limiting step with the free energy barrier of about 20 kcal/mol, and the substrate protonation less influences the activity of CU-NH, completely different from other nucleoside hydrolases. The combined random acceleration molecular dynamics and MD simulations (RAMD MD) have been used to explore possible transportation channels and the most favorable pathways for the delivery of the naphthalimide-polyamine complex to two drug sites of BSA, as well as their thermodynamic and dynamic properties. MM MD and RAMD MD simulations show that the selective specificity of these two drug sites strongly depends on the electron-withdrawing and electron-donating substituents of the substrate, which can be used for the design of new naphthalimide drugs. Such site-selectivity-dependence design strategy may open up an avenue in drug discovery.With the fast development of low-scaling QM methodologies and the application of artificial intelligence neural network approaches to theoretical and computational chemistry, it can be expected to reduce the computational costs of QM treatment and improve the efficiency of QM/MM, which makes the direct QM/MM study of larger and more complex systems become possible. In addition, the improvement and development of coarse-grained models can surmount the spatial and temporal limitation of multiscale calculations and simulations. Accordingly, global simulations of extremely large and complicated biosystems can be realized in the future, which may provide a theoretical basis for further expanding the application of enzyme engineering.

Is inhibition process better described with MD(QM/MM) simulations? The case of urokinase type plasminogen activator inhibitors

Evaluation of Interactions Between Urokinase Plasminogen and Inhibitors Using Molecular Dynamic Simulation and Free-Energy Calculation.

In Silico Screening of Natural Products As Upar Inhibitors Via Multiple Structure-Based Docking and Molecular Dynamics Simulations.

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations

Understanding the Selectivity of Inhibitors Toward PI4KIIIα and PI4KIIIβ Based Molecular Modeling

Crystal Structures of Urokinase-Type Plasminogen Activator in Complex with 4-(Aminomethyl) Benzoic Acid and 4-(Aminomethyl-phenyl)-methanol

Challenges for drug discovery - a case study of urokinase receptor inhibition.

Understanding the Molecular Mechanism of Binding Modes of Aurora A Inhibitors by Long Time Scale Gpu Dynamics

Structural basis for therapeutic intervention of uPA/uPAR system

Distinctive Binding Modes and Inhibitory Mechanisms of Two Peptidic Inhibitors of Urokinase-Type Plasminogen Activator with Isomeric P1 Residues.

Probing the interaction mechanism of small molecule inhibitors with matriptase based on molecular dynamics simulation and free energy calculations.

Qm/Mm Study Of The Enzymatic Catalysis

Small Molecules Engage Hot Spots Through Cooperative Binding to Inhibit a Tight Protein-Protein Interaction.

Structure-based and multiple potential three-dimensional quantitative structure-activity relationship (SB-MP-3D-QSAR) for inhibitor design.

Structure-Based Engineering Of Species Selectivity In The Interaction Between Urokinase And Its Receptor Implication For Preclinical Cancer Therapy

Comparative Molecular Field Analysis Of Non-Peptidic Inhibitors Of Thrombin

Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation

Binding Affinity and Dissociation Pathway Predictions for a Series of USP7 Inhibitors with Pyrimidinone Scaffold by Multiple Computational Methods.

Probing the Structural Requirements of A-type Aurora Kinase Inhibitors Using 3D-QSAR and Molecular Docking Analysis

Molecular Modeling Study of the Induced-Fit Effect on Kinase Inhibition: the Case of Fibroblast Growth Factor Receptor 3 (FGFR3).

Structure-based QSAR, Molecule Design and Bioassays of Protease-Activated Receptor 1 Inhibitors