Interligand Overhauser effects in type II dihydrofolate reductase

NMR and computational studies reveal novel aspects in molecular recognition of unsaturated fatty acids with non-labelled serum albumin

An approach based on the combined use of saturation transfer difference (STD), Tr-NOESY and Inter-ligand NOEs for PHArmacophore Mapping (INPHARMA) NMR techniques and docking calculations is reported, for the first time, for mapping interactions and specific binding sites of caproleic acid (10 : 1 cis-9), oleic acid (18 : 1 cis-9), linoleic acid (18 : 2 cis-9,12) and linolenic (18 : 3, cis-9,12,15) free fatty acids (FFAs) with non-labelled serum albumin (BSA/HSA). Significant negative inter-ligand NOEs between the FFAs and the drugs ibuprofen and warfarin, through competition experiments, were observed. The inter-ligand NOEs and docking calculations were interpreted in terms of competitive binding mode, the significant folding of the bis allylic region and the presence of two orientations of the FFAs in the warfarin binding site (FA7), due to two potential distinctive anchoring polar groups of amino acids. This conformational flexibility is the reason that, the location and conformational states of the FFAs in the binding site of warfarin could not be determined accurately, despite numerous available X-ray structural studies. α-Linolenic acid competes favourably with warfarin at the binding site FA7. Isothermal titration calorimetry experiments of the preformed HSA/α-linolenic acid complex upon titration with warfarin show a significant reduction in the binding constant of warfarin, in very good agreement with NMR and computational data. The combined use, therefore, of STD, Tr-NOESY and INPHARMA NMR, ITC and docking calculations may find promising applications in the field of protein-lipid recognition research.

Computational Development of Inhibitors of Plasmid-Borne Bacterial Dihydrofolate Reductase

Resistance to trimethoprim and other antibiotics targeting dihydrofolate reductase may arise in bacteria harboring an atypical, plasmid-encoded, homotetrameric dihydrofolate reductase, called R67 DHFR. Although developing inhibitors to this enzyme may be expected to be promising drugs to fight trimethoprim-resistant strains, there is a paucity of reports describing the development of such molecules. In this manuscript, we describe the design of promising lead compounds to target R67 DHFR. Density-functional calculations were first used to identify the modifications of the pterin core that yielded derivatives likely to bind the enzyme and not susceptible to being acted upon by it. These unreactive molecules were then docked to the active site, and the stability of the docking poses of the best candidates was analyzed through triplicate molecular dynamics simulations, and compared to the binding stability of the enzyme–substrate complex. Molecule 32 ([6-(methoxymethyl)-4-oxo-3,7-dihydro-4H-pyrano[2,3-d]pyrimidin-2-yl]methyl-guanidinium) was shown by this methodology to afford extremely stable binding towards R67 DHFR and to prevent simultaneous binding to the substrate. Additional docking and molecular dynamics simulations further showed that this candidate also binds strongly to the canonical prokaryotic dihydrofolate reductase and to human DHFR, and is therefore likely to be useful to the development of chemotherapeutic agents and of dual-acting antibiotics that target the two types of bacterial dihydrofolate reductase.

Process of Fragment-Based Lead Discovery-A Perspective from NMR

Fragment-based lead discovery (FBLD) has proven fruitful during the past two decades for a variety of targets, even challenging protein–protein interaction (PPI) systems. Nuclear magnetic resonance (NMR) spectroscopy plays a vital role, from initial fragment-based screening to lead generation, because of its power to probe the intrinsically weak interactions between targets and low-molecular-weight fragments. Here, we review the NMR FBLD process from initial library construction to lead generation. We describe technical aspects regarding fragment library design, ligand- and protein-observed screening, and protein–ligand structure model generation. For weak binders, the initial hit-to-lead evolution can be guided by structural information retrieved from NMR spectroscopy, including chemical shift perturbation, transferred pseudocontact shifts, and paramagnetic relaxation enhancement. This perspective examines structure-guided optimization from weak fragment screening hits to potent leads for challenging PPI targets.

Tales of Dihydrofolate Binding to R67 Dihydrofolate Reductase

