The crystal structure of dihydrofolate reductase from Thermotoga maritima: molecular features of thermostability

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

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

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

Integron mediated antimicrobial resistance in diarrheagenic Escherichia coli in children: in vitro and in silico analysis

The exchange of genes between bacterial chromosome and plasmid(s) and their integration into integrons are mainly responsible for acquisition and dissemination of antibiotic resistance. We investigated the role of integrons and their underlying molecular mechanisms leading to development of adaptability in E. coli and eventual resistance to antimicrobials. Escherichia coli isolates (n = 120); including 40 diarrheagenic isolates, an even number of isolates from cases other than diarrhea, and equal number of isolates from healthy children recovered from fresh stool samples were used for identification of integron genes and gene cassettes. The association of integrons with antibiotic resistance was assayed before phylogenetic analysis. DNA sequence analysis revealed class 1 and 2 integrons in 55.83% and 21.66% isolates, respectively. The integron presence was found significantly associated with the probability of antibiotic resistance in E. coli; the association being highest with class 1 integron. Modelling and molecular docking along with molecular dynamics simulation analyses found ceftriaxone and amoxicillin as potential inhibitors of dihydrofolate reductase (DHFR). The class 1 integrons of these pathogenic isolates can serve as prospective therapeutic targets using specific silencing strategies and combinational antimicrobial therapy. The findings may be useful for the development of a potent and versatile drug for DHFR inhibition.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

Sequence and Structure Properties Uncover the Natural Classification of Protein Complexes Formed by Intrinsically Disordered Proteins via Mutual Synergistic Folding

Intrinsically disordered proteins mediate crucial biological functions through their interactions with other proteins. Mutual synergistic folding (MSF) occurs when all interacting proteins are disordered, folding into a stable structure in the course of the complex formation. In these cases, the folding and binding processes occur in parallel, lending the resulting structures uniquely heterogeneous features. Currently there are no dedicated classification approaches that take into account the particular biological and biophysical properties of MSF complexes. Here, we present a scalable clustering-based classification scheme, built on redundancy-filtered features that describe the sequence and structure properties of the complexes and the role of the interaction, which is directly responsible for structure formation. Using this approach, we define six major types of MSF complexes, corresponding to biologically meaningful groups. Hence, the presented method also shows that differences in binding strength, subcellular localization, and regulation are encoded in the sequence and structural properties of proteins. While current protein structure classification methods can also handle complex structures, we show that the developed scheme is fundamentally different, and since it takes into account defining features of MSF complexes, it serves as a better representation of structures arising through this specific interaction mode.

Why are some Enzymes Dimers? Flexibility and Catalysis in Thermotoga Maritima Dihydrofolate Reductase

Dihydrofolate reductase from Thermotoga maritima (TmDFHFR) is a dimeric thermophilic enzyme that catalyzes the hydride transfer from the cofactor NADPH to dihydrofolate less efficiently than other DHFR enzymes, such as the mesophilic analogue Escherichia coli DHFR (EcDHFR). Using QM/MM potentials, we show that the reduced catalytic efficiency of TmDHFR is most likely due to differences in the amino acid sequence that stabilize the M20 loop in an open conformation, which prevents the formation of some interactions in the transition state and increases the number of water molecules in the active site. However, dimerization provides two advantages to the thermophilic enzyme: it protects its structure against denaturation by reducing thermal fluctuations and it provides a less negative activation entropy, toning down the increase of the activation free energy with temperature. Our molecular picture is confirmed by the analysis of the temperature dependence of enzyme kinetic isotope effects in different DHFR enzymes.

High-Order Epistasis in Catalytic Power of Dihydrofolate Reductase Gives Rise to a Rugged Fitness Landscape in the Presence of Trimethoprim Selection

Evolutionary fitness landscapes of several antibiotic target proteins have been comprehensively mapped showing strong high-order epistasis between mutations, but understanding these effects at the biochemical and structural levels remained open. Here, we carried out an extensive experimental and computational study to quantitatively understand the evolutionary dynamics of Escherichia coli dihydrofolate reductase (DHFR) enzyme in the presence of trimethoprim-induced selection. To facilitate this, we developed a new in vitro assay for rapidly characterizing DHFR steady-state kinetics. Biochemical and structural characterization of resistance-conferring mutations targeting a total of ten residues spanning the substrate binding pocket of DHFR revealed distinct changes in the catalytic efficiencies of mutated DHFR enzymes. Next, we measured biochemical parameters (Km, Ki, and kcat) for a mutant library carrying all possible combinations of six resistance-conferring DHFR mutations and quantified epistatic interactions between them. We found that the high-order epistasis in catalytic power of DHFR (kcat and Km) creates a rugged fitness landscape under trimethoprim selection. Taken together, our data provide a concrete illustration of how epistatic coupling at the level of biochemical parameters can give rise to complex fitness landscapes, and suggest new strategies for developing mutant specific inhibitors.

Increased substrate affinity in the Escherichia coli L28R dihydrofolate reductase mutant causes trimethoprim resistance

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.

A Rapid Analysis of Variations in Conformational Behavior during Dihydrofolate Reductase Catalysis

