Tales of Dihydrofolate Binding to R67 Dihydrofolate Reductase

From a binding module to essential catalytic activity: how nature stumbled on a good thing

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Dual-Target Inhibitors of the Folate Pathway Inhibit Intrinsically Trimethoprim-Resistant DfrB Dihydrofolate Reductases

Trimethoprim (TMP) is widely used to treat infections in humans and in livestock, accelerating the incidence of TMP resistance. The emergent and largely untracked type II dihydrofolate reductases (DfrBs) are intrinsically TMP-resistant plasmid-borne Dfrs that are structurally and evolutionarily unrelated to chromosomal Dfrs. We report kinetic characterization of the known DfrB family members. Their kinetic constants are conserved and all are poorly inhibited by TMP, consistent with TMP resistance. We investigate their inhibition with known and novel bisubstrate inhibitors of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK). Importantly, all are inhibited by the HPPK inhibitors, making these molecules dual-target inhibitors of two folate pathway enzymes that are strictly microbial.

Structure-Based Design of Dimeric Bisbenzimidazole Inhibitors to an Emergent Trimethoprim-Resistant Type II Dihydrofolate Reductase Guides the Design of Monomeric Analogues

The worldwide use of the broad-spectrum antimicrobial trimethoprim (TMP) has induced the rise of TMP-resistant microorganisms. In addition to resistance-causing mutations of the microbial chromosomal dihydrofolate reductase (Dfr), the evolutionarily and structurally unrelated type II Dfrs (DfrBs) have been identified in TMP-resistant microorganisms. DfrBs are intrinsically TMP-resistant and allow bacterial proliferation when the microbial chromosomal Dfr is TMP-inhibited, making these enzymes important targets for inhibitor development. Furthermore, DfrBs occur in multiresistance plasmids, potentially accelerating their dissemination. We previously reported symmetrical bisbenzimidazoles that are the first selective inhibitors of the only well-characterized DfrB, DfrB1. Here, their diversification provides a new series of inhibitors (K _i = 1.7-12.0 μM). Our results reveal two prominent features: terminal carboxylates and inhibitor length allow the establishment of essential interactions with DfrB1. Two crystal structures demonstrate the simultaneous binding of two inhibitor molecules in the symmetrical active site. Observations of those dimeric inhibitors inspired the design of monomeric analogues, binding in a single copy yet offering similar inhibition potency (K _i = 1.1 and 7.4 μM). Inhibition of a second member of the DfrB family, DfrB4, suggests the generality of these inhibitors. These results provide key insights into inhibition of the highly TMP-resistant DfrBs, opening avenues to downstream development of antibiotics for combatting this emergent source of resistance.

Small Angle Neutron Scattering Studies of R67 Dihydrofolate Reductase, a Tetrameric Protein with Intrinsically Disordered N-Termini

R67 dihydrofolate reductase (DHFR) is a homotetramer with a single active site pore and no sequence or structural homology with chromosomal DHFRs. The R67 enzyme provides resistance to trimethoprim, an active site-directed inhibitor of Escherichia coli DHFR. Sixteen to twenty N-terminal amino acids are intrinsically disordered in the R67 dimer crystal structure. Chymotrypsin cleavage of 16 N-terminal residues results in an active enzyme with a decreased stability. The space sampled by the disordered N-termini of R67 DHFR was investigated using small angle neutron scattering. From a combined analysis using molecular dynamics and the program SASSIE ( http://www.smallangles.net/sassie/SASSIE_HOME.html ), the apoenzyme displays a radius of gyration (Rg) of 21.46 ± 0.50 Å. Addition of glycine betaine, an osmolyte, does not result in folding of the termini as the Rg increases slightly to 22.78 ± 0.87 Å. SASSIE fits of the latter SANS data indicate that the disordered N-termini sample larger regions of space and remain disordered, suggesting they might function as entropic bristles. Pressure perturbation calorimetry also indicated that the volume of R67 DHFR increases upon addition of 10% betaine and decreased at 20% betaine because of the dehydration of the protein. Studies of the hydration of full-length R67 DHFR in the presence of the osmolytes betaine and dimethyl sulfoxide find around 1250 water molecules hydrating the protein. Similar studies with truncated R67 DHFR yield around 400 water molecules hydrating the protein in the presence of betaine. The difference of ∼900 waters indicates the N-termini are well-hydrated.

Integron-Associated DfrB4, a Previously Uncharacterized Member of the Trimethoprim-Resistant Dihydrofolate Reductase B Family, Is a Clinically Identified Emergent Source of Antibiotic Resistance

Whole-genome sequencing of trimethoprim-resistant Escherichia coli clinical isolates identified a member of the trimethoprim-resistant type II dihydrofolate reductase gene family (dfrB). The dfrB4 gene was located within a class I integron flanked by multiple resistance genes. This arrangement was previously reported in a 130.6-kb multiresistance plasmid. The DfrB4 protein conferred a >2,000-fold increased trimethoprim resistance on overexpression in E. coli Our results are consistent with the finding that dfrB4 contributes to clinical trimethoprim resistance.

Aspects of Weak Interactions between Folate and Glycine Betaine

Folate, or vitamin B9, is an important compound in one-carbon metabolism. Previous studies have found weaker binding of dihydrofolate to dihydrofolate reductase in the presence of osmolytes. In other words, osmolytes are more difficult to remove from the dihydrofolate solvation shell than water; this shifts the equilibrium toward the free ligand and protein species. This study uses vapor-pressure osmometry to explore the interaction of folate with the model osmolyte, glycine betaine. This method yields a preferential interaction potential (μ23/RT value). This value is concentration-dependent as folate dimerizes. The μ23/RT value also tracks the deprotonation of folate's N3-O4 keto-enol group, yielding a pKa of 8.1. To determine which folate atoms interact most strongly with betaine, the interaction of heterocyclic aromatic compounds (as well as other small molecules) with betaine was monitored. Using an accessible surface area approach coupled with osmometry measurements, deconvolution of the μ23/RT values into α values for atom types was achieved. This allows prediction of μ23/RT values for larger molecules such as folate. Molecular dynamics simulations of folate show a variety of structures from extended to L-shaped. These conformers possess μ23/RT values from -0.18 to 0.09 m-1, where a negative value indicates a preference for solvation by betaine and a positive value indicates a preference for water. This range of values is consistent with values observed in osmometry and solubility experiments. As the average predicted folate μ23/RT value is near zero, this indicates folate interacts almost equally well with betaine and water. Specifically, the glutamate tail prefers to interact with water, while the aromatic rings prefer betaine. In general, the more protonated species in our small molecule survey interact better with betaine as they provide a source of hydrogens (betaine is not a hydrogen bond donor). Upon deprotonation of the small molecule, the preference swings toward water interaction because of its hydrogen bond donating capacities.