Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases

Metagenomic exploration of cold-active enzymes for detergent applications: Characterization of a novel, cold-active and alkali-stable GH8 endoglucanase from ikaite columns in SW Greenland

Microbial communities from extreme environments are largely understudied, but are essential as producers of metabolites, including enzymes, for industrial processes. As cultivation of most microorganisms remains a challenge, culture-independent approaches for enzyme discovery in the form of metagenomics to analyse the genetic potential of a community are rapidly becoming the way forward. This study focused on analysing a metagenome from the cold and alkaline ikaite columns in Greenland, identifying 282 open reading frames (ORFs) that encoded putative carbohydrate-modifying enzymes with potential applications in, for example detergents and other processes where activity at low temperature and high pH is desired. Seventeen selected ORFs, representing eight enzyme families were synthesized and expressed in two host organisms, Escherichia coli and Aliivibrio wodanis. Aliivibrio wodanis demonstrated expression of a more diverse range of enzyme classes compared to E. coli, emphasizing the importance of alternative expression systems for enzymes from extremophilic microorganisms. To demonstrate the validity of the screening strategy, we chose a recombinantly expressed cellulolytic enzyme from the metagenome for further characterization. The enzyme, Cel240, exhibited close to 40% of its relative activity at low temperatures (4°C) and demonstrated endoglucanase characteristics, with a preference for cellulose substrates. Despite low sequence similarity with known enzymes, computational analysis and structural modelling confirmed its cellulase-family affiliation. Cel240 displayed activity at low temperatures and good stability at 25°C, activity at alkaline pH and increased activity in the presence of CaCl2, making it a promising candidate for detergent and washing industry applications.

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

A Bibliometric Analysis: Current Perspectives and Potential Trends of Enzyme Thermostability from 1991-2022

Thermostability is considered a crucial parameter to evaluate the viability of enzymes in industrial applications. Over the past 31 years, many studies have been reported on the thermostability of enzymes. However, there is no systematic bibliometric analysis of publications on the thermostability of enzymes. In this study, 16,035 publications related to the thermostability of enzymes were searched and collected, showing an increasing annual trend. China contributed the most publications, while the United States had the highest citation count. International Journal of Biological Macromolecules is the most productive journal in the research field. Moreover, Chinese acad sci and Khosro Khajeh are the most active institutions and prolific authors in the field, respectively. Analysis of references with the strongest citation bursts and keyword co-occurrences, magnetic nanoparticles, metal-organic frameworks, molecular dynamics, and rational design are current hot spots and significant future research directions. This study is the first comprehensive bibliometric analysis summarizing trends and developments in enzyme thermostability research. Our findings could provide scholars with an understanding of the fundamental knowledge framework of the field and identify recent potential hotspots and research trends that could facilitate the discovery of collaboration opportunities.

Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

Charge reversal mutations in mesophilic-thermophilic orthologous protein pairs and their role in enhancing coulombic interaction energy

Proteins from thermophilic organisms are a matter of immense interest for decades because of its application in fields like de-novo protein design, thermostable variants of biocatalysts etc. Previous studies have found several sequence and structural adaptations related to thermal stability, while charge reversal study remains ignored. Here we address whether charge reversal mutations naturally occur in mesophilic-thermophilic/hyperthermophilic orthologous proteins. Do they contribute to thermal stability? Our systematic study on 1550 mesophilic-thermophilic/hyperthermophilic orthologous protein pairs with remarkable structural and topological similarity, shows gain in coulombic interaction energy in thermophilic/hyperthermophilic proteins at short range associated with partially exposed and buried charge reversal mutations, which may enhance thermostability. Our findings call forth its application in future protein engineering studies. Communicated by Ramaswamy H. Sarma.

Versatility of subtilisin: A review on structure, characteristics, and applications

Due to its thermostability and high pH compatibility, subtilisin is most known for its role as an additive for detergents in which it is categorized as a serine protease according to MEROPS database. Subtilisin is typically isolated from various bacterial species of the Bacillus genus such as Bacillus subtilis, B. amyloliquefaciens, B. licheniformis, and various other organisms. It is composed of 268-275 amino acid residues and is initially secreted in the precursor form, preprosubtilisin, which is composed of 29-residues signal peptide, 77-residues propeptide, and 275-residues active subtilisin. Subtilisin is known for the presence of high and low affinity calcium binding sites in its structure. Native subtilisin has general properties of thermostability, tolerance to neutral to high pH, broad specificity, and calcium-dependent stability, which contribute to the versatility of subtilisin applicability. Through protein engineering and immobilization technologies, many variants of subtilisin have been generated, which increase the applicability of subtilisin in various industries including detergent, food processing and packaging, synthesis of inhibitory peptides, therapeutic, and waste management applications.

The First Comprehensive Biodiversity Study of Culturable Fungal Communities Inhabiting Cryoconite Holes in the Werenskiold Glacier on Spitsbergen (Svalbard Archipelago, Arctic)

Cryoconite holes on glacier surfaces are a source of cold-adapted microorganisms, but little is known about their fungal inhabitants. Here, we provide the first report of distinctive fungal communities in cryoconite holes in the Werenskiold Glacier on Spitsbergen (Svalbard Archipelago, Arctic). Due to a combination of two incubation temperatures (7 °C and 24 ± 0.5 °C) and two media during isolation (PDA, YPG), as well as classical and molecular identification approaches, we were able to identify 23 different fungi (21 species and 2 unassigned species). Most of the fungi cultured from cryoconite sediment were ascomycetous filamentous micromycetes. However, four representatives of macromycetes were also identified (Bjerkandera adusta, Holwaya mucida, Orbiliaceae sp., and Trametes versicolor). Some of the described fungi possess biotechnological potential (Aspergillus pseudoglaucus, A. sydowii, Penicillium expansum, P. velutinum, B. adusta, and T. versicolor), thus, we propose the Arctic region as a source of new strains for industrial applications. In addition, two phytopathogenic representatives were present (P. sumatraense, Botrytis cinerea), as well as one potentially harmful to humans (Cladosporium cladosporioides). To the best of our knowledge, we are the first to report the occurrence of A. pseudoglaucus, C. allicinum, C. ramotenellum, P. sumatraense, P. velutinum, P. cumulodentata, B. adusta, and T. versicolor in polar regions. In all likelihood, two unassigned fungus species (Orbiliaceae and Dothideomycetes spp.) might also be newly described in such environments. Additionally, due to experimenting with 10 sampling sites located at different latitudes, we were able to conclude that the number of fungal spores decreases as one moves down the glacier. Considering the prevalence and endangerment of glacial environments worldwide, such findings suggest their potential as reservoirs of fungal diversity, which should not be overlooked.

Thermostability engineering of industrial enzymes through structure modification

Thermostability is an essential requirement of enzymes in the industrial processes to catalyze the reactions at high temperatures; thus, enzyme engineering through directed evolution, semi-rational design and rational design are commonly employed to construct desired thermostable mutants. Several strategies are implemented to fulfill enzymes' thermostability demand including decreasing the entropy of the unfolded state through substitutions Gly → Xxx or Xxx → Pro, hydrogen bond, salt bridge, introducing two different simultaneous interactions through single mutant, hydrophobic interaction, filling the hydrophobic cavity core, decreasing surface hydrophobicity, truncating loop, aromatic-aromatic interaction and introducing positively charged residues to enzyme surface. In the current review, horizons about compatibility between secondary structures and substitutions at preferable structural positions to generate the most desirable thermostability in industrial enzymes are broadened. KEY POINTS: • Protein engineering is a powerful tool for generating thermostable industrial enzymes. • Directed evolution and rational design are practical approaches in enzyme engineering. • Substitutions in preferable structural positions can increase thermostability.

Identification of lactoferrin-derived peptides as potential inhibitors against the main protease of SARS-CoV-2

