Role of electrostatics in cold adaptation: A comparative study of eury- and stenopsychrophilic triose phosphate isomerase

Psychrophilic (cold-active) organisms have developed enzymes that facilitate sufficient metabolic activity at low temperatures to sustain life. This occurs through molecular adaptations that tend to increase protein flexibility at the expense of stability. However, psychrophiles also vary in their growth conditions. Eurypsychrophiles thrive over a wide temperature range and often prefer temperatures above 20 °C, while stenopsychrophiles grow optimally below 15 °C and are more narrowly adapted to cold temperatures. To elucidate differences between these two classes of enzymes, we here compare the stability and unfolding kinetics of two orthologues of the basal household enzyme triose phosphate isomerase, one from the stenopsychrophilic Antarctic permafrost bacterium Rhodonellum psychrophilum (sTPI) and the other from the eurypsychrophilic Greenland ikaite column bacterium Rhodococcus sp. JG-3 (eTPI). Remarkably, sTPI proved significantly more thermostable and resistant to chemical denaturation than its eurypsychrophilic counterpart, eTPI, in the absence of ionic components in solution, whereas inclusion of electrostatic screening agents in the form of sodium chloride or the charged denaturant guanidinium chloride largely cancelled out this difference. Thus, electrostatics play a prominent role in stabilizing the stenopsychrophilic sTPI, and a mandatory low-temperature growth environment does not preclude the development of considerable thermotolerance for individual enzymes. We were able to increase the thermostability of sTPI using an evolutionary machine learning model, which transferred several sTPI residues into the eTPI active site. While the stabilizing effect was modest, the combination of individual mutations was additive, underscoring the potential of combining multiple beneficial mutations to achieve enhanced enzyme properties.

Deciphering the folding code of collagens

Collagen proteins contain a characteristic structural motif called a triple helix. During the self-assembly of this motif, three polypeptides form a folding nucleus at the C-termini and then propagate towards the N-termini like a zip-chain. While polypeptides from human collagens contain up to a 1000 amino acids, those found in bacteria can contain up to 6000 amino acids. Additionally, the collagen polypeptides are also frequently interrupted by non-helical sequences that disrupt folding and reduce stability. Given the length of polypeptides and the disruptive interruptions, compensating mechanisms that stabilize against local unfolding during propagation and offset the entropic cost of folding are not fully understood. Here, we show that the information for the correct folding of collagen triple helices is encoded in their sequence as interchain electrostatic interactions, which likely act as molecular clamps that prevent local unfolding. In the case of humans, disrupting these electrostatic interactions is associated with severe to lethal diseases.

Stable monomers in the ancestral sequence reconstruction of the last opisthokont common ancestor of dimeric triosephosphate isomerase

Function and structure are strongly coupled in obligated oligomers such as Triosephosphate isomerase (TIM). In animals and fungi, TIM monomers are inactive and unstable. Previously, we used ancestral sequence reconstruction to study TIM evolution and found that before these lineages diverged, the last opisthokonta common ancestor of TIM (LOCATIM) was an obligated oligomer that resembles those of extant TIMs. Notably, calorimetric evidence indicated that ancestral TIM monomers are more structured than extant ones. To further increase confidence about the function, structure, and stability of the LOCATIM, in this work, we applied two different inference methodologies and the worst plausible case scenario for both of them, to infer four sequences of this ancestor and test the robustness of their physicochemical properties. The extensive biophysical characterization of the four reconstructed sequences of LOCATIM showed very similar hydrodynamic and spectroscopic properties, as well as ligand-binding energetics and catalytic parameters. Their 3D structures were also conserved. Although differences were observed in melting temperature, all LOCATIMs showed reversible urea-induced unfolding transitions, and for those that reached equilibrium, high conformational stability was estimated (ΔGTot = 40.6-46.2 kcal/mol). The stability of the inactive monomeric intermediates was also high (ΔGunf = 12.6-18.4 kcal/mol), resembling some protozoan TIMs rather than the unstable monomer observed in extant opisthokonts. A comparative analysis of the 3D structure of ancestral and extant TIMs shows a correlation between the higher stability of the ancestral monomers with the presence of several hydrogen bonds located in the "bottom" part of the barrel.

Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

DomainMapper: Accurate domain structure annotation including those with non-contiguous topologies

Automated domain annotation is an important tool for structural informatics. These pipelines typically involve searching query sequences against hidden Markov model (HMM) profiles, yielding matches to profiles for various domains. However, domain annotation can be ambiguous or inaccurate when proteins contain domains with non-contiguous residue ranges, and especially when insertional domains are hosted within them. Here, we present DomainMapper, an algorithm that accurately assigns a unique domain structure annotation to a query sequence, including those with complex topologies. We validate our domain assignments using the AlphaFold database and confirm that non-contiguity is pervasive (10.74% of all domains in yeast and 4.52% in human). Using this resource, we find that certain folds have strong propensities to be non-contiguous or insertional across the Tree of Life. DomainMapper is freely available and can be ran as a single command-line function.

Synergistic effect of compounds directed to triosephosphate isomerase, a combination to develop drug against trichomoniasis

The development of new drugs is continuous in the world; currently, saving resources (both economic ones and time) and preventing secondary effects have become a necessity for drug developers. Trichomoniasis is the most common nonviral sexually transmitted infection affecting more than 270 million people around the world. In our research group, we focussed on developing a selective and more effective drug against Trichomonas vaginalis, and we previously reported on a compound, called A4, which had a trichomonacidal effect. Later, we determined another compound, called D4, which also had a trichomonacidal effect together with favorable toxicity results. Both A4 and D4 are directed at the enzyme triosephosphate isomerase. Thus, we made combinations between the two compounds, in which we determined a synergistic effect against T. vaginalis, determining the IC₅₀ and the toxicity of the best relationship to obtain the trichomonacidal effect. With these results, we can propose a combination of compounds that represents a promising alternative for the development of a new therapeutic strategy against trichomoniasis.

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

•

The TIM barrel is a versatile fold to understand structure-stability relationships.

•

A collection of de novo TIM barrels with improved hydrophobic cores was designed.

•

DeNovoTIMs are reversible in chemical and thermal unfolding, which is uncommon in TIM barrels.

•

Epistatic effects play a central role in DeNovoTIMs stabilization.

•

DeNovoTIMs navigate a previously uncharted region of the stability landscape.

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.

Sortase mutants with improved protein thermostability and enzymatic activity obtained by consensus design

Staphylococcus aureus sortase A (SaSrtA) is an enzyme that anchors proteins to the cell surface of Gram-positive bacteria. During the transpeptidation reaction performed by SaSrtA, proteins containing an N-terminal glycine can be covalently linked to another protein with a C-terminal LPXTG motif (X being any amino acid). Since the sortase reaction can be performed in vitro as well, it has found many applications in biotechnology. Although sortase-mediated ligation has many advantages, SaSrtA is limited by its low enzymatic activity and dependence on Ca2+. In our study, we evaluated the thermodynamic stability of the SaSrtA wild type and found the enzyme to be stable. We applied consensus analysis to further improve the enzyme's stability while at the same time enhancing the enzyme's activity. As a result, we found thermodynamically improved, more active and Ca2+-independent mutants. We envision that these new variants can be applied in conjugation reactions in low Ca2+ environments.

Native aggregation is a common feature among triosephosphate isomerases of different species

Triosephosphate isomerase (TIM) is an enzyme of the glycolysis pathway which exists in almost all types of cells. Its structure is the prototype of a motif called TIM-barrel or (α/β)₈ barrel, which is the most common fold of all known enzyme structures. The simplest form in which TIM is catalytically active is a homodimer, in many species of bacteria and eukaryotes, or a homotetramer in some archaea. Here we show that the purified homodimeric TIMs from nine different species of eukaryotes and one of an extremophile bacterium spontaneously form higher order aggregates that can range from 3 to 21 dimers per macromolecular complex. We analysed these aggregates with clear native electrophoresis with normal and inverse polarity, blue native polyacrylamide gel electrophoresis, liquid chromatography, dynamic light scattering, thermal shift assay and transmission electron and fluorescence microscopies, we also performed bioinformatic analysis of the sequences of all enzymes to identify and predict regions that are prone to aggregation. Additionally, the capacity of TIM from Trypanosoma brucei to form fibrillar aggregates was characterized. Our results indicate that all the TIMs we studied are capable of forming oligomers of different sizes. This is significant because aggregation of TIM may be important in some of its non-catalytic moonlighting functions, like being a potent food allergen, or in its role associated with Alzheimer’s disease.

Structural basis for the modulation of plant cytosolic triosephosphate isomerase activity by mimicry of redox-based modifications

Reactive oxidative species (ROS) and S-glutathionylation modulate the activity of plant cytosolic triosephosphate isomerases (cTPI). Arabidopsis thaliana cTPI (AtcTPI) is subject of redox regulation at two reactive cysteines that function as thiol switches. Here we investigate the role of these residues, AtcTPI-Cys13 and At-Cys218, by substituting them with aspartic acid that mimics the irreversible oxidation of cysteine to sulfinic acid and with amino acids that mimic thiol conjugation. Crystallographic studies show that mimicking AtcTPI-Cys13 oxidation promotes the formation of inactive monomers by reposition residue Phe75 of the neighboring subunit, into a conformation that destabilizes the dimer interface. Mutations in residue AtcTPI-Cys218 to Asp, Lys, or Tyr generate TPI variants with a decreased enzymatic activity by creating structural modifications in two loops (loop 7 and loop 6) whose integrity is necessary to assemble the active site. In contrast with mutations in residue AtcTPI-Cys13, mutations in AtcTPI-Cys218 do not alter the dimeric nature of AtcTPI. Therefore, modifications of residues AtcTPI-Cys13 and AtcTPI-Cys218 modulate AtcTPI activity by inducing the formation of inactive monomers and by altering the active site of the dimeric enzyme, respectively. The identity of residue AtcTPI-Cys218 is conserved in the majority of plant cytosolic TPIs, this conservation and its solvent-exposed localization make it the most probable target for TPI regulation upon oxidative damage by reactive oxygen species. Our data reveal the structural mechanisms by which S-glutathionylation protects AtcTPI from irreversible chemical modifications and re-routes carbon metabolism to the pentose phosphate pathway to decrease oxidative stress.

