Search

Dihydrofolate reductase (DHFR), due to its universality and the depth with which it has been studied, is a model system in the study of protein dynamics. Myriad previous works have identified networks of residues in positions near to and remote from the active site that are involved in the dynamics. For example, specific mutations on the Met20 loop in Escherichia coli DHFR (N23PP/S148A) are known to disrupt millisecond-time scale motions as well as reduce catalytic activity. However, how and if networks of dynamically coupled residues influence the evolution of DHFR is still an unanswered question. In this study, we first identify, by statistical coupling analysis and molecular dynamic simulations, a network of coevolving residues that possesses increased correlated motions. We then go on to show that allosteric communication in this network is knocked down in N23PP/S148A mutant E. coli DHFR. We also identify two sites in the human DHFR sector which may accommodate the Met20 loop double proline motif. Finally, we demonstrate a concerted evolutionary change in the human DHFR allosteric networks, which maintains dynamic communication. These findings strongly implicate protein dynamics as a driving force for evolution.

A new superior bacteria complementation model was achieved for testing antifolate compounds and investigating antifolate resistance in the dihydrofolate reductase (DHFR) enzyme of the malaria parasite. Earlier models depended on the addition of trimethoprim (TMP) to chemically suppress the host Escherichia coli (Ec) DHFR function. However, incomplete suppression of EcDHFR and potential interference of antibiotics needed to maintain plasmids for complementary gene expression can complicate the interpretations. To overcome such limitations, the folA (F) and thyA (T) genes were genetically knocked out (Δ) in E. coli BL21(DE3). The resulting EcΔFΔT cells were thymidine auxotroph where thymidine supplementation or functional complementation with heterologous DHFR-thymidylate synthase (TS) is needed to restore the loss of gene functions. When tested against pyrimethamine (PYR) and its analogs designed to target Plasmodium falciparum (Pf) DHFR-TS, the 50 % inhibitory concentration values obtained from EcΔFΔT surrogates expressing wildtype (PfTM4) or double mutant (PfK1) DHFR-TS showed strong correlations to the results obtained from the standard in vitro P. falciparum growth inhibition assay. Interestingly, while TMP had little effect on the susceptibility to PYR and analogs in EcΔFΔT expressing PfDHFR-TS, it hypersensitized the chemically knockdown E. coli BL21(DE3) expressing PfTM4 DHFR-TS but desensitized the one carrying PfK1 DHFR-TS. The low intrinsic expression level of PfTM4 in E. coli BL21(DE3) by western blot analysis may explain the hypersensitive to antifolates of chemical knockdown bacteria surrogate. These results demonstrated the usefulness of EcΔFΔT surrogate as a new tool for antifolate antimalarial screening with potential application for investigation of antifolate resistance mechanism., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier B.V.)

The Escherichia coli GroEL/ES chaperonin system facilitates protein folding in an ATP-driven manner. There are <100 obligate clients of this system in E. coli although GroEL can interact and assist the folding of a multitude of proteins in vitro. It has remained unclear, however, which features distinguish obligate clients from all the other proteins in an E. coli cell. To address this question, we established a system for selecting mutations in mouse dihydrofolate reductase (mDHFR), a GroEL interactor, that diminish its dependence on GroEL for folding. Strikingly, both synonymous and non-synonymous codon substitutions were found to reduce mDHFR's dependence on GroEL. The non-synonymous substitutions increase the rate of spontaneous folding whereas computational analysis indicates that the synonymous substitutions appear to affect translation rates at specific sites., (© 2024 The Author(s). Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

