Next Article in Journal
Gene Variants of the OAS/RNase L Pathway and Their Association with Severity of Symptoms and Outcome of SARS-CoV-2 Infection
Previous Article in Journal
Bilateral Scaphoid Fractures: A Systematic Literature Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Analysis of Deleterious nsSNPs in INS Gene Associated with Permanent Neonatal Diabetes Mellitus

by
Elsadig Mohamed Ahmed
*,
Mohamed E. Elangeeb
,
Khalid Mohamed Adam
,
Hytham Ahmed Abuagla
,
Abubakr Ali Elamin MohamedAhmed
,
Elshazali Widaa Ali
,
Elmoiz Idris Eltieb
,
Ali M. Edris
,
Hiba Mahgoub Ali Osman
,
Ebtehal Saleh Idris
and
Khalil A. A. Khalil
Department of Medical Laboratory Sciences, College of Applied Medical Sciences, University of Bisha, P.O. Box 551, Bisha 61922, Saudi Arabia
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2024, 14(4), 425; https://doi.org/10.3390/jpm14040425
Submission received: 20 March 2024 / Revised: 6 April 2024 / Accepted: 10 April 2024 / Published: 17 April 2024
(This article belongs to the Section Omics/Informatics)

Abstract

:
Insulin gene mutations affect the structure of insulin and are considered a leading cause of neonatal diabetes and permanent neonatal diabetes mellitus PNDM. These mutations can affect the production and secretion of insulin, resulting in inadequate insulin levels and subsequent hyperglycemia. Early discovery or prediction of PNDM can aid in better management and treatment. The current study identified potential deleterious non-synonymous single nucleotide polymorphisms nsSNPs in the INS gene. The analysis of the nsSNPs in the INS gene was conducted using bioinformatics tools by implementing computational algorithms including SIFT, PolyPhen2, SNAP2, SNPs & GO, PhD-SNP, MutPred2, I-Mutant, MuPro, and HOPE tools to investigate the prediction of the potential association between nsSNPs in the INS gene and PNDM. Three mutations, C96Y, P52R, and C96R, were shown to potentially reduce the stability and function of the INS protein. These mutants were subjected to MDSs for structural analysis. Results suggested that these three potential pathogenic mutations may affect the stability and functionality of the insulin protein encoded by the INS gene. Therefore, these changes may influence the development of PNDM. Further researches are required to fully understand the various effects of mutations in the INS gene on insulin synthesis and function. These data can aid in genetic testing for PNDM to evaluate its risk and create treatment and prevention strategies in personalized medicine.

1. Introduction

Diabetes mellitus encompasses various metabolic disorders characterized by elevated blood sugar levels due to impaired insulin secretion or increased insulin resistance [1]. The dysfunction in the INS gene may contribute to the onset of diabetes mellitus [2]. The INS gene is located on chromosome 11 in humans. Insulin is the primary name for the gene [3].
The synthesis and production of insulin hormone take place in beta cells in the Langerhans islet of the pancreas [4,5]. These two processes begin with the synthesis of a longer precursor molecule called preproinsulin, the initial result of a translational product from the INS gene [2,4]. Preproinsulin is a 110-amino-acid polypeptide with a 24-amino-acid signal peptide. It is then cleaved to form proinsulin, which consists of an A chain, a B chain, and a connecting peptide (C-peptide). Proinsulin undergoes further cleavage to form mature insulin hormone and C-peptide. The C-peptide is released when insulin hormone is secreted and serves as a marker for its production [2,6]. The INS gene is transcribed into mRNA, which serves as a template for synthesizing insulin hormone. The mRNA is translated in the endoplasmic reticulum, and the nascent preproinsulin undergoes post-translational modifications, including signal peptide cleavage [5]. Insulin hormone is packaged into secretory granules within the pancreatic beta cells. It is released from beta cells in response to elevated blood glucose levels in the circulation to regulate and maintain the elevated blood glucose levels [7,8].
Neonatal diabetes typically manifests within the initial six months of life. Infants who experience the early onset of diabetes frequently possess mutations in one allele of the gene responsible for producing the insulin hormone [9]. Two types of neonatal diabetes mellitus have been reported: permanent neonatal diabetes mellitus (PNDM) and transient neonatal diabetes mellitus (TNDM) [10]. Mutations in the INS gene influence different phases of the biosynthesis of insulin hormone. These mutations are thought to interfere with either the cleavage process of the proinsulin chain or the assembly of the (A and B) chains to produce insulin, resulting in disturbed glucose levels in the blood [11,12]. INS gene mutations were described as a cause of PNDM [13,14].
As INS gene mutations were considered to affect the structure of insulin hormone and prescribed as a leading cause of PNDM, prediction of these mutations can lead to the early discovery or prediction of PNDM and better management and treatment of the disease. This study was conducted to identify potentially deleterious nsSNPs in the INS gene.

2. Methodology

2.1. Work Plan

This study utilized a computational approach consistent with other previous studies [15]. Different bioinformatics tools were used in this study to assess the deleterious effects of non-synonymous single-nucleotide polymorphisms (nsSNPs) in the INS gene. The nsSNPs considered to have high-risk profiles were analyzed to evaluate their stability, structural impact, and conservation. The methods used to identify and categorize putative functional nsSNPs in the INS gene are shown in Figure 1.

2.2. Data Collection

Data collection involved retrieving the nucleotide and amino acid sequences of the insulin protein, identified by accession numbers NP_001278826.1 and NG_050578.1, respectively, from the NCBI database in FASTA format. Information regarding single-nucleotide polymorphisms (SNPs) within the INS gene was obtained from the SNP database of NCBI via the following URL: http://www.ncbi.nlm.nih.gov/snp/ (accessed on 26 September 2023). Data of the INS gene and insulin protein were sourced from the PDB and UniProt databases, accessible using the web address https://www.uniprot.org/uniprotkb/P01308/entry (accessed on 26 September 2023).

2.3. Investigation of the Impact of nsSNPs on Protein Function

The impact of genetic variation on protein function was investigated using several tools, including Sorts Intolerant from Tolerant (SIFT), PolyPhen2, and SNAP2. The SIFT sequences were performed through the website https://sift.bii.a-star.edu.sg/www/SIFT_seq_submit2.html (accessed on 26 September 2023). SIFT predicts changes in protein function resulting from amino acid sequence alterations [16]. The website http://genetics.bwh.harvard.edu/pph2 (accessed on 26 September 2023) was accessed for PolyPhen2. The PolyPhen2 predicts the potential impact of amino acid substitutions on protein function and structure [17]. SNAP2 was accessed using https://www.rostlab.org/servces/snap/ (accessed on 26 September 2023). It predicts changes in secondary protein structure due to nsSNPs [18].
The SIFT, PolyPhen2, and SNAP2 tools were employed to assess the potential effects of SNPs; each offer unique approaches for classifying nsSNPs and determining confidence scores using various input criteria. These tools were employed to assess the functional consequences of non-synonymous single amino acid substitutions, differentiating between substitutions that are tolerable and those that are detrimental to protein function. SIFT predicts the impact of an amino acid change on protein function by analyzing sequence homology and the physical properties of amino acids [16]. A tolerance index score is assigned by SIFT, with scores below 0.05 indicating deleterious variations [16]. The PolyPhen2 assigns a numerical value ranging from 0 to 1, to forecast the potential effect of an amino acid substitution on the structure and function of a protein. A value of 0 indicates no impact, while a value of 1 indicates the most harmful consequence [17]. In contrast, SNAP2 conducts a comparative analysis of genomes to predict potential functional consequences of single amino acid substitutions at the protein level [18]. Deleterious SNPs identified through all mentioned servers underwent further computational analysis.

2.4. Prediction of SNP Disease Association

The correlation between filtered SNPs and diseases was assessed using SNPS & GO and Predictor of human Deleterious Single Nucleotide Polymorphisms (PhD-SNP) [19,20], accessible at http://snps-and-go.biocomp.unibo.it/snps-and-go/ (accessed on 26 September 2023) and http://snps.biofold.org/phd-snp/phd-snp.html (accessed on 26 September 2023), respectively. SNPS & GO utilizes Support Vector Machine (SVM) computational methods to detect harmful single amino acid substitutions, incorporating Gene Ontology (GO) annotations for classification. Variants scoring above 0.5 are deemed disease-causing [19]. The PhD-SNP utilizes a support vector machine (SVM) classifier to classify disease-associated SNPs. During the classification method, the categorization of amino acid changes as either neutral or disease-associated relies heavily on the exploitation of sequence and profile data [20]. The sequence profile calculation entails utilizing an input vector derived from the amino acid frequencies in both the wild-type and mutant sequences, the number of matched sequences, and the conservation score at the replacement site. If the PhD-SNP score exceeds 0.5, it indicates the existence of a mutation that is accountable for the onset of a disease [20]. The SNPs that were found to be responsible for causing disease and came across all mentioned servers were subjected to further computational analysis.

