# Prediction of the Impact of Deleterious Nonsynonymous Single Nucleotide Polymorphisms on the HumanRRM2B Gene: A Molecular Modeling Study

Authors: Chaimaa Ait El Cadi; Al Mehdi Krami; Hicham Charoute; Zouhair Elkarhat; Najat Sifeddine; Hamid Lakhiari; Hassan Rouba; Abdelhamid Barakat; Halima Nahili
Journal: BioMed Research International
 (2020)

---

## Abstract

RRM2B gene encodes ribonucleoside-diphosphate reductase subunit M2 B, the p53-inducible small subunit (p53R2) of ribonucleotide reductase (RNR), an enzyme catalyzing dNTP synthesis for mitochondrial DNA. Defects in this gene may cause severe mitochondrial disease affecting mainly the nervous system. This study is aimed at examining the effect of deleterious nonsynonymous SNP (nsSNP) on the structure of the RRM2B protein, using a variety of prediction tools followed by a molecular modeling analysis. After using 13 algorithms, 19 nsSNPs were predicted deleterious. Among these variants, 18 decreased the protein stability and 16 were localized in very highly conserved regions. Protein 3D structure analysis showed that 18 variants changed amino acid interactions. These results concur with what has been found in experimental trials; 7 deleterious nsSNPs were previously reported in patients suffering from genetic disorders affecting the nervous system. Thus, our study will provide useful information to design more efficient and fast genetic tests to find RRM2B gene mutations.


---

## Body

## 1. Introduction

Mitochondria are semiautonomous, self-reproducing organelles that occur in the cytoplasm of all cells of most eukaryotes. In that manner, a mitochondrial function can be provided by either mitochondrial or nuclear DNA; therefore, any disruption of both genetic materials can lead to severe mitochondrial diseases [1]. Human ribonucleotide reductase (RNR) is a heterotetramer enzyme that catalyzes the synthesis of dNTPs by direct reduction of ribonucleoside diphosphates to deoxyribonucleoside diphosphates required for DNA replication; it consists of two subunits: a large catalytic subunit called R1 and a small subunit called R2 [2, 3]. Actually, cells have two forms of the R2 subunit, a form regulated by the cell cycle expressed abundantly at the S phase and a second form regulated by the p53 tumor suppressor protein called p53R2.This form is necessary for the synthesis and repair of mitochondrial DNA (mtDNA) inside nonproliferative cells [3, 4]. Indeed, this form is encoded by a gene known as ribonucleoside diphosphate reductase subunit M2 B or RRM2B located on chromosome 8 in position 8q22.3 [5].Mutations within theRRM2B gene result in a variety of inherited mitochondrial diseases either recessive or dominant and can be divided into 2 major types: diseases characterized by mtDNA depletion and diseases characterized by multiple deletions of the mtDNA. mtDNA depletion, where inheritance is autosomal recessive, is generally characterized by severe multisystemic manifestations such as encephalomyopathy with proximal renal tubulopathy usually fatal in early age [5–10]. However, multiple mtDNA deletions cause tissue-specific cytochrome c oxidase (COX) deficiency. Inheritance can be either autosomal recessive with progressive external ophthalmoplegia (PEO) and multisystem involvement manifesting during early childhood/adulthood or autosomal dominant with less severe manifestations appearing during late adulthood [11–15]. We can also mention Kearns-Sayre syndrome (KSS) as one of the rare phenotypes that may occur due to mutations in the RRM2B gene. RRM2B gene variants can be found within all exons, but exon 9 is specifically considered as a mutation hotspot, since variations in this region are responsible for a truncated protein presumed to cause a dominant-negative or gain-of-function effect on the heterotetrameric structure of the RNR enzyme [11, 12].There is currently a large number ofin silico tools used to predict the structural impact caused by amino acid (aa) changes in a given molecule. These tools are not very accurate; however, they can still be used as an initial filter of potentially deleterious changes. Subsequently, we can first resort to molecular modeling which allows us to observe the impact of these changes on the protein 3D structure. Then, advanced analyses such as molecular dynamics simulation can be used to better assess the impact caused by these variants. However, the use of such a method requires a large computing power.In this study, we will investigate the impact of nonsynonymous single nucleotide polymorphisms (nsSNPs) on theRRM2B protein, followed by molecular modeling analysis. This will allow us to evaluate the effect of potentially deleterious variations on the protein structure.


## 2. Data and Methods

### 2.1. Data and Sequence Extraction