Homotetrameric R67 dihydrofolate reductase possesses 222 symmetry and a single active site pore. This situation results in a promiscuous binding site that accommodates either the substrate, dihydrofolate (DHF), or the cofactor, NADPH. NADPH interacts more directly with the protein as it is larger than the substrate. In contrast, the p-aminobenzoyl-glutamate tail of DHF, as monitored by nuclear magnetic resonance and crystallography, is disordered when bound. To explore whether smaller active site volumes (which should decrease the level of tail disorder by confinement effects) alter steady state rates, asymmetric mutations that decreased the half-pore volume by ∼35% were constructed. Only minor effects on k(cat) were observed. To continue exploring the role of tail disorder in catalysis, 1-ethyl-3-[3-(dimethylamino)propyl]carbodiimide-mediated cross-linking between R67 DHFR and folate was performed. A two-folate, one-tetramer complex results in the loss of enzyme activity where two symmetry-related K32 residues in the protein are cross-linked to the carboxylates of two bound folates. The tethered folate could be reduced, although with a ≤30-fold decreased rate, suggesting decreased dynamics and/or suboptimal positioning of the cross-linked folate for catalysis. Computer simulations that restrain the dihydrofolate tail near K32 indicate that cross-linking still allows movement of the p-aminobenzoyl ring, which allows the reaction to occur. Finally, a bis-ethylene-diamine-α,γ-amide folate adduct was synthesized; both negatively charged carboxylates in the glutamate tail were replaced with positively charged amines. The K(i) for this adduct was ∼9-fold higher than for folate. These various results indicate a balance between folate tail disorder, which helps the enzyme bind substrate while dynamics facilitates catalysis.

NMR-based analysis of protein-ligand interactions

Physiological processes are mainly controlled by intermolecular recognition mechanisms involving protein-protein and protein-ligand (low molecular weight molecules) interactions. One of the most important tools for probing these interactions is high-field solution nuclear magnetic resonance (NMR) through protein-observed and ligand-observed experiments, where the protein receptor or the organic compounds are selectively detected. NMR binding experiments rely on comparison of NMR parameters of the free and bound states of the molecules. Ligand-observed methods are not limited by the protein molecular size and therefore have great applicability for analysing protein-ligand interactions. The use of these NMR techniques has considerably expanded in recent years, both in chemical biology and in drug discovery. We review here three major ligand-observed NMR methods that depend on the nuclear Overhauser effect-transferred nuclear Overhauser effect spectroscopy, saturation transfer difference spectroscopy and water-ligand interactions observed via gradient spectroscopy experiments-with the aim of reporting recent developments and applications for the characterization of protein-ligand complexes, including affinity measurements and structural determination.

Further studies on the role of water in R67 dihydrofolate reductase

Previous osmotic pressure studies of two nonhomologous dihydrofolate reductase (DHFR) enzymes found tighter binding of the nicotinamide adenine dinucleotide phosphate cofactor upon addition of neutral osmolytes. This result is consistent with water release accompanying binding. In contrast, osmotic stress studies found weaker binding of the dihydrofolate (DHF) substrate for both type I and type II DHFRs in the presence of osmolytes; this observation can be explained if dihydrofolate interacts with osmolytes and shifts the equilibrium from the enzyme-bound state toward the unbound substrate. Nuclear magnetic resonance experiments support this hypothesis, finding that osmolytes interact with dihydrofolate. To consider binding without added osmolytes, a high-pressure approach was used. In this study, the type II enzyme, R67 DHFR, was subjected to high hydrostatic pressure (HHP). Both enzyme activity and fluorescence measurements find the protein tolerates pressures up to 200 MPa. Binding of the cofactor to R67 DHFR weakens with increasing pressure, and a positive association volume of 11.4 ± 0.5 cm(3)/mol was measured. Additionally, an activation volume of 3.3 ± 0.5 cm(3)/mol describing k(cat)/K(m(DHF)) was determined from progress curve analysis. Results from these HHP experiments suggest water release accompanies binding of both the cofactor and DHF to R67 DHFR. In an additional set of experiments, isothermal titration calorimetry studies in H2O and D2O find that water reorganization dominates the enthalpy associated with binding of DHF to R67 DHFR·NADP(+), while no obvious effects occur for cofactor binding. The combined results indicate that water plays an active role in ligand binding to R67 DHFR.

Fragment-based discovery of novel thymidylate synthase leads by NMR screening and group epitope mapping

Solution-state nuclear magnetic resonance (NMR) is a versatile tool for the study of binding interactions between small molecules and macromolecular targets. We applied ligand-based NMR techniques to the study of human thymidylate synthase (hTS) using known nanomolar inhibitors and a library of small molecule fragments. Screening by NMR led to the rapid identification of ligand pairs that bind in proximal sites within the cofactor-binding pocket of hTS. Screening hits were used as search criteria within commercially available sources, and a subset of catalog analogs were tested for potency by in vitro assay and binding affinity by quantitative saturation transfer difference (STD)-NMR titration. Two compounds identified by this approach possess low micromolar affinity and potency, as well as excellent binding efficiency against hTS. Relative binding orientations for both leads were modeled using AutoDock, and the most likely bound conformations were validated using experimentally derived STD-NMR binding epitope data. These ligands represent novel starting points for fragment-based drug design of non-canonical TS inhibitors, and their binding epitopes highlight important and previously unexploited interactions with conserved residues in the cofactor-binding site.

