Lessons in stability from thermophilic proteins

Deep generative modeling of temperature-dependent structural ensembles of proteins

Deep learning has revolutionized protein structure prediction, but capturing conformational ensembles and structural variability remains an open challenge. While molecular dynamics (MD) is the foundation method for simulating biomolecular dynamics, it is computationally expensive. Recently, deep learning models trained on MD have made progress in generating structural ensembles at reduced cost. However, they remain limited in modeling atomistic details and, crucially, incorporating the effect of environmental factors. Here, we present aSAM (atomistic structural autoencoder model), a latent diffusion model trained on MD to generate heavy atom protein ensembles. Unlike most methods, aSAM models atoms in a latent space, greatly facilitating accurate sampling of side chain and backbone torsion angle distributions. Additionally, we extended aSAM into the first reported transferable generator conditioned on temperature, named aSAMt. Trained on the large and open mdCATH dataset, aSAMt captures temperature-dependent ensemble properties and demonstrates generalization beyond training temperatures. By comparing aSAMt ensembles to long MD simulations of fast folding proteins, we find that high-temperature training enhances the ability of deep generators to explore energy landscapes. Finally, we also show that our MD-based aSAMt can already capture experimentally observed thermal behavior of proteins. Our work is a step towards generalizable ensemble generation to complement physics-based approaches.

Conformational Stability at Low Temperatures Using Single Protein Nanoaperture Optical Tweezers

Nanoaperture optical tweezers allow for trapping single proteins and detecting their conformational changes without modifying the protein, i.e., being free from labels or tethers. While past works have used laser heating as a way to vary the local temperature, this does not allow for probing of lower temperature values. Here we investigate the lower temperature dynamics of individual Bovine Serum Albumin (BSA) proteins with the help of a custom Peltier cooling stage. The BSA transitions between the normal (N) and fast (F) states. The normal form of BSA has a maximum occupancy at 21 ± 1 °C, which is interpreted as its maximum stability point for the compact N form with respect to the F form. In this way, it is possible to find the relative thermodynamic parameters of single proteins without requiring any modifications to the intrinsic structure.

Thermostability Improvement of the Chitinase from Bacillus circulans for Efficient Chitin Oligosaccharide Production via Computational Design

The chitinase A1 from Bacillus circulans WL-12 (BcChiA1) exhibits promising potential for producing chitin oligosaccharides (CHOs), while its application is limited by its poor thermal stability. In this study, a set of thermostable variants were obtained by modifying BcChiA1 using a comprehensive strategy based on a computer-aided design. A combination of five beneficial single-point mutations (S67G/K177R/A220V/N257Y/N271E) to BcChiA1 generated a markedly improved variant, Mu5. Mu5 exhibited a half-life of 295 min at 60 °C, which was 59 times higher than that of BcChiA1. Furthermore, Mu5 was reused for chitin conversion, releasing 86.14 ± 3.73 mM of CHOs after five reaction cycles. Molecular dynamics simulation and structural analysis revealed that these enhancements were driven by increased structural rigidity and compactness, resulting in a protein conformation that was less prone to thermal denaturation. This combined approach through computational design yielded a thermostable BcChiA1 variant, potentially facilitating CHOs production in economical way.

Amino acid sequence encodes protein abundance shaped by protein stability at reduced synthesis cost

Understanding what drives protein abundance is essential to biology, medicine, and biotechnology. Driven by evolutionary selection, an amino acid sequence is tailored to meet the required abundance of a proteome, underscoring the intricate relationship between sequence and functional demand. Yet, the specific role of amino acid sequences in determining proteome abundance remains elusive. Here we show that the amino acid sequence alone encodes over half of protein abundance variation across all domains of life, ranging from bacteria to mouse and human. With an attempt to go beyond predictions, we trained a manageable-size Transformer model to interpret latent factors predictive of protein abundances. Intuitively, the model's attention focused on the protein's structural features linked to stability and metabolic costs related to protein synthesis. To probe these relationships, we introduce MGEM (Mutation Guided by an Embedded Manifold), a methodology for guiding protein abundance through sequence modifications. We find that mutations which increase predicted abundance have significantly altered protein polarity and hydrophobicity, underscoring a connection between protein structural features and abundance. Through molecular dynamics simulations we revealed that abundance-enhancing mutations possibly contribute to protein thermostability by increasing rigidity, which occurs at a lower synthesis cost.

The Structural Stability of Membrane Proteins Revisited: Combined Thermodynamic and Spectral Phasor Analysis of SDS-induced Denaturation of a Thermophilic Cu(I)-transport ATPase

Assessing membrane protein stability is among the major challenges in protein science due to their inherent complexity, which complicates the application of conventional biophysical tools. In this work, sodium dodecyl sulfate-induced denaturation of AfCopA, a Cu(I)-transport ATPase from Archaeoglobus fulgidus, was explored using a combined model-free spectral phasor analysis and a model-dependent thermodynamic analysis. Decrease in tryptophan and 1-anilino-naphthalene-8-sulfonate fluorescence intensity, displacements in the spectral phasor space, and the loss of ATPase activity were reversibly induced by this detergent. Refolding from the SDS-induced denatured state yields an active enzyme that is functionally and spectroscopically indistinguishable from the native state of the protein. Phasor analysis of Trp spectra allowed us to identify two intermediate states in the SDS-induced denaturation of AfCopA, a result further supported by principal component analysis. In contrast, traditional thermodynamic analysis detected only one intermediate state, and including the second one led to overparameterization. Additionally, ANS fluorescence spectral analysis detected one more intermediate and a gradual change at the level of the hydrophobic transmembrane surface of the protein. Based on this evidence, a model for acquiring the native structure of AfCopA in a membrane-like environment is proposed.

Deciphering Structural Traits for Thermal and Kinetic Stability across Protein Family Evolution through Ancestral Sequence Reconstruction

