Home SNP based analysis depicts phenotypic variability in heme oxygenase-1 protein
Article Open Access

SNP based analysis depicts phenotypic variability in heme oxygenase-1 protein

  • Pratichi Singh ORCID logo EMAIL logo , Syed Habeeb Ahmed , Irfan Ahmad and Mohammad Mahtab Alam
Published/Copyright: February 16, 2023

Abstract

Objectives

The heme oxygenase-1 (HMOX1) gene is a very critical player in cell homeostasis and takes part in heme catabolism. The HMOX1 gene possesses antioxidant, antiapoptotic anti-inflammatory, and antithrombotic properties. This study aimed to identify the deleterious SNPs which may alter the functional and structural attributes of the HMOX1 protein.

Methods

Deleterious SNPs were predicted using Polyphen-1, SIFT Blink, Polyphen-2, I-MUTANT 2.0, PROVEAN, PANTHER, MAPP, SNAP, and PhD-SNP. The 3D structure of the native protein was modelled using ITASSER and validated using PROCHECK. Mutant structures were created through SWISS PDB Viewer. All structures were energy-minimized using GROMACS. Protein-protein interaction (PPI) network was done using STRING.

Results

Three deleterious SNPs (rs146227657, rs373577583, and rs138349040 corresponding to A88D, A131V, and A206E respectively) in the HMOX1 gene were predicted. The structural analysis revealed notable differences in the structural attributes of wild-type and mutant structures. Furthermore, a PPI network was generated for the HMOX1 gene which predicts its interaction with other crucial cellular signaling molecules.

Conclusions

Three reported SNPs (A88D, A131V, and A206E) were identified as highly deleterious in the HMOX1 coding region that can alter the protein secondary structure, protein stability, and its conservation. This paves a new way to study the phenotype of the HMOX1 gene and its correlating diseases.

Introduction

Heme oxygenase-1 (HMOX1 or HO1, OMIM: 141250), also known as HSP32 (32-kDa heat shock protein) is a key rate-limiting enzyme that mediates the degradation of the heme into three products including biliverdin, free ferrous ion, and CO (carbon monoxide). The generated biliverdin is further reduced to bilirubin with the help of biliverdin reductase [1]. These metabolic products participate in several cellular processes that embrace inflammation, proliferation, apoptosis, and oxidative stress [2], [3], [4], [5], [6]. Heme oxygenase (HO) occurs mainly in three isoforms HO-I (Heme oxygenase-1), HO-2 (Heme oxygenase-2), and HO-3 (Heme oxygenase-3). These isoforms are the products of different genes and are differently regulated. HMOX1 is highly inducible which can be induced by various stresses or stimuli, such as heme or hemoglobin, lipopolysaccharide, heat shock, heavy metals, UV radiation, growth factors, phorbol esters, nitric oxide, hydrogen peroxide, hyperoxia, hypoxia, endotoxins, and cytokines [7], [8], [9], [10], [11]. In contrast, HMOX2 is constitutively articulated in most tissues. HMOX3 has the same protein structure as HMOX2 but with lower enzymatic activity than HMOX2 and is not well-characterized [12]. Despite heme catabolic activity, HMOX1 also engages in several disease progressions that include anti-oxidation and immune-modulatory properties and some skin diseases like psoriasis and vitiligo [12]. The deep-routed significance of HMOX1 in several diseases makes its functional investigation of missense mutation a unique method in establishing an enhanced and better therapeutic or diagnostic approach.

In recent times, a plethora of bioinformatics investigations are performed to screen missense mutation inside the protein-coding gene sequence of any gene via sequence or structural based information. The in-silico approach categorizes these SNPs as deleterious or lethal by considering numerous features like physio-chemical properties of the protein, sequence-conservation among several species, and structural attributes [13]. Structural variants in protein occur due to their biochemical properties, such as acidic, basic, hydrophilic or hydrophobic, etc., amino acid substitutions, and the substitution position [13].

In this study, the two diverse approaches including support vector-based and empirical-based methods were employed for the screening of non-synonymous SNPs. The objective of this research was to identify the deleterious or harmful SNPs which may tend to alter the functional and structural attributes of the HMOX1 protein.

