MultiClassMetabo: A Superior Classification Model Constructed Using Metabolic Markers in Multiclass Metabolomics.

Multiclass metabolomics has become a popular technique for revealing the mechanisms underlying certain physiological processes, different tumor types, or different therapeutic responses. In multiclass metabolomics, it is highly important to uncover the underlying biological information on biosamples by identifying the metabolic markers with the most associations and classifying the different sample classes. The classification problem of multiclass metabolomics is more difficult than that of the binary problem. To date, various methods exist for constructing classification models and identifying metabolic markers consisting of well-established techniques and newly emerging machine learning algorithms. However, how to construct a superior classification model using these methods remains unclear for a given multiclass metabolomic data set. Herein, MultiClassMetabo has been developed for constructing a superior classification model using metabolic markers identified in multiclass metabolomics. MultiClassMetabo can enable online services, including (a) identifying metabolic markers by marker identification methods, (b) constructing classification models by classification methods, and (c) performing a comprehensive assessment from multiple perspectives to construct a superior classification model for multiclass metabolomics. In summary, MultiClassMetabo is distinguished for its capability to construct a superior classification model using the most appropriate method through a comprehensive assessment, which makes it an important complement to other available tools in multiclass metabolomics. MultiClassMetabo can be accessed at http://idrblab.cn/multiclassmetabo/.

Mechanistic insights into SARS-CoV-2 epidemic via revealing the features of SARS-CoV-2 coding proteins and host responses upon its infection.

SUMMARY: There are seven known coronaviruses that infect humans: four mild coronaviruses, including HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1, only cause mild respiratory diseases, and three severe coronaviruses, including SARS-CoV, MERS-CoV and SARS-CoV-2, can cause severe respiratory diseases even death of infected patients. Both infection and death caused by SARS-CoV-2 are still rapidly increasing worldwide. In this study, we demonstrate that viral coding proteins of SARS-CoV-2 have distinct features and are most, medium and least conserved with SARS-CoV, MERS-CoV and the rest four mild coronaviruses (HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1), respectively. Moreover, expression of host responsive genes (HRG), HRG-enriched biological processes and HRG-enriched KEGG pathways upon infection of SARS-CoV-2 shows slightly overlapping with SARS-CoV and MERS-CoV but distinctive to the four mild coronaviruses. Interestingly, enrichment of overactivation of neutrophil by HRGs is only and commonly found in infections of severe SARS-CoV-2, SARS-CoV and MERS-CoV but not in the other four mild coronaviruses, and the related gene networks show different patterns. Clinical data support that overactivation of neutrophil for severe patients can be one major factor for the similar clinical symptoms observed in SARS-CoV-2 infection compared to infections of the other two severe coronavirus (SARS-CoV and MERS-CoV). Taken together, our study provides a mechanistic insight into SARS-CoV-2 epidemic via revealing the conserved and distinct features of SARS-CoV-2, raising the critical role of dysregulation of neutrophil for SARS-CoV-2 infection. AVAILABILITY AND IMPLEMENTATION: All data sources and analysis methods related to this manuscript are available in the methods, supplementary materials and GEO database (https://www.ncbi.nlm.nih.gov/geo/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Digital gene expression profiling analysis of DNA repair pathways in colon cancer stem population of HT29 cells.

Cancer stem cells (CSCs) contribute to the relapse and development of new neoplasm lesions. While most available clinical approaches, such as chemical and radiation therapies, will kill the majority of cancer cells, they do not kill them all. Some resisting cells, like CSCs, are able to survive due to their excellent self-maintaining capabilities, even in challenging environments. In the present study, we investigated the mRNA level of DNA repair genes of colon CSCs from the HT29 cell line in response to single-strand damage and double-strand breaks, as well as the evident upregulation of key genes in base excision repair, mismatch repair, non-homologous end-joining, and homologous recombination pathways in these cells. Digital gene expression analysis identified upregulated genes in CD44+ HT29 cells that may play important roles in DNA repair. Our results reveal that colon CSCs bear efficient DNA repair abilities, which might explain the survival of colon CSCs after repeated chemical and radiation therapy.

Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K.

To obtain detailed information about the effect of the solvent temperatures on protein dynamics, multiple long molecular dynamics (MD) simulations of serine protease proteinase K with the solute and solvent coupled to different temperatures (either 300 or 180 K) have been performed. Comparative analyses demonstrate that the internal flexibility and mobility of proteinase K are strongly dependent on the solvent temperatures but weakly on the protein temperatures. The constructed free energy landscapes (FELs) at the high solvent temperatures exhibit a more rugged surface, broader spanning range, and higher minimum free energy level than do those at the low solvent temperatures. Comparison between the dynamic hydrogen bond (HB) numbers reveals that the high solvent temperatures intensify the competitive HB interactions between water molecules and protein surface atoms, and this in turn exacerbates the competitive HB interactions between protein internal atoms, thus enhancing the conformational flexibility and facilitating the collective motions of the protein. A refined FEL model was proposed to explain the role of the solvent mobility in facilitating the cascade amplification of microscopic motions of atoms and atomic groups into the global collective motions of the protein.

Insights into Protein-Ligand Interactions: Mechanisms, Models, and Methods.

Molecular recognition, which is the process of biological macromolecules interacting with each other or various small molecules with a high specificity and affinity to form a specific complex, constitutes the basis of all processes in living organisms. Proteins, an important class of biological macromolecules, realize their functions through binding to themselves or other molecules. A detailed understanding of the protein-ligand interactions is therefore central to understanding biology at the molecular level. Moreover, knowledge of the mechanisms responsible for the protein-ligand recognition and binding will also facilitate the discovery, design, and development of drugs. In the present review, first, the physicochemical mechanisms underlying protein-ligand binding, including the binding kinetics, thermodynamic concepts and relationships, and binding driving forces, are introduced and rationalized. Next, three currently existing protein-ligand binding models--the "lock-and-key", "induced fit", and "conformational selection"--are described and their underlying thermodynamic mechanisms are discussed. Finally, the methods available for investigating protein-ligand binding affinity, including experimental and theoretical/computational approaches, are introduced, and their advantages, disadvantages, and challenges are discussed.

A proposed adhesin AoMad1 helps nematode-trapping fungus Arthrobotrys oligospora recognizing host signals for life-style switching.

The nematode-trapping fungus Arthrobotrys oligospora is an important natural enemy of nematodes. It can capture nematodes by producing a special mycelial structure called adhesive network or trap. The trap is also a signature of the fungus switching from the saprophytic lifestyle to the predacious lifestyle. At present, little is known about the mechanism of lifestyle switch in nematode-trapping fungi. Here we describe the effect of a cell wall protein called AoMad1 on lifestyle switch. The disruption of the AoMad1-encoding gene in A. oligospora resulted in the formation of more traps in the presence of nematodes. Interestingly, the mutant strain was more sensitive to certain nitrogen sources as trap inducers than the wild type strain. The microscopic examinations revealed that the AoMad1-deletion mutant lacked cell surface adhesive materials and the cell wall structures were more porous than wild-type strains. A great of genes were differentially expressed by transcriptomic analysis when trap formation was induced by sodium nitrate compared to the wild type strain, many of them were related to nitrogen metabolism, host-pathogen interaction and mycelia development. The results suggest that AoMad1 plays an important role in life style switching in A. oligospora.

The DAF-16/FOXO transcription factor functions as a regulator of epidermal innate immunity.

The Caenorhabditis elegans DAF-16 transcription factor is critical for diverse biological processes, particularly longevity and stress resistance. Disruption of the DAF-2 signaling cascade promotes DAF-16 activation, and confers resistance to killing by pathogenic bacteria, such as Pseudomonas aeruginosa, Staphylococcus aureus, and Enterococcus faecalis. However, daf-16 mutants exhibit similar sensitivity to these bacteria as wild-type animals, suggesting that DAF-16 is not normally activated by these bacterial pathogens. In this report, we demonstrate that DAF-16 can be directly activated by fungal infection and wounding in wild-type animals, which is independent of the DAF-2 pathway. Fungal infection and wounding initiate the Gαq signaling cascade, leading to Ca(2+) release. Ca(2+) mediates the activation of BLI-3, a dual-oxidase, resulting in the production of reactive oxygen species (ROS). ROS then activate DAF-16 through a Ste20-like kinase-1/CST-1. Our results indicate that DAF-16 in the epidermis is required for survival after fungal infection and wounding. Thus, the EGL-30-Ca(2+)-BLI-3-CST-1-DAF-16 signaling represents a previously unknown pathway to regulate epidermal damage response.

TWO INNEXINS OF Spodoptera litura INFLUENCES HEMICHANNEL AND GAP JUNCTION FUNCTIONS IN CELLULAR IMMUNE RESPONSES.

Insect cellular immune responses include encapsulation, nodule formation, and phagocytosis. Hemichannels and gap junctions are involved in these cellular actions. Innexins (Inxs: analogous to the vertebrate connexins) form hemichannels and gap junctions, but the molecular mechanisms underlying their biology is still unclear. In this article, we reported a steady-state level of Inxs (SpliInxs) in hemocytes of Spodoptera litura, which formed nonfunctional hemichannels on the cell surface to maintain normal metabolism. We also reported that two innnexins (SpliInx2 and SpliInx3) were expressed significantly higher in hemocytes compared to other tissues, suggesting that they play important roles in hemocytes. Amino acid analysis found that two cysteine residues in two extracellular loops provided the capability for SpliInx2 and SpliInx3 hemichannels to dock into gap junctions. Western blotting demonstrated that both extracellular and intracellular loops of SpliInx3 and the extracellular loops of SpliInx2 might undergo posttranslational modification during the formation of a steady-state hemichannel. During hemichannel formation, SpliInx2 presented as one isoform, while SpliInx3 presented as three isoforms. These results provide fundamental knowledge for further study of how steady-state levels of SpliInxs are dynamically adjusted to perform cellular immune responses under immune challenge.

Associated bacteria of different life stages of Meloidogyne incognita using pyrosequencing-based analysis.

The root knot nematode (RKN), Meloidogyne incognita, belongs to the most damaging plant pathogens worldwide, and is able to infect almost all cultivated plants, like tomato. Recent research supports the hypothesis that bacteria often associated with plant-parasitic nematodes, function as nematode parasites, symbionts, or commensal organisms etc. In this study, we explored the bacterial consortia associated with M. incognita at different developmental stages, including egg mass, adult female and second-stage juvenile using the pyrosequencing approach. The results showed that Proteobacteria, with a proportion of 71-84%, is the most abundant phylum associated with M. incognita in infected tomato roots, followed by Actinobacteria, Bacteroidetes, Firmicutes etc. Egg mass, female and second-stage juvenile of M. incognita harbored a core microbiome with minor difference in communities and diversities. Several bacteria genera identified in M. incognita are recognized cellulosic microorganisms, pathogenic bacteria, nitrogen-fixing bacteria and antagonists to M. incognita. Some genera previously identified in other plant-parasitic nematodes were also found in tomato RKNs. The potential biological control microorganisms, including the known bacterial pathogens and nematode antagonists, such as Actinomycetes and Pseudomonas, showed the largest diversity and proportion in egg mass, and dramatically decreased in second-stage juvenile and female of M. incognita. This is the first comprehensive report of bacterial flora associated with the RKN identified by pyrosequencing-based analysis. The results provide valuable information for understanding nematode-microbiota interactions and may be helpful in the development of novel nematode-control strategies.

DNA methylation plays important roles in lifestyle transition of Arthrobotrys oligospora.

Trap formation is the key indicator of carnivorous lifestyle transition of nematode-trapping fungi (NTF). Here, the DNA methylation profile was explored during trap induction of Arthrobotrys oligospora, a typical NTF that captures nematodes by developing adhesive networks. Whole-genome bisulfite sequencing identified 871 methylation sites and 1979 differentially methylated regions (DMRs). This first-of-its-kind investigation unveiled the widespread presence of methylation systems in NTF, and suggested potential regulation of ribosomal RNAs through DNA methylation. Functional analysis indicated DNA methylation's involvement in complex gene regulations during trap induction, impacting multiple biological processes like response to stimulus, transporter activity, cell reproduction and molecular function regulator. These findings provide a glimpse into the important roles of DNA methylation in trap induction and offer new insights for understanding the molecular mechanisms driving carnivorous lifestyle transition of NTF.

Genomic and proteomic analyses of the fungus Arthrobotrys oligospora provide insights into nematode-trap formation.

Nematode-trapping fungi are "carnivorous" and attack their hosts using specialized trapping devices. The morphological development of these traps is the key indicator of their switch from saprophytic to predacious lifestyles. Here, the genome of the nematode-trapping fungus Arthrobotrys oligospora Fres. (ATCC24927) was reported. The genome contains 40.07 Mb assembled sequence with 11,479 predicted genes. Comparative analysis showed that A. oligospora shared many more genes with pathogenic fungi than with non-pathogenic fungi. Specifically, compared to several sequenced ascomycete fungi, the A. oligospora genome has a larger number of pathogenicity-related genes in the subtilisin, cellulase, cellobiohydrolase, and pectinesterase gene families. Searching against the pathogen-host interaction gene database identified 398 homologous genes involved in pathogenicity in other fungi. The analysis of repetitive sequences provided evidence for repeat-induced point mutations in A. oligospora. Proteomic and quantitative PCR (qPCR) analyses revealed that 90 genes were significantly up-regulated at the early stage of trap-formation by nematode extracts and most of these genes were involved in translation, amino acid metabolism, carbohydrate metabolism, cell wall and membrane biogenesis. Based on the combined genomic, proteomic and qPCR data, a model for the formation of nematode trapping device in this fungus was proposed. In this model, multiple fungal signal transduction pathways are activated by its nematode prey to further regulate downstream genes associated with diverse cellular processes such as energy metabolism, biosynthesis of the cell wall and adhesive proteins, cell division, glycerol accumulation and peroxisome biogenesis. This study will facilitate the identification of pathogenicity-related genes and provide a broad foundation for understanding the molecular and evolutionary mechanisms underlying fungi-nematodes interactions.

Phylogenetic relationship among East Asian species of the Stegana genus group (Diptera, Drosophilidae).

The phylogenetic relationship among 27 East Asian species of the Stegana genus group was reconstructed using DNA sequences of mitochondrial (COI and ND2) and nuclear (28S) genes. The results lent support to the current generic/subgeneric taxonomic classification in the genus group with the exceptions of the paraphyly of the genus Parastegana and the subgenus Oxyphortica in the genus Stegana. The ancestral areas and divergence times in the genus group were reconstructed/estimated, and accordingly, the biogeographical history of this important clade was discussed. It was proposed that, the evolution of the plant family Fagaceae, especially Quercus, may have played a certain role in facilitating the diversification of the Stegana genus group.

A Deep Learning Approach for Predicting Antigenic Variation of Influenza A H3N2.

Modeling antigenic variation in influenza (flu) virus A H3N2 using amino acid sequences is a promising approach for improving the prediction accuracy of immune efficacy of vaccines and increasing the efficiency of vaccine screening. Antigenic drift and antigenic jump/shift, which arise from the accumulation of mutations with small or moderate effects and from a major, abrupt change with large effects on the surface antigen hemagglutinin (HA), respectively, are two types of antigenic variation that facilitate immune evasion of flu virus A and make it challenging to predict the antigenic properties of new viral strains. Despite considerable progress in modeling antigenic variation based on the amino acid sequences, few studies focus on the deep learning framework which could be most suitable to be applied to this task. Here, we propose a novel deep learning approach that incorporates a convolutional neural network (CNN) and bidirectional long-short-term memory (BLSTM) neural network to predict antigenic variation. In this approach, CNN extracts the complex local contexts of amino acids while the BLSTM neural network captures the long-distance sequence information. When compared to the existing methods, our deep learning approach achieves the overall highest prediction performance on the validation dataset, and more encouragingly, it achieves prediction agreements of 99.20% and 96.46% for the strains in the forthcoming year and in the next two years included in an existing set of chronological amino acid sequences, respectively. These results indicate that our deep learning approach is promising to be applied to antigenic variation prediction of flu virus A H3N2.

Thermus thermophilus proteins that are differentially expressed in response to growth temperature and their implication in thermoadaptation.

As a kind of important extremophiles to realize the adaptation of life at high temperatures, thermophiles have attracted extensive studies. However, the pathways of thermophile proteins related to thermoadaptation remain to be addressed. Our study showed that there existed two types of protein profiles for the thermophile Thermus thermophilus wl in response to temperature change. One of them came from cultures growing below 65 degrees C, which was close to the optimal growth temperature, and another from cultures at or above 65 degrees C. These protein profiles were confirmed by Northern blots. On the basis of the proteomic and computational analyses, it was found that the thermophile proteins related to thermoadaptation might be involved in metabolic pathways as well as the stabilities and modifications of DNA and proteins. Interestingly, for the basic metabolism glycolysis, the phosphoglucomutase was up-regulated at below-optimum temperature, while the glyceraldehyde-3-phosphate dehydrogenase was up-regulated at above-optimum temperature, suggesting that different regulations might be used for basic metabolism at different temperatures. To characterize the proteins in response to high temperatures, superoxide dismutase (SOD), an important enzyme in organism to remove free radical produced in stress environment such as high temperature, was selected as a target protein for this investigation. SOD was inactivated to construct a SOD mutant. The results showed that the SOD protein was essential in thermoadaptation of T. thermophilus. Our study, therefore, presented the thermophile proteins required for thermoadaptation and their possible pathways in thermoadaptation.

Distinguishing between cancer driver and passenger gene alteration candidates via cross-species comparison: a pilot study.

BACKGROUND: We are developing a cross-species comparison strategy to distinguish between cancer driver- and passenger gene alteration candidates, by utilizing the difference in genomic location of orthologous genes between the human and other mammals. As an initial test of this strategy, we conducted a pilot study with human colorectal cancer (CRC) and its mouse model C57BL/6J ApcMin/+, focusing on human 5q22.2 and 18q21.1-q21.2. METHODS: We first performed bioinformatics analysis on the evolution of 5q22.2 and 18q21.1-q21.2 regions. Then, we performed exon-targeted sequencing, real time quantitative polymerase chain reaction (qPCR), and real time quantitative reverse transcriptase PCR (qRT-PCR) analyses on a number of genes of both regions with both human and mouse colon tumors. RESULTS: These two regions (5q22.2 and 18q21.1-q21.2) are frequently deleted in human CRCs and encode genuine colorectal tumor suppressors APC and SMAD4. They also encode genes such as MCC (mutated in colorectal cancer) with their role in CRC etiology unknown. We have discovered that both regions are evolutionarily unstable, resulting in genes that are clustered in each human region being found scattered at several distinct loci in the genome of many other species. For instance, APC and MCC are within 200 kb apart in human 5q22.2 but are 10 Mb apart in the mouse genome. Importantly, our analyses revealed that, while known CRC driver genes APC and SMAD4 were disrupted in both human colorectal tumors and tumors from ApcMin/+ mice, the questionable MCC gene was disrupted in human tumors but appeared to be intact in mouse tumors. CONCLUSIONS: These results indicate that MCC may not actually play any causative role in early colorectal tumorigenesis. We also hypothesize that its disruption in human CRCs is likely a mere result of its close proximity to APC in the human genome. Expanding this pilot study to the entire genome may identify more questionable genes like MCC, facilitating the discovery of new CRC driver gene candidates.

Expansion of Adhesion Genes Drives Pathogenic Adaptation of Nematode-Trapping Fungi.

Understanding how fungi interact with other organisms has significant medical, environmental, and agricultural implications. Nematode-trapping fungi (NTF) can switch to pathogens by producing various trapping devices to capture nematodes. Here we perform comparative genomic analysis of the NTF with four representative trapping devices. Phylogenomic reconstruction of these NTF suggested an evolutionary trend of trapping device simplification in morphology. Interestingly, trapping device simplification was accompanied by expansion of gene families encoding adhesion proteins and their increasing adhesiveness on trap surfaces. Gene expression analysis revealed a consistent up-regulation of the adhesion genes during their lifestyle transition from saprophytic to nematophagous stages. Our results suggest that the expansion of adhesion genes in NTF genomes and consequential increase in trap surface adhesiveness are likely the key drivers of fungal adaptation in trapping nematodes, providing new insights into understanding mechanisms underlying infection and adaptation of pathogenic fungi.

The lifestyle transition of Arthrobotrys oligospora is mediated by microRNA-like RNAs.

The lifestyle transition of fungi, defined as switching from taking organic material as nutrients to pathogens, is a fundamental phenomenon in nature. However, the mechanisms of such transition remain largely unknown. Here we show microRNA-like RNAs (milRNAs) play a key role in fungal lifestyle transition for the first time. We identified milRNAs by small RNA sequencing in Arthrobotrys oligospora, a known nematode-trapping fungus. Among them, 7 highly expressed milRNAs were confirmed by northern-blot analysis. Knocking out two milRNAs significantly decreased A. oligospora's ability to switch lifestyles. We further identified that two of these milRNAs were associated with argonaute protein QDE-2 by RNA-immunoprecipitation (RIP) analysis. Three of the predicted target genes of milRNAs were found in immunoprecipitation (IP) products of QDE-2. Disruption of argonaute gene qde-2 also led to serious defects in lifestyle transition. Interestingly, knocking out individual milRNAs or qde-2 lead to diverse responses under different conditions, and qde-2 itself may be targeted by the milRNAs. Collectively, it indicates the lifestyle transition of fungi is mediated by milRNAs through RNA interference (RNAi) machinery, revealing the wide existence of miRNAs in fungi kingdom and providing new insights into understanding the adaptation of fungi from scavengers to predators and the mechanisms underlying fungal infections.

DA and Xiao-two giant and composite LTR-retrotransposon-like elements identified in the human genome.

We discovered two new complex elements while studying large genomic rearrangements and segmental duplications in the human genome. Both resemble bacterial composite DNA transposon Tn9, consisting of a core flanked by mobile elements, except that the flanking element is not a DNA transposon but instead is long terminal repeat retrotransposon-like with human endogenous retrovirus and satellite sequences. Based on the core size, we named them Xiao ( approximately 30 kb) and DA ( approximately 280 kb), meaning small and big, respectively, in Chinese. Xiao originated from a 19p region encoding olfactory receptor 7E members after the human/ape divergence from Old World monkeys, while DA likely evolved from a Xiao by inserting approximately 200 kb of chimeric sequence from 16p and 21q into the Xiao core, resulting in a target site duplication of 3.4 kb. DA/Xiao was identified in 30 loci on 12 chromosomes, and only DAs mediated intrachromosomal rearrangements, based on our reconstructed human-mouse-rat ancestral genome and the rhesus macaque genome.

The cAMP-PKA pathway-mediated fat mobilization is required for cold tolerance in C. elegans.

Low temperature has a great impact on animal life. Homoiotherms such as mammals increase their energy expenditure to produce heat by activating the cAMP-protein kinase A (PKA)-hormone-sensitive lipase (HSL) pathway under cold stress. Although poikilothermic animals do not have the ability to regulate body temperature, whether this pathway is required for cold tolerance remains unknown. We have now achieved this using the genetically tractable model animal Caenorhabditis elegans. We demonstrate that cold stress activates PKA signaling, which in turn up-regulates the expression of a hormone-sensitive lipase hosl-1. The lipase induces fat mobilization, leading to glycerol accumulation, thereby protecting worms against cold stress. Our findings provide an example of an evolutionarily conserved mechanism for cold tolerance that has persisted in both poikilothermic and homoeothermic animals.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations.

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.