Crystal structure of a type II dihydrofolate reductase catalytic ternary complex

Type II dihydrofolate reductase (DHFR) is a plasmid-encoded enzyme that confers resistance to bacterial DHFR-targeted antifolate drugs. It forms a symmetric homotetramer with a central pore which functions as the active site. Its unusual structure, which results in a promiscuous binding surface that accommodates either the dihydrofolate (DHF) substrate or the NADPH cofactor, has constituted a significant limitation to efforts to understand its substrate specificity and reaction mechanism. We describe here the first structure of a ternary R67 DHFR.DHF.NADP+ catalytic complex, resolved to 1.26 A. This structure provides the first clear picture of how this enzyme, which lacks the active site carboxyl residue that is ubiquitous in Type I DHFRs, is able to function. In the catalytic complex, the polar backbone atoms of two symmetry-related I68 residues provide recognition motifs that interact with the carboxamide on the nicotinamide ring, and the N3-O4 amide function on the pteridine ring. This set of interactions orients the aromatic rings of substrate and cofactor in a relative endo geometry in which the reactive centers are held in close proximity. Additionally, a central, hydrogen-bonded network consisting of two pairs of Y69-Q67-Q67'-Y69' residues provides an unusually tight interface, which appears to serve as a "molecular clamp" holding the substrates in place in an orientation conducive to hydride transfer. In addition to providing the first clear insight regarding how this extremely unusual enzyme is able to function, the structure of the ternary complex provides general insights into how a mutationally challenged enzyme, i.e., an enzyme whose evolution is restricted to four-residues-at-a-time active site mutations, overcomes this fundamental limitation.

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

Structure of the Q67H mutant of R67 dihydrofolate reductase-NADP+ complex reveals a novel cofactor binding mode

Plasmid-encoded bacterial R67 dihydrofolate reductase (DHFR) is a NADPH-dependent enzyme unrelated to chromosomal DHFR in amino acid sequence and structure. R67 DHFR is insensitive to the bacterial drug trimethoprim in contrast to chromosomal DHFR. The crystal structure of Q67H mutant of R67 DHFR bound to NADP(+) has been determined at 1.15 angstroms resolution. The cofactor assumes an extended conformation with the nicotinamide ring bound near the center of the active site pore, the ribose and pyrophosphate group (PP(i)) extending toward the outer pore. The ribonicotinamide exhibits anti conformation as in chromosomal DHFR complexes. The relative orientation between the PP(i) and the nicotinamide ribose differs from that observed in chromosomal DHFR-NADP(+) complexes. The coenzyme displays symmetrical binding mode with several water-mediated hydrogen bonds with the protein besides ionic, stacking, and van der Waals interactions. The structure provides a molecular basis for the observed stoichiometry and cooperativity in ligand binding. The ternary model based on the present structure and the previous R67 DHFR-folate complex provides insight into the catalytic mechanism and indicates that the relative orientation of the reactants in plasmid DHFR is different from that seen in chromosomal DHFRs.

Enhanced signal dispersion in saturation transfer difference experiments by conversion to a 1D-STD-homodecoupled spectrum

The saturation transfer difference (STD) experiment is a rich source of information on topological aspects of ligand binding to a receptor. The epitope mapping is based on a magnetization transfer after signal saturation from the receptor to the ligand, where interproton distances permit this process. Signal overlap in the STD spectrum can cause difficulties to correctly assign and/or quantitate the measured enhancements. To address this issue we report here a modified version of the routine experiment and a processing scheme that provides a 1D-STD homodecoupled spectrum (i.e. an experiment in which all STD signals appear as singlets) with line widths similar to those in original STD spectrum. These refinements contribute to alleviate problems of signal overlap. The experiment is based on 2D-J-resolved spectroscopy, one of the fastest 2D experiments under conventional data sampling in the indirect dimension, and provides excellent sensitivity, a key factor for the difference experiments.

Integron-sequestered dihydrofolate reductase: a recently redeployed enzyme