Inorganic materials depleted of heavy stable isotopes are known to deviate strongly in some physicochemical properties from their isotopically natural counterparts. Here we explored for the first time the effect of simultaneous depletion of the heavy carbon, hydrogen, oxygen and nitrogen isotopes on the bacterium E. coli and the enzymes expressed in it. Bacteria showed faster growth, with most proteins exhibiting higher thermal stability, while for recombinant enzymes expressed in depleted media, faster kinetics was discovered. At room temperature, luciferase, thioredoxin and dihydrofolate reductase and Pfu DNA polymerase showed up to a 250 % increase in activity compared to the native counterparts, with an additional ∼50 % increase at 10 °C. Diminished conformational and vibrational entropy is hypothesized to be the cause of the accelerated kinetics. Ultralight enzymes may find an application where extreme reaction rates are required., (© 2023 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Protein space is a rich analogy for genotype-phenotype maps, where amino acid sequence is organized into a high-dimensional space that highlights the connectivity between protein variants. It is a useful abstraction for understanding the process of evolution, and for efforts to engineer proteins towards desirable phenotypes. Few mentions of protein space consider how protein phenotypes can be described in terms of their biophysical components, nor do they rigorously interrogate how forces like epistasis-describing the nonlinear interaction between mutations and their phenotypic consequences-manifest across these components. In this study, we deconstruct a low-dimensional protein space of a bacterial enzyme (dihydrofolate reductase; DHFR) into "subspaces" corresponding to a set of kinetic and thermodynamic traits [k&#95;{cat}, K&#95;{M}, K&#95;{i}, and T&#95;{m} (melting temperature)]. We then examine how combinations of three mutations (eight alleles in total) display pleiotropy, or unique effects on individual subspace traits. We examine protein spaces across three orthologous DHFR enzymes (Escherichia coli, Listeria grayi, and Chlamydia muridarum), adding a genotypic context dimension through which epistasis occurs across subspaces. In doing so, we reveal that protein space is a deceptively complex notion, and that future applications to bioengineering should consider how interactions between amino acid substitutions manifest across different phenotypic subspaces.

Nutritive symbiosis between bacteria and ticks is observed across a range of ecological contexts; however, little characterization on the molecular components responsible for this symbiosis has been done. Previous studies in our lab demonstrated that Rickettsia monacensis str. Humboldt (strain Humboldt) can synthesize folate de novo via the folate biosynthesis pathway involving folA, folC, folE, folKP, and ptpS genes. In this study, expression of the strain Humboldt folA gene within a folA mutant Escherichia coli construct was used to functionally characterize the strain Humboldt folA folate gene in vivo. The strain Humboldt folA folate gene was subcloned into a TransBac vector and transformed into a folA mutant E. coli construct. The mutant containing strain Humboldt folA subclone and a pFE604 clone of the knocked-out folA gene was cured of pFE604. Curing of the folA mutant E. coli construct was successful using acridine orange and 43.5 °C incubation temperature. The plasmid curing assay showed curing efficiency of the folA mutant at 100%. Functional complementation was assessed by growth phenotype on minimal media with and without IPTG between strain Humboldt folA and E. coli folA. Large and homogenous wild-type colony growth was observed for both strain Humboldt and E. coli folA on minimal media with 0.1 mM IPTG, wild-type growth for strain Humboldt folA and pin-point growth for E. coli folA on 0.01 mM IPTG, and pin-point growth without IPTG for both strain Humboldt and E. coli folA. This study provides evidence substantiating the in vivo functionality of strain Humboldt folA in producing functional gene products for folate biosynthesis., Competing Interests: Declaration of Competing Interest The authors declare that there is no conflict of interest in this study., (Copyright © 2023 The Author(s). Published by Elsevier GmbH.. All rights reserved.)

AAA+ proteolytic machines play essential roles in maintaining and rebalancing the cellular proteome in response to stress, developmental cues, and environmental changes. Of the five AAA+ proteases in Escherichia coli, FtsH is unique in its attachment to the inner membrane and its function in degrading both membrane and cytosolic proteins. E. coli dihydrofolate reductase (DHFR) is a stable and biophysically well-characterized protein, which a previous study found resisted FtsH degradation despite the presence of an ssrA degron. By contrast, we find that FtsH degrades DHFR fused to a long peptide linker and ssrA tag. Surprisingly, we also find that FtsH degrades DHFR with shorter linkers and ssrA tag, and without any linker or tag. Thus, FtsH must be able to recognize a sequence element or elements within DHFR. We find that FtsH degradation of DHFR is noncanonical in the sense that it does not rely upon recognition of an unstructured polypeptide at or near the N-terminus or C-terminus of the substrate. Results using peptide-array experiments, mutant DHFR proteins, and fusion proteins suggest that FtsH recognizes an internal sequence in a species of DHFR that is partially unfolded. Overall, our findings provide insight into substrate recognition by FtsH and indicate that its degradation capacity is broader than previously reported., (© 2022 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

To what degree are individual structural elements within proteins modular such that similar structures from unrelated proteins can be interchanged? We study subdomain modularity by creating 20 chimeras of an enzyme, Escherichia coli dihydrofolate reductase (DHFR), in which a catalytically important, 10-residue α-helical sequence is replaced by α-helical sequences from a diverse set of proteins. The chimeras stably fold but have a range of diminished thermal stabilities and catalytic activities. Evolutionary coupling analysis indicates that the residues of this α-helix are under selection pressure to maintain catalytic activity in DHFR. Reversion to phenylalanine at key position 31 was found to partially restore catalytic activity, which could be explained by evolutionary coupling values. We performed molecular dynamics simulations using replica exchange with solute tempering. Chimeras with low catalytic activity exhibit nonhelical conformations that block the binding site and disrupt the positioning of the catalytically essential residue D27. Simulation observables and in vitro measurements of thermal stability and substrate-binding affinity are strongly correlated. Several E. coli strains with chromosomally integrated chimeric DHFRs can grow, with growth rates that follow predictions from a kinetic flux model that depends on the intracellular abundance and catalytic activity of DHFR. Our findings show that although α-helices are not universally substitutable, the molecular and fitness effects of modular segments can be predicted by the biophysical compatibility of the replacement segment., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a "sample and select" method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls ( 1 H, 15 N, and 13 C'). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution., (Copyright © 2020 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Many enzymes are known to change conformations during their catalytic cycle, but the role of these protein motions is not well understood. Escherichia coli dihydrofolate reductase (DHFR) is a small, flexible enzyme that is often used as a model system for understanding enzyme dynamics. Recently, native tryptophan fluorescence was used as a probe to study micro- to millisecond dynamics of DHFR. Yet, because DHFR has five native tryptophans, the origin of the observed conformational changes could not be assigned to a specific region within the enzyme. Here, we use DHFR mutants, each with a single tryptophan as a probe for temperature jump fluorescence spectroscopy, to further inform our understanding of DHFR dynamics. The equilibrium tryptophan fluorescence of the mutants shows that each tryptophan is in a different environment and that wild-type DHFR fluorescence is not a simple summation of all the individual tryptophan fluorescence signatures due to tryptophan-tryptophan interactions. Additionally, each mutant exhibits a two-phase relaxation profile corresponding to ligand association/dissociation convolved with associated conformational changes and a slow conformational change that is independent of ligand association and dissociation, similar to the wild-type enzyme. However, the relaxation rate of the slow phase depends on the location of the tryptophan within the enzyme, supporting the conclusion that the individual tryptophan fluorescence dynamics do not originate from a single collective motion, but instead report on local motions throughout the enzyme.

Evolutionary trajectories and mutational landscapes of drug-resistant bacteria are influenced by cell-intrinsic and extrinsic factors. In this study, I demonstrated that loss of the Lon protease altered susceptibility of Escherichia coli to trimethoprim and that these effects were strongly contingent on the drug concentration and genetic background. Lon, an AAA + ATPase, is a bacterial master regulator protease involved in cytokinesis, suppression of transposition events, and clearance of misfolded proteins. I show that Lon deficiency enhances intrinsic drug tolerance at sub-MIC levels of trimethoprim. As a result, loss of Lon, though disadvantageous under drug-free conditions, has a selective advantage at low concentrations of trimethoprim. At high drug concentrations, however, Lon deficiency is detrimental for E. coli I show that the former is explained by suppression of drug efflux by Lon, while the latter can be attributed to SulA-dependent hyperfilamentation. On the other hand, deletion of lon in a trimethoprim-resistant mutant E. coli strain (harboring the Trp30Gly dihydrofolate reductase [DHFR] allele) directly potentiates resistance by enhancing the in vivo stability of mutant DHFR. Using extensive mutational analysis at 3 hot spots of resistance, I show that many resistance-conferring mutations render DHFR prone to proteolysis. This trade-off between gaining resistance and losing in vivo stability limits the number of mutations in DHFR that can confer trimethoprim resistance. Loss of Lon expands the mutational capacity for acquisition of trimethoprim resistance. This paper identifies the multipronged action of Lon in trimethoprim resistance in E. coli and provides mechanistic insight into how genetic backgrounds and drug concentrations may alter the potential for antimicrobial resistance evolution. IMPORTANCE Understanding the evolutionary dynamics of antimicrobial resistance is vital to curb its emergence and spread. Being fundamentally similar to natural selection, the fitness of resistant mutants is a key parameter to consider in the evolutionary dynamics of antimicrobial resistance (AMR). Various intrinsic and extrinsic factors modulate the fitness of resistant bacteria. This study demonstrated that Lon, a bacterial master regulator protease, influences drug tolerance and resistance. Lon is a key regulator of several fundamental processes in bacteria, including cytokinesis. I demonstrated that Lon deficiency produces highly contingent phenotypes in E. coli challenged with trimethoprim and can expand the mutational repertoire available to E. coli to evolve resistance. This multipronged influence of Lon on drug resistance provides an illustrative instance of how master regulators shape the response of bacteria to antibiotics., (Copyright © 2020 American Society for Microbiology.)

Cell-based therapeutics have considerable promise across diverse medical specialties; however, reliable human imaging of the distribution and trafficking of genetically engineered cells remains a challenge. We developed positron emission tomography (PET) probes based on the small-molecule antibiotic trimethoprim (TMP) that can be used to image the expression of the Escherichia coli dihydrofolate reductase enzyme (eDHFR) and tested the ability of [ 18 F]-TMP, a fluorine-18 probe, to image primary human chimeric antigen receptor (CAR) T cells expressing the PET reporter gene eDHFR, yellow fluorescent protein (YFP), and Renilla luciferase (rLuc). Engineered T cells showed an approximately 50-fold increased bioluminescent imaging signal and 10-fold increased [ 18 F]-TMP uptake compared to controls in vitro. eDHFR-expressing anti-GD2 CAR T cells were then injected into mice bearing control GD2 - and GD2 + tumors. PET/computed tomography (CT) images acquired on days 7 and 13 demonstrated early residency of CAR T cells in the spleen followed by on-target redistribution to the GD2 + tumors. This was corroborated by autoradiography and anti-human CD8 immunohistochemistry. We found a high sensitivity of detection for identifying tumor-infiltrating CD8 CAR T cells, ∼11,000 cells per mm 3 . These data suggest that the [ 18 F]-TMP/eDHFR PET pair offers important advantages that could better allow investigators to monitor immune cell trafficking to tumors in patients., (Copyright © 2019 The American Society of Gene and Cell Therapy. Published by Elsevier Inc. All rights reserved.)

The mechanisms of adaptation to inactivation of essential genes remain unknown. Here we inactivate E. coli dihydrofolate reductase (DHFR) by introducing D27G,N,F chromosomal mutations in a key catalytic residue with subsequent adaptation by an automated serial transfer protocol. The partial reversal G27- > C occurred in three evolutionary trajectories. Conversely, in one trajectory for D27G and in all trajectories for D27F,N strains adapted to grow at very low metabolic supplement (folAmix) concentrations but did not escape entirely from supplement auxotrophy. Major global shifts in metabolome and proteome occurred upon DHFR inactivation, which were partially reversed in adapted strains. Loss-of-function mutations in two genes, thyA and deoB , ensured adaptation to low folAmix by rerouting the 2-Deoxy-D-ribose-phosphate metabolism from glycolysis towards synthesis of dTMP. Multiple evolutionary pathways of adaptation converged to a suboptimal solution due to the high accessibility to loss-of-function mutations that block the path to the highest, yet least accessible, fitness peak., Competing Interests: JR, ES No competing interests declared, (© 2019, Rodrigues and Shakhnovich.)

Background: In Gram-negative bacteria, over 30 different genes are known to encode a trimethoprim-insensitive dihydrofolate reductase that confers resistance to trimethoprim., Objectives: To determine whether a gene encoding a putative dihydrofolate reductase found in type 2 IncC plasmids isolated between 2002 and 2013 in healthcare facilities in Melbourne, Australia, confers trimethoprim resistance., Methods: Conjugation was used to transfer plasmids into a laboratory Escherichia coli. A PCR-amplified fragment was cloned into pUC19 using Gibson Assembly and transformed into E. coli. The level of resistance to trimethoprim was determined using broth microdilution. MEGA (7.0.26) and Geneious Prime (7.0.9) were used to examine the relationship to known Dfr proteins., Results: The conjugative IncC plasmid pEc158 from a 2002 Melbourne clinical E. coli isolate was shown to transfer trimethoprim resistance. The putative DfrA protein encoded by a dfrA gene in pEc158 shares <40% amino acid identity with any previously identified DfrA protein. This gene was cloned and found to confer trimethoprim resistance. The gene and protein were named dfrA35/DfrA35. In pEc158 the dfrA35 gene is located near the ori end of a partial copy of the CR1 element, within a complex resistance island. It is found in the same location in further closely-related type 2 IncC plasmids from Klebsiella pneumoniae (Melbourne, 2013), which were not transfer proficient., Conclusions: Resistance determinants continue to be found and will be missed using website-associated databases to infer phenotypes from genome sequences rather than direct phenotypic testing., (© The Author(s) 2019. Published by Oxford University Press on behalf of the British Society for Antimicrobial Chemotherapy. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Whole-genome sequencing of trimethoprim-resistant Escherichia coli strains MF2165 and PF9285 from healthy Swiss fattening calves revealed a so far uncharacterized dihydrofolate reductase gene, dfrA35 Functionality and association with trimethoprim resistance were demonstrated by cloning and expressing dfrA35 in E. coli The DfrA35 protein showed the closest amino acid identity (49.4%) to DfrA20 from Pasteurella multocida and to the Dfr determinants DfrG (41.2%), DfrD (40.8%), and DfrK (40.0%) found in Gram-positive bacteria. The dfrA35 gene was integrated within a florfenicol/chloramphenicol-sulfonamide resistance IS CR2 element ( floR -IS CR2 - dfrA35 - sul2 ) next to a Tn 21 -like transposon that contained genes with resistance to sulfonamides ( sul1 ), streptomycin ( aadA1 ), gentamicin/tobramycin/kanamycin ( aadB ), and quaternary ammonium compounds ( qacE Δ 1 ). A search of GenBank databases revealed that dfrA35 was present in 26 other E. coli strains from different origins as well as in Acinetobacter IMPORTANCE The presence of dfrA35 associated with IS CR2 in Escherichia coli from animals, as well as its presence in other E. coli strains from different sources and countries and in Acinetobacter , highlights the global spread of this gene and its potential for further dissemination. The genetic link of IS CR2 - dfrA35 with other antibiotic and disinfectant resistance genes showed that multidrug-resistant E. coli may be selected and maintained by the use of either one of several antimicrobials., (Copyright © 2019 Wüthrich et al.)

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. IMPORTANCE E. 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., (Copyright © 2018 American Society for Microbiology.)

Dihydrofolate reductase (DHFR) and 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) transformylase/IMP cyclohydrolase (PurH) play key roles in maintaining folate pools in cells, and are targets of antimicrobial and anticancer drugs. While the activities of bacterial DHFR and PurH on their classical substrates (DHF and 10-CHO-THF, respectively) are known, their activities and kinetic properties of utilisation of 10-CHO-DHF are unknown. We have determined the kinetic properties ( k cat / K m ) of conversion of 10-CHO-DHF to 10-CHO-THF by DHFR, and to DHF by PurH. We show that DHFR utilises 10-CHO-DHF about one third as efficiently as it utilises DHF. The 10-CHO-DHF is also utilised (as a formyl group donor) by PurH albeit slightly less efficiently than 10-CHO-THF. The utilisation of 10-CHO-DHF by DHFR is ~50 fold more efficient than its utilisation by PurH. A folate deficient Escherichia coli (∆ pabA ) grows well when supplemented with adenine, glycine, thymine and methionine, the metabolites that arise from the one-carbon metabolic pathway. Notably, when the ∆ pabA strain harboured a folate transporter, it grew in the presence of 10-CHO-DHF alone, suggesting that it (10-CHO-DHF) can enter one-carbon metabolic pathway to provide the required metabolites. Thus, our studies reveal that both DHFR and PurH could utilise 10-CHO-DHF for folate homeostasis in E. coli .

Ribosomes are present inside bacterial cells at micromolar concentrations and occupy up to 20% of the cell volume. Under these conditions, even weak quinary interactions between ribosomes and cytosolic proteins can affect protein activity. By using in-cell and in vitro NMR spectroscopy, and biophysical techniques, we show that the enzymes, adenylate kinase and dihydrofolate reductase, and the respective coenzymes, ATP and NADPH, bind to ribosomes with micromolar affinity, and that this interaction suppresses the enzymatic activities of both enzymes. Conversely, thymidylate synthase, which works together with dihydrofolate reductase in the thymidylate synthetic pathway, is activated by ribosomes. We also show that ribosomes impede diffusion of green fluorescent protein in vitro and contribute to the decrease in diffusion in vivo. These results strongly suggest that ribosome-mediated quinary interactions contribute to the differences between in vitro and in vivo protein activities and that ribosomes play a previously under-appreciated nontranslational role in regulating cellular biochemistry.

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.

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., (Copyright © 2017 American Society for Microbiology.)

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three-dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N-termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native-state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild-type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild-type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild-type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins., (© 2016 The Protein Society.)

Occasionally, a mutation in an exposed loop region causes a significant change in protein function and/or stability. A single mutation Gly67Val of E. coli dihydrofolate reductase (DHFR) in the exposed CD loop is such an example. We have carried out the chemical shift assignments for H(N), N(H), C(α) and C(β) atoms of the Gly67Val mutant of E. coli DHFR complexed with folate at pH 7.0, 35 °C, and then evaluated the H(N), N(H), C(α) and C(β) chemical shift changes caused by the mutation. The result indicates that, while the overall secondary structure remains the same, the single mutation Gly67Val causes site-specific conformational changes of the polypeptide backbone restricted around the adenosine-binding subdomain (residues 38-88) and not in the distant catalytic domain.

Fitness landscapes of drug resistance constitute powerful tools to elucidate mutational pathways of antibiotic escape. Here, we developed a predictive biophysics-based fitness landscape of trimethoprim (TMP) resistance for Escherichia coli dihydrofolate reductase (DHFR). We investigated the activity, binding, folding stability, and intracellular abundance for a complete set of combinatorial DHFR mutants made out of three key resistance mutations and extended this analysis to DHFR originated from Chlamydia muridarum and Listeria grayi We found that the acquisition of TMP resistance via decreased drug affinity is limited by a trade-off in catalytic efficiency. Protein stability is concurrently affected by the resistant mutants, which precludes a precise description of fitness from a single molecular trait. Application of the kinetic flux theory provided an accurate model to predict resistance phenotypes (IC50) quantitatively from a unique combination of the in vitro protein molecular properties. Further, we found that a controlled modulation of the GroEL/ES chaperonins and Lon protease levels affects the intracellular steady-state concentration of DHFR in a mutation-specific manner, whereas IC50 is changed proportionally, as indeed predicted by the model. This unveils a molecular rationale for the pleiotropic role of the protein quality control machinery on the evolution of antibiotic resistance, which, as we illustrate here, may drastically confound the evolutionary outcome. These results provide a comprehensive quantitative genotype-phenotype map for the essential enzyme that serves as an important target of antibiotic and anticancer therapies.

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.

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.

Vacuum-ultraviolet (VUV) circular dichroism (CD) spectroscopy has recently been used for secondary structure analysis of proteins; however, the contribution of aromatic side chains to protein VUV CD spectra is unresolved. In this report, VUV CD spectra of 10 Escherichia coli dihydrofolate reductase (DHFR) mutants, in which each phenylalanine or tyrosine residue was mutated to leucine, were measured down to 175 nm at 25 °C and pH 8.0 to elucidate the contributions of these aromatic side chains to the high-energy transitions of peptide bonds. The VUV CD spectra of these mutants were different from the spectrum of the wild-type protein, indicating that the contribution of the phenylalanine and tyrosine side chains of DHFR extends to the VUV region. Furthermore, the VUV CD spectrum and the folate- or NADP(+)-induced spectral change of F103L mutant DHFR indicated a modification and regeneration of exciton coupling between the Trp47 and Trp74 side chains, respectively, suggesting that exciton coupling may also contribute to the CD spectrum of DHFR in the VUV region. These results should be useful for theoretically characterizing the contribution of aromatic side chains to protein CD spectra, leading to the improvement of protein secondary-structure analysis by VUV CD spectroscopy.

Interaction between antimalarial drugs is important in determining the outcome of chemotherapy using drug combinations. Inhibitors of dihydrofolate reductase (DHFR) such as pyrimethamine and of dihydropteroate synthase (DHPS) such as sulfa drugs are known to have synergistic interactions. However, studies of the synergism are complicated by the fact that the malaria parasite can also salvage exogenous folates, and the salvage may also be affected by the drugs. It is desirable to have a convenient system to study interaction of DHFR and DHPS inhibitors without such complications. Here, we describe the use of Escherichia coli transformed with malarial DHFR and DHPS, while its own corresponding genes have been inactivated by optimal concentration of trimethoprim and genetic knockout, respectively, to study the interaction of the inhibitors. Marked synergistic effects are observed for all combinations of pyrimethamine and sulfa inhibitors in the presence of trimethoprim. At 0.05μM trimethoprim, sum of fractional inhibitory concentrations, ΣFIC of pyrimethamine with sulfadoxine, pyrimethamine with sulfathiazole, pyrimethamine with sulfamethoxazole, and pyrimethamine with dapsone are in the range of 0.24-0.41. These results show synergism between inhibitors of the two enzymes even in the absence of folate transport and uptake. This bacterial surrogate system should be useful as a tool for assessing the interactions of drug combinations between the DHFR and DHPS inhibitors., (Copyright © 2015 Elsevier B.V. All rights reserved.)

The enzyme dihydrofolate reductase (DHFR, E) from Escherichia coli is a paradigm for the role of protein dynamics in enzyme catalysis. Previous studies have shown that the enzyme progresses through the kinetic cycle by modulating the dynamic conformational landscape in the presence of substrate dihydrofolate (DHF), product tetrahydrofolate (THF), and cofactor (NADPH or NADP(+)). This study focuses on the quantitative description of the relationship between protein fluctuations and product release, the rate-limiting step of DHFR catalysis. NMR relaxation dispersion measurements of millisecond time scale motions for the E:THF:NADP(+) and E:THF:NADPH complexes of wild-type and the Leu28Phe (L28F) point mutant reveal conformational exchange between an occluded ground state and a low population of a closed state. The backbone structures of the occluded ground states of the wild-type and mutant proteins are very similar, but the rates of exchange with the closed excited states are very different. Integrated analysis of relaxation dispersion data and THF dissociation rates measured by stopped-flow spectroscopy shows that product release can occur by two pathways. The intrinsic pathway consists of spontaneous product dissociation and occurs for all THF-bound complexes of DHFR. The allosteric pathway features cofactor-assisted product release from the closed excited state and is utilized only in the E:THF:NADPH complexes. The L28F mutation alters the partitioning between the pathways and results in increased flux through the intrinsic pathway relative to the wild-type enzyme. This repartitioning could represent a general mechanism to explain changes in product release rates in other E. coli DHFR mutants.

A protein, differing in origin, may exhibit variable physicochemical behaviour, difference in sequence homology, fold and function. Thus studying structure-function relationship of proteins from altered sources is meaningful in the sense that it may give rise to comparative aspects of their sequence-structure-function relationship. Dihydrofolate reductase is an enzyme involved in cell cycle regulation. It is a significant enzyme as.a target for developing anticancer drugs. Hence, detailed understanding of structure-function relationships of wide variants of the enzyme dihydrofolate reductase would be important for developing an inhibitor or an antagonist against the enzyme involved in the cellular developmental processes. In this communication, we have reported the comparative structure-function relationship between E. coli and human dihydrofolate reductase. The differences in the unfolding behaviour of these two proteins have been investigated to understand various properties of these two proteins like relative' stability differences and variation in conformational changes under identical denaturing conditions. The equilibrium unfolding mechanism of dihydrofolate reductase proteins using guanidine hydrochloride as a denaturant in the presence of various types of osmolytes has been monitored using loss in enzymatic activity, intrinsic tryptophan fluorescence and an extrinsic fluorophore 8-anilino-1-naphthalene-sulfonic acid as probes. It has been observed that osmolytes, such as 1M sucrose, and 30% glycerol, provided enhanced stability to both variants of dihydrofolate reductase. Their level of stabilisation has been observed to be dependent on intrinsic protein stability. It was observed that 100 mM proline does not show any 'significant stabilisation to either of dihydrofolate reductases. In the present study, it has been observed that the human protein is relatively less stable than the E.coli counterpart.

Linking the molecular effects of mutations to fitness is central to understanding evolutionary dynamics. Here, we establish a quantitative relation between the global effect of mutations on the E. coli proteome and bacterial fitness. We created E. coli strains with specific destabilizing mutations in the chromosomal folA gene encoding dihydrofolate reductase (DHFR) and quantified the ensuing changes in the abundances of 2,000+ E. coli proteins in mutant strains using tandem mass tags with subsequent LC-MS/MS. mRNA abundances in the same E. coli strains were also quantified. The proteomic effects of mutations in DHFR are quantitatively linked to phenotype: the SDs of the distributions of logarithms of relative (to WT) protein abundances anticorrelate with bacterial growth rates. Proteomes hierarchically cluster first by media conditions, and within each condition, by the severity of the perturbation to DHFR function. These results highlight the importance of a systems-level layer in the genotype-phenotype relationship., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

In this study, of the hundred Escherichia coli strains isolated from feral Pigeon faeces, eighty five strains were resistant to one or more antibiotics and fifteen sensitive to all the antibiotics tested. The only strain (among all antibiotic-resistant E. coli isolates) that possessed class 1 integron was PGB01. The dihydrofolate reductase gene of the said integron was cloned, sequenced and expressed in E. coli JM109. Since PGB01 was native to pigeon's gut, we have compared the growth of PGB01 at two different temperatures, 42°C (normal body temperature of pigeon) and 37°C (optimal growth temperature of E. coli; also the human body temperature), with E. coli K12. It was found that PGB01 grew better than the laboratory strain E. coli K12 at 37°C as well as at 42°C. In the thermal fitness assay, it was observed that the cells of PGB01 were better adapted to 42°C, resembling the average body temperature of pigeon. The strain PGB01 also sustained more microwave mediated thermal stress than E. coli K12 cells. The NMR spectra of the whole cells of PGB01 varied from E. coli K12 in several spectral peaks relating some metabolic adaptation to thermotolerance. On elevating the growth temperature from 37°C to 42°C, susceptibility to kanamycin (both strains were sensitive to it) of E. coli K12 was increased, but in case of PGB01 no change in susceptibility took place. We have also attempted to reveal the basis of trimethoprim resistance phenotype conferred by the dfrA7 gene homologue of PGB01. Molecular Dynamics (MD) simulation study of docked complexes, PGB01-DfrA7 and E. coli TMP-sensitive-Dfr with trimethoprim (TMP) showed loss of some of the hydrogen and hydrophobic interaction between TMP and mutated residues in PGB01-DfrA7-TMP complex compared to TMP-sensitive-Dfr-TMP complex. This loss of interaction entails decrease in affinity of TMP for PGB01-DfrA7 compared to TMP-sensitive-Dfr.

To investigate the contribution of solvent environments to the enzymatic function of Escherichia coli dihydrofolate reductase (DHFR), the salt-, pH-, and pressure-dependence of the enzymatic function of the wild-type protein were compared with those of the active-site mutant D27E in relation to their structure and stability. The salt concentration-dependence of enzymatic activity indicated that inorganic cations bound to and inhibited the activity of wild-type DHFR at neutral pH. The BaCl2 concentration-dependence of the (1)H-(15)N HSQC spectra of the wild-type DHFR-folate binary complex showed that the cation-binding site was located adjacent to the Met20 loop. The insensitivity of the D27E mutant to univalent cations, the decreased optimal pH for its enzymatic activity, and the increased Km and Kd values for its substrate dihydrofolate suggested that the substrate-binding cleft of the mutant was slightly opened to expose the active-site side chain to the solvent. The marginally increased fluorescence intensity and decreased volume change due to unfolding of the mutant also supported this structural change or the modified cavity and hydration. Surprisingly, the enzymatic activity of the mutant increased with pressurization up to 250MPa together with negative activation volumes of -4.0 or -4.8mL/mol, depending on the solvent system, while that of the wild-type was decreased and had positive activation volumes of 6.1 or 7.7mL/mol. These results clearly indicate that the insertion of a single methylene at the active site could substantially change the enzymatic reaction mechanism of DHFR, and solvent environments play important roles in the function of this enzyme., (© 2013.)

Although allosteric effector antibodies are used widely as modulators of receptors and enzymes, experimental analysis of their mechanism remains highly challenging. Here, we investigate the molecular mechanisms of allosteric and non-allosteric effector antibodies in an experimentally tractable system, consisting of single-domain antibodies (nanobodies) that target the model enzyme dihydrofolate reductase (DHFR) from Escherichia coli. A panel of thirty-five nanobodies was isolated using several strategies to increase nanobody diversity. The nanobodies exhibit a variety of effector properties, including partial inhibition, strong inhibition and stimulation of DHFR activity. Despite these diverse effector properties, chemical shift perturbation NMR epitope mapping identified only two epitope regions: epitope α is a new allosteric site that is over 10Å from the active site, while epitope β is located in the region of the Met20 loop. The structural basis for DHFR allosteric inhibition or activation upon nanobody binding to the α epitope was examined by solving the crystal structures of DHFR in complex with Nb113 (an allosteric inhibitor) and Nb179 (an allosteric activator). The structures suggest roles for conformational constraint and altered protein dynamics, but not epitope distortion, in the observed allosteric effects. The crystal structure of a β epitope region binder (ca1698) in complex with DHFR is also reported. Although CDR3 of ca1698 occupies the substrate binding site, ca1698 displays linear mixed inhibition kinetics instead of simple competitive inhibition kinetics. Two mechanisms are proposed to account for this apparent anomaly. Evidence for structural convergence of ca1698 and Nb216 during affinity maturation is also presented., (© 2013.)

With the rapidly growing wealth of genomic data, experimental inquiries on the functional significance of important divergence sites in protein evolution are becoming more accessible. Here we trace the evolution of dihydrofolate reductase (DHFR) and identify multiple key divergence sites among 233 species between humans and bacteria. We connect these sites, experimentally and computationally, to changes in the enzyme's binding properties and catalytic efficiency. One of the identified evolutionarily important sites is the N23PP modification (∼mid-Devonian, 415-385 Mya), which alters the conformational states of the active site loop in Escherichia coli dihydrofolate reductase and negatively impacts catalysis. This enzyme activity was restored with the inclusion of an evolutionarily significant lid domain (G51PEKN in E. coli enzyme; ∼2.4 Gya). Guided by this evolutionary genomic analysis, we generated a human-like E. coli dihydrofolate reductase variant through three simple mutations despite only 26% sequence identity between native human and E. coli DHFRs. Molecular dynamics simulations indicate that the overall conformational motions of the protein within a common scaffold are retained throughout evolution, although subtle changes to the equilibrium conformational sampling altered the free energy barrier of the enzymatic reaction in some cases. The data presented here provide a glimpse into the evolutionary trajectory of functional DHFR through its protein sequence space that lead to the diverged binding and catalytic properties of the E. coli and human enzymes.

What are the molecular properties of proteins that fall on the radar of protein quality control (PQC)? Here we mutate the E. coli's gene encoding dihydrofolate reductase (DHFR) and replace it with bacterial orthologous genes to determine how components of PQC modulate fitness effects of these genetic changes. We find that chaperonins GroEL/ES and protease Lon compete for binding to molten globule intermediate of DHFR, resulting in a peculiar symmetry in their action: overexpression of GroEL/ES and deletion of Lon both restore growth of deleterious DHFR mutants and most of the slow-growing orthologous DHFR strains. Kinetic steady-state modeling predicts and experimentation verifies that mutations affect fitness by shifting the flux balance in cellular milieu between protein production, folding, and degradation orchestrated by PQC through the interaction with folding intermediates., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Low protein solubility is a problem in many areas of protein science. Although chemical methods have been developed to solubilize proteins these are not always effective and add to the cost of producing the protein. One way of overcoming these difficulties is to evolve the protein to be more soluble. A major hurdle in this process is the ability to select mutant proteins with enhanced solubility from a large library of randomly mutated proteins. In this article, we describe such a method. The method relies on the fact that increasing the expression of dihydrofolate reductase (DHFR) makes Escherichia coli resistant to Trimethoprim (TMP). Proteins fused to DHFR will produce chimeras with altered levels of resistance to TMP. This variation in TMP resistance can be used to identify mutant proteins with enhanced solubility.

The micro-viscosity and molecular crowding experienced by specific proteins can regulate their dynamics and function within live cells. Taking advantage of the emerging TMP-tag technology, we present the design, synthesis and application of a hybrid genetic-chemical molecular rotor probe whose fluorescence lifetime can report protein-specific micro-environments in live cells.

Mutations create the genetic diversity on which selective pressures can act, yet also create structural instability in proteins. How, then, is it possible for organisms to ameliorate mutation-induced perturbations of protein stability while maintaining biological fitness and gaining a selective advantage? Here we used site-specific chromosomal mutagenesis to introduce a selected set of mostly destabilizing mutations into folA--an essential chromosomal gene of Escherichia coli encoding dihydrofolate reductase (DHFR)--to determine how changes in protein stability, activity, and abundance affect fitness. In total, 27 E. coli strains carrying mutant DHFR were created. We found no significant correlation between protein stability and its catalytic activity nor between catalytic activity and fitness in a limited range of variation of catalytic activity observed in mutants. The stability of these mutants is strongly correlated with their intracellular abundance, suggesting that protein homeostatic machinery plays an active role in maintaining intracellular concentrations of proteins. Fitness also shows a significant correlation with intracellular abundance of soluble DHFR in cells growing at 30 °C. At 42 °C, the picture was mixed, yet remarkable: A few strains carrying mutant DHFR proteins aggregated, rendering them nonviable, but, intriguingly, the majority exhibited fitness higher than wild type. We found that mutational destabilization of DHFR proteins in E. coli is counterbalanced at 42 °C by their soluble oligomerization, thereby restoring structural stability and protecting against aggregation.

The question of whether protein motions play a role in the chemical step of enzymatic catalysis has generated much controversy in recent years. Debate has recently reignited over possible dynamic contributions to catalysis in dihydrofolate reductase, following conflicting conclusions from studies of the N23PP/S148A variant of the Escherichia coli enzyme. By investigating the temperature dependence of kinetic isotope effects, we present evidence that the reduction in the hydride transfer rate constants in this variant is not a direct result of impairment of conformational fluctuations. Instead, the conformational state of the enzyme immediately before hydride transfer, which determines the electrostatic environment of the active site, affects the rate constant for the reaction. Although protein motions are clearly important for binding and release of substrates and products, there appears to be no detectable dynamic coupling of protein motions to the hydride transfer step itself.

A significant contemporary question in enzymology involves the role of protein dynamics and hydrogen tunneling in enhancing enzyme catalyzed reactions. Here, we report a correlation between the donor-acceptor distance (DAD) distribution and intrinsic kinetic isotope effects (KIEs) for the dihydrofolate reductase (DHFR) catalyzed reaction. This study compares the nature of the hydride-transfer step for a series of active-site mutants, where the size of a side chain that modulates the DAD (I14 in E. coli DHFR) is systematically reduced (I14V, I14A, and I14G). The contributions of the DAD and its dynamics to the hydride-transfer step were examined by the temperature dependence of intrinsic KIEs, hydride-transfer rates, activation parameters, and classical molecular dynamics (MD) simulations. Results are interpreted within the framework of the Marcus-like model where the increase in the temperature dependence of KIEs arises as a direct consequence of the deviation of the DAD from its distribution in the wild type enzyme. Classical MD simulations suggest new populations with larger average DADs, as well as broader distributions, and a reduction in the population of the reactive conformers correlated with the decrease in the size of the hydrophobic residue. The more flexible active site in the mutants required more substantial thermally activated motions for effective H-tunneling, consistent with the hypothesis that the role of the hydrophobic side chain of I14 is to restrict the distribution and dynamics of the DAD and thus assist the hydride-transfer. These studies establish relationships between the distribution of DADs, the hydride-transfer rates, and the DAD's rearrangement toward tunneling-ready states. This structure-function correlation shall assist in the interpretation of the temperature dependence of KIEs caused by mutants far from the active site in this and other enzymes, and may apply generally to C-H→C transfer reactions., (© 2011 American Chemical Society)

Recent work indicates a general architecture for proteins in which sparse networks of physically contiguous and coevolving amino acids underlie basic aspects of structure and function. These networks, termed sectors, are spatially organized such that active sites are linked to many surface sites distributed throughout the structure. Using the metabolic enzyme dihydrofolate reductase as a model system, we show that: (1) the sector is strongly correlated to a network of residues undergoing millisecond conformational fluctuations associated with enzyme catalysis, and (2) sector-connected surface sites are statistically preferred locations for the emergence of allosteric control in vivo. Thus, sectors represent an evolutionarily conserved "wiring" mechanism that can enable perturbations at specific surface positions to rapidly initiate conformational control over protein function. These findings suggest that sectors enable the evolution of intermolecular communication and regulation., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Dihydrofolate reductase (DHFR) is a well-studied drug target and a paradigm for understanding enzyme catalysis. Preparation of pure DHFR samples, in defined ligand-bound states, is a prerequisite for in vitro studies and drug discovery efforts. We use NMR spectroscopy to monitor ligand content of human and Escherichia coli DHFR (ecDHFR), which bind different co-purifying ligands during expression in bacteria. An alternate purification strategy yields highly pure DHFR complexes, containing only the desired ligands, in the quantities required for structural studies. Interestingly, ecDHFR is bound to endogenous THF while human DHFR is bound to NADP. Consistent with these findings, a designed "humanized" mutant of ecDHFR switches binding specificity in the cell., (Copyright © 2011 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

The characterization of microsecond dynamics in the folding of multisubdomain proteins has been a major challenge in understanding their often complex folding mechanisms. Using a continuous-flow mixing device coupled with fluorescence lifetime detection, we report the microsecond folding dynamics of dihydrofolate reductase (DHFR), a two-subdomain α/β/α sandwich protein known to begin folding in this time range. The global dimensions of early intermediates were monitored by Förster resonance energy transfer, and the dynamic properties of the local Trp environments were monitored by fluorescence lifetime detection. We found that substantial collapse occurs in both the locally connected adenosine binding subdomain and the discontinuous loop subdomain within 35 μs of initiation of folding from the urea unfolded state. During the fastest observable ∼550 μs phase, the discontinuous loop subdomain further contracts, concomitant with the burial of Trp residue(s), as both subdomains achieve a similar degree of compactness. Taken together with previous studies in the millisecond time range, a hierarchical assembly of DHFR--in which each subdomain independently folds, subsequently docks, and then anneals into the native conformation after an initial heterogeneous global collapse--emerges. The progressive acquisition of structure, beginning with a continuously connected subdomain and spreading to distal regions, shows that chain entropy is a significant organizing principle in the folding of multisubdomain proteins and single-domain proteins. Subdomain folding also provides a rationale for the complex kinetics often observed., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Conformational dynamics play a key role in enzyme catalysis. Although protein motions have clear implications for ligand flux, a role for dynamics in the chemical step of enzyme catalysis has not been clearly established. We generated a mutant of Escherichia coli dihydrofolate reductase that abrogates millisecond-time-scale fluctuations in the enzyme active site without perturbing its structural and electrostatic preorganization. This dynamic knockout severely impairs hydride transfer. Thus, we have found a link between conformational fluctuations on the millisecond time scale and the chemical step of an enzymatic reaction, with broad implications for our understanding of enzyme mechanisms and for design of novel protein catalysts.

Post-translational glycosylation is one of the most abundant forms of covalent protein modification in eukaryotic cells. It plays an important role in determining the properties of proteins, and stabilizes many proteins against thermal denaturation. Protein glycosylation may establish a surface microenvironment that resembles that of unglycosylated proteins in concentrated solutions of sugars and other polyols. We have used site-directed mutagenesis to introduce a series of unique cysteine residues into a cysteine-free double mutant (DM, C85A/C152S) of dihydrofolate reductase from Escherichia coli (EcDHFR). The resulting triple mutants, DM-N18C, DM-R52C, DM-D87C and DM-D132C EcDHFR, were alkylated with glucose, N-acetylglucosamine, lactose and maltotriose iodoacetamides. We found little effect on catalysis or stability in three cases. However, when DM-D87C EcDHFR is glycosylated, stability is increased by between 1.5 and 2.6 kcal.mol(-1) in a sugar-dependent manner. D87 is found in a hinge region of EcDHFR that loses structure early in the thermal denaturation process, whereas the other glycosylation sites are found in regions involved in the later stages of temperature-induced unfolding. Glycosylation at this site may improve the stability of EcDHFR by protecting a region of the enzyme that is particularly prone to denaturation.

Residues distal from the active site in dihydrofolate reductase (DHFR) have regulatory roles in catalytic reaction and also folding stability. The couplings of the distal residues to the ones in the active site have been analyzed using site-directed mutants. To expand our understanding of the structural and functional influences of distal residue mutation, we explored the structural stability and enzymatic activity of deletion mutants. Deletion has greater structural and dynamical impacts on the corresponding part than site-directed mutation does. Thus, deletion amplifies the effects caused by distal mutations, which should make the mutual couplings among the distant residues more apparent. We focused on residues 52, 67, 121, and 145 in the four distinct loops of DHFR. All the single-residue deletion mutants showed marked reduction in stability, except for Delta52 in an alphaC-betaC loop. Double deletion mutants showed that the loop alphaC-betaC has nonadditive couplings with the betaF-betaG and betaG-betaH loops regarding stability. Single deletion to the loops alphaC-betaC or betaC-betaD resulted in considerable activity reduction, demonstrating that the loops couple to the residues near the active site. The four loops were shown to be functionally interdependent from the double deletion experiments., (Copyright 2009 Elsevier B.V. All rights reserved.)

Background: Trimethoprim resistance is increasing in Enterobacteriaceae. In 2004-2006 an intervention on trimethoprim use was conducted in Kronoberg County, Sweden, resulting in 85% reduction in trimethoprim prescriptions. We investigated the distribution of dihydrofolate reductase (dfr)-genes and integrons in Escherichia coli and Klebsiella pneumoniae and the effect of the intervention on this distribution., Methodology/principal Findings: Consecutively isolated E. coli (n = 320) and K. pneumoniae (n = 54) isolates phenotypically resistant to trimethoprim were studied. All were investigated for the presence of dfrA1, dfrA5, dfrA7, dfrA8, dfrA12, dfrA14, dfrA17 and integrons class I and II. Isolates negative for the seven dfr-genes (n = 12) were also screened for dfr2d, dfrA3, dfrA9, dfrA10, dfrA24 and dfrA26. These genes accounted for 96% of trimethoprim resistance in E. coli and 69% in K. pneumoniae. The most prevalent was dfrA1 in both species. This was followed by dfrA17 in E. coli which was only found in one K. pneumoniae isolate. Class I and II Integrons were more common in E. coli (85%) than in K. pneumoniae (57%). The distribution of dfr-genes did not change during the course of the 2-year intervention., Conclusions/significance: The differences observed between the studied species in terms of dfr-gene and integron prevalence indicated a low rate of dfr-gene transfer between these two species and highlighted the possible role of narrow host range plasmids in the spread of trimethoprim resistance. The stability of dfr-genes, despite large changes in the selective pressure, indirectly suggests a low fitness cost of dfr-gene carriage.

Sequence perturbation analysis is a powerful method to reveal roles of an amino acid residue in function and stability of a protein. By using and improving this method, we studied roles of highly conserved Arg44 of Escherichia coli dihydrofolate reductase (DHFR) in its function and stability. Here, we introduced systematic amino acid substitutions at this position and found that all 19 kinds of amino acid substitutions were tolerated, but the mutations significantly reduced the enzymatic activity and the binding affinity toward the cofactor NADPH. Moreover, the mutational effects on the cofactor binding affinity were well correlated with those on the catalytic activity, indicating that the R44X mutations affect the catalytic activity mainly by modulating the cofactor binding affinity. On the other hand, thermal denaturation measurements showed that most mutations stabilized the protein. Comparison between the mutational effects and various amino acid indices taken from the AAindex database indicated that hydrophobicity and polarity are key determinants of amino acids favorable at this position. These results suggest that through electrostatic interactions Arg44 plays a functional role in retaining the cofactor binding affinity at the cost of the protein stability., (Copyright 2009 Elsevier Inc. All rights reserved.)

Tetrahydromonapterin is a major pterin in Escherichia coli and is hypothesized to be the cofactor for phenylalanine hydroxylase (PhhA) in Pseudomonas aeruginosa, but neither its biosynthetic origin nor its cofactor role has been clearly demonstrated. A comparative genomics analysis implicated the enigmatic folX and folM genes in tetrahydromonapterin synthesis via their phyletic distribution and chromosomal clustering patterns. folX encodes dihydroneopterin triphosphate epimerase, which interconverts dihydroneopterin triphosphate and dihydromonapterin triphosphate. folM encodes an unusual short-chain dehydrogenase/reductase known to have dihydrofolate and dihydrobiopterin reductase activity. The roles of FolX and FolM were tested experimentally first in E. coli, which lacks PhhA and in which the expression of P. aeruginosa PhhA plus the recycling enzyme pterin 4a-carbinolamine dehydratase, PhhB, rescues tyrosine auxotrophy. This rescue was abrogated by deleting folX or folM and restored by expressing the deleted gene from a plasmid. The folX deletion selectively eliminated tetrahydromonapterin production, which far exceeded folate production. Purified FolM showed high, NADPH-dependent dihydromonapterin reductase activity. These results were substantiated in P. aeruginosa by deleting tyrA (making PhhA the sole source of tyrosine) and folX. The DeltatyrA strain was, as expected, prototrophic for tyrosine, whereas the DeltatyrA DeltafolX strain was auxotrophic. As in E. coli, the folX deletant lacked tetrahydromonapterin. Collectively, these data establish that tetrahydromonapterin formation requires both FolX and FolM, that tetrahydromonapterin is the physiological cofactor for PhhA, and that tetrahydromonapterin can outrank folate as an end product of pterin biosynthesis.

Many of the antibiotic resistance elements used in molecular biology have idiosyncratic limitations. For example, beta-lactam selections rely on antibiotics that are unstable to hydrolysis and allow satellite colonies to form upon extended incubation, and tetracycline selections typically give rise to widely varying colony sizes and lower transformation efficiencies. Although prokaryotic Type II dihydrofolate reductase (dfr) genes have long been considered to have potential utility for the selection of plasmids and mobile elements in bacteria, practical limitations to the quality of those selections, mostly relating to background and inefficiency, have led for the most part to their underuse. I describe here the construction of a Type IIa dfr prokaryotic expression cassette that confers strong resistance against trimethoprim (Tmp), a bactericidal dfr inhibiting antibiotic. The Tmp-resistance cassette provides consistent and efficient selections and plasmid transformation frequencies equivalent to those encountered with beta-lactamases.