Materials and methods

Data mining

The variant genomic region of the HMOX1 gene is harbored from the 1,000 genome project and NCBI database on 11th February 2021. Our research was narrowed down to Homo sapiens. A total of 1954 (HMOX1) non-reductant SNPs were curated. Among the 1954 HMOX1 gene variants, the regional SNPs comprise 65 synonymous SNPs (ss-SNP), 161 missenses non-synonymous (ns-SNP), 84 3′ UTR, 21 5′ UTR, and 1,334 intronic SNPs. The amino acid sequence of the HMOX1 gene (UniProt ID P09601) was recovered from the UniProt database. UniProt database is a collective resource of proteins where all the protein sequence data are accessible. HMOX1 gene contains 288 amino acids sequences downloaded in FASTA format for H. sapiens.

Tools for identifying deleterious SNPs

All the nsSNPs of the HMOX1 gene were exposed to various computational tools like SIFT Blink, PolyPhen-2, PROVEAN, I-MUTANT 2.0, PANTHER, MAPP, SNAP, and PhD-SNP to find its deleterious effect [1415]. The entire nine tools follow different algorithm-based on the homology of the sequences and evolutionary-based conservation methods. The impact of each variation in amino acid over the native residues was prioritized based on combinatorial scoring by the above-mentioned tools.

SIFT blink

SIFT blink (sorting intolerant from tolerant) predicts whether an amino acid substitution may affect the physical and functional properties of the amino acid sequence and ultimately the encoded protein. SIFT evaluates the pathogenicity of a missense mutation based on the degree of conservation of the residues [13, 14, 17]. SIFT is functional to laboratory-induced missense mutations and naturally occurring non-synonymous polymorphisms. SIFT has developed as one of the standard tools for portraying missense variation [16, 17]. SIFT analyzes a single protein using pre-computed BLAST from NCBI. The mutation can be refereed based on the scores; when the probability score is ≥0.05 or <0.05 it is designated as tolerant or deleterious, respectively.

PolyPhen

PolyPhen (Polymorphism Phenotyping v 1&2) is software that forecasts the possible impact of an amino acid substitution on a protein’s structure and its function, PolyPhen score is the indication of possibly damaging substitution. The FASTA format of the native sequence of a protein and the substitution of native and mutant amino acids is given as an input. The output results as score, sensitivity, specificity, and calculate the PSIC (position-specific independent count) and elucidate the effects as benign, probably or possibly damaging [13, 18].

The PolyPhen-2 score ranges between 0.0 and 1.0 (tolerated to deleterious respectively). The variants having a score between 0.0 and 0.15 are considered to be benign, the scores between 0.15 and 1.0 are considered to be damaging and the scores between 0.85 and 1.0 are considered most probably damaging.

PROVEAN

PROVEAN is another tool that works for both SNPs and indels [19]. PROVEAN (Protein Variation Effect Analyzer) envisages whether an amino acid substitution impacts the biological function of any protein [19]. It requires alignment-based scoring metrics to segregate detrimental mutations from deletions, insertions, point mutations, and multiple amino acid substitutions. The resulting score of ≤ 2.5 specifies a disease-causing or damaging variation, and ≥ 2.5 scores indicate tolerated or neutral variants [19].

I-MUTANT

I-MUTANT 2.0 is a sustain vector machine-based web server that performs the automatic guess of protein stability fluctuates upon single-site mutations [20]. I-MUTANT 2.0 can express the stability change of the mutated protein or mutant structure and the stability of the folded protein. I-MUTANT 2.0 allows the assortment of the prediction of protein stability changes at different series of temperatures and pH. FASTA sequences are provided as an input and output protein stability is obtained in the form of energy calculated as DDG in kcal/Mol [20].

PANTHER

PANTHER database defines as Protein Analysis Through Evolutionary Relationship. It uses HMM (hidden markov models) algorithm. The subPSEC scores (Substitution Position Specific Evolutionary Conservation) are calculated by the evolutionary relationship of proteins. Panther requires the protein sequence in FASTA format as input and validates the impact of protein function. The protein sequence with subPSEC score 0 is considered neutral and the score ≤−3 is deemed to be detrimental mutation [21].