2.5. Prediction of SNPs on INS Protein Function Related to Pathogenicity

The MutPred2 web server, available at http://mutpred.mutdb.org (accessed on 26 September 2023), is used to classify mutations as either neutral or associated with disease [21,22]. This tool applies a machine-learning approach to predict the molecular mechanisms underlying the pathogenicity of amino acid substitutions. By analyzing a wide array of fifty different protein properties, it evaluates the effects of these substitutions. A score greater than 0.5 indicates a mutation that may be pathogenic [21].

2.6. Impact of SNPs on Protein Stability

The I-Mutant algorithm, leveraging Support Vector Machine (SVM) technology, forecasts changes in protein stability induced by missense nsSNPs [23]. Its output, known as the reliability index (RI), is scaled from 0 to 10, where a score of 10 represents the utmost reliability. I-Mutant is available at http://folding.biofold.org/i-mutant/imutant2.0.html (accessed on 26 September 2023). It requires the amino acid sequence of the INS protein and details about the mutating residues, such as their locations, for its analysis. Additionally, the MuPro tool, accessible at https://mupro.proteomics.ics.uci.edu (accessed on 26 September 2023), is employed for evaluating protein stability [24]. MuPro utilizes machine learning approaches, including SVM and Neural Networks, to predict how individual amino acid alterations affect protein stability. These methods were honed on an extensive database of mutations and validated through a 20-fold cross-validation process, achieving an accuracy rate of over 84%. Scores range between −1.0 and 1.0 to signify prediction reliability. A significant benefit of these approaches is their capacity to foresee alterations in protein stability without needing information on the tertiary structure.

2.7. Analysis of Protein Properties Utilizing the HOPE Tool

The HOPE software tool, which can be accessed at the URL https://www.cmbi.umcn.nl/hope (accessed on 27 September 2023), is designed to streamline the analysis of protein mutations [25]. It harnesses data from the UniProt database to assess the effects of single-nucleotide variations on the structural and functional aspects of proteins. HOPE produces in-depth reports that encompass textual explanations, graphical representations, and interactive visualizations, offering a detailed view of the mutation’s consequences. Mutants identified as deleterious nsSNPs by the previously mentioned tools were chosen for further scrutiny. These selected mutants underwent Molecular Dynamics Simulation (MDS) for an enhanced level of analysis.

2.8. Performing of Molecular Dynamics Simulations (MDSs)

To study changes in the structure over time in wild-type and mutant structures, Molecular Dynamics Simulations (MDSs) were carried out using GROMACS version 2020.6 on a Google Collaboration Pro notebook with substantial RAM. The GROMACS-OPLS-AA force field was used for initial calculations. Both systems were immersed in cubic boxes filled halfway with water molecules, with a minimal radius of 1 nm. To balance the system, 10 sodium ions (Na+) were introduced using the GROMACS genion tool. Energy minimization was achieved through the steepest descent algorithm, with an energy step size of 0.01 and a maximum of 50,000 iterations. The system’s stability was preserved at a Parrinello–Rahman pressure of 1 bar and a Berendsen temperature of 300 K. Electrostatic interactions were computed using the Partial Mesh Ewald (PME) technique, with a short-range cut-off for electrostatic and van der Waals interactions set at 1.0 nm. The neighbor list was updated every 10 ps. The LINCS algorithm, with a time step of 0.002 ps, maintained all bond constraints, including heavy atom-H bonds. An isothermal compressibility of 4.5 × 10−5 was used, with coupling constants for temperature and pressure set to 0.1 ps and 2.0 ps, respectively. The system underwent equilibration for 100 ps in both the NPT (constant number of particles, pressure, and temperature) and NVT (constant number of particles, volume, and temperature) ensembles, with a Berendsen temperature of 300 K and a Parrinello–Rahman pressure of 1 bar. Trajectories were recorded every 1 ps during 10 ns of molecular dynamics simulations for both wild-type and mutant structures. Various analysis tools like g_rms, g_rmsf, g_sasa, g_Rg, and g_density were utilized to assess parameters such as the RMSD (root-mean-square deviation), RMSF (root-mean-square fluctuation), SASA (solvent accessible surface area), and Rg (radius of gyration) to compare structural differences between the wild-type and mutant structures.

2.9. Displaying Three-Dimensional Structural Change Using PyMol Software

A 3D simulation of the wild-type INS protein and putative mutations was performed to display structural changes using the PyMol software Version 2.0.

3. Results

This study aimed to identify the potential deleterious nsSNPs in the INS gene by implementing computational algorithms including SIFT, PolyPhen2, SNAP2, SNPs & GO, PhD-SNP, MutPred2, MuPro, and I-Mutant. These algorithms analyze the influence of nsSNPs on Insulin protein function, stability, and disease association. Subsequently, the HOPE software tool was utilized to evaluate the effect of nsSNPs on the structural and functional properties of the INS gene. The nsSNPs predicted by all utilized algorithms were analyzed for structural changes over time, in both wild-type and mutant structures, through MDSs using the GROMACS tool.

3.1. Data Collection

SNPs associated with the INS gene and its corresponding Insulin protein sequence, identified by UniProt ID P01308, were obtained from the NCBI dbSNP and UniProt Knowledgebase databases. A total of 8,147,365 SNPs were mapped to the INS gene sequence. Subsequently, 130 nsSNPs were detected in the coding region. These nsSNPs are suspected to have the potential to cause missense mutations and significantly impact protein structure and function. The nsSNPs present in the coding sequence of the INS protein were considered for our analysis.

3.2. Exploring the Influence of SNPs on Protein Function

The analysis in this study focused on predicting nsSNPs deleteriousness on the INS gene to determine their potential impact on structural and functional properties. Out of 130 nsSNPs, 86 were classified as deleterious by the SIFT tool. Additionally, 67 nsSNPs were identified as probably damaging and potentially harmful by the PolyPhen2 tool. Confirmation of the harmful effect of these nsSNPs was conducted using SNAP2. All identified nsSNPs using the PolyPhen2 tool were confirmed as deleterious by SNAP2. Thus, a total of 67 nsSNPs were identified as the most deleterious based on the above-mentioned software (Table 1). These 67 nsSNPs were used for further analysis in the next predicting tools.

3.3. Prediction of SNP–Disease Association

A combined approach utilizing the SNP & GO and PhD-SNP tools was used to analyze the 67 nsSNPs for disease susceptibility. Results of SNP & GO suggested that all 67 nsSNPs were considered to have the potential to be associated with disease, with scores ranging from 9 to 10. However, PhD-SNP analysis identified only 12 nsSNPs out of the 67 as significantly correlated with disease. The remaining 55 nsSNPs showed neutral effects (Table 2). These 12 variants were subjected to additional analysis to obtain their potential impact.

3.4. Prediction of SNPs on INS Protein Function Related to Pathogenicity

This study examines the genetic pathways that cause disease and presents the probability ratings obtained from analyzing 12 different nsSNPs using the MutPred2 server. MutPred2 analysis categorized 8 out of the 12 examined nsSNPs as variants associated with pathological conditions, indicating alterations in the protein’s structure (Table 3). Changes included a loss of helix, disulfide linkage, and loop, and altered transmembrane protein, signal peptide, ordered interface, and metal binding. However, MutPred2 predicted a neutral impact on protein function for four mutants.

3.5. Predicting the Impact of SNPs on Protein Stability

I-Mutant2.0 and MuPro software were used to investigate the impact of predicted deleterious mutations on INS protein stability. Findings derived from I-Mutant2.0 exhibited coherence, indicating that seven nsSNPs were anticipated to diminish the protein’s stability in a decreasing manner. However, only one nsSNP was found to increase protein stability. However, findings derived from MuPro indicated that four nsSNPs out of the seven were anticipated to diminish the protein’s stability. Three nsSNPs were found to increase protein stability (Table 4).

3.6. Analysis of Protein Properties

The HOPE tool analysis revealed that three nsSNPs were predicted to damage protein structure. However, L14Q was found not to affect the protein structure (Table 5). From the HOPE report, the mutations (P52R) Proline into Arginine at position 52, (C96Y) Cysteine into Tyrosine at position 96, and (C96R) Cysteine into Arginine at position 96 showed that the mutant residue was bigger than the wild-type residue in all mutations. And the Wild-type residue was more hydrophobic than the mutant residue. However, the Wild-type residue charge was NEUTRAL and the Mutant residue charge was POSITIVE in mutations (P52R and C96R), with no charge change in mutation (C96Y). The Mutant and the wild-type residues were not very similar in all mutants (P52R, C96Y, and C96R). All three mutations were considered probably damaging to the protein.