The introduction and wide use of antibacterial drugs has resulted in the emergence of resistant organisms. DfrB dihydrofolate reductase (DHFR) is a bacterial enzyme that is uniquely associated with mobile gene cassettes within integrons, and confers resistance to the drug trimethoprim. This enzyme has intrigued microbiologists since it was discovered more than thirty years ago because of its simple structure, enzymatic inefficiency and its virtual insensitivity to trimethoprim. Here, for the first time, a comprehensive discussion of genetic, evolutionary, structural and functional studies of this enzyme is presented together. This information supports the ideas that DfrB DHFR is a poorly adapted catalyst and has recently been recruited to perform a novel enzymatic activity in response to selective pressure.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Searching sequence space: two different approaches to dihydrofolate reductase catalysis

There are numerous examples of proteins that catalyze the same reaction while possessing different structures. This review focuses on two dihydrofolate reductases (DHFRs) that have disparate structures and discusses how the catalytic strategies of these two DHFRs are driven by their respective scaffolds. The two proteins are E. coli chromosomal DHFR (Ec DHFR) and a type II R-plasmid-encoded DHFR, typified by R67 DHFR. The former has been described as a very well evolved enzyme with an efficiency of 0.15, while the latter has been suggested to be a model for a "primitive" enzyme that has not yet been optimized by evolution. This comparison underlines what is important to catalysis in these two enzymes and concurrently highlights fundamental issues in enzyme catalysis.

Targeting apoptosis via chemical design: inhibition of bid-induced cell death by small organic molecules

Bid is a key member of the Bcl-2 family proteins involved in the control of the apoptotic cascade in cells, leading to cell death. Uncontrolled cell death is associated with several human pathologies, such as neurodegenerative diseases and ischemic injuries. Therefore, Bid represents a potential yet unexplored and challenging target for strategies aimed at the development of therapeutic agents. Here we show that a multidisciplinary NMR-based approach that we named SAR by ILOEs (structure activity relationships by interligand nuclear Overhauser effect) allowed us to rationally design a series of 4-phenylsulfanyl-phenylamine derivatives that are capable of occupying a deep hydrophobic crevice on the surface of Bid. These compounds represent the first antiapoptotic small molecules targeting a Bcl-2 protein as shown by their ability to inhibit tBid-induced SMAC release, caspase-3 activation, and cell death.

"Catch 222," the effects of symmetry on ligand binding and catalysis in R67 dihydrofolate reductase as determined by mutations at Tyr-69

R67 dihydrofolate reductase (R67 DHFR) catalyzes the transfer of a hydride ion from NADPH to dihydrofolate, generating tetrahydrofolate. The homotetrameric enzyme provides a unique environment for catalysis as both ligands bind within a single active site pore possessing 222 symmetry. Mutation of one active site residue results in concurrent mutation of three additional symmetry-related residues, causing large effects on binding of both ligands as well as catalysis. For example, mutation of symmetry-related tyrosine 69 residues to phenylalanine (Y69F), results in large increases in Km values for both ligands and a 2-fold rise in the kcat value for the reaction (Strader, M. B., Smiley, R. D., Stinnett, L. G., VerBerkmoes, N. C., and Howell, E. E. (2001) Biochemistry 40, 11344-11352). To understand the interactions between specific Tyr-69 residues and each ligand, asymmetric Y69F mutants were generated that contain one to four Y69F mutations. A general trend observed from isothermal titration calorimetry and steady-state kinetic studies of these asymmetric mutants is that increasing the number of Y69F mutations results in an increase in the Kd and Km values. In addition, a comparison of steady-state kinetic values suggests that two Tyr-69 residues in one half of the active site pore are necessary for NADPH to exhibit a wild-type Km value. A tyrosine 69 to leucine mutant was also generated to approach the type(s) of interaction(s) occurring between Tyr-69 residues and the ligands. These studies suggest that the hydroxyl group of Tyr-69 is important for interactions with NADPH, whereas both the hydroxyl group and hydrophobic ring atoms of the Tyr-69 residues are necessary for proper interactions with dihydrofolate.

Role of Lys-32 residues in R67 dihydrofolate reductase probed by asymmetric mutations