MAPP

MAPP defines as a Multivariate analysis of protein polymorphism. Six physicochemical properties are used to evaluate missense variants compromising of charge, polarity, hydropathy, free energy in alpha-helical and beta-sheet conformation, and side-chain volume. The quantification of the physicochemical properties determines the high predictive score. The strong relationship between missense variant and physicochemical characteristics of protein prolongs human disease. It also gives the confidence level of the predicted mutations [22].

SNAP

SNAP defines screening for non-acceptable polymorphisms. This tool can potentially classify the non-synonymous SNPs in all the proteins into neutral and non-neutral (deleterious). It uses the sequence-based algorithm. For each case, SNAP gives a relatability index, which tells the level of confidence of each prediction. SNAP predicts the functional effects of nsSNPs by integrating predicted features of protein structure (including the secondary structure and solvent accessibility), evolutionary information, residue conservation, and other significant information. The sequence of the protein is required as an input. The deleterious mutations with low confidence are reported by SNAP [23].

PhD-SNP

PhD-SNP defines as a Predictor of human Deleterious Single Nucleotide Polymorphisms. It uses a support vector machine (SVM) based algorithm to classify non-synonymous SNPs into human genetic disease-causing (deleterious) or benign mutations [24].

Modeling of HMOX1 protein structure

The presence of SNP can significantly transform protein stability; for better insight on SNPs and their role in protein structural, functional, and stability, it is required to gather information about the 3D structure. The 3D structures allow the comparative analysis of the protein models that checks the structures before and after model refinement and differentiate between the residue geometry and overall structural geometry. The 3D structure of the native protein was modeled using ITASSER [25]. I-TASSER (Iterative Threading Assembly Refinement) is a hierarchical approach to predict protein structure and function. This server provides the most precise structural and functional calculations using state-of-the-art algorithms. After validating the 3D structure using PROCHECK [26] mutants were created through SWISS PDB Viewer [27]. Both native and mutant structures were energy minimized using Gromacs, (GROMOS 43B1 force field) for energy minimization [28].

In this study, the validation of the HMOX1 mutants was done using SAVES (Structure Analysis and Verification Server). This meta-server runs for validating and checking the protein structures throughout and even after model refinement. Moreover, various physical and chemical properties of the HMOX1 protein such as root mean square deviation (RMSD), hydrogen bond, total energy, and conserved residues were calculated using amino acid sequences and its PDB structure.

Protein-protein interaction (PPI) network

STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) is a biological database and a web (universal) resource of known and predicted protein-protein interactions. The STRING database contains information from several sources and computational prediction methods, experimental data, and public text collections. This database provides information of about 9.8 million proteins from >2,000 organisms [29]. STRING has been developed by a consortium of academic institutions, including KU, SIB, CPR, TUD, EMBL, and UZH. STRING imports data resulting from experimentally derived protein interactions through literature review. Connections with proteins in STRING are provided with a confidence score, and with some additional information such as protein domains and their 3D structures, all in a stable and reliable identifier space [29].

Results and discussion

In this study, nsSNPs were selected for HMOX1 using various parameters of coding synonymous (sSNPs) region, coding non-synonymous (nsSNPs) region, mRNA untranslated regions (5′UTR and 3′UTR), and introns. Out of 1954 SNPs, most SNPs were from the intronic region, and there were 161 nsSNPs selected for screening. The functional effect of these 161 sorted nsSNPs was further evaluated by different bioinformatics tools.

Prediction of deleterious nsSNPs

SNP analysis was executed using various bioinformatic tools like SIFT Blink, Polyphen-1, PolyPhen-2, I-MUTANT 2.0, PROVEAN, PANTHER, MAPP, SNAP, and PhD-SNP. Each tool works on a specific algorithm; therefore, by merging different computational tools, prediction of deleterious nsSNPs can be achieved more precisely [13, 14]. Hence, to ease the false positive predictions, a combinatorial approach was determined. Moreover, only those nsSNPs were noted which were generally predicted to be deleterious by all the above-mentioned tools. By comparing the results of SIFT Blink, Polyphen-1, PolyPhen-2, Provean I-MUTANT 2.0, PANTHER, MAPP, SNAP, and PhD-SNP tools out of 161 SNPs, 17 SNPs were depicted deleterious mentioned in Table 1.