3.7. Molecular Dynamics Simulation (MDS) Results

Figure 2 displays a Root-Mean-Square Deviation (RMSD) analysis of proteins throughout molecular dynamics simulations (MDSs). The RMSD graphic demonstrates the conformational stability of a wild-type protein in comparison to its three variations, C96R, C96Y, and P52R, throughout a 1000 ps simulation period. Each variety is represented by a distinct colored line, whereas the wild type is displayed in black. The path of the wild-type protein is used as a standard to evaluate how mutations affect protein stability. Upon initial examination, all protein variations, including the wild-type, exhibited a rise in RMSD values as the simulation advances, a common occurrence as the protein navigates through different conformations. The wild-type stabilized at about 0.6 nm, showing a plateau after initial fluctuations. The C96R mutation, highlighted in red, showed a comparable level of stabilization, indicating that this mutation does not greatly affect the protein’s structural stability. The C96Y variant, shown in green, displayed a slightly greater RMSD, reaching approximately 0.7 nm, suggesting a potential minor increase in structural flexibility or a departure from the original structure. The P52R mutation, highlighted in blue, exhibited a consistent rise in RMSD, peaking at about 1.5nm, indicating significant structural alterations or instability caused by this mutation. The RMSD study shows that the P52 R mutation has a large effect on the protein’s structural stability, whereas the C96R mutation has a minimal influence, and the C96Y mutation causes a minor increase in flexibility. These results suggest that the mutation at position 52 causes a greater departure from the protein’s original structure compared to mutations at position 96. None of the variants exhibited symptoms of reaching a plateau within the 1000 ps duration. The wild-type protein remained stable for around 200 ps, serving as a reference point for studying the characteristics of its many forms. The RMSD analysis data are crucial for comprehending the structural consequences of these mutations and their possible influence on the protein’s function.
Concerning the results presented in Figure 3, the graph displays the Root-Mean-Square fluctuation (RMSF) of individual atoms in a protein throughout the MD simulation. The RMSF quantifies the degree of atomic movement from their mean position, indicating the variability in flexibility across various regions of the protein structure. The RMSF graph illustrates the flexibility of a wild-type protein and its variations, C96R, C96Y, and P52R, spanning approximately 2000 atoms. A line graph illustrates the flexibility of each protein variant, with distinct colors representing each variant, including the wild-type one. Fluctuations are quantified in nanometers (nm), with the vertical axis indicating the extent of fluctuation and the horizontal axis indicating the atom’s position in the protein. Upon analysis, it is clear that the wild-type protein displayed a fluctuation pattern that acts as a reference for comparing with the mutations. All variations exhibited regions with high and low fluctuation, indicating sections of the protein that are either more flexible or stiffer. The C96R variant (in red) typically aligned with the behavior of the wild type, indicating that the mutation does not significantly alter the protein’s flexibility. The C96Y variant (in green) and P52R variant (in blue) exhibited higher fluctuation in certain places, suggesting probable areas of enhanced flexibility or disorder perhaps caused by the mutations. The P52R variant displayed prominent peaks, indicating potential structural changes or increased mobility in certain protein areas relative to the wild-type and other variants. The data show that the P52R mutation caused significant variation in atomic locations, which could impact the structural integrity and function of the protein. The C96Y mutation enhanced flexibility to a lesser extent than P52R. The C96R mutation had little effect on flexibility, closely approaching the behavior of the wild-type protein. Fluctuations in RMSF are important for comprehending the functional impacts of mutations. Areas with increased fluctuations may indicate active or binding sites, or structural domains crucial for protein stability.
In Figure 4, the radius of gyration (Rg) of a wild-type protein and its three variations, C96R, C96Y, and P52R, as a function of time throughout a simulation duration of 1000 ps is shown. The Rg quantifies the compactness of a protein and the distribution of its mass relative to an axis. The figure illustrates the difference in compactness between the wild-type protein and its variations during the MD simulation. The radius is quantified in nanometers (nm) and monitors alterations in the protein’s morphology and dimensions over time. The graph displays the average distance of the protein’s atoms from its center of mass, with each line representing a different version indicated by different colors. A larger radius of gyration indicates a less compact protein structure, whereas a smaller number suggests a more compact structure.
During the simulation, all protein structures displayed variations in their radius of gyration, suggesting alterations in compactness and tertiary structure as time progressed. The wild-type protein exhibited varying levels of compactness, with the radius of gyration values oscillating between around 2.5 nm and 2.9 nm. The C96R variation (red line) displayed comparable fluctuations, indicating that this mutation does not substantially change the overall compactness of the protein. The C96Y variant (green line) showed a somewhat larger radius of gyration across most of the simulation, suggesting a less dense structure compared to the wild type. The P52R variant showed a unique pattern where its radius of gyration stayed consistently lower than the others during the initial half of the simulation but eventually increased to reach the highest level of extension. This indicates a notable structural change caused by the mutation. The data from the radius of gyration research indicate that the P52R mutation has the most significant impact on the protein’s compactness, potentially influencing its stability and function. The C96Y mutation slightly reduced the protein’s compactness; however, the C96R mutation had less impact on the protein’s compactness relative to the wild type. The differences in the radius of gyration, especially in the P52R variant, may indicate modified functional dynamics or alterations in the protein’s capacity to interact with other molecules. These findings can impact the comprehension of the structural and functional characteristics of the protein being analyzed.
Figure 5 shows the Solvent Accessible Surface Area (SASA) of a wild-type protein and three variations, C96R, C96Y, and P52R, during a 1000 ps molecular dynamics simulation. SASA is a crucial structural metric that represents the surface area of a protein available for interaction with a solvent, usually water. This parameter is pertinent to the stability, folding, and interactions of proteins. The SASA figure displays how the surface area exposed to solvent changed over time for both the wild-type and mutant proteins during the simulation. Values are expressed in square nanometers (nm2), with time shown in ps. The wild-type protein is shown by the black line, and the mutants are illustrated in red (C96R), green (C96Y), and blue (P52R). Initial observations demonstrate that all protein variations exhibited a reduction in SASA over time during the simulation, indicating a possible shift towards a more condensed structure or a modification in conformation dynamics, leading to decreased solvent exposure. The C96R variant initially exhibited SASA similar to the wild type but deviated towards the end, suggesting minor variations in exposure or folding patterns that are absent in the wild type. Throughout the simulation, the C96Y and P52R variations consistently displayed higher SASA values than the wild type, with the C96Y variant displaying the most notable difference. This may suggest a more open conformation or changed folding pattern that enhances the solvent accessibility of specific residues. The P52R variation had higher SASA values than the wild type and displayed a complicated pattern with frequent fluctuations, indicating a potentially more dynamic or flexible protein structure. This study indicates that mutations at positions 96 and 52 impact the protein’s interaction with the solvent, potentially influencing protein stability and function. The heightened solvent exposure in the C96Y and P52R variations could impact the protein’s surface chemistry and its ability to interact with other biomolecules. The reduced SASA observed in the later stages of the simulation for the C96R variant compared to the wild type may suggest changes in the tertiary structure or dynamics that could impact the protein’s biological function. The results emphasize the significance of SASA in elucidating the impact of mutations on protein structure and function.
Figure 6 quantifies the hydrogen bonds present in a wild-type protein and its three variations, C96R, C96Y, and P52R, over a 1000 ps molecular dynamics simulation. Hydrogen bonds play a crucial role in maintaining the stability of protein secondary and tertiary structures and can impact protein dynamics. The graph measures the hydrogen bonding arrangements in the molecular structures of the wild-type and mutant proteins as time progresses. The number of hydrogen bonds in the wild-type protein is represented by a black line that varies during the simulation, demonstrating the dynamic nature of intramolecular interactions. The mutations C96R, C96Y, and P52R are depicted by red, green, and blue lines, respectively, each displaying a unique hydrogen bonding arrangement in comparison to the original type. The comparison of hydrogen bond dynamics shows that the wild-type protein typically maintained a specific range of hydrogen bonds, mostly between 35 and 55. The C96R variant closely resembled the wild type in the amount of hydrogen bonds, indicating that this mutation does not significantly change the hydrogen bonding network. The C96Y variation had a wider range and a slightly larger quantity of hydrogen bonds, suggesting potential changes in intramolecular interactions or the enhanced stability of the secondary structure. The P52R variation exhibited the highest variety in hydrogen bond numbers, potentially indicating alterations in the protein’s structural stability and folding pattern as a result of the mutation. Proteins’ structural integrity and function are frequently dependent on their hydrogen bonding networks. The data indicate that the P52R mutation had the most significant impact on hydrogen bond formation during the simulation, potentially leading to substantial changes in protein structure and dynamics. The C96Y mutation exhibited enhanced hydrogen bonding but did not display the same level of variability as the P52R mutation, indicating a less significant effect on the protein’s stability. The similarity in hydrogen bond patterns between the C96R variation and the wild type suggests that this mutation has a limited impact on the protein’s structural stability. The observations could impact the protein’s heat stability, folding dynamics, and interaction with other molecules.