R67 dihydrofolate reductase (R67 DHFR) is a novel protein encoded by an R-plasmid that confers resistance to the antibiotic, trimethoprim. This homotetrameric enzyme possesses 222 symmetry, which imposes numerous constraints on the single active site pore, including a "one-site-fits-both" strategy for binding its ligands, dihydrofolate (DHF) and NADPH. Previous studies uncovered salt effects on binding and catalysis (Hicks, S. N., Smiley, R. D., Hamilton, J. B., and Howell, E. E. (2003) Biochemistry 42, 10569-10578), however the one or more residues that participate in ionic contacts with the negatively charged tail of DHF as well as the phosphate groups in NADPH were not identified. Several studies predict that Lys-32 residues were involved, however mutations at this residue destabilize the R67 DHFR homotetramer. To study the role of Lys-32 in binding and catalysis, asymmetric K32M mutations have been utilized. To create asymmetry, individual mutations were added to a tandem array of four in-frame gene copies. These studies show one K32M mutation is tolerated quite well, whereas addition of two mutations has variable effects. Two double mutants, K32M:1+2 and K32M: 1+4, which place the mutations on opposite sides of the pore, reduce kcat. However a third double mutant, K32M: 1+3, that places two mutations on the same half pore, enhances kcat 4- to 5-fold compared with the parent enzyme, albeit at the expense of weaker binding of ligands. Because the kcat/Km values for this double mutant series are similar, these mutations appear to have uncovered some degree of non-productive binding. This non-productive binding mode likely arises from formation of an ionic interaction that must be broken to allow access to the transition state. The K32M:1+3 mutant data suggest this interaction is an ionic interaction between Lys-32 and the charged tail of dihydrofolate. This unusual catalytic scenario arises from the 222 symmetry imposed on the single active site pore.

The effect of relaxation on the epitope mapping by saturation transfer difference NMR

The effect of longitudinal relaxation of ligand protons on saturation transfer difference (STD) was investigated by using a known binding system, dihydrofolate reductase and trimethoprim. The results indicate that T1 relaxation of ligand protons has a severe interference on the epitope map derived from a STD measurement. When the T1s of individual ligand protons are distinctly different, STD experiments may not give an accurate epitope map for the ligand-target interactions. Measuring the relaxation times prior to mapping is strongly advised. A saturation time shorter than T1s is suggested for improving the potential epitope map. Reduction in temperature was seen to enhance the saturation efficiency in small to medium size targets.

Vibrational structure of dihydrofolate bound to R67 dihydrofolate reductase

R67 is a Type II dihydrofolate reductase (DHFR) that catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate by facilitating the addition of a proton to N5 of DHF and the transfer of a hydride ion from NADPH to C6. Because this enzyme is a plasmid-encoded DHFR from trimethoprim-resistant bacteria, extensive studies on R67 with various methods have been performed to elucidate its reaction mechanism. Here, Raman difference measurements, conducted on the ternary complex of R67.NADP(+).DHF believed to be an accurate mimic of the productive DHFR.NADPH.DHF complex, show that the pK(a) of N5 in the complex is less than 4. This is in clear contrast to the behavior observed in Escherichia coli DHFR, a substantially more efficient enzyme, where the pK(a) of bound DHF at N5 is increased to 6.5 compared with its solution value of 2.6. A comparison of the ternary complexes in R67 and E. coli DHFRs suggests that enzymic raising of the pK(a) at N5 can significantly increase the catalytic efficiency of the hydride transfer step. However, R67 shows that even without such a strategy an effective DHFR can still be designed.

One site fits both: a model for the ternary complex of folate + NADPH in R67 dihydrofolate reductase, a D2 symmetric enzyme

R67 dihydrofolate reductase (DHFR) is a novel enzyme that confers resistance to the antibiotic trimethoprim. The crystal structure of R67 DHFR displays a toroidal structure with a central active-site pore. This homotetrameric protein exhibits 222 symmetry, with only a few residues from each chain contributing to the active site, so related sites must be used to bind both substrate (dihydrofolate) and cofactor (NADPH) in the productive R67 DHFR.NADPH.dihydrofolate complex. Whereas the site of folate binding has been partially resolved crystallographically, an interesting question remains: how can the highly symmetrical active site also bind and orient NADPH for catalysis? To model this ternary complex, we employed DOCK and SLIDE, two methods for docking flexible ligands into proteins using quite different algorithms. The bound pteridine ring of folate (Fol I) from the crystal structure of R67 DHFR was used as the basis for docking the nicotinamide-ribose-Pi (NMN) moiety of NADPH. NMN was positioned by both DOCK and SLIDE on the opposite side of the pore from Fol I, where it interacts with Fol I at the pore's center. Numerous residues serve dual roles in binding. For example, Gln 67 from both the B and D subunits has several contacts with the pteridine ring, while the same residue from the A and C subunits has several contacts with the nicotinamide ring. The residues involved in dual roles are generally amphipathic, allowing them to make both hydrophobic and hydrophilic contacts with the ligands. The result is a 'hot spot' binding surface allowing the same residues to co-optimize the binding of two ligands, and orient them for catalysis.