Solvent environments significantly affect the enzymatic function of Escherichia coli dihydrofolate reductase: comparison of wild-type protein and active-site mutant D27E

Pressure Adaptations in Deep-Sea Moritella Dihydrofolate Reductases: Compressibility versus Stability

Simple Summary

Deep-sea organisms must have proteins that function under high hydrostatic pressure to survive. Adaptations used in proteins from “pressure-loving” piezophiles may include greater compressibility or greater stability against pressure-induced destabilization. However, while greater compressibility can be accomplished by greater void volume, larger cavities in a protein have been associated with greater destabilization and even unfolding as pressure is increased. Here, computer simulations of dihydrofolate reductase from a moderate piezophile and a hyperpiezophile were performed to understand the balance between adaptations for greater compressibility and those against pressure destabilization and unfolding. The results indicate that while compressibility appears to be important for deep-sea microbes, adaptation for the greatest depths may be to prevent water penetration into the interior.

Abstract

Proteins from “pressure-loving” piezophiles appear to adapt by greater compressibility via larger total cavity volume. However, larger cavities in proteins have been associated with lower unfolding pressures. Here, dihydrofolate reductase (DHFR) from a moderate piezophile Moritella profunda (Mp) isolated at ~2.9 km in depth and from a hyperpiezophile Moritella yayanosii (My) isolated at ~11 km in depth were compared using molecular dynamics simulations. Although previous simulations indicate that MpDHFR is more compressible than a mesophile DHFR, here the average properties and a quasiharmonic analysis indicate that MpDHFR and MyDHFR have similar compressibilities. A cavity analysis also indicates that the three unique mutations in MyDHFR are near cavities, although the cavities are generally similar in size in both. However, while a cleft overlaps an internal cavity, thus forming a pathway from the surface to the interior in MpDHFR, the unique residue Tyr103 found in MyDHFR forms a hydrogen bond with Leu78, and the sidechain separates the cleft from the cavity. Thus, while Moritella DHFR may generally be well suited to high-pressure environments because of their greater compressibility, adaptation for greater depths may be to prevent water entry into the interior cavities.

Adaptations for Pressure and Temperature in Dihydrofolate Reductases

Enzymes from extremophilic microbes that live in extreme conditions are generally adapted so that they function under those conditions, although adaptations for extreme temperatures and pressures can be difficult to unravel. Previous studies have shown mutation of Asp27 in Escherichia coli dihydrofolate reductase (DHFR) to Glu27 in Moritella profunda (Mp). DHFR enhances activity at higher pressures, although this may be an adaptation for cold. Interestingly, MpDHFR unfolds at ~70 MPa, while Moritella yayanosii (My) was isolated at depths corresponding to ~110 MPa, indicating that MyDHFR might be adapted for higher pressures. Here, these adaptations are examined using molecular dynamics simulations of DHFR from different microbes in the context of not only experimental studies of activity and stability of the protein but also the evolutionary history of the microbe. Results suggest Tyr103 of MyDHFR may be an adaptation for high pressure since Cys103 in helix F of MpDHFR forms an intra-helix hydrogen bond with Ile99 while Tyr103 in helix F of MyDHFR forms a hydrogen bond with Leu78 in helix E. This suggests the hydrogen bond between helices F and E in MyDHFR might prevent distortion at higher pressures.

How adding a single methylene to dihydrofolate reductase can change its conformational dynamics

Studies of the effects of pressure on proteins from piezophilic (pressure-loving) microbes compared with homologous proteins from mesophilic microbes have been relatively rare. Interestingly, such studies of dihydrofolate reductase show that a single-site mutation from an aspartic acid to a glutamic acid can reverse the pressure-dependent monotonic decrease in activity to that in a monotonic pressure-dependent activation. This residue is near the active site but is not thought to directly participate in the catalytic mechanism. Here, the ways that addition of one carbon to the entire protein could lead to such a profound difference in pressure effects are explored using molecular dynamics simulations. The results indicate that the glutamate changes the coupling between a helix and the β-sheet due to the extra flexibility of the side chain, which further changes correlated motions of other regions of the protein.

Adaptations for Pressure and Temperature Effects on Loop Motion in Escherichia coli and Moritella profunda Dihydrofolate Reductase

Determining how enzymes in piezophilic microbes function at high pressure can give insights into how life adapts to living at high pressure. Here, the effects of pressure and temperature on loop motions are compared Escherichia coli (Ec) and Moritella profunda (Mp) dihydrofolate reductase (DHFR) via molecular dynamics simulations at combinations of the growth temperature and pressure of the two organisms. Analysis indicates that a flexible CD loop in MpDHFR is an adaptation for cold because it makes the adenosine binding subdomain more flexible. Also, analysis indicates that the Thr113-Glu27 hydrogen bond in MpDHFR is an adaptation for high pressure because it provides flexibility within the loop subdomain compared to the very strong Thr113-Asp27 hydrogen bond in EcDHFR, and affects the correlation of the Met20 and GH loops. In addition, the results suggest that temperature might affect external loops more strongly while pressure might affect motion between elements within the protein.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions

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

