Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Structural insights into the effects of glycerol on ligand binding to cytochrome P450

New antitubercular drugs are vital due to the spread of resistant strains. Carbethoxyhexyl imidazole (CHImi) inhibits cytochrome P450 CYP124, which is a steroid-metabolizing enzyme that is important for the survival of Mycobacterium tuberculosis in macrophages. The available crystal structure of the CYP124-CHImi complex reveals two glycerol molecules in the active site. A 1.15 Å resolution crystal structure of the glycerol-free CYP124-CHimi complex reported here shows multiple conformations of CHImi and the CYP124 active site which were previously restricted by glycerol. Complementary molecular dynamics simulations show coherence of the ligand and enzyme conformations. Spectrophotometric titration confirmed the influence of glycerol on CHImi binding: the affinity decreases more than tenfold in glycerol-containing buffer. In addition, it also showed that glycerol has a similar effect on other azole and triazole CYP124 ligands. Together, these data show that glycerol may compromise structural-functional studies and impede rational drug-design campaigns.

Effects of Osmolytes on Ligand Binding to Dihydropteroate Synthase from Bacillus anthracis

Osmolyte interactions with ligands can affect their affinity for proteins and are dependent upon the cosolute and the functional groups of the ligand. Here, we explored ligand binding to Bacillus anthracis dihydropteroate synthase (BaDHPS) under osmotic stress conditions. Osmolyte effects were specific to the cosolute and ligand, suggesting interaction of the osmolytes with the free ligands in solution. The association rates of pterin pyrophosphate were mostly unaffected by the osmolytes, except for a 2-fold decrease in the presence of 1 M trehalose, while the dissociation rates decreased in most osmolyte solutions. The viscosity and dielectric constant of the solution did not correlate with the effects of the osmolytes. Experimental results were compared with predicted preferential interaction coefficients (Δμ₂₃/RT) between the osmolytes and ligands. The Δμ₂₃/RT were able to predict the experimental data for most of the osmolytes. Trehalose and proline effects did not correlate with the predicted values, indicating that these two osmolytes may affect binding in more complex ways than simple preferential interactions. Additionally, osmolytes weakly interacted with the sulfa drug sulfathiazole, which altered its affinity for BaDHPS, suggesting that these types of weak interactions can also impact drug binding. As osmolytes affect ligands binding to two different folate cycle enzymes (DHFRs and DHPS), we predicted how ligand binding to other folate cycle enzymes will be altered by the presence of osmolytes.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP⁺ binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These nonspecific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

In Vivo Titration of Folate Pathway Enzymes

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

Modulating Enzyme Activity by Altering Protein Dynamics with Solvent

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized; however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ∼3 fold at 20% (v/v) IPA concentration. Stopped-flow kinetic measurements find that the pH-independent k_hydride rate decreases by 2.2 fold. X-ray crystallographic enzyme structures show no noticeable differences, while computational studies indicate that the transition state and electrostatic effects were identical for water and mixed solvent conditions; quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzyme's ability to access its functionally relevant conformational substates, explaining the decreased k_hydride rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

Investigation of Classical Organic and Ionic Liquid Cosolvents for Early-Stage Screening in Fragment-Based Inhibitor Design with Unrelated Bacterial and Human Dihydrofolate Reductases

Drug design by methods such as fragment screening requires effective solubilization of millimolar concentrations of small organic compounds while maintaining the properties of the biological target. We investigate four organic solvents and three 1-butyl-3-methylimidazolium (BMIm)-based ionic liquids (ILs) as cosolvents to establish conditions for screening two structurally unrelated dihydrofolate reductases (DHFRs) that are prime drug targets. Moderate concentrations (10%-15%) of cosolvents had little effect on inhibition of the microbial type II R67 DHFR and of human DHFR (hDHFR), while higher concentrations of organic cosolvents generally decreased activity of both DHFRs. In contrast, a specific IL conserved the activity of one DHFR, while severely reducing the activity of the other, and vice versa, illustrating the differing effect of ILs on distinct protein folds. Most of the cosolvents investigated preserved the fold of R67 DHFR and had little effect on binding of the cofactor NADPH, but reduced the productive affinity for its substrate. In contrast, cosolvents resulted in modest structural destabilization of hDHFR with little effect on productive affinity. We conclude that the organic cosolvents, methanol, dimethylformamide, and dimethylsulfoxide, offer the most balanced conditions for early-stage compound screening as they maintain sufficient biological activity of both DHFRs while allowing for compound dissolution in the millimolar range. However, IL cosolvents showed poor capacity to solubilize organic compounds at millimolar concentrations, mitigating their utility in early-stage screening. Nonetheless, ILs could provide an alternative to classical organic cosolvents when low concentrations of inhibitors are used, as when characterizing higher affinity inhibitors.

