Combined quantum and molecular mechanics (QM/MM) study of the ionization state of 8-methylpterin substrate bound to dihydrofolate reductase

Simple item page
dc.contributor.authorCummins, Peter
dc.contributor.authorGready, Jill
dc.date.accessioned2015-12-13T23:18:29Z
dc.date.issued2000
dc.date.updated2015-12-12T08:57:01Z
dc.description.abstractThe prediction of stabilities of ion pairs relative to the corresponding neutral pairs is a challenging computational problem but very important in biological systems where the occurrence of such pairs is often mechanistically significant. Here we have computed the relative free energy for the protonated and neutral forms of the mechanism-based substrate 8-methylpterin bound to dihydrofolate reductase (DHFR) in a ternary complex with cofactor nicotinamide adenine dinucleotide phosphate (NADPH). The free energy components required to calculate relative affinities of the protonated substrate hydrogen (H) bonded to the unprotonated form of the conserved active-site carboxylate residue (Glu30 for chicken DHFR) and neutral substrate H bonded to the neutral (carboxylic acid) form of Glu30 were obtained using both ab initio QM methods and combined semiempirical (AMI) QM and MM (QM/MM) methods. The energy difference for the H-bonded systems was first calculated in a vacuum using ab initio QM methods for a model system consisting of the Glu30 side chain and the substrate. The free energy of interaction of this system with the protein/solvent surroundings was then computed using a coordinate-coupled free energy perturbation (FEP) method implemented within the molecular dynamics (MD) simulation scheme. The ab initio QM results show the neutral form is more stable in a vacuum, but that this relative stability is reduced in the enzyme due to the more favorable electrostatic interactions made by the ion pair. The neutral pair is calculated to be approximately 5 ± 2 kcal/mol (ΔΔGbind) more stable than the ion pair, but the results are suggestive of a small free energy gap and the possibility of a low energy barrier for the proton transfer. The results are in qualitative agreement with other theoretical studies on similar types of H-bonded systems but appear to contradict our own experimental data for these DHFR ligands. Reasons for the apparent discrepancy are discussed.
dc.identifier.issn1520-6106
dc.identifier.urihttp://hdl.handle.net/1885/90201
dc.publisherAmerican Chemical Society
dc.sourceJournal of Physical Chemistry B
dc.titleCombined quantum and molecular mechanics (QM/MM) study of the ionization state of 8-methylpterin substrate bound to dihydrofolate reductase
dc.typeJournal article
local.bibliographicCitation.lastpage4510
local.bibliographicCitation.startpage4503
local.contributor.affiliationCummins, Peter, College of Medicine, Biology and Environment, ANU
local.contributor.affiliationGready, Jill, College of Medicine, Biology and Environment, ANU
local.contributor.authoremailu9508363@anu.edu.au
local.contributor.authoruidCummins, Peter, u9508363
local.contributor.authoruidGready, Jill, u9508375
local.description.embargo2037-12-31
local.description.notesImported from ARIES
local.description.refereedYes
local.identifier.absfor030799 - Theoretical and Computational Chemistry not elsewhere classified
local.identifier.absfor060107 - Enzymes
local.identifier.ariespublicationMigratedxPub20500
local.identifier.citationvolume104
local.identifier.scopusID2-s2.0-0001685622
local.identifier.uidSubmittedByMigrated
local.type.statusPublished Version

Downloads

Original bundle
Now showing 1 - 1 of 1
No Thumbnail Available
Name:
01_Cummins_Combined_quantum_and_molecular_2000.pdf
Size:
110.41 KB
Format:
Adobe Portable Document Format
Download

Collections

ANU Research Publications

DSpace software copyright © 2002-2024 LYRASIS