Exploring the Catalytic Mechanism of Cas9 Using Information Inferred from Endonuclease VII.

Elucidating the nature of the gene editing mechanism of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) is an important task in view of the role of this breakthrough to the advancement of human medicine. In particular, it is crucial to understand the catalytic mechanism of Cas9 (one of the CRISPR associated proteins) and its role in confirming accurate editing. Thus, we focus in this work on an attempt to analyze the catalytic mechanism of Cas9. Considering the absence of detailed structural information on the active form of Cas9, we use an empirical valence bond (EVB) which is calibrated on the closely related mechanism of T4 endonuclease VII. The calibrated EVB is then used in studying the reaction of Cas9, while trying several structural models. It is found that the catalytic activation requires a large conformational change, where K848 or other positively charged group moves from a relatively large distance toward the scissile phosphate. This conformational change leads to the change in position of the Mg2+ ion and to a major reduction in the activation barrier for the catalytic reaction. Our finding provides an important clue on the nature of the catalytic activation of CAS9 and thus should help in elucidating a key aspect of the gene editing process. For example, the approach used here should be effective in exploring the nature of off target activation and its relationship to the energetics of the unwinding process. This strategy may offer ways to improve the selectivity of Cas9.

Validating the Water Flooding Approach by Comparing It to Grand Canonical Monte Carlo Simulations.

The study of the function of proteins on a quantitative level requires consideration of the water molecules in and around the protein. This requirement presents a major computational challenge due to the fact that the insertion of water molecules can have a very high activation barrier and would require a long simulation time. Recently, we developed a water flooding (WF) approach which is based on a postprocessing Monte Carlo ranking of possible water configurations. This approach appears to provide a very effective way for assessing the insertion free energies and determining the most likely configurations of the internal water molecules. Although the WF approach was used effectively in modeling challenging systems that have not been addressed reliably by other microscopic approaches, it was not validated by a comparison to the more rigorous grand canonical Monte Carlo (GCMC) method. Here we validate the WF approach by comparing its performance to that of the GCMC method. It is found that the WF approach reproduces the GCMC results in well-defined test cases but does so much faster. This established the WF approach as a useful strategy for finding correct water configurations in proteins and thus to provide a powerful way for studies of the functions of proteins.

Simulating the fidelity and the three Mg mechanism of pol η and clarifying the validity of transition state theory in enzyme catalysis.

Pol η belongs to the important Y family of DNA polymerases that can catalyze translesion synthesis across sites of damaged DNA. This activity involves the reduced fidelity of Pol η for 8-oxo-7,8-dhyedro-2'-deoxoguanosin(8-oxoG). The fundamental interest in Pol η has grown recently with the demonstration of the importance of a 3rd Mg2+ ion. The current work explores both the fidelity of Pol η and the role of the 3rd metal ion, by using empirical valence bond (EVB) simulations. The simulations reproduce the observed trend in fidelity and shed a new light on the role of the 3rd metal ion. It is found that this ion does not lead to a major catalytic effect, but most probably plays an important role in reducing the product release barrier. Furthermore, it is concluded, in contrast to some implications, that the effect of this metal does not violate transition state theory, and the evaluation of the catalytic effect must conserve the molecular composition upon moving from the reactant to the transition state. Proteins 2017; 85:1446-1453. © 2017 Wiley Periodicals, Inc.

The control of the discrimination between dNTP and rNTP in DNA and RNA polymerase.

Understanding the origin of discrimination between rNTP and dNTP by DNA/RNA polymerases is important both for gaining fundamental knowledge on the corresponding systems and for advancing the design of specific drugs. This work explores the nature of this discrimination by systematic calculations of the transition state (TS) binding energy in RB69 DNA polymerase (gp43) and T7 RNA polymerase. The calculations reproduce the observed trend, in particular when they included the water contribution obtained by the water flooding approach. Our detailed study confirms the idea that the discrimination is due to the steric interaction between the 2'OH and Tyr416 in DNA polymerase, while the electrostatic interaction is the source of the discrimination in RNA polymerase. Proteins 2016; 84:1616-1624. © 2016 Wiley Periodicals, Inc.

Exploring the mechanism of DNA polymerases by analyzing the effect of mutations of active site acidic groups in Polymerase ß.

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase ß as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc.