Extension of QM/MM docking and its applications to metalloproteins

To overcome the limitation of conventional docking methods which assume fixed charge model from force field parameters, combined quantum mechanics/molecular mechanics (QM/MM) method has been applied to docking as a variable charge model and shown to exhibit improvement on the docking accuracy over fixed charge based methods. However, it has also been shown that there are a number of examples for which adoption of variable‐charge model fails to reproduce the native binding modes. In particular, for metalloproteins, previously implemented method of QM/MM docking failed most often. This class of proteins has highly polarized binding sites at which high‐coordinate‐numbered metal ions reside. We extend the QM/MM docking method so that protein atoms surrounding the binding site along with metal ions are included as quantum region, as opposed to only ligand atoms. This extension facilitates the required scaling of partial charges on metal ions leading to prediction of correct binding modes in metalloproteins. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Hybrid QM/MM studies on energetics of malonaldehyde in condensed phase

Energetics of eight enol isomers of a malonaldehyde (MA) at the ground state in CCl₄ solvent environment have been investigated using a hybrid quantum mechanical (QM)/molecular mechanics (MM) method. It is found that relative energies of the isomers slightly change due to interactions between MA and surrounding atoms. In an isolated environment all eight isomers have stable planar structure. On the other hand, most of the isomers have nonplanar structure in CCl₄ whose interactions with the solute molecules are, however, weak. Mainly, structural changes are found in the H? O? C?C dihedral angle ?_{H? O? C?C}, i.e., a hydrogen atom that is connected directly with an oxygen atom, is located at a nonplanar position, and other atoms remain almost planar. Vertical excitation energies of low‐lying excited states at the resultant optimized structure of each isomer are evaluated. The vertical excitation energies in CCl₄ are almost the same as those in the isolated environment, but some changes were found in the triplet excitation states. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

QM/MM method for metal–organic interfaces

We present a QM/MM method for modeling metal/organic interfaces, which incorporates contributions from long‐range electron correlation, characteristic to metals and non‐bonded interactions in organic systems. This method can be used to study structurally irregular systems. We apply the method to model finite size domains of self‐assembled monolayers on the gold (111) surface and discuss the influence of boundary effects on the electrostatic and electronic properties of these systems. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Pseudobond ab initio QM/MM approach and its applications to enzyme reactions

This perspective article mainly focuses on the development and applications of a pseudobond ab initio QM/MM approach to study enzyme reactions. The following aspects of methodology development are discussed: the approaches for the QM/MM covalent boundary problem, an efficient iterative optimization procedure, the methods to determine enzyme reaction paths, and the approaches to calculate free energy change in enzyme reactions. Several applications are described to illustrate the capability of the methods. Finally, future directions are discussed.     

应用QM与ABEEM/MM结合的方法研究NAMI-A的结构性质

应用量子力学(QM)与ABEEM浮动电荷力场(ABEEM/MM)相结合的方法研究了抗癌药物NAMI-A在水溶液中的结构性质. 所有的结构优化都是在DFT的B3LYP方法下采用6-31G(d,p)和LanL2DZ基组完成的, 没有加入任何限制性条件. 结果表明, 优化得到的NAMI-A构型受不同环境及方法的影响均有变化. 与气相中得到的构型相比, QM/MM迭代优化得到构型要比PCM的构型变化更明显. QM/MM (ABEEM/MM)迭代优化得到的NAMI-A构型比QM/MM (OPLS-AA)的变化要小. 总之, 溶剂通过极化效应对NAMI-A结构、电荷分布及径向分布函数等性质均有影响, 客观地处理极化效应才能正确地反映QM区的性质.     

Interfacing Q-Chem and CHARMM to perform QM/MM reaction path calculations

A hybrid quantum mechanical/molecular mechanical (QM/MM) potential energy function with Hartree-Fock, density functional theory (DFT), and post-HF (RIMP2, MP2, CCSD) capability has been implemented in the CHARMM and Q-Chem software packages. In addition, we have modified CHARMM and Q-Chem to take advantage of the newly introduced replica path and the nudged elastic band methods, which are powerful techniques for studying reaction pathways in a highly parallel (i.e., parallel/parallel) fashion, with each pathway point being distributed to a different node of a large cluster. To test our implementation, a series of systems were studied and comparisons were made to both full QM calculations and previous QM/MM studies and experiments. For instance, the differences between HF, DFT, MP2, and CCSD QM/MM calculations of H2O...H2O, H2O...Na+, and H2O...Cl- complexes have been explored. Furthermore, the recently implemented polarizable Drude water model was used to make comparisons to the popular TIP3P and TIP4P water models for doing QM/MM calculations. We have also computed the energetic profile of the chorismate mutase catalyzed Claisen rearrangement at various QM/MM levels of theory and have compared the results with previous studies. Our best estimate for the activation energy is 8.20 kcal/mol and for the reaction energy is -23.1 kcal/mol, both calculated at the MP2/6-31+G(d)//MP2/6-31+G(d)/C22 level of theory.     

