Nuclear magnetic resonance study of the state of protonation of inhibitors bound to mutant dihydrofolate reductase lacking the active-site carboxyl

Kinetic Barrier to Enzyme Inhibition Is Manipulated by Dynamical Local Interactions in E. coli DHFR

Dihydrofolate reductase (DHFR) is an important drug target and a highly studied model protein for understanding enzyme dynamics. DHFR's crucial role in folate synthesis renders it an ideal candidate to understand protein function and protein evolution mechanisms. In this study, to understand how a newly proposed DHFR inhibitor, 4'-deoxy methyl trimethoprim (4'-DTMP), alters evolutionary trajectories, we studied interactions that lead to its superior performance over that of trimethoprim (TMP). To elucidate the inhibition mechanism of 4'-DTMP, we first confirmed, both computationally and experimentally, that the relative binding free energy cost for the mutation of TMP and 4'-DTMP is the same, pointing the origin of the characteristic differences to be kinetic rather than thermodynamic. We then employed an interaction-based analysis by focusing first on the active site and then on the whole enzyme. We confirmed that the polar modification in 4'-DTMP induces additional local interactions with the enzyme, particularly, the M20 loop. These changes are propagated to the whole enzyme as shifts in the hydrogen bond networks. To shed light on the allosteric interactions, we support our analysis with network-based community analysis and show that segmentation of the loop domain of inhibitor-bound DHFR must be avoided by a successful inhibitor.

Increased substrate affinity in the Escherichia coli L28R dihydrofolate reductase mutant causes trimethoprim resistance

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme with an essential role in cell metabolism. DHFR catalyzes the reduction of dihydrofolate to tetrahydrofolate, which is a precursor for purine and thymidylate synthesis. Several DHFR targeting antifolate drugs including trimethoprim, a competitive antibacterial inhibitor, have therefore been developed and are clinically used. Evolution of resistance against antifolates is a common public health problem rendering these drugs ineffective. To combat the resistance problem, it is important to understand resistance-conferring changes in the DHFR structure and accordingly develop alternative strategies. Here, we structurally and dynamically characterize Escherichia coli DHFR in its wild type (WT) and trimethoprim resistant L28R mutant forms in the presence of the substrate and its inhibitor trimethoprim. We use molecular dynamics simulations to determine the conformational space, loop dynamics and hydrogen bond distributions at the active site of DHFR for the WT and the L28R mutant. We also report their experimental kcat, Km, and Ki values, accompanied by isothermal titration calorimetry measurements of DHFR that distinguish enthalpic and entropic contributions to trimethoprim binding. Although mutations that confer resistance to competitive inhibitors typically make enzymes more promiscuous and decrease affinity to both the substrate and the inhibitor, strikingly, we find that the L28R mutant has a unique resistance mechanism. While the binding affinity differences between the WT and the mutant for the inhibitor and the substrate are small, the newly formed extra hydrogen bonds with the aminobenzoyl glutamate tail of DHF in the L28R mutant leads to increased barriers for the dissociation of the substrate and the product. Therefore, the L28R mutant indirectly gains resistance by enjoying prolonged binding times in the enzyme-substrate complex. While this also leads to slower product release and decreases the catalytic rate of the L28R mutant, the overall effect is the maintenance of a sufficient product formation rate. Finally, the experimental and computational analyses together reveal the changes that occur in the energetic landscape of DHFR upon the resistance-conferring L28R mutation. We show that the negative entropy associated with the binding of trimethoprim in WT DHFR is due to water organization at the binding interface. Our study lays the framework to study structural changes in other trimethoprim resistant DHFR mutants.

Development, validation, and application of adapted PEOE charges to estimate pKa values of functional groups in protein-ligand complexes

For routine pK(a) calculations of protein-ligand complexes in drug design, the PEOE method to compute partial charges was modified. The new method is applicable to a large scope of proteins and ligands. The adapted charges were parameterized using experimental free energies of solvation of amino acids and small organic ligands. For a data set of 80 small organic molecules, a correlation coefficient of r(2) = 0.78 between calculated and experimental solvation free energies was obtained. Continuum electrostatics pK(a) calculations based on the Poisson-Boltzmann equation were carried out on a validation set of nine proteins for which 132 experimental pK(a) values are known. In total, an overall RMSD of 0.88 log units between calculated and experimentally determined data is achieved. In particular, the predictions of significantly shifted pK(a) values are satisfactory, and reasonable estimates of protonation states in the active sites of lysozyme and xylanase could be obtained. Application of the charge-assignment and pK(a)-calculation procedure to protein-ligand complexes provides clear structural interpretations of experimentally observed changes of protonation states of functional groups upon complex formation. This information is essential for the interpretation of thermodynamic data of protein-ligand complex formation and provides the basis for the reliable factorization of the free energy of binding in enthalpic and entropic contributions. The modified charge-assignment procedure forms the basis for future automated pK(a) calculations of protein-ligand complexes.

Neutron diffraction studies of Escherichia coli dihydrofolate reductase complexed with methotrexate

Hydrogen atoms play a central role in many biochemical processes yet are difficult to visualize by x-ray crystallography. Spallation neutron sources provide a new arena for protein crystallography with TOF measurements enhancing data collection efficiency and allowing hydrogen atoms to be located in smaller crystals of larger biological macromolecules. Here we report a 2.2-A resolution neutron structure of Escherichia coli dihydrofolate reductase (DHFR) in complex with methotrexate (MTX). Neutron data were collected on a 0.3-mm(3) D(2)O-soaked crystal at the Los Alamos Neutron Scattering Center. This study provides an example of using spallation neutrons to study protein dynamics, to identify protonation states directly from nuclear density maps, and to analyze solvent structure. Our structure reveals that the occluded loop conformation [monomer (mon.) A] of the DHFR.MTX complex undergoes greater H/D exchange compared with the closed-loop conformer (mon. B), partly because the Met-20 and beta(F-G) loops readily exchange in mon. A. The eight-stranded beta sheet of both DHFR molecules resists H/D exchange more than the helices and loops. However, the C-terminal strand, betaH, in mon. A is almost fully exchanged. Several D(2)Os form hydrogen bonds with exchanged amides. At the active site, the N1 atom of MTX is protonated and thus charged when bound to DHFR. Several D(2)Os are observed at hydrophobic surfaces, including two pockets near the MTX-binding site. A previously unidentified D(2)O hydrogen bonds with the catalytic D27 in mon. B, stabilizing its negative charge.

Lipophilic conjugates of methotrexate with short-chain alkylamino acids as DHFR inhibitors. Synthesis, biological evaluation, and molecular modeling

Pursuing previous researches on lipophilic conjugates of methotrexate, aimed at over-crossing a form of transport resistance shown by some tumor cell lines toward the drug, a new series of derivatives is described in which the drug alpha- and gamma-carboxyl groups have been linked through amide bonds to short-chain alpha-alkylamino acids (4-6 carbon atoms). A specific NMR study was performed to delineate the stereochemistry of the conjugates. The inhibitory activity of these compounds against the target enzyme, (bovine liver) dihydrofolate reductase, and a sensitive (CCRF-CEM) and a transport-resistant tumor cell subline (CEM-MTX) were assessed. The conjugates showed the ability of retaining the same inhibitory activity also against the resistant cell subline, against which the parent drug was much less active than against the wild one; the alpha,gamma-bis(hexyl) derivative was the most active term of the series. Docking studies are in agreement with the proposed mode of interaction of these conjugates with the human DHFR.

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.

Dynamics of trimethoprim bound to dihydrofolate reductase--a deuterium NMR study

We have employed deuterium NMR techniques to determine the dynamics of trimethoprim (TMP) in a binary complex with dihydrofolate reductase (DHFR) or in a ternary complex with DHFR and cofactor NADP+ in the fully hydrated state. TMP was deuterated at the following positions: (2',6'-D2)TMP, (3'-Ome-D3)TMP and (3',4'-Ome-D6)TMP. Dynamics of TMP were deduced from lineshape simulation and relaxation measurements of the deuterium NMR powder spectra of the three samples obtained at various temperatures. The results showed that in the polycrystalline state the TMP molecule is very rigid. The only detectable motion is the methyl group rotation at a rate of 10(10) s-1 at 25 degrees C, as determined from simulation of the partially relaxed powder patterns. When bound to DHFR a residual deuterium quadrupole splitting of 140 kHz was observed for (2',6'-D2)TMP at temperatures up to 30 degrees C, suggesting that the benzyl ring in the bound state is also very rigid. In contrast, in the binary complex with DHFR the methoxyl groups of TMP undergo librational motion of 10(7) s-1 about the C3-O bond at an amplitude of 54 degrees for the meta methoxyl group and about the C4-O bond at an amplitude of 70 degrees and similar rate for the para methoxyl group at 30 degrees C. The presence of the cofactor, NADP+, appears to tighten up the binding pocket such that the motion freedom of TMP is more restricted. The rigidity of TMP in a protein complex as revealed by our deuterium NMR results is in accord with the tight binding of TMP to DHFR.

The conformations of trimethoprim/E. coli dihydrofolate reductase complexes. A 15N and 31P NMR study

We have employed 15N and 31P NMR techniques to characterize the conformations of trimethoprim (TMP)/E. coli dihydrofolate reductase (DHFR) complexes in the presence and absence of NADPH and NADP+. A single conformation was observed for TMP/DHFR, NADP+/DHFR, NADPH/DHFR, and TMP/NADPH/DHFR complexes. In the ternary complex of TMP/NADP+/DHFR both the 15N and 31P spectra revealed the presence of two conformations. However, the conformations of TMP and NADP+ in the ternary complex may not be correlated, resulting in the possible existence of four conformations for the protein ternary complex.

Crystal structures of recombinant human dihydrofolate reductase complexed with folate and 5-deazafolate

The 2.3-A crystal structure of recombinant human dihydrofolate reductase (EC 1.5.1.3, DHFR) has been solved as a binary complex with folate (a poor substrate at neutral pH) and also as a binary complex with an inhibitor, 5-deazafolate. The inhibitor appears to be protonated at N8 on binding, whereas folate is not. Rotation of the peptide plane joining I7 and V8 from its position in the folate complex permits hydrogen bonding of 5-deazafolate's protonated N8 to the backbone carbonyl of I7, thus contributing to the enzyme's greater affinity for 5-deazafolate than for folate. In this respect it is likely that bound 5-deazafolate furnishes a model for 7,8-dihydrofolate binding and, in addition, resembles the transition state for folate reduction. A hypothetical transition-state model for folate reduction, generated by superposition of the DHFR binary complexes human.5-deazafolate and chicken liver.NADPH, reveals a 1-A overlap of the binding sites for folate's pteridine ring and the dihydronicotinamide ring of NADPH. It is proposed that this binding-site overlap accelerates the reduction of both folate and 7,8-dihydrofolate by simultaneously binding substrate and cofactor with a sub van der Waals separation that is optimal for hydride transfer.

NMR studies of interactions of ligands with dihydrofolate reductase

NMR spectroscopy is a useful technique for studying interactions, conformations and dynamic processes within ligand-protein complexes. Several examples of the application of the method to studies of complexes of anti-folate drugs with their target enzyme, dihydrofolate reductase, are discussed.

Probing the salt bridge in the dihydrofolate reductase-methotrexate complex by using the coordinate-coupled free-energy perturbation method

The importance of the ionic interaction due to the formation of the salt bridge between the Asp-27 and the pteridine ring in Escherichia coli dihydrofolate reductase-methotrexate complex has been studied by using the free-energy perturbation method. The calculation suggests that the ion-pair contribution to the binding energy is insignificant, as the enzyme surroundings do not stabilize the salt bridge to the extent of the desolvation of the charged groups. The activation barrier for the proton exchange between the pteridine ring and the Asp-27 is calculated to be 20.1 kcal/mol (1 cal = 4.184 J) by using the coordinate-coupled perturbation method, implying that this may be a channel to the proton exchange from the pteridine ring to the solvent. The Gibbs-energy difference of binding between the Asn-27 and Ser-27 is calculated to be 3.2 kcal/mol and is mainly due to the electrostatic interactions.