The first branching point in porphyrin biosynthesis: a systematic docking, molecular dynamics and quantum mechanical/molecular mechanical study of substrate binding and mechanism of uroporphyrinogen-III decarboxylase

In humans, uroporphyrinogen decarboxylase is intimately involved in the synthesis of heme, where the decarboxylation of the uroporphyrinogen-III occurs in a single catalytic site. Several variants of the mechanistic proposal exist; however, the exact mechanism is still debated. Thus, using an ONIOM quantum mechanical/molecular mechanical approach, the mechanism by which uroporphyrinogen decarboxylase decarboxylates ring D of uroporphyrinogen-III has been investigated. From the study performed, it was found that both Arg37 and Arg50 are essential in the decarboxylation of ring D, where experimentally both have been shown to be critical to the catalytic behavior of the enzyme. Overall, the reaction was found to have a barrier of 10.3 kcal mol(-1) at 298.15 K. The rate-limiting step was found to be the initial proton transfer from Arg37 to the substrate before the decarboxylation. In addition, it has been found that several key interactions exist between the substrate carboxylate groups and backbone amides of various active site residues as well as several other functional groups.     

The MM2QM tool for combining docking,molecular dynamics,molecular mechanics,and quantum mechanics†

The use of the MM2QM tool in a combined docking + molecular dynamics (MD) + molecular mechanics (MM) + quantum mechanical (QM) binding affinity prediction study is presented, and the tool itself is discussed. The system of interest is Mycobacterium tuberculosis (MTB) pantothenate synthetase in complexes with three highly similar sulfonamide inhibitors, for which crystal structures are available. Starting from the structure of MTB pantothenate synthetase in the “open” conformation and following the combined docking + MD + MM + QM procedure, we were able to capture the closing of the enzyme binding pocket and to reproduce the position of the ligands with an average root mean square deviation of 1.6 Å. Protein–ligand interaction energies were reproduced with an average error lower than 10%. The discussion on the MD part and a protein flexibility importance is carried out. The presented approach may be useful especially for finding analog inhibitors or improving drug candidates. © 2012 Wiley Periodicals, Inc.     

Calculation of wave‐functions with frozen orbitals in mixed quantum mechanics/molecular mechanics methods. II. Application of the local basis equation

The application of the local basis equation (Ferenczy and Adams, J. Chem. Phys. 2009 , 130, 134108) in mixed quantum mechanics/molecular mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) methods is investigated. This equation is suitable to derive local basis nonorthogonal orbitals that minimize the energy of the system and it exhibits good convergence properties in a self‐consistent field solution. These features make the equation appropriate to be used in mixed QM/MM and QM/QM methods to optimize orbitals in the field of frozen localized orbitals connecting the subsystems. Calculations performed for several properties in divers systems show that the method is robust with various choices of the frozen orbitals and frontier atom properties. With appropriate basis set assignment, it gives results equivalent with those of a related approach [G. G. Ferenczy previous paper in this issue] using the Huzinaga equation. Thus, the local basis equation can be used in mixed QM/MM methods with small size quantum subsystems to calculate properties in good agreement with reference Hartree–Fock–Roothaan results. It is shown that bond charges are not necessary when the local basis equation is applied, although they are required for the self‐consistent field solution of the Huzinaga equation based method. Conversely, the deformation of the wave‐function near to the boundary is observed without bond charges and this has a significant effect on deprotonation energies but a less pronounced effect when the total charge of the system is conserved. The local basis equation can also be used to define a two layer quantum system with nonorthogonal localized orbitals surrounding the central delocalized quantum subsystem. © 2013 Wiley Periodicals, Inc.     

Convergence in the QM‐only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase

We report systematic quantum mechanics‐only (QM‐only) and QM/molecular mechanics (MM) calculations on an enzyme‐catalyzed reaction to assess the convergence behavior of QM‐only and QM/MM energies with respect to the size of the chosen QM region. The QM and MM parts are described by density functional theory (typically B3LYP/def2‐SVP) and the CHARMM force field, respectively. Extending our previous work on acetylene hydratase with QM regions up to 157 atoms (Liao and Thiel, J. Chem. Theory Comput. 2012, 8, 3793), we performed QM/MM geometry optimizations with a QM region M4 composed of 408 atoms, as well as further QM/MM single‐point calculations with even larger QM regions up to 657 atoms. A charge deletion analysis was conducted for the previously used QM/MM model ( M3a , with a QM region of 157 atoms) to identify all MM residues with strong electrostatic contributions to the reaction energetics (typically more than 2 kcal/mol), which were then included in M4 . QM/MM calculations with this large QM region M4 lead to the same overall mechanism as the previous QM/MM calculations with M3a , but there are some variations in the relative energies of the stationary points, with a mean absolute deviation (MAD) of 2.7 kcal/mol. The energies of the two relevant transition states are close to each other at all levels applied (typically within 2 kcal/mol), with the first (second) one being rate‐limiting in the QM/MM calculations with M3a ( M4 ). QM‐only gas‐phase calculations give a very similar energy profile for QM region M4 (MAD of 1.7 kcal/mol), contrary to the situation for M3a where we had previously found significant discrepancies between the QM‐only and QM/MM results (MAD of 7.9 kcal/mol). Extension of the QM region beyond M4 up to M7 (657 atoms) leads to only rather small variations in the relative energies from single‐point QM‐only and QM/MM calculations (MAD typically about 1–2 kcal/mol). In the case of acetylene hydratase, a model with 408 QM atoms thus seems sufficient to achieve convergence in the computed relative energies to within 1–2 kcal/mol.Copyright © 2013 Wiley Periodicals, Inc.     

The adaptive buffered force QM/MM method in the CP2K and AMBER software packages

The implementation and validation of the adaptive buffered force (AdBF) quantum‐mechanics/molecular‐mechanics (QM/MM) method in two popular packages, CP2K and AMBER are presented. The implementations build on the existing QM/MM functionality in each code, extending it to allow for redefinition of the QM and MM regions during the simulation and reducing QM‐MM interface errors by discarding forces near the boundary according to the buffered force‐mixing approach. New adaptive thermostats, needed by force‐mixing methods, are also implemented. Different variants of the method are benchmarked by simulating the structure of bulk water, water autoprotolysis in the presence of zinc and dimethyl‐phosphate hydrolysis using various semiempirical Hamiltonians and density functional theory as the QM model. It is shown that with suitable parameters, based on force convergence tests, the AdBF QM/MM scheme can provide an accurate approximation of the structure in the dynamical QM region matching the corresponding fully QM simulations, as well as reproducing the correct energetics in all cases. Adaptive unbuffered force‐mixing and adaptive conventional QM/MM methods also provide reasonable results for some systems, but are more likely to suffer from instabilities and inaccuracies. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Integrated QMMM and Monte Carlo methods for analysis of adsorptive interactions between goethite cluster,carbon nanotubes,and arsenate

Iron‐immobilized nanoporous carbon is a well‐known adsorbent used in treating arsenic‐contaminated waters. In this contribution, we present findings on the adsorptive interactions and dynamics of arsenate–goethite cluster ([FeO(OH)]₆) with carbon nanotubes (CNTs) using hybridized quantum mechanics/molecular mechanics (QMMM) calculations. The CNTs adsorption mechanism is of interest since a better understanding of the fundamental interactions between arsenate, goethite, and carbon surfaces would translate to advances in CNT‐based adsorbent production and utilization. Novel applications of general amber force field (GAFF) and isobaric‐isothermal Gibbs ensemble Monte Carlo (NpT‐GEMC) methods are described. By the abovementioned methods, we postulate that the [FeO(OH)]₆/CNT‐2.3 (diameter 2.3 nm ‐ mesoporous) system enhances the qualitative (i.e., improved chemisorption) rather than the quantitative adsorptive aspect (i.e., total ions adsorbed) in comparison to the [FeO(OH)]₆/CNT‐1.6 (diameter 1.6 nm ‐ microporous) system.     

From the Micro- to the Macro-: Managing the Conservation of the Warship,Vasa

Sweden's famous warship, Vasa, sank on her maiden voyage in August 1628, and remained on the bottom of Stockholm harbour for 333 years. Raised in 1961, she became the first large-scale wooden object to be treated with polyethylene glycol (PEG). In the summer of 2000 a number of acidic salt precipitations were noticed on the surface of the ship and on wooden artefacts in the storerooms. An international research project has been established to look into the causes of this problem and suggest possible re-treatments. Meanwhile projects are underway to monitor movements in the ship, to build a better support system, and to replace the thousands of iron bolts holding the structure together, while a sophisticated new climate system has recently been installed in the museum.     

Computation of internal coordinates,derivatives, and gradient expressions: torsion and improper torsion

Laplacians and gradient dot products are required for the recently developed internal coordinate quantum Monte Carlo method. New formulas are presented for these quantities for torsion and improper torsion angles. The Laplacians can also be used to economize calculation of sets of second derivatives used in molecular mechanics and other methods. Formulas for torsion angle gradient dot products and Laplacians, and completely new formulas for improper torsion, are presented. In addition, calculations of cos τ and sin τ, some suitable for energy subroutines and others for force subroutines, are shown. Finally, in a related development, several sets of conditions for three atom linearity or four atom planarity involving internal coordinate derivatives are reported. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 553–561, 2000     