Halophilic mechanism of the enzymatic function of a moderately halophilic dihydrofolate reductase from Haloarcula japonica strain TR-1

Dihydrofolate (DHF) reductase coded by a plasmid of the extremely halophilic archaeon Haloarcula japonica strain TR-1 (HjDHFR P1) shows moderate halophilicity on enzymatic activity at pH 6.0, although there is no significant effect of NaCl on its secondary structure. To elucidate the salt-activation and -inactivation mechanisms of this enzyme, we investigated the effects of pH and salt concentration, deuterium isotope effect, steady-state kinetics, and rapid-phase ligand-binding kinetics. Enzyme activity was increased eightfold by the addition of 500 mM NaCl at pH 6.0, fourfold by 250 mM at pH 8.0, and became independent of salt concentration at pH 10.0. Full isotope effects observed at pH 10.0 under 0-1000 mM NaCl indicated that the rate of hydride transfer, which was the rate-determining step at the basic pH region, was independent of salt concentration. Conversely, rapid-phase ligand-binding experiments showed that the amplitude of the DHF-binding reaction increased and the tetrahydrofolate (THF)-releasing rate decreased with increasing NaCl concentration. These results suggested that the salt-activation mechanism of HjDHFR P1 is via the population change of the anion-unbound and anion-bound conformers, which are binding-incompetent and -competent conformations for DHF, respectively, while that of salt inactivation is via deceleration of the THF-releasing rate, which is the rate-determining step at the neutral pH region.

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.

Ligand-Dependent Conformational Dynamics of Dihydrofolate Reductase

Enzymes are known to change among several conformational states during turnover. The role of such dynamic structural changes in catalysis is not fully understood. The influence of dynamics in catalysis can be inferred, but not proven, by comparison of equilibrium structures of protein variants and protein-ligand complexes. A more direct way to establish connections between protein dynamics and the catalytic cycle is to probe the kinetics of specific protein motions in comparison to progress along the reaction coordinate. We have examined the enzyme model system dihydrofolate reductase (DHFR) from Escherichia coli with tryptophan fluorescence-probed temperature-jump spectroscopy. We aimed to observe the kinetics of the ligand binding and ligand-induced conformational changes of three DHFR complexes to establish the relationship among these catalytic steps. Surprisingly, in all three complexes, the observed kinetics do not match a simple sequential two-step process. Through analysis of the relationship between ligand concentration and observed rate, we conclude that the observed kinetics correspond to the ligand binding step of the reaction and a noncoupled enzyme conformational change. The kinetics of the conformational change vary with the ligand's identity and presence but do not appear to be directly related to progress along the reaction coordinate. These results emphasize the need for kinetic studies of DHFR with highly specific spectroscopic probes to determine which dynamic events are coupled to the catalytic cycle and which are not.

Network of remote and local protein dynamics in dihydrofolate reductase catalysis