Protein flexibility is central to enzyme catalysis, yet it remains challenging both to predict conformational behavior on the basis of analysis of amino acid sequence and protein structure and to provide the necessary breadth of experimental support to any such predictions. Here a generic and rapid procedure for identifying conformational changes during dihydrofolate reductase (DHFR) catalysis is described. Using DHFR from Escherichia coli (EcDHFR), selective side-chain ¹³C labeling of methionine and tryptophan residues is shown to be sufficient to detect the closed-to-occluded conformational transition that follows the chemical step in the catalytic cycle, with clear chemical shift perturbations found for both methionine methyl and tryptophan indole groups. In contrast, no such perturbations are seen for the DHFR from the psychrophile Moritella profunda, where the equivalent conformational change is absent. Like EcDHFR, Salmonella enterica DHFR shows experimental evidence of a large-scale conformational change following hydride transfer that relies on conservation of a key hydrogen bonding interaction between the M20 and GH loops, directly comparable to the closed-to-occluded conformational change observed in EcDHFR. For the hyperthermophile Thermotoga maritima, no chemical shift perturbations were observed, suggesting that no major conformational change occurs during the catalytic cycle. In spite of their conserved tertiary structures, DHFRs display variations in conformational sampling that occurs concurrently with catalysis.

Reduction of Folate by Dihydrofolate Reductase from Thermotoga maritima

Mammalian dihydrofolate reductases (DHFRs) catalyze the reduction of folate more efficiently than the equivalent bacterial enzymes do, despite typically having similar efficiencies for the reduction of their natural substrate, dihydrofolate. In contrast, we show here that DHFR from the hyperthermophilic bacterium Thermotoga maritima can catalyze reduction of folate to tetrahydrofolate with an efficiency similar to that of reduction of dihydrofolate under saturating conditions. Nuclear magnetic resonance and mass spectrometry experiments showed no evidence of the production of free dihydrofolate during either the EcDHFR- or TmDHFR-catalyzed reductions of folate, suggesting that both enzymes perform the two reduction steps without release of the partially reduced substrate. Our results imply that the reaction proceeds more efficiently in TmDHFR than in EcDHFR because the more open active site of TmDHFR facilitates protonation of folate. Because T. maritima lives under extreme conditions where tetrahydrofolate is particularly prone to oxidation, this ability to salvage folate may impart an advantage to the bacterium by minimizing the squandering of a valuable cofactor.

Protein motions and dynamic effects in enzyme 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.

Minimization of dynamic effects in the evolution of dihydrofolate reductase

Protein isotope labeling is a powerful technique to probe functionally important motions in enzyme catalysis and can be applied to investigate the conformational dynamics of proteins. Previous investigations have indicated that dynamic coupling is detrimental to catalysis by dihydrofolate reductase (DHFR) from the mesophile Escherichia coli (EcDHFR). Comparison of DHFRs from organisms adapted to survive at a wide range of temperatures suggests that dynamic coupling in DHFR catalysis has been minimized during evolution; it arises from reorganizational motions needed to facilitate charge transfer events. Contrary to the behaviour observed for the DHFR from the moderate thermophile Geobacillus stearothermophilus (BsDHFR), the chemical transformation catalyzed by the cold-adapted bacterium Moritella profunda (MpDHFR) is only weakly affected by protein isotope substitutions at low temperatures, but the isotopically substituted enzyme is a substantially inferior catalyst at higher, non-physiological temperatures. QM/MM studies revealed that this behaviour is caused by the enzyme's structural sensitivity to temperature changes, which enhances unfavorable dynamic coupling at higher temperatures by promoting additional recrossing trajectories on the transition state dividing surface. We postulate that these motions are minimized by fine-tuning DHFR flexibility through optimization of the free energy surface of the reaction, such that a nearly static reaction-ready configuration with optimal electrostatic properties is maintained under physiological conditions.

Structure and Stability of the Dimeric Triosephosphate Isomerase from the Thermophilic Archaeon Thermoplasma acidophilum

Thermoplasma acidophilum is a thermophilic archaeon that uses both non-phosphorylative Entner-Doudoroff (ED) pathway and Embden-Meyerhof-Parnas (EMP) pathway for glucose degradation. While triosephosphate isomerase (TPI), a well-known glycolytic enzyme, is not involved in the ED pathway in T. acidophilum, it has been considered to play an important role in the EMP pathway. Here, we report crystal structures of apo- and glycerol-3-phosphate-bound TPI from T. acidophilum (TaTPI). TaTPI adopts the canonical TIM-barrel fold with eight α-helices and parallel eight β-strands. Although TaTPI shares ~30% sequence identity to other TPIs from thermophilic species that adopt tetrameric conformation for enzymatic activity in their harsh physiological environments, TaTPI exists as a dimer in solution. We confirmed the dimeric conformation of TaTPI by analytical ultracentrifugation and size-exclusion chromatography. Helix 5 as well as helix 4 of thermostable tetrameric TPIs have been known to play crucial roles in oligomerization, forming a hydrophobic interface. However, TaTPI contains unique charged-amino acid residues in the helix 5 and adopts dimer conformation. TaTPI exhibits the apparent T_d value of 74.6°C and maintains its overall structure with some changes in the secondary structure contents at extremely acidic conditions (pH 1–2). Based on our structural and biophysical analyses of TaTPI, more compact structure of the protomer with reduced length of loops and certain patches on the surface could account for the robust nature of Thermoplasma acidophilum TPI.