Table 1:

List of deleterious nsSNPs computed by various bioinformatics tools.

Rs-id Nucleotide change Mutation PROVEAN I-MUTANT MAPPa PhD-SNPa PolyPhen-1a PolyPhen-2a SIFTa SNAPa PANTHERa Predict-SNPa
rs150399064 A/T K18M Deleterious 0.16 62% 51% 74% 68% 79% 50% 61%
rs77672261 G/T E29D Deleterious −0.51 57% 61% 67% 43% 79% 61% 61%
rs377208959 C/T I65T Deleterious −3.09 66% 58% 67% 47% 79% 61% 57% 63%
rs200856037 C/G N68K Deleterious −0.29 65% 68% 59% 63% 53% 58% 49% 60%
rs369019087 A/G Y78C Deleterious 0.63 75% 68% 74% 63% 79% 62% 72%
rs146227657 A/C A88D Deleterious −1.39 59% 88% 74% 59% 79% 62% 65% 87%
rs373577583 C/T A131V Deleterious 1.1 75% 88% 74% 81% 79% 72% 65% 87%
rs201596112 A/G R136H Deleterious −0.67 62% 88% 59% 81% 79% 81% 87%
rs374687228 A/G G139S Deleterious −1.33 57% 88% 74% 81% 79% 72% 87%
rs368454530 A/G G144D Deleterious −0.76 75% 88% 74% 81% 79% 85% 87%
rs5755713 C/G Q152H Deleterious −0.32 70% 59% 67% 70% 79% 71% 63%
rs200966095 A/T K153I Deleterious 0.67 59% 88% 74% 45% 79% 50% 76%
rs200255845 A/C K153N Deleterious −0.11 64% 68% 67% 47% 45% 67% 60%
rs202031269 A/C K153Q Deleterious −0.07 73% 55% 67% 40% 53% 61% 65%
rs138349040 A/C/T A206E Deleterious −1.05 88% 88% 74% 81% 79% 81% 66% 87%
rs202094347 G/T Q212H Deleterious −1.6 70% 78% 59% 43% 79% 58% 48% 60%
rs35980144 C/T P266L Deleterious −0.95 59% 74% 65% 79% 50% 66% 76%
  1. Commonly deleterious mutations are highlighted in green by all the computational tools. a% gives the expected accuracy, level of confidence [15].

Among these 17 SNPs, three SNPs with rs-id rs146227657, rs373577583, and rs138349040 were reported damaging/harmful by all the nine computational tools mentioned above. The mutation arises due to change in single nucleotide positions at 88, 131, and 206 from A to D, A to V, and A to E respectively [30, 31] highlighted in Table 1. Prior studies have proposed that identifying detrimental ns-SNPs using various algorithmic tools improves prediction accuracy [32]. Therefore, here these reported deleterious mutations (A88D, A131V, and A206E) may hinder the function and expression of HMOX1 protein.

Modeling of mutant HMOX1 protein structures

Protein stability can be reformed by the presence of SNP, for a better understanding of single nucleotide changes and their role in protein’s functional and structural stability, it is obligatory to gather information about the 3D structure. The 3D structure allows the comparative analysis of the protein models that checks the structures before and after model refinement and differentiates between the residue geometry and overall structural geometry. The 3D structure of the native protein was modeled using automated homology modeling via ITASSER and is validated using PROCHECK and Ramachandran plot. The mutants (A88D, A131V, and A206E) were created through SWISS PDB Viewer, native and mutant structures were energy minimized using GROMACS by default for minimization [28]. The single nucleotide polymorphism that occurred in native and mutant structures of HMOX1 protein is revealed in Figure 1. These mutations disrupt the structure and function of the protein. Physiological parameters for instance total energy, root mean square deviation (RMSD), hydrogen bond, Tm-scores were also calculated for these reported substitutions. The mutant A88D has a total energy of −13,271.421 KJ/Mol, RMSD 0.04 Å, and a hydrogen bond distance of 2.18. Mutant A131V has a total energy of −13,490.983 KJ/Mol, RMSD 0.038 Å, and a hydrogen bond distance of 1.98 whereas, mutant A206E has a total energy of −12,420.129 KJ/Mol, RMSD 0.104 Å, and hydrogen bond distance of 2.09. The evolutionary conservation profile of HMOX1 protein of every single amino acid position was also calculated using ConSurf [33] which gives the Tm-score of 1.0000. Supplementary Material Figure 1 demonstrates the conserved residues of HMOX1 protein. Solvent accessibility was also analyzed for HMOX1 protein using SABLE server representing secondary structure shown in Figure 2.