The sequence of theRRM2B protein has been extracted in Fasta format from UniProt (ID: Q7LG56) (https://www.uniprot.org/uniprot/Q7LG56). Meanwhile, the SNPs were generated from NCBI’s dbSNP database (https://www.ncbi.nlm.nih.gov/snp/).


### 2.2. Most Damaging Variant Identification

In order to have an idea about the impact of nsSNPs on the structure of theRRM2B protein, we used the following algorithms: SIFT [16], POLYPHEN [17], SNAP [18], PhD-SNP [19], Condel [20], PROVEAN [21], M-Cap [22], LRT [23], Mutation Assessor [24], Mutation Taster [25], PredictSNP [26], MAPP [27], and PANTHER [28]. Similarly to other studies [29], we maintained the standard cutoff or the threshold scores of all prediction tools to select the most deleterious amino acid changes on the protein function.


### 2.3. Stability and Conservation Analysis

I-Mutant (http://folding.biofold.org/i-mutant/) is a stability prediction software that measures the degree of protein destabilization and gives the predicted free energy change value DDG, which is the Gibbs free energy value from the mutated protein minus the Gibbs free energy value from the wild type expressed in Kcal/mol. The software judges an SNP “Increasing” the stability of the protein if the DDG value is superior to 0 or “Decreasing” it if the DDG value is below 0.Mupro (http://mupro.proteomics.ics.uci.edu/) is also another stability prediction web server that uses support vector machines and neural networks, to predict the impact of a single site amino acid mutation on protein stability. This tool provides a prediction reliability score between -1 and 1.Conservation analysis was conducted by ConSurf (http://consurf.tau.ac.il/), a web server identifying the functional regions of a protein by evaluating the degree of its conservation.We also attempted to examine the exact number of conserved residues upon theRRM2B protein sequence, using the Structurally Conserved Region (SCR) prediction web server (http://prodata.swmed.edu/scr_prediction/index.php).


### 2.4. Molecular Modeling of Native and Variant Forms of the RRM2B Protein

For the purpose of having a clear picture of the effect of each variant on the structure ofRRM2B, we need the three-dimensional structure (3D) of native protein and its different mutated forms. This was done by homology modeling using the SWISS-MODEL server [30, 31].


### 2.5. Sequence Visualization and Interaction Change Observation

After obtaining the PDB files of the 3D models, YASARA software [32] was used to visualize and evaluate possible changes in the interactions between amino acids comparing to the native model, using the various features and tools available.


## 2.1. Data and Sequence Extraction

The sequence of theRRM2B protein has been extracted in Fasta format from UniProt (ID: Q7LG56) (https://www.uniprot.org/uniprot/Q7LG56). Meanwhile, the SNPs were generated from NCBI’s dbSNP database (https://www.ncbi.nlm.nih.gov/snp/).


## 2.2. Most Damaging Variant Identification

In order to have an idea about the impact of nsSNPs on the structure of theRRM2B protein, we used the following algorithms: SIFT [16], POLYPHEN [17], SNAP [18], PhD-SNP [19], Condel [20], PROVEAN [21], M-Cap [22], LRT [23], Mutation Assessor [24], Mutation Taster [25], PredictSNP [26], MAPP [27], and PANTHER [28]. Similarly to other studies [29], we maintained the standard cutoff or the threshold scores of all prediction tools to select the most deleterious amino acid changes on the protein function.


## 2.3. Stability and Conservation Analysis

I-Mutant (http://folding.biofold.org/i-mutant/) is a stability prediction software that measures the degree of protein destabilization and gives the predicted free energy change value DDG, which is the Gibbs free energy value from the mutated protein minus the Gibbs free energy value from the wild type expressed in Kcal/mol. The software judges an SNP “Increasing” the stability of the protein if the DDG value is superior to 0 or “Decreasing” it if the DDG value is below 0.Mupro (http://mupro.proteomics.ics.uci.edu/) is also another stability prediction web server that uses support vector machines and neural networks, to predict the impact of a single site amino acid mutation on protein stability. This tool provides a prediction reliability score between -1 and 1.Conservation analysis was conducted by ConSurf (http://consurf.tau.ac.il/), a web server identifying the functional regions of a protein by evaluating the degree of its conservation.We also attempted to examine the exact number of conserved residues upon theRRM2B protein sequence, using the Structurally Conserved Region (SCR) prediction web server (http://prodata.swmed.edu/scr_prediction/index.php).


## 2.4. Molecular Modeling of Native and Variant Forms of the RRM2B Protein

For the purpose of having a clear picture of the effect of each variant on the structure ofRRM2B, we need the three-dimensional structure (3D) of native protein and its different mutated forms. This was done by homology modeling using the SWISS-MODEL server [30, 31].


## 2.5. Sequence Visualization and Interaction Change Observation

After obtaining the PDB files of the 3D models, YASARA software [32] was used to visualize and evaluate possible changes in the interactions between amino acids comparing to the native model, using the various features and tools available.


## 3. Results

### 3.1. SNP’s Distribution

During this study, we managed to extract a total of 8464 SNPs of theRRM2B gene using the dbSNP database of NCBI. 6621 of these SNPs are in the intron region, 720 in the 3′ UTR region, 552 upstream, 184 nonsynonymous, 125 in the 5′ UTR region, 114 downstream, 74 synonymous, and 74 others (Figure 1). This study focuses only on nonsynonymous SNPs.Figure 1
Distribution of SNPs present in theRRM2B gene.


### 3.2. Deleterious nsSNPs of the RRM2B Gene

The results of the 13 prediction algorithms showed that 19 among 184 nsSNPs were deleterious (Table1). For the majority of prediction software, 1 is the higher score and 0 is the lower score, except for the SIFT tool where the higher score is 0 and the lower score is 1.Table 1
Software prediction and scores for the 19 deleterious nsSNP of theRRM2B gene. For most prediction software: D: deleterious. Mutation Assessor: H: high. Mutation Taster: D: disease causing and A: disease causing automatic.
dbSNP IDSNPSIFTPOLYPHENELRTM-CAPPhD-SNPSNAPMutation AssessorMutation tasterPROVEANCondelPredictSNPMAPPPANTHERrs515726181R41WD (0)D (0.999)DDD (0.875)D (0.805)HDDD (0.935)D (0.869)D (0.426)D (0.780)rs200273673R41PD (0)D (0.997)DDD (0.885)D (0.848)HDDD (0.911)D (0.869)D (0.560)D (0.718)rs267607025R110PD (0)D (0.992)DDD (0.858)D (0.622)HDDD (0.892)D (0.869)D (0.588)D (0.718)rs515726187F123SD (0)D (1)DDD (0.817)D (0.622)HDDD (0.945)D (0.869)D (0.409)D (0.714)rs1211197364Y124CD (0)D (1)DDD (0.858)D (0.805)HDDD (0.945)D (0.869)D (0.588)—rs1185904532E131AD (0)D (1)DDD (0.875)D (0.848)HDDD (0.945)D (0.869)D (0.633)D (0.661)rs515726188E131KD (0,01)D (0.999)DDD (0.885)D (0.848)HDDD (0.896)D (0.869)D (0.483)D (0.483)rs1313813193S139RD (0)D (1)DDD (0.858)D (0.848)HDDD (0.945)D (0.869)D (0.571)—rs1169870960I142TD (0,02)D (0.992)DDD (0.608)D (0.555)HDDD (0.831)D (0.869)D (0.657)D (0.661)rs515726190R186GD (0)D (0.997)DDD (0.875)D (0.720)HDDD (0.911)D (0.869)D (0.573)—rs515726191E194GD (0)D (0.994)DDD (0.817)D (0.805)HDDD (0.897)D (0.869)D (0.571)—rs121918308E194KD (0)D (0.996)DDD (0.875)D (0.848)HADD (0.906)D (0.869)D (0.462)—rs515726192G195RD (0,02)D (0.998)DDD (0.875)D (0.805)HDDD (0.857)D (0.869)D (0.633)D (0.742)rs863224192G200ED (0)D (0.998)DDD (0.885)D (0.848)HDDD (0.919)D (0.869)D (0.877)—rs1283028277F206VD (0)D (0.912)DDD (0.817)D (0.6221)HDDD (0.808)D (0.869)D (0.462)D (0.687)rs575109470G216RD (0)D (1)DDD (0.885)D (0.720)HDDD (0.945)D (0.869)D (0.761)—rs1422185855D227VD (0)D (1)DDD (0.875)D (0.805)HDDD (0.945)D (0.869)D (0.775)—rs772868983H231QD (0)D (0.998)DDD (0.858)D (0.848)HDDD (0.919)D (0.869)D (0.509)D (0.687)rs1215473912A235TD (0)D (0.96)DDD (0.676)D (0.555)HDDD (0.848)D (0.869)D (0.633)D (0.611)


### 3.3. Stability Analysis

To investigate the effect of each of the 19 nsSNPs on the stability of theRRM2B protein, two web tools were used: I-Mutant and Mupro. According to the DDG values, only 1 nsSNP (Y124C) increased the stability of the RRM2B protein according to I-Mutant prediction, while for Mupro, all the 19 nsSNPs decreased the protein stability (Table 2).Table 2
Stability software predictions and DDG scores for the 19 deleterious nsSNPs of theRRM2B gene.
dbSNP IDSNPI-MutantMuproPredictionDDG valuePredictionDDG valuers515726181R41WDecrease-0.57Decrease-1.809rs200273673R41PDecrease-1.06Decrease-2.19rs267607025R110PDecrease-1.47Decrease-1.312rs515726187F123SDecrease-2.49Decrease-1.688rs1211197364Y124CIncrease0.56Decrease-0.935rs1185904532E131ADecrease-0.04Decrease-1.24rs515726188E131KDecrease-0.63Decrease-1.644rs1313813193S139RDecrease-0.75Decrease-0.849rs1169870960I142TDecrease-2.42Decrease-3.633rs515726190R186GDecrease-1.91Decrease-1.33rs515726191E194GDecrease-2.61Decrease-1.121rs121918308E194KDecrease-1.5Decrease-0.523rs515726192G195RDecrease-0.82Decrease-0.393rs863224192G200EDecrease-0.94Decrease-0.399rs1283028277F206VDecrease-2.2Decrease-1.157rs575109470G216RDecrease-1.08Decrease-0.503rs1422185855D227VDecrease-0.92Decrease-0.507rs772868983H231QDecrease-0.62Decrease-0.199rs1215473912A235TDecrease-1.63Decrease-1.261


### 3.4. Conservation Analysis

According to the conservation scale color code, we can see that theRRM2B protein sequence contains a high proportion of conserved regions, which make this protein more vulnerable for structural property alterations. These observations were confirmed through the examination of the exact number of conserved residues upon the RRM2B protein sequence, using the SCR prediction web server. The results showed that the RRM2B protein contains a total of 225 (out of 351) conserved residue representing 64% of the whole protein sequence (Table 3).Table 3
Number of conserved residues in the RRM2B protein provided by the SCR prediction server.
Sequence length351Structural conservation cutoff0.74Residues predicted to be conserved225Using the ConSurf server, all the 19 nsSNP residues are predicted to be located in conserved regions of theRRM2B protein, representing exposed or buried and/or structural or functional residues. It is to mention that a structural residue is an amino acid that provides a structural framework to the protein, while a functional residue is an amino acid that mediates interactions of the protein with other biomolecules.According to the conservation results, 16 nsSNPs were found to be localized in very highly conserved regions (9 on the conservation scale); 9 residues were predicted to be both exposed and functional, whereas 7 amino acids were buried and structural.Two buried nsSNPs were located in a conserved region (8 on the conservation scale). The remaining SNP was found in a moderately conserved region (7 on the conservation scale) (Figure2).Figure 2
Conservation results of theRRM2B protein provided by the ConSurf web server.


### 3.5. Molecular Modeling

The SWISS-MODEL server generated 3D structures of the wild-type and mutated proteins based on the 4DJN template. This template was the perfect candidate for the homology modeling, thanks to its high resolution of 2.2 Å and other important variables (Qmean, coverage and sequence identity). Protein PDB files were downloaded from this web server and visualized using the YASARA software (Figure3).Figure 3
The RRM2B protein structure predicted using the SWISS-MODEL server.


### 3.6. Structural Analysis

After visualizing the structure of the wild-type and mutated proteins using YASARA, we have noticed several differences in the interactions between amino acids (hydrogen or hydrophobic bonds).As illustrated in Figure4(a), the Arg amino acid in position 41 of the native form of RRM2B protein that has a hydrophobic bond with Gln 48 and hydrogen bonds with Ile 44 and Glu 119 lost interactions after being replaced by Pro (Figure 4(b)) or by Trp (Figure 4(c)).Figure 4
Comparison of the native RRM2B protein structure and two mutant forms. (a) The structural model of the wild-type protein (Arg 41). (b) The structural model of the first mutated form (Pro 41). (c) The structural model of the second mutated form (Trp 41). Discontinuous cylinders represent hydrogen bonds; continuous lines represent hydrophobic bonds. Red residues are the main residues where the nsSNPs appeared; purple residues are those who that have a hydrogen bond with the main residue; green residues are those who that have a hydrophobic bond with the main residue; magenta residues mark the loss of a bond between an amino acid and the main residue in the variant form that existed in the wild-type form.
(a)(b)(c)In the wild-type protein, the Phe 123 residue kept his 2 hydrogen bonds with Gln 127 and Glu 119 (Figure5(a)) and his hydrophobic bond with Phe 234 when replaced by Ser (Figure 5(b)). Amino acid interaction analysis for all studied nsSNPs is summarized in Table 4, and their corresponding illustrations can be found in the supplementary materials.Figure 5
Comparison of the native RRM2B protein structure and a mutant form. (a) The structural model of the wild-type protein (Phe 123). (b) The structural model of the mutated form (Ser 123).
(a)(b)Table 4
Analysis of nsSNP effect on hydrogen bonds and hydrophobic interactions.
dbSNP IDSNPHydrogen bondsHydrophobic bondsWild typeVariantWild typeVariantrs515726181R41PIle 44 and Glu 119 (2)—Gln 48—rs200273673R41WIle 44 and Glu 119 (2)—Gln 48—rs267607025R110PGlu 114, Asp 178, Ile 176, and Asn 106Glu 114Glu 109Ile 176 and Phe 111rs515726187F123SGlu 119 and Gln 127Glu 119 and Gln 127Phe 234Phe 234rs1211197364Y124CIle 128 and Ala 120Ile 128 and Ala 120Ser 112, Leu 107, Phe 111, and Phe 234Ser 112rs1185904532E131AGln 127 (2) and Ser 135Gln 127 and Ser 135Gln 127—rs515726188E131KGln 127 (2) and Ser 135Gln 127 and Ser 135Gln 127Asp 100, Val 103, and His 231rs1313813193S139RSer 135, Asp 143, and Phe 45Ser 135 and Asp 143Phe 45Phe 45 and Ile 44rs1169870960I142TTyr 138 and Ile 146Tyr 138 and Ile 146Phe 156, Ala 97, and Leu 93Phe 156rs515726190R186GGlu 114 (2), Phe 190, Thr 182, and Ser 181Phe 190 and Thr 182Phe 111—rs515726191E194GPhe 190Phe 190—Phe 198rs121918308E194KPhe 190Phe 190 (2)—His 231rs515726192G195RAla 191Ala 191 and Cys 232—Val 289rs863224192G200EAla 204Ala 204—Ile 286rs1283028277F206VPhe 202 and Lys 210Phe 202 and Lys 210Phe 311, Met 312, and Glu 222Phe 311rs575109470G216RAsp 74 and Ser 220 (2)Asp74 and Ser 220 (2)Leu 71—rs1422185855D227VLeu 223, Arg 226 (2), Trp 64, and His 231Leu 223 and His 231—His 231rs772868983H231QAsp 227 and Ala 235Asp 227 and Ala 235Gln 127—rs1215473912A235THis 231 and Phe 239His 231 and Phe 239Ala 191 and Val 187Ala 191(2): two interactions with the same residue.Mutated proteins were superimposed with the wild type to calculate the root-mean-square deviation (RMSD). For all structures, the RMSD (Å) values did not deviate significantly from the native protein (Table5).Table 5
Root-mean-square deviation (RMSD) of deleterious nsSNPs.
dbSNP IDSNPRMSD (Å)rs515726181R41W0.0274rs200273673R41P0.0236rs267607025R110P0.0446rs515726187F123S0.0032rs1211197364Y124C0.0019rs1185904532E131A0.0165rs515726188E131K0.0275rs1313813193S139R0.0964rs1169870960I142T0.0025rs515726190R186G0.0114rs515726191E194G0.0160rs121918308E194K0.0239rs515726192G195R0.0378rs863224192G200E0.0537rs1283028277F206V0.0059rs575109470G216R0.0684rs1422185855D227V0.0101rs772868983H231Q0.0181rs1215473912A235T0.0251


## 3.1. SNP’s Distribution

During this study, we managed to extract a total of 8464 SNPs of theRRM2B gene using the dbSNP database of NCBI. 6621 of these SNPs are in the intron region, 720 in the 3′ UTR region, 552 upstream, 184 nonsynonymous, 125 in the 5′ UTR region, 114 downstream, 74 synonymous, and 74 others (Figure 1). This study focuses only on nonsynonymous SNPs.Figure 1
Distribution of SNPs present in theRRM2B gene.


## 3.2. Deleterious nsSNPs of the RRM2B Gene

The results of the 13 prediction algorithms showed that 19 among 184 nsSNPs were deleterious (Table1). For the majority of prediction software, 1 is the higher score and 0 is the lower score, except for the SIFT tool where the higher score is 0 and the lower score is 1.Table 1
Software prediction and scores for the 19 deleterious nsSNP of theRRM2B gene. For most prediction software: D: deleterious. Mutation Assessor: H: high. Mutation Taster: D: disease causing and A: disease causing automatic.
dbSNP IDSNPSIFTPOLYPHENELRTM-CAPPhD-SNPSNAPMutation AssessorMutation tasterPROVEANCondelPredictSNPMAPPPANTHERrs515726181R41WD (0)D (0.999)DDD (0.875)D (0.805)HDDD (0.935)D (0.869)D (0.426)D (0.780)rs200273673R41PD (0)D (0.997)DDD (0.885)D (0.848)HDDD (0.911)D (0.869)D (0.560)D (0.718)rs267607025R110PD (0)D (0.992)DDD (0.858)D (0.622)HDDD (0.892)D (0.869)D (0.588)D (0.718)rs515726187F123SD (0)D (1)DDD (0.817)D (0.622)HDDD (0.945)D (0.869)D (0.409)D (0.714)rs1211197364Y124CD (0)D (1)DDD (0.858)D (0.805)HDDD (0.945)D (0.869)D (0.588)—rs1185904532E131AD (0)D (1)DDD (0.875)D (0.848)HDDD (0.945)D (0.869)D (0.633)D (0.661)rs515726188E131KD (0,01)D (0.999)DDD (0.885)D (0.848)HDDD (0.896)D (0.869)D (0.483)D (0.483)rs1313813193S139RD (0)D (1)DDD (0.858)D (0.848)HDDD (0.945)D (0.869)D (0.571)—rs1169870960I142TD (0,02)D (0.992)DDD (0.608)D (0.555)HDDD (0.831)D (0.869)D (0.657)D (0.661)rs515726190R186GD (0)D (0.997)DDD (0.875)D (0.720)HDDD (0.911)D (0.869)D (0.573)—rs515726191E194GD (0)D (0.994)DDD (0.817)D (0.805)HDDD (0.897)D (0.869)D (0.571)—rs121918308E194KD (0)D (0.996)DDD (0.875)D (0.848)HADD (0.906)D (0.869)D (0.462)—rs515726192G195RD (0,02)D (0.998)DDD (0.875)D (0.805)HDDD (0.857)D (0.869)D (0.633)D (0.742)rs863224192G200ED (0)D (0.998)DDD (0.885)D (0.848)HDDD (0.919)D (0.869)D (0.877)—rs1283028277F206VD (0)D (0.912)DDD (0.817)D (0.6221)HDDD (0.808)D (0.869)D (0.462)D (0.687)rs575109470G216RD (0)D (1)DDD (0.885)D (0.720)HDDD (0.945)D (0.869)D (0.761)—rs1422185855D227VD (0)D (1)DDD (0.875)D (0.805)HDDD (0.945)D (0.869)D (0.775)—rs772868983H231QD (0)D (0.998)DDD (0.858)D (0.848)HDDD (0.919)D (0.869)D (0.509)D (0.687)rs1215473912A235TD (0)D (0.96)DDD (0.676)D (0.555)HDDD (0.848)D (0.869)D (0.633)D (0.611)


## 3.3. Stability Analysis

To investigate the effect of each of the 19 nsSNPs on the stability of theRRM2B protein, two web tools were used: I-Mutant and Mupro. According to the DDG values, only 1 nsSNP (Y124C) increased the stability of the RRM2B protein according to I-Mutant prediction, while for Mupro, all the 19 nsSNPs decreased the protein stability (Table 2).Table 2
Stability software predictions and DDG scores for the 19 deleterious nsSNPs of theRRM2B gene.
dbSNP IDSNPI-MutantMuproPredictionDDG valuePredictionDDG valuers515726181R41WDecrease-0.57Decrease-1.809rs200273673R41PDecrease-1.06Decrease-2.19rs267607025R110PDecrease-1.47Decrease-1.312rs515726187F123SDecrease-2.49Decrease-1.688rs1211197364Y124CIncrease0.56Decrease-0.935rs1185904532E131ADecrease-0.04Decrease-1.24rs515726188E131KDecrease-0.63Decrease-1.644rs1313813193S139RDecrease-0.75Decrease-0.849rs1169870960I142TDecrease-2.42Decrease-3.633rs515726190R186GDecrease-1.91Decrease-1.33rs515726191E194GDecrease-2.61Decrease-1.121rs121918308E194KDecrease-1.5Decrease-0.523rs515726192G195RDecrease-0.82Decrease-0.393rs863224192G200EDecrease-0.94Decrease-0.399rs1283028277F206VDecrease-2.2Decrease-1.157rs575109470G216RDecrease-1.08Decrease-0.503rs1422185855D227VDecrease-0.92Decrease-0.507rs772868983H231QDecrease-0.62Decrease-0.199rs1215473912A235TDecrease-1.63Decrease-1.261


## 3.4. Conservation Analysis

According to the conservation scale color code, we can see that theRRM2B protein sequence contains a high proportion of conserved regions, which make this protein more vulnerable for structural property alterations. These observations were confirmed through the examination of the exact number of conserved residues upon the RRM2B protein sequence, using the SCR prediction web server. The results showed that the RRM2B protein contains a total of 225 (out of 351) conserved residue representing 64% of the whole protein sequence (Table 3).Table 3
Number of conserved residues in the RRM2B protein provided by the SCR prediction server.
Sequence length351Structural conservation cutoff0.74Residues predicted to be conserved225Using the ConSurf server, all the 19 nsSNP residues are predicted to be located in conserved regions of theRRM2B protein, representing exposed or buried and/or structural or functional residues. It is to mention that a structural residue is an amino acid that provides a structural framework to the protein, while a functional residue is an amino acid that mediates interactions of the protein with other biomolecules.According to the conservation results, 16 nsSNPs were found to be localized in very highly conserved regions (9 on the conservation scale); 9 residues were predicted to be both exposed and functional, whereas 7 amino acids were buried and structural.Two buried nsSNPs were located in a conserved region (8 on the conservation scale). The remaining SNP was found in a moderately conserved region (7 on the conservation scale) (Figure2).Figure 2
Conservation results of theRRM2B protein provided by the ConSurf web server.


## 3.5. Molecular Modeling

The SWISS-MODEL server generated 3D structures of the wild-type and mutated proteins based on the 4DJN template. This template was the perfect candidate for the homology modeling, thanks to its high resolution of 2.2 Å and other important variables (Qmean, coverage and sequence identity). Protein PDB files were downloaded from this web server and visualized using the YASARA software (Figure3).Figure 3
The RRM2B protein structure predicted using the SWISS-MODEL server.


## 3.6. Structural Analysis

After visualizing the structure of the wild-type and mutated proteins using YASARA, we have noticed several differences in the interactions between amino acids (hydrogen or hydrophobic bonds).As illustrated in Figure4(a), the Arg amino acid in position 41 of the native form of RRM2B protein that has a hydrophobic bond with Gln 48 and hydrogen bonds with Ile 44 and Glu 119 lost interactions after being replaced by Pro (Figure 4(b)) or by Trp (Figure 4(c)).Figure 4
Comparison of the native RRM2B protein structure and two mutant forms. (a) The structural model of the wild-type protein (Arg 41). (b) The structural model of the first mutated form (Pro 41). (c) The structural model of the second mutated form (Trp 41). Discontinuous cylinders represent hydrogen bonds; continuous lines represent hydrophobic bonds. Red residues are the main residues where the nsSNPs appeared; purple residues are those who that have a hydrogen bond with the main residue; green residues are those who that have a hydrophobic bond with the main residue; magenta residues mark the loss of a bond between an amino acid and the main residue in the variant form that existed in the wild-type form.
(a)(b)(c)In the wild-type protein, the Phe 123 residue kept his 2 hydrogen bonds with Gln 127 and Glu 119 (Figure5(a)) and his hydrophobic bond with Phe 234 when replaced by Ser (Figure 5(b)). Amino acid interaction analysis for all studied nsSNPs is summarized in Table 4, and their corresponding illustrations can be found in the supplementary materials.Figure 5
Comparison of the native RRM2B protein structure and a mutant form. (a) The structural model of the wild-type protein (Phe 123). (b) The structural model of the mutated form (Ser 123).
(a)(b)Table 4
Analysis of nsSNP effect on hydrogen bonds and hydrophobic interactions.
dbSNP IDSNPHydrogen bondsHydrophobic bondsWild typeVariantWild typeVariantrs515726181R41PIle 44 and Glu 119 (2)—Gln 48—rs200273673R41WIle 44 and Glu 119 (2)—Gln 48—rs267607025R110PGlu 114, Asp 178, Ile 176, and Asn 106Glu 114Glu 109Ile 176 and Phe 111rs515726187F123SGlu 119 and Gln 127Glu 119 and Gln 127Phe 234Phe 234rs1211197364Y124CIle 128 and Ala 120Ile 128 and Ala 120Ser 112, Leu 107, Phe 111, and Phe 234Ser 112rs1185904532E131AGln 127 (2) and Ser 135Gln 127 and Ser 135Gln 127—rs515726188E131KGln 127 (2) and Ser 135Gln 127 and Ser 135Gln 127Asp 100, Val 103, and His 231rs1313813193S139RSer 135, Asp 143, and Phe 45Ser 135 and Asp 143Phe 45Phe 45 and Ile 44rs1169870960I142TTyr 138 and Ile 146Tyr 138 and Ile 146Phe 156, Ala 97, and Leu 93Phe 156rs515726190R186GGlu 114 (2), Phe 190, Thr 182, and Ser 181Phe 190 and Thr 182Phe 111—rs515726191E194GPhe 190Phe 190—Phe 198rs121918308E194KPhe 190Phe 190 (2)—His 231rs515726192G195RAla 191Ala 191 and Cys 232—Val 289rs863224192G200EAla 204Ala 204—Ile 286rs1283028277F206VPhe 202 and Lys 210Phe 202 and Lys 210Phe 311, Met 312, and Glu 222Phe 311rs575109470G216RAsp 74 and Ser 220 (2)Asp74 and Ser 220 (2)Leu 71—rs1422185855D227VLeu 223, Arg 226 (2), Trp 64, and His 231Leu 223 and His 231—His 231rs772868983H231QAsp 227 and Ala 235Asp 227 and Ala 235Gln 127—rs1215473912A235THis 231 and Phe 239His 231 and Phe 239Ala 191 and Val 187Ala 191(2): two interactions with the same residue.Mutated proteins were superimposed with the wild type to calculate the root-mean-square deviation (RMSD). For all structures, the RMSD (Å) values did not deviate significantly from the native protein (Table5).Table 5
Root-mean-square deviation (RMSD) of deleterious nsSNPs.
dbSNP IDSNPRMSD (Å)rs515726181R41W0.0274rs200273673R41P0.0236rs267607025R110P0.0446rs515726187F123S0.0032rs1211197364Y124C0.0019rs1185904532E131A0.0165rs515726188E131K0.0275rs1313813193S139R0.0964rs1169870960I142T0.0025rs515726190R186G0.0114rs515726191E194G0.0160rs121918308E194K0.0239rs515726192G195R0.0378rs863224192G200E0.0537rs1283028277F206V0.0059rs575109470G216R0.0684rs1422185855D227V0.0101rs772868983H231Q0.0181rs1215473912A235T0.0251


## 4. Discussion

Pathogenic variants affectingRRM2B protein structural properties can lead to remarkable mtDNA disruption whether on the qualitative (accumulation of multiple mtDNA deletions) or quantitative (depletion of mtDNA copy number) level, causing the appearance of serious phenotypes.In order to understand the molecular origin of these phenotypes, experimental validation seems to be crucial, although these methods require both time and resource investment. In this regard, computational tools may be an intelligent way to do this kind of studies, through the easy and open access to different software and servers dedicated to the extraction of useful information about the impact of mutations, these tools can orientate the experimental studies in the future to be more efficient and precise. Through this study, we attempted to investigate the effect of deleterious nsSNPs on the structure of theRRM2B protein, via different bioinformatic tools. After using a total of 13 algorithms, only nsSNPs jugged deleterious were taken into consideration allowing us to move from 8464 SNPs to 19 deleterious nsSNPs: eight SNPs in exon 6 (R186G, E194G, E194K, G195R, G200E, F206V, G216R, and D227V), 7 in exon 4 (R110P, F123S, Y124C, E131A, E131K, S139R, and I142T), 2 in exon 7 (H231Q and A235T), and 2 in exon 2 (R41W and R41P).A stability analysis was then conducted showing that only 1 nsSNP (Y124C) increased the stability of the RRM2B protein according to I-Mutant prediction, while for the Mupro tool, all the 19 nsSNPs decreased the stability. For the conservation analysis, the results have shown that the RRM2B protein is predicted to have a total of 225 conserved residues and that 16 out of the analyzed nsSNPs are localized in very highly conserved regions; 9 SNPs represent exposed and functional residues, while 7 SNPs represent buried and structural residues. These results may lead us to suspect that most of these variants may cause important alterations of the RRM2B protein.Our results have shown that the majority of variants are responsible for important impacts on the protein structure and amino acid interactions. Only one nsSNPs (F123S) had no impact on the amino acid interactions; these results cannot exclude the fact that this variant may have an important effect on the protein structure, since it was predicted to be deleterious, stability decreasing, and located in a highly conserved region.A variety of experimental trials concerning theRRM2B gene have shown that some of the 19 nsSNPs concerned by this study are responsible for serious manifestations within subjects of various ages and phenotypes. Mutation p.Arg41Trp was found in a patient with PEO, ptosis, migraine, neck flexion weakness (mild), and proximal muscle weakness, the age of onset was 60, and the patient died at the age of 66. While the heterozygous mutation p.Arg186Gly was found along with the p.Thr218Ile mutation in a 43-year-old patient with Severe PEO, asymmetrical ptosis, proximal and distal muscle weakness, ataxia, SNHL (sensorineural hearing loss), facial weakness, low BMI (Body Mass Index), leukoencephalopathy, and depression, with an age of onset of 11 years old, another patient with Severe PEO, ptosis, SNHL, dysphonia, proximal muscle weakness, ataxia, cataracts, glaucoma, and IHD (ischemic heart disease) had the p.Gly195Arg mutation, with an age of onset of 50 years old [15].On the other hand, Acham-Roschitz et al. reported a case of a baby girl of 8 weeks that suffered from muscular hypotonia with lack of head control as well as congenital deafness, nephrocalcinosis, and Fanconi syndrome. The genetic analysis revealed the presence of the p.Phe123Ser mutation [7]. The p.Glu194Lys mutation was reported in two siblings, a boy that showed trunk hypotonia and tubulopathy shortly after birth and died at the age of 2 months after status epilepticus and his sister that died at the same age and showed neonatal hyperlactatemia in addition to some similar complications as her brother. This mutation was also found in a baby girl that showed hypotonia soon after birth and developed respiratory distress, hyperlactatemia, and tubulopathy at an age of 3 months then died at the age of 4 months [5]. Finally, according to the study of Pitceathly et al., the p.Glu131Lys mutation was found along with p.Arg41Gln mutation in a patient with KSS syndrome [13].Many studies are aimed at explaining the mechanism of these diseases as well as the functional consequence of their corresponding variants. Actually, the p.Arg41Trp mutation is predicted to prevent salt bridge formation. This bridge has an important role in conformation change that controls iron binding; this prediction was based on previous molecular modeling studies showing that p.Arg41Gln prevents salt bridge formation [13, 33].As for mutations in the Arg 186 residue, it is demonstrated that they might be responsible for protein folding efficiency reduction leading to autosomal recessive disease, due to the presence of the residue at the end of theα-helix [15].Likewise, Gly 195 substitution impacts some of its neighboring amino acids that contribute to iron atom coordination; this disruption is associated with adult-onset chronic progressive external ophthalmoplegia (CPEO)/multiple mtDNA deletion disorders [13].All this experimental evidence correlates strongly with the results obtained from the bioinformatic study that we have conducted, thus proving and explaining the effect of deleterious nsSNPs on theRRM2B gene.


## 5. Conclusion

In the current study, we identified 19 nonsynonymous single nucleotide polymorphisms based on different computational prediction tools. To investigate the structural impact of these nsSNPs, molecular modeling analyses were done. A potential impact on protein stability was observed in 18 mutated structures in comparison with native protein. In addition, 18 variants may affect hydrogen bond and hydrophobic interactions. Our results allow for a better understanding of the effect ofRRM2B gene nsSNPs. Therefore, these SNPs should be considered important candidates in the genetic screening of neurological diseases, such as PEO and encephalomyopathy.



---

*Source: 7614634-2020-07-26.xml*