Natural proteins are frequently marginally stable, and an increase in environmental temperature can easily lead to unfolding. As a result, protein engineering to improve protein stability is an area of intensive research. Nonetheless, since there is usually a high degree of structural homology between proteins from thermophilic organisms and their mesophilic counterparts, the identification of structural determinants for thermoadaptation is challenging. Moreover, in many cases, it has become clear that the success of stabilization strategies is often dependent on the evolutionary history of a protein family. In the last few years, the use of ancestral sequence reconstruction (ASR) as a tool for elucidation of the evolutionary history of functional traits of a protein family has gained strength. Here, we used ASR to trace the evolutionary pathways between mesophilic and thermophilic kinases that participate in the biosynthetic pathway of vitamin B1 in bacteria. By combining biophysics approaches, X-ray crystallography, and molecular dynamics simulations, we found that the thermal stability of these enzymes correlates with their kinetic stability, where the highest thermal/kinetic stability is given by an increase in small hydrophobic amino acids that allow a higher number of interatomic hydrophobic contacts, making this type of interaction the main support for stability in this protein architecture. The results highlight the potential benefits of using ASR to explore the evolutionary history of protein sequence and structure to identify traits responsible for the kinetic and thermal stability of any protein architecture.

Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

Differentiating stable and unstable protein using convolution neural network and molecular dynamics simulations

Protein stability is a critical aspect of molecular biology and biochemistry, hinges on an intricate balance of thermodynamic and structural factors. Determining protein stability is crucial for understanding and manipulating biological machineries, as it directly correlated with the protein function. Thus, this study delves into the intricacies of protein stability, highlighting its dependence on various factors, including thermodynamics, thermal conditions, and structural properties. Moreover, a notable focus is placed on the free energy change of unfolding (ΔGunfolding), change in heat capacity (ΔCp) with protein structural transition, melting temperature (Tm) and number of disulfide bonds, which are critical parameters in understanding protein stability. In this study, a machine learning (ML) predictive model was developed to estimate these four parameters using the primary sequence of the protein. The shortfall of available tools for protein stability prediction based on multiple parameters propelled the completion of this study. Convolutional Neural Network (CNN) with multiple layers was adopted to develop a more reliable ML model. Individual predictive models were prepared for each property, and all the prepared models showed results with high accuracy. The R2 (coefficient of determination) of these models were 0.79, 0.78, 0.92 and 0.92, respectively, for ΔG, ΔCp, Tm and disulfide bonds. A case study on stability analysis of two homologous proteins was presented to validate the results predicted through the developed model. The case study included in silico analysis of protein stability using molecular docking and molecular dynamic simulations. This validation study assured the accuracy of each model in predicting the stability associated properties. The alignment of physics-based principles with ML models has provided an opportunity to develop a fast machine learning solution to replace the computationally demanding physics-based calculations used to determine protein stability. Furthermore, this work provided valuable insights into the impact of mutation on protein stability, which has implications for the field of protein engineering. The source codes are available at https://github.com/Growdeatechnology.

Extremozyme-Based Biosensors for Environmental Pollution Monitoring: Recent Developments

Extremozymes combine high specificity and sensitivity with the ability to withstand extreme operational conditions. This work presents an overview of extremozymes that show potential for environmental monitoring devices and outlines the latest advances in biosensors utilizing these unique molecules. The characteristics of various extremozymes described so far are presented, underlining their stability and operational conditions that make them attractive for biosensing. The biosensor design is discussed based on the detection of photosynthesis-inhibiting herbicides as a case study. Several biosensors for the detection of pesticides, heavy metals, and phenols are presented in more detail to highlight interesting substrate specificity, applications or immobilization methods. Compared to mesophilic enzymes, the integration of extremozymes in biosensors faces additional challenges related to lower availability and high production costs. The use of extremozymes in biosensing does not parallel their success in industrial applications. In recent years, the "collection" of recognition elements was enriched by extremozymes with interesting selectivity and by thermostable chimeras. The perspectives for biosensor development are exciting, considering also the progress in genetic editing for the oriented immobilization of enzymes, efficient folding, and better electron transport. Stability, production costs and immobilization at sensing interfaces must be improved to encourage wider applications of extremozymes in biosensors.

Contribution of a surface salt bridge to the protein stability of deep-sea Shewanella benthica cytochrome c'

Two homologous cytochromes c', SBCP and SVCP, from deep-sea Shewanella benthica and Shewanella violacea respectively exhibit only nine surface amino acid substitutions, along with one at the N-terminus. Despite the small sequence difference, SBCP is thermally more stable than SVCP. Here, we examined the thermal stability of SBCP variants, each containing one of the nine substituted residues in SVCP, and found that the SBCP K87V variant was the most destabilized. We then determined the X-ray crystal structure of the SBCP K87V variant at a resolution of 2.1 Å. The variant retains a four-helix bundle structure similar to the wild-type, but notable differences are observed in the hydration structure around the mutation site. Instead of forming of the intrahelical salt bridge between Lys-87 and Asp-91 in the wild-type, a clathrate-like hydration around Val-87 through a hydrogen bond network with the nearby amino acid residues is observed. This network potentially enhances the ordering of surrounding water molecules, leading to an entropic destabilization of the protein. These results suggest that the unfavorable hydrophobic hydration environment around Val-87 and the inability to form the Asp-91-mediated salt bridge contribute to the observed difference in stability between SBCP and SVCP. These findings will be useful in future protein engineering for controlling protein stability through the manipulation of surface intrahelical salt bridges.

Rational Redesign of Chitosanase to Enhance Thermostability and Catalytic Activity to Produce Chitooligosaccharides with a Relatively High Degree of Polymerization

Chitooligosaccharides (hdpCOS) with a high degree of polymerization (hdp, DP 4-10) generally have greater biological activities than those of low-DP (ldp, DP 2-3) COS. Chitosanase from Bacillus amyloliquefaciens KCP2 (Csn46) can degrade chitosan to more hdpCOS at high temperature (70 °C), but low thermal stability at this temperature makes it unsuitable for industrial application; the wild-type enzyme can only produce COS (DP 2-4) at lower temperatures. Several thermostable mutants were obtained by modifying chitosanase using a comprehensive strategy based on a computer-aided mutant design. A combination of four beneficial single-point mutations (A129L/T175 V/K70T/D34G) to Csn46 was selected to obtain a markedly improved mutant, Mut4, with a half-life at 60 °C extended from 34.31 to 690.80 min, and the specific activity increased from 1671.73 to 3528.77 U/mg. Mut4 produced COS with DPs of 2-4 and 2-7 at 60 and 70 °C, respectively. Therefore, Mut4 has the potential to be applied to the industrial-scale preparation of hdpCOS with high biological activity.

Susceptibility of Cu(I) transport ATPases to sodium dodecyl sulfate. Relevance of the composition of the micellar phase

Sodium dodecyl sulfate (SDS) is a well-known protein denaturing agent. A less known property of this detergent is that it can activate or inactivate some enzymes at sub-denaturing concentrations. In this work we explore the effect of SDS on the ATPase activity of a hyper-thermophilic and a mesophilic Cu(I) ATPases reconstituted in mixed micelles of phospholipids and a non-denaturing detergent. An iterative procedure was used to evaluate the partition of SDS between the aqueous and the micellar phases, allowing to determine the composition of micelles prepared from phospholipid/detergent mixtures. The incubation of enzymes with SDS in the presence of different amounts of phospholipids reveals that higher SDS concentrations are required to obtain the same degree of inactivation when the initial concentration of phospholipids is increased. Remarkably, we found that, if represented as a function of the mole fraction of SDS in the micelle, the degree of inactivation obtained at different amounts of amphiphiles converges to a single inactivation curve. To interpret this result, we propose a simple model involving active and inactive enzyme molecules in equilibrium. This model allowed us to estimate the Gibbs free energy change for the inactivation process and its derivative with respect to the mole fraction of SDS in the micellar phase, the latter being a measure of the susceptibility of the enzyme to SDS. Our results showed that the inactivation free energy changes are similar for both proteins. Conversely, susceptibility to SDS is significantly lower for the hyperthermophilic ATPase, suggesting an inverse relation between thermophilicity and susceptibility to SDS.

Molecular Dynamics Simulations of HPr Proteins from a Thermophilic and a Mesophilic Organism: A Comparative Thermal Study

The histidine-containing phosphocarrier (HPr) is a monomeric protein conserved in Gram-positive bacteria, which may be of mesophilic or thermophilic nature. In particular, the HPr protein from the thermophilic organism B. stearothermophilus is a good model system for thermostability studies, since experimental data, such as crystal structure and thermal stability curves, are available. However, its unfolding mechanism at higher temperatures is yet unclear at a molecular level. Therefore, in this work, we researched the thermal stability of this protein using molecular dynamics simulations, subjecting it to five different temperatures during a time span of 1 μs. The analyses of the structural parameters and molecular interactions were compared with those of the mesophilic homologue HPr protein from B. subtilis. Each simulation was run in triplicate using identical conditions for both proteins. The results showed that the two proteins lose stability as the temperature increases, but the mesophilic structure is more affected. We found that the salt bridge network formed by the triad of Glu3-Lys62-Glu36 residues and the salt bridge made up of Asp79-Lys83 ion pair are key factors to keep stable the thermophilic protein, maintaining the hydrophobic core protected and the structure packed. In addition, these molecular interactions neutralize the negative surface charge, acting as "natural molecular staples".

Computational Insights and In Silico Characterization of a Novel Mini-Lipoxygenase from Nostoc Sphaeroides and Its Application in the Quality Improvement of Steamed Bread

Lipoxygenase (EC1.13.11.12, LOX) has been potentially used in the food industry for food quality improvement. However, the low activity, poor thermal stability, narrow range of pH stability, as well as undesirable isoenzymes and off-flavors, have hampered the application of current commercial LOX. In this study, a putative mini-lipoxygenase gene from cyanobacteria, Nostoc sphaeroides (NsLOX), was cloned and expressed in E. coli BL21. NsLOX displayed only 26.62% structural identity with the reported LOX from Cyanothece sp., indicating it as a novel LOX. The purified NsLOX showed the maximum activity at pH 8.0 and 15 °C, with superior stability at a pH range from 6.0 to 13.0, retaining about 40% activity at 40 °C for 90 min. Notably, NsLOX exhibited the highest specific activity of 78,080 U/mg towards linoleic acid (LA), and the kinetic parameters-K_m, k_cat, and k_cat/K_m-attain values of 19.46 μM, 9199.75 s^-1, and 473.85 μM^-1 s^-1, respectively. Moreover, the activity of NsLOX was obviously activated by Ca²⁺, but it was completely inhibited by Zn²⁺ and Cu²⁺. Finally, NsLOX was supplied in steamed bread and contributed even better improved bread quality than the commercial LOX. These results suggest NsLOX as a promising substitute of current commercial LOX for application in the food industry.

Enhancing the Conformational Stability of the cl-Par-4 Tumor Suppressor via Site-Directed Mutagenesis

Intrinsically disordered proteins play important roles in cell signaling, and dysregulation of these proteins is associated with several diseases. Prostate apoptosis response-4 (Par-4), an approximately 40 kilodalton proapoptotic tumor suppressor, is a predominantly intrinsically disordered protein whose downregulation has been observed in various cancers. The caspase-cleaved fragment of Par-4 (cl-Par-4) is active and plays a role in tumor suppression by inhibiting cell survival pathways. Here, we employed site-directed mutagenesis to create a cl-Par-4 point mutant (D313K). The expressed and purified D313K protein was characterized using biophysical techniques, and the results were compared to that of the wild-type (WT). We have previously demonstrated that WT cl-Par-4 attains a stable, compact, and helical conformation in the presence of a high level of salt at physiological pH. Here, we show that the D313K protein attains a similar conformation as the WT in the presence of salt, but at an approximately two times lower salt concentration. This establishes that the substitution of a basic residue for an acidic residue at position 313 alleviates inter-helical charge repulsion between dimer partners and helps to stabilize the structural conformation.