Environmental Adaptation of Dihydrofolate Reductase from Deep-Sea Bacteria

In order to elucidate the molecular adaptation mechanisms of enzymes to the high hydrostatic pressure of the deep sea, we cloned, purified, and characterized more than ten dihydrofolate reductases (DHFRs) from bacteria living in deep-sea and ambient atmospheric pressure environments. The nucleotide and amino acid sequences of these DHFRs indicate the deep-sea bacteria are adapted to their environments after the differentiation of their genus from ancestors inhabiting atmospheric pressure environments. In particular, the backbone structure of the deep-sea DHFR from Moritella profunda (mpDHFR) almost overlapped with the normal homolog from Escherichia coli (ecDHFR). Thus, those of other DHFRs would also overlap on the basis of their sequence similarities. However, the structural stability of both DHFRs was quite different: compared to ecDHFR, mpDHFR was more thermally stable but less stable against urea and pressure unfolding. The smaller volume changes due to unfolding suggest that the native structure of mpDHFR has a smaller cavity and/or enhanced hydration compared to ecDHFR. High hydrostatic pressure reduced the enzymatic activity of many DHFRs, but three deep-sea DHFRs and the D27E mutant of ecDHFR exhibited pressure-dependent activation. The inverted activation volumes from positive to negative values indicate the modification of their structural dynamics, conversion of the rate-determining step of the enzymatic reaction, and different contributions of the cavity and hydration to the transition-state structure. Since the cavity and hydration depend on amino acid side chains, DHFRs would adapt to the deep-sea environment by regulating the cavity and hydration by substituting their amino acid side chains without altering their backbone structure. The results of this study clearly indicate that the cavity and hydration play important roles in the adaptation of enzymes to the deep-sea environment.

Protein isotope effects in dihydrofolate reductase from Geobacillus stearothermophilus show entropic-enthalpic compensatory effects on the rate constant

Catalysis by dihydrofolate reductase from the moderately thermophilic bacterium Geobacillus stearothermophilus (BsDHFR) was investigated by isotope substitution of the enzyme. The enzyme kinetic isotope effect for hydride transfer was close to unity at physiological temperatures but increased with decreasing temperatures to a value of 1.65 at 5 °C. This behavior is opposite to that observed for DHFR from Escherichia coli (EcDHFR), where the enzyme kinetic isotope effect increased slightly with increasing temperature. These experimental results were reproduced in the framework of variational transition-state theory that includes a dynamical recrossing coefficient that varies with the mass of the protein. Our simulations indicate that BsDHFR has greater flexibility than EcDHFR on the ps-ns time scale, which affects the coupling of the environmental motions of the protein to the chemical coordinate and consequently to the recrossing trajectories on the reaction barrier. The intensity of the dynamic coupling in DHFRs is influenced by compensatory temperature-dependent factors, namely the enthalpic barrier needed to achieve an ideal transition-state configuration with minimal nonproductive trajectories and the protein disorder that disrupts the electrostatic preorganization required to stabilize the transition state. Together with our previous studies of other DHFRs, the results presented here provide a general explanation why protein dynamic effects vary between enzymes. Our theoretical treatment demonstrates that these effects can be satisfactorily reproduced by including a transmission coefficient in the rate constant calculation, whose dependence on temperature is affected by the protein flexibility.

Structure and evolution of the spliceosomal peptidyl-prolyl cis-trans isomerase Cwc27

Cwc27 is a spliceosomal cyclophilin-type peptidyl-prolyl cis-trans isomerase (PPIase). Here, the crystal structure of a relatively protease-resistant N-terminal fragment of human Cwc27 containing the PPIase domain was determined at 2.0 Å resolution. The fragment exhibits a C-terminal appendix and resides in a reduced state compared with the previous oxidized structure of a similar fragment. By combining multiple sequence alignments spanning the eukaryotic tree of life and secondary-structure prediction, Cwc27 proteins across the entire eukaryotic kingdom were identified. This analysis revealed the specific loss of a crucial active-site residue in higher eukaryotic Cwc27 proteins, suggesting that the protein evolved from a prolyl isomerase to a pure proline binder. Noting a fungus-specific insertion in the PPIase domain, the 1.3 Å resolution crystal structure of the PPIase domain of Cwc27 from Chaetomium thermophilum was also determined. Although structurally highly similar in the core domain, the C. thermophilum protein displayed a higher thermal stability than its human counterpart, presumably owing to the combined effect of several amino-acid exchanges that reduce the number of long side chains with strained conformations and create new intramolecular interactions, in particular increased hydrogen-bond networks.

Lean forward: Genetic analysis of temperature-sensitive mutants unfolds the secrets of oligomeric protein complex assembly