Active-Site Glu165 Activation in Triosephosphate Isomerase and Its Deprotonation Kinetics

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) via an enediol(ate) intermediate. The active-site residue Glu165 serves as the catalytic base during catalysis. It abstracts a proton from C1 carbon of DHAP to form the reaction intermediate and donates a proton to C2 carbon of the intermediate to form product GAP. Our difference Fourier transform infrared spectroscopy studies on the yeast TIM (YeTIM)/phosphate complex revealed a C═O stretch band at 1706 cm^-1 from the protonated Glu165 carboxyl group at pH 7.5, indicating that the p K_a of the catalytic base is increased by >3.0 pH units upon phosphate binding, and that the Glu165 carboxyl environment in the complex is still hydrophilic in spite of the increased p K_a. Hence, the results show that the binding of the phosphodianion group is part of the activation mechanism which involves the p K_a elevation of the catalytic base Glu165. The deprotonation kinetics of Glu165 in the μs to ms time range were determined via infrared (IR) T-jump studies on the YeTIM/phosphate and ("heavy enzyme") [U-¹³C,-¹⁵N]YeTIM/phosphate complexes. The slower deprotonation kinetics in the ms time scale is due to phosphate dissociation modulated by the loop motion, which slows down by enzyme mass increase to show a normal heavy enzyme kinetic isotope effect (KIE) ∼1.2 (i.e., slower rate in the heavy enzyme). The faster deprotonation kinetics in the tens of μs time scale is assigned to temperature-induced p K_a decrease, while phosphate is still bound, and it shows an inverse heavy enzyme KIE ∼0.89 (faster rate in the heavy enzyme). The IR static and T-jump spectroscopy provides atomic-level resolution of the catalytic mechanism because of its ability to directly observe the bond breaking/forming process.

Structure and conformational stability of the triosephosphate isomerase from Zea mays. Comparison with the chemical unfolding pathways of other eukaryotic TIMs

We studied the structure, function and thermodynamic properties for the unfolding of the Triosephosphate isomerase (TIM) from Zea mays (ZmTIM). ZmTIM shows a catalytic efficiency close to the diffusion limit. Native ZmTIM is a dimer that dissociates upon dilution into inactive and unfolded monomers. Its thermal unfolding is irreversible with a T_m of 61.6 ± 1.4 °C and an activation energy of 383.4 ± 11.5 kJ mol^-1. The urea-induced unfolding of ZmTIM is reversible. Transitions followed by catalytic activity and spectroscopic properties are monophasic and superimposable, indicating that ZmTIM unfolds/refolds in a two-state behavior with an unfolding ΔG°_(H20) = 99.8 ± 5.3 kJ mol^-1. This contrasts with most other studied TIMs, where folding intermediates are common. The three-dimensional structure of ZmTIM was solved at 1.8 Å. A structural comparison with other eukaryotic TIMs shows a similar number of intramolecular and intermolecular interactions. Interestingly the number of interfacial water molecules found in ZmTIM is lower than those observed in most TIMs that show folding intermediates. Although with the available data, there is no clear correlation between structural properties and the number of equilibrium intermediates in the unfolding of TIM, the identification of such structural properties should increase our understanding of folding mechanisms.

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.

De novo design of a four-fold symmetric TIM-barrel protein with atomic-level accuracy

Despite efforts for over 25 years, de novo protein design has not succeeded in achieving the TIM-barrel fold. Here we describe the computational design of four-fold symmetrical (β/α)8 barrels guided by geometrical and chemical principles. Experimental characterization of 33 designs revealed the importance of side chain-backbone hydrogen bonds for defining the strand register between repeat units. The X-ray crystal structure of a designed thermostable 184-residue protein is nearly identical to that of the designed TIM-barrel model. PSI-BLAST searches do not identify sequence similarities to known TIM-barrel proteins, and sensitive profile-profile searches indicate that the design sequence is distant from other naturally occurring TIM-barrel superfamilies, suggesting that Nature has sampled only a subset of the sequence space available to the TIM-barrel fold. The ability to design TIM barrels de novo opens new possibilities for custom-made enzymes.