Next Article in Journal
Administration of Tamoxifen Can Regulate Changes in Gene Expression during the Acute Phase of Traumatic Spinal Cord Injury
Previous Article in Journal
MicroRNAs as Potential Regulators of GSK-3β in Renal Cell Carcinoma
Previous Article in Special Issue
A Comparative Study on the Carcass and Meat Chemical Composition, and Lipid-Metabolism-Related Gene Expression in Korean Hanwoo and Brindle Chikso Cattle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Investigation of AKT2 Gene and Protein Abnormalities Reveals Potential Association with Insulin Resistance and Type 2 Diabetes

by
M. E. Elangeeb
1,*,
Imadeldin Elfaki
2,
M. A. Elkhalifa
3,
Khalid M. Adam
4,
A. O. Alameen
5,
Ahmed Kamaleldin Elfadl
6,7,
Ibrahim Altedlawi Albalawi
8,
Kholoud S. Almasoudi
9,
Reema Almotairi
9,
Basim S. O. Alsaedi
10,
Marwan H. Alhelali
10,
Mohammad Muzaffar Mir
11,
Dnyanesh Amle
12 and
Rashid Mir
9,*
1
Department of Basic Medical Sciences, College of Applied Medical Sciences, University of Bisha, Bisha 61922, Saudi Arabia
2
Department of Biochemistry, Faculty of Science, University of Tabuk, Tabuk 47512, Saudi Arabia
3
Department of Anatomy, Faculty of Medicine and Health Sciences, University of Bisha, Bisha 61922, Saudi Arabia
4
Department of Medical Laboratory Sciences, College of Applied Medical Sciences, University of Bisha, Bisha 61922, Saudi Arabia
5
Department of Biomedical Science, Faculty of Veterinary Medicine, King Faisal University, Alahssa 31982, Saudi Arabia
6
Veterinary Research Section, Ministry of Municipality, Doha P.O. Box 35081, Qatar
7
Department of Pathology, Faculty of Veterinary Medicine, University of Khartoum, Khartoum 11115, Sudan
8
Department of Surgical Oncology, Faculty of Medicine, University of Tabuk, Tabuk 47512, Saudi Arabia
9
Department of Medical Lab Technology, Prince Fahad Bin Sultan Chair for Biomedical Research, Faculty of Applied Medical Sciences, University of Tabuk, Tabuk 71491, Saudi Arabia
10
Department of Statistics, University of Tabuk, Tabuk 47512, Saudi Arabia
11
Department of Basic Medical Sciences, College of Medicine, University of Bisha, Bisha 61922, Saudi Arabia
12
Department of Biochemistry, All India Institute of Medical Sciences, Nagpur 441108, India
*
Authors to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2023, 45(9), 7449-7475; https://doi.org/10.3390/cimb45090471
Submission received: 18 July 2023 / Revised: 12 August 2023 / Accepted: 24 August 2023 / Published: 12 September 2023
(This article belongs to the Special Issue Understanding Adipose Tissue: A Molecular Perspective)

Abstract

:
Type 2 diabetes (T2D) develops from insulin resistance (IR) and the dysfunction of pancreatic beta cells. The AKT2 protein is very important for the protein signaling pathway, and the non-synonymous SNP (nsSNPs) in AKT2 gene may be associated with T2D. nsSNPs can result in alterations in protein stability, enzymatic activity, or binding specificity. The objective of this study was to investigate the effect of nsSNPs on the AKT2 protein structure and function that may result in the induction of IR and T2D. The study identified 20 variants that were considered to be the most deleterious based on a range of analytical tools included (SIFT, PolyPhen2, Mut-pred, SNAP2, PANTHER, PhD-SNP, SNP&Go, MUpro, Cosurf, and I-Mut). Two mutations, p.A179T and p.L183Q, were selected for further investigation based on their location within the protein as determined by PyMol. The results indicated that mutations, p.A179T and p.L183Q alter the protein stability and functional characteristics, which could potentially affect its function. In order to conduct a more in-depth analysis of these effects, a molecular dynamics simulation was performed for wildtype AKT2 and the two mutants (p.A179T and p.L183Q). The simulation evaluated various parameters, including temperature, pressure, density, RMSD, RMSF, SASA, and Region, over a period of 100 ps. According to the simulation results, the wildtype AKT2 protein demonstrated higher stability in comparison to the mutant variants. The mutations p.A179T and p.L183Q were found to cause a reduction in both protein stability and functionality. These findings underscore the significance of the effects of nsSNPs (mutations p.A179T and p.L183Q) on the structure and function of AKT2 that may lead to IR and T2D. Nevertheless, they require further verifications in future protein functional, protein–protein interaction, and large-scale case–control studies. When verified, these results will help in the identification and stratification of individuals who are at risk of IR and T2D for the purpose of prevention and treatment.

1. Introduction

Diabetes mellitus (DM) is a metabolic disease characterized by elevated blood sugar. The general types of DM include type 1 DM, which is caused by the absence of insulin due to pancreatic beta cell destruction mediated by autoimmune mechanisms [1]. Type 2 diabetes mellitus (T2D) is characterized by insulin resistance (IR) and pancreatic beta cell dysfunction [2]. The third type of DM is called gestational DM (GDM), which develops during pregnancy (in the late second trimester or early in the third trimester) [3]. GDM is caused by any degree of glucose intolerance and defective pancreatic beta cells [3,4]. In addition, there are other types of DM such as maturity-onset diabetes of the young (MODY) and neonatal DM [4]. T2D is induced by the interaction of genetic and environmental or lifestyle risk factors [5]. Lifestyle or environmental risk factors include poor diet, reduced physical activity, depression, stress, and smoking [5]. In the last decade, genome-wide association studies (GWAs) have revealed the association of specific loci with diseases such as cancers, T2D, and cardiovascular disease [6,7,8,9,10,11]. For example, loci involved in the insulin signaling pathway may lead to IR and T2D [12,13]. The AKT protein is composed of 480 amino acid residues and comprises three domains: the pleckstrin homology (PH), middle catalytic (kinase) domain, and regulatory C-terminal disordered tail (C-tail) [14,15]. The PH domain anchors AKT to the membrane [15]. AKTs are classified into three isoforms (AKT1, AKT2, and AKT3) based on the differences in serine/threonine amino acid residues [15]. AKT1 is found ubiquitously, and AKT2 is found in insulin-sensitive tissues (liver, muscles, and adipose tissues) [15], whereas AKT3 is expressed in the testes and brain [15].
The serine/threonine kinase AKT2 is an effector protein in the insulin-signaling pathway and is connected to insulin metabolic actions [16]. It is encoded by the AKT2 gene (located on chromosome 19q13.1). The phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that phosphorylate phosphatidylinositol, which is a component of eukaryotic cell membranes [17]. The AKT2 protein is expressed in various cell types within the human body, such as adipocytes, skeletal muscle cells, hepatocytes, and the pancreas [15,16]. The expression of AKT2 is stimulated by different stimuli such as cell stress, movement, hormones, and drugs that influence cellular metabolism [16]. There are two signaling pathways stimulated by insulin: the PI3K regulates metabolic function, and the mitogen-activated protein kinase (MAPK) signaling pathway is implicated in cell survival and mitogenesis [18].
PI3K/AKT signaling is crucial for regulating insulin responsiveness and other physiological and pathological processes, including growth, glucose homeostasis, lipid metabolism, protein synthesis, cell proliferation, survival, and angiogenesis [7,15,19,20,21]. In muscles, insulin regulates metabolism through the promotion of sugar uptake, glycogenesis, and the synthesis of protein via the PI3K/AKT signaling pathway [15]. It has been reported that about 90% of glucose metabolism is stimulated by insulin in muscles, which has important roles in energy homeostasis [15]. In adipose tissues, insulin–AKT signaling enhances the utilization of sugar and the synthesis of proteins and lipids [15,22]. In hepatocytes during the fed state, the PI3K/AKT signaling pathway decreases gluconeogenesis and glycogen breakdown and enhances glycogenesis and lipogenesis for peripheral utilization and storage [15,23].
IR and abnormalities in lipid metabolism have been associated with the dysfunction of AKT2 and other elements of the insulin signaling pathway [16,24]. Gene variations (p.G972R, and p.T608R) of insulin receptor substrate 1 (IRS-1), PI3K and (R274H, R208 K and R467W) and in AKT2 were reported to be associated with T2D, fasting hyperglycemia and postprandial hyperinsulinemia [24]. In addition, the p.Pro50Thr AKT2 mutation was associated with T2D risk in the Finnish population [25]. Moreover, gene polymorphisms in AKT1 and AKT2 were associated with T2D, polycystic ovarian syndrome (PCOS) and cancer in the Chinese population [26,27,28]. Furthermore, mutation (p.E17K) in AKT2 gene can result in defective insulin signaling pathway, which causes reduced levels of plasma insulin, ketone bodies, fatty acids and glucose [29]. It has been suggested that autonomous activation of the downstream insulin signaling pathway as a treatment strategy for patients carrying the mutation (p.E17K) should be considered [29], since the ATK2 protein has crucial functions in the insulin signaling pathway and defective insulin signaling pathway leads to metabolic syndrome, IR and T2D [30,31,32]. In the present study, we used bioinformatics tools to examine the non-synonymous SNP (nsSNPs) in the AKT2 protein that may influence its structure and function.

2. Materials and Methods

Figure 1 Schematic illustration of the steps of this study.

2.1. Dataset Download

The FASTA-formatted nucleotide sequence (Accession number: NG 012038.2) and amino acid sequence (NP 001229956.1) of the AKT2 protein were obtained from the National Center for Biotechnology Information (NCBI) database, which can be accessed at http://www.ncbi.nlm.nih.gov (accessed on 23 January 2023). The single nucleotide polymorphisms (SNPs) databases for the AKT2 (GLU17LYS) gene can be found in the NCBI’s SNP database (SNP), which is accessible at http://www.ncbi.nlm.nih.gov/snp/ (accessed on 23 January 2023). Similar data were acquired from the Online Mendelian Inheritance in Man (OMIM) database (http://www.omim.org (accessed on 23 January 2023)) regarding the AKT2 gene and protein.

