Theoretical studies on pteridines. 2. Geometries, tautomer, ionization and reduction energies of substrates and inhibitors of dihydrofolate reductase

Energetically most likely substrate and active-site protonation sites and pathways in the catalytic mechanism of dihydrofolate reductase

Despite much experimental and computational study, key aspects of the mechanism of reduction of dihydrofolate (DHF) by dihydrofolate reductase (DHFR) remain unresolved, while the secondary DHFR-catalyzed reduction of folate has been little studied. Major differences between proposed DHF mechanisms are whether the carboxylate group of the conserved active-site Asp or Glu residue is protonated or ionized during the reaction, and whether there is direct protonation of N5 or a proton shuttle from an initially protonated carboxylate group via O4. We have addressed these questions for both reduction steps with a comprehensive set of ab initio quantum chemical calculations on active-site fragment complexes, including the carboxyl side chain and, progressively, all other polar active-site residue groups including conserved water molecules. Addition of two protons in two steps was considered. The polarization effects of the remainder of the enzyme system were approximated by a dielectric continuum self-consistent reaction field (SCRF) model using an effective dielectric constant (epsilon) of 2. Optimized geometries were calculated using the density functional (B3LYP) method and Onsager SCRF model with the 6-31G basis. Single-point energy calculations were then carried out at the B3LYP/6-311+G level with either the Onsager or dielectric polarizable continuum model. Additional checking calculations at MP2 and HF levels, or with other basis sets or values of epsilon, were also done. From the results, the conserved water molecule, corresponding to W206 in the E. coli DHFR complexes, that is H-bonded to both the OD2 oxygen atom of the carboxyl (Asp) side chain and O4 of the pterin/dihydropterin ring, appears critically important and may determine the protonation site for the enzyme-bound substrates. In the absence of W206, the most stable monoprotonated species are the neutral-pair 4-enol forms of substrates with the carboxyl group OD2 oxygen protonated and H-bonded to N3. If W206 is included, then the most stable forms are still the neutral-pair complexes but now for the N3-H keto forms with the protonated OD2 atom H-bonding with W206. A second proton addition to these complexes gives protonations at N8 (folate) or N5 (DHF). Calculated H-bond distances correlate well with those for the conserved W206 observed in many X-ray structures. For all structures with occluded M20 loop conformations (closed active site), OD2-N3 distances are less than OD2-NA2 distances, which is consistent with those calculated for protonated OD2 complexes. Thus, the results (B3LYP; epsilon = 2 calculations) support a mechanism for both folate and DHF reduction in which the OD2 carboxyl oxygen is first protonated, followed by a direct protonation at N8 (folate) and N5 (DHF) to obtain the active cation complexes, i.e., doubly protonated. The results do not support a proposed protonated carboxyl with DHF in the enol form for the Michaelis complex, nor an ionized carboxyl with protonated enol-DHF as a catalytic intermediate. However, as additional calculations for the monoprotonated complete complexes show a reduction in the energy differences between the neutral-pair keto and ion-pair keto (N8- or N5-protonated) forms, we are extending the treatment using combined quantum mechanics and molecular mechanics (QM/MM) and molecular dynamics simulation methods to refine the description of the protein/solvent environment and prediction of the relative stabilization free energies of the various (OD2, O4, N5, and N8) protonation sites.

Computer-aided drug design: a free energy perturbation study on the binding of methyl-substituted pterins and N5-deazapterins to dihydrofolate reductase

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of 8-methylpterins and 8-methyl-N5-deazapterins to dihydrofolate reductase (DHFR). Methyl-substitution at the 5, 6 and 7 positions in the N-heterocyclic ring gives rise to a variety of ring substituent patterns and biological activity: several of these methyl derivatives of the 8-methyl parent compounds (8-methylpterin and 8-methyl-N5-deazapterin) have been identified as substrates or inhibitors of vertebrate DHFR in previous work. The calculated free energy differences reveal that the methyl-substituted compounds are thermodynamically more stable than the primary compounds (8-methylpterin and 8-methyl-N5-deazapterin) when bound to the enzyme, due largely to hydrophobic hydration phenomena. Methyl substitution at the 5 and/or 7 positions in the 6-methyl-substituted compounds has only a small effect on the stability of ligand binding. Furthermore, repulsive interactions between the 6-methyl substituent and DHFR are minimal, suggesting that the 6-methyl position is optimal for binding. The results also show that similarly substituted 8-methylpterins and 8-methyl-N5-deazapterins have very similar affinities for binding to DHFR. The computer simulation predictions are in broad agreement with experimental data obtained from kinetic studies, i.e. 6,8-dimethylpterin is a more efficient substrate than 8-methylpterin and 6,8-dimethyl-N5-deazapterin is a better inhibitor than 8-methyl-N5-deazapterin.

Computer-aided design of mechanism-based pterin analogues and MD/FEP simulations of their binding to dihydrofolate reductase

Using a combined theoretical and experimental approach we have been able to predict several chemical properties and the contributions of the many factors which determine the macroscopic binding behaviour of these new mechanism-based compounds with DHFR, and also analyse experimental data to develop structure-activity relationships.

Computer studies on the stereostructure and quantum chemical properties of 6-pyruvoyl tetrahydropterin, the key intermediate of tetrahydrobiopterin biosynthesis

The optimized geometry of the conformation of atoms constituting the 6-pyruvoyl tetrahydropterin molecule, the labile key intermediate of tetrahydrobiopterin biosynthesis, was obtained by molecular orbital calculations within the MINDO/3 framework. The stereostructure of the molecule showing the preferred mode for binding to sepiapterin reductase or pyruvoyl tetrahydropterin reductase was drawn in perspective. The resulting structure with the equatorial staggered configuration of the 6-1',2'-dioxopropyl (pyruvoyl) side chain indicated that O(1') and H(6) were located in the trans position around the C(6)-C(1') bond and that the two vicinal carbonyls in the side chain were fixed in the incomplete trans form. The calculation of atomic charges and LUMO coefficients of these carbonyls suggests that the C2'-carbonyl may be more reactive toward NADPH than the C1'-carbonyl in the enzymatic reaction.

Computational studies on pterins and speculations on the mechanism of action of dihydrofolate reductase

The significance of the enol form of the pterin ring in enzymatic reduction of dihydrofolate by DHFR is discussed on the basis of the results of ab initio calculations carried out on the keto/enol tautomers of 6-methyl-7, 8-dihydropterin as the model compound for the natural substrate, dihydrofolate.