3.8. Displaying Three-Dimensional Structural Changes Using PyMol Software

Figure 7 displays the three-dimensional structure of the INS protein (PDB ID: 3i40), the wild type, and its mutations P52R, C96R, and C96Y. The sites of the mutations are shown in the red color. Figure 8, Figure 9, Figure 10 and Figure 11 display the three-dimensional structure of the wild type and its mutations, P52R, C96R, and C96Y, respectively.

4. Discussion

In this study, we investigated the impact of nsSNPs in the INS gene on PNDM using computational methods. In this research, we sought to enhance the understanding of genetic factors influencing PNDM and their potential implications for personalized medicine. Our findings are expected to highlight the molecular basis of the illness, identify new diagnostic indicators, and uncover potential treatment targets. This study aligns with research that identified the association of nsSNPs with diabetes mellitus, and research emphasized the relevance of genetic changes in the INS gene in the development of diabetes mellitus [26].
Many computational methods were used to analyze the impact of nsSNPs on the INS gene protein structure and function that are potentially contributing to the onset of PNDM. The approach outlined in this study for predicting the deleteriousness of nsSNPs within the INS gene is comprehensive and systematic, utilizing multiple tools to assess the potential impact on protein structure and function. The use of SIFT, PolyPhen2, and SNAP2 tools provides a multifaceted analysis, enhancing confidence in identifying harmful variants (Table 1). Comparing our study with other similar research using computational tools can provide insights into the consistency and reliability of our findings. A study by Kumar et al. (2009) conducted a comprehensive evaluation of the SIFT and PolyPhen2 algorithms for predicting the impact of amino acid substitutions on protein function. Kumar et al. (2009) found that while both tools generally perform well, there can be discrepancies in their predictions, highlighting the importance of using multiple prediction methods for increased accuracy [27]. Furthermore, the incorporation of SNAP2 as a confirmation tool adds another layer of validation to the predictions made by PolyPhen2, enhancing the robustness of the findings. Additionally, a study performed by Ioannidis et al. (2016) emphasized the need for caution when interpreting computational predictions of nsSNP deleteriousness, highlighting the potential for false positives and false negatives. They underscored the importance of experimental validation to confirm the functional impact of predicted variants [28]. Our study employed a rigorous approach to predict the deleteriousness of nsSNPs within the INS gene, incorporating multiple computational tools and validation steps. This approach enhanced confidence in the identified deleterious variants.
A combined approach was utilized using SNP & GO and PhD-SNP tools to analyze the disease susceptibility of nsSNPs within the INS gene to provide valuable insights into their potential association with diseases. SNP & GO analysis provides a score range from 9 to 10, suggesting a potential association of 67 nsSNPs with the disease (Table 2). This method likely considers functional annotations and gene ontology terms associated with the SNPs to infer their potential disease relevance. On the other hand, the PhD-SNP methodology identifies only 12 out of the 67 nsSNPs as significantly correlated with disease, with the remaining 55 nsSNPs predicted to have neutral effects (Table 2). A study conducted by Martelotto et al. (2014) evaluated the performance of SNP & GO in predicting the disease association of nsSNPs. They found that while SNP & GO can provide valuable insights into the potential functional impact of nsSNPs, it may also yield false positive predictions due to the reliance on gene ontology annotations, which can be incomplete or inaccurate [29]. In contrast, the PhD-SNP methodology employs a different approach, likely incorporating features such as protein structure, evolutionary conservation, and functional annotations to predict disease association. A study conducted by Bendl et al. (2014) demonstrated the effectiveness of PhD-SNP in identifying disease-associated nsSNPs, highlighting its utility in prioritizing variants for further analysis [30]. The identification of 12 nsSNPs as significantly correlated with disease by PhD-SNP suggests a subset of variants with a higher likelihood of impacting disease susceptibility. While SNP & GO and PhD-SNP provide valuable insights into the potential disease association of nsSNPs, the differences in their predictions highlight the need for the careful consideration and validation of computational findings.
Then, MutPred2 was applied to analyze nsSNPs within the INS gene to provide valuable insights into nsSNPs’ potential impact on protein functionality related to pathogenicity. The identification of eight nsSNPs associated with pathological conditions suggests a subset of variants with a higher likelihood of disrupting normal protein function (Table 3). The specific alterations predicted by MutPred2, such as the loss of helix, disulfide linkage, loop, transmembrane protein, signal peptide, ordered interface, and metal binding, provide mechanistic insights into how these nsSNPs may lead to pathogenicity. Comparing these results with findings from other analyses, such as results obtained from SIFT, PolyPhen2, SNAP2, SNP & GO, and PhD-SNP, can help corroborate and contextualize the predictions made by MutPred2. For instance, if an nsSNP is predicted to be deleterious by multiple algorithms and also associated with disease by SNP & GO and PhD-SNP, it strengthens the confidence in its potential pathogenicity, as indicated by MutPred2.
A study by Pejaver et al. (2017) evaluated the performance of MutPred2 in predicting the pathogenicity of nsSNPs and found it to be one of the top-performing tools in terms of sensitivity and specificity. However, they also emphasized the importance of integrating predictions from multiple tools and experimental validation to improve the accuracy of pathogenicity predictions [21]. The MutPred2 analysis provides valuable insights into the potential impact of nsSNPs on INS protein functionality related to pathogenicity [21].
Assessment of predicted deleterious mutations on INS protein stability using both I-Mutant2.0 and MuPro tools provides valuable insights into the potential impact of nsSNPs on protein stability. However, the discrepancies between the predictions obtained from these two tools warrant careful consideration and comparison with findings from the other analyses. In this study, results from I-Mutant2.0 suggested that seven out of eight of the analyzed nsSNPs are expected to diminish protein stability, with only one nsSNP predicted to increase stability. Conversely, the findings from MuPro indicated that four out of seven nsSNPs are anticipated to diminish protein stability, while three are predicted to increase stability (Table 4). Comparing these results with findings from other computational tools, such as SIFT, PolyPhen2, SNAP2, and MutPred2, can help assess the consistency of predictions regarding the impact of nsSNPs on protein stability. If an nsSNP is consistently predicted to be deleterious by multiple tools and associated with disease susceptibility, it strengthens the confidence in its potential to disrupt protein stability and function.
Khan, et al. (2010) evaluated the performance of I-Mutant2.0 and MuPro in predicting the impact of nsSNPs on protein stability. They found that while these tools can provide useful predictions, there can be discrepancies between their results due to differences in underlying algorithms and features considered. While the predictions from I-Mutant2.0 and MuPro provide insights into the potential impact of nsSNPs on INS protein stability, discrepancies between the results highlight the importance of considering multiple computational tools [31].
The use of the HOPE tool to evaluate the effects of single-nucleotide variations (SNVs) on the structural and functional properties of proteins, specifically focusing on the four nsSNPs predicted to diminish protein stability by the MuPro tool, provides additional insights into their potential impact. The identification of three nsSNPs (C96Y, P52R, and C96R) as potentially damaging to the protein’s structure by the HOPE tool aligns with findings from other computational analyses, such as MutPred2, which predicted alterations in protein structure for these variants. This consistency across multiple prediction methods strengthens the confidence in the detrimental effects of these nsSNPs on protein structure and function. Conversely, our study finding that the L14Q variant was not predicted to affect protein structure by the HOPE tool (Table 5) contrasts with its classification as deleterious by other computational tools, such as MutPred2, I-Mutant, and MuPro. This discrepancy highlights the importance of considering multiple prediction methods and experimental validation to accurately assess the impact of nsSNPs on protein structure and function.
Venselaar et al. (2010) evaluated the performance of the HOPE tool in predicting the effects of mutations on protein structure and function. They found that HOPE can aid in understanding the potential impact of gene mutations on protein stability and function [25]. Therefore, using the HOPE tool to evaluate the effects of nsSNPs on protein structure and function complements findings from other computational analyses.
Our investigations pinpointed four variations as the most harmful using the combined results of different analytical methods. Three possible mutants, named C96Y, P52R, and C96R, were chosen among the four nsSNPs, identified using the HOPE tool (Table 5) as the most harmful mutants, for further study using molecular dynamic simulations (MDSs). The detected mutations were shown to alter the protein’s stability and functional characteristics, potentially affecting its biological activity. A molecular dynamics simulation for these three nsSNPs was conducted to assess the differences in protein dynamics between the wild-type and mutant INS proteins (Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6).
In the analysis of nsSNPs in the INS gene, molecular dynamics simulations can be employed to understand the structural dynamics and stability of the INS protein affected by the variants.
RMSD analysis evaluated the overall structural stability of proteins during the MD simulations. Our findings suggest that while the C96R mutation has minimal impact on stability, the C96Y mutation causes a slight increase in flexibility, and the P52R mutation induces significant structural alterations, as evidenced by a higher RMSD value. These observations align with computational predictions and highlight the structural consequences of these mutations. RMSF analysis provides insights into the flexibility of individual atoms within the protein structure. Results indicate that the C96R mutation closely resembles the wild type, while the C96Y and P52R mutations exhibit increased fluctuations in certain regions, suggesting potential alterations in flexibility and structural dynamics caused by these mutations. The radius of gyration reflects the compactness of the protein structure. Data showed that while the C96R mutation has minimal impact on compactness, the C96Y mutation slightly reduces it, and the P52R mutation significantly increases it, indicating substantial structural changes induced by this mutation. SASA analysis revealed changes in the protein’s interaction with the surrounding solvent. Results suggest that the C96Y and P52R mutations increase solvent accessibility, potentially affecting surface chemistry and molecular interactions, whereas the C96R mutation has a limited impact on solvent exposure. Hydrogen bonds are crucial for maintaining protein stability and structure. Data indicate that while the C96R mutation maintains a similar hydrogen bonding pattern to the wild type, the C96Y mutation exhibits enhanced hydrogen bonding, and the P52R mutation displays variability, suggesting significant alterations in protein stability and folding dynamics.
These computational analyses provide valuable insights into the structural and functional consequences of specific mutations in the insulin protein. The detected mutations were shown to alter the protein’s stability and function, potentially affecting its biological activity. A molecular dynamics simulation was conducted to assess the differences in protein dynamics between the wild-type and mutant INS proteins. The simulation lasted 1000 picoseconds and involved analyzing various metrics such as the RMSD, RMSF, SASA, and RG. Comparisons were conducted to analyze the functional traits and robustness of the protein structures. After minimizing, marginal variances in the RMSD plots of backbone atoms showed that the wild-type structure was more stable than the mutations. The RMSF plots showed that the native and mutant structures displayed the highest degree of flexibility. The RMSF estimates averaged over residual elements and trajectories showed that atom fluctuation within the structures was significantly higher in the second third. The inquiry analyzed the radius of gyration and SASA to determine how the unfolding of the mutant protein structure affects the solubility of surface area and protein activity. A discrepancy in RG values suggests decreased stability. The current results indicate the original structure had a more limited range of RG values than the mutated structures. The solvent-accessible surface area of the wild-type INS differed from that of the mutations of the INS protein. These three mutations, C96Y, P52R, and C96R, are suggested to cause a decrease in INS stability and functionality. Figure 7 displays the three-dimensional structure of the selected INS gene mutations, that are likely suggested to introduce structural disruptions potentially affecting the stability and function of the protein.
From the literature, Laurenzano et al., in 2019, detected a P9R missense variant of the INS gene in a patient, which occurs in the signal peptide of the preproinsulin molecule in an amino acid residue that is highly conserved [32]. Also, the results of Balboa et al., in 2018, demonstrated that neonatal diabetes-associated INS mutations lead to defective beta-cell mass expansion, leading to diabetes development [9]. Moreover, Liu et al., in 2015, reported that a preproinsulin-L13R mutation was identified experimentally in a baby girl who developed severe diabetes on the second day after birth [8]. Our results differ from other results in the literature. Hence, we compare the results of our study with the results of other published studies on INS gene mutations associated with PNDM in Table 6.
However, additional research may be required to link these computational results with experimental data to understand the effects of these mutations on protein behavior fully.