Quantifying Electronic Effects in QM and QM/MM Biomolecular Modeling with the Fukui Function

Multi-scale quantum-mechanical/molecular-mechanical(QM/MM) and large-scale QM simulation provide valuable insight into enzyme mechanism and structure-property relationships. Analysis of the electron density afforded through these methods can enhance our understanding of how the enzyme environment modulates reactivity at the enzyme active site. From this perspective, tools from conceptual density functional theory to interrogate electron densities can provide added insight into enzyme function. We recently introduced the highly parallelizable Fukui shift analysis(FSA) method, which identifies how frontier states of an active site are altered by the presence of an additional QM residue to identify when QM treatment of a residue is essential as a result of quantum-mechanically affecting the behavior of the active site. We now demonstrate and analyze distance and residue dependence of Fukui function shifts in pairs of residues representing different non-covalent interactions. We also show how the interpretation of the Fukui function as a measure of relative nucleophilicity provides insight into enzymes that carry out S_N2 methyl transfer. The FSA method represents a promising approach for the systematic, unbiased determination of quantum mechanical effects in enzymes and for other complex systems that necessitate multi-scale modeling.     

Spectral Features of Canthaxanthin in HCP2. A QM/MM Approach

The increased interest in sequencing cyanobacterial genomes has allowed the identification of new homologs to both the N-terminal domain (NTD) and C-terminal domain (CTD) of the Orange Carotenoid Protein (OCP). The N-terminal domain homologs are known as Helical Carotenoid Proteins (HCPs). Although some of these paralogs have been reported to act as singlet oxygen quenchers, their distinct functional roles remain unclear. One of these paralogs (HCP2) exclusively binds canthaxanthin (CAN) and its crystal structure has been recently characterized. Its absorption spectrum is significantly red-shifted, in comparison to the protein in solution, due to a dimerization where the two carotenoids are closely placed, favoring an electronic coupling interaction. Both the crystal and solution spectra are red-shifted by more than 50 nm when compared to canthaxanthin in solution. Using molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) studies of HCP2, we aim to simulate these shifts as well as obtain insight into the environmental and coupling effects of carotenoid–protein interactions.     

QM/MM study of the taxadiene synthase mechanism

Combined quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the reaction mechanism of taxadiene synthase (TXS). TXS catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to taxadiene (T) and four minor cyclic products. All these products originate from the deprotonation of carbocation intermediates. The reaction profiles for the conversion of GGPP to T as well as to minor products were calculated for different configurations of relevant TXS carbocation complexes. The QM region was treated at the M06-2X/TZVP level, while the CHARMM27 force field was used to describe the MM region. The QM/MM calculations suggest a reaction pathway for the conversion of GGPP to T, which slightly differs from previous proposals regarding the number of reaction steps and the conformation of the carbocations. The QM/MM results also indicate that the formation of minor products via water-assisted deprotonation of the carbocations is highly exothermic, by about −7 to −23 kcal/mol. Curiously, however, the computed barriers and reaction energies indicate that the formation of some of the minor products is more facile than the formation of T. Thus, the present QM/MM calculations provide detailed insights into possible reaction pathways and into the origin of the promiscuity of TXS, but they do not reproduce the product distribution observed experimentally. © 2019 Wiley Periodicals, Inc.     

MM and QM/MM Modeling of Threonyl-tRNA Synthetase: Model Testing and Simulations

Aminoacyl-tRNA synthetases are centrally important enzymes in protein synthesis. We have investigated threonyl-tRNA synthetase from E. coli, complexed with reactants, using molecular mechanics and combined quantum mechanical/molecular mechanical (QM/MM) techniques. These modeling methods have the potential to provide molecular level understanding of enzyme catalytic processes. Modeling of this enzyme presents a number of challenges. The procedure of system preparation and testing is described in detail. For example, the number of metal ions at the active site, and their positions, were investigated. Molecular dynamics simulations suggest that the system is most stable when it contains only one magnesium ion, and the zinc ion is removed. Two different QM/MM methods were tested in models based on the findings of MM molecular dynamics simulations. AM1/CHARMM calculations resulted in unrealistic structures for the phosphates in this system. This is apparently due to an error of AM1. PM3/CHARMM calculations proved to be more suitable for this enzyme system. These results will provide a useful basis for future modeling investigations of the enzyme mechanism and dynamics.     

On the effect of a variation of the force field, spatial boundary condition and size of the QM region in QM/MM MD simulations