Multisubunit protein complexes are essential for cellular function. Genetic analysis of essential processes requires special tools, among which temperature-sensitive (Ts) mutants have historically been crucial. Many researchers assume that the effect of temperature on such mutants is to drive their proteolytic destruction. In fact, degradation-mediated elimination of mutant proteins likely explains only a fraction of the phenotypes associated with Ts mutants. Here I discuss insights gained from analysis of Ts mutants in oligomeric proteins, with particular focus on the study of septins, GTP-binding subunits of cytoskeletal filaments whose structures and functions are the subject of current investigation in my and many other labs. I argue that the kinds of unbiased forward genetic approaches that generate Ts mutants provide information that is largely inaccessible to modern reverse genetic methodologies, and will continue to drive our understanding of higher-order assembly by septins and other oligomeric proteins.

Different dynamical effects in mesophilic and hyperthermophilic dihydrofolate reductases

The role of protein dynamics in the reaction catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) has been examined by enzyme isotope substitution ((15)N, (13)C, (2)H). In contrast to all other enzyme reactions investigated previously, including DHFR from Escherichia coli (EcDHFR), for which isotopic substitution led to decreased reactivity, the rate constant for the hydride transfer step is not affected by isotopic substitution of TmDHFR. TmDHFR therefore appears to lack the coupling of protein motions to the reaction coordinate that have been identified for EcDHFR catalysis. Clearly, dynamical coupling is not a universal phenomenon that affects the efficiency of enzyme catalysis.

Thermal adaptation of dihydrofolate reductase from the moderate thermophile Geobacillus stearothermophilus

The thermal melting temperature of dihydrofolate reductase from Geobacillus stearothermophilus (BsDHFR) is ~30 °C higher than that of its homologue from the psychrophile Moritella profunda. Additional proline residues in the loop regions of BsDHFR have been proposed to enhance the thermostability of BsDHFR, but site-directed mutagenesis studies reveal that these proline residues contribute only minimally. Instead, the high thermal stability of BsDHFR is partly due to removal of water-accessible thermolabile residues such as glutamine and methionine, which are prone to hydrolysis or oxidation at high temperatures. The extra thermostability of BsDHFR can be obtained by ligand binding, or in the presence of salts or cosolvents such as glycerol and sucrose. The sum of all these incremental factors allows BsDHFR to function efficiently in the natural habitat of G. stearothermophilus, which is characterized by temperatures that can reach 75 °C.

Structure of diaminohydroxyphosphoribosylaminopyrimidine deaminase/5-amino-6-(5-phosphoribosylamino)uracil reductase from Acinetobacter baumannii

The bifunctional diaminohydroxyphosphoribosylaminopyrimidine deaminase/5-amino-6-(5-phosphoribosylamino)uracil reductase (RibD) represents a potential antibacterial drug target. The structure of recombinant Acinetobacter baumannii RibD is reported in orthorhombic and tetragonal crystal forms at 2.2 and 2.0 Å resolution, respectively. Comparisons with orthologous structures in the Protein Data Bank indicated close similarities. The tetragonal crystal form was obtained in the presence of guanosine monophosphate, which surprisingly was observed to occupy the adenine-binding site of the reductase domain.

Interaction of dihydrofolate reductase and aminoglycoside adenyltransferase enzyme from Klebsiella pneumoniae multidrug resistant strain DF12SA with clindamycin: a molecular modelling and docking study

Klebsiella pneumoniae strain DF12SA (HQ114261) was isolated from diabetic foot wounds. The strain showed resistance against ampicillin, kanamycin, gentamicin, streptomycin, spectinomycin, trimethoprim, tetracycline, meropenem, amikacin, piperacillin/tazobactam, augmentin, co-trimoxazole, carbapenems, penicillins and cefoperazone, and was sensitive to clindamycin. Molecular characterization of the multidrug-resistance phenotype revealed the presence of a class 1 integron containing two genes, a dihydrofolate reductase (DHFR) (PF00186), which confers resistance to trimethoprim; and aminoglycoside adenyltransferase (AadA) (PF01909), which confers resistance to streptomycin and spectinomycin. A class 1 integron in K. pneumoniae containing these two genes was present in eight (18.18%) out of 44 different diabetic foot ulcer (DFU) patients. Hence, there is a need to develop therapeutics that inhibit growth of multidrug resistant K. pneumoniae in DFU patients and still achieve amputation control. Am attempt was made to create a 3D model and find a suitable inhibitor using an in silico study. Rational drug design/testing requires crystal structures for DHFR and AadA. However, the structures of DHFR and AadA from K. pneumoniae are not available. Modelling was performed using Swiss Model Server and Discovery Studio 3.1. The PDBSum server was used to check stereo chemical properties using Ramachandran plot analysis of modeled structures. Clindamycin was found to be suitable inhibitor of DHFR and AadA. A DockingServer based on Autodock & Mopac was used for docking calculations. The amino acid residues Ser(32), Ile(46), Glu(53), Gln(54), Phe(57), Thr(72), Met(76), Val(78), Leu(79), Ser(122), Tyr(128), Ile(151) in case of DHFR and Phe(34), Asp(60), Arg(63), Gln(64), Leu(68), Glu(87), Thr(89), Val(90) for AadA were found to be responsible for positioning clindamycin into the active site. The study identifies amino acid residues crucial to 'DHFR and AadA -drug' and 'DHFR and AadA -inhibitor' interactions that might be useful in the ongoing search for a versatile DHFR and AadA -inhibitor.

Soluble oligomerization provides a beneficial fitness effect on destabilizing mutations

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.

Characterization of an oxidoreductase from the arylamine N-acetyltransferase operon in Mycobacterium smegmatis

Mycobacterium tuberculosis, the most successful bacterial pathogen, causes tuberculosis, a disease that still causes more than 2 million deaths per year. Arylamine N-acetyltransferase is an enzyme that is conserved in most Mycobacterium spp. The nat gene belongs to an operon that is important for the intracellular survival of M. tuberculosis within macrophages. The nat operon in Mycobacterium smegmatis and other fast-growing mycobacterial species has a unique organization containing genes with uncharacterized function. Here, we describe the biochemical, biophysical and structural characterization of the MSMEG_0308 gene product (MS0308) of the M. smegmatis nat operon. While characterizing the function of MS0308, we validated the oxidoreductase property; however, we found that the enzyme was not utilizing dihydrofolate as its substrate, hence we first report that MS0308 is not a dihydrofolate reductase, as annotated in the genome. The structure of this oxidoreductase was solved at 2.0 Å in complex with the cofactor NADPH and has revealed the hydrophobic pocket where the endogenous substrate binds.

Protein dynamics and enzyme catalysis: insights from simulations

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Mechanisms of protein oligomerization, the critical role of insertions and deletions in maintaining different oligomeric states

The main principles of protein-protein recognition are elucidated by the studies of homooligomers which in turn mediate and regulate gene expression, activity of enzymes, ion channels, receptors, and cell-cell adhesion processes. Here we explore oligomeric states of homologous proteins in various organisms to better understand the functional roles and evolutionary mechanisms of homooligomerization. We observe a great diversity in mechanisms controlling oligomerization and focus in our study on insertions and deletions in homologous proteins and how they enable or disable complex formation. We show that insertions and deletions which differentiate monomers and dimers have a significant tendency to be located on the interaction interfaces and about a quarter of all proteins studied and forty percent of enzymes have regions which mediate or disrupt the formation of oligomers. We suggest that relatively small insertions or deletions may have a profound effect on complex stability and/or specificity. Indeed removal of complex enabling regions from protein structures in many cases resulted in the complete or partial loss of stability. Moreover, we find that insertions and deletions modulating oligomerization have a lower aggregation propensity and contain a larger fraction of polar, charged residues, glycine and proline compared to conventional interfaces and protein surface. Most likely, these regions may mediate specific interactions, prevent nonspecific dysfunctional aggregation and preclude undesired interactions between close paralogs therefore separating their functional pathways. Last, we show how the presence or absence of insertions and deletions on interfaces might be of practical value in annotating protein oligomeric states.

The temperature dependence of the kinetic isotope effects of dihydrofolate reductase from Thermotoga maritima is influenced by intersubunit interactions

Dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) is unique among structurally characterized chromosomal DHFRs in that it forms a stable homodimer. Dimerization is believed to play a key role in the high thermal stability of TmDHFR, which is reflected in a melting temperature in excess of 85 degrees C. The dimer interface of TmDHFR is composed of a hydrophobic core with charged residues around the periphery. In particular, Lys129 of each subunit forms three-membered salt bridges with Glu136 and Glu138 of the other subunit. To probe the role of these salt bridges in the dimerization and thermal stability of TmDHFR, we generated a series of variants (TmDHFR-K129E, TmDHFR-E136K, TmDHFR-E138K, and TmDHFR-E136K/E138K) in which these residues were exchanged for residues whose side chains bear the opposite charge. Our results indicate that these salt bridges are key for the high thermal stability of TmDHFR but are not a requirement for dimerization. Although the rate of dihydrofolate reduction by TmDHFR is not significantly affected by the loss of the K129-E136-E138 salt bridges, changes to the temperature dependence of the kinetic isotope effect on hydride transfer are observed. These changes are in agreement with the proposal that DHFR catalysis may be affected by changes to the conformational ensemble of the enzyme rather than only to the coupling of protein motions to the reaction coordinate.

Amino acid substitutions at protein-protein interfaces that modulate the oligomeric state

Despite similarities in their sequence and structure, there are a number of homologous proteins that adopt various oligomeric states. Comparisons of these homologous protein pairs, in terms of residue substitutions at the protein-protein interfaces, have provided fundamental characteristics that describe how proteins interact with each other. We have prepared a dataset composed of pairs of related proteins with different homo-oligomeric states. Using the protein complexes, the interface residues were identified, and using structural alignments, the shadow-interface residues have been defined as the surface residues that align with the interface residues. Subsequently, we investigated residue substitutions between the interfaces and the shadow interfaces. Based on the degree of the contributions to the interactions, the aligned sites of the interfaces and shadow interfaces were divided into primary and secondary sites; the primary sites are the focus of this work. The primary sites were further classified into two groups (i.e. exposed and buried) based on the degree to which the residue is buried within the shadow interfaces. Using these classifications, two simple mechanisms that mediate the oligomeric states were identified. In the primary-exposed sites, the residues on the shadow interfaces are replaced by more hydrophobic or aromatic residues, which are physicochemically favored at protein-protein interfaces. In the primary-buried sites, the residues on the shadow interfaces are replaced by larger residues that protrude into other proteins. These simple rules are satisfied in 23 out of 25 Structural Classification of Proteins (SCOP) families with a different-oligomeric-state pair, and thus represent a basic strategy for modulating protein associations and dissociations.