5. Conclusions

This study delved into the assessment of the potential correlation between non-synonymous single-nucleotide polymorphisms (nsSNPs) within the INS gene and PNDM utilizing an array of computational methodologies. Techniques such as SIFT, PolyPhen2, Mutpred2, SNAP2, SNP & Go, PhD-SNP, MuPro, I-Mutant, HOPE, and MDSs were employed. Among the findings, three mutations, C96Y, P52R, and C96R, were identified to diminish the stability and functionality of the INS protein. These results suggest that these mutations could impact the stability and function of the insulin protein encoded by the INS gene, potentially influencing the onset of PNDM. Further investigations are warranted to comprehensively grasp the diverse ramifications of INS gene mutations on insulin production and function. Upon validation, these findings could be instrumental in genetic screenings for PNDM, aiding in risk assessment and the formulation of therapeutic or preventive measures.

Author Contributions

Conceptualization, E.M.A., M.E.E., A.A.E.M., E.W.A., E.I.E., H.A.A., A.M.E., H.M.A.O., E.S.I., K.A.A.K. and K.M.A.; Methodology, E.M.A., M.E.E. and K.M.A.; Software, E.M.A., M.E.E. and K.M.A.; Validation, E.M.A., M.E.E., K.M.A., H.A.A., A.A.E.M., E.W.A., E.I.E., A.M.E., H.M.A.O., E.S.I. and K.A.A.K.; Formal analysis, E.M.A., M.E.E. and K.M.A.; Data curation, E.M.A., M.E.E. and K.M.A.; Writing—original draft preparation, E.M.A., M.E.E. and K.M.A.; Writing—review and editing, E.M.A., M.E.E., K.M.A., A.A.E.M., E.W.A., E.I.E., H.A.A., A.M.E., H.M.A.O., E.S.I. and K.M.A.; Funding acquisition, E.M.A., M.E.E. and K.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data concerning SNPs associated with the INS gene, Insulin protein sequence identified by UniProt ID P01308, and INS nsSNPs are available at the NCBI database, via the URL: http://www.ncbi.nlm.nih.gov/snp/ (accessed on 26 September 2023) and UniProt Knowledgebase database.

Acknowledgments

The authors are thankful to the University of Bisha, Saudi Arabia for the financial support through the Research Support Program for the Publication in the Ranked Journals.

Conflicts of Interest

The authors declare that no conflicts of interest exist.