Modeling the role of charged residues in thermophilic proteins by rotamer and dynamic cross correlation analysis

Discerning the determinants of protein thermostability is very important both from the theoretical and applied perspective. Different lines of evidence seem to indicate that a dynamical network of salt bridges/charged residues plays a fundamental role in the thermostability of enzymes. In this work, we applied measures of dynamic variance, like the Gini coefficients, Kullback-Leibler (KL) divergence and dynamic cross correlation (DCC) coefficients to compare the behavior of 3 pairs of homologous proteins from the thermophilic bacterium Thermus thermophilus and mesophilic Escherichia coli. Molecular dynamic (MD) simulations of these proteins were performed at 303 K and 363 K. In the characterization of their side chain rotamer distributions, the corresponding Gini coefficients and KL-divergence both revealed significant correlations with temperature. Similarly, a DCC analysis revealed a higher trend to de-correlate the movement of charged residues at higher temperatures in the thermophilic proteins, when compared with their mesophilic homologues. These results highlight the importance of dynamic electrostatic network interactions for the thermostability of enzymes.

Genome-wide signatures of adaptation to extreme environments in red algae

The high temperature, acidity, and heavy metal-rich environments associated with hot springs have a major impact on biological processes in resident cells. One group of photosynthetic eukaryotes, the Cyanidiophyceae (Rhodophyta), has successfully thrived in hot springs and associated sites worldwide for more than 1 billion years. Here, we analyze chromosome-level assemblies from three representative Cyanidiophyceae species to study environmental adaptation at the genomic level. We find that subtelomeric gene duplication of functional genes and loss of canonical eukaryotic traits played a major role in environmental adaptation, in addition to horizontal gene transfer events. Shared responses to environmental stress exist in Cyanidiales and Galdieriales, however, most of the adaptive genes (e.g., for arsenic detoxification) evolved independently in these lineages. Our results underline the power of local selection to shape eukaryotic genomes that may face vastly different stresses in adjacent, extreme microhabitats.

Sequence-specific capture and concentration of viral RNA by type III CRISPR system enhances diagnostic

Type-III CRISPR-Cas systems have recently been adopted for sequence-specific detection of SARS-CoV-2. Here, we repurpose the type III-A CRISPR complex from Thermus thermophilus (TtCsm) for programmable capture and concentration of specific RNAs from complex mixtures. The target bound TtCsm complex generates two cyclic oligoadenylates (i.e., cA3 and cA4) that allosterically activate ancillary nucleases. We show that both Can1 and Can2 nucleases cleave single-stranded RNA, single-stranded DNA, and double-stranded DNA in the presence of cA4. We integrate the Can2 nuclease with type III-A RNA capture and concentration for direct detection of SARS-CoV-2 RNA in nasopharyngeal swabs with 15 fM sensitivity. Collectively, this work demonstrates how type-III CRISPR-based RNA capture and concentration simultaneously increases sensitivity, limits time to result, lowers cost of the assay, eliminates solvents used for RNA extraction, and reduces sample handling.

Recent advances in the improvement of cyanobacterial enzymes for bioalkane production

The use of biologically produced alkanes has attracted considerable attention as an alternative energy source to petroleum. In 2010, the alkane synthesis pathway in cyanobacteria was found to include two small globular proteins, acyl-(acyl carrier protein [ACP]) reductase (AAR) and aldehyde deformylating oxygenase (ADO). AAR produces fatty aldehydes from acyl-ACPs/CoAs, which are then converted by ADO to alkanes/alkenes equivalent to diesel oil. This discovery has paved the way for alkane production by genetically modified organisms. Since then, many studies have investigated the reactions catalyzed by AAR and ADO. In this review, we first summarize recent findings on structures and catalytic mechanisms of AAR and ADO. We then outline the mechanism by which AAR and ADO form a complex and efficiently transfer the insoluble aldehyde produced by AAR to ADO. Furthermore, we describe recent advances in protein engineering studies on AAR and ADO to improve the efficiency of alkane production in genetically engineered microorganisms such as Escherichia coli and cyanobacteria. Finally, the role of alkanes in cyanobacteria and future perspectives for bioalkane production using AAR and ADO are discussed. This review provides strategies for improving the production of bioalkanes using AAR and ADO in cyanobacteria for enabling the production of carbon-neutral fuels.

Cryo-Electron Microscopy Snapshots of Eukaryotic Membrane Proteins in Native Lipid-Bilayer Nanodiscs

New technologies for purifying membrane-bound protein complexes in combination with cryo-electron microscopy (EM) have recently allowed the exploration of such complexes under near-native conditions. In particular, polymer-encapsulated nanodiscs enable the study of membrane proteins at high resolution while retaining protein-protein and protein-lipid interactions within a lipid bilayer. However, this powerful technology has not been exploited to address the important question of how endogenous─as opposed to overexpressed─membrane proteins are organized within a lipid environment. In this work, we demonstrate that biochemical enrichment protocols for native membrane-protein complexes from Chaetomium thermophilum in combination with polymer-based lipid-bilayer nanodiscs provide a substantial improvement in the quality of recovered endogenous membrane-protein complexes. Mass spectrometry results revealed ∼1123 proteins, while multiple 2D class averages and two 3D reconstructions from cryo-EM data furnished prominent structural signatures. This integrated methodological approach to enriching endogenous membrane-protein complexes provides unprecedented opportunities for a deeper understanding of eukaryotic membrane proteomes.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.

Analyzing structural features of proteins from deep-sea organisms