Figure 1: 
Native and mutant structure of HMOX1 protein; (A) Native structure displaying alanine (A) at 88th position; (B) mutant structure displaying aspartic acid at 88th position (protein is in yellow color whereas amino acid residue is in magenta color); (C) native structure displaying Alanine (A) at 131th position; (D) mutant structure displaying valine (V) at 131th position (protein is in green color whereas amino acid residue is in yellow color); (E) native structure displaying alanine (A) at 206th position; (F) mutant structure displaying glutamic acid; (E) at 206th position (protein is in cyan color whereas amino acid residue is in orange color).
Figure 1:

Native and mutant structure of HMOX1 protein; (A) Native structure displaying alanine (A) at 88th position; (B) mutant structure displaying aspartic acid at 88th position (protein is in yellow color whereas amino acid residue is in magenta color); (C) native structure displaying Alanine (A) at 131th position; (D) mutant structure displaying valine (V) at 131th position (protein is in green color whereas amino acid residue is in yellow color); (E) native structure displaying alanine (A) at 206th position; (F) mutant structure displaying glutamic acid; (E) at 206th position (protein is in cyan color whereas amino acid residue is in orange color).

Figure 2: 
Solvent accessibility of every single amino acid of HMOX1 protein.
Figure 2:

Solvent accessibility of every single amino acid of HMOX1 protein.

Network assembly of HMOX1 gene

The protein-protein interaction (PPI) network of the HMOX1 gene was constructed by the STRING database. This includes descriptive statistics such as the number of nodes, expected numbers of edges, average node degree, PPI enrichment p-value, and average clustering coefficient. PPI network of HMOX1 gene in humans was pulled out utilizing predicted functional neighborhood, co-expression, gene-fusion, co-occurrence, experiments, databases, text-mining, homology, and score. The network gained contains 11 nodes and 23 edges, with the average node degree of 4.18, the expected number of edges was 12, and 0.00274 of PPI enrichment p-value as shown in Figure 3.

Figure 3: 
Network assembly of HMOX1 protein.
Figure 3:

Network assembly of HMOX1 protein.

The 10 interacted partners of Heme Oxygenase −1 (HMOX1) were Biliverdin reductase A (BLVRA), Biliverdin reductase B (BLVRB), Ferrochelatase (FECH), Hephaestin (HEPH), Frataxin (FXN), Ceruloplasmin (CP), Nuclear factor erythroid 2-related factor 2 (NFE2L2), NAD(P)H dehydrogenase [quinone] 1 (NQO1), Transcription factor AP-1 (JUN), and Mitogen-activated protein kinase 14 (MAPK14) showing the highest connectivity in the network. These all reported genes are involved in various activities like heme metabolic process, oxidation-reduction process, response to hydrogen peroxide, response to an antibiotic, heme catabolic process, homeostasis, response to a toxic substance, cellular response to oxidative stress, response to an inorganic substance, drug, aging, apoptotis etc. The molecular, biological and cellular functional enrichment of the network is illustrated in Supplementary Material Figure 2.

This report suggests that heme oxygenase (HMOX1) plays a fundamental role in homeostasis and cellular hormetic stress response. The mutation (A88D, A131V, and A206E) in this gene may affect the gene expression and it may concern the pathophysiology of a wide spectrum of human disorders like diabetes, obesity, cardiovascular disease, pulmonary disease, gastrointestinal ailments, dermatitis, kidney dysfunction, cancer, etc., [34]. Consequently, heme oxygenase has been well-thought-out as one of the most projecting targets for therapeutic involvement in the management of various human diseases [35, 36].

