Lactate dehydrogenase from the hyperthermophilic bacterium thermotoga maritima: the crystal structure at 2.1 A resolution reveals strategies for intrinsic protein stabilization

Biochemical, structural and dynamical characterizations of the lactate dehydrogenase from Selenomonas ruminantium provide information about an intermediate evolutionary step prior to complete allosteric regulation acquisition in the super family of lactate and malate dehydrogenases

In this work, we investigated the lactate dehydrogenase (LDH) from Selenomonas ruminantium (S. rum), an enzyme that differs at key amino acid positions from canonical allosteric LDHs. The wild type (Wt) of this enzyme recognises pyuvate as all LDHs. However, introducing a single point mutation in the active site loop (I85R) allows S. Rum LDH to recognize the oxaloacetate substrate as a typical malate dehydrogenase (MalDH), whilst maintaining homotropic activation as an LDH. We report the tertiary structure of the Wt and I85RLDH mutant. The Wt S. rum enzyme structure binds NADH and malonate, whilst also resembling the typical compact R-active state of canonical LDHs. The structure of the mutant with I85R was solved in the Apo State (without ligand), and shows no large conformational reorganization such as that observed with canonical allosteric LDHs in Apo state. This is due to a local structural feature typical of S. rum LDH that prevents large-scale conformational reorganization. The S. rum LDH was also studied using Molecular Dynamics simulations, probing specific local deformations of the active site that allow the S. rum LDH to sample the T-inactive state. We propose that, with respect to the LDH/MalDH superfamily, the S. rum enzyme possesses a specificstructural and dynamical way to ensure homotropic activation.

Deciphering Evolutionary Trajectories of Lactate Dehydrogenases Provides New Insights into Allostery

Lactate dehydrogenase (LDH, EC.1.1.127) is an important enzyme engaged in the anaerobic metabolism of cells, catalyzing the conversion of pyruvate to lactate and NADH to NAD+. LDH is a relevant enzyme to investigate structure-function relationships. The present work provides the missing link in our understanding of the evolution of LDHs. This allows to explain (i) the various evolutionary origins of LDHs in eukaryotic cells and their further diversification and (ii) subtle phenotypic modifications with respect to their regulation capacity. We identified a group of cyanobacterial LDHs displaying eukaryotic-like LDH sequence features. The biochemical and structural characterization of Cyanobacterium aponinum LDH, taken as representative, unexpectedly revealed that it displays homotropic and heterotropic activation, typical of an allosteric enzyme, whereas it harbors a long N-terminal extension, a structural feature considered responsible for the lack of allosteric capacity in eukaryotic LDHs. Its crystallographic structure was solved in 2 different configurations typical of the R-active and T-inactive states encountered in allosteric LDHs. Structural comparisons coupled with our evolutionary analyses helped to identify 2 amino acid positions that could have had a major role in the attenuation and extinction of the allosteric activation in eukaryotic LDHs rather than the presence of the N-terminal extension. We tested this hypothesis by site-directed mutagenesis. The resulting C. aponinum LDH mutants displayed reduced allosteric capacity mimicking those encountered in plants and human LDHs. This study provides a new evolutionary scenario of LDHs that unifies descriptions of regulatory properties with structural and mutational patterns of these important enzymes.

Protein Conformational Space at the Edge of Allostery: Turning a Non-allosteric Malate Dehydrogenase into an "Allosterized" Enzyme using Evolution Guided Punctual Mutations

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.

Unravelling the Adaptation Mechanisms to High Pressure in Proteins

Life is thought to have appeared in the depth of the sea under high hydrostatic pressure. Nowadays, it is known that the deep biosphere hosts a myriad of life forms thriving under high-pressure conditions. However, the evolutionary mechanisms leading to their adaptation are still not known. Here, we show the molecular bases of these mechanisms through a joint structural and dynamical study of two orthologous proteins. We observed that pressure adaptation involves the decoupling of protein–water dynamics and the elimination of cavities in the protein core. This is achieved by rearranging the charged residues on the protein surface and using bulkier hydrophobic residues in the core. These findings will be the starting point in the search for a complete genomic model explaining high-pressure adaptation.

Structure, Function, and Thermodynamics of Lactate Dehydrogenases from Humans and the Malaria Parasite P. falciparum

Temperature adaptation is ubiquitous among all living organisms, yet the molecular basis for this process remains poorly understood. It can be assumed that for parasite-host systems, the same enzymes found in both organisms respond to the same selection factor (human body temperature) with similar structural changes. Herein, we report the existence of a reversible temperature-dependent structural transition for the glycolytic enzyme lactate dehydrogenase (LDH) from the malaria parasite Plasmodium falciparum (pfLDH) and human heart (hhLDH) occurring in the temperature range of human fever. This transition is observed for LDHs from psychrophiles, mesophiles, and moderate thermophiles in their operating temperature range. Thermodynamic analysis reveals unique thermodynamic signatures of the LDH-substrate complexes defining a specific temperature range to which human LDH is adapted and parasite LDH is not, despite their common mesophilic nature. The results of spectroscopic analysis combined with the available crystallographic data reveal the existence of an active center within pfLDH that imparts psychrophilic structural properties to the enzyme. This center consists of two pockets, one formed by the five amino acids (5AA insert) within the substrate specificity loop and the other by the active site, that mutually regulate one another in response to temperature and induce structural and functional changes in the Michaelis complex. Our findings pave the way toward a new strategy for malaria treatments and drug design using therapeutic agents that inactivate malarial LDH selectively at a specific temperature range of the cyclic malaria paroxysm.