Protein adaptations to extreme environmental conditions are drivers in biotechnological process optimization and essential to unravel the molecular limits of life. Most proteins with such desirable adaptations are found in extremophilic organisms inhabiting extreme environments. The deep sea is such an environment and a promising resource that poses multiple extremes on its inhabitants. Conditions like high hydrostatic pressure and high or low temperature are prevalent and many deep-sea organisms tolerate multiple of these extremes. While molecular adaptations to high temperature are comparatively good described, adaptations to other extremes like high pressure are not well-understood yet. To fully unravel the molecular mechanisms of individual adaptations it is probably necessary to disentangle multifactorial adaptations. In this study, we evaluate differences of protein structures from deep-sea organisms and their respective related proteins from nondeep-sea organisms. We created a data collection of 1281 experimental protein structures from 25 deep-sea organisms and paired them with orthologous proteins. We exhaustively evaluate differences between the protein pairs with machine learning and Shapley values to determine characteristic differences in sequence and structure. The results show a reasonable discrimination of deep-sea and nondeep-sea proteins from which we distinguish correlations previously attributed to thermal stability from other signals potentially describing adaptions to high pressure. While some distinct correlations can be observed the overall picture appears intricate.

Essentiality of core hydrophobicity to the structure and function of archaeal chromatin protein Cren7

Studies on the structure-function relationship of protein greatly help to understand not only the principles of protein folding but also the rationales of protein engineering. Crenarchaeal chromatin protein Cren7 provides an excellent research model for this issue. The small protein adopts a 'β-barrel' fold, formed by the double-stranded antiparallel β-sheet 1 tightly packing with the triple-stranded antiparallel β-sheet 2. The simple structure of Cren7 is stabilized by the hydrophobic core between the β-sheets, consisting of the side chains of V8, V10, L20, V25, F41 and F50. In the present work, mutation analyses by alanine substitution of each of the residues in the hydrophobic core were performed. Circular dichroism spectra and nuclear magnetic resonance analyses showed that mutation of F41 led to a significant misfolding of Cren7 through disruption of the β-sheets. Meanwhile, the mutant F41A showed a reduced thermostatility (Tm of 53.2 °C), as compared with the wild-type Cren7 (Tm > 80 °C). Biolayer interferometry and nick-closure assays showed the largely unchanged activities in DNA binding and supercoiling of F41A, indicating the DNA interface of Cren7 was generally retained in F41A. However, F41A was unable to mediate DNA bridging, probably due to the impairment in forming oligomers/polymers on DNA. Atomic force microscopic images of the F41A-DNA complexes also revealed that F41A nearly completely lost the ability to compact DNA into highly condensed structures. Our results not only reveal the critical role of F41 in protein folding of Cren7 but also provide new insights into the structure-function relationships of thermostable proteins.

Crystal structure of thermally stable homodimeric cytochrome c'-β from Thermus thermophilus

Cytochrome c'-β is a heme protein that belongs to the cytochrome P460 family and consists of homodimeric subunits with a predominantly antiparallel β-sheet fold. Here, the crystal structure of cytochrome c'-β from the thermophilic Thermus thermophilus (TTCP-β) is reported at 1.74 Å resolution. TTCP-β has a typical antiparallel β-sheet fold similar to that of cytochrome c'-β from the moderately thermophilic Methylococcus capsulatus (MCCP-β). The phenylalanine cap structure around the distal side of the heme is also similar in TTCP-β and MCCP-β, indicating that both proteins similarly bind nitric oxide and carbon monoxide, as observed spectroscopically. Notably, TTCP-β exhibits a denaturation temperature of 117°C, which is higher than that of MCCP-β. Mutational analysis reveals that the increased homodimeric interface area of TTCP-β contributes to its high thermal stability. Furthermore, 14 proline residues, which are mostly located in the TTCP-β loop regions, possibly contribute to the rigid loop structure compared with MCCP-β, which has only six proline residues. These findings, together with those from phylogenetic analysis, suggest that the structures of Thermus cytochromes c'-β, including TTCP-β, are optimized for function under the high-temperature conditions in which the source organisms live.

Insights on Rigidity and Flexibility at the Global and Local Levels of Protein Structures and Their Roles in Homologous Psychrophilic, Mesophilic, and Thermophilic Proteins: A Computational Study

The rigidity and flexibility of homologous psychrophilic (P), mesophilic (M), and thermophilic (T) proteins have been investigated at the global and local levels in terms of "packing factors" and "atomic fluctuations" obtained from B-factors. For comparison of atomic fluctuations, correction of errors by considering errors in B-factors from all sources in a consolidated manner and conversion of the fluctuations to the same temperature have been suggested and validated. The results indicate no differences in the global values like the average packing factor among the three classes of protein homologues, but at local levels there are differences. A comparison of homologous protein triplets show that the average atomic fluctuations at a given temperature mainly obey the order P > M > T. Packing factors and the atomic fluctuations are anti-correlated, suggesting that altering the rigidity of the active site might be a potential strategy to make tailor-made psychrophilic or thermophilic proteins from their mesophilic homologues. The computer codes developed and used in this work are available at the link https://github.com/Munna-Sarkar/proteins-rigidity-flexibility.git.

Viral protein instability enhances host-range evolvability

Viruses are highly evolvable, but what traits endow this property? The high mutation rates of viruses certainly play a role, but factors that act above the genetic code, like protein thermostability, are also expected to contribute. We studied how the thermostability of a model virus, bacteriophage λ, affects its ability to evolve to use a new receptor, a key evolutionary transition that can cause host-range evolution. Using directed evolution and synthetic biology techniques we generated a library of host-recognition protein variants with altered stabilities and then tested their capacity to evolve to use a new receptor. Variants fell within three stability classes: stable, unstable, and catastrophically unstable. The most evolvable were the two unstable variants, whereas seven of eight stable variants were significantly less evolvable, and the two catastrophically unstable variants could not grow. The slowly evolving stable variants were delayed because they required an additional destabilizing mutation. These results are particularly noteworthy because they contradict a widely supported contention that thermostabilizing mutations enhance evolvability of proteins by increasing mutational robustness. Our work suggests that the relationship between thermostability and evolvability is more complex than previously thought, provides evidence for a new molecular model of host-range expansion evolution, and identifies instability as a potential predictor of viral host-range evolution.