Solvent effects on catalysis by Escherichia coli dihydrofolate reductase

Hydride transfer catalyzed by dihydrofolate reductase (DHFR) has been described previously within an environmentally coupled model of hydrogen tunneling, where protein motions control binding of substrate and cofactor to generate a tunneling ready conformation and modulate the width of the activation barrier and hence the reaction rate. Changes to the composition of the reaction medium are known to perturb protein motions. We have measured kinetic parameters of the reaction catalyzed by DHFR from Escherichia coli in the presence of various cosolvents and cosolutes and show that the dielectric constant, but not the viscosity, of the reaction medium affects the rate of reaction. Neither the primary kinetic isotope effect on the reaction nor its temperature dependence were affected by changes to the bulk solvent properties. These results are in agreement with our previous report on the effect of solvent composition on catalysis by DHFR from the hyperthermophile Thermotoga maritima. However, the effect of solvent on the temperature dependence of the kinetic isotope effect on hydride transfer catalyzed by E. coli DHFR is difficult to explain within a model, in which long-range motions couple to the chemical step of the reaction, but may indicate the existence of a short-range promoting vibration or the presence of multiple nearly isoenergetic conformational substates of enzymes with similar but distinct catalytic properties.

Structural ordering of disordered ligand-binding loops of biotin protein ligase into active conformations as a consequence of dehydration

Mycobacterium tuberculosis (Mtb), a dreaded pathogen, has a unique cell envelope composed of high fatty acid content that plays a crucial role in its pathogenesis. Acetyl Coenzyme A Carboxylase (ACC), an important enzyme that catalyzes the first reaction of fatty acid biosynthesis, is biotinylated by biotin acetyl-CoA carboxylase ligase (BirA). The ligand-binding loops in all known apo BirAs to date are disordered and attain an ordered structure only after undergoing a conformational change upon ligand-binding. Here, we report that dehydration of Mtb-BirA crystals traps both the apo and active conformations in its asymmetric unit, and for the first time provides structural evidence of such transformation. Recombinant Mtb-BirA was crystallized at room temperature, and diffraction data was collected at 295 K as well as at 120 K. Transfer of crystals to paraffin and paratone-N oil (cryoprotectants) prior to flash-freezing induced lattice shrinkage and enhancement in the resolution of the X-ray diffraction data. Intriguingly, the crystal lattice rearrangement due to shrinkage in the dehydrated Mtb-BirA crystals ensued structural order of otherwise flexible ligand-binding loops L4 and L8 in apo BirA. In addition, crystal dehydration resulted in a shift of approximately 3.5 A in the flexible loop L6, a proline-rich loop unique to Mtb complex as well as around the L11 region. The shift in loop L11 in the C-terminal domain on dehydration emulates the action responsible for the complex formation with its protein ligand biotin carboxyl carrier protein (BCCP) domain of ACCA3. This is contrary to the involvement of loop L14 observed in Pyrococcus horikoshii BirA-BCCP complex. Another interesting feature that emerges from this dehydrated structure is that the two subunits A and B, though related by a noncrystallographic twofold symmetry, assemble into an asymmetric dimer representing the ligand-bound and ligand-free states of the protein, respectively. In-depth analyses of the sequence and the structure also provide answers to the reported lower affinities of Mtb-BirA toward ATP and biotin substrates. This dehydrated crystal structure not only provides key leads to the understanding of the structure/function relationships in the protein in the absence of any ligand-bound structure, but also demonstrates the merit of dehydration of crystals as an inimitable technique to have a glance at proteins in action.

A single intersubunit salt bridge affects oligomerization and catalytic activity in a bacterial quinone reductase

YhdA, a thermostable NADPH:FMN oxidoreductase from Bacillus subtilis, reduces quinones via a ping-pong bi-bi mechanism with a pronounced preference for NADPH. The enzyme occurs as a stable tetramer in solution. The two extended dimer surfaces are packed against each other by a 90 rotation of one dimer with respect to the other. This assembly is stabilized by the formation of four salt bridges between K109 and D137 of the neighbouring protomers. To investigate the importance of the ion pair contacts, the K109L and D137L single replacement variants, as well as the K109L/D137L and K109D/D137K double replacement variants, were generated, expressed, purified, crystallized and biochemically characterized. The K109L and D137L variants form dimers instead of tetramers, whereas the K109L/D137L and K109D/D137K variants appear to exist in a dimer-tetramer equilibrium in solution. The crystal structures of the K109L and D137L variants confirm the dimeric state, with the K109L/D137L and K109D/D137K variants adopting a tetrameric assembly. Interestingly, all protein variants show a drastically reduced quinone reductase activity in steady-state kinetics. Detailed analysis of the two half reactions revealed that the oxidative half reaction is not affected, whereas reduction of the bound FMN cofactor by NADPH is virtually abolished. Inspection of the crystal structures indicates that the side chain of K109 plays a dual role by forming a salt bridge to D137, as well as stabilizing a glycine-rich loop in the vicinity of the FMN cofactor. In all protein variants, this glycine-rich loop exhibits a much higher mobility, compared to the wild-type. This appears to be incompatible with NADPH binding and thus leads to abrogation of flavin reduction.