Molecular dynamics calculations and bionformatic studies of dihydrofolate reductase (DHFR) have suggested a network of coupled motions across the whole protein that is correlated to the reaction coordinate. Experimental studies demonstrated that distal residues G121, M42 and F125 in E. coli DHFR participate in that network. The missing link in our understanding of DHFR catalysis is the lack of a mechanism by which such remote residues can affect the catalyzed chemistry at the active site. Here, we present a study of the temperature dependence of intrinsic kinetic isotope effects (KIEs) that indicates synergism between a remote residue in that dynamic network, G121, and the active site's residue I14. The intrinsic KIEs for the I14A-G121V double mutant showed steeper temperature dependence (ΔE_a(T-H)) than expected from comparison of the wild type and two single mutants. That effect was non-additive, i.e., ΔE_a(T-H)G121V +ΔE_{a(T-H) I14A} < ΔE_{a(T-H) double mutant}, which indicates a synergism between the two residues. This finding links the remote residues in the network under investigation to the enzyme's active site, providing a mechanism by which these residues can be coupled to the catalyzed chemistry. This experimental evidence validates calculations proposing that both remote and active site residues constitute a network of coupled promoting motions correlated to the bond activation step (C-H→C hydride transfer in this case). Additionally, the effect of I14A and G121V mutations on single turnover rates was additive rather than synergistic. Although single turnover rate measurements are more readily available and thus more popular than assessing intrinsic kinetic isotope effects, the current finding demonstrates that for these rates, which in DHFR reflect several microscopic rate constants, can fall short of revealing the nature of the C-H bond activation per se.

SEDPHAT--a platform for global ITC analysis and global multi-method analysis of molecular interactions

Isothermal titration calorimetry experiments can provide significantly more detailed information about molecular interactions when combined in global analysis. For example, global analysis can improve the precision of binding affinity and enthalpy, and of possible linkage parameters, even for simple bimolecular interactions, and greatly facilitate the study of multi-site and multi-component systems with competition or cooperativity. A pre-requisite for global analysis is the departure from the traditional binding model, including an 'n'-value describing unphysical, non-integral numbers of sites. Instead, concentration correction factors can be introduced to account for either errors in the concentration determination or for the presence of inactive fractions of material. SEDPHAT is a computer program that embeds these ideas and provides a graphical user interface for the seamless combination of biophysical experiments to be globally modeled with a large number of different binding models. It offers statistical tools for the rigorous determination of parameter errors, correlations, as well as advanced statistical functions for global ITC (gITC) and global multi-method analysis (GMMA). SEDPHAT will also take full advantage of error bars of individual titration data points determined with the unbiased integration software NITPIC. The present communication reviews principles and strategies of global analysis for ITC and its extension to GMMA in SEDPHAT. We will also introduce a new graphical tool for aiding experimental design by surveying the concentration space and generating simulated data sets, which can be subsequently statistically examined for their information content. This procedure can replace the 'c'-value as an experimental design parameter, which ceases to be helpful for multi-site systems and in the context of gITC.

Thermodynamics and solvent linkage of macromolecule-ligand interactions

Binding involves two steps, desolvation and association. While water is ubiquitous and occurs at high concentration, it is typically ignored. In vitro experiments typically use infinite dilution conditions, while in vivo, the concentration of water is decreased due to the presence of high concentrations of molecules in the cellular milieu. This review discusses isothermal titration calorimetry approaches that address the role of water in binding. For example, use of D2O allows the contribution of solvent reorganization to the enthalpy component to be assessed. Further, the addition of osmolytes will decrease the water activity of a solution and allow effects on Ka to be determined. In most cases, binding becomes tighter in the presence of osmolytes as the desolvation penalty associated with binding is minimized. In other cases, the osmolytes prefer to interact with the ligand or protein, and if their removal is more difficult than shedding water, then binding can be weakened. These complicating layers can be discerned by different slopes in ln(Ka) vs osmolality plots and by differential scanning calorimetry in the presence of the osmolyte.

Investigation of osmolyte effects on FolM: comparison with other dihydrofolate reductases