Biochemical, structural and dynamical studies reveal strong differences in the thermal-dependent allosteric behavior of two extremophilic lactate dehydrogenases

In this work, we combined biochemical and structural investigations with molecular dynamics (MD) simulations to analyze the very different thermal-dependent allosteric behavior of two lactate dehydrogenases (LDH) from thermophilic bacteria. We found that the enzyme from Petrotoga mobilis (P. mob) necessitates an absolute requirement of the allosteric effector (fructose 1, 6-bisphosphate) to ensure functionality. In contrast, even without allosteric effector, the LDH from Thermus thermophilus (T. the) is functional when the temperature is raised. We report the crystal structure of P. mob LDH in the Apo state solved at 1.9 Å resolution. We used this structure and the one from T. the, obtained previously, as a starting point for MD simulations at various temperatures. We found clear differences between the thermal dynamics, which accounts for the behavior of the two enzymes. Our work demonstrates that, within an allosteric enzyme, some areas act as local gatekeepers of signal transmission, allowing the enzyme to populate either the T-inactive or the R-active states with different degrees of stringency.

Phylogenetic, functional and structural characterization of a GH10 xylanase active at extreme conditions of temperature and alkalinity

Endoxylanases active under extreme conditions of temperature and alkalinity can replace the use of highly pollutant chemicals in the pulp and paper industry. Searching for enzymes with these properties, we carried out a comprehensive bioinformatics study of the GH10 family. The phylogenetic analysis allowed the construction of a radial cladogram in which protein sequences putatively ascribed as thermophilic and alkaliphilic appeared grouped in a well-defined region of the cladogram, designated TAK Cluster. One among five TAK sequences selected for experimental analysis (Xyn11) showed extraordinary xylanolytic activity under simultaneous conditions of high temperature (90 °C) and alkalinity (pH 10.5). Addition of a carbohydrate binding domain (CBM2) at the C-terminus of the protein sequence further improved the activity of the enzyme at high pH. Xyn11 structure, which has been solved at 1.8 Å resolution by X-ray crystallography, reveals an unusually high number of hydrophobic, ionic and hydrogen bond atomic interactions that could account for the enzyme’s extremophilic nature.

Electrostimulation of hyperthermophile Thermotoga neapolitana cultures

Hyperthermophile bioelectrochemical systems are seldom investigated although their superior control of microbial consortium and thermodynamic advantages. Hyperthermophilic Thermotogales, for instance, are able to produce hydrogen and lactic acid from wastes better than mesophilic bacteria. Here, the electrostimulation of Thermotoga neapolitana in single-chamber electrochemical bioreactors is studied. The glucose fermentation under CO₂ pressure, as model metabolism, was tested at 80 °C. Results show that a dynamic polarization (±0.8 to ±1.2 V) drives glucose fermentation and biofilm stasis on electrodes. Under this condition, production of lactic acid (33 vs 12 mM) and yields of acetate and hydrogen (with lactic/acetic acid ratio of 1.18) were higher than those achieved with static polarization or open-circuit. Dynamic polarization is possibly exploitable to stimulate T. neapolitana in a hyperthermophile electrochemical system for various applications including control of power-to-gas processes or production of value-added products (hydrogen and lactic acid) from sugary wastes.

Evolution of Protein Structure and Stability in Global Warming

This review focuses on the molecular signatures of protein structures in relation to evolution and survival in global warming. It is based on the premise that the power of evolutionary selection may lead to thermotolerant organisms that will repopulate the planet and continue life in general, but perhaps with different kinds of flora and fauna. Our focus is on molecular mechanisms, whereby known examples of thermoresistance and their physicochemical characteristics were noted. A comparison of interactions of diverse residues in proteins from thermophilic and mesophilic organisms, as well as reverse genetic studies, revealed a set of imprecise molecular signatures that pointed to major roles of hydrophobicity, solvent accessibility, disulfide bonds, hydrogen bonds, ionic and π-electron interactions, and an overall condensed packing of the higher-order structure, especially in the hydrophobic regions. Regardless of mutations, specialized protein chaperones may play a cardinal role. In evolutionary terms, thermoresistance to global warming will likely occur in stepwise mutational changes, conforming to the molecular signatures, such that each "intermediate" fits a temporary niche through punctuated equilibrium, while maintaining protein functionality. Finally, the population response of different species to global warming may vary substantially, and, as such, some may evolve while others will undergo catastrophic mass extinction.

Heat induces end to end repetitive association in P. furiosus L-asparaginase which enables its thermophilic property

It remains undeciphered how thermophilic enzymes display enhanced stability at elevated temperatures. Taking l-asparaginase from P. furiosus (PfA) as an example, we combined scattering shapes deduced from small-angle X-ray scattering (SAXS) data at increased temperatures with symmetry mates from crystallographic structures to find that heating caused end-to-end association. The small contact point of self-binding appeared to be enabled by a terminal short β-strand in N-terminal domain, Leu¹⁷⁹-Val-Val-Asn¹⁸² (LVVN). Interestingly, deletion of this strand led to a defunct enzyme, whereas suplementation of the peptide LVVN to the defunct enzyme restored structural frameworkwith mesophile-type functionality. Crystal structure of the peptide-bound defunct enzyme showed that one peptide ispresent in the same coordinates as in original enzyme, explaining gain-of lost function. A second peptide was seen bound to the protein at a different location suggesting its possible role in substrate-free molecular-association. Overall, we show that the heating induced self-assembly of native shapes of PfA led to an apparent super-stable assembly.

Aggregation-promoting conditions necessary to create the complexes by acylphosphatase from the hyperthermophile Sulfolobus solfataricus

The structural transition from the globular to the amyloid form of proteins requires aggregation-promoting conditions. The protein example of this category is acylphosphatase from the hyperthermophile Sulfolobus solfataricus. This protein represents a structure with a well-defined hydrophobic core. This is why the complexation (including oligomerization) of this protein is of low probability. The chain fragment participating in aggregation in comparison to the status with respect to the fuzzy oil drop model is discussed in this paper.

Biosurfactant-Producing Capability and Prediction of Functional Genes Potentially Beneficial to Microbial Enhanced Oil Recovery in Indigenous Bacterial Communities of an Onshore Oil Reservoir

Microbial enhanced oil recovery (MEOR) is a bio-based technology with economic and environmental benefits. The success of MEOR depends greatly on the types and characteristics of indigenous microbes. The aim of this study was to evaluate the feasibility of applying MEOR at Mae Soon Reservoir, an onshore oil reservoir experiencing a decline in its production rate. We investigated the capability of the reservoir's bacteria to produce biosurfactants, and evaluated the potentials of uncultured indigenous bacteria to support MEOR by means of prediction of MEOR-related functional genes, based on a set of metagenomic 16s rRNA gene data. The biosurfactant-producing bacteria isolated from the oil-bearing sandstones from the reservoir belonged to one species: Bacillus licheniformis, with one having the ability to decrease surface tension from 72 to 32 mN/m. Gene sequences responsible for biosurfactant (licA3), lipase (lipP1) and catechol 2,3-dioxygenase (C23O) were detected in these isolates. The latter two, and other genes encoding MEOR-related functional proteins such as enoyl-CoA hydratase and alkane 1-monooxygenase, were predicted in the bacterial communities residing the reservoir's sandstones. Exposure of these sandstones to nutrients, consisting of KNO₃ and NaH₂PO₄, resulted in an increase in the proportions of some predicted functional genes. These results indicated the potentials of MEOR application at Mae Soon site. Using the approaches demonstrated in this study would also assist evaluation of the feasibility of applying MEOR in oil reservoirs, which may be enhanced by an appropriate nutrient treatment.

Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations

Proteins are sensitive to temperature, and abrupt changes in the normal temperature conditions can have a profound impact on both structure and function, leading to protein unfolding. However, the adaptation of certain organisms to extreme conditions raises questions about the structural features that permit the structure and function of proteins to be preserved under these adverse conditions. To gain insight into the molecular basis of protein thermostability in the globin family, we have examined three representative examples: human neuroglobin, horse heart myoglobin, and Drosophila hemoglobin, which differ in their melting temperatures and coordination states of the heme iron in the absence of external ligands. In order to elucidate the possible mechanisms that govern the thermostability of these proteins, microsecond-scale classical molecular dynamics simulations were performed at different temperatures. Structural fluctuations and essential dynamics were analyzed, indicating that the flexibility of the CD region, which includes the two short C and D helixes and the connecting CD loop, is directly related to the thermostability. We observed that a larger inherent flexibility of the protein produces higher thermostability, probably concentrating the thermal fluctuations observed at high temperature in flexible regions, preventing unfolding. Globally, the results of this work improve our understanding of thermostability in the globin family.

Thermodynamic and Structural Adaptation Differences between the Mesophilic and Psychrophilic Lactate Dehydrogenases

The thermodynamics of substrate binding and enzymatic activity of a glycolytic enzyme, lactate dehydrogenase (LDH), from both porcine heart, phLDH (Sus scrofa; a mesophile), and mackerel icefish, cgLDH (Chamapsocephalus gunnari; a psychrophile), were investigated. Using a novel and quite sensitive fluorescence assay that can distinguish protein conformational changes close to and distal from the substrate binding pocket, a reversible global protein structural transition preceding the high-temperature transition (denaturation) was surprisingly found to coincide with a marked change in enzymatic activity for both LDHs. A similar reversible structural transition of the active site structure was observed for phLDH but not for cgLDH. An observed lower substrate binding affinity for cgLDH compared to that for phLDH was accompanied by a larger contribution of entropy to ΔG, which reflects a higher functional plasticity of the psychrophilic cgLDH compared to that of the mesophilic phLDH. The natural osmolyte, trimethylamine N-oxide (TMAO), increases stability and shifts all structural transitions to higher temperatures for both orthologs while simultaneously reducing catalytic activity. The presence of TMAO causes cgLDH to adopt catalytic parameters like those of phLDH in the absence of the osmolyte. Our results are most naturally understood within a model of enzyme dynamics whereby different conformations of the enzyme that have varied catalytic parameters (i.e., binding and catalytic proclivity) and whose population profiles are temperature-dependent and influenced by osmolytes interconvert among themselves. Our results also show that adaptation can be achieved by means other than gene mutations and complements the synchronic evolution of the cellular milieu.

Is the hydrophobic core a universal structural element in proteins?

The hydrophobic core, when subjected to analysis based on the fuzzy oil drop model, appears to be a universal structural component of proteins irrespective of their secondary, supersecondary, and tertiary conformations. A study has been performed on a set of nonhomologous proteins representing a variety of CATH categories. The presence of a well-ordered hydrophobic core has been confirmed in each case, regardless of the protein’s biological function, chain length or source organism. In light of fuzzy oil drop (FOD) analysis, various supersecondary forms seem to share a common structural factor in the form of a hydrophobic core, emerging either as part of the whole protein or a specific domain. The variable status of individual folds with respect to the FOD model reflects their propensity for conformational changes, frequently associated with biological function. Such flexibility is expressed as variable stability of the hydrophobic core, along with specific encoding of potential conformational changes which depend on the properties of helices and β-folds.

Differential effects of metal ions on TCDD-induced cytotoxicity and cytochrome P4501A1 gene expression in a zebrafish liver (ZFL) cell-line

Trace metal ions and trace organic compounds are common co-contaminants in the environment that pose risks to human health. We evaluated the effects of four metal ions (As(3+), Cu(2+), Hg(2+), and Zn(2+)) on 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) induced cytotoxicity and the expression of the cytochrome P4501A1 gene (cyp1a1) in the zebrafish liver (ZFL) cell line. A metal accumulation study showed that Cu and Zn did not accumulate in ZFL cells. However, As and Hg did accumulate, which resulted in the inhibition of TCDD-mediated induction of cyp1a1 mRNA and protein expression, and 7-ethoxyresorufin O-deethylase activity. A luciferase assay showed that both As(3+) and Hg(2+) inhibited the TCDD-induced activity of gene constructs containing either synthetic 3XRE or a distal cyp1a1 promoter region, implying that the decreased levels of TCDD-induced cyp1a1 were due to transcriptional effects. A proteomic study showed that the toxic effects of As(3+) might be due to changes in cellular metabolic processes, the cellular stimulation response and the cellular redox state in ZFL cells.

Group 3 LEA protein model peptides protect enzymes against desiccation stress

We tested whether model peptides for group 3 late embryogenesis abundant (G3LEA) proteins, which we developed previously, are capable of maintaining the catalytic activities of enzymes dried in their presence. Three different peptides were compared: 1) PvLEA-22, which consists of two tandem repeats of the 11-mer motif found in G3LEA proteins from an African sleeping chironomid; 2) PvLEA-44, which is made of four tandem repeats of the same 11-mer motif; and 3) a peptide whose amino acid composition is the same as that of PvLEA-22, but whose sequence is scrambled. We selected two enzymes, lactate dehydrogenase (LDH) and β-d-galactosidase (BDG), as targets because they have different isoelectric point (pI) values, in the alkaline and acidic range, respectively. While these enzymes were almost inactivated when dried alone, their catalytic activity was preserved at ≥70% of native levels in the presence of any of the above three peptides. This degree of protection is comparable to that conferred by several full-length G3LEA proteins, as reported previously for LDH. Interestingly, the protective activity of the peptides was enhanced slightly when they were mixed with trehalose, especially when the molar content of the peptides was low. On the basis of these results, the G3LEA model peptides show promise as protectants for the dry preservation of enzymes/proteins with a wide range of pI values.

A Streamlined, Automated Protocol for the Production of Milligram Quantities of Untagged Recombinant Rat Lactate Dehydrogenase A Using ÄKTAxpressTM

We developed an efficient, automated 2-step purification protocol for the production of milligram quantities of untagged recombinant rat lactate dehydrogenase A (rLDHA) from E. coli, using the ÄKTAxpress™ chromatography system. Cation exchange followed by size exclusion results in average final purity in excess of 93% and yields ~ 14 milligrams per 50 ml of original cell culture in EnPresso B media, in under 8 hrs, including all primary sample processing and column equilibration steps. The protein is highly active and coherent biophysically and a viable alternative to the more problematic human homolog for structural and ligand-binding studies; an apo structure of untagged rLDHA was solved to a resolution 2.29 Å (PDB ID 5ES3). Our automated methodology uses generic commercially available pre-packed columns and simple buffers, and represents a robust standard method for the production of milligram amounts of untagged rLDHA, facilitating a novel fragment screening approach for new inhibitors.

The D-Lactate Dehydrogenase from Sporolactobacillus inulinus Also Possessing Reversible Deamination Activity

Hydroxyacid dehydrogenases are responsible for the conversion of 2-keto acids to 2-hydroxyacids and have a wide range of biotechnological applications. In this study, a D-lactate dehydrogenase (D-LDH) from a Sporolactobacillus inulinus strain was experimentally verified to have both the D-LDH and glutamate dehydrogenase (GDH) activities (reversible deamination). The catalytic mechanism was demonstrated by identification of key residues from the crystal structure analysis and site-directed mutagenesis. The Arg²³⁴ and Gly⁷⁹ residues of this enzyme play a significant role in both D-LDH and GDH activities. His²⁹⁵ and Phe²⁹⁸ in DLDH744 were identified to be key residues for lactate dehydrogenase (LDH) activity only whereas Tyr¹⁰¹ is a unique residue that is critical for GDH activity. Characterization of the biochemical properties contributes to understanding of the catalytic mechanism of this novel D-lactate dehydrogenase enzyme.

Structural and catalytic effects of proline substitution and surface loop deletion in the extended active site of human carbonic anhydrase II

Bioengineering of a thermophilic enzyme starting from a mesophilic scaffold has proven to be a significant challenge, as several stabilizing elements have been proposed to be the foundation of thermal stability, including disulfide bridges, surface loop reduction, ionic pair networks, proline substitutions and aromatic clusters. This study emphasizes the effect of increasing the rigidity of human carbonic anhydrase II (HCA II; EC 4.2.1.1) via incorporation of proline residues at positions 170 and 234, which are located in surface loops that are able to accommodate restrictive main-chain conformations without rearrangement of the surrounding peptide backbone. Additionally, the effect of the compactness of HCA II was examined by deletion of a surface loop (residues 230-240) that had been previously identified as a possible source of thermal stability for the hyperthermophilic carbonic anhydrase isolated from the bacterium Sulfurihydrogenibium yellowstonense YO3AOP1. Differential scanning calorimetry analysis of these HCA II variants revealed that these structural modifications had a minimum effect on the thermal stability of the enzyme, while kinetic studies showed unexpected effects on the catalytic efficiency and proton transfer rates. X-ray crystallographic analysis of these HCA II variants showed that the electrostatic potential and configuration of the highly acidic loop (residues 230-240) play an important role in its high catalytic activity. Based on these observations and previous studies, a picture is emerging of the various components within the general structural architecture of HCA II that are key to stability. These elements may provide blueprints for rational thermal stability engineering of other enzymes.

DATABASE

Structural data have been submitted to the Protein Data Bank under accession numbers 4QK1 (K170P), 4QK2 (E234P) and 4QK3 (Δ230-240).

Backbone circularization of Bacillus subtilis family 11 xylanase increases its thermostability and its resistance against aggregation

While using a serine (S) as linker for circularization increases the thermostability, a longer linker (RGKCWE) leads to reduced aggregation after heat shock at elevated temperatures.

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.

Thermostable artificial enzyme isolated by in vitro selection

Artificial enzymes hold the potential to catalyze valuable reactions not observed in nature. One approach to build artificial enzymes introduces mutations into an existing protein scaffold to enable a new catalytic activity. This process commonly results in a simultaneous reduction of protein stability as an undesired side effect. While protein stability can be increased through techniques like directed evolution, care needs to be taken that added stability, conversely, does not sacrifice the desired activity of the enzyme. Ideally, enzymatic activity and protein stability are engineered simultaneously to ensure that stable enzymes with the desired catalytic properties are isolated. Here, we present the use of the in vitro selection technique mRNA display to isolate enzymes with improved stability and activity in a single step. Starting with a library of artificial RNA ligase enzymes that were previously isolated at ambient temperature and were therefore mostly mesophilic, we selected for thermostable active enzyme variants by performing the selection step at 65°C. The most efficient enzyme, ligase 10C, was not only active at 65°C, but was also an order of magnitude more active at room temperature compared to related enzymes previously isolated at ambient temperature. Concurrently, the melting temperature of ligase 10C increased by 35 degrees compared to these related enzymes. While low stability and solubility of the previously selected enzymes prevented a structural characterization, the improved properties of the heat-stable ligase 10C finally allowed us to solve the three-dimensional structure by NMR. This artificial enzyme adopted an entirely novel fold that has not been seen in nature, which was published elsewhere. These results highlight the versatility of the in vitro selection technique mRNA display as a powerful method for the isolation of thermostable novel enzymes.

The core of allosteric motion in Thermus caldophilus L-lactate dehydrogenase

For Thermus caldophilus L-lactate dehydrogenase (TcLDH), fructose 1,6-bisphosphate (FBP) reduced the pyruvate S(0.5) value 10(3)-fold and increased the V(max) value 4-fold at 30 °C and pH 7.0, indicating that TcLDH has a much more T state-sided allosteric equilibrium than Thermus thermophilus L-lactate dehydrogenase, which has only two amino acid replacements, A154G and H179Y. The inactive (T) and active (R) state structures of TcLDH were determined at 1.8 and 2.0 Å resolution, respectively. The structures indicated that two mobile regions, MR1 (positions 172-185) and MR2 (positions 211-221), form a compact core for allosteric motion, and His(179) of MR1 forms constitutive hydrogen bonds with MR2. The Q4(R) mutation, which comprises the L67E, H68D, E178K, and A235R replacements, increased V(max) 4-fold but reduced pyruvate S(0.5) only 5-fold in the reaction without FBP. In contrast, the P2 mutation, comprising the R173Q and R216L replacements, did not markedly increase V(max), but 10(2)-reduced pyruvate S(0.5), and additively increased the FBP-independent activity of the Q4(R) enzyme. The two types of mutation consistently increased the thermal stability of the enzyme. The MR1-MR2 area is a positively charged cluster, and its center approaches another positively charged cluster (N domain cluster) across the Q-axis subunit interface by 5 Å, when the enzyme undergoes the T to R transition. Structural and kinetic analyses thus revealed the simple and unique allosteric machinery of TcLDH, where the MR1-MR2 area pivotally moves during the allosteric motion and mediates the allosteric equilibrium through electrostatic repulsion within the protein molecule.

Molecular basis of thermal stability in truncated (2/2) hemoglobins

BACKGROUND

Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures.

METHODS

We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position.

RESULTS

The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible.

CONCLUSIONS

This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability.

GENERAL SIGNIFICANCE

These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

Purification and properties of a monomeric lactate dehydrogenase from yak Hypoderma sinense larva

The objective of the present study was to study the characteristics of lactate dehydrogenase (LDH) from Hypoderma sinense larva. H. sinense larvae were collected from yak (Bos grunniens) and identified by a PCR-RFLP method. Analysis of LDH activity showed that the total LDH activity in H. sinense larva was negatively correlated with the length of larva. Polyacrylamide gel electrophoresis of the extracts of H. sinense larvae revealed one band of LDH, which was then purified by affinity chromatography and gel filtration. This enzyme showed an approximately 36 kDa band on SDS-gel under both reducing and non-reducing conditions, in addition, size exclusion chromatography analysis showed that its molecular weight was smaller than bovine serum albumin (67 kDa), indicating that it contains only one subunit. Michaelis constants (Km) values assay revealed that LDH from H. sinense larva showed significantly lower Km for lactate than other animals. LDH of H. sinense larva was stable at 60 °C for 15 min, and also exhibited high catalytic efficiency in a wide range of pH. HgCl₂ at the concentration of 0.1mM significantly decreased the activity of LDH from H. sinense larva but not at the concentration of 0.01 mM. The results of the present study demonstrate that LDH from H. sinense larva is a thermal stable and pH insensitive enzyme suitable for catalyzing both forward and reverse reactions.

Molecular dynamics of a thermostable multicopper oxidase from Thermus thermophilus HB27: structural differences between the apo and holo forms

Molecular dynamic (MD) simulations have been performed on Tth-MCO, a hyperthermophilic multicopper oxidase from thermus thermophilus HB27, in the apo as well as the holo form, with the aim of exploring the structural dynamic properties common to the two conformational states. According to structural comparison between this enzyme and other MCOs, the substrate in process to electron transfer in an outer-sphere event seems to transiently occupy a shallow and overall hydrophobic cavity near the Cu type 1 (T1Cu). The linker connecting the β-strands 21 and 24 of the second domain (loop (β21–β24)_D2) has the same conformation in both states, forming a flexible lid at the entrance of the electron-transfer cavity. Loop (β21–β24)_D2 has been tentatively assigned a role occluding the access to the electron-transfer site. The dynamic of the loop (β21–β24)_D2 has been investigated by MD simulation, and results show that the structures of both species have the same secondary and tertiary structure during almost all the MD simulations. In the simulation, loop (β21–β24)_D2 of the holo form undergoes a higher mobility than in the apo form. In fact, loop (β21–β24)_D2 of the holo form experiences a conformational change which enables exposure to the electron-transfer site (open conformation), while in the apo form the opposite effect takes place (closed conformation). To confirm the hypothesis that the open conformation might facilitate the transient electron-donor molecule occupation of the site, the simulation was extended another 40 ns with the electron-donor molecule docked into the protein cavity. Upon electron-donor molecule stabilization, loops near the cavity reduce their mobility. These findings show that coordination between the copper and the protein might play an important role in the general mobility of the enzyme, and that the open conformation seems to be required for the electron transfer process to T1Cu.

A molecular design that stabilizes active state in bacterial allosteric L-lactate dehydrogenases

l-Lactate dehydrogenase (l-LDH) of Lactobacillus casei (LCLDH) is a typical bacterial allosteric l-LDH that requires fructose 1,6-bisphosphate (FBP) for its enzyme activity. A mutant LCLDH was designed to introduce an inter-subunit salt bridge network at the Q-axis subunit interface, mimicking Lactobacillus pentosus non-allosteric l-LDH (LPLDH). The mutant LCLDH exhibited high catalytic activity with hyperbolic pyruvate saturation curves independently of FBP, and virtually the equivalent K(m) and V(m) values at pH 5.0 to those of the fully activated wild-type enzyme with FBP, although the K(m) value was slightly improved with FBP or Mn(2+) at pH 7.0. The mutant enzyme exhibited a markedly higher apparent denaturating temperature (T(1/2)) than the wild-type enzyme in the presence of FBP, but showed an even lower T(1/2) without FBP, where it exhibited higher activation enthalpy of inactivation (ΔH(‡)). This result is consistent with the fact that the active state is more unstable than the inactive state in allosteric equilibrium of LCLDH. The LPLDH-like network appears to be conserved in many bacterial non-allosteric l-LDHs and dimeric l-malate dehydrogenases, and thus to be a key for the functional divergence of bacterial l-LDHs during evolution.

PDBsum new things

PDBsum (http://www.ebi.ac.uk/pdbsum) provides summary information about each experimentally determined structural model in the Protein Data Bank (PDB). Here we describe some of its most recent features, including figures from the structure's key reference, citation data, Pfam domain diagrams, topology diagrams and protein-protein interactions. Furthermore, it now accepts users' own PDB format files and generates a private set of analyses for each uploaded structure.

Evidence for subdomain flexibility in Drosophila melanogaster acetylcholinesterase

The catalytic domain of the acetylcholinesterases is composed of a single polypeptide chain, the folding of which determines two subdomains. We have linked these two subdomains by mutating two residues, I327 and D375, to cysteines, to form a disulfide bridge. As a consequence, the hydrodynamic radius of the protein was reduced, suggesting that there is some flexibility in the subdomain connection. In addition to the smaller size, the mutated protein is more stable than the wild-type protein. Therefore, the flexibility between the two domains is a weak point in terms of protein stability. As expected from the location of the disulfide bond at the rim of the active site, the kinetic studies show that it affects interactions with peripheral ligands and the entrance of some of the bulkier substrates, like o-nitrophenyl acetate. In addition, the mutations affect the catalytic step for o-nitrophenyl acetate and phosphorylation by organophosphates, suggesting that this movement between the two subdomains is connected with the cooperativity between the peripheral and catalytic sites.

Activity, stability and structural studies of lactate dehydrogenases adapted to extreme thermal environments

Lactate dehydrogenase (LDH) catalyzes the conversion of pyruvate to lactate with concomitant oxidation of NADH during the last step in anaerobic glycolysis. In the present study, we present a comparative biochemical and structural analysis of various LDHs adapted to function over a large temperature range. The enzymes were from Champsocephalus gunnari (an Antarctic fish), Deinococcus radiodurans (a mesophilic bacterium) and Thermus thermophilus (a hyperthermophilic bacterium). The thermodynamic activation parameters of these LDHs indicated that temperature adaptation from hot to cold conditions was due to a decrease in the activation enthalpy and an increase in activation entropy. The crystal structures of these LDHs have been solved. Pairwise comparisons at the structural level, between hyperthermophilic versus mesophilic LDHs and mesophilic versus psychrophilic LDHs, have revealed that temperature adaptation is due to a few amino acid substitutions that are localized in critical regions of the enzyme. These substitutions, each having accumulating effects, play a role in either the conformational stability or the local flexibility or in both. Going from hot- to cold-adapted LDHs, the various substitutions have decreased the number of ion pairs, reduced the size of ionic networks, created unfavorable interactions involving charged residues and induced strong local disorder. The analysis of the LDHs adapted to extreme temperatures shed light on how evolutionary processes shift the subtle balance between overall stability and flexibility of an enzyme.

Mechanisms for stabilisation and the maintenance of solubility in proteins from thermophiles

Background

The database of protein structures contains representatives from organisms with a range of growth temperatures. Various properties have been studied in a search for the molecular basis of protein adaptation to higher growth temperature. Charged groups have emerged as key distinguishing factors for proteins from thermophiles and mesophiles.

Results

A dataset of 291 thermophile-derived protein structures is compared with mesophile proteins. Calculations of electrostatic interactions support the importance of charges, but indicate that increases in charge contribution to folded state stabilisation do not generally correlate with the numbers of charged groups. Relative propensities of charged groups vary, such as the substitution of glutamic for aspartic acid sidechains. Calculations suggest an energetic basis, with less dehydration for longer sidechains. Most other properties studied show weak or insignificant separation of proteins from moderate thermophiles or hyperthermophiles and mesophiles, including an estimate of the difference in sidechain rotameric entropy upon protein folding. An exception is increased burial of alanine and proline residues and decreased burial of phenylalanine, methionine, tyrosine and tryptophan in hyperthermophile proteins compared to those from mesophiles.

Conclusion

Since an increase in the number of charged groups for hyperthermophile proteins is separable from charged group contribution to folded state stability, we hypothesise that charged group propensity is important in the context of protein solubility and the prevention of aggregation. Accordingly we find some separation between mesophile and hyperthermophile proteins when looking at the largest surface patch that does not contain a charged sidechain. With regard to our observation that aromatic sidechains are less buried in hyperthermophile proteins, further analysis indicates that the placement of some of these groups may facilitate the reduction of folding fluctuations in proteins of the higher growth temperature organisms.

The crystal structure of a hyperthermostable subfamily II isocitrate dehydrogenase from Thermotoga maritima

Isocitrate dehydrogenase (IDH) from the hyperthermophile Thermotoga maritima (TmIDH) catalyses NADP+- and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. It belongs to the beta-decarboxylating dehydrogenase family and is the only hyperthermostable IDH identified within subfamily II. Furthermore, it is the only IDH that has been characterized as both dimeric and tetrameric in solution. We solved the crystal structure of the dimeric apo form of TmIDH at 2.2 A. The R-factor of the refined model was 18.5% (R(free) 22.4%). The conformation of the TmIDH structure was open and showed a domain rotation of 25-30 degrees compared with closed IDHs. The separate domains were found to be homologous to those of the mesophilic mammalian IDHs of subfamily II and were subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that stabilization of the N- and C-termini via long-range electrostatic interactions were important for the higher thermostability of TmIDH. Moreover, the number of intra- and intersubunit ion pairs was higher and the ionic networks were larger compared with the mesophilic IDHs. Other factors likely to confer higher stability in TmIDH were a less hydrophobic and more charged accessible surface, a more hydrophobic subunit interface, more hydrogen bonds per residue and a few loop deletions. The residues responsible for the binding of isocitrate and NADP+ were found to be highly conserved between TmIDH and the mammalian IDHs and it is likely that the reaction mechanism is the same.

Structural basis of the substrate subsite and the highly thermal stability of xylanase 10B from Thermotoga maritima MSB8

The crystal structure of xylanase 10B from Thermotoga maritima MSB8 (TmxB), a hyperthermostable xylanase, has been solved in its native form and in complex with xylobiose or xylotriose at 1.8 A resolution. In order to gain insight into the substrate subsite and the molecular features for thermal stability, we compared TmxB with family 10 xylanase structures from nine microorganisms. As expected, TmxB folds into a (beta/alpha)8-barrel structure, which is common among the glycoside hydrolase family 10. The enzyme active site and the environment surrounding the xylooligosaccharide of TmxB are highly similar to those of family 10 xylanases. However, only two xylose moieties were found in its binding pocket from the TmxB-xylotriose complex structure. This finding suggests that TmxB could be a potential biocatalyst for the large-scale production of xylobiose. The result of structural analyses also indicated that TmxB possesses some additional features that account for its thermostability. In particular, clusters of aromatic residues together with a lack of exposed hydrophobic residues are characteristic of the TmxB structure. TmxB has also a significant number of ion pairs on the protein surface that are not found in other thermophilic family 10 xylanases.

Protein molecular function prediction by Bayesian phylogenomics

We present a statistical graphical model to infer specific molecular function for unannotated protein sequences using homology. Based on phylogenomic principles, SIFTER (Statistical Inference of Function Through Evolutionary Relationships) accurately predicts molecular function for members of a protein family given a reconciled phylogeny and available function annotations, even when the data are sparse or noisy. Our method produced specific and consistent molecular function predictions across 100 Pfam families in comparison to the Gene Ontology annotation database, BLAST, GOtcha, and Orthostrapper. We performed a more detailed exploration of functional predictions on the adenosine-5'-monophosphate/adenosine deaminase family and the lactate/malate dehydrogenase family, in the former case comparing the predictions against a gold standard set of published functional characterizations. Given function annotations for 3% of the proteins in the deaminase family, SIFTER achieves 96% accuracy in predicting molecular function for experimentally characterized proteins as reported in the literature. The accuracy of SIFTER on this dataset is a significant improvement over other currently available methods such as BLAST (75%), GeneQuiz (64%), GOtcha (89%), and Orthostrapper (11%). We also experimentally characterized the adenosine deaminase from Plasmodium falciparum, confirming SIFTER's prediction. The results illustrate the predictive power of exploiting a statistical model of function evolution in phylogenomic problems. A software implementation of SIFTER is available from the authors.

Three amino acids that are critical to formation and stability of the P22 tailspike trimer

The P22 tailspike protein folds by forming a folding competent monomer species that forms a dimeric, then a non-native trimeric (protrimer) species by addition of folding competent monomers. We have found three residues, R549, R563, and D572, which play a critical role in both the stability of the native tailspike protein and assembly and maturation of the protrimer. King and colleagues reported previously that substitution of R563 to glutamine inhibited protrimer formation. We now show that the R549Q and R563K variants significantly delay the protrimer-to-trimer transition both in vivo and in vitro. Previously, variants that destabilize intermediates have shown wild-type chemical stability. Interestingly, both the R549Q and R563K variants destabilize the tailspike trimer in guanidine denaturation studies, indicating that they represent a new class of tailspike folding variants. R549Q has a midpoint of unfolding at 3.2M guanidine, compared to 5.6M for the wild-type tailspike protein, while R563K has a midpoint of unfolding of 1.8 M. R549Q and R563K also denature over a broader pH range than the wild-type tailspike protein and both proteins have increased sensitivity to pH during refolding, suggesting that both residues are involved in ionic interactions. Our model is that R563 and D572 interact to stabilize the adjacent turn, aiding the assembly of the dimer and protrimer species. We believe that the interaction between R563 and D572 is also critical following assembly of the protrimer to properly orient D572 in order to form a salt bridge with R549 during protrimer maturation.

Structural Basis for the Thermostability of Ferredoxin from the Cyanobacterium Mastigocladus laminosus

Plant-type ferredoxins (Fds) carry a single [2Fe-2S] cluster and serve as electron acceptors of photosystem I (PSI). The ferredoxin from the thermophilic cyanobacterium Mastigocladus laminosus displays optimal activity at 65 degrees C. In order to reveal the molecular factors that confer thermostability, the crystal structure of M.laminosus Fd (mFd) was determined to 1.25 A resolution and subsequently analyzed in comparison with four similar plant-type mesophilic ferredoxins. The topologies of the plant-type ferredoxins are similar, yet two structural determinants were identified that may account for differences in thermostability, a salt bridge network in the C-terminal region, and the flexible L1,2 loop that increases hydrophobic accessible surface area. These conclusions were verified by three mutations, i.e. substitution of L1,2 into a rigid beta-turn ((Delta)L1,2) and two point mutations (E90S and E96S) that disrupt the salt bridge network at the C-terminal region. All three mutants have shown reduced electron transfer (ET) capabilities and [2Fe-2S] stability at high temperatures in comparison to the wild-type mFd. The results have also provided new insights into the involvement of the L1,2 loop in the Fd interactions with its electron donor, the PSI complex.

Isocitrate dehydrogenase from the hyperthermophile Aeropyrum pernix: X-ray structure analysis of a ternary enzyme-substrate complex and thermal stability

Isocitrate dehydrogenase from Aeropyrum pernix (ApIDH) is a homodimeric enzyme that belongs to the beta-decarboxylating dehydrogenase family and is the most thermostable IDH identified. It catalyzes the NADP+ and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. We have solved the crystal structures of a native ApIDH at 2.2 A, a pseudo-native ApIDH at 2.1 A, and of ApIDH in complex with NADP+, Ca2+ and d-isocitrate at 2.3 A. The pseudo-native ApIDH is in complex with etheno-NADP+ which was located at the surface instead of in the active site revealing a novel adenine-nucleotide binding site in ApIDH. The native and the pseudo-native ApIDHs were found in an open conformation, whereas one of the subunits of the ternary complex was closed upon substrate binding. The closed subunit showed a domain rotation of 19 degrees compared to the open subunit. The binding of isocitrate in the closed subunit was identical with that of the binary complex of porcine mitochondrial IDH, whereas the binding of NADP+ was similar to that of the ternary complex of IDH from Escherichiacoli. The reaction mechanism is likely to be conserved in the different IDHs. A proton relay chain involving at least five solvent molecules, the 5'-phosphate group of the nicotinamide-ribose and a coupled lysine-tyrosine pair in the active site, is postulated as essential in both the initial and the final steps of the catalytic reaction of IDH. ApIDH was found to be highly homologous to the mesophilic IDHs and was subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that a disulfide bond at the N terminus and a seven-membered inter-domain ionic network at the surface are major determinants for the higher thermostability of ApIDH compared to EcIDH. Furthermore, the total number of ion pairs was dramatically higher in ApIDH compared to the mesophilic IDHs if a cutoff of 4.2 A was used. A calculated net charge of only +1 compared to -19 and -25 in EcIDH and BsIDH, respectively, suggested a high degree of electrostatic optimization, which is known to be an important determinant for increased thermostability.