Letif Mones

On the Divalent Metal Ion Dependence of DNA Cleavage by Restriction Endonucleases of the EcoRI Family

Restriction endonucleases of the PD...D/EXK family need Mg(2+) for DNA cleavage. Whereas Mg(2+) (or Mn(2+)) promotes catalysis, Ca(2+) (without Mg(2+)) only supports DNA binding. The role of Mg(2+) in DNA cleavage by restriction endonucleases has elicited many hypotheses, differing mainly in the number of Mg(2+) involved in catalysis. To address this problem, we measured the Mg(2+) and Mn(2+) concentration dependence of DNA cleavage by BamHI, BglII, Cfr10I, EcoRI, EcoRII (catalytic domain), MboI, NgoMIV, PspGI, and SsoII, which were reported in co-crystal structure analyses to bind one (BglII and EcoRI) or two (BamHI and NgoMIV) Me(2+) per active site. DNA cleavage experiments were carried out at various Mg(2+) and Mn(2+) concentrations at constant ionic strength. All enzymes show a qualitatively similar Mg(2+) and Mn(2+) concentration dependence. In general, the Mg(2+) concentration optimum (between approximately 1 and 10 mM) is higher than the Mn(2+) concentration optimum (between approximately 0.1 and 1 mM). At still higher Mg(2+) or Mn(2+) concentrations, the activities of all enzymes tested are reduced but can be reactivated by Ca(2+). Based on these results, we propose that one Mg(2+) or Mn(2+) is critical for restriction enzyme activation, and binding of a second Me(2+) plays a role in modulating the activity. Steady-state kinetics carried out with EcoRI and BamHI suggest that binding of a second Mg(2+) or Mn(2+) mainly leads to an increase in K(m), such that the inhibitory effect of excess Mg(2+) or Mn(2+) can be overcome by increasing the substrate concentration. Our conclusions are supported by molecular dynamics simulations and are consistent with the structural observations of both one and two Me(2+) binding to these enzymes.

The Energy Gap as a Universal Reaction Coordinate for the Simulation of Chemical Reactions

The selection of a proper reaction coordinate is a major bottleneck in simulations of chemical reactions in complex systems. Increasing the number of variables that are used to bias the reaction largely affects the convergence and leads to an unbearable increase in computational price. This problem can be overcome by employing a complex reaction coordinate that depends on many geometrical variables of the system, such as the energy gap (EGAP) in the empirical valence bond (EVB) method. EGAP depends on all of the coordinates of the system, and its robustness has been demonstrated for a variety of enzymatic reactions. In this work, we demonstrate that EGAP, derived from a classical representation, can be used as a reaction coordinate in systems described with any quantum chemistry Hamiltonian. Benefits of using EGAP as a reaction coordinate as compared to a traditional geometrical variable are illustrated in the case of a symmetric nucleophilic substitution reaction in water solution. EGAP is shown to provide a significantly more efficient sampling and allows a better localization of the transition state as compared to a geometrical reaction coordinate.

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.

Tests of an adaptive QM/MM calculation on free energy profiles of chemical reactions in solution.

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Optimization of reorganization energy drives evolution of the designed Kemp eliminase KE07

Understanding enzymatic evolution is essential to engineer enzymes with improved activities or to generate enzymes with tailor-made activities. The computationally designed Kemp eliminase KE07 carries out an unnatural reaction by converting of 5-nitrobenzisoxazole to cyanophenol, but its catalytic efficiency is significantly lower than those of natural enzymes. Three series of designed Kemp eliminases (KE07, KE70, KE59) were shown to be evolvable with considerable improvement in catalytic efficiency. Here we use the KE07 enzyme as a model system to reveal those forces, which govern enzymatic evolution and elucidate the key factors for improving activity. We applied the Empirical Valence Bond (EVB) method to construct the free energy pathway of the reaction in the original KE07 design and the evolved R7 1/3H variant. We analyzed catalytic effect of residues and demonstrated that not all mutations in evolution are favorable for activity. In contrast to the small decrease in the activation barrier, in vitro evolution significantly reduced the reorganization energy. We developed an algorithm to evaluate group contributions to the reorganization energy and used this approach to screen for KE07 variants with potential for improvement. We aimed to identify those mutations that facilitate enzymatic evolution, but might not directly increase catalytic efficiency. Computational results in accord with experimental data show that all mutations, which appear during in vitro evolution were either neutral or favorable for the reorganization energy. These results underscore that distant mutations can also play role in optimizing efficiency via their contribution to the reorganization energy. Exploiting this principle could be a promising strategy for computer-aided enzyme design. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