Alteration of Coaxial Heme Ligands Reveals the Role of Heme in Bacterioferritin from Mycobacterium tuberculosis

The uptake and utilization of iron remains critical for the survival/virulence of the host/pathogens in spite of the limitations (low bioavailability/high toxicity) associated with this nutrient. Both the host and pathogens manage to overcome these problems by utilizing the iron repository protein nanocages, ferritins, which not only sequester and detoxify the free Fe(II) ions but also decrease the iron solubility gap by synthesizing/encapsulating the Fe(III)-oxyhydroxide biomineral in its central hollow nanocavity. Bacterial pathogens including Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, encode a distinct subclass of ferritins called bacterioferritin (BfrA), which binds heme, the versatile redox cofactor, via coaxial, conserved methionine (M52) residues at its subunit-dimer interfaces. However, the exact role of heme in Mtb BfrA remains yet to be established. Therefore, its coaxial ligands were altered via site-directed mutagenesis, which resulted in both heme-bound (M52C; ∼1 heme per cage) and heme-free (M52H and M52L) variants, indicating the importance of M52 residues as preferential heme binding axial ligands in Mtb BfrA. All these variants formed intact nanocages of similar size and iron-loading ability as that of wild-type (WT) Mtb BfrA. However, the as-isolated heme-bound variants (WT and M52C) exhibited enhanced protein stability and reductive iron mobilization as compared to their heme-free analogues (M52H and M52L). Further, increasing the heme content in BfrA variants by reconstitution not only enhanced the cage stability but also facilitated the iron mobilization, suggesting the role of heme. In contrary, heme altered the ferroxidase activity to a lesser extent despite facilitating the accumulation of the reactive intermediates formed during the course of the reaction. The current study suggests that heme in Mtb BfrA enhances the overall stability of the protein and possibly acts as an intrinsic electron relay station to influence the iron mineral dissolution and thus may be associated with Mtb's pathogenicity.

Surface residues and non-additive interactions stabilize a consensus homeodomain protein

Despite the widely reported success of consensus design in producing highly stabilized proteins, little is known about the physical mechanisms underlying this stabilization. Here we explore the potential sources of stabilization by performing a systematic analysis of the 29 substitutions that we previously found to collectively stabilize a consensus homeodomain compared to an extant homeodomain. By separately introducing groups of consensus substitutions that alter or preserve charge state, occur at varying degrees of residue burial, and occur at positions of varying degrees of conservation, we determine the extent to which these three features contribute to the consensus stability enhancement. Surprisingly, we find that the largest total contribution to stability comes from consensus substitutions on the protein surface and that the largest per-substitution contributions come from substitutions that maintain charge state. This finding suggests that although consensus proteins are often enriched in charged residues, consensus stabilization does not result primarily from interactions involving charged residues. Although consensus substitutions at strongly conserved positions also contribute disproportionately to stabilization, significant stabilization is also contributed from substitutions at weakly conserved positions. Furthermore, we find that identical consensus substitutions show larger stabilizing effects when introduced into the consensus background than when introduced into an extant homeodomain, indicating that synergistic, stabilizing interactions among the consensus residues contribute to consensus stability enhancement of the homeodomain. By measuring DNA binding affinity for the same set of variants, we find that although consensus design of the homeodomain increases both affinity and folding stability, it does so using a largely non-overlapping set of substitutions.

Abundance and Stability as Common Properties of Allergens

There have been many attempts to identify common biophysical properties which differentiate allergens from their non-immunogenic counterparts. This review will focus on recent studies which examine two such factors: abundance and stability. Anecdotal accounts have speculated that the elevated abundance of potential allergens would increase the likelihood of human exposure and thus the probability of sensitization. Similarly, the stability of potential allergens dictates its ability to remain a viable immunogen during the transfer from the source to humans. This stability could also increase the resilience of potential allergens to both gastric and endosomal degradation, further skewing the immune system toward allergy. Statistical analyses confirm both abundance and stability as common properties of allergens, while epidemiological surveys show a correlation between exposure levels (abundance) and allergic disease. Additional studies show that changes in protein stability can predictably alter gastric/endosomal processing and immunogenicity, providing a mechanistic link between stability and allergenicity. However, notable exceptions exist to both hypotheses which highlight the multifaceted nature of immunological sensitization, and further inform our understanding of some of these other factors and their contribution to allergic disease.

Novel Enzymes From the Red Sea Brine Pools: Current State and Potential

The Red Sea is a marine environment with unique chemical characteristics and physical topographies. Among the various habitats offered by the Red Sea, the deep-sea brine pools are the most extreme in terms of salinity, temperature and metal contents. Nonetheless, the brine pools host rich polyextremophilic bacterial and archaeal communities. These microbial communities are promising sources for various classes of enzymes adapted to harsh environments - extremozymes. Extremozymes are emerging as novel biocatalysts for biotechnological applications due to their ability to perform catalytic reactions under harsh biophysical conditions, such as those used in many industrial processes. In this review, we provide an overview of the extremozymes from different Red Sea brine pools and discuss the overall biotechnological potential of the Red Sea proteome.

Comprehensive mutagenesis to identify amino acid residues contributing to the difference in thermostability between two originally thermostable ancestral proteins

Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.

The anatomy of unfolding of Yfh1 is revealed by site-specific fold stability analysis measured by 2D NMR spectroscopy

Most techniques allow detection of protein unfolding either by following the behaviour of single reporters or as an averaged all-or-none process. We recently added 2D NMR spectroscopy to the well-established techniques able to obtain information on the process of unfolding using resonances of residues in the hydrophobic core of a protein. Here, we questioned whether an analysis of the individual stability curves from each resonance could provide additional site-specific information. We used the Yfh1 protein that has the unique feature to undergo both cold and heat denaturation at temperatures above water freezing at low ionic strength. We show that stability curves inconsistent with the average NMR curve from hydrophobic core residues mainly comprise exposed outliers that do nevertheless provide precious information. By monitoring both cold and heat denaturation of individual residues we gain knowledge on the process of cold denaturation and convincingly demonstrate that the two unfolding processes are intrinsically different.

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.