Enzymes from piezophiles

The discovery of microbial communities in extreme conditions that would seem hostile to life leads to the question of how the molecules making up these microbes can maintain their structure and function. While microbes that live under extremes of temperature have been heavily studied, those that live under extremes of pressure, or "piezophiles", are now increasingly being studied because of advances in sample collection and high-pressure cells for biochemical and biophysical measurements. Here, adaptations of enzymes in piezophiles against the effects of pressure are discussed in light of recent experimental and computational studies. However, while concepts from studies of enzymes from temperature extremophiles can provide frameworks for understanding adaptations by piezophile enzymes, the effects of temperature and pressure on proteins differ in significant ways. Thus, the state of the knowledge of adaptation in piezophile enzymes is still in its infancy and many more experiments and computational studies on different enzymes from a variety of piezophiles are needed.

Quasiharmonic Analysis of the Energy Landscapes of Dihydrofolate Reductase from Piezophiles and Mesophiles

A quasiharmonic analysis (QHA) method is used to compare the potential energy landscapes of dihydrofolate reductase (DHFR) from a piezophile (pressure-loving organism), Moritella profunda (Mp), and a mesophile, Escherichia coli (Ec). The QHA method considers atomic fluctuations of the protein as motions of an atom in a local effective potential created by neighboring atoms and quantitates it in terms of effective force constants, isothermal compressibilities, and thermal expansivities. The analysis indicates that the underlying potential energy surface of MpDHFR is inherently softer than that of EcDHFR. In addition, on picosecond time scales, the energy surfaces become more similar under the growth conditions of Mp and Ec. On these time scales, DHFR behaves as expected; namely, increasing temperature makes the effective energy minimum less steep because thermal fluctuations increase the available volume, whereas increasing pressure steepens it because compression reduces the available volume. Our longer simulations show that, on nanosecond time scales, increasing temperature has a similar effect as on picosecond time scales because thermal fluctuations increase the volume even more on a longer time scale. However, these simulations also indicate that, on nanosecond time scales, pressure makes the local potential less steep, contrary to picosecond time scales. Further examination of the QHA indicates the nanosecond pressure response may originate at picosecond time scales at the exterior of the protein, which suggests that protein-water interactions may be involved. The results may lead to understanding adaptations in enzymes made by piezophiles that enable them to function at higher pressures.

Similar structural stabilities of 3-isopropylmalate dehydrogenases from the obligatory piezophilic bacterium Shewanella benthica strain DB21MT-2 and its atmospheric congener S. oneidensis strain MR-1

We previously found that the enzymatic activity of 3-isopropylmalate dehydrogenase from the obligatory piezophilic bacterium Shewanella benthica strain DB21MT-2 (SbIPMDH) was pressure-tolerant up to 100 MPa, but that from its atmospheric congener S. oneidensis strain MR-1 (SoIPMDH) was pressure-sensitive. Such characteristics were determined by only one amino acid residue at position 266, serine (SoIPMDH) or alanine (SbIPMDH) [Y. Hamajima et al. Extremophiles 20: 177, 2016]. In this study, we investigated the structural stability of these enzymes. At pH 7.6, SoIPMDH was slightly more stable against hydrostatic pressure than SbIPMDH, contrary to the physiological pressures of their normal environments. Pressure unfolding of these IPMDHs followed a two-state unfolding model between a native dimer and two unfolded monomers, and the dimer structure was pressure-tolerant up to 200 MPa, employing a midpoint pressure of 245.3 ± 0.1 MPa and a volume change of -225 ± 24 mL mol^-1 for the most unstable mutant, SbIPMDH A266S. Thus, their pressure-dependent activity did not originate from structural perturbations such as unfolding or dimer dissociation. Conversely, urea-induced unfolding of these IPMDHs followed a three-state unfolding model, including a dimer intermediate. Interestingly, the first transition was strongly pH-dependent but pressure-independent; however, the second transition showed the opposite pattern. Obtained volume changes due to urea-induced unfolding were almost equal for both IPMDHs, approximately +10 and -30 mL mol^-1 for intermediate formation and dimer dissociation, respectively. These results indicated that both IPMDHs have similar structural stability, and a pressure-adaptation mechanism was provided for only the enzymatic activity of SbIPMDH.

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

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

Pressure adaptation of 3-isopropylmalate dehydrogenase from an extremely piezophilic bacterium is attributed to a single amino acid substitution