References

  1. Dilworth, L.; Facey, A.; Omoruyi, F. Diabetes Mellitus and Its Metabolic Complications: The Role of Adipose Tissues. Int. J. Mol. Sci. 2021, 22, 7644. [Google Scholar] [CrossRef]
  2. Wang, H.; Saint-Martin, C.; Xu, J.; Ding, L.; Wang, R.; Feng, W.; Liu, M.; Shu, H.; Fan, Z.; Haataja, L.; et al. Biological behaviors of mutant proinsulin contribute to the phenotypic spectrum of diabetes associated with insulin gene mutations. Mol. Cell. Endocrinol. 2020, 518, 111025. [Google Scholar] [CrossRef] [PubMed]
  3. Owerbach, D.; Bell, G.I.; Rutter, W.J.; Brown, J.A.; Shows, T.B. The insulin gene is located on the short arm of chromosome 11 in humans. Diabetes 1981, 30, 267–270. [Google Scholar] [CrossRef] [PubMed]
  4. Murthy, M.H.; Jasbi, P.; Lowe, W.; Kumar, L.; Olaosebikan, M.; Roger, L.; Yang, J.; Lewinski, N.; Daniels, N.; Cowen, L.; et al. Insulin signaling and pharmacology in humans and in corals. PeerJ 2024, 12, e16804. [Google Scholar] [CrossRef] [PubMed]
  5. Ataie-Ashtiani, S.; Forbes, B. A review of the biosynthesis and structural implications of insulin gene mutations linked to human disease. Cells 2023, 12, 1008. [Google Scholar] [CrossRef] [PubMed]
  6. Jamiołkowska-Sztabkowska, M.; Głowińska-Olszewska, B.; Bossowski, A. C-peptide and residual β-cell function in pediatric diabetes—State of the art. Pediatr. Endocrinol. Diabetes Metab. 2021, 27, 123–133. [Google Scholar] [CrossRef]
  7. Rorsman, P.; Ashcroft, F.M. Pancreatic β-Cell Electrical Activity and Insulin Secretion: Of Mice and Men. Physiol. Rev. 2018, 98, 117–214. [Google Scholar] [CrossRef]
  8. Liu, M.; Sun, J.; Cui, J.; Chen, W.; Guo, H.; Barbetti, F.; Arvan, P. INS-gene mutations: From genetics and beta cell biology to clinical disease. Mol. Asp. Med. 2015, 42, 3–18. [Google Scholar] [CrossRef] [PubMed]
  9. Balboa, D.; Saarimäki-Vire, J.; Borshagovski, D.; Survila, M.; Lindholm, P.; Galli, E.; Eurola, S.; Ustinov, J.; Grym, H.; Huopio, H.; et al. Insulin mutations impair beta-cell development in a patient-derived iPSC model of neonatal diabetes. eLife 2018, 7, e38519. [Google Scholar] [CrossRef] [PubMed]
  10. Haris, B.; Mohammed, I.; Al-Khawaga, S.; Hussain, K. Homozygous Insulin Promotor Gene Mutation Causing Permanent Neonatal Diabetes Mellitus and Childhood Onset Autoantibody Negative Diabetes in the Same Family. Int. Med. Case Rep. J. 2022, 15, 35–41. [Google Scholar] [CrossRef] [PubMed]
  11. Tikhonovich, Y.V.; Petryaykina, E.E.; Timofeev, A.V.; Zubkova, N.A.; Kolodkina, A.A.; Sorkina, E.L.; Vasiliev, E.V.; Petrov, V.M.; Andrianova, E.A.; Zilberman, L.I.; et al. Clinical, hormonal, and molecular-genetic characteristics of monogenic diabetes mellitus associated with the mutations in the INS gene. Diabetes Mellit. 2021, 24, 414–421. [Google Scholar] [CrossRef]
  12. Can, N.T.; Vu, D.C.; Bui, T.P.; Nguyen, K.N.; Nguyen, D.P.; Craig, M.; Ellard, S.; Nguyen, H.T. AB125. Neonatal diabetes mellitus due to insulin gene mutation. Ann. Transl. Med. 2015, 3 (Suppl. S2), AB125. [Google Scholar] [CrossRef]
  13. Bee, Y.M.; Zhao, Y.; Ellard, S.; Hattersley, A.T.; Yap, F. Permanent neonatal diabetes in siblings with novel C109Y INS mutation transmitted by an unaffected parent with somatic mosaicism. Pediatr. Diabetes 2014, 15, 324–328. [Google Scholar] [CrossRef] [PubMed]
  14. Al Shaikh, A.; Shirah, B.; Alzelaye, S. A homozygous mutation in the insulin gene (INS) causing autosomal recessive neonatal diabetes in Saudi families. Ann. Pediatr. Endocrinol. Metab. 2020, 25, 42. [Google Scholar] [CrossRef] [PubMed]
  15. Islam, M.J.; Parves, M.R.; Mahmud, S.; Tithi, F.A.; Reza, M.A. Assessment of structurally and functionally high-risk nsSNPs impacts on human bone morphogenetic protein receptor type IA (BMPR1A) by computational approach. Comput. Biol. Chem. 2019, 80, 31–45. [Google Scholar] [CrossRef] [PubMed]
  16. Sim, N.L.; Kumar, P.; Hu, J.; Henikoff, S.; Schneider, G.; Ng, P.C. SIFT web server: Predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012, 40, W452–W457. [Google Scholar] [CrossRef] [PubMed]
  17. Adzhubei, I.; Jordan, D.M.; Sunyaev, S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 2013, 76, 7–20. [Google Scholar] [CrossRef] [PubMed]
  18. Hecht, M.; Bromberg, Y.; Rost, B. Better prediction of functional effects for sequence variants. BMC Genom. 2015, 16, S1. [Google Scholar] [CrossRef] [PubMed]
  19. Manfredi, M.; Savojardo, C.; Martelli, P.L.; Casadio, R. E-SNPs&GO: Embedding of protein sequence and function improves the annotation of human pathogenic variants. Bioinformatics 2022, 38, 5168–5174. [Google Scholar] [PubMed]
  20. 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–2734. [Google Scholar] [CrossRef] [PubMed]
  21. Pejaver, V.; Urresti, J.; Lugo-Martinez, J.; Pagel, K.A.; Lin, G.N.; Nam, H.J.; Mort, M.; Cooper, D.N.; Sebat, J.; Iakoucheva, L.M.; et al. Inferring the molecular and phenotypic impact of amino acid variants with MutPred2. Nat. Commun. 2020, 11, 5918. [Google Scholar] [CrossRef] [PubMed]
  22. Stenson, P.D.; Mort, M.; Ball, E.V.; Evans, K.; Hayden, M.; Heywood, S.; Hussain, M.; Phillips, A.D.; Cooper, D.N. The Human Gene Mutation Database: Towards a comprehensive repository of inherited mutation data for medical research, genetic diagnosis, and next-generation sequencing studies. Hum. Genet. 2017, 136, 665–677. [Google Scholar] [CrossRef] [PubMed]
  23. Bava, K.A.; Gromiha, M.M.; Uedaira, H.; Kitajima, K.; Sarai, A. ProTherm, version 4.0: Thermodynamic database for proteins and mutants. Nucleic Acids Res. 2004, 32, D120–D121. [Google Scholar] [CrossRef] [PubMed]
  24. Cheng, J.; Randall, A.; Baldi, P. Prediction of protein stability changes for single-site mutations using support vector machines. Proteins Struct. Funct. Bioinform. 2006, 62, 1125–1132. [Google Scholar] [CrossRef] [PubMed]
  25. Venselaar, H.; Te Beek, T.A.; Kuipers, R.K.; Hekkelman, M.L.; Vriend, G. Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces. BMC Bioinform. 2010, 11, 548. [Google Scholar] [CrossRef] [PubMed]
  26. Khetan, S.; Kursawe, R.; Youn, A.; Lawlor, N.; Jillette, A.; Marquez, E.J.; Ucar, D.; Stitzel, M.L. Type 2 Diabetes-Associated Genetic Variants Regulate Chromatin Accessibility in Human Islets. Diabetes 2018, 67, 2466–2477. [Google Scholar] [CrossRef] [PubMed]
  27. Kumar, P.; Henikoff, S.; Ng, P.C. A comprehensive evaluation of the SIFT and PolyPhen2 algorithms for predicting the impact of amino acid substitutions on protein function. Bioinformatics 2009, 25, 2843–2844. [Google Scholar]
  28. Ioannidis, N.M.; Rothstein, J.H.; Pejaver, V.; Middha, S.; McDonnell, S.K.; Baheti, S.; Musolf, A.; Li, Q.; Holzinger, E.; Karyadi, D.; et al. REVEL: An ensemble method for predicting the pathogenicity of rare missense variants. Am. J. Hum. Genet. 2016, 99, 877–885. [Google Scholar] [CrossRef] [PubMed]
  29. Martelotto, L.G.; Ng, C.K.; Piscuoglio, S.; Weigelt, B.; Reis-Filho, J.S. Breast cancer intra-tumor heterogeneity. Breast Cancer Res. 2014, 16, 210. [Google Scholar] [CrossRef] [PubMed]
  30. Bendl, J.; Stourac, J.; Salanda, O.; Pavelka, A.; Wieben, E.D.; Zendulka, J.; Brezovsky, J.; Damborsky, J. PredictSNP: Robust and accurate consensus classifier for prediction of disease-related mutations. PLoS Comput. Biol. 2014, 10, e1003440. [Google Scholar] [CrossRef] [PubMed]
  31. Khan, S.; Vihinen, M. Performance of protein stability predictors. Hum. Mutat. 2010, 31, 675–684. [Google Scholar] [CrossRef] [PubMed]
  32. Laurenzano, S.E.; McFall, C.; Nguyen, L.; Savla, D.; Coufal, N.G.; Wright, M.S.; Tokita, M.; Dimmock, D.; Kingsmore, S.F.; Newfield, R.S. Neonatal diabetes mellitus due to a novel variant in the INS gene. Mol. Case Stud. 2019, 5, a004085. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flowchart to identify and categorize nsSNPs in the INS gene; each step indicates the tools used. If one nsSNP is considered deleterious in each step using certain tool it will be directly analyzed in the following tool or in the next step using another computational tool. In steps 8 and 9, tools are used to analyze and display structural changes.