COVID-19 is a global health emergency that causes serious concerns. A global effort is underway to identify drugs for the treatment of COVID-19. One possible solution to the present problem is to develop drugs that can inhibit SARS-CoV-2 main protease (M^pro), a coronavirus protein that been considered as one among many drug targets. In this work, lactoferrin from Bos taurus L. was in silico hydrolyzed. The bioactivity, water solubility, and ADMET properties of the generated peptides were predicted using various online tools. The molecular interactions between M^pro and the peptides were studied using molecular docking and molecular dynamic simulation. The results demonstrated that peptide GSRY was predicted to have better physicochemical properties, and the value of '-C DOCKER interaction energy' between peptide GSRY and M^pro was 80.8505 kcal/mol. The interaction between the peptide GSRY and the native ligand N3 co-crystallized with M^pro had overlapped amino acids, i.e., HIS163, GlY143, GLU166, GLN189 and MET165. Molecular dynamic simulation revealed that M^pro/GSRY complexes were stable. Collectively, the peptide GSRY may be a potential candidate drug against M^pro of SARS-CoV-2.

Experimental diagnostic of sequence-variant dynamic perturbations revealed by broadband dielectric spectroscopy

Genetic diversity leads to protein robustness, adaptability, and failure. Some sequence variants are structurally robust but functionally disturbed because mutations bring the protein onto unfolding/refolding routes resulting in misfolding diseases (e.g., Parkinson). We assume dynamic perturbations introduced by mutations foster the alternative unfolding routes and test this possibility by comparing the unfolding dynamics of the heat-labile enterotoxin B pentamers and the cholera toxin B pentamers, two pentamers structurally and functionally related and robust to 17 sequence variations. The B-subunit thermal unfolding dynamics are monitored by broadband dielectric spectroscopy in nanoconfined and weakly hydrated conditions. Distinct dielectric signals reveal the different B-subunits unfolding dynamics. Combined with network analyses, the experiments pinpoint the role of three mutations A1T, E7D, and E102A, in diverting LTB₅ to alternative unfolding routes that protect LTB₅ from dissociation. Altogether, the methodology diagnoses dynamics faults that may underlie functional disorder, drug resistance, or higher virulence of sequence variants.

Exploring the Cold-Adaptation Mechanism of Serine Hydroxymethyltransferase by Comparative Molecular Dynamics Simulations

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to μs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.

Increased surface charge in the protein chaperone Spy enhances its anti-aggregation activity

Chaperones are essential components of the protein homeostasis network. There is a growing interest in optimizing chaperone function, but exactly how to achieve this aim is unclear. Here, using a model chaperone, the bacterial protein Spy, we demonstrate that substitutions that alter the electrostatic potential of Spy's concave, client-binding surface enhance Spy's anti-aggregation activity. We show that this strategy is more efficient than one that enhances the hydrophobicity of Spy's surface. Our findings thus challenge the traditional notion that hydrophobic interactions are the major driving forces that guide chaperone-substrate binding. Kinetic data revealed that both charge- and hydrophobicity-enhanced Spy variants release clients more slowly, resulting in a greater "holdase" activity. However, increasing short-range hydrophobic interactions deleteriously affected Spy's ability to capture substrates, thus reducing its in vitro chaperone activity toward fast-aggregating substrates. Our strategy in chaperone surface engineering therefore sought to fine-tune the different molecular forces involved in chaperone-substrate interactions rather than focusing on enhancing hydrophobic interactions. These results improve our understanding of the mechanistic basis of chaperone-client interactions and illustrate how protein surface-based mutational strategies can facilitate the rational improvement of molecular chaperones.

The Energetic Origin of Different Catalytic Activities in Temperature-Adapted Trypsins