Conclusions

In conclusion, three SNPs A88D, A131V, and A206E were identified as highly deleterious in the HMOX1 coding region via multiple computational platforms. These mutations may be of high concern in the HMOX1 gene and its associated diseases as they can alter the protein secondary structure, protein stability, and its conservation. Additionally, the protein-protein interaction network gives an insight into the functional enrichment of HMOX1 and its allied genes. Hence, this study highlights the functional SNPs that affect the phenotypic variability of the protein. Likewise, it paves a way for researchers associated with wet-lab analysis to characterize the impact of each nsSNPs on protein function is difficult, expensive, and time-consuming to design experiments.


Corresponding author: Pratichi Singh, Department of Biosciences, Galgotias University, Greater Noida, Uttar Pradesh, 203201, India, E-mail:

Funding source: Scientific Research Deanship at King Khalid University, Abha, Saudi Arabia

Award Identifier / Grant number: RGP. 01-48-42

Acknowledgments

The authors are grateful to Scientific Research Deanship at King Khalid University, Abha, Saudi Arabia for their financial support through the Small Research Group Project under grant number (RGP. 01-48-42).

  1. Research funding: No specific funding received.

  2. Author contributions: PS designed and conceived the idea for the study. SHA performed the experiments and write the first manuscript draft. PS, IA and MHA analyzed the data and edited the manuscript final version.

  3. Competing interests: The author declares no conflict of interest.

  4. Informed consent: Not applicable.

  5. Ethical approval: Not applicable.

References

1. Yoshida, T, Migita, CT. Mechanism of heme degradation by heme oxygenase. J Inorg Biochem 2000;82:33–41. https://doi.org/10.1016/s0162-0134(00)00156-2.Search in Google Scholar PubMed

2. Chen, S, Wang, X, Nisar, MF, Lin, M, Zhong, JL. Heme oxygenases: cellular multifunctional and protective molecules against UV-Induced oxidative stress. Oxid Med Cell Longev 2019;2019:1–17. https://doi.org/10.1155/2019/5416728. 31885801.Search in Google Scholar PubMed PubMed Central

3. Dennery, PA. Regulation and role of heme oxygenase in oxidative injury. Curr Top Cell Regul 2001;36:181–99.10.1016/S0070-2137(01)80008-XSearch in Google Scholar PubMed

4. Kim, YM, Choi, BM, Kim, YS, Kwon, YG, Kibbe, MR, Billiar, TR, et al.. Protective effect of p53 in vascular smooth muscle cells against nitric oxide-induced apoptosis is mediated by up-regulation of heme oxygenase-2. BMB Rep 2008;41:164–9. https://doi.org/10.5483/bmbrep.2008.41.2.164.Search in Google Scholar PubMed

5. Schipper, HM. Heme oxygenase-1: role in brain aging and neuro degeneration. Exp Gerontol 2000;35:821–30. https://doi.org/10.1016/s0531-5565(00)00148-0.Search in Google Scholar PubMed

6. Zhang, Y, Furuyama, K, Kaneko, K, Ding, Y, Ogawa, K, Yoshizawa, M, et al.. Hypoxia reduces the expression of heme oxygenase-2 in various types of human cell lines: a possible strategy for the maintenance of intracellular heme level. FEBS J 2006;273:3136–47. https://doi.org/10.1111/j.1742-4658.2006.05319.x.Search in Google Scholar PubMed

7. Yoshida, T, Biro, P, Cohen, T, Müller, RM, Shibahara, S. Human heme oxygenase cDNA and induction of its mRNA by hemin. Eur J Biochem 1988;171:457–61. https://doi.org/10.1111/j.1432-1033.1988.tb13811.x.Search in Google Scholar PubMed

8. Alam, J, Shibahara, S, Smith, A. Transcriptional activation of the heme oxygenase gene by heme and cadmium in mouse hepatoma cells. J Biol Chem 1989;264:6371–5. https://doi.org/10.1016/s0021-9258(18)83358-0.Search in Google Scholar