CD8+ T cells specific for conserved coronavirus epitopes correlate with milder disease in COVID-19 patients

A central feature of the SARS-CoV-2 pandemic is that some individuals become severely ill or die, whereas others have only a mild disease course or are asymptomatic. Here we report development of an improved multimeric αβ T cell staining reagent platform, with each maxi-ferritin "spheromer" displaying 12 peptide-MHC complexes. Spheromers stain specific T cells more efficiently than peptide-MHC tetramers and capture a broader portion of the sequence repertoire for a given peptide-MHC. Analyzing the response in unexposed individuals, we find that T cells recognizing peptides conserved amongst coronaviruses are more abundant and tend to have a "memory" phenotype, compared to those unique to SARS-CoV-2. Significantly, CD8+ T cells with these conserved specificities are much more abundant in COVID-19 patients with mild disease versus those with a more severe illness, suggesting a protective role.

Improving the enzyme property of D-allulose 3-epimerase from a thermophilic organism of Halanaerobium congolense through rational design

The rare sugar d-allulose is an attractive sucrose substitute due to its sweetness and ultra-low caloric value. It can be produced from D-fructose using d-allulose 3-epimerase (DAE) as the biocatalyst. However, most of the reported DAEs show low catalytic efficiency and poor thermostability, which limited their further use in food industrial. Here, a putative d-allulose 3-epimerase from a thermophilic organism of Halanaerobium congolense (HcDAE) was characterized, showing optimal activity at pH 8.0 and 70 °C in the presence of Mg²⁺. Saturation mutagenesis of Y7, C66, and I108, the putative residues responsible for substrate recognition at the O-4, -5, and -6 atoms of D-fructose was performed, and it yielded the triple mutant Y7H/C66L/I108A with improved activity toward D-fructose (345 % of wild-type enzyme). The combined mutant Y7H/C66L/I108A/R156C/K260C exhibited a half-half (t_1/2) of 5.2 h at 70 °C and an increase of the T_m value by 6.5 °C due to the introduction of disulfide bridges between intersubunit with increased interface interactions. The results indicate that mutants could be used as industrial biocatalysts for d-allulose production.

Mini-Intein Structures from Extremophiles Suggest a Strategy for Finding Novel Robust Inteins

Inteins are prevalent among extremophiles. Mini-inteins with robust splicing properties are of particular interest for biotechnological applications due to their small size. However, biochemical and structural characterization has still been limited to a small number of inteins, and only a few serve as widely used tools in protein engineering. We determined the crystal structure of a naturally occurring Pol-II mini-intein from Pyrococcus horikoshii and compared all three mini-inteins found in the genome of P. horikoshii. Despite their similar sizes, the comparison revealed distinct differences in the insertions and deletions, implying specific evolutionary pathways from distinct ancestral origins. Our studies suggest that sporadically distributed mini-inteins might be more promising for further protein engineering applications than highly conserved mini-inteins. Structural investigations of additional inteins could guide the shortest path to finding novel robust mini-inteins suitable for various protein engineering purposes.

A Comparative Study to Decipher the Structural and Dynamics Determinants Underlying the Activity and Thermal Stability of GH-11 Xylanases

With the growing need for renewable sources of energy, the interest for enzymes capable of biomass degradation has been increasing. In this paper, we consider two different xylanases from the GH-11 family: the particularly active GH-11 xylanase from Neocallimastix patriciarum, NpXyn11A, and the hyper-thermostable mutant of the environmentally isolated GH-11 xylanase, EvXyn11^TS. Our aim is to identify the molecular determinants underlying the enhanced capacities of these two enzymes to ultimately graft the abilities of one on the other. Molecular dynamics simulations of the respective free-enzymes and enzyme–xylohexaose complexes were carried out at temperatures of 300, 340, and 500 K. An in-depth analysis of these MD simulations showed how differences in dynamics influence the activity and stability of these two enzymes and allowed us to study and understand in greater depth the molecular and structural basis of these two systems. In light of the results presented in this paper, the thumb region and the larger substrate binding cleft of NpXyn11A seem to play a major role on the activity of this enzyme. Its lower thermal stability may instead be caused by the higher flexibility of certain regions located further from the active site. Regions such as the N-ter, the loops located in the fingers region, the palm loop, and the helix loop seem to be less stable than in the hyper-thermostable EvXyn11^TS. By identifying molecular regions that are critical for the stability of these enzymes, this study allowed us to identify promising targets for engineering GH-11 xylanases. Eventually, we identify NpXyn11A as the ideal host for grafting the thermostabilizing traits of EvXyn11^TS.

General method to stabilize mesophilic proteins in hyperthermal water

The stability of protein structures and biological functions at normal temperature is closely linked with the universal aqueous environment of organisms. Preserving bioactivities of proteins in hyperthermia water would expand their functional capabilities beyond those in native environments. However, only a limited number of proteins derived from hyperthermophiles are thermostable at elevated temperatures. Triggered by this, here we describe a general method to stabilize mesophilic proteins in hyperthermia water. The mesophilic proteins, protected by amphiphilic polymers with multiple binding sites, maintain their secondary and tertiary structures after incubation even in boiling water. This approach, outside the conventional environment for bioactivities of mesophilic proteins, provides a general strategy to dramatically increase the T_m (melting temperature) of mesophilic proteins without any changes to amino sequences of the native proteins. Current work offers a new insight with protein stability engineering for potential application, including vaccine storage and enzyme engineering.

•

Preserving bioactivities of proteins in hyperthermia water is promising.

•

Amphiphilic polymers could protect mesophilic proteins even in boiling water.

•

Mesophilic proteins protected by amphiphilic polymers show dramatically increased T_m.

•