Psychrophilic enzymes were always observed to have higher catalytic activity (k cat) than their mesophilic homologs at room temperature, while the origin of this phenomenon remains obscure. Here, we used two different temperature-adapted trypsins, the psychrophilic Atlantic cod trypsin (ACT) and the mesophilic bovine trypsin (BT), as a model system to explore the energetic origin of their different catalytic activities using computational methods. The results reproduce the characteristic changing trends in the activation free energy, activation enthalpy, and activation entropy between the psychrophilic and mesophilic enzymes, where, in particular, the slightly decreased activation free energy of ACT is determined by its considerably reduced activation enthalpy rather than by its more negative activation entropy compared to BT. The calculated electrostatic contributions to the solvation free energies in the reactant state/ground sate (RS/GS) and transition state (TS) show that, going from BT to ACT, the TS stabilization has a predominant effect over the RS stabilization on lowering the activation enthalpy of ACT. Comparison between the solvation energy components reveals a more optimized electrostatic preorganization to the TS in ACT, which provides a larger stabilization to the TS through reducing the reorganization energy, thus resulting in the lower activation enthalpy and hence lower activation free energy of ACT. Thus, it can be concluded that it is the difference in the protein electrostatic environment, and hence its different stabilizing effects on the TS, that brings about the different catalytic activities of different temperature-adapted trypsins.

New Insight into Mechanisms of Protein Adaptation to High Temperatures: A Comparative Molecular Dynamics Simulation Study of Thermophilic and Mesophilic Subtilisin-Like Serine Proteases

In high-temperature environments, thermophilic proteins must possess enhanced thermal stability in order to maintain their normal biological functions. However, the physicochemical basis of the structural stability of thermophilic proteins at high temperatures remains elusive. In this study, we performed comparative molecular dynamics simulations on thermophilic serine protease (THM) and its homologous mesophilic counterpart (PRK). The comparative analyses of dynamic structural and geometrical properties suggested that THM adopted a more compact conformation and exhibited more intramolecular interactions and lower global flexibility than PRK, which could be in favor of its thermal stability in high-temperature environments. Comparison between protein solvent interactions and the hydrophobicity of these two forms of serine proteases showed that THM had more burial of nonpolar areas, and less protein solvent hydrogen bonds (HBs), indicating that solvent entropy maximization and mobility may play a significant role in THM's adaption to high temperature environments. The constructed funnel-like free energy landscape (FEL) revealed that, in comparison to PRK, THM had a relatively flat and narrow free energy surface, and a lower minimum free energy level, suggesting that the thermophilic form had lower conformational diversity and flexibility. Combining the FEL theory and our simulation results, we conclude that the solvent (entropy force) plays a significant role in protein adaption at high temperatures.

Molecular dynamics simulations reveal distinct differences in conformational dynamics and thermodynamics between the unliganded and CD4-bound states of HIV-1 gp120

The entry of human immunodeficiency virus type I (HIV-1) into host cells is initiated by binding to the cell-surface receptor CD4, which induces a conformational transition of the envelope (Env) glycoprotein gp120 from the closed, unliganded state to the open, CD4-bound state. Despite many available structures in these two states, detailed aspects on the dynamics and thermodynamics of gp120 remain elusive. Here, we performed microsecond-scale (μs-scale) multiple-replica molecular dynamics (MD) simulations to explore the differences in the conformational dynamics, protein motions, and thermodynamics between the unliganded and CD4-bound/complexed forms of gp120. Comparative analyses of MD trajectories reveal that CD4 binding promotes the structural deviations/changes and conformational flexibility, loosens the structural packing, and complicates the molecular motions of gp120. Comparison of the constructed free energy landscapes (FELs) reveals that the CD4-complexed gp120 has more conformational substates, larger conformational entropy, and lower thermostability than the unliganded form. Therefore, the unliganded conformation represents a structurally and energetically stable "ground state" for the full-length gp120. The observed great increase in the mobility of V1/V2 and V3 along with their more versatile movement directions in the CD4-bound gp120 compared to the unliganded form suggests that their orientations with respect to each other and to the structural core determine the differences in the conformational dynamics and thermodynamics between the two gp120 forms. The results presented here provide a basis by which to better understand the functional and immunological properties of gp120 and, furthermore, to deploy appropriate strategies for the development of anti-HIV-1 drugs or vaccines.