In pursuit of virtual lead optimization: pruning ensembles of receptor structures for increased efficiency and accuracy during docking

Representing receptors as ensembles of protein conformations during docking is a powerful method to approximate protein flexibility and increase the accuracy of the resulting ranked list of compounds. Unfortunately, docking compounds against a large number of ensemble members can increase computational cost and time investment. In this article, we present an efficient method to evaluate and select the most contributive ensemble members prior to docking for targets with a conserved core of residues that bind a ligand moiety. We observed that ensemble members that preserve the geometry of the active site core are most likely to place ligands in the active site with a conserved orientation, generally rank ligands correctly and increase interactions with the receptor. A relative distance approach is used to quantify the preservation of the three-dimensional interatomic distances of the conserved ligand-binding atoms and prune large ensembles quickly. In this study, we investigate dihydrofolate reductase as an example of a protein with a conserved core; however, this method for accurately selecting relevant ensemble members a priori can be applied to any system with a conserved ligand-binding core, including HIV-1 protease, kinases, and acetylcholinesterase. Representing a drug target as a pruned ensemble during in silico screening should increase the accuracy and efficiency of high-throughput analyses of lead analogs.

Probing coupled motions in enzymatic hydrogen tunnelling reactions

Much work has gone into understanding the physical basis of the enormous catalytic power of enzymes over the last 50 years or so. Nevertheless, the detailed mechanism used by Nature's catalysts to speed chemical transformations remains elusive. DHFR (dihydrofolate reductase) has served as a paradigm to study the relationship between the structure, function and dynamics of enzymatic transformations. A complex reaction cascade, which involves rearrangements and movements of loops and domains of the enzyme, is used to orientate cofactor and substrate in a reactive configuration from which hydride is transferred by quantum mechanical tunnelling. In the present paper, we review results from experiments that probe the influence of protein dynamics on the chemical step of the reaction catalysed by TmDHFR (DHFR from Thermotoga maritima). This enzyme appears to have evolved an optimal structure that can maintain a catalytically competent conformation under extreme conditions.

Solvent effects on environmentally coupled hydrogen tunnelling during catalysis by dihydrofolate reductase from Thermotoga maritima

Protein motions may be perturbed by altering the properties of the reaction medium. Here we show that dielectric constant, but not viscosity, affects the rate of the hydride-transfer reaction catalysed by dihydrofolate reductase from Thermotoga maritima (TmDHFR), in which quantum-mechanical tunnelling has previously been shown to be driven by protein motions. Neither dielectric constant nor viscosity directly alters the kinetic isotope effect of the reaction or the mechanism of coupling of protein motions to tunnelling. Glycerol and sucrose cause a significant increase in the rate of hydride transfer, but lead to a reduction in the magnitude of the kinetic isotope effect as well as an extension of the temperature range over which "passive" protein dynamics (rather than "active" gating motions) dominate the reaction. Our results are in agreement with the proposal that non-equilibrium dynamical processes (promoting motions) drive the hydride-transfer reaction in TmDHFR.

In pursuit of virtual lead optimization: the role of the receptor structure and ensembles in accurate docking

Accurate ranking during in silico lead optimization is critical to drive the generation of new ligands with higher affinity, yet it is especially difficult because of the subtle changes between analogs. In order to assess the role of the structure of the receptor in delivering accurate lead ranking results, we docked a set of forty related inhibitors to structures of one species of dihydrofolate reductase (DHFR) derived from crystallographic, NMR solution data, and homology models. In this study, the crystal structures yielded the superior results: the compounds were placed in the active site in the conserved orientation and the docking scores for 80% percent of the compounds clustered into the same bins as the measured affinity. Single receptor structures derived from NMR data or homology models did not serve as accurate docking receptors. To our knowledge, these are the first experiments that assess ranking of homologous lead compounds using a variety of receptor structures. We then extended the study to investigate whether ensembles, either computationally or experimentally derived, of all of the single starting structures aid, hinder or have no effect on the performance of the starting template. Impressively, when ensembles of receptor structures derived from NMR data or homology models were employed, docking accuracy improved to a level equal to that of the high resolution crystal structures. The same experiments using a second species of DHFR and set of ligands confirm the results. A comparison of the structures of the individual ensemble members to the starting structures shows that the effect of the ensembles can be ascribed to protein flexibility in addition to absorption of computational error.

Built-in loops allow versatility in domain-domain interactions: lessons from self-interacting domains

Compilations of domain-domain interactions based on solved structures suggest there are distinct domain pairs that are used repeatedly in different protein contexts to mediate protein-protein interactions. However, not all protein pairs with the corresponding domains that can potentially mediate interaction do interact, even when they are colocalized and coexpressed. It is conceivable that there are structural and sequence features, below the domain level, that play a role in determining the potential of domains to mediate protein-protein interactions. Here, we discover such features by comparing domains that, on the one hand, mediate homodimerization of proteins and, on the other, occur in different proteins that are documented as monomers. Intriguingly, this comparison uncovered surface loops that can be considered as determinants of the interactions. There are enabling loops, which mediate the domain interactions, and disabling loops that prevent the interactions. The presence of the enabling/disabling loops is consistent with the fulfillment/prevention of the interaction and is highly preserved in evolution. This suggests that, along with the preservation of structural elements that enable interaction, evolution maintains elements intended to prevent unwanted interactions. The enabling and disabling loops discovered in this study have implications in prediction of protein-protein interactions, by pointing to the protein regions that determine the interaction. Our results extend the hierarchy of attributes that collectively establish the modularity of domain-mediated protein-protein interactions.

