1
|
Abbat S, Jaladanki CK, Bharatam PV. Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis. J Mol Graph Model 2018; 87:76-88. [PMID: 30508692 DOI: 10.1016/j.jmgm.2018.11.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2018] [Revised: 11/16/2018] [Accepted: 11/19/2018] [Indexed: 11/18/2022]
Abstract
The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP+ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP+ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.
Collapse
Affiliation(s)
- Sheenu Abbat
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India
| | - Chaitanya K Jaladanki
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India
| | - Prasad V Bharatam
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India; Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India.
| |
Collapse
|
2
|
Doron D, Major DT, Kohen A, Thiel W, Wu X. Hybrid Quantum and Classical Simulations of the Dihydrofolate Reductase Catalyzed Hydride Transfer Reaction on an Accurate Semi-Empirical Potential Energy Surface. J Chem Theory Comput 2011; 7:3420-37. [PMID: 26598171 DOI: 10.1021/ct2004808] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dihydrofolate reductase (DHFR) catalyzes the reduction of 7,8-dihydrofolate by nicotinamide adenine dinucleotide phosphate hydride (NADPH) to form 5,6,7,8-tetrahydrofolate and oxidized nicotinamide. DHFR is a small, flexible, monomeric protein with no metals or SS bonds and serves as one of the enzymes commonly used to examine basic aspects in enzymology. In the current work, we present extensive benchmark calculations for several model reactions in the gas phase that are relevant to the DHFR catalyzed hydride transfer. To this end, we employ G4MP2 and CBS-QB3 ab initio calculations as well as numerous density functional theory methods. Using these results, we develop two specific reaction parameter (SRP) Hamiltonians based on the semiempirical AM1 method. The first generation SRP Hamiltonian does not account for dispersion, while the second generation SRP accounts for dispersion implicitly via the AM1 core-repulsion functions. These SRP semiempirical Hamiltonians are subsequently used in hybrid quantum mechanics/molecular mechanics simulations of the DHFR catalyzed reaction. Finally, kinetic isotope effects are computed using a mass-perturbation-based path-integral approach.
Collapse
Affiliation(s)
- Dvir Doron
- Department of Chemistry, The Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Dan Thomas Major
- Department of Chemistry, The Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Amnon Kohen
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Xin Wu
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
3
|
Alonso H, Cummins PL, Gready JE. Methyltetrahydrofolate:corrinoid/iron−sulfur Protein Methyltransferase (MeTr): Protonation State of the Ligand and Active-Site Residues. J Phys Chem B 2009; 113:14787-96. [DOI: 10.1021/jp900181g] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hernán Alonso
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Peter L. Cummins
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Jill E. Gready
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| |
Collapse
|
4
|
Kumarasiri M, Baker GA, Soudackov AV, Hammes-Schiffer S. Computational approach for ranking mutant enzymes according to catalytic reaction rates. J Phys Chem B 2009; 113:3579-83. [PMID: 19235997 DOI: 10.1021/jp810363k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
A computationally efficient approach for ranking mutant enzymes according to the catalytic reaction rates is presented. This procedure requires the generation and equilibration of the mutant structures, followed by the calculation of partial free energy curves using an empirical valence bond potential in conjunction with biased molecular dynamics simulations and umbrella integration. The individual steps are automated and optimized for computational efficiency. This approach is used to rank a series of 15 dihydrofolate reductase mutants according to the hydride transfer reaction rate. The agreement between the calculated and experimental changes in the free energy barrier upon mutation is encouraging. The computational approach predicts the correct direction of the change in free energy barrier for all mutants, and the correlation coefficient between the calculated and experimental data is 0.82. This general approach for ranking protein designs has implications for protein engineering and drug design.
Collapse
Affiliation(s)
- Malika Kumarasiri
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | | | | | |
Collapse
|
5
|
Khavrutskii IV, Price DJ, Lee J, Brooks CL. Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate. Protein Sci 2007; 16:1087-100. [PMID: 17473015 PMCID: PMC2206655 DOI: 10.1110/ps.062724307] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2006] [Revised: 03/05/2007] [Accepted: 03/06/2007] [Indexed: 10/23/2022]
Abstract
We evaluate the pK(a) of dihydrofolate (H(2)F) at the N(5) position in three ternary complexes with Escherichia coli dihydrofolate reductase (ecDHFR), namely ecDHFR(NADP(+):H(2)F) in the closed form (1), and the Michaelis complexes ecDHFR(NADPH:H(2)F) in the closed (2) and occluded (3) forms, by performing free energy perturbation with molecular dynamics simulations (FEP/MD). Our simulations suggest that in the Michaelis complex the pK(a) is modulated by the Met20 loop fluctuations, providing the largest pK(a) shift in substates with a "tightly closed" loop conformation; in the "partially closed/open" substates, the pK(a) is similar to that in the occluded complex. Conducive to the protonation, tightly closing the Met20 loop enhances the interactions of the cofactor and the substrate with the Met20 side chain and aligns the nicotinamide ring of the cofactor coplanar with the pterin ring of the substrate. Overall, the present study favors the hypothesis that N(5) is protonated directly from solution and provides further insights into the mechanism of the substrate protonation.
Collapse
Affiliation(s)
- Ilja V Khavrutskii
- The Scripps Research Institute, Department of Molecular Biology, TPC6, La Jolla, California 92037, USA
| | | | | | | |
Collapse
|
6
|
Cummins PL, Rostov IV, Gready JE. Calculation of a Complete Enzymic Reaction Surface: Reaction and Activation Free Energies for Hydride-Ion Transfer in Dihydrofolate Reductase. J Chem Theory Comput 2007; 3:1203-11. [DOI: 10.1021/ct600313b] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Peter L. Cummins
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Ivan V. Rostov
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Jill E. Gready
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| |
Collapse
|
7
|
Cavalli A, Carloni P, Recanatini M. Target-Related Applications of First Principles Quantum Chemical Methods in Drug Design. Chem Rev 2006; 106:3497-519. [PMID: 16967914 DOI: 10.1021/cr050579p] [Citation(s) in RCA: 96] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Andrea Cavalli
- Department of Pharmaceutical Sciences, University of Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | | | | |
Collapse
|
8
|
Gao J, Ma S, Major DT, Nam K, Pu J, Truhlar DG. Mechanisms and free energies of enzymatic reactions. Chem Rev 2006; 106:3188-209. [PMID: 16895324 PMCID: PMC4477011 DOI: 10.1021/cr050293k] [Citation(s) in RCA: 317] [Impact Index Per Article: 17.6] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jiali Gao
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Shuhua Ma
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Dan T. Major
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Kwangho Nam
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Jingzhi Pu
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, Digital Technology Center, University of Minnesota, Minneapolis, Minnesota 55455
| |
Collapse
|
9
|
Abstract
This review examines the linkage between protein conformational motions and enzyme catalysis. The fundamental issues related to this linkage are probed in the context of two enzymes that catalyze hydride transfer, namely dihydrofolate reductase and liver alcohol dehydrogenase. The extensive experimental and theoretical studies addressing the role of protein conformational changes in these enzyme reactions are summarized. Evidence is presented for a network of coupled motions throughout the protein fold that facilitate the chemical reaction. This network is comprised of fast thermal motions that are in equilibrium as the reaction progresses along the reaction coordinate and that lead to slower equilibrium conformational changes conducive to the chemical reaction.
Collapse
Affiliation(s)
- Sharon Hammes-Schiffer
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, USA.
| | | |
Collapse
|
10
|
Abstract
Theoretical perspectives on hydrogen transfer reactions in enzymes are presented. The proton-coupled electron transfer reaction catalyzed by soybean lipoxygenase and the hydride transfer reaction catalyzed by dihydrofolate reductase are discussed. The first reaction is nonadiabatic and involves two distinct electronic states, while the second reaction is predominantly adiabatic and occurs on the electronic ground state. Theoretical studies indicate that hydrogen tunneling and protein motion play significant roles in both reactions. In both cases, the proton donor-acceptor distance decreases relative to its equilibrium value to enable efficient hydrogen tunneling. Equilibrium thermal motions of the protein lead to conformational changes that facilitate hydrogen transfer, but the nonequilibrium dynamical aspects of these motions have negligible impact.
Collapse
Affiliation(s)
- Sharon Hammes-Schiffer
- Department of Chemistry, 104 Chemistry Building, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
11
|
Hammes-Schiffer S. Quantum-classical simulation methods for hydrogen transfer in enzymes: a case study of dihydrofolate reductase. Curr Opin Struct Biol 2005; 14:192-201. [PMID: 15093834 DOI: 10.1016/j.sbi.2004.03.008] [Citation(s) in RCA: 56] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
A variety of theoretical approaches have been used to investigate hydrogen transfer in enzymatic reactions. The free energy barriers for hydrogen transfer in enzymes have been calculated using classical molecular dynamics simulations in conjunction with quantum mechanical/molecular mechanical and empirical valence bond potentials. Nuclear quantum effects have been included with several different approaches. Applications of these approaches to hydride transfer in dihydrofolate reductase are consistent with experimental measurements and provide significant insight into the protein conformational changes that facilitate the hydride transfer reaction.
Collapse
Affiliation(s)
- Sharon Hammes-Schiffer
- Department of Chemistry, 152 Davey Laboratory, Pennsylvania State University, University Park, Pennsylvania 16802, USA.
| |
Collapse
|
12
|
Cummins PL, Gready JE. Computational methods for the study of enzymic reaction mechanisms III: A perturbation plus QM/MM approach for calculating relative free energies of protonation. J Comput Chem 2005; 26:561-8. [PMID: 15726569 DOI: 10.1002/jcc.20192] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We describe a coupling parameter, that is, perturbation, approach to effectively create and annihilate atoms in the quantum mechanical Hamiltonian within the closed shell restricted Hartree-Fock formalism. This perturbed quantum mechanical atom (PQA) method is combined with molecular mechanics (MM) methods (PQA/MM) within a molecular dynamics simulation, to model the protein environment (MM region) effects that also make a contribution to the overall free energy change. Using the semiempirical PM3 method to model the QM region, the application of this PQA/MM method is illustrated by calculation of the relative protonation free energy of the conserved OD2 (Asp27) and the N5 (dihydrofolate) proton acceptor sites in the active site of Escherichia coli dihydrofolate reductase (DHFR) with the bound nicotinamide adenine dinucleotide phosphate (NADPH) cofactor. For a number of choices for the QM region, the relative protonation free energy was calculated as the sum of contributions from the QM region and the interaction between the QM and MM regions via the thermodynamic integration (TI) method. The results demonstrate the importance of including the whole substrate molecule in the QM region, and the overall protein (MM) environment in determining the relative stabilities of protonation sites in the enzyme active site. The PQA/MM free energies obtained by TI were also compared with those estimated by a less computationally demanding nonperturbative method based on the linear response approximation (LRA). For some choices of QM region, the total free energies calculated using the LRA method were in very close agreement with the PQA/MM values. However, the QM and QM/MM component free energies were found to differ significantly between the two methods.
Collapse
Affiliation(s)
- Peter L Cummins
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | | |
Collapse
|
13
|
Wong KF, Watney JB, Hammes-Schiffer S. Analysis of Electrostatics and Correlated Motions for Hydride Transfer in Dihydrofolate Reductase. J Phys Chem B 2004. [DOI: 10.1021/jp048565v] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Kim F. Wong
- Department of Chemistry, 152 Davey Laboratory, The Pennsylvania State University, University Park, Pennsylvania 16802
| | - James B. Watney
- Department of Chemistry, 152 Davey Laboratory, The Pennsylvania State University, University Park, Pennsylvania 16802
| | - Sharon Hammes-Schiffer
- Department of Chemistry, 152 Davey Laboratory, The Pennsylvania State University, University Park, Pennsylvania 16802
| |
Collapse
|
14
|
Thorpe IF, Brooks CL. The coupling of structural fluctuations to hydride transfer in dihydrofolate reductase. Proteins 2004; 57:444-57. [PMID: 15382243 DOI: 10.1002/prot.20219] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The energy barrier for hydride transfer in wild-type G121V and G121S variants of Escherichia coli dihydrofolate reductase (DHFR) fluctuates in a time-dependent manner. This fluctuation may be attributed to structural changes in the protein that modulate the site of chemistry. Despite being far from the active site, mutations at position 121 of DHFR reduce the hydride transfer rate of the enzyme. This occurrence has been suggested to arise from modifications to the conformational ensemble of the protein. We elucidate the effects of the G121S and G121V mutations on the hydride transfer barrier by identifying structural changes in the protein that correlate with lowered barriers. The effect of these structural parameters on the hydride transfer barrier may be rationalized by simple considerations of the geometric constraints of the hydride transfer reaction. Fluctuations of these properties are associated with specific backbone dihedral angles of residues within the Methione-20 (M20) loop. The dihedral angle preferences are mediated by interactions with the region of the enzyme in the vicinity of residue 121 and are translated into distinct ligand conformations. We predict mutations within the M20 loop that may alter the conformational space explored by DHFR. Such mutational changes are anticipated to adjust the hydride transfer efficacy of DHFR by modifying equilibrium distributions of hydride transfer barriers found in the enzyme.
Collapse
Affiliation(s)
- Ian F Thorpe
- Department of Molecular Biology (TPC-06), The Scripps Research Institute, La Jolla, California 92037, USA
| | | |
Collapse
|
15
|
Swanwick RS, Shrimpton PJ, Allemann RK. Pivotal Role of Gly 121 in Dihydrofolate Reductase from Escherichia coli: The Altered Structure of a Mutant Enzyme May Form the Basis of Its Diminished Catalytic Performance. Biochemistry 2004; 43:4119-27. [PMID: 15065854 DOI: 10.1021/bi036164k] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The structure and folding of dihydrofolate reductase (DHFR) from Escherichia coli and the mutant G121V-DHFR, in which glycine 121 in the exterior FG loop was replaced with valine, were studied by molecular dynamics simulations and CD and fluorescence spectroscopy. The importance of residue 121 for the chemical step during DHFR catalysis had been demonstrated previously. High-temperature MD simulations indicated that while DHFR and G121V-DHFR followed similar unfolding pathways, the strong contacts between the M20 loop and the FG loop in DHFR were less stable in the mutant. These contacts have been proposed to be involved in a coupled network of interactions that influence the protein dynamics and promote catalysis [Benkovic, S. J., and Hammes-Schiffer, S. (2003) Science 301, 1196-1202]. CD spectroscopy of DHFR and G121V-DHFR indicated that the two proteins existed in different conformations at room temperature. While the thermally induced unfolding of DHFR was highly cooperative with a midpoint at 51.6 +/- 0.7 degrees C, G121V-DHFR exhibited a gradual decrease in its level of secondary structure without a clear melting temperature. Temperature-induced unfolding and renaturation from the urea-denatured state revealed that both proteins folded via highly fluorescent intermediates. The formation of these intermediates occurred with relaxation times of 149 +/- 4.5 and 256 +/- 13 ms for DHFR and G121V-DHFR, respectively. The fluorescence intensity for the intermediates formed during refolding of G121V-DHFR was approximately twice that of the wild-type. While the fluorescence intensity then slowly decayed for DHFR toward a state representing the native protein, G121V-DHFR appeared to be trapped in a highly fluorescent state. These results suggest that the reduced catalytic activity of G121V-DHFR is the consequence of nonlocal structural effects that may result in a perturbation of the network of promoting motions.
Collapse
|
16
|
Ferrer S, Silla E, Tuñón I, Martí S, Moliner V. Catalytic Mechanism of Dihydrofolate Reductase Enzyme. A Combined Quantum-Mechanical/Molecular-Mechanical Characterization of the N5 Protonation Step. J Phys Chem B 2003. [DOI: 10.1021/jp0354898] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
17
|
Thorpe IF, Brooks CL. Barriers to Hydride Transfer in Wild Type and Mutant Dihydrofolate Reductase from E. coli. J Phys Chem B 2003. [DOI: 10.1021/jp035734n] [Citation(s) in RCA: 53] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ian F. Thorpe
- Department of Molecular Biology (TPC6), Center for Theoretical Biological Physics, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037
| | - Charles L. Brooks
- Department of Molecular Biology (TPC6), Center for Theoretical Biological Physics, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037
| |
Collapse
|
18
|
Kedzierski P, Sokalski W, Cheng H, Mitchell J, Leszczynski J. DFT study of the reaction proceeding in the cytidine deaminase. Chem Phys Lett 2003. [DOI: 10.1016/j.cplett.2003.10.042] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
19
|
Computational methods for the study of enzymic reaction mechanisms. II. An overlapping mechanically embedded method for hybrid semi-empirical-QM/MM calculations. ACTA ACUST UNITED AC 2003. [DOI: 10.1016/s0166-1280(03)00303-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
20
|
Abstract
Dihydrofolate Reductase (DHFR) catalyzes the reduction of dihydrofolate (H2F) to tetrahydrofolate. On the basis of 10-12.5 ns molecular dynamics simulations of two conformations (closed and occluded) of the ternary DHFR/NADPH/H2F complex from Escherichia coli and a free energy perturbation approach, we have calculated the pKa value for the N5 atom in H2F. Our results suggest that the N5 atom in H2F is responsible for the pH dependency of the catalyzed reaction, meaning that DHFR facilitates protonation of H2F by approximately 4 pKa units. The mechanism behind this increase is due to favorable electrostatic interactions between the Asp27 residue and a proton at the N5 atom. The electrostatic interactions are enhanced by a hydrophobic active site, which to a large extent is made hydrophobic by the M20 loop in DHFR. Moreover, we find that the conformation imposed on H2F by DHFR to some extent also favors protonation of the N5 atom. Our results add support to previous findings and suggestions by Callender and co-workers [e.g., Deng, J.; Callender, R. J. Am. Chem. Soc. 1998, 120, 7730-7737] and explain why mutation of Asp27 may lead to severely reduced activity at neutral pH.
Collapse
Affiliation(s)
- Thomas H Rod
- Department of Molecular Biology, TPC6, 10550 North Torrey Pines Road, La Jolla, California 92037, USA
| | | |
Collapse
|