3-Isopropylmalate dehydrogenase (IPMDH) from the extreme piezophile Shewanella benthica (SbIPMDH) is more pressure-tolerant than that from the atmospheric pressure-adapted Shewanella oneidensis (SoIPMDH). To understand the molecular mechanisms of this pressure tolerance, we analyzed mutated enzymes. The results indicate that only a single mutation at position 266, corresponding to Ala (SbIPMDH) and Ser (SoIPMDH), essentially affects activity under higher-pressure conditions. Structural analyses of SoIPMDH suggests that penetration of three water molecules into the cleft around Ser266 under high-pressure conditions could reduce the activity of the wild-type enzyme; however, no water molecule is observed in the Ala266 mutant.

Surface water retardation around single-chain polymeric nanoparticles: critical for catalytic function?

A library of water-soluble dynamic single-chain polymeric nanoparticles (SCPN) was prepared using a controlled radical polymerisation technique followed by the introduction of functional groups, including probes at targeted positions. The combined tools of electron paramagnetic resonance (EPR) and Overhauser dynamic nuclear polarization (ODNP) reveal that these SCPNs have structural and surface hydration properties resembling that of enzymes.

Vacuum-Ultraviolet Circular Dichroism Spectra of Escherichia coli Dihydrofolate Reductase and Its Mutants: Contributions of Phenylalanine and Tyrosine Side Chains and Exciton Coupling of Two Tryptophan Side Chains

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.

Cofactor-Mediated Conformational Dynamics Promote Product Release From Escherichia coli Dihydrofolate Reductase via an Allosteric Pathway

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.

Thermal stabilization of dihydrofolate reductase using monte carlo unfolding simulations and its functional consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.

All-atom molecular simulations have provided valuable insight into the workings of molecular machines and the folding and unfolding of proteins. However, commonly employed molecular dynamics simulations suffer from a limitation in accessible time scale, making it difficult to model large-scale unfolding events in a realistic amount of simulation time without employing unrealistically high temperatures. Here, we describe a rapid all-atom Monte Carlo simulation approach to simulate unfolding of the essential bacterial enzyme Dihydrofolate Reductase (DHFR) and all possible single point-mutants. We use these simulations to predict which mutants will be more thermodynamically stable (i.e., reside more often in the native folded state vs. the unfolded state) than the wild-type protein, and we confirm our predictions experimentally, creating several highly stable and catalytically active mutants. Thermally stable active engineered proteins can be used as a starting point in directed evolution experiments to evolve new functions on the background of this additional “reservoir of stability.” The stabilized enzyme may be able to accumulate a greater number of destabilizing yet functionally important mutations before unfolding, protease digestion, and aggregation abolish its activity.

Structure, Activity and Stereoselectivity of NADPH-Dependent Oxidoreductases Catalysing the S-Selective Reduction of the Imine Substrate 2-Methylpyrroline

Oxidoreductases from Streptomyces sp. GF3546 [3546-IRED], Bacillus cereus BAG3X2 (BcIRED) and Nocardiopsis halophila (NhIRED) each reduce prochiral 2-methylpyrroline (2MPN) to (S)-2-methylpyrrolidine with >95 % ee and also a number of other imine substrates with good selectivity. Structures of BcIRED and NhIRED have helped to identify conserved active site residues within this subgroup of imine reductases that have S selectivity towards 2MPN, including a tyrosine residue that has a possible role in catalysis and superimposes with an aspartate in related enzymes that display R selectivity towards the same substrate. Mutation of this tyrosine residue-Tyr169-in 3546-IRED to Phe resulted in a mutant of negligible activity. The data together provide structural evidence for the location and significance of the Tyr residue in this group of imine reductases, and permit a comparison of the active sites of enzymes that reduce 2MPN with either R or S selectivity.

Effects of salt on the structure, stability, and function of a halophilic dihydrofolate reductase from a hyperhalophilic archaeon, Haloarcula japonica strain TR-1

The effects of salt on the structure, stability, and enzymatic function of a novel dihydrofolate reductase (HjDHFR P1) from a hyperhalophilic archaeon, Haloarcula japonica strain TR-1 living in a Japanese saltern, were studied using ultraviolet absorption, circular dichroism (CD), and fluorescence spectroscopy. HjDHFR P1 had a partial structure at pH 8.0 in the absence of NaCl, and the addition of NaCl (0-500 mM concentration) induced significant structural formation to HjDHFR P1. The addition of NADPH, which is a coenzyme for its catalytic reaction, and lowering the pH from 8 to 6 also induced the same CD change, indicating the formation of the NADPH-binding site in HjDHFR P1. The NaCl dependence of thermal and urea-induced unfolding measurements suggested that protein stability increased depending on NaCl concentration regardless of structural formation, and HjDHFR P1 achieved the same stability as Escherichia coli DHFR at 750 mM NaCl. Halophilic characteristics were also observed for enzymatic function, although its structure had already formed under the conditions that enzymatic activity was measured at due to the presence of NADPH. These results suggest that the halophilic mechanism on structural stability and function was caused by factors other than structural formation, which are suggested to be the contributions of preferential interactions between the protein and salt ions and the specific binding of salt ions.

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.