Ab Initio QM/MM Study Shows a Highly Dissociated SN2 Hydrolysis Mechanism for the cGMP-Specific Phosphodiesterase-5

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Dynamical Basis of Allosteric Activation for the Plasmodium falciparum Protein Kinase G

The Plasmodium falciparum cGMP-dependent protein kinase (PfPKG) is required for the progression of the Plasmodium's life cycle and is therefore a promising malaria drug target. PfPKG includes four cGMP-binding domains (CBD-A to CBD-D). CBD-D plays a crucial role in PfPKG regulation as it is the primary determinant for the inhibition and cGMP-dependent activation of the catalytic domain. Hence, it is critical to understand how CBD-D is allosterically regulated by cGMP. Although the apo versus holo conformational changes of CBD-D have been reported, information on the intermediates of the activation pathway is currently lacking. Here, we employed molecular dynamics simulations to model four key states along the thermodynamic cycle for the cGMP-dependent activation of the PfPKG CBD-D domain. The simulations were compared to NMR data, and they revealed that the PfPKG CBD-D activation pathway samples a compact intermediate in which the N- and C-terminal helices approach the central β-barrel. In addition, by comparing the cGMP-bound active and inactive states, the essential binding interactions that differentiate these states were identified. The identification of structural and dynamical features unique to the cGMP-bound inactive state provides a promising basis to design PfPKG-selective allosteric inhibitors as a viable treatment for malaria.

Discovery and Optimization of Chromone Derivatives as Novel Selective Phosphodiesterase 10 Inhibitors

Phosphodiesterase 10 (PDE10) inhibitors have received much attention as promising therapeutic agents for central nervous system (CNS) disorders such as schizophrenia and Huntington's disease. Recently, a hit compound 1 with a novel chromone scaffold has shown moderate inhibitory activity against PDE10A (IC₅₀ = 500 nM). Hit-to-lead optimization has resulted in compound 3e with an improved inhibitory activity (IC₅₀ = 6.5 nM), remarkable selectivity (>95-fold over other PDEs), and good metabolic stability (RLM t_1/2 = 105 min) by using an integrated strategy (molecular modeling, chemical synthesis, bioassay, and cocrystal structure). The cocrystal structural information provides insights into the binding pattern of 3e in the PDE10A catalytic domain to highlight the key role of the halogen and hydrogen bonds toward Tyr524 and Tyr693, respectively, thereby resulting in high selectivity against other PDEs. These new observations are of benefit for the rational design of the next generation PDE10 inhibitors for CNS disorders.

Catalytic Roles of Coenzyme Pyridoxal-5'-phosphate (PLP) in PLP-dependent Enzymes: Reaction Pathway for Methionine-γ-lyase-catalyzed L-methionine Depletion

Pyridoxal-5'-phosphate (PLP), the active form of vitamin B₆, is an important and versatile coenzyme involved in a variety of enzymatic reactions, accounting for about 4% of all classified activities. However, the detailed catalytic reaction pathways for PLP-dependent enzymes remain to be explored. Methionine-γ-lyase (MGL), a promising alternative anti-tumor agent to conventional chemotherapies whose catalytic mechanism is highly desired for guiding further development of re-engineered enzymes, was used as a representative PLP-dependent enzyme, and the catalytic mechanism for L-Met elimination by MGL was explored at the first-principles quantum mechanical/molecular mechanical (QM/MM) level with umbrella sampling. The QM/MM calculations revealed that the enzymatic reaction pathway consists of 4 stages for a total of 19 reaction steps with five intermediates captured in available crystal structures. Furthermore, the more comprehensive role of PLP was revealed. Besides the commonly known role of "electron sink", coenzyme PLP can also assist proton transfer and temporarily store the excess proton generated in some intermediate states by using its hydroxyl group and phosphate group. Thus, PLP is participated in most of the 19 steps. This study not only provided a theoretical basis for further development and re-engineering MGL as a potential anti-tumor agent, but also revealed the comprehensive role of PLP which could be used to explore the mechanisms of other PLP-dependent enzymes.

Analysis of proton wires in the enzyme active site suggests a mechanism of c-di-GMP hydrolysis by the EAL domain phosphodiesterases

We report for the first time a hydrolysis mechanism of the cyclic dimeric guanosine monophosphate (c-di-GMP) by the EAL domain phosphodiesterases as revealed by molecular simulations. A model system for the enzyme-substrate complex was prepared on the base of the crystal structure of the EAL domain from the BlrP1 protein complexed with c-di-GMP. The nucleophilic hydroxide generated from the bridging water molecule appeared in a favorable position for attack on the phosphorus atom of c-di-GMP. The most difficult task was to find a pathway for a proton transfer to the O3' atom of c-di-GMP to promote the O3'P bond cleavage. We show that the hydrogen bond network extended over the chain of water molecules in the enzyme active site and the Glu359 and Asp303 side chains provides the relevant proton wires. The suggested mechanism is consistent with the structural, mutagenesis, and kinetic experimental studies on the EAL domain phosphodiesterases. Proteins 2016; 84:1670-1680. © 2016 Wiley Periodicals, Inc.