Figure 1. Flowchart to identify and categorize nsSNPs in the INS gene; each step indicates the tools used. If one nsSNP is considered deleterious in each step using certain tool it will be directly analyzed in the following tool or in the next step using another computational tool. In steps 8 and 9, tools are used to analyze and display structural changes.
Jpm 14 00425 g001
Figure 2. A Root-Mean-Square Deviation (RMSD) analysis of proteins throughout molecular dynamics (MD) simulations.
Figure 2. A Root-Mean-Square Deviation (RMSD) analysis of proteins throughout molecular dynamics (MD) simulations.
Jpm 14 00425 g002
Figure 3. The Root-Mean-Square fluctuation (RMSF) of individual atoms in a protein throughout the MD simulation.
Figure 3. The Root-Mean-Square fluctuation (RMSF) of individual atoms in a protein throughout the MD simulation.
Jpm 14 00425 g003
Figure 4. The radius of gyration of a wild-type protein and its three variations (C96R, C96Y, and P52R) as a function of time throughout a simulation duration of 1000 ps.
Figure 4. The radius of gyration of a wild-type protein and its three variations (C96R, C96Y, and P52R) as a function of time throughout a simulation duration of 1000 ps.
Jpm 14 00425 g004
Figure 5. The Solvent Accessible Surface Area (SASA) of a wild-type protein and the three variants (C96R, C96Y, and P52R) during a 1000 ps molecular dynamics simulation.
Figure 5. The Solvent Accessible Surface Area (SASA) of a wild-type protein and the three variants (C96R, C96Y, and P52R) during a 1000 ps molecular dynamics simulation.
Jpm 14 00425 g005
Figure 6. Quantification of the hydrogen bonds present in a wild-type protein and its three variations (C96R, C96Y, and P52R) over a 1000 ps molecular dynamics simulation.
Figure 6. Quantification of the hydrogen bonds present in a wild-type protein and its three variations (C96R, C96Y, and P52R) over a 1000 ps molecular dynamics simulation.
Jpm 14 00425 g006
Figure 7. Three-dimensional structure of the wild-type INS protein (PDB ID: 4i30), and its mutations, P52R, C96R, and C96Y. The sites of the mutations are shown in the surface presentation and red color. This figure was prepared using PyMol Version 2.0.
Figure 7. Three-dimensional structure of the wild-type INS protein (PDB ID: 4i30), and its mutations, P52R, C96R, and C96Y. The sites of the mutations are shown in the surface presentation and red color. This figure was prepared using PyMol Version 2.0.
Jpm 14 00425 g007
Figure 8. The three-dimensional structure of the wild-type INS protein (PDB ID: 3i40). This figure was prepared using PyMol software.
Figure 8. The three-dimensional structure of the wild-type INS protein (PDB ID: 3i40). This figure was prepared using PyMol software.
Jpm 14 00425 g008
Figure 9. The three-dimensional structure of the INS protein (PDB ID: 3i40) for the mutation P52R; the site of the mutation is red. This figure was prepared using PyMol software.
Figure 9. The three-dimensional structure of the INS protein (PDB ID: 3i40) for the mutation P52R; the site of the mutation is red. This figure was prepared using PyMol software.
Jpm 14 00425 g009
Figure 10. The three-dimensional structure of the INS protein (PDB ID: 3i40) for the mutation C96R; the site of the mutation is red. This figure was prepared using PyMol software.
Figure 10. The three-dimensional structure of the INS protein (PDB ID: 3i40) for the mutation C96R; the site of the mutation is red. This figure was prepared using PyMol software.
Jpm 14 00425 g010
Figure 11. The three-dimensional structure of the INS protein (PDB ID: 3i40) for the mutation C96Y; the site of the mutation is red. This figure was prepared using PyMol software.
Figure 11. The three-dimensional structure of the INS protein (PDB ID: 3i40) for the mutation C96Y; the site of the mutation is red. This figure was prepared using PyMol software.
Jpm 14 00425 g011
Table 1. Functional analysis results using SIFT, PolyPhen-2, and SNAP2 tools. Transcription ID is ENST00000381330.5.
Table 1. Functional analysis results using SIFT, PolyPhen-2, and SNAP2 tools. Transcription ID is ENST00000381330.5.
NCBI DATABASESIFTPolyPhen-2SNAP2
Variant IDAlleleMAFMutationPredictionScorePredictionScorePredictionAccuracy %
rs28933985C|GT = 0./0R89Pdeleterious0.00Probably1.0074/8585%
rs28933985C|TT = 0./0R89Hdeleterious0.00Probably1.0068/8080%
rs80356664C|GT = 0./0G32Rdeleterious0.00Probably1.0095/9595%
rs80356664C|TT = 0./0G32Sdeleterious0.00Probably1.0088//9191%
rs80356666A|CC = 0./0C43Gdeleterious0.00Probably1.0086/9191%
rs80356667C|AA = 0./0G47Vdeleterious0.00Probably1.0085/9191%
rs80356668A|CC = 0./0F48Cdeleterious0.00Probably1.0081/9191%
rs80356668A|GC = 0./0F48Sdeleterious0.00Probably1.0089/9191%
rs80356669G|ANAR89Cdeleterious0.00Probably1.0062/8080%
rs80356670C|AA = 0./0G90Cdeleterious0.00Probably1.0065//8080%
rs80356671C|GT = 0./0C96Sdeleterious0.01Probably0.99581/9191%
rs80356671C|TT = 0./0C96Ydeleterious0.00Probably1.0087/9191%
rs80356672T|CC = 0./0Y108Cdeleterious0.00Probably1.0078/8585%
rs121908261G|AA = 0./0R55Cdeleterious0.00Probably1.0061/8080%
rs121908272G|CNAH29Ddeleterious0.01Probably1.0078/8585%
rs121908273A|GT = 0./0L35Pdeleterious0.00Probably1.0096/9595%
rs121908273A|TT = 0./0L35Qdeleterious0.00Probably1.0091/9595%
rs121908276G|CNAS101Cdeleterious0.00Probably1.0054/7575%
rs121908277T|CNAY103Cdeleterious0.02Probably1.0055/7575%
rs121918102C|ANAV92Ldeleterious0.00Probably1.0077/8585%
rs145038693G|AC = 0./0P52Ldeleterious0.00Probably1.0073/8585%
rs145038693G|CC = 0./0P52Rdeleterious0.00Probably0.99877/8585%
rs145038693G|TC = 0./0P52Hdeleterious0.00Probably1.0075/8585%
rs148685531G|CA = 0.000142/2F49Ldeleterious0.00Probably0.99481/9191%
rs397515521C|ANAM1Ideleterious0.00Probably0.98177/8585%
rs397515521C|TNAM1Ideleterious0.00Probably0.98177/8585%
rs765512575C|TT = 0.000105/2G44Rdeleterious0.00Probably0.99686/9191%
rs1057524907T|CNAE93Gdeleterious0.01Probably0.97977/8585%
rs1252051752T|ANAS98Cdeleterious0.04Probably0.99637/6666%
rs1278232284G|AA = 0.000004/1L35Vdeleterious0.00Probably1.0075/8585%
rs1278232284G|CA = 0.000004/1L35Mdeleterious0.00Probably1.0077/8585%
rs1564912274T|CC = 0.000004/1H34Rdeleterious0.00Probably0.99764/8080%
rs1564912274T|GC = 0.000004/1H34Pdeleterious0.00Probably0.99982/9191%
rs1845839718A|GNAC96Rdeleterious0.00Probably1.0090/9595%
rs2133672742C|ANAC109Fdeleterious0.00Probably1.0083//9191%
rs2133672778A|CNAY108Ddeleterious0.00Probably1.0090/9595%
rs2133676660G|ANAL39Fdeleterious0.00Probably1.0072/8585%
rs2133676747C|ANAG32Vdeleterious0.00Probably1.0088/9191%
rs2133676771A|CNAC31Gdeleterious0.00Probably1.0090//9595%
rs28933985C|TT = 0./0R89Hdeleterious0.00Probably1.0068/8080%
rs80356664C|GT = 0./0G32Rdeleterious0.00Probably1.0095/9595%
rs80356664C|TT = 0./0G32Sdeleterious0.00Probably1.0088/9191%
rs80356666A|CC = 0./0C43Gdeleterious0.00Probably1.0086//9191%
rs80356667C|AA = 0./0G47Vdeleterious0.00Probably1.0085/9191%
rs80356668A|CC = 0./0F48Cdeleterious0.00Probably1.0078/8585%
rs80356668A|GC = 0./0F48Sdeleterious0.00Probably1.0077/8585%
rs80356669G|ANAR89Cdeleterious0.00Probably1.0062/8080%
rs80356670C|AA = 0./0G90Cdeleterious0.00Probably1.0065/8080%
rs80356671C|GT = 0./0C96Sdeleterious0.01Probably0.99581/9191%
rs80356672T|CC = 0./0Y108Cdeleterious0.00Probably1.0078/8585%
rs11557614G|AA = 0./0A38Vdeleterious0.01Probably0.99170/8585%
rs144093133C|AT = 0.000061/6G75Vdeleterious0.00Probably0.99938/6666%
rs202244834T|GG = 0./0K53Tdeleterious0.02Probably0.96354/7575%
rs760425445C|TT = 0.000004/1G90Ddeleterious0.00Probably0.98780/9191%
rs781016664G|AA = 0.000012/3R56Wdeleterious0.00Probably1.0076/8585%
rs983508038C|TT = 0.000004/1R56Qdeleterious0.00Probably1.0057/7575%
rs1182567488T|GG = 0.000071/1Q28Hdeleterious0.00Probably0.99558/7575%
rs1184417816C|TT = 0.000111/1R55Hdeleterious0.02Probably1.0060/8080%
rs1213888316A|GG = 0./0L11Pdeleterious0.00Probably0.99461/8080%
rs1313322794C|GG = 0.000004/1G47Rdeleterious0.00Probably1.0091/9595%
rs1313490068A|GG = 0.000224/1L13Pdeleterious0.00Probably0.99783/9191%
rs1460766978G|AA = 0./0P52Sdeleterious0.00Probably1.0064/8080%
rs1460766978G|TA = 0./0P52Tdeleterious0.00Probably1.0067/8080%
rs1845838687T|CC = 0./0Q104Rdeleterious0.02Probably0.97659/7575%
rs1845873763C|TT = 0.00005/1E59Kdeleterious0.02Probably0.98635/6666%
rs1845879949A|TT = 0./0L14Qdeleterious0.00Probably1.0061/8080%
rs2133677126A|TT = 0.00007/2M1Kdeleterious0.00Probably0.99594/9595%
NA = not available.
Table 2. Pathogenicity analysis results using SNP & GO and PhD-SNP tools.
Table 2. Pathogenicity analysis results using SNP & GO and PhD-SNP tools.
NCBI DATABASESNPs & GOPhD-SNP.
Variant IDAlleleMutationPredictionSCOREPredictionSCORE
rs28933985C|GR89PDISESASE10DISESASE2
rs80356671C|TC96YDISESASE10DISESASE1
rs121908273A|GL35PDISESASE10DISESASE3
rs145038693G|AP52LDISESASE10DISESASE1
rs145038693G|CP52RDISESASE10DISESASE2
rs145038693G|TP52HDISESASE9DISESASE1
rs765512575C|TG44RDISESASE10DISESASE0
rs1564912274T|CH34RDISESASE10DISESASE2
rs1564912274T|GH34PDISESASE10DISESASE1
rs1845839718A|GC96RDISESASE10DISESASE5
rs1460766978G|AP52SDISESASE10DISESASE0
rs1845879949A|TL14QDISESASE9DISESASE5
Table 3. Pathogenicity analysis results using MutPred2 tools.
Table 3. Pathogenicity analysis results using MutPred2 tools.
NCBI DATABASEMUTPRED2
Variant IDAlleleMutationEFFECTSCORE
rs28933985C|GR89PAffect0.791
rs80356671C|TC96YAffect0.860
rs121908273A|GL35PAffect0.925
rs145038693G|CP52RAffect0.570
rs765512575C|TG44RAffect0.716
rs1564912274T|GH34PAffect0.663
rs1845839718A|GC96RAffect0.910
rs1845879949A|TL14QAffect0.627
Table 4. Results of the impact of SNPs on protein stability analysis using I-MUTANT and MuPro tool.
Table 4. Results of the impact of SNPs on protein stability analysis using I-MUTANT and MuPro tool.
NCBI DATABASEI-MUTANTMuPro
Variant IDAlleleMutationStabilityReliability IndexEFFECTConfidence Score
rs28933985C|GR89PDECREASE5Increase0.078515012
rs80356671C|TC96YDECREASE4Decrease−0.080856117
rs121908273A|GL35PDECREASE4Increase0.065423903
rs145038693G|CP52RDECREASE0Decrease−1.0
rs765512575C|TG44RDECREASE1Increase0.49104088
rs1845839718A|GC96RDECREASE8Decrease−0.574683230494305
rs1845879949A|TL14QDECREASE8Decrease−0.62610731178958
Table 5. Summary of HOPE results for the different amino acid mutations of the INS gene.
Table 5. Summary of HOPE results for the different amino acid mutations of the INS gene.
Variant IDrs145038693rs80356671rs1845839718rs1845879949
MutationsP52RC96YC96RL14Q
Proline into Arginine at position 52Cysteine into Tyrosine at position 96Cysteine into Arginine at position 96Leucine into Glutamine at position 14
AlleleG|CC|TA|GA|T
Amino Acids Properties
SizeThe mutant residue is bigger than the wild-type residue.The mutant residue is bigger than the wild-type residue.The mutant residue is bigger than the wild-type residueThe mutant residue is bigger than the wild-type residue.
Hydrophobicity valueWild-type residue is more hydrophobic than the mutant residue.Wild-type residue is more hydrophobic than the mutant residue.Wild-type residue is more hydrophobic than the mutant residueWild-type residue is more hydrophobic than the mutant residue.
ChargeWild-type residue charge was NEUTRALNo changeWild-type residue charge was NEUTRALNo change
Mutant residue charge was POSITIVE---Mutant residue charge was POSITIVE.----
Variant’s score0.903225660.99475780.99398980.8815106
ConservationMutant and wild-type residues were not very similarWild-type residue was very conservedWild-type residue was very conservedWild-type residue was very conserved
The mutant residue was located near a highly conserved positionThe mutant residue was located near a highly conserved positionMutant residue was located near a highly conserved position
Effect of mutation on proteinProbably damaging to the protein.Probably damaging to the protein.Probably damaging to the protein.Mutation occurs without damaging the protein
--- = No change (mutant residue charge); ---- = blank (conservation).
Table 6. Comparison of the results of our study with the results of other published studies on INS gene mutations associated with PNDM.
Table 6. Comparison of the results of our study with the results of other published studies on INS gene mutations associated with PNDM.
Studies on INS Gene Mutations Associated with PNDM
This StudyLaurenzano et al. 2019
[32]
Balboa et al. 2018
[9]
Liu et al. 2015
[8]
MutationsP52RP9RC96RL13R
C96R C109Y
C96Y
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ahmed, E.M.; Elangeeb, M.E.; Adam, K.M.; Abuagla, H.A.; MohamedAhmed, A.A.E.; Ali, E.W.; Eltieb, E.I.; Edris, A.M.; Ali Osman, H.M.; Idris, E.S.; et al. Computational Analysis of Deleterious nsSNPs in INS Gene Associated with Permanent Neonatal Diabetes Mellitus. J. Pers. Med. 2024, 14, 425. https://doi.org/10.3390/jpm14040425