9. Keyse, SM, Tyrrell, RM. Induction of the heme oxygenase gene in human skin fibroblasts by hydrogen peroxide and UVA (365 nm) radiation: evidence for the involvement of the hydroxyl radical. Carcinogenesis 1990;11:787–91. https://doi.org/10.1093/carcin/11.5.787.Search in Google Scholar PubMed

10. Cantoni, L, Rossi, C, Rizzardini, M, Gadina, M, Ghezzi, P. Interleukin-1 and tumour necrosis factor induce hepatic haem oxygenase. Feedback regulation by glucocorticoids. Biochem J 1991;279:891–4. https://doi.org/10.1042/bj2790891.Search in Google Scholar PubMed PubMed Central

11. Rizzardini, M, Carelli, M, Cabello Porras, MR, Cantoni, L. Mechanisms of endotoxin-induced haem oxygenase mRNA accumulation in mouse liver: synergism by glutathione depletion and protection by N-acetylcysteine. Biochem J 1994;304:477–83. https://doi.org/10.1042/bj3040477.Search in Google Scholar PubMed PubMed Central

12. Kishimoto, Y, Kondo, K, Momiyama, Y. The protective role of heme oxygenase-1 in atherosclerotic diseases. Int J Mol Sci 2019;20:3628–43. https://doi.org/10.3390/ijms20153628.Search in Google Scholar PubMed PubMed Central

13. Khan, I, Ansari, IA, Singh, P, Dass, JFP. Prediction of functionally significant single nucleotide polymorphisms in PTEN tumor suppressor gene: an in silico approach. Biotechnol Appl Biochem 2017;64:657–66. https://doi.org/10.1002/bab.1483.Search in Google Scholar PubMed

14. Singh, P, Dass, JF. A multifaceted computational report on the variants effect on KIR2DL3 and IFNL3 candidate gene in HCV clearance. Mol Biol Rep 2016;43:1101–17. https://doi.org/10.1007/s11033-016-4044-5.Search in Google Scholar PubMed

15. Bendl, J, Stourac, J, Salanda, O, Pavelka, A, Wieben, ED, Zendulka, J, et al.. PredictSNP: robust and accurate consensus classifier for prediction of disease-related mutations. PLoS Comput Biol 2014;10:e1003440-1–11. https://doi.org/10.1371/journal.pcbi.1003440. 24453961.Search in Google Scholar PubMed PubMed Central

16. Kumar, P, Henikoff, S, Ng, PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc 2009;4:1073. https://doi.org/10.1038/nprot.2009.86.Search in Google Scholar PubMed

17. Ng, PC, Henikoff, S. Predicting the effects of amino acid substitutions on protein function. Annu Rev Genom Hum Genet 2006;7:61–80. https://doi.org/10.1146/annurev.genom.7.080505.115630.Search in Google Scholar PubMed

18. Adzhubei, IA, Schmidt, S, Peshkin, L, Ramensky, VE, Gerasimova, A, Bork, P, et al.. A method and server for predicting damaging missense mutations. Nat Methods 2010;7:248–9. https://doi.org/10.1038/nmeth0410-248.Search in Google Scholar PubMed PubMed Central

19. Choi, Y, Chan, AP. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 2015;31:2745–7. https://doi.org/10.1093/bioinformatics/btv195.Search in Google Scholar PubMed PubMed Central

20. Capriotti, E, Fariselli, P, Casadio, R. I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res 2005;33(2 Suppl):W306–10. https://doi.org/10.1093/nar/gki375.Search in Google Scholar PubMed PubMed Central

21. Thomas, PD, Campbell, MJ, Kejariwal, A, Mi, H, Karlak, B, Daverman, R, et al.. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res 2003;13:2129–41. https://doi.org/10.1101/gr.772403.Search in Google Scholar PubMed PubMed Central

22. Stone, EA, Sidow, A. Physicochemical constraint violation by missense substitutions mediates impairment of protein function and disease severity. Genome Res 2005;15:978–86. https://doi.org/10.1101/gr.3804205.Search in Google Scholar PubMed PubMed Central

23. Bromberg, Y, Rost, B. SNAP: predict effect of non-synonymous polymorphisms on function. Nucleic Acids Res 2007;35:3823–35. https://doi.org/10.1093/nar/gkm238.Search in Google Scholar PubMed PubMed Central

24. Capriotti, E, Calabrese, R, Casadio, R. Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information. Bioinformatics 2006;22:2729–34. https://doi.org/10.1093/bioinformatics/btl423.Search in Google Scholar PubMed

25. Yang, J, Yan, R, Roy, A, Xu, D, Poisson, J, Zhang, Y. The I-TASSER suite: protein structure and function prediction. Nat Methods 2015;12:7. https://doi.org/10.1038/nmeth.3213.Search in Google Scholar PubMed PubMed Central

26. Laskowski, RA, MacArthur, MW, Moss, DS, Thornton, JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr 1993;26:283–91. https://doi.org/10.1107/s0021889892009944.Search in Google Scholar

27. Johansson, MU, Zoete, V, Michielin, O, Guex, N. Defining and searching for structural motifs using DeepView/Swiss-PdbViewer. BMC Bioinf 2012;13:173. https://doi.org/10.1186/1471-2105-13-173.Search in Google Scholar PubMed PubMed Central

28. Abraham, MJ, Murtola, T, Schulz, R, Páll, S, Smith, JC, Hess, B, et al.. GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015;1:19–25. https://doi.org/10.1016/j.softx.2015.06.001.Search in Google Scholar

29. Szklarczyk, D, Morris, JH, Cook, H, Kuhn, M, Wyder, S, Simonovic, M, et al.. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res 2017;45:362–8. https://doi.org/10.1093/nar/gkw937. 27924014.Search in Google Scholar PubMed PubMed Central

30. Lad, L, Schuller, DJ, Shimizu, H, Friedman, J, Li, H, de Montellano, PR, et al.. Comparison of the heme-free and-bound crystal structures of human heme oxygenase-1. J Biol Chem 2003;278:7834–43. https://doi.org/10.1074/jbc.m211450200.Search in Google Scholar PubMed

31. Tian, S, Yang, X, Zhao, Q, Zheng, J, Huang, H, Chen, Y, et al.. Association between a heme oxygenase-2 genetic variant and risk of Parkinson’s disease in Han Chinese. Neurosci Lett 2017;642:119–22. https://doi.org/10.1016/j.neulet.2017.02.008.Search in Google Scholar PubMed

32. Agúndez, JA, García-Martín, E, Martínez, C, Benito-León, J, Millán-Pascual, J, Díaz-Sánchez, M, et al.. Heme oxygenase-1 and 2 common genetic variants and risk for multiple sclerosis. Sci Rep 2016;6:1–7. https://doi.org/10.1038/srep20830.Search in Google Scholar PubMed PubMed Central

33. Ashkenazy, H, Abadi, S, Martz, E, Chay, O, Mayrose, I, Pupko, T, et al.. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res 2016;44:W344–50. https://doi.org/10.1093/nar/gkw408.Search in Google Scholar PubMed PubMed Central

34. Surh, YJ, Chung, HT, Na, HK, Dulak, J, Stec, DE. Progress in heme oxygenase research. Arch Biochem Biophys 2020;685:108321. https://doi.org/10.1016/j.abb.2020.108321.Search in Google Scholar PubMed PubMed Central

35. Facchinetti, MM. Heme-oxygenase-1. Antioxidants Redox Signal 2020;32:1239–42. https://doi.org/10.1089/ars.2020.8065.Search in Google Scholar PubMed

36. Medina, MV, Sapochnik, D, Garcia Solá, M, Coso, O. Regulation of the expression of heme oxygenase-1: signal transduction, gene promoter activation, and beyond. Antioxidants Redox Signal 2020;32:1033–44. https://doi.org/10.1089/ars.2019.7991.Search in Google Scholar PubMed PubMed Central


Supplementary Material

This article contains supplementary material (https://doi.org/10.1515/tjb-2021-0250).


Received: 2021-10-25
Accepted: 2022-01-07
Published Online: 2023-02-16

© 2023 the author(s), published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 29.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/tjb-2021-0250/html
Scroll to top button