A weak association between osmolytes and dihydrofolate (DHF) decreases the affinity of the substrate for the Escherichia coli chromosomal and R67 plasmid dihydrofolate reductase (DHFR) enzymes. To test whether the osmolyte-DHF association also interferes with binding of DHF to FolM, an E. coli enzyme that possesses weak DHFR activity, ligand binding was monitored in the presence of osmolytes. The affinity of FolM for DHF, measured by kcat/Km(DHF), was decreased by the addition of an osmolyte. Additionally, binding of the antifolate drug, methotrexate, to FolM was weakened by the addition of an osmolyte. The changes in ligand binding with water activity were unique for each osmolyte, indicating preferential interaction between the osmolyte and folate and its derivatives; however, additional evidence provided support for further interactions between FolM and osmolytes. Binding of the reduced nicotinamide adenine dinucleotide phosphate (NADPH) cofactor to FolM was monitored by isothermal titration calorimetry as a control for protein-osmolyte association. In the presence of betaine (proposed to be the osmolyte most excluded from protein surfaces), the NADPH Kd decreased, consistent with dehydration effects. However, other osmolytes did not tighten binding to the cofactor. Rather, dimethyl sulfoxide (DMSO) had no effect on the NADPH Kd, while ethylene glycol and polyethylene glycol 400 weakened cofactor binding. Differential scanning calorimetry of FolM in the presence of osmolytes showed that both DMSO and ethylene glycol decreased the stability of FolM, while betaine increased the stability of the protein. These results suggest that some osmolytes can destabilize FolM by preferentially interacting with the protein. Further, these weak attractions can impede ligand binding. These various contributions have to be considered when interpreting osmotic pressure results.

Quantifying the molecular origins of opposite solvent effects on protein-protein interactions

Although the nature of solvent-protein interactions is generally weak and non-specific, addition of cosolvents such as denaturants and osmolytes strengthens protein-protein interactions for some proteins, whereas it weakens protein-protein interactions for others. This is exemplified by the puzzling observation that addition of glycerol oppositely affects the association constants of two antibodies, D1.3 and D44.1, with lysozyme. To resolve this conundrum, we develop a methodology based on the thermodynamic principles of preferential interaction theory and the quantitative characterization of local protein solvation from molecular dynamics simulations. We find that changes of preferential solvent interactions at the protein-protein interface quantitatively account for the opposite effects of glycerol on the antibody-antigen association constants. Detailed characterization of local protein solvation in the free and associated protein states reveals how opposite solvent effects on protein-protein interactions depend on the extent of dewetting of the protein-protein contact region and on structural changes that alter cooperative solvent-protein interactions at the periphery of the protein-protein interface. These results demonstrate the direct relationship between macroscopic solvent effects on protein-protein interactions and atom-scale solvent-protein interactions, and establish a general methodology for predicting and understanding solvent effects on protein-protein interactions in diverse biological environments.

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.

Weak interactions between folate and osmolytes in solution

Previous osmotic stress studies on the role of solvent in two structurally unrelated dihydrofolate reductases (DHFRs) found weaker binding of dihydrofolate (DHF) to either enzyme in the presence of osmolytes. To explain these unusual results, weak interactions between DHF and osmolytes were proposed, with a competition between osmolyte and DHFR for DHF. High osmolyte concentrations will inhibit binding of the cognate pair. To evaluate this hypothesis, we devised a small molecule approach. Dimerization of folate, monitored by nuclear magnetic resonance, was weakened 2-3-fold upon addition of betaine or dimethyl sulfoxide (DMSO), supporting preferential interaction of either osmolyte with the monomer (as it possesses a larger surface area). Nuclear Overhauser effect (NOE) spectroscopy experiments found a positive NOE for the interaction of the C3'/C5' benzoyl ring protons with the C9 proton in buffer; however, a negative NOE was observed upon addition of betaine or DMSO. This change indicated a decreased tumbling rate, consistent with osmolyte interaction. Osmotic stress experiments also showed that betaine, DMSO, and sucrose preferentially interact with folate. Further, studies with the folate fragments, p-aminobenzoic acid and pterin 6-carboxylate, revealed interactions for both model compounds with betaine and sucrose. In contrast, DMSO was strongly excluded from the pterin ring but preferentially interacted with the p-aminobenzoyl moiety. These interactions are likely to be important in vivo because of the crowded conditions of the cell where weak contacts can more readily compete with specific binding interactions.