2.2. Predicting the Effects of nsSNPs on AKT2 Protein Function

The assessment of potential effects of SNPs involved the utilization of several online tools and servers, namely the Sorting Intolerant From Tolerant (SIFT) Sequence, Polymorphism Phenotyping (poly-Phen2), Single Amino Acid Predication (SNAP2), and Protein Analysis Through Evolutionary Relationships (PANTHER). These tools and servers were employed to evaluate the impact of SNPs. Four servers were utilized to classify nsSNPs and assess their confidence ratings based on multiple inputs. The impact on protein function, specifically for non-synonymous single amino acid substitutions, can be determined by analyzing the predictions provided by SIFT, PolyPhen2, SNAP2, and PANTHER. These predictions classify the impact as either tolerable or damaging. The SIFT algorithm employs a tolerance index score to evaluate SNPs, whereby variations with scores below 0.05 are deemed deleterious. In contrast, the PolyPhen algorithm assigns a numerical score ranging from 0 to 1, where a score of 0 indicates a neutral effect and a score of 1 indicates the most severe negative effect. Another tool utilized in this study was SNAP2, which facilitates the comparison of genomes and provides predictions regarding functional consequences at the amino acid level. The SNAP2 algorithm utilizes a protein sequence and a set of Single Amino Acid Variants (SAVs) as input in order to forecast the impact of each substitution on the molecular function of the protein. The prediction scores span from a value of −100, indicating complete neutrality, to +100, indicating a significant effect. The final tool, known as PANTHER, employs a scoring system that spans from −1 to 1 in order to assess the projected influence of the input single nucleotide polymorphisms (SNPs) on protein functionality. The functionality of a single nucleotide polymorphism (SNP) can be deduced by evaluating its score. A positive score indicates that the SNP is functional and has an effect on protein function, while a negative score suggests that the SNP is non-functional. Greater absolute values are indicative of a more pronounced predicted impact on protein function. The single nucleotide polymorphisms (SNPs) of the AKT2 gene, which were found to be common among the fourth group of servers, were subsequently subjected to computational analysis.

2.3. The Prediction of Associations between nsSNPs and Disease