Experimental Protein Molecular Dynamics: Broadband Dielectric Spectroscopy coupled with nanoconfinement

Protein dynamics covers multiple spatiotemporal scale processes, among which slow motions, not much understood even though they are underlying protein folding and protein functions. Protein slow motions are associated with structural heterogeneity, short-lived and poorly populated conformations, hard to detect individually. In addition, they involve collective motions of many atoms, not easily tracked by simulation and experimental devices. Here we propose a biophysical approach, coupling geometrical nanoconfinement and broadband dielectric spectroscopy (BDS), which distinguishes protein conformations by their respective molecular dynamics. In particular, protein-unfolding intermediates, usually poorly populated in macroscopic solutions are detected. The protein dynamics is observed under unusual conditions (sample nanoconfinement and dehydration) highlighting the robustness of protein structure and protein dynamics to a variety of conditions consistent with protein sustainability. The protein dielectric signals evolve with the temperature of thermal treatments indicating sensitivity to atomic and molecular interaction changes triggered by the protein thermal unfolding. As dipole fluctuations depend on both collective large-scale motions and local motions, the approach offers a prospect to track in-depth unfolding events.

Effect of Temperature and pH on the Secondary Structure and Denaturation Process of Jumbo Squid Hepatopancreas Cathepsin D

BACKGROUND

Cathepsin D is a lysosomal enzyme that is found in all organisms acting in protein turnover, in humans it is present in some types of carcinomas, and it has a high activity in Parkinson's disease and a low activity in Alzheimer disease. In marine organisms, most of the research has been limited to corroborate the presence of this enzyme. It is known that cathepsin D of some marine organisms has a low thermostability and that it has the ability to have activity at very acidic pH. Cathepsin D of the Jumbo squid (Dosidicus gigas) hepatopancreas was purified and partially characterized. The secondary structure of these enzymes is highly conserved so the role of temperature and pH in the secondary structure and in protein denaturation is of great importance in the study of enzymes. The secondary structure of cathepsin D from jumbo squid hepatopancreas was determined by means of circular dichroism spectroscopy.

OBJECTIVE

In this article, our purpose was to determine the secondary structure of the enzyme and how it is affected by subjecting it to different temperature and pH conditions.

METHODS

Circular dichroism technique was used to measure the modifications of the secondary structure of cathepsin D when subjected to different treatments. The methodology consisted in dissecting the hepatopancreas of squid and freeze drying it. Then a crude extract was prepared by mixing 1: 1 hepatopancreas with assay buffer, the purification was in two steps; the first step consisted of using an ultrafiltration membrane with a molecular cut of 50 kDa, and the second step, a pepstatin agarose resin was used to purification the enzyme. Once the enzyme was purified, the purity was corroborated with SDS PAGE electrophoresis, isoelectric point and zymogram. Circular dichroism is carried out by placing the sample with a concentration of 0.125 mg / mL in a 3 mL quartz cell. The results were obtained in mdeg (millidegrees) and transformed to mean ellipticity per residue, using 111 g/mol molecular weight/residue as average. Secondary-structure estimation from the far-UV CD spectra was calculated using K2D Dichroweb software.

RESULTS

It was found that α helix decreases at temperatures above 50 °C and above pH 4. Heating the enzyme above 70°C maintains a low percentage of α helix and increases β sheet. Far-UV CD measurements of cathepsin D showed irreversible thermal denaturation. The process was strongly dependent on the heating rate, accompanied by a process of oligomerization of the protein that appears when the sample is heated, and maintained a certain time at this temperature. An amount typically between 3 and 4% α helix of their secondary structure remains unchanged. It is consistent with an unfolding process kinetically controlled due to the presence of an irreversible reaction. The secondary structure depends on pH, and a pH above 4 causes α helix structures to be modified.

CONCLUSION

In conclusion, cathepsin D from jumbo squid hepatopancreas showed retaining up to 4% α helix at 80°C. The thermal denaturation of cathepsin D at pH 3.5 is under kinetic control and follows an irreversible model.