The method offers application prospect for vaccine storage and enzyme engineering.

The concept of protein folding/unfolding and its impacts on human health

Proteins have evolved in specific 3D structures and play different functions in cells and determine various reactions and pathways. The newly synthesized amino acid chains once depart ribosome must crumple into three-dimensional structures so can be biologically active. This process of protein that makes a functional molecule is called protein folding. The protein folding is both a biological and a physicochemical process that depends on the sequence of it. In fact, this process occurs more complicated and in some cases and in exposure to some molecules like glucose (glycation), mistaken folding leads to amyloid structures and fatal disorders called conformational diseases. Such conditions are detected by the quality control system of the cell and these abnormal proteins undergo renovation or degradation. This scenario takes place by the chaperones, chaperonins, and Ubiquitin-proteasome complex. Understanding of protein folding mechanisms from different views including experimental and computational approaches has revealed some intermediate ensembles such as molten globule and has been subjected to biophysical and molecular biology attempts to know more about prevalent conformational diseases.

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.

Structural changes induced by ligand binding drastically increase the thermostability of the Ser/Thr protein kinase TpkD from Thermus thermophilus HB8

Thermophilic proteins maintain their structure at high temperatures through a combination of various factors. Here, we report the ligand-induced stabilization of a thermophilic Ser/Thr protein kinase. Thermus thermophilus TpkD unfolds completely at 55 °C despite the optimum growth temperature of 75 °C. Unexpectedly, we found that the TpkD structure is drastically stabilized by its natural ligands ATP and ADP, as evidenced by the increase in the melting temperature to 80 °C. Such a striking effect of a substrate on thermostability has not been reported for other protein kinases. Conformational changes upon ATP binding were observed in fluorescence quenching and limited proteolysis experiments. Urea denaturation of Trp mutants suggested that ATP binding affects not only the ATP-binding site, but also the remote regions. Our findings shed light on thermoadaptation of thermophilic proteins.

Statistical and molecular dynamics (MD) simulation approach to investigate the role of intrinsically disordered regions of shikimate dehydrogenase in microorganisms surviving at different temperatures

Hyperthermophiles, a subset of prokaryotes that thrive in adverse temperatures, potentially utilize the protein molecular biosystem for maintaining thermostability in a wide range of temperatures. Recent studies revealed that these organisms have smaller proportions of intrinsically disordered proteins. In this study, we performed sequence and structural analysis to investigate the maintenance of protein conformation and their stability at different temperatures. The sequence analysis reveals the higher proportion of charged amino acids are responsible for preventing the helix formation and, hence, become disordered regions. For structural analysis, we chose shikimate dehydrogenase from four species, namely Listeria monocytogenes, Escherichia coli, Thermus thermophilus, and Methanopyrus kandleri, and evaluated the protein adaptation at 283 K, 300 K, and 395 K temperatures. From this investigation, we found more residues of shikimate dehydrogenase prefer an order-to-disorder transition at 395 K only for hyperthermophilic species. The solvent-accessible surface area (SASA) and hydrogen-bond analysis revealed that the tertiary conformation and the number of hydrogen bonds for hyperthermophilic shikimate dehydrogenase are highly preserved at 395 K, compared to 300 K. Our simulation results conjointly provide shikimate dehydrogenase of hyperthermophile which resists high temperatures through stronger protein tertiary conformations.

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.

Hierarchical Conformational Dynamics Confers Thermal Adaptability to preQ1 RNA Riboswitches

Single-stranded noncoding regulatory RNAs, as exemplified by bacterial riboswitches, are highly dynamic. The conformational dynamics allow the riboswitch to reach maximum switching efficiency under appropriate conditions. Here we characterize the conformational dynamics of preQ₁ riboswitches from mesophilic and thermophilic bacterial species at various temperatures. With the integrative use of small-angle X-ray scattering, NMR, and molecular dynamics simulations, we model the ensemble-structures of the preQ₁ riboswitch aptamers without or with a ligand bound. We show that the preQ₁ riboswitch is sufficiently dynamic and fluctuating among multiple folding intermediates only near the physiological temperature of the microorganism. The hierarchical folding dynamics of the RNA involves the docking of 3'-tail to form a second RNA helix and the helical stacking to form an H-type pseudoknot structure. Further, we show that RNA secondary and tertiary dynamics can be modulated by temperature and by the length of an internal loop. The coupled equilibria between RNA folding intermediates are essential for preQ₁ binding, and a four-state exchange model can account for the change of ligand-triggered switching efficiency with temperature. Together, we have established a relationship between the hierarchical dynamics and riboswitch function, and illustrated how the RNA adapts to high temperature.

The use of consensus sequence information to engineer stability and activity in proteins

The goal of protein design is to create proteins that are stable, soluble, and active. Here we focus on one approach to protein design in which sequence information is used to create a "consensus" sequence. Such consensus sequences comprise the most common residue at each position in a multiple sequence alignment (MSA). After describing some general ideas that relate MSA and consensus sequences and presenting a statistical thermodynamic framework that relates consensus and non-consensus sequences to stability, we detail the process of designing a consensus sequence and survey reports of consensus design and characterization from the literature. Many of these consensus proteins retain native biological activities including ligand binding and enzyme activity. Remarkably, in most cases the consensus protein shows significantly higher stability than extant versions of the protein, as measured by thermal or chemical denaturation, consistent with the statistical thermodynamic model. To understand this stability increase, we compare various features of consensus sequences with the extant MSA sequences from which they were derived. Consensus sequences show enrichment in charged residues (most notably glutamate and lysine) and depletion of uncharged polar residues (glutamine, serine, and asparagine). Surprisingly, a survey of stability changes resulting from point substitutions show little correlation with residue frequencies at the corresponding positions within the MSA, suggesting that the high stability of consensus proteins may result from interactions among residue pairs or higher-order clusters. Whatever the source, the large number of reported successes demonstrates that consensus design is a viable route to generating active and in many cases highly stabilized proteins.