On the relationship between thermal stability and catalytic power of enzymes

The possible relationship between the thermal stability and the catalytic power of enzymes is of great current interest. In particular, it has been suggested that thermophilic or hyperthermophilic (Tm) enzymes have lower catalytic power at a given temperature than the corresponding mesophilic (Ms) enzymes, because the thermophilic enzymes are less flexible (assuming that flexibility and catalysis are directly correlated). These suggestions presume that the reduced dynamics of the thermophilic enzymes is the reason for their reduced catalytic power. The present paper takes the specific case of dihydrofolate reductase (DHFR) and explores the validity of the above argument by simulation approaches. It is found that the Tm enzymes have restricted motions in the direction of the folding coordinate, but this is not relevant to the chemical process, since the motions along the reaction coordinate are perpendicular to the folding motions. Moreover, it is shown that the rate of the chemical reaction is determined by the activation barrier and the corresponding reorganization energy, rather than by dynamics or flexibility in the ground state. In fact, as far as flexibility is concerned, we conclude that the displacement along the reaction coordinate is larger in the Tm enzyme than in the Ms enzyme and that the general trend in enzyme catalysis is that the best catalyst involves less motion during the reaction than the less optimal catalyst. The relationship between thermal stability and catalysis appears to reflect the fact that to obtain small electrostatic reorganization energy it is necessary to invest some folding energy in the overall preorganization process. Thus, the optimized catalysts are less stable. This trend is clearly observed in the DHFR case.

Prediction of the structure of symmetrical protein assemblies

Biological supramolecular systems are commonly built up by the self-assembly of identical protein subunits to produce symmetrical oligomers with cyclical, icosahedral, or helical symmetry that play roles in processes ranging from allosteric control and molecular transport to motor action. The large size of these systems often makes them difficult to structurally characterize using experimental techniques. We have developed a computational protocol to predict the structure of symmetrical protein assemblies based on the structure of a single subunit. The method carries out simultaneous optimization of backbone, side chain, and rigid-body degrees of freedom, while restricting the search space to symmetrical conformations. Using this protocol, we can reconstruct, starting from the structure of a single subunit, the structure of cyclic oligomers and the icosahedral virus capsid of satellite panicum virus using a rigid backbone approximation. We predict the oligomeric state of EscJ from the type III secretion system both in its proposed cyclical and crystallized helical form. Finally, we show that the method can recapitulate the structure of an amyloid-like fibril formed by the peptide NNQQNY from the yeast prion protein Sup35 starting from the amino acid sequence alone and searching the complete space of backbone, side chain, and rigid-body degrees of freedom.

Implication for buried polar contacts and ion pairs in hyperthermostable enzymes

Understanding the structural basis of thermostability and thermoactivity, and their interdependence, is central to the successful future exploitation of extremophilic enzymes in biotechnology. However, the structural basis of thermostability is still not fully characterized. Ionizable residues play essential roles in proteins, modulating protein stability, folding and function. The dominant roles of the buried polar contacts and ion pairs have been reviewed by distinguishing between the inside polar contacts and the total intramolecular polar contacts, and by evaluating their contribution as molecular determinants for protein stability using various protein structures from hyperthermophiles, thermophiles and mesophilic organisms. The analysis revealed that the remarkably increased number of internal polar contacts in a monomeric structure probably play a central role in enhancing the melting temperature value up to 120 degrees C for hyperthermophilic enzymes from the genus Pyrococcus. These results provide a promising contribution for improving the thermostability of enzymes by modulating buried polar contacts and ion pairs.

Nanoscale dewetting transition in protein complex folding

In a previous study, a surprising drying transition was observed to take place inside the nanoscale hydrophobic channel in the tetramer of the protein melittin. The goal of this paper is to determine if there are other protein complexes capable of displaying a dewetting transition during their final stage of folding. We searched the entire protein data bank (PDB) for all possible candidates, including protein tetramers, dimers, and two-domain proteins, and then performed the molecular dynamics (MD) simulations on the top candidates identified by a simple hydrophobic scoring function based on aligned hydrophobic surface areas. Our large scale MD simulations found several more proteins, including three tetramers, six dimers, and two two-domain proteins, which display a nanoscale dewetting transition in their final stage of folding. Even though the scoring function alone is not sufficient (i.e., a high score is necessary but not sufficient) in identifying the dewetting candidates, it does provide useful insights into the features of complex interfaces needed for dewetting. All top candidates have two features in common: (1) large aligned (matched) hydrophobic areas between two corresponding surfaces, and (2) large connected hydrophobic areas on the same surface. We have also studied the effect on dewetting of different water models and different treatments of the long-range electrostatic interactions (cutoff vs PME), and found the dewetting phenomena is fairly robust. This work presents a few proteins other than melittin tetramer for further experimental studies of the role of dewetting in the end stages of protein folding.