Predictor of human deleterious single nucleotide polymorphisms (PhD-SNP) (http://snps.biofold.org/phd-snp/phd-snp.html (accessed on 23 January 2023)) and SNPs&GO (http://snps-and-go.biocomp.unibo.it/snps-and-go/ (accessed on 23 January 2023)) were used for determining association of filtered SNPs with disease. The PhD-SNP online bioinformatics tool predicts the association between SNPs and diseases using machine learning methods. It assigns SNPs a score between 0 and 9 based on the probability of them resulting in a disease, classifying them as either disease-associated or neutral. The proportion of PhD accuracy SNP is 78% [33]. SNPs&GO is a precise tool that predicts disease-related amino acid changes at a particular position within a protein and offers functional classifications, with an overall prediction accuracy of 82% [34]. The UniProt accession number (P31751) of the AKT2 protein and the locations of the original and modified amino acids were the data needed by SNPs&GO.

2.4. Predicting the Impact of SNPs on Protein Stability

Single nucleotide polymorphisms (SNPs) can modulate protein strength by either reducing or enhancing protein stability. In order to predict these effects with maximum confidence, a pair of tools was used. The tool known as I-Mutant, accessible at (http://folding.biofold.org/i-mutant/i-mutant2.0.html (accessed on 23 January 2023)), was utilized to forecast the effect of SNPs on the alteration of a protein’s stable state. It has been reported to possess an accuracy of up to 77% [35]. The I-Mutant algorithm was provided with the amino acid sequence of the AKT2 protein, along with information regarding the specific residues that had undergone mutations, including their respective positions. MUpro (https://mupro.proteomics.ics.uci.edu/ (accessed on 23 January 2023)) is a set of machine learning algorithms that detect alterations in protein states and stability due to amino acid mutations [36]. The inputs utilized by MUpro were identical to those employed by I-Mutant. However, MUpro additionally considers the locations of substitutions, along with the initial and altered amino acid residues.

2.5. Prediction of SNP Pathogencicty and Function Alteration on AKT2 Protein Using Mut-Pred Server

The Mut-pred web server (http://mutpred.mutdb.org/#qform (accessed on 23 January 2023)) is a valuable resource for forecasting the prospective influence of missense mutations on the functionality of proteins. The Mut-pred tool utilizes a comprehensive set of fifty distinct protein properties in order to determine the impact of substitutions. The objective is to forecast the probability of a particular missense mutation being either deleterious or benign. The requisite data for Mut-pred comprises the amino acid sequence of the protein that encompasses the mutation, in conjunction with the identity and position of the mutated residue. The resultant output furnishes an anticipated pathogenicity score, along with particulars regarding the precise characteristics of the protein function that are projected to be impacted by the mutation. A MutPred2 score exceeding 0.5 indicates the presence of a pathogenic mutation.

Analyzing Protein Sequence Conservation Using ConSurf Server

The examination of the conservation of the AKT2 protein was conducted utilizing the web-based programme ConSurf, accessible at http://ConSurf.tau.ac.il/ (accessed on 23 January 2023). This tool enables the high-throughput prediction of functional characteristics pertaining to specific protein regions. A conservation analysis value is assigned to each protein residue on a scale ranging from 1 to 9. Residues with scores of 1–3 indicate variabilities, scores of 4–6 indicate average conservation, and scores of 7–9 indicate highly conserved regions. In order to utilize the application, it is necessary to input the protein sequence in FASTA format.

2.6. Determination of the Active Binding Site Using PyMol Software

The software PyMol Version 2.0. was utilized to locate the binding site where alterations in amino acid composition took place. The utilization of pyMOL entails the examination of residues located in the active site. The aforementioned software is a tool that is open-source in nature and is utilized for the purpose of molecular visualization. It has the ability to produce three-dimensional images of small molecules with a superior level of quality. Thus, two potential mutants (A179T, L183Q) were selected from the 20 nsSNPs deemed to be the most deleterious mutants in order to conduct further analysis through molecular dynamics (MD) simulation.

2.7. Performing Molecular Dynamics Simulations (MDS)

GROMACS version 2020.6 was used on a Google Collaboration pro notebook with high RAM and python 3.8 to run molecular dynamics (MD) calculations on both wildtype and mutant structures to look into the structural changes over time. For the first computation, the GROMACS-OPLS-AA force field was used. The systems were submerged in a cubic box filled with water molecules that had marginal radii of 1 nm (Å10), and the simulation box was neutralized by introducing 10 sodium Na+ ions using the GROMACS genion tool. The energy minimization was performed using the steepest descent algorithm with an energy step size of 0.01 and a maximum of 50000 iterations. To maintain a stable system, the Parrinello–Rahman (pcouple) pressure of 1 bar and a Berendsen temperature (tcouple) of 300 K were used. For temperature and pressure, the coupling constant was changed to 0.1 and 2.0 ps, respectively.
The electrostatic interactions in the system were computed using the Partial Mesh Ewald (PME) algorithm. The short-range cut-off for both Van der Waals (rvdw) and electrostatic (rcoulomb) interactions was set at 1.0 nm. The neighbor list (nstlist) was updated every 10 ps. All bond constraints, including heavy atom-H bonds, were maintained using the LINCS algorithm with a time step of 0.002 ps. The isothermal compressibility was set to 4.5 × 10−5. The system was equilibrated for 100 ps each in NPT (constant number of particles, pressure, and temperature) and NVT (constant number of particles, volume, and temperature) ensembles with the Berendsen temperature (tcouple) at 300 K and the Parrinello–Rahman (pcouple) pressure at 1 bar with coupling constants set to 0.1 ps and 2.0 ps for temperature and pressure, respectively. Subsequently, the native and mutant structures were subjected to 10 ns of molecular dynamics simulations, with trajectories recorded every 1 ps. The comparative analysis of structural deviations between the wildtype and mutant structures was conducted using various tools, such as g_rms, g_rmsf, g_sasa, g_Rg, and g_density, to analyze RMSD (root mean square deviation), RMSF (root mean square fluctuation), SASA (solvent accessible surface area), Rg (radius of gyration), pressure, temperature, and density plots. The application XMGRACE was used to create each plot.

3. Results

The present study examined the association of nsSNPs in the AKT2 gene and IR and T2D. Understanding the effects of on protein function is crucial for unraveling disease genetics and finding potential therapies. Computational methods utilize evolutionary conservation, protein structure, and function to predict nsSNP impacts on stability and structure. Investigation of the nsSNP AKT2’s effects offers insights into its role in IR and T2D development, aiding in understanding of the underlying disease mechanisms and identifying therapeutic targets. In addition, predicting nsSNP effects in the AKT2 gene can pinpoint high-risk individuals for tailored interventions and advance comprehension of genetic variation role in IR and T2D. Deciphering nsSNP impact on AKT2 could lead to personalized treatments targeting genetic variants, improving IR and glucose regulation for those at T2D risk. We present our results on AKT2 nsSNPs and their relevance to IR and T2D susceptibility, aiming to contribute to disease understanding and therapeutic development.

3.1. The Anticipation of the Characteristics of SNPs

A summary of the results of the four different web servers, SIFT, PolyPhen2, SNAP2 and PANTHER, regarding the effects of different SNPs on the AKT2 protein are shown in Table 1.
The aforementioned table presents data pertaining to the SIFT score predictions for individual amino acid substitutions, which are utilized in the assessment of the effect on protein functionality. The data presented in the Table 1 suggest that the majority of substitutions are anticipated to be harmful, as evidenced by their SIFT score of less than 0.05, which implies a significant probability of impacting protein functionality.
Nonetheless, certain substitutions are anticipated to be innocuous or slightly harmful, exhibiting elevated SIFT scores. In general, the SIFT score presents a valuable tool for forecasting the impact of SNPs on protein functionality. According to Table 1, the impact of the substitution is measured in terms of PolyPhen and SNAP2 scores, which are calculated by the respective programs. The PolyPhen score predicts the impact of the substitution as either “probably damaging”, “possibly damaging”, or “benign”. The SNAP2 score predicts the effect of the substitution as either “effect” or “no effect”.
To predict the SNPs’ nature with high confidence, SIFT, PolyPhen2, SNAP2, and PANTHER servers were employed concurrently. SIFT and PolyPhen2 received 533 missense nsSNPs in total. After consulting the databases it is integrated with, it offered each reported SNP an index score. A total of 142 (27%) of the nsSNPs were indicated by SIFT indexing to be harmful with a score of ≤0.05. Moreover, 99 out of 533 nsSNPs were reported to be highly harmful with a confidence level of 0.00.
A total of 99 of the nsSNPs were classified as “probably damaging” by the PolyPhen indexing, which is based on structural data and multiple sequence alignment, and 43 as “possibly damaging”. As a result, 99 (18.6%) of nsSNPs had negative PolyPhen predictions. Combining data from SIFT and PolyPhen allowed the removal of 142 nsSNPs in order to produce the most accurate results. After that, SNAP2 and PANTHER servers received these 142 nsSNPs for additional verification. The 99 nsSNPs were then sent to the SNAP2 server to determine their effect. The findings showed that 75 nsSNPs were found to affect the protein functionally. A summary of the results related the nsSNP functional effects is shown in Table 1.

3.2. Identifying Disease-Linked SNPs

We employed SNP&GO and PhD-SNP tools to identify disease-associated SNPs and predict the connection between nsSNPs in AKT2 and IR and T2D, exploring their impact on AKT2 structure and function and IR and T2D development. The limited set of polymorphisms that have undergone filtration and reduction were subsequently evaluated to determine their association with diseases. Subsequently, 75 nsSNPs with a high level of confidence were submitted to SNP&GO and PhD-SNP. Both individuals computed a numerical value to classify every nsSNP as either benign or pathogenic. A single nsSNP was deemed neutral if its prediction score was less than 0.5, whereas a score greater than 0.5 indicated that the SNP was disease-causing. Out of a total of 75 nsSNPs, 47 variants were determined to have a neutral effect, while the remaining 28 nsSNPs were identified as having a causal association with disease. To reduce the number of nsSNPs, a cut off value of p ≥ 0.80 was utilized for PhD-SNP scores, resulting in a total of 20 nsSNPs. Table 2 depicts the outcomes that were subjected to filtration.

3.3. The Prognostication of the Impact of Protein Stability

We used the I-Mutant and MuPro servers to try to figure out what would happen to the structure of 20 possible nsSNPs. Based on the results, 20 nsSNPs were thought to stop proteins from working by making them less stable. Moreover, as per the Mut-pro server, it was indicated that only 3 nsSNP out of 20 showed an increase on the protein’s functionality due to structure impact. The findings are presented in Table 3.

3.4. Prediction of the Pathogenicity and Underlying Functional Alterations of nsSNP Effects on AKT2 Protein

We utilized The Mut-pred web server to predict the pathogenicity and underlying functional changes caused by nsSNPs on the AKT2 protein. This approach aids in understanding the potential impact of these genetic variations on protein function and their association with disease, such as T2D related to insulin resistance. The Mut-pred web server offers valuable insights into the molecular mechanisms linking nsSNPs in AKT2 to disease susceptibility and provides a foundation for further investigation. The Mut-pred web server generates outcomes that aid in predicting the probable influence of missense mutations on protein function that is linked with diseases. Furthermore, Mut-pred provides details regarding the particular characteristics of the protein that are anticipated to be impacted by the mutation, encompassing solvent accessibility, secondary structure, and sequence conservation. According to the results presented in Table 4 of Mut-pred, all 20 nsSNP were found to have an impact on protein function. However, it was observed that 10 out of the 20 nsSNPs had a solely negative effect on the AKT2 protein, attributable to alterations in the protein’s structure and characteristics.

3.5. Study of Sequence Conservation

The ConSurf server was utilized to investigate the preservation of 20 nsSNPs that are associated with T2D. As per the findings, a set of three nsSNPs exhibited variability, whereas four nsSNPs demonstrated a modest conservation level. A set of nsSNPs that were determined to be conserved was analyzed, revealing that seven of them displayed both structural and spatial conservation, and they were located in buried regions. Six nsSNPs were predicted to be functionally conserved and exposed. The ConSurf tool employs chromatic schemes to represent the extent of sequence conservation and offers predictions about the nature and location of modifications on a numerical scale ranging from 1 to 9. Table 5 and Figure 2 provide a comprehensive summary of the 20 nsSNPs that were detected, including their functional and structural implications, as well as their phylogenetic conservation scores.
Figure 2: Analyses of the AKT2 structure showing alanine (A) 179 and leucine (L) 183 are conserved amino residues.

3.6. Molecular Dynamics Simulation

Through these simulations, we gained valuable insights into the dynamic behavior of the protein, its interactions, and conformational changes that may influence AKT2 structure and function that may lead IR and T2D.

3.6.1. Temperature, Pressure and Density

Solvents and ions were then balanced around the protein structures after energy minimization. NVT (thermostat) and NPT (bariostat) equilibrations were carried out to make sure that the molecules were oriented correctly within the structures. In an isothermal–isobaric ensemble, these two phases were used to stabilize the systems’ temperature and pressure, respectively. The computational tool known as “gms/grompp” was utilized to perform calculations for the NVT and NPT ensembles of the system. The equilibrated process involves temperature, pressure, and density, leading to the establishment of a density plot from NPT. Figure 3 depicts the temperature profiles of the protein systems in their native and mutant forms at a temperature of 300 Kelvin. The plot illustrates significant temperature fluctuations in the conformations during the 100 ps equilibrium phase. Nonetheless, this conduct was not entirely unforeseen, given that the mean temperature of the majority of the conformations fell within a favorable range, with the native structure exhibiting an average temperature of 301 K, and the mutant A179T and L183Q displaying temperatures of 301.5 and 302.9 K, respectively (Figure 3). Likewise, a pressure–time graph (Figure 4) and a density–time graph (illustrated in Figure 5) were generated.
Figure 3: A graph illustrating the variations in temperature over time for both the wildtype and mutant forms of AKT2, as simulated using the GROMACS software. https://www.gromacs.org/. The wildtype is visually represented by the colour purple, while the A179T mutant is depicted in green, and the L183Q mutant is represented by the colour blue.
Figure 4 contains a graph illustrating the variations in pressure over time for both the wildtype and mutant forms of AKT2, as simulated using the GROMACS software. The wildtype is visually represented by the colour purple, while the A179T mutant is depicted in green, and the L183Q mutant is represented by the color blue.
The pressure–time graph (Figure 4) indicated that the wildtype structure exhibited peaks above 110 bars, while the A179T and L183Q mutant structures displayed higher peaks of 260 bars and 180 bars, respectively. The plot for mutant A179T exhibits a local minimum of −500 bar, which deviates from the native structure’s value of −350 bars and that of L183Q is −650 bars. The graph depicting pressure over time indicated that the mean pressure value for all conformations in the native structure was 95 bars. However, for A179T and L183Q, the average pressure values were 105 bars and 100 bar, respectively. The present study examines the density values of the native akt2 protein and its mutant variant A179T over a period of 100 ps (Figure 5).
Figure 5 contains a graph illustrating the variations in density over time for both the wildtype and mutant forms of AKT2, as simulated using the GROMACS software. The wildtype is visually represented by the purple color, while the A179T mutant is depicted in green, and the L183Q mutant is represented by the blue color.
The average density value observed for both the native protein and the mutant variant was 1015 kg/m−3, with minor fluctuations. However, it is noteworthy that the mutant A179T exhibited a lower initial density value of 985 kg/m−3. The L183Q mutant displays analogous fluctuations in pattern to the native AKT2 and mutant A179T; however, it maintains a higher value of 1020.5 Kg/m−3.

3.6.2. Root-Mean-Square Deviation (RMSD) and Root-Mean-Square-Fluctuation (RMSF)

Comparative analysis of the RMSF and RMSD of AKT2 wildtype and its A179T and L183Q mutations was performed. To compare and depict the structural and functional behavior of the wildtype AKT2 gene and its two mutant forms, A179T and L183Q, the current work used a molecular dynamics simulation method. The RMSD was used to evaluate the trajectories of the native and mutant forms of the AKT2 protein. The trajectory file for the C-alpha backbone least-squares fit model and g_rms were used to construct the RMSD graphs shown in Figure 6A.
Figure 6 shows the root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) for the wildtype and the mutants of the AKT2 protein.
The initial conformations of the three structures exhibited close proximity; however, a sudden deviation was observed at 10 ps, resulting in an RMSD value of 0.11 nm for the two mutants.
The wildtype structure achieved its highest RMSD value at 8 ps, while the A179T and L183Q mutants reached their maximum values at 14 ps and 15 ps, respectively. The RMSD graph comparing the wildtype AKT2 gene to its A179T and L183Q mutants, which were aligned to the backbone, indicated that the wildtype AKT2 gene exhibited greater stability than the A179T and L183Q mutants. Subsequently, the aforementioned mutants have been observed to impact the protein dynamics of the wildtype AKT2. This is evidenced by alterations in the RMSD graph of the backbone of the mutants, thereby providing a suitable basis for subsequent analysis. The analysis of RMSF values of native and mutant residues was conducted to ascertain the impact of mutations on the dynamic behavior at the atomic level (Figure 6B). The RMSF of the native residues was found to range from 0.041 nm to 0.38 nm, with an average value of 0.21 nm. The RMSF values for the A179T and L183Q residues were determined to range from 0.032 to 0.38 and 0.04 to 0.41, respectively, with average values of 0.22 and 0.23.

3.6.3. Radius of Gyration (Rg) and Solvent Accessibility Surface Area (SASA)

Analyzing the Radius of Gyration (Rg) and Solvent Accessibility Surface Area (SASA) of both the wildtype and the A179T and L183Q mutations in AKT2 is vital in understanding their impact on AKT2 structure and function. These structural metrics offer insights into the compactness and exposure of amino acids, shedding light on potential conformational changes that contribute to the complex relationship between these mutations and IR and T2D. The concepts of radius of gyration and SASA are of interest in various fields of study. The concept of Rg pertains to the root-mean-square distance that all atoms occupy from their shared center, as determined by their mass and weight. This makes it easier to understand all of the components of protein analysis. Figure 7 compares the normal and mutant protein structures at a temperature of 300 K by plotting the radius of gyration for alpha carbon atoms over time. Upon completion of the simulation, it was observed that the structure of the two mutants, A179T (Figure 7B) and L183Q (Figure 7C), exhibited normal behavior like that of the wildtype (Figure 7A), characterized by consistently low Rg values (<2.63 nm and 2.68, respectively). In contrast, the Rg values for the native structures remaining higher ranged from 3.09 nm to 3.05 nm. SASA refers to the degree of accessibility of the biomolecular surface area to solvent molecules. A reduction in the SASA value signifies a contraction in the structure and a decrease in the surface area. Figure 8 illustrates that the native structure exhibited a significantly elevated SASA value, which ranged from 305 nm2 to 317.5 nm2, with an average of 315 nm2. The SASA values for A179T (Figure 8B) and L183Q (Figure 8B) were observed within the range of 253 nm2 to 260 nm2, and 293.5 nm2 to 298.5 nm2, with an average of 256 nm2 and 296 nm2, respectively.
Figure 7 shows the radius of gyration (Rg) of Cα atoms of native and mutant-type AKT2 protein.
Figure 8: Solvent accessible surface area of wildtype and mutant AKT2 protein.
Figure 9 shows the three-dimensional structure of AKT2 with the sites of mutations indicated with arrows.

4. Discussion

The prevalence of DM is increasing worldwide and it has a serious socioeconomic burden [37]. T2D represents about 90% of all cases of DM [38]. T2D is induced by insulin resistance (IR) or impaired insulin action and pancreatic beta cell demise [39]. IR is a reduced cellular metabolic responsiveness to insulin and a defective lowering effect of circulating or injected insulin on blood sugar [40]. Impaired insulin signaling, for example, because of a defect of AKT2 or a dysfunction of PI3K/AKT signal transduction, can lead to IR [41]. IR is associated with the risk factors of the metabolic syndrome [42]. These risk factor include obesity and dyslipidemia, hyperinsulinemia, high blood pressure, and elevated blood sugar [42,43]. IR and metabolic syndrome are associated with increased susceptibility to T2D and cerebrocardiovascular disease [42,43].
The serine/threonine kinase AKT2, also called protein kinase B (PKB), is implicated in the regulation of insulin signaling pathway. Insulin has important roles in the metabolism of glucose, proteins, and fats [44]. AKT2 is an important mediator in the PI3K pathway, and it binds different downstream proteins that influence cellular metabolism and other process [18]. Gene variations in insulin signaling genes or genes involved in carbohydrate metabolism were reported to be associated with T2D [7,8,45,46]. Insulin enhances the glucose uptake from the circulation in muscles and fat tissues via enhancing the movement of glucose transporter (GLUT4) from their intracellular compartment to the plasma membrane [47]. This movement of GLUT4 from their pool and infusion in the plasma membrane is mediated by AKT2 and its downstream effectors [48]. In the liver, insulin signaling mediated by the AKT2 protein is very important for the inhibition of hepatic glucose synthesis and stimulating lipogenesis [16,49]. Defective AKT2 protein in the pancreatic beta cells causes decreased insulin secretion, which has negative effects on insulin levels in the peripheral tissues [16]. Moreover, it induces insulin insensitively in muscles, adipocytes, and hepatocytes, causing decreased glycogen absorption in muscle cells, reduces glycogenesis in hepatocytes, enhances the breakdown of lipids and induces IR [16]. IR, or the impaired insulin action, that is, irresponsiveness to the insulin in adipose tissues, muscles, liver and other tissues, represents an important component of the pathophysiology of metabolic syndrome [50,51]. Previous studies reported that impaired insulin action is a step preceding the induction of T2D and comorbidities such as abdominal obesity, hyperlipidemia, hepatic steatosis, PCOS, increased blood pressure, and cardiovascular disease [16,50,52]. AKT2 gene variations were associated with PCOS [27], T2D [25], and cancer [53]. The AKT2 variant with alternative splicing is proposed to code the AKT2 protein with a defective C-terminal regulatory region without the hydrophobic motif regulatory site [54]. In the present study, a computational investigation was conducted to assess the effect of nsSNPs on the AKT2 protein structure and function that may lead to the development of IR and T2D.
Our study identified 20 variants as the most detrimental effect based on the collective outcomes of various analytical tools. Then, two potential mutants (A179T, L183Q) were selected from the 20 nsSNPs deemed to be most pathogenic mutants for further investigation through molecular dynamics (MD) simulation. The identified mutations were observed to modify the protein’s stability and function properties, which may have an impact on its biological activity. A molecular dynamics simulation was performed to evaluate the disparities in protein dynamics between the wildtype and mutant AKT2 proteins. The simulation was conducted over a period of 100 ps, during which a range of parameters, including temperature, pressure, density, RMSD, RMSF, SASA, and Rg, were evaluated and scrutinized. Comparisons were performed to investigate the functional characteristics and stability of the protein structures. Upon conducting minimization, marginal variances in temperature, pressure, and density were observed between the native and mutant conformations over the course of the simulation, although the mutant p.L183Q revealed slight increases in these parameters. The RMSD plots of backbone atoms indicated that the wildtype structure demonstrated a higher degree of stability in comparison to the mutants. The RMSF plots indicated that the native and mutant structures exhibited the greatest level of flexibility. The RMSF calculations, which were averaged across the residual elements and trajectories, indicated that the fluctuation of atoms within the structures was notably greater in the first third. The investigation encompassed an analysis of the radius of gyration and SASA to ascertain the impact of the mutant protein structure’s unfolding on the solubility of surface area and the activity of the protein. A higher Rg value is indicative of reduced stability. The present findings suggest that the native structure exhibited a narrower range of Rg values in comparison to the mutants. Moreover, the solvent accessible surface area of the wildtype AKT2 is different from the mutants of the AKT2 protein. This may indicate that the two mutations (p.A179T and p.L183Q) lead to a reduction in AKT2 stability and functionality.
Alanine (A) is the second simplest amino acid; it is hydrophobic. Meanwhile, threonine contains a hydroxyl group; it is a hydrophilic [55]. In p.A179T, the substitution of alanine by threonine can have a profound effect on the AKT2 protein structure and function [56]. Moreover, p.A179T may lead to a loss of helix according to the Mut-pred web server result, which may affect the AKT2 protein. In addition, in p.L183Q, the substitution of leucine (L), a hydrophobic amino acid, by glutamine (Q), a hydrophilic amino acid [55], may exert an impact on the AKT2 protein structure and function. Furthermore, according to the Mut-pred web server, p.L183Q may lead to the loss of acetylation at lysine (K) 181 residue, which may influence the AKT2 protein. This result may be consistent with that of a study which reported that a missense mutation (p.R274H) in the AKT2 generates dysfunctional AKT2 and results in impaired insulin action, hyperinsulinemia, and T2D [29]. Our results are also in agreement with those of a study that reported that a partial loss of function in p.Pro50Thr mutation in AKT2 causes reduced glucose uptake in insulin-sensitive tissues, and it is associated with elevated fasting blood insulin and risk of T2D [31]. To the best of our knowledge, this is the first study highlighting the importance of these mutations (p.A179T and p.L183Q) on the AKT2 protein structure and function. These results need further verification before introduction in the clinic setting via well-designed large-scale case–control, protein–protein interaction (PPI) and protein functional studies [57,58,59,60] to examine the effect of these mutations on the AKT2 protein structure, function and the insulin signaling pathway.

5. Conclusions

In this study, we examined the potential association of non-synonymous SNP (nsSNPs) in the AKT2 gene with IR and T2D. Results indicated that two mutations, p.A179T and p.L183Q, may cause a reduction in AKT2 protein stability and functionality which can be a cause for IR, metabolic syndrome and T2D. These results require further verification in protein functional, PPI and large case–control studies. After being verified, these results can be used in genetic testing for T2D in risk stratification for treatment or prevention strategies.

Author Contributions

M.E.E., I.E. and R.M.; Methodology: M.E.E.; analyses: M.E.E., I.E., R.M., M.A.E., K.M.A., A.O.A., A.K.E., I.A.A., K.S.A., R.A., B.S.O.A., M.H.A., M.M.M. and D.A. Data curation: M.E.E., I.E., R.M., M.A.E., K.M.A., A.O.A., A.K.E., I.A.A., K.S.A., R.A., B.S.O.A., M.H.A., M.M.M. and D.A.; Manuscript drafting M.E.E.; Manuscript Writing Editing: M.E.E. and I.E. 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.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest exist.

References

  1. Roep, B.O.; Thomaidou, S.; van Tienhoven, R.; Zaldumbide, A. Type 1 diabetes mellitus as a disease of the beta-cell (do not blame the immune system?). Nat. Rev. Endocrinol. 2021, 17, 150–161. [Google Scholar] [CrossRef] [PubMed]
  2. Wysham, C.; Shubrook, J. Beta-cell failure in type 2 diabetes: Mechanisms, markers, and clinical implications. Postgrad. Med. 2020, 132, 676–686. [Google Scholar] [CrossRef] [PubMed]
  3. McIntyre, H.D.; Catalano, P.; Zhang, C.; Desoye, G.; Mathiesen, E.R.; Damm, P. Gestational diabetes mellitus. Nat. Rev. Dis. Primers 2019, 5, 47. [Google Scholar] [CrossRef]
  4. Banday, M.Z.; Sameer, A.S.; Nissar, S. Pathophysiology of diabetes: An overview. Avicenna J. Med. 2020, 10, 174–188. [Google Scholar] [CrossRef] [PubMed]
  5. Kolb, H.; Martin, S. Environmental/lifestyle factors in the pathogenesis and prevention of type 2 diabetes. BMC Med. 2017, 15, 131. [Google Scholar] [CrossRef] [PubMed]
  6. Alzahrani, O.R.; Alatwi, H.E.; Alharbi, A.A.; Alessa, A.H.; Al-Amer, O.M.; Alanazi, A.F.R.; Shams, A.M.; Alomari, E.; Naser, A.Y.; Alzahrani, F.A.; et al. Identification and Characterization of Novel Mutations in Chronic Kidney Disease (CKD) and Autosomal Dominant Polycystic Kidney Disease (ADPKD) in Saudi Subjects by Whole-Exome Sequencing. Medicina 2022, 58, 1657. [Google Scholar] [CrossRef]
  7. Elfaki, I.; Mir, R.; Abu-Duhier, F.M.; Khan, R.; Sakran, M. Phosphatidylinositol 3-kinase Glu545Lys and His1047Tyr Mutations are not Associated with T2D. Curr. Diabetes Rev. 2020, 16, 881–888. [Google Scholar] [CrossRef] [PubMed]
  8. Elfaki, I.; Mir, R.; Duhier, F.M.A.; Alotaibi, M.A.; Alalawy, A.I.; Barnawi, J.; Babakr, A.T.; Mir, M.M.; Altayeb, F.; Mirghani, H.; et al. Clinical Implications of MiR128, Angiotensin I Converting Enzyme and Vascular Endothelial Growth Factor Gene Abnormalities and Their Association with T2D. Curr. Issues Mol. Biol. 2021, 43, 1859–1875. [Google Scholar] [CrossRef]
  9. Jha, C.K.; Mir, R.; Elfaki, I.; Javid, J.; Babakr, A.T.; Banu, S.; Chahal, S.M.S. Evaluation of the Association of Omentin 1 rs2274907 A > T and rs2274908 G > A Gene Polymorphisms with Coronary Artery Disease in Indian Population: A Case Control Study. J. Pers. Med. 2019, 9, 30. [Google Scholar] [CrossRef] [PubMed]
  10. Mir, R.; Elfaki, I.; Javid, J.; Barnawi, J.; Altayar, M.A.; Albalawi, S.O.; Jalal, M.M.; Tayeb, F.J.; Yousif, A.; Ullah, M.F.; et al. Genetic Determinants of Cardiovascular Disease: The Endothelial Nitric Oxide Synthase 3 (eNOS3), Kruppel-Like Factor-14 (KLF-14), Methylenetetrahydrofolate Reductase (MTHFR), MiRNAs27a and Their Association with the Predisposition and Susceptibility to Coronary Artery Disease. Life 2022, 12, 1905. [Google Scholar] [PubMed]
  11. Cai, L.; Wheeler, E.; Kerrison, N.D.; Luan, J.; Deloukas, P.; Franks, P.W.; Amiano, P.; Ardanaz, E.; Bonet, C.; Fagherazzi, G.; et al. Genome-wide association analysis of type 2 diabetes in the EPIC-InterAct study. Sci. Data 2020, 7, 393. [Google Scholar] [CrossRef] [PubMed]
  12. Pei, J.; Wang, B.; Wang, D. Current Studies on Molecular Mechanisms of Insulin Resistance. J. Diabetes Res. 2022, 2022, 1863429. [Google Scholar] [CrossRef] [PubMed]
  13. Batista, T.M.; Haider, N.; Kahn, C.R. Correction to: Defining the underlying defect in insulin action in type 2 diabetes. Diabetologia 2022, 65, 1064. [Google Scholar] [CrossRef] [PubMed]
  14. Chu, N.; Viennet, T.; Bae, H.; Salguero, A.; Boeszoermenyi, A.; Arthanari, H.; Cole, P.A. The structural determinants of PH domain-mediated regulation of Akt revealed by segmental labeling. Elife 2020, 9, e59151. [Google Scholar] [CrossRef]
  15. Huang, X.; Liu, G.; Guo, J.; Su, Z. The PI3K/AKT pathway in obesity and type 2 diabetes. Int. J. Biol. Sci. 2018, 14, 1483–1496. [Google Scholar] [CrossRef]
  16. Miao, R.; Fang, X.; Wei, J.; Wu, H.; Wang, X.; Tian, J. Akt: A Potential Drug Target for Metabolic Syndrome. Front. Physiol. 2022, 13, 822333. [Google Scholar] [CrossRef]
  17. Yang, J.; Nie, J.; Ma, X.; Wei, Y.; Peng, Y.; Wei, X. Targeting PI3K in cancer: Mechanisms and advances in clinical trials. Mol. Cancer 2019, 18, 26. [Google Scholar] [CrossRef]
  18. Caruso, M.; Zhang, X.; Ma, D.; Yang, Z.; Qi, Y.; Yi, Z. Novel Endogenous, Insulin-Stimulated Akt2 Protein Interaction Partners in L6 Myoblasts. PLoS ONE 2015, 10, e0140255. [Google Scholar] [CrossRef]
  19. Hoxhaj, G.; Manning, B.D. The PI3K-AKT network at the interface of oncogenic signalling and cancer metabolism. Nat. Rev. Cancer 2020, 20, 74–88. [Google Scholar] [CrossRef]
  20. Alzahrani, O.R.; Mir, R.; Alatwi, H.E.; Hawsawi, Y.M.; Alharbi, A.A.; Alessa, A.H.; Albalawi, E.S.; Elfaki, I.; Alalawi, Y.; Moharam, L.; et al. Potential Impact of PI3K-AKT Signaling Pathway Genes, KLF-14, MDM4, miRNAs 27a, miRNA-196a Genetic Alterations in the Predisposition and Progression of Breast Cancer Patients. Cancers 2023, 15, 1281. [Google Scholar] [CrossRef]
  21. Teng, Y.; Fan, Y.; Ma, J.; Lu, W.; Liu, N.; Chen, Y.; Pan, W.; Tao, X. The PI3K/Akt Pathway: Emerging Roles in Skin Homeostasis and a Group of Non-Malignant Skin Disorders. Cells 2021, 10, 1219. [Google Scholar] [CrossRef]
  22. Zhang, J.; Liu, F. Tissue-specific insulin signaling in the regulation of metabolism and aging. IUBMB Life 2014, 66, 485–495. [Google Scholar] [CrossRef] [PubMed]
  23. Titchenell, P.M.; Quinn, W.J.; Lu, M.; Chu, Q.; Lu, W.; Li, C.; Chen, H.; Monks, B.R.; Chen, J.; Rabinowitz, J.D.; et al. Direct Hepatocyte Insulin Signaling Is Required for Lipogenesis but Is Dispensable for the Suppression of Glucose Production. Cell Metab. 2016, 23, 1154–1166. [Google Scholar] [CrossRef] [PubMed]
  24. Boucher, J.; Kleinridders, A.; Kahn, C.R. Insulin receptor signaling in normal and insulin-resistant states. Cold Spring Harb. Perspect. Biol. 2014, 6, a009191. [Google Scholar] [CrossRef] [PubMed]
  25. Manning, A.; Highland, H.M.; Gasser, J.; Sim, X.; Tukiainen, T.; Fontanillas, P.; Grarup, N.; Rivas, M.A.; Mahajan, A.; Locke, A.E.; et al. A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk. Diabetes 2017, 66, 2019–2032. [Google Scholar] [CrossRef]
  26. Zhu, J.; Wang, M.; He, J.; Zhu, M.; Wang, J.C.; Jin, L.; Wang, X.F.; Yang, Y.J.; Xiang, J.Q.; Wei, Q. Polymorphisms in the AKT1 and AKT2 genes and oesophageal squamous cell carcinoma risk in an Eastern Chinese population. J. Cell. Mol. Med. 2016, 20, 666–677. [Google Scholar] [CrossRef]
  27. Li, J.; Mo, S.; Sun, Y.; Huang, H.; Wang, S. Association between AKT2 gene polymorphism and polycystic ovary syndrome: A case-control study. Gynecol. Endocrinol. 2021, 37, 635–639. [Google Scholar] [CrossRef]
  28. Yin, X.; Xu, Z.; Zhang, Z.; Li, L.; Pan, Q.; Zheng, F.; Li, H. Association of PI3K/AKT/mTOR pathway genetic variants with type 2 diabetes mellitus in Chinese. Diabetes Res. Clin. Pract. 2017, 128, 127–135. [Google Scholar] [CrossRef]
  29. Arya, V.B.; Flanagan, S.E.; Schober, E.; Rami-Merhar, B.; Ellard, S.; Hussain, K. Activating AKT2 mutation: Hypoinsulinemic hypoketotic hypoglycemia. J. Clin. Endocrinol. Metab. 2014, 99, 391–394. [Google Scholar] [CrossRef]
  30. Zhao, X.; An, X.; Yang, C.; Sun, W.; Ji, H.; Lian, F. The crucial role and mechanism of insulin resistance in metabolic disease. Front. Endocrinol. 2023, 14, 1149239. [Google Scholar] [CrossRef]
  31. Latva-Rasku, A.; Honka, M.J.; Stancakova, A.; Koistinen, H.A.; Kuusisto, J.; Guan, L.; Manning, A.K.; Stringham, H.; Gloyn, A.L.; Lindgren, C.M.; et al. A Partial Loss-of-Function Variant in AKT2 Is Associated With Reduced Insulin-Mediated Glucose Uptake in Multiple Insulin-Sensitive Tissues: A Genotype-Based Callback Positron Emission Tomography Study. Diabetes 2018, 67, 334–342. [Google Scholar] [CrossRef] [PubMed]
  32. Zhou, Q.; Sun, W.W.; Chen, J.C.; Zhang, H.L.; Liu, J.; Lin, Y.; Lin, P.C.; Wu, B.X.; An, Y.P.; Huang, L.; et al. Phenylalanine impairs insulin signaling and inhibits glucose uptake through modification of IRbeta. Nat. Commun. 2022, 13, 4291. [Google Scholar] [CrossRef] [PubMed]
  33. Johnson, A.D.; Handsaker, R.E.; Pulit, S.L.; Nizzari, M.M.; O’Donnell, C.J.; de Bakker, P.I. SNAP: A web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 2008, 24, 2938–2939. [Google Scholar] [CrossRef] [PubMed]
  34. Calabrese, R.; Capriotti, E.; Fariselli, P.; Martelli, P.L.; Casadio, R. Functional annotations improve the predictive score of human disease-related mutations in proteins. Hum. Mutat. 2009, 30, 1237–1244. [Google Scholar] [CrossRef]
  35. 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]
  36. Armon, A.; Graur, D.; Ben-Tal, N. ConSurf: An algorithmic tool for the identification of functional regions in proteins by surface mapping of phylogenetic information. J. Mol. Biol. 2001, 307, 447–463. [Google Scholar] [CrossRef]
  37. Al Mansour, M.A. The Prevalence and Risk Factors of Type 2 Diabetes Mellitus (DMT2) in a Semi-Urban Saudi Population. Int. J. Environ. Res. Public Health 2019, 17, 7. [Google Scholar] [CrossRef]
  38. Reed, J.; Bain, S.; Kanamarlapudi, V. A Review of Current Trends with Type 2 Diabetes Epidemiology, Aetiology, Pathogenesis, Treatments and Future Perspectives. Diabetes Metab. Syndr. Obes. 2021, 14, 3567–3602. [Google Scholar] [CrossRef]
  39. Cerf, M.E. Beta cell dysfunction and insulin resistance. Front. Endocrinol. 2013, 4, 37. [Google Scholar] [CrossRef]
  40. Czech, M.P. Insulin action and resistance in obesity and type 2 diabetes. Nat. Med. 2017, 23, 804–814. [Google Scholar] [CrossRef]
  41. Li, M.; Chi, X.; Wang, Y.; Setrerrahmane, S.; Xie, W.; Xu, H. Trends in insulin resistance: Insights into mechanisms and therapeutic strategy. Signal Transduct. Target. Ther. 2022, 7, 216. [Google Scholar] [CrossRef] [PubMed]
  42. Hayden, M.R. Overview and New Insights into the Metabolic Syndrome: Risk Factors and Emerging Variables in the Development of Type 2 Diabetes and Cerebrocardiovascular Disease. Medicina 2023, 59, 561. [Google Scholar] [CrossRef] [PubMed]
  43. Chen, W.; Pang, Y. Metabolic Syndrome and PCOS: Pathogenesis and the Role of Metabolites. Metabolites 2021, 11, 869. [Google Scholar] [CrossRef] [PubMed]
  44. Jaiswal, N.; Gavin, M.G.; Quinn, W.J., 3rd; Luongo, T.S.; Gelfer, R.G.; Baur, J.A.; Titchenell, P.M. The role of skeletal muscle Akt in the regulation of muscle mass and glucose homeostasis. Mol. Metab. 2019, 28, 1–13. [Google Scholar] [CrossRef]
  45. Elfaki, I.; Mir, R.; Mir, M.M.; AbuDuhier, F.M.; Babakr, A.T.; Barnawi, J. Potential Impact of MicroRNA Gene Polymorphisms in the Pathogenesis of Diabetes and Atherosclerotic Cardiovascular Disease. J. Pers. Med. 2019, 9, 51. [Google Scholar] [CrossRef]
  46. Elfaki, I.; Mir, R.; Almutairi, F.M.; Duhier, F.M.A. Cytochrome P450: Polymorphisms and Roles in Cancer, Diabetes and Atherosclerosis. Asian Pac. J. Cancer Prev. 2018, 19, 2057–2070. [Google Scholar]
  47. Satoh, T. Molecular mechanisms for the regulation of insulin-stimulated glucose uptake by small guanosine triphosphatases in skeletal muscle and adipocytes. Int. J. Mol. Sci. 2014, 15, 18677–18692. [Google Scholar] [CrossRef]
  48. Lim, C.Y.; Bi, X.; Wu, D.; Kim, J.B.; Gunning, P.W.; Hong, W.; Han, W. Tropomodulin3 is a novel Akt2 effector regulating insulin-stimulated GLUT4 exocytosis through cortical actin remodeling. Nat. Commun. 2015, 6, 5951. [Google Scholar] [CrossRef]
  49. Eshaghi, F.S.; Ghazizadeh, H.; Kazami-Nooreini, S.; Timar, A.; Esmaeily, H.; Mehramiz, M.; Avan, A.; Ghayour-Mobarhan, M. Association of a genetic variant in AKT1 gene with features of the metabolic syndrome. Genes. Dis. 2019, 6, 290–295. [Google Scholar] [CrossRef]
  50. Saltiel, A.R. Insulin signaling in health and disease. J. Clin. Investig. 2021, 131, 1710–1711. [Google Scholar] [CrossRef]
  51. Fazakerley, D.J.; Krycer, J.R.; Kearney, A.L.; Hocking, S.L.; James, D.E. Muscle and adipose tissue insulin resistance: Malady without mechanism? J. Lipid Res. 2019, 60, 1720–1732. [Google Scholar] [CrossRef] [PubMed]
  52. Piche, M.E.; Tchernof, A.; Despres, J.P. Obesity Phenotypes, Diabetes, and Cardiovascular Diseases. Circ. Res. 2020, 126, 1477–1500. [Google Scholar] [CrossRef] [PubMed]
  53. Sung, J.S.; Park, K.H.; Kim, S.T.; Kim, Y.H. Discovery and Evaluation of Polymorphisms in the AKT2 and AKT3 Promoter Regions for Risk of Korean Lung Cancer. Genom. Inform. 2012, 10, 167–174. [Google Scholar] [CrossRef]
  54. Plotz, G.; Lopez-Garcia, L.A.; Brieger, A.; Zeuzem, S.; Biondi, R.M. Alternative AKT2 splicing produces protein lacking the hydrophobic motif regulatory region. PLoS ONE 2020, 15, e0242819. [Google Scholar] [CrossRef]
  55. Aftabuddin, M.; Kundu, S. Hydrophobic, hydrophilic, and charged amino acid networks within protein. Biophys. J. 2007, 93, 225–231. [Google Scholar] [CrossRef] [PubMed]
  56. Podoly, E.; Hanin, G.; Soreq, H. Alanine-to-threonine substitutions and amyloid diseases: Butyrylcholinesterase as a case study. Chem. Biol. Interact. 2010, 187, 64–71. [Google Scholar] [CrossRef]
  57. Elfaki, I.; Bayer, P.; Mueller, J.W. A potential transcriptional regulator is out-of-frame translated from the metallothionein 2A messenger RNA. Anal. Biochem. 2011, 409, 159–161. [Google Scholar] [CrossRef] [PubMed]
  58. Elfaki, I.; Knitsch, A.; Matena, A.; Bayer, P. Identification and characterization of peptides that bind the PPIase domain of Parvulin17. J. Pept. Sci. 2013, 19, 362–369. [Google Scholar] [CrossRef]
  59. Rodrigues, C.H.M.; Myung, Y.; Pires, D.E.V.; Ascher, D.B. mCSM-PPI2: Predicting the effects of mutations on protein-protein interactions. Nucleic Acids Res. 2019, 47, W338–W344. [Google Scholar] [CrossRef]
  60. Jubb, H.C.; Pandurangan, A.P.; Turner, M.A.; Ochoa-Montano, B.; Blundell, T.L.; Ascher, D.B. Mutations at protein-protein interfaces: Small changes over big surfaces have large impacts on human health. Prog. Biophys. Mol. Biol. 2017, 128, 3–13. [Google Scholar] [CrossRef]
Figure 1. Workflow of the nsSNP detection and analysis. The nucleotide sequences are retrieved in FASTA format, and SNPs are assessed whether synonymous, nonsynonymous, or nonsense. Nonsynonymous SNPs (nsSNPs) are assessed on the subject of whether they have deleterious effects on the protein encoded by the gene by molecular dynamic simulation.
Figure 1. Workflow of the nsSNP detection and analysis. The nucleotide sequences are retrieved in FASTA format, and SNPs are assessed whether synonymous, nonsynonymous, or nonsense. Nonsynonymous SNPs (nsSNPs) are assessed on the subject of whether they have deleterious effects on the protein encoded by the gene by molecular dynamic simulation.
Cimb 45 00471 g001
Figure 2. The ConSurf analysis yielded results pertaining to the conservation of residues. The ConSurf results display a range of colours that correspond to the level of confidence regarding sequence conservation. In this colour scheme, the sky-blue colour represents variable residues, while the dark purple colour represents highly conserved residues.
Figure 2. The ConSurf analysis yielded results pertaining to the conservation of residues. The ConSurf results display a range of colours that correspond to the level of confidence regarding sequence conservation. In this colour scheme, the sky-blue colour represents variable residues, while the dark purple colour represents highly conserved residues.
Cimb 45 00471 g002
Figure 3. The variations in temperature over time for both the wildtype and mutant forms of AKT2.
Figure 3. The variations in temperature over time for both the wildtype and mutant forms of AKT2.
Cimb 45 00471 g003
Figure 4. The variations in pressure over time for both the wildtype and mutant forms of AKT2.
Figure 4. The variations in pressure over time for both the wildtype and mutant forms of AKT2.
Cimb 45 00471 g004
Figure 5. The variations in density over time for the wildtype and mutant AKT2 protein.
Figure 5. The variations in density over time for the wildtype and mutant AKT2 protein.
Cimb 45 00471 g005
Figure 6. The root-mean-square deviation (RMSD) (A) and root-mean-square fluctuation (RMSF) (B) graphs were generated for the backbone atoms of both the wildtype and mutant variants of the AKT2 protein. These simulations were conducted using GROMACS version 5.1.2. The wildtype is visually represented by the colour purple, while the A179T mutant is depicted in green and the L183Q mutant is depicted in blue.
Figure 6. The root-mean-square deviation (RMSD) (A) and root-mean-square fluctuation (RMSF) (B) graphs were generated for the backbone atoms of both the wildtype and mutant variants of the AKT2 protein. These simulations were conducted using GROMACS version 5.1.2. The wildtype is visually represented by the colour purple, while the A179T mutant is depicted in green and the L183Q mutant is depicted in blue.
Cimb 45 00471 g006
Figure 7. The radius of gyration (Rg) of Cα atoms of native and mutant-type AKT2 protein. Wildtype (A) is shown in black, A179T (B) mutant is shown in blue and L183Q (C) is shown in green.
Figure 7. The radius of gyration (Rg) of Cα atoms of native and mutant-type AKT2 protein. Wildtype (A) is shown in black, A179T (B) mutant is shown in blue and L183Q (C) is shown in green.
Cimb 45 00471 g007
Figure 8. Solvent accessible surface area of wildtype (A) and mutant-type AKT protein. Wildtype is shown in black, A179T (B) mutant is shown in Blue and L183Q (C) is shown in green.
Figure 8. Solvent accessible surface area of wildtype (A) and mutant-type AKT protein. Wildtype is shown in black, A179T (B) mutant is shown in Blue and L183Q (C) is shown in green.
Cimb 45 00471 g008
Figure 9. The cartoon structure of the AKT2 protein (PDB ID: 3D0E) showing the sites of the mutations p.A179T and p.L183Q with surface presentations indicated with arrows. This figure is prepared using YASARA view.
Figure 9. The cartoon structure of the AKT2 protein (PDB ID: 3D0E) showing the sites of the mutations p.A179T and p.L183Q with surface presentations indicated with arrows. This figure is prepared using YASARA view.
Cimb 45 00471 g009
Table 1. The effects of amino acid substitution on the AKT2 gene as predicted by SIFT (deleterious or tolerated), PolyPhen (damaging or benign), SNAP2 and PANTHER.
Table 1. The effects of amino acid substitution on the AKT2 gene as predicted by SIFT (deleterious or tolerated), PolyPhen (damaging or benign), SNAP2 and PANTHER.
SIFTPolyphenSNAP2PANTHER
Variant IDAlleleProteinMutationScorePredictionScorePredictionEffectsExpected Accuracy %pdelPrediction
rs1017033490C|TENST00000311278.10R480H0.00deleterious0.988probably damagingEffect65%0.57probably damaging
rs769806554G|AENST00000311278.10R480C0.00deleterious0.992probably damagingEffect63%0.57probably damaging
rs769806554G|CENST00000311278.10R480G0.00deleterious0.986probably damagingEffect75%0.57probably damaging
rs1447335111T|CENST00000311278.10Y475C0.00deleterious1probably damagingEffect80%0.85probably damaging
rs1447335111T|GENST00000311278.10Y475S0.00deleterious0.991probably damagingEffect91%0.85probably damaging
rs1368133745A|TENST00000311278.10F470I0.00deleterious1probably damagingEffect80%0.85probably damaging
rs1251045508G|CENST00000311278.10S432C0.01deleterious0.999probably damagingEffect61%0.85probably damaging
rs757239082A|GENST00000311278.10M404T0.09deleterious0.644possible damagingEffect82%0.85probably damaging
rs757239082A|CENST00000311278.10V403G0.00deleterious1probably damagingEffect80%0.57probably damaging
rs1278736092C|GENST00000311278.10G396R0.01deleterious0.999probably damagingEffect53%0.86probably damaging
rs1207581157C|TENST00000311278.10G396R0.01deleterious0.999probably damagingEffect53%0.86probably damaging
rs1316408057C|AENST00000311278.10G395W0.01deleterious1probably damagingEffect66%0.74probably damaging
rs1371711287C|AENST00000311278.10G394C0.00deleterious1probably damagingEffect63%0.89probably damaging
rs867381867A|TENST00000311278.10L385H0.00deleterious1probably damagingEffect91%0.89probably damaging
rs767039609G|AENST00000311278.10A377V0.00deleterious1probably damagingEffect63%0.86probably damaging
rs752215436G|AENST00000311278.10P370L0.00deleterious1probably damagingEffect53%0.89probably damaging
rs754942500G|AENST00000311278.10R357C0.00deleterious1possible damagingEffect72%0.57probably damaging
rs749701848G|AENST00000311278.10P349L0.00deleterious0.998probably damagingEffect53%0.89probably damaging
rs1470884234G|AENST00000311278.10R347C0.01deleterious0.985probably damagingEffect85%0.85probably damaging
rs1198321649C|TENST00000311278.10G346S0.01deleterious1probably damagingEffect66%0.89probably damaging
rs1485158704A|TENST00000311278.10M344K0.01deleterious1probably damagingEffect71%0.85probably damaging
rs765026585G|AENST00000311278.10R329W0.00deleterious1possible damagingEffect78%0.86probably damaging
rs765026585C|AENST00000311278.10D326Y0.03deleterious1probably damagingEffect80%0.85probably damaging
rs1257580103G|AENST00000311278.10A318V0.00deleterious1probably damagingEffect71%0.89probably damaging
rs1458251691T|CENST00000311278.10E315G0.00deleterious0.995probably damagingEffect75%0.86probably damaging
rs1437512007G|AENST00000311278.10P314L0.00deleterious0.998probably damagingEffect80%0.86probably damaging
rs753989242G|CENST00000311278.10P314A0.00deleterious1probably damagingEffect80%0.86probably damaging
rs529742537T|CENST00000311278.10M307V0.02deleterious0.998probably damagingEffect71%0.74probably damaging
rs200395003C|TENST00000311278.10D284N0.04deleterious0.994possible damagingEffect53%0.86probably damaging
rs777146223C|TENST00000311278.10R274H0.00deleterious0.994probably damagingEffect85%0.89probably damaging
rs139506765G|AENST00000311278.10R274C0.00deleterious0.999probably damagingEffect95%0.89probably damaging
rs778561687C|TENST00000311278.10V271M0.00deleterious1probably damagingEffect95%0.74probably damaging
rs778561687G|AENST00000311278.10R269W0.02deleterious0.948probably damagingEffect80%0.57probably damaging
rs745550650A|TENST00000311278.10Y265N0.01deleterious0.985probably damagingEffect85%0.86probably damaging
rs1568518681G|CENST00000311278.10A252G0.01deleterious1probably damagingEffect85%0.78probably damaging
rs930976777G|AENST00000311278.10R251W0.00deleterious0.9probably damagingEffect66%0.86probably damaging
rs1568518786G|AENST00000311278.10R245C0.01deleterious1probably damagingEffect75%0.89probably damaging
rs753149827G|AENST00000311278.10R243W0.00deleterious0.976probably damagingEffect82%0.85probably damaging
rs756750912T|CENST00000311278.10N233D0.01deleterious0.992probably damagingEffect66%0.86probably damaging
rs1020231672A|GENST00000311278.10Y231H0.00deleterious0.999probably damagingEffect63%0.86probably damaging
rs145305228T|GENST00000311278.10T221P0.00deleterious1probably damagingEffect85%0.89probably damaging
rs1169786300G|CENST00000311278.10L215V0.00deleterious0.997probably damagingEffect80%0.86probably damaging
rs372625056A|GENST00000311278.10L212P0.00deleterious0.997probably damagingEffect71%0.86probably damaging
rs1166756826G|CENST00000311278.10P210R0.01deleterious1probably damagingEffect95%0.89probably damaging
rs1422114933G|AENST00000311278.10R202W0.00deleterious1probably damagingEffect71%0.85probably damaging
rs1004075798T|CENST00000311278.10T197A0.01deleterious0.875probably damagingEffect85%0.86probably damaging
rs776561628G|AENST00000311278.10R184W0.00deleterious0.76possible damagingEffect71%0.57probably damaging
rs1568519097A|TENST00000311278.10L183Q0.00deleterious0.996probably damagingEffect71%0.86probably damaging
rs1323113005C|TENST00000311278.10A179T0.01deleterious0.995probably damagingEffect66%0.89probably damaging
rs763266151T|GENST00000311278.10T174P0.01deleterious0.99probably damagingEffect91%0.85probably damaging
rs144395843G|AENST00000311278.10R170W0.00deleterious0.994probably damagingEffect85%0.57probably damaging
rs980537301C|TENST00000311278.10V169M0.01deleterious0.997probably damagingEffect80%0.85probably damaging
rs140987550A|GENST00000311278.10M149T0.02deleterious0.907probably damagingEffect71%0.74probably damaging
rs762379584G|AENST00000311278.10S126F0.00deleterious0.99probably damagingEffect66%0.74probably damaging
rs121434593A|GENST00000311278.10I103T0.02deleterious0.823probably damagingEffect59%0.85probably damaging
rs774959037T|AENST00000311278.10I103F0.02deleterious0.999probably damagingEffect66%0.85probably damaging
rs1304063423C|AENST00000311278.10R101W0.02deleterious0.913probably damagingEffect59%0.5possibly damaging
rs1347953319C|TENST00000311278.10E97G0.02deleterious0.99probably damagingEffect71%0.85probably damaging
rs757245351C/AENST00000311278.10D91Y0.00deleterious0.732possible damagingEffect66%0.74probably damaging
rs758367800A/GENST00000311278.10F88L0.01deleterious0.999probably damagingEffect66%0.86probably damaging
rs747207035C/TENST00000311278.10E85K0.02deleterious0.945probably damagingEffect66%0.86probably damaging
rs1207387979C/A/TENST00000311278.10V74F0.04deleterious1probably damagingEffect59%0.5possibly damaging
rs145907048A/CENST00000311278.10F73V0.00deleterious0.93probably damagingEffect75%0.85probably damaging
rs1428621686G/AENST00000311278.10P68L0.01deleterious1probably damagingEffect53%0.85probably damaging
rs773553776G/A/CENST00000311278.10P51R0.03deleterious0.825probably damagingEffect75%0.85probably damaging
rs771461640G/CENST00000311278.10P51A0.03deleterious0.999probably damagingEffect63%0.85probably damaging
rs1367814764T/GENST00000311278.10Y26S0.00deleterious0.791possible damagingEffect91%0.85probably damaging
rs775105991C/TENST00000311278.10R25Q0.00deleterious1probably damagingEffect91%0.86probably damaging
rs1363221582A/GENST00000311278.10I19T0.00deleterious0.987probably damagingEffect80%0.86probably damaging
rs763896696C/A/TENST00000311278.10E17K0.05deleterious1probably damagingEffect85%0.86probably damaging
rs988297708G/AENST00000311278.10R15C0.00deleterious1probably damagingEffect80%0.86probably damaging
rs1227145174T/GENST00000311278.10K14T0.00deleterious1probably damagingEffect91%0.86probably damaging
rs754286414C/TENST00000311278.10E9K0.00deleterious0.998probably damagingEffect91%0.85probably damaging
rs1373751808T/A/CENST00000311278.10I7F0.02deleterious0.93probably damagingEffect71%0.27probably benign
rs746407532C/AENST00000311278.10V6F0.01deleterious0.752possible damagingEffect66%0.57probably damaging
Table 2. Disease association of the non-synonymous SNP (nsSNPs) in AKT2 protein as predicted by SNP&GO and PhD-SN.
Table 2. Disease association of the non-synonymous SNP (nsSNPs) in AKT2 protein as predicted by SNP&GO and PhD-SN.
PhD-SNPSNP&GO
Variant IDAlleleMutationScorePredictionScorePrediction
rs1447335111T|CY475C3Disease1Disease
rs1447335111T|GY475S6Disease2Disease
rs1368133745A|TF470I8Disease5Disease
rs1371711287C|AG394C4Disease4Disease
rs749701848G/AP349L4Disease6Disease
rs1198321649C/TG346S0Disease7Disease
rs1485158704A/TM344K9Disease8Disease
rs765026585C/AD326Y4Disease5Disease
rs1257580103G/AA318V7Disease6Disease
rs930976777G/AR251W2Disease5Disease
rs753149827G/AR243W2Disease2Disease
rs145305228T/GT221P5Disease3Disease
rs1166756826G/CP210R4Disease6Disease
rs1422114933G/AR202W1Disease3Disease
rs1568519097A/TL183Q0Disease4Disease
rs1323113005C/TA179T2Disease7Disease
rs763266151T/GT174P7Disease4Disease
rs758367800A/GF88L8Disease1Disease
rs145907048A/CF73V2Disease4Disease
rs775105991C/TR25Q7Disease3Disease
Table 3. Predication of the effects of non-synonymous SNP (nsSNPs) in the AKT2 protein stability (increasing or decreasing).
Table 3. Predication of the effects of non-synonymous SNP (nsSNPs) in the AKT2 protein stability (increasing or decreasing).
Variant IDAlleleMutationI-MutantMUpro
rs1447335111T|CY475CDecreaseDecrease
rs1447335111T|GY475SDecreaseDecrease
rs1368133745A|TF470IDecreaseDecrease
rs1371711287C|AG394CDecreaseDecrease
rs749701848G/AP349LDecreaseDecrease
rs1198321649C|TG346SDecreaseDecrease
rs1485158704A|TM344KDecreaseDecrease
rs765026585C|AD326YDecreaseDecrease
rs1257580103G|AA318VDecreaseIncrease
rs930976777G|AR251WDecreaseDecrease
rs753149827G|AR243WDecreaseDecrease
rs145305228T|GT221PDecreaseDecrease
rs1166756826G|CP210RDecreaseIncrease
rs1422114933G|AR202WDecreaseIncrease
rs1568519097A|TL183QDecreaseDecrease
rs1323113005C|TA179TDecreaseDecrease
rs763266151T|GT174PDecreaseDecrease
rs758367800A|GF88LDecreaseDecrease
rs145907048A|CF73VDecreaseDecrease
rs775105991C|TR25QDecreaseDecrease
Table 4. Predict the pathogenicity and underlying functional alterations of non-synonymous SNP (nsSNPs) in AKT2 protein.
Table 4. Predict the pathogenicity and underlying functional alterations of non-synonymous SNP (nsSNPs) in AKT2 protein.
Mut-Pred
Variant IDAlleleMutationScoreEffectsFunction Affected
rs1447335111T|CY475C0.9(−)Loss of Relative solvent accessibility
rs1371711287C|AG394C0.9(−)Loss of Intrinsic disorder; Loss of B-factor; Loss of Methylation at K390
rs765026585C|AD326Y0.9(−)Loss of Relative solvent accessibility
rs930976777G|AR251W0.8(−)Loss of Helix; Loss of Allosteric site at Y255
rs1422114933G|AR202W0.6(−)Loss of Helix
rs1568519097A|TL183Q0.9(−)Loss of Acetylation at K181
rs1323113005C|TA179T0.9(−)Loss of Helix
rs763266151T|GT174P0.8(−)Loss of Strand
rs758367800A|GF88L0.8(−)Loss of Strand
rs145907048A|CF73V0.9(−)Loss of Strand
(−): negative impact on protein function.
Table 5. Detailed summary of important non-synonymous SNP (nsSNPs) identified.
Table 5. Detailed summary of important non-synonymous SNP (nsSNPs) identified.
Functional ConsequenceDisease AssociationPhylogeneticsProtein StabilityFunction
nsSNPMutationSIFTPolyphenSNAP2PANTHERSNP&GOPhD-SNPConSurfI-MutMUproMut-Pred
rs1447335111Y475C********5, b, average******
rs1447335111Y475S********4, b, average******+
rs1368133745F470I********9, b,s, conserved******+
rs1371711287G394C********9, e,f, conserved******
rs749701848P349L********9, e,f, conserved******+
rs1198321649G346S********9, e,f, conserved******+
rs1485158704M344K********7, b, conserved******+
rs765026585D326Y********4, e, average******
rs1257580103A318V********9, b,s, conserved***Increase+
rs930976777R251W********5, e, average******
rs753149827R243W********5, e, average******+
rs145305228T221P********9, e,f, conserved******+
rs1166756826P210R********9, e,f, conserved***Increase+
rs1422114933R202W********5, e, average***Increase
rs1568519097L183Q********6, b, conserved******
rs1323113005A179T********9, b,s, conserved******
rs763266151T174P********5, e, average******
rs758367800F88L********7, b, conserved******
rs145907048F73V********6, b, conserved******
rs775105991R25Q********8, b,f, conserved******+
*: deleterious; **: disease-associated; ***: decrease protein stability. e—An exposed residue according to the neural network algorithm. f—A predicted functional residue (highly conserved and exposed). b—A buried residue according to the neural network algorithm. s—A predicted structural residue (highly conserved and buried); (−) negative effects; (+) positive effects.
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

Elangeeb, M.E.; Elfaki, I.; Elkhalifa, M.A.; Adam, K.M.; Alameen, A.O.; Elfadl, A.K.; Albalawi, I.A.; Almasoudi, K.S.; Almotairi, R.; Alsaedi, B.S.O.; et al. In Silico Investigation of AKT2 Gene and Protein Abnormalities Reveals Potential Association with Insulin Resistance and Type 2 Diabetes. Curr. Issues Mol. Biol. 2023, 45, 7449-7475. https://doi.org/10.3390/cimb45090471

AMA Style

Elangeeb ME, Elfaki I, Elkhalifa MA, Adam KM, Alameen AO, Elfadl AK, Albalawi IA, Almasoudi KS, Almotairi R, Alsaedi BSO, et al. In Silico Investigation of AKT2 Gene and Protein Abnormalities Reveals Potential Association with Insulin Resistance and Type 2 Diabetes. Current Issues in Molecular Biology. 2023; 45(9):7449-7475. https://doi.org/10.3390/cimb45090471

Chicago/Turabian Style

Elangeeb, M. E., Imadeldin Elfaki, M. A. Elkhalifa, Khalid M. Adam, A. O. Alameen, Ahmed Kamaleldin Elfadl, Ibrahim Altedlawi Albalawi, Kholoud S. Almasoudi, Reema Almotairi, Basim S. O. Alsaedi, and et al. 2023. "In Silico Investigation of AKT2 Gene and Protein Abnormalities Reveals Potential Association with Insulin Resistance and Type 2 Diabetes" Current Issues in Molecular Biology 45, no. 9: 7449-7475. https://doi.org/10.3390/cimb45090471

APA Style

Elangeeb, M. E., Elfaki, I., Elkhalifa, M. A., Adam, K. M., Alameen, A. O., Elfadl, A. K., Albalawi, I. A., Almasoudi, K. S., Almotairi, R., Alsaedi, B. S. O., Alhelali, M. H., Mir, M. M., Amle, D., & Mir, R. (2023). In Silico Investigation of AKT2 Gene and Protein Abnormalities Reveals Potential Association with Insulin Resistance and Type 2 Diabetes. Current Issues in Molecular Biology, 45(9), 7449-7475. https://doi.org/10.3390/cimb45090471

Article Metrics

Back to TopTop