AMA Style

Ahmed EM, Elangeeb ME, Adam KM, Abuagla HA, MohamedAhmed AAE, Ali EW, Eltieb EI, Edris AM, Ali Osman HM, Idris ES, et al. Computational Analysis of Deleterious nsSNPs in INS Gene Associated with Permanent Neonatal Diabetes Mellitus. Journal of Personalized Medicine. 2024; 14(4):425. https://doi.org/10.3390/jpm14040425

Chicago/Turabian Style

Ahmed, Elsadig Mohamed, Mohamed E. Elangeeb, Khalid Mohamed Adam, Hytham Ahmed Abuagla, Abubakr Ali Elamin MohamedAhmed, Elshazali Widaa Ali, Elmoiz Idris Eltieb, Ali M. Edris, Hiba Mahgoub Ali Osman, Ebtehal Saleh Idris, and et al. 2024. "Computational Analysis of Deleterious nsSNPs in INS Gene Associated with Permanent Neonatal Diabetes Mellitus" Journal of Personalized Medicine 14, no. 4: 425. https://doi.org/10.3390/jpm14040425

APA Style

Ahmed, E. M., Elangeeb, M. E., Adam, K. M., Abuagla, H. A., MohamedAhmed, A. A. E., Ali, E. W., Eltieb, E. I., Edris, A. M., Ali Osman, H. M., Idris, E. S., & Khalil, K. A. A. (2024). Computational Analysis of Deleterious nsSNPs in INS Gene Associated with Permanent Neonatal Diabetes Mellitus. Journal of Personalized Medicine, 14(4), 425. https://doi.org/10.3390/jpm14040425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop