Consideration of the pH-dependent inhibition of dihydrofolate reductase by methotrexate.

Poisson-Boltzmann calculations were used to determine the pKa of protein functional groups in the unliganded dihydrofolate reductase enzyme, and the pKa of protein and ligand groups in methotrexate-enzyme complexes. The results reported here are in conflict with two fundamental tenets of dihydrofolate reductase inhibition by methotrexate: (1) Asp27 is not expected to be protonated near pH 6.5 in the apoenzyme as previously proposed based on fitting of empirical equations to binding data, and (2) the calculated pKa for the pteridine N1 of the inhibitor while bound to the protein is significantly lower than that estimated for this group from interpretation of NMR data (>10). In fact, the electrostatic calculations and complementary quantum chemical calculations indicate that Asp27 is likely protonated when methotrexate is bound, resulting in a neutral dipole-dipole interaction rather than a salt-bridge between the enzyme and the inhibitor. Reasons for this discrepancy with the experimental data are discussed. Furthermore, His45 and Glu17 in the Escherichia coli enzyme are proposed to be in part responsible for the pH dependence of the conformational degeneracy in the inhibitor-enzyme complex.

A perspective on biological catalysis.

We have analysed enzyme catalysis through a re-examination of the reaction coordinate. The ground state of the enzyme-substrate complex is shown to be related to the transition state through the mean force acting along the reaction path; as such, catalytic strategies cannot be resolved into ground state destabilization versus transition state stabilization. We compare the role of active-site residues in the chemical step with the analogous role played by solvent molecules in the environment of the noncatalysed reaction. We conclude that enzyme catalysis is significantly enhanced by the ability of the enzyme to preorganize the reaction environment. This complementation of the enzyme to the substrate's transition state geometry acts to eliminate the slow components of solvent reorganization required for reactions in aqueous solution. Dramatically strong binding of the transition state geometry is not required.

Conservative and nonconservative mutations in proteins: anomalous mutations in a transport receptor analyzed by free energy and quantum chemical calculations.

Experimental studies on a bacterial sulfate receptor have indicated anomalous relative binding affinities for the mutations Ser130-->Cys,Ser130-->Gly, and Ser130-->Ala. The loss of affinity for sulfate in the former mutation was previously attributed to a greater steric effect on the part of the Cys side chain relative to the Ser side chain, whereas the relatively small loss of binding affinity for the latter two mutations was attributed to the loss of a single hydrogen bond. In this report we present quantum chemical and statistical thermodynamic studies of these mutations. Qualitative results from these studies indicate that for the Ser130-->Cys mutation the large decrease in binding affinity is in part caused by steric effects, but also significantly by the differential work required to polarize the Cys thiol group relative to the Ser hydroxyl group. The Gly mutant cobinds a water molecule in the same location as the Ser side chain resulting in a relatively small decrease in binding affinity. Results for the Ala mutant are in disagreement with experimental results but are likely to be limited by insufficient sampling of configuration space due to physical constraints applied during the simulation.