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

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.

Unexplored Excipients in Biotherapeutic Formulations: Natural Osmolytes as Potential Stabilizers Against Thermally Induced Aggregation of IgG1 Biotherapeutics

Monoclonal antibodies (mAbs), while incredibly successful, are prone to a variety of degradation pathways, the most significant of which is aggregation. One of the most commonly used strategy to overcome protein aggregation is addition of excipients to the formulation. Osmolytes such as trehalose, sucrose, and glycine are widely used. In this paper, we explore potential use of naturally occurring osmolytes such as betaine, sarcosine, ectoine, and hydroxyectoine for reducing aggregation of mAb therapeutics. Experimentation has been performed on two IgG1 mAbs via accelerated stability studies. A variety of analytical tools have been used for monitoring the impact, dynamic light scattering (DLS) for colloidal stability, Fourier transform infrared (FTIR) spectroscopy and fluorescence spectroscopy for conformational stability and the higher order structure (HOS), and differential scanning calorimetry (DSC) for thermal stability. No significant impact of osmolyte addition was observed on protein structure, on comparative Fc receptor (FcRn) binding, and on biocompatibility as per our hemolytic assay. Our results rank the osmolytes’ stabilizing trend to be sarcosine > betaine > hydroxyectoine > ectoine. Sarcosine emerged as the most successful osmolyte rendering highest degree of protection against aggregation. Our data support the prospect of using these osmolytes as successful excipients for mAb formulations.

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.

Changing relations between proteins and osmolytes: a choice of nature

The stabilization and destabilization of the protein in the presence of any additive is mainly attributed to its preferential exclusion from protein surface and its preferential binding to the protein surface, respectively.

Large cosolutes, small cosolutes, and dihydrofolate reductase activity

Protein enzymes are the main catalysts in the crowded and complex cellular interior, but their activity is almost always studied in dilute buffered solutions. Studies that attempt to recreate the cellular interior in vitro often utilize synthetic polymers as crowding agents. Here, we report the effects of the synthetic polymer cosolutes Ficoll, dextran, and polyvinylpyrrolidone, and their respective monomers, sucrose, glucose, and 1-ethyl-2-pyrrolidone, on the activity of the 18-kDa monomeric enzyme, Escherichia coli dihydrofolate reductase. At low concentrations, reductase activity increases relative to buffer and monomers, suggesting a macromolecular effect. However, the effect decreases at higher concentrations, approaching, and, in some cases, falling below buffer values. We also assessed activity in terms of volume occupancy, viscosity, and the overlap concentration (where polymers form an interwoven mesh). The trends vary with polymer family, but changes in activity are within threefold of buffer values. We also compiled and analyzed results from previous studies and conclude that alterations of steady-state enzyme kinetics in solutions crowded with synthetic polymers are idiosyncratic with respect to the crowding agent and enzyme.

Observed crowding effects on Mycobacterium tuberculosis 2-trans-enoyl-ACP (CoA) reductase enzyme activity are not due to excluded volume only

The cellular milieu is a complex and crowded aqueous solution. Macromolecular crowding effects are commonly studied in vitro using crowding agents. The aim of the present study was to evaluate the effects, if any, of macromolecular synthetic crowding agents on the apparent steady-state kinetic parameters (K m , k cat , and k cat /K m ) of Mycobacterium tuberculosis 2-trans-enoyl-ACP (CoA) reductase (InhA). Negligible effects on InhA activity were observed for ficoll 70, ficoll 400 and dextran 70. A complex effect was observed for PEG 6000. Glucose and sucrose showed, respectively, no effect on InhA activity and decreased k cat /K m for NADH and k cat for 2-trans-dodecenoyl-CoA. Molecular dynamics results suggest that InhA adopts a more compact conformer in sucrose solution. The effects of the crowding agents on the energy (E a and E η ), enthalpy (∆H # ), entropy (∆S # ), and Gibbs free energy (∆G # ) of activation were determined. The ∆G # values for all crowding agents were similar to buffer, suggesting that excluded volume effects did not facilitate stable activated ES # complex formation. Nonlinear Arrhenius plot for PEG 6000 suggests that "soft" interactions play a role in crowding effects. The results on InhA do not unequivocally meet the criteria for crowding effect due to exclude volume only.

Pterin function in bacteria

Pterins are widely conserved biomolecules that play essential roles in diverse organisms. First described as enzymatic cofactors in eukaryotic systems, bacterial pterins were discovered in cyanobacteria soon after. Several pterin structures unique to bacteria have been described, with conjugation to glycosides and nucleotides commonly observed. Despite this significant structural diversity, relatively few biological functions have been elucidated. Molybdopterin, the best studied bacterial pterin, plays an essential role in the function of the Moco cofactor. Moco is an essential component of molybdoenzymes such as sulfite oxidase, nitrate reductase, and dimethyl sulfoxide reductase, all of which play important roles in bacterial metabolism and global nutrient cycles. Outside of the molybdoenzymes, pterin cofactors play important roles in bacterial cyanide utilization and aromatic amino acid metabolism. Less is known about the roles of pterins in nonenzymatic processes. Cyanobacterial pterins have been implicated in phenotypes related to UV protection and phototaxis. Research describing the pterin-mediated control of cyclic nucleotide metabolism, and their influence on virulence and attachment, points to a possible role for pterins in regulation of bacterial behavior. In this review, we describe the variety of pterin functions in bacteria, compare and contrast structural and mechanistic differences, and illuminate promising avenues of future research.

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.

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.