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.

Genomic analysis of Paenibacillus sp. MDMC362 from the Merzouga desert leads to the identification of a potentially thermostable catalase

Microorganisms in hot deserts face heat and other environmental conditions, such as desiccation, UV radiation, or low nutrient availability. Therefore, this hostile environment harbour microorganisms with acquired characteristics related to survival in their habitat, which can be exploited in biotechnology. In this work, the genome of Paenibacillus sp. MDMC362 isolated from the Merzouga desert in Morocco was sequenced to understand its survival strategy's genetic basis; and to evaluate the thermostability of a catalase extracted from genomic annotation files using molecular dynamics. Paenibacillus sp. MDMC362 genome was rich in genetic elements involved in the fight against different stresses, notably temperature stress, UV radiations, osmotic stress, carbon starvation, and oxidative stress. Indeed, we could identify genes of the operons groES-groEL and hrcA-grpE-dnaK and those involved in the different stages of sporulation, which can help the bacteria to survive the high temperatures imposed by a desertic environment. We also observed the genetic components of the UvrABC system and additional mechanisms involved in DNA repair, which help overcome UV radiation damage. Other genes have been identified in the genome, like those coding for ectoine and proline, that aids fight osmotic stress and desiccation. Catalase thermostability investigation using molecular dynamics showed that the protein reached stability and conserved its compactness at temperatures up to 373.15 K. These results suggest a potential thermostability of the enzyme. Since the studied protein is a core protein, thermostability could be conserved among Paenibacillus sp. MDMC362 closely related strains; however, bacteria from harsh environments may have a slight advantage regarding protein stability.

Glycosyl hydrolase 11 (xynA) gene with xylanase activity from thermophilic bacteria isolated from thermal springs

Background

Hemicellulose is one of the copious polymer in lignocellulosic biomass (LCB). It is primarily composed of xylan linked by β-1,4 glycosidic bonds. Xylanase preferentially cleaves the β-1,4-glycosidic bonds in the xylan backbone resulting in complete hydrolysis of the biomass. Thermostable variants of glycoside hydrolases act as robust catalysts, not only in degradation but also during processing, to obtain specific carbohydrate-containing chemicals and materials (Ramasamy et al. in Madras Agric J 107(special):1. 10.29321/MAJ.2020.000382, 2020).

Results

The xylanase production by two thermophilic bacteria isolated from thermal springs was evaluated. In addition, the gene encoding this industrially vital enzyme was isolated and characterized, and its protein structure was analyzed. The thermophilic bacteria producing xylanases were isolated from augmented sawdust and banana fiber biomass from hot springs of Himachal Pradesh and identified as Bacillus subtilis VSDB5 and Bacillus licheniformis KBFB4 using 16S rRNA gene sequencing. The persistent xylanase activity revealed that the enzyme is secreted extracellularly with the maximum activity of 0.76 IU mL−1 and 1.0 IU mL−1 at 6 h and 12 h of growth by KBFB4 and VSDB5, respectively, under submerged fermentation. Both the strains exhibited the maximum activity at pH 6 and a temperature of 50 °C. The xylanases of KBFB4 and VSDB5 were thermostable and retained 40% of their activity at 60 °C after incubation for 30 min. Xylanase of VSDB5 had wide thermotolerance and retained 20% of its activity from 60 to 80 °C, whereas xylanase of KBFB4 showed wide alkali tolerance and retained 80% of its activity until pH 10. The xylanase (xynA)-encoding gene (650 bp) cloned from both the strains using specific primers showed 98 to 99% homology to β-1,4-endoxylanase gene. Further in silico analysis predicted that the xylanase protein, with a molecular weight of 23 kDa, had a high pI (9.44–9.65), which explained the alkaline nature of the enzyme and greater aliphatic index (56.29). This finding suggested that the protein is thermostable. Multiple sequence alignment and homology modeling of the protein sequence revealed that the gene product belonged to the GH11 family, indicating its possible application in bioconversion.

Conclusion

The strains B. subtilis VSDB5 and B. licheniformis KBFB4 obtained from hot springs of Himachal Pradesh produced potent and alkali-tolerant thermostable xylanases, which belong to the GH11 family. The enzyme can be supplemented in industrial applications for biomass conversion at high temperatures and pH (or in processes involving alkali treatment).

Graphical Abstract