Probing the range of applicability of structure‐ and energy‐adjusted QM/MM link bonds

The hydrogen‐capping method is one of the most popular and widely used coupling‐schemes for quantum mechanics/molecular mechanics (QM/MM)‐molecular dynamics simulations of macromolecular systems. This is mostly due to the fact that it is fairly convenient to implement and parametrize, thus providing an excellent compromise between accuracy and computational effort. In this work, a viable and straight‐forward approach to optimize the placing of the link atom on a suitable distance ratio between the frontier atoms is discussed. To further increase the accuracy, instead of global parameters for all amino acids, different parameter sets for each type of amino acid are derived. The dependency of the link bond parameters on the chemical environment and the used QM‐method is probed to assess the range of applicability of the parametrization. Suitable sets of parameters for RI‐MP2, B3LYP, (RI)‐B3LYP‐D3, and RI‐BLYP‐D3 at triple‐zeta level for all relevant proteinogenic amino acids are presented. Furthermore, the scope and range of the perturbation, stemming from the introduction of link bonds is evaluated through application of the presented QM/MM scheme in calculations of the active site of 15S‐lipoxygenase. © 2015 Wiley Periodicals, Inc.     

Machine learning and signed particles,an alternative and efficient way to simulate quantum systems

Recently neural networks have been applied in the context of the signed particle formulation of quantum mechanics to rapidly and reliably compute the Wigner kernel of any provided potential. Important advantages were introduced, such as the reduction of the amount of memory required for the simulation of a quantum system by avoiding the storage of the kernel in a multi-dimensional array, as well as attainment of consistent speedup by the ability to realize the computation only on the cells occupied by signed particles. An inherent limitation was the number of hidden neurons to be equal to the number of cells of the discretized real space. In this work, anew network architecture is presented, decreasing the number of neurons in its hidden layer, thereby reducing the complexity of the network and achieving an additional speedup. The approach is validated on a onedimensional quantum system consisting of a Gaussian wave packet interacting with a potential barrier.     

The Catalytic Mechanism of the Marine‐Derived Macrocyclase PatGmac

Cyclic peptides are a class of compounds with high therapeutic potential, possessing bioactivities including antitumor and antiviral (including anti‐HIV). Despite their desirability, efficient design and production of these compounds has not been achieved to date. The catalytic mechanism of patellamide macrocyclization by the PatG macrocyclase domain has been computationally investigated by using quantum mechanics/molecular mechanics methodology, specifically ONIOM(M06/6‐311++G(2d,2p):ff94//B3LYP/6‐31G(d):ff94). The mechanism proposed herein begins with a proton transfer from Ser783 to His 618 and from the latter to Asp548. Nucleophilic attack of Ser783 on the substrate leads to the formation of an acyl–enzyme covalent complex. The leaving group Ala‐Tyr‐Asp‐Gly (AYDG) of the substrate is protonated by the substrate's N terminus, leading to the breakage of the P1?P1′ bond. Finally, the substrate's N terminus attacks the P1 residue, decomposing the acyl–enzyme complex forming the macrocycle. The formation and decomposition of the acyl–enzyme complex have the highest activation free energies (21.1 kcal mol^?1 and 19.8 kcal mol^?1 respectively), typical of serine proteases. Understanding the mechanism behind the macrocyclization of patellamides will be important to the application of the enzymes in the pharmaceutical and biotechnological industries.     

共轭烯烃的分子几何构型、电子结构和生成热的DMM和AM1、PM3的比较研究

运用Ｄｅｌｆｔ分子力学（ＤＭＭ）力场和程序以及半经验分子轨道ＡＭ１和ＰＭ３方法计算研究了丁二烯、苯、甲苯、联苯、苯乙烯、富烯、、环辛四烯、［２，２］对环烷和菲等１０个共轭烯烃分子的几何构型、电子结构和生成热．ＤＭＭ计算的几何构型和生成热与实验结果相吻合，电荷分布结果与从头计算结果较接近．ＡＭ１和ＰＭ３计算的几何构型较好，但计算的生成热与实验结果偏差较大．ＰＭ３计算值比ＡＭ１的稍好．     

Integrating quantum and molecular mechanics

A computer algorithm is developed for integrating density functional quantum mechanics into a molecular mechanics program. The computationally infeasible aspects of the standard LCAO-MO approach (the iterative calculation of eigenvectors and the requirement of orthogonal expansions for the orbitals) are replaced with an efficient use of optimization via the trace theorem of linear algebra. The construction of a basis is also described for expanding the electron density that transforms with the molecular geometry. The combination of the trace method and the basis allow the solution for one configuration of atoms and electrons to be tracked over a wide range of internal conformations. The approach is readily adaptable to being used in the context of an imposed classical field that allows it to be used on part of a macromolecular complex. The initial implementation in the program AMMP is described. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1618–1633, 1999     

Molecular modeling of intercalation complexes of antitumor active 9-aminoacridine and a [d,e]-anellated isoquinoline derivative with base paired deoxytetranucleotides

Summary Intercalators are molecules capable of sliding between DNA base pairs without breaking up the hydrogen bonds between the DNA bases. On the basis of molecular mechanics calculations structural, models of B-DNA tetranucleotide intercalation complexes of some cytostatic active 9-aminoacridines and of a [d, e]-anellated isoquinoline derivative are presented. The drug complexes are stabilized by energetically favouredvan der Waals interactions and by selective hydrogen bonds between the side chains of the drugs and the DNA bases. Semiempirical quantum chemistry calculations revealed that the chromophoric system of the intercalators is able to form ,-charge-transfer interactions with the purine bases of the base paired deoxytetranucleotides. The theoretical findings are of interest for a more specific drug design of cytostatically active agents.

---

Molecular Modeling von Interkalationskomplexen antitumoraktiver 9-Aminoacridine sowie eines [d, e]-anellierten Isochinolinderivates mit basengepaarten Desoxytetranukleotiden

Zusammenfassung Interkalatoren sind Moleküle, die in der Lage sind, sich zwischen DNA-Basenpaare einzulagern, ohne die Wasserstoffbrücken zwischen den DNA-Basen aufzubrechen. Auf der Basis von molekülmechanischen Rechnungen werden Tetranukleotid-Interkalationskomplexe von verschiedenen zytostatisch aktiven 9-Aminoacridinen und von einem [d, e]-anellierten Isochinolinderivat präsentiert. Die Komplexe werden durch energetisch günstigevan der Waals-Interaktionen sowie durch selektive Wasserstoffbrückenbindungen zwischen den Seitenketten der Wirkstoffe und den DNA-Basen stabilisiert. Semiempirische quantenchemische Rechnungen ergaben, daß der Chromophor der Interkalatoren in der Lage ist, ,-charge-transfer Wechselwirkungen mit den Purinbasen der basengepaarten Desoxytetranukleotide auszubilden. Die theoretischen Ergebnisse sind für ein spezifischeres Wirkstoffdesign zytostatisch aktiver Verbindungen von Interesse.

     

Comparison of radii sets,entropy, QM methods,and sampling on MM‐PBSA,MM‐GBSA,and QM/MM‐GBSA ligand binding energies of F. tularensis enoyl‐ACP reductase (FabI)

To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM‐PBSA, MM‐GBSA, and QM/MM‐GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM‐GBSA using 6 ns MD simulation trajectories together with GB^neck2, PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r² = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called “multiple independent sampling”), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB^HCT and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r² = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive “multiple independent sampling method” may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM‐P(G)BSA were limited to inhibitors’ relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program. © 2015 Wiley Periodicals, Inc.     

Building capacity for undergraduate education and training in computational molecular science: A collaboration between the MERCURY consortium and the Molecular Sciences Software Institute

The Molecular Sciences Software Institute (MolSSI) is an National Science Foundation (NSF) funded institute that focuses on improving software, education, and training in the computational molecular sciences. Through a collaboration with the Molecular Education and Research Consortium in Undergraduate computational chemistRY (MERCURY), the MolSSI has developed resources for undergraduate and other early career students to lay an educational foundation for the next generation of computational molecular scientists. The resources focus on introducing best practices in software engineering to students from the very start to make their software more useable, maintainable, and reproducible.     

Solvents effects on the mechanism of cellulose hydrolysis: A QM/MM study

This article reports a combined quantum mechanics/molecular mechanics (QM/MM) investigation on the acid hydrolysis of cellulose in water using two different models, cellobiose and a 40‐unit cellulose chain. The explicitly treated solvent molecules strongly influence the conformations, intramolecular hydrogen bonds, and exoanomeric effects in these models. As these features are largely responsible for the barrier to cellulose hydrolysis, the present QM/MM results for the pathways and reaction intermediates in water are expected to be more realistic than those from a former density functional theory (DFT) study with implicit solvent (CPCM). However, in a qualitative sense, there is reasonable agreement between the DFT/CPCM and QM/MM predictions for the reaction mechanism. Differences arise mainly from specific solute–solvent hydrogen bonds that are only captured by QM/MM and not by DFT/CPCM. © 2015 Wiley Periodicals, Inc.