The role of reorganization energy in rational enzyme design

Computational design is becoming an integral component in developing novel enzymatic activities. Catalytic efficiencies of man-made enzymes however are far behind their natural counterparts. The discrepancy between laboratory and naturally evolved enzymes suggests that a major catalytic factor is still missing in the computational process. Reorganization energy, which is the origin of catalytic power of natural enzymes, has not been exploited yet for design. As exemplified in case of KE07 Kemp eliminase, this quantity is optimized by directed evolution. Mutations beneficial for evolution, but without direct impact on catalysis can be identified based on contributions to reorganization energy. We propose to incorporate the reorganization energy in scaffold selection to provide highly evolvable initial designs.

Metal-binding sites at the active site of restriction endonuclease BamHI can conform to a one-ion mechanism.

Abstract The number of metal ions required for phosphoryl transfer in restriction endonucleases is still an unresolved question in molecular biology. The two Ca2+ and Mn2+ ions observed in the pre- and post-reactive complexes of BamHI conform to the classical two-metal ion choreography. We probed the Mg2+ cofactor positions at the active site of BamHI by molecular dynamics simulations with one and two metal ions present and identified several catalytically relevant sites. These can mark the pathway of a single ion during catalysis, suggesting its critical role, while a regulatory function is proposed for a possible second ion.

Empirical valence bond simulations of the chemical mechanism of ATP to cAMP conversion by anthrax edema factor.

The two-metal catalysis by the adenylyl cyclase domain of the anthrax edema factor toxin was simulated using the empirical valence bond (EVB) quantum mechanical/molecular mechanical approach. These calculations considered the energetics of the nucleophile deprotonation and the formation of a new P-O bond in aqueous solution and in the enzyme-substrate complex present in the crystal structure models of the reactant and product states of the reaction. Our calculations support a reaction pathway that involves metal-assisted transfer of a proton from the nucleophile to the bulk aqueous solution followed by subsequent formation of an unstable pentavalent intermediate that decomposes into cAMP and pyrophosphate (PPi). This pathway involves ligand exchange in the first solvation sphere of the catalytic metal. At 12.9 kcal/mol, the barrier for the last step of the reaction, the cleavage of the P-O bond to PPi, corresponds to the highest point on the free energy profile for this reaction pathway. However, this energy is too close to the value of 11.4 kcal/mol calculated for the barrier of the nucleophilic attack step to reach a definitive conclusion about the rate-limiting step. The calculated reaction mechanism is supported by reasonable agreement between the experimental and calculated catalytic rate constant decrease caused by the mutation of the active site lysine 346 to arginine.

A universal preconditioner for simulating condensed phase materials

We introduce a universal sparse preconditioner that accelerates geometry optimisation and saddle point search tasks that are common in the atomic scale simulation of materials. Our preconditioner is based on the neighbourhood structure and we demonstrate the gain in computational efficiency in a wide range of materials that include metals, insulators, and molecular solids. The simple structure of the preconditioner means that the gains can be realised in practice not only when using expensive electronic structure models but also for fast empirical potentials. Even for relatively small systems of a few hundred atoms, we observe speedups of a factor of two or more, and the gain grows with system size. An open source Python implementation within the Atomic Simulation Environment is available, offering interfaces to a wide range of atomistic codes.

Preconditioners for the geometry optimisation and saddle point search of molecular systems

A class of preconditioners is introduced to enhance geometry optimisation and transition state search of molecular systems. We start from the Hessian of molecular mechanical terms, decompose it and retain only its positive definite part to construct a sparse preconditioner matrix. The construction requires only the computation of the gradient of the corresponding molecular mechanical terms that are already available in popular force field software packages. For molecular crystals, the preconditioner can be combined straightforwardly with the exponential preconditioner recently introduced for periodic systems. The efficiency is demonstrated on several systems using empirical, semiempirical and ab initio potential energy surfaces.