During the past years, the use of combined quantum-classical, QM/MM, methods for the study of complex biomolecular processes, such as enzymatic reactions and photocycles, has increased considerably. The quality of the results obtained from QM/MM calculations is largely dependent on five aspects to be considered when setting up a molecular model: the QM Hamiltonian, the MM Hamiltonian or force field, the boundary and coupling between the QM and MM regions, the size of the QM region and the boundary condition for the MM region. In this study, we systematically investigate the influence of a variation of the molecular mechanics force field and the size of the QM region in QM/MM MD simulations on properties of the photoactive part of the blue light photoreceptor protein AppA. For comparison, we additionally performed classical MD simulations and studied the effect of a variation of the type of spatial boundary condition. The classical boundary conditions and the force field used in a QM/MM MD simulation are shown to have non-neglegible effects upon the structural and energetic properties of the protein which makes it advisable to minimize computational artifacts in QM/MM MD simulations by application of periodic boundary conditions and a thermodynamically calibrated force field. A comparison of the structural and energetic properties of MD simulations starting from two alternative, different X-ray structures for the blue light utilizing flavin protein in its dark state indicates a slight preference of the two force fields used for the so-called Anderson structure over the Jung structure.     

LICHEM: A QM/MM program for simulations with multipolar and polarizable force fields

We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER–HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single‐point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo‐bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization. © 2016 Wiley Periodicals, Inc.     

Infrared spectroscopy of N-methylacetamide in water from high-level QM/MM calculations简

The infrared(IR) spectra of the N-methylacetamide molecule in water are calculated by using the MD simulation with high-level QM/MM corrections. The B3LYP and MP2 levels with 6-311++G** basis set are used for the QM region, respectively. Our results show all IR spectra at the B3LYP level are well consistent with the corresponding MP2 results. A dynamical charge fluctuation is observed for each atom along the simulation trajectories due to the electrostatic polarization(EP) effects from surrounding solvent environment. We find that the QM/MM corrected IR spectra satisfactorily reproduce the experimental vibrational features of amide I–III modes.     

A QM/MM Study on the Initiation Reaction of Firefly Bioluminescence—Enzymatic Oxidation of Luciferin

Among all bioluminescent organisms, the firefly is the most famous, with a high luminescent efficiency of 41%, which is widely used in the fields of biotechnology, biomedicine and so on. The entire bioluminescence (BL) process involves a series of complicated in-vivo chemical reactions. The BL is initiated by the enzymatic oxidation of luciferin (LH₂). However, the mechanism of the efficient spin-forbidden oxygenation is far from being totally understood. Via MD simulation and QM/MM calculations, this article describes the complete process of oxygenation in real protein. The oxygenation of luciferin is initiated by a single electron transfer from the trivalent anionic LH₂ (L³⁻) to O₂ to form ¹[L^•2−…O₂^•−]; the entire reaction is carried out along the ground-state potential energy surface to produce the dioxetanone (FDO⁻) via three transition states and two intermediates. The low energy barriers of the oxygenation reaction and biradical annihilation involved in the reaction explain this spin-forbidden reaction with high efficiency. This study is helpful for understanding the BL initiation of fireflies and the other oxygen-dependent bioluminescent organisms.     

QM/MM calculation of solvent effects on absorption spectra of guanine

Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited‐state calculations were extracted from ground‐state molecular dynamics (MD) simulations using the self‐consistent‐charge density functional tight binding (SCC‐DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited‐state calculations used time‐dependent density functional theory (TDDFT) and the DFT‐based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto‐N7H and keto‐N9H guanine, with particular focus on solvent effects in the low‐energy spectrum of the keto‐N9H tautomer. When compared with the vertical excitation energies of gas‐phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC‐DFTB‐based rather than B3LYP‐based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas‐phase and solution spectra, typically ca. 0.1–0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent‐induced structural changes and from electrostatic solute–solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Molecular Dynamics and Combined QM/MM Studies on the Deactivation of anti-Tubercular Drug Isoniazid by Arylamine N-Acetyltransferases

Both a molecule dynamic study and a combined quantum mechanics and mole-cule mechanics(QM/MM) study on the acetylating deactivation mechanism of isoniazid were presented.This type of reaction was catalyzed by arylamine N-acetyltransferases(NATs) and the results strongly support a direct acetyl group transfer process rather than a stepwise one.The isoniazid was strictly restrained in proper relative position to accept the acetyl group by a Hydrogen-bond network formed by the residues at the active center.The residues,His110 and Cys70,would be functioned as 'general base' rather than 'general acid'.If all the residues(including H2O molecules) were removed from the system,the activation energy will be increased from 145.1 to 243.3 kJ/mol.The calculations met the experimental data with good agreement.