A Simulation Analysis and Screening of Deleterious Nonsynonymous Single Nucleotide Polymorphisms (nsSNPs) in Sheep <i>LEP</i> Gene

Girmay, Shishay; Ahmad, Hafiz Ishfaq; Zahra, Quratul Ain

doi:https://doi.org/10.1155/2022/7736485

BioMed Research International

On this page

Abstract Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 7736485 | https://doi.org/10.1155/2022/7736485

A Simulation Analysis and Screening of Deleterious Nonsynonymous Single Nucleotide Polymorphisms (nsSNPs) in Sheep LEP Gene

Shishay Girmay,¹Hafiz Ishfaq Ahmad,²and Quratul Ain Zahra³

Academic Editor: Amjad Bashir

Received24 May 2022

Accepted22 Jul 2022

Published08 Aug 2022

Abstract

Leptin is a polypeptide hormone produced in the adipose tissue and governs many processes in the body. Recently, polymorphisms in the LEP gene revealed a significant change in body weight regulation, energy balance, food intake, and reproductive hormone secretion. This study considers its crucial role in the regulation of the economically important traits of sheep. Several computational tools, including SIFT, Predict SNP2, SNAP2, and PROVEAN, have been used to screen out the deleterious nsSNPs. Following the screening of 11 nsSNPs in the sheep genome, 5 nsSNPs, T86M (C → T), D98N (G → A), N136T (A → C), R142Q (G → A), and P157Q (C → A), were predicted to have a significant deleterious effect on the LEP protein function, leading to phenotypic difference. The analysis of proteins’ stability change due to amino acid substitution using the I-stable, SDM, and DynaMut consistently confirmed that three nsSNPs (T86M (C → T), D98N (G → A), and P157Q (C → A)) increased protein stability. It is suggested that these three nsSNPs may enhance the evolvability of LEP protein, which is vital for the evolutionary adaptation of sheep. Our findings demonstrate that the five nsSNPs reported in this study might be responsible for sheep’s structural and functional modifications of LEP protein. This is the first comprehensive report on the sheep LEP gene. It narrow downs the candidate nsSNPs for in vitro experiments to facilitate the development of reliable molecular markers for associated traits.

1. Introductions

Leptin is one of the major hormones secreted by adipocytes. It is the primary function of regulating homeostatic control of energy balance, metabolism, neuroendocrine system, and other functions through its effects on the central nervous system [1, 2].

The short negative feedback maintains the body leptin proportion through pituitary cells, or via the long feedback performed by neurosecretory cells, interacting with corresponding receptors in the hypothalamus (paraventricular, lateral, ventromedial, and dorsomedial nuclei), leads to repression of orexigenic peptide production and stimulation of anorexigenic factors [3]. Then, the brain regulates the balance of energy expenditure from the body and the amount of energy stored in the body of the organism [4]. Usually, blood leptin levels are proportionate to the mass of adipose tissue, so the fatter individuals have, the more leptin they have circulating in their blood. Nonetheless, people and animals with a higher level of adipose tissue seldom show leptin synthesis or leptin production [5]. The reduction in the hormone’s effectiveness may be associated with a fault in drawing hormonal signals or the ability of leptin to penetrate the blood-brain barrier, resulting in leptin resistance [6]. In turn, leptin deficiency or polymorphisms in the leptin pathways increase appetite and energy intake, ultimately leading to obesity [7].

Single nucleotide polymorphisms are essential and the most common types of mutations associated with complex traits [8]. In particular, the missense mutations located in the coding region, which alter the amino acid configuration, may cause a significant change in the structure and function of the native protein [9]. In vitro, functional comparison of mutated proteins with their corresponding wild-type proteins associated with specific traits of animals is a common practice that enables us to identify the impact of each nsSNP in the corresponding protein [10, 11]. However, the experimental design for each mutational change is laborious and time-consuming. Therefore, it is reasonable and economical to carry out the basic work required for mutation data mining and predict its effect on protein properties through computational analysis (Nailwal and Chauhan, 2017). Screening out the deleterious nsSNPs with a significant consequence on phenotype from the tolerant nsSNPs (without phenotypic changes) has a pivotal role in understanding the molecular basis of the polygenetic traits in sheep. Therefore, considering the role played by the LEP gene in regulating body weight, energy balance, and feed intake traits in sheep, we retrieved all of the mutations found in sheep LEP gene deposited in the Ensembl sheep genome browser until March 2022. The objective of this study was to characterize the potential variations in the ovine LEP gene using in silico analysis.

2. Materials and Methods

We used various computational tools to screen out the functional effects of the nsSNPs in the LEP gene. Details of each bioinformatics tool are presented below and summarized in Figure 1.

2.1. Collecting SNPs and Protein’s Sequence from the Databases

Sheep leptin (LEP) gene SNPs were collected from the Ovis aries Ensembl genome browser 106 (https://uswest.ensembl.org/Ovis_aries/Info/Index) [12]. The transcript sequence and the protein encoded by the sheep LEP gene were retrieved from the Ensemble database (https://ensembl.org/) [13]. Then, the UniProt ID for the amino acid (UniProtKB - W5NWK7) was obtained from UniProt Protein Database (http://www.uniprot.org).

2.2. Nonsynonymous SNP Functional Analysis for LEP

Four tools, SIFT, PredictSNP2, PROVEAN, and SNAP2, were used to predict the practical context of missense mutations in the LEP gene.

Firstly, the amino acid substitutions that affect protein structure and function lead to phenotypic changes were screened by SIFT (Sorting Intolerant From Tolerant; http://sift.jcvi.org/) tool [14]. The identification numbers (rsIDs) of each nsSNP of the LEP gene obtained from Ensembl were submitted as a query to SIFT for homology searching. Results were obtained as SIFT scores which were classified as damaging (0.00–0.05), potentially damaging (0.051–0.10), borderline (0.101–0.20), or tolerant (0.201–1.00) [15]. Then, PredictSNP2 (http://loschmidt.chemi.muni.cz/predictsnp2) was used to identify the deleterious nsSNPs. PredictSNP2 is a consensus classifier for predicting the effect of amino acid substitutions based on the output of eight different amino acid-based predictors http://loschmidt.chemi.muni.cz/predictsnp [16].

PROVEAN (Protein Variation Effect Analyzer: http://provean.jcvi.org/index.php) is used to predict the level of impact on protein structure and biological function. A protein FASTA sequence with amino acid alterations is used as the input query. It classifies nsSNPs as deleterious or neutral based on whether the final score falls below a threshold of − 2.5. Scores over this threshold are considered neutral [17].

SNAP2 predicts the functional impact of mutations [18]. SNAP2 predicts the functional impact of mutations (Turkson, 2004). SNAP2 is a taught classifier based on the “neural network” machine learning device. It uses a number of sequence and variant properties to discriminate between effect and neutral variants/nonsynonymous SNPs. The SNAP2 achieved sustained two-state accuracy (effect/neutral) of 82 percent in cross-validation across 100,000 experimentally annotated variations (at an AUC of 0.9). In other words, this represents a major improvement over existing methods [19] (the website https://rostlab.org/services/snap2web/ is where you may find it).

2.3. Analysis of Protein’s Stability Change upon Amino Acid Substitution

Protein evolution is mainly governed by protein stability [20]. After folding, we analyzed the relationship between LEP mutations and protein stability based on a smaller free energy change (ΔG or dG). In contrast, the difference in folding free energy change between wild-type and mutant protein (ΔΔG or ddG) is often considered an impact factor on protein stability changes [21]. The protein structure stability and the molecular and structural effects of protein-coding variants were predicted using the integrated predictor for protein stability change upon a single mutation called iStable (http://predictor.nchu.edu.tw/iStable) [22]. Furthermore, the change in protein stability upon mutation was estimated using a site-directed mutator (SDM) (http://marid.bioc.cam.ac.uk/sdm2/prediction) a server [23]. Finally, the DynaMut (http://biosig.unimelb.edu.au/dynamut/) was used to verify the impact of mutations on protein conformation, flexibility, and stability and to visualize the protein dynamics [24].

2.4. Structural Conformation and Conservation Analysis of LEP

Highly conserved functional regions of the protein coded by the LEP gene was identified by ConSurf tool (http://consurf.tau.ac.il/) [25–27].

2.5. Prediction of Secondary Structure

The secondary structure of LEP was predicted using the PSIPRED server available at (http://bioinf.cs.ucl.ac.uk/psipred/) [28]. It is based on a two-stage neural network with the implementation of position-specific scoring matrices constructed from PSI-BLAST to predict the available secondary structures of a protein [29].

2.6. Homology Modelling

The 3D structures of the protein encoded by LEP were constructed using four different homology modeling tools. No crystal structure with an appropriate length of this protein was available in the protein data bank.

2.6.1. Homology Modelling by Swiss-Model Server

Swiss-Model workplace available at (https://swissmodel.expasy.org/), a web-based tool for the homology modeling of protein, was used to avail the three-dimensional structure of LEP [30, 31]. The CAMEO system calculates the accuracy of the constructed model. Swiss-Model is an automated tool based on evolutionary information which searches for the best sequence-template alignment from its high-throughput template library (SMTL) to build the model [32].

2.6.2. Homology Modelling by Phyre-2 Server

Phyre-2 tool (http://www.sbg.bio.ic.ac.uk/phyre2/), based on the Hidden Markov method, was used to predict the homology-based three-dimensional structure of the query amino acid sequence. It combines five steps to build a model: (1) collecting homologous sequence, (2) screening of fold library, (3) loop modelling, (4) multiple template modeling by ab initio folding simulation Poing, and (5) side chain placement [33, 34].

2.6.3. Homology Modeling by RaptorX Server

RaptorX server available at (http://raptorx.uchicago.edu/StructurePrediction/predict/), using RaptorX-Boost and RaptorX-MSA to construct the three-dimensional structure of a protein, predicted the 3D model of the protein coded by CDKN1A. It combines a nonlinear scoring function and a probabilistic consistency algorithm to predict the model structure [35].

2.6.4. Verification of Three-Dimensional Model of Protein

RAMPAGE Ramachandran plot analysis was used for the verification of 3D structures. It provides the number of residues in the favored, allowed, and outlier regions [36]. If a good proportion of residues lie in the favored and allowed region, then the model is predicted to be good.

3. Results

3.1. Data Mining

Data mining of the sheep genome browser from Ensembl database on March 12, 2022, revealed 543 SNPs of the LEP gene in sheep, of which 11 were missense mutations. Detailed information about these SNPs is shown in Table 1.

3.2. Nonsynonymous SNPs Functional Analysis for LEP

Out of the 11 nsSNPs subjected to SIFT tool to predict their impact level on the LEP protein function as deleterious or tolerated, the nine nsSNPs have shown a damaging effect with an average prediction score of 0.00–0.02, and the remaining two nsSNPs were tolerant (T86M and R196Q) Table 2.

PredictSNP2 further evaluated all the 11nsSNPs, including the nine predicted as deleterious by SIFT. The PredictSNP2 is a consensus classifier tool for predicting the effect of amino acid substitutions based on the output of eight computational prediction tools. The PredictSNP2 tool prediction outcomes revealed that eight nsSNPs were predicted as deleterious (Table S1). The SNAP2 was used to confirm further the sequence variants’ anticipated effects on the LEP protein function.

3.3. Mutant Protein Stability Prediction for LEP

The iStable Meta server tool demonstrated that the amino acid changes in R2C, L35P, D98N, N136T, R142Q, P157Q, V181L, and R196Q cause loss and decreased stability of the LEP protein. However, the other three mutations (T75M, T86M, and V94I) increase the protein stability Table 3. These results were cross-checked with a site-directed mutator (SDM), and DynaMut servers found that eight nsSNPs significantly change the stability compared to the wild-type proteins.

Double-checking with site-directed mutator (SDM) and DynaMut servers consistently verified the molecular consequences of 5 nsSNPs. Out of these, three nsSNPs (T86M, D98N, and P157Q) increased stability by decreasing the molecular flexibility of the wild-type proteins. However, the other two nsSNPs (N136T and R142Q) revealed a decrease in stability and molecular flexibility Table S4.

Furthermore, the interatomic interactions prediction outcomes between wild type (left side) and after single point mutations (right side) are presented in Figures 2(a)–2(e). The wild-type and mutant residues are colored in light green and represented as sticks alongside the surrounding residues that take part in any interactions.

(a)

(b)

(c)

(d)

(e)

3.4. Structural Conformation and Conservation Analysis by ConSurf Server

Our ConSurf results indicated that nsSNPs at position L35P, T75M, V94I, D98N, and R142Q are in the highly conserved region, with a conservation score 9. Likewise, nsSNPs R2C, N136T, and P157Q have shown a conservation score of 8 (Figure S1). From this, we can speculate that these nsSNPs have a potential effect on LEP protein.

3.5. LEP Protein Secondary Structure Prediction by PSIPRED

The alpha helix and beta sheet distribution and the coil are exposed by PSIPRED (Figure S2). Among the secondary structures, the highest in percentage was coils (51.5%) followed by alpha helix (48.5%) and no beta-sheet (0.0%).

3.6. Homology Modelling

Structure refined and energy minimized homology models of LEP by Swiss-Model, Phyre-2, ConSurf, and RaptorX servers are illustrated using computational program QMEAN (Qualitative Model Energy ANalysis) (Figure S3). These models were checked for validation through Ramachandran plot analysis. Some residues in the model may lie in favorable (blue), allowed (green), or disallowed (red) regions of the Ramachandran plot. This coloring indicates residues that may have problems with the backbone phi/psi angles.

3.7. Ramachandran Plot Analysis

Ramachandran plot is an x-y plot of phi/psi dihedral angles between NC-alpha and Calpha-C bonds to evaluate a protein’s backbone conformation. The Ramachandran plot of the wild-type protein in the Swiss-Model showed 132 residues (93.6%) in the favored region, 8 residues in allowing region (5.7%), and one residue (0.7%) in the outer region (Figure S4A). The LEP hypothetical protein built by Phyre2 revealed 135 residues (95.7%) in favored region, 4 residues (2.8%) in allowing region, and 2 residues (1.4) in the outer part. Similarly, the Ramachandran plot of the wild LEP protein model constructed by the ConSurf server demonstrates 158 residues (95%) in the favored region, 4 residues (2.4%) in the allowing part, and 3 residues (1.8%) in the outer region. The last server we have been used for homology modelling was RaptorX. Here, the Ramachandran plot analysis of the model produced with this tool shows that the number of residues in the favored region was 52 (96.3%), the number of residues in allowed region 2 (3.7%), and the number of residues in outlier region was 0 (0.0%). Overall, the homology models evidenced the models of LEP protein were good. They can be used for further experiments and a better understanding of the protein’s biological activity.

4. Discussion

Missense mutations in genes have a profound impact on protein function. Extensive computational analysis of the phenotypic features attributed to nsSNPs may reveal the susceptible variants by interrupting the original function [37, 38]. In humans, mutations of the LEP gene have been associated with obesity in different populations [39]. The sheep LEP gene, located on chromosome 4, encodes 204 amino acids annotated with 11 domains and is currently associated with 543 variations [12].

Leptin is an endocrine hormone member of the long-chain helical cytokine family. It has multiple effects on regulating food intake, energy expenditure, body weight, and immune responses [40]. Polymorphism in the leptin gene has been suggested to be connected to cattle carcass composition differences [38, 41, 42]. For example, a nucleotide switched from cytosine (C) to thymine (T) causes an amino acid change in the plasma leptin circulation and a higher 12^th rib fat and marbling score [41]. In addition, the T-allele in cattle has been reported to be associated with fatter carcasses, whereas the C-allele is associated with leaner carcasses [43].

Numerous research studies have reported polymorphisms in the sheep leptin (LEP) gene [44, 45] and its association with food intake [46], growth traits [47, 48], and carcass and mutton quality traits [49, 50]. Consequently, amino acid substitutions in the LEP gene have been recommended as predictors of relative differences among individuals for such economic traits [51].

In this study, following the screening of nsSNPs as deleterious or neutral upon their effect on protein function using the SIFT, predictSNP2, PROVEAN, and SNAP2, five nsSNPs T86M (C → T, rs593507294), D98N (G → A, rs426762318), N136T (A → C, rs429690456), R142Q (G → A, rs409584889), and P157Q (C → A, rs1093355763) were predicted to have a significant impact on the protein structure, stability, and function by the majority of the tools. As described [37] in the human CDKN1A gene and ADIPOR2 gene by [52], and the CSN3 gene in cattle by [53], dbSNP-based studies have been reported. Our study employed similar tools to predict the impact of the nsSNPs in the LEP gene and found important results.

The sequence of a protein determines the protein’s physicochemical properties, such as protein structure, protein thermodynamic stability, the ability to interact with other molecules, and catalytic capacity [54, 55]. Thus, these properties, in turn, determine protein function [56].

Further analysis of these nsSNPs on protein functionality was predicted indirectly based on the effects exerted on protein stability. First, the I-stable server predicted that all the 11 nsSNPs impacted protein stability. Then, the report from I-stable was crosschecked using SDM and DynaMut. The stability analysis verified by SDM and DynaMut servers consistently narrowed to 5 nsSNPs. Three nsSNPs (T86M, D98N, and P157Q) had increased protein stability, whereas the other two nsSNPs (N136T and R142Q) showed a destabilizing effect on the wild-type proteins. DynaMut predicts the molecular consequences of the mutations in-depth through changes in protein dynamics and stability from vibrational entropy change [24]. The DynaMut-prediction outcomes of the Δ vibrational entropy energy between wild-type and mutants (ΔΔSVib ENCoM) for the 3 nsSNPs were − 0.046 kcal.mol⁻¹.K⁻¹, − 0.056 kcal.mol⁻¹.K⁻¹, and − 0.027 kcal.mol⁻¹.K⁻¹, respectively. Our data suggest that the vibrational entropy change of protein upon mutation of the amino acid decreases the molecular flexibility of the protein. Therefore, the three nsSNPs that increased the protein stability may increase evolvability by allowing a protein to accept a wider range of beneficial mutations while still folding to its innate structure [57]. Our findings are compatible with [58], who reported that the marginal protein thermostability is a consequence of the mutation-selection balance. Similarly, [57, 59] explained that mutants derived from highly stable variants of a given protein are more likely to retain their fold and are consequently more likely to develop novel or improved protein functions.

In contrast to the other two nsSNPs, N136T (A → C, rs429690456) and R142Q (G → A, rs409584889) were evidenced by reduced stability. Table 3 suggests that these missense mutations might lead to new functions of the leptin gene in sheep. This finding is analogous to [21] stating that the evolution of new enzymatic specificities is accompanied by a loss of the protein’s thermodynamic stability (DDG), implying a trade-off between acquiring new function and stability. Most likely, the mutations confer new functions that destabilize [57]. Thus, in this study, the mutations that cause a destabilizing effect on the thermodynamic stability of the protein through increasing the molecular flexibility may constrain the native evolution of leptin protein in sheep. Taken all together, the effect of the deleterious nsSNPs in the LEP gene may have a diverse impact on the performance of the sheep for the particular traits. The results from this study provide promising insight for future in vitro studies to specify the impact of these sorted nsSNPs on sheep performance.

5. Conclusion

A simulation-based study to detect the nsSNPs using many computational programs employed in this study has revealed 5 amino acid substitutions: T86M (C → T, rs593507294), D98N (G → A, rs426762318), N136T (A → C, rs429690456), R142Q (G → A, rs409584889), and P157Q (C → A, rs1093355763) were found to be functionally deleterious. These nsSNPs in the sheep LEP gene may contribute to the functional discrepancies compared to the wild LEP gene, which significantly regulates energy homeostasis, body weight control, and reproduction traits. Most of the variants identified in this study are located in the conserved domain of leptin. We also calculated the free energy changes for mutants and wild-type LEP proteins to evaluate their stability. Our data provide evidence for the functional role of the 5nsSNPs, which is helpful for further in vitro study to facilitate the development of reliable molecular markers for future practical application in sheep breeding.

Data Availability

All datasets generated/analyzed for this study are included in the manuscript and submitted as supplementary.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

S. G. and H.I. A conceived and designed the methodology, Q. A. Z. performed data analysis, and S. G. wrote the manuscript.

Acknowledgments

The authors are thankful to Samara University for providing encouragement and facilities and colleagues in the Department of Animal Science for their valuable comments and suggestions. The research was funded by Samara University, Ethiopia.

Supplementary Materials

Supplementary Materials Table S1. Deleterious SNPs predicted by predictSNP2. Table S2. nsSNP analysis by SNAP². Table S 3. nsSNP analysis by PROVEAN. Table S4. SDM and DynaMut tools analysis results of the effect of missense mutations on protein stability. Figure S1. Prediction of evolutionary conserved amino acid residues by ConSurf server. Conservation score is represented as the color coding bars. Figure S2. Protein secondary structure predictions by PSIPRED tool. The graphical output of PSIPRED prediction of secondary structure of the sheep LEP protein shows 6 α-helices extends from 41^th to 55^th, 66^thto 84^th, 110^th to 125^th, 133^th to 150^th, 164^th to 176^th, and 178^th to 198^th residue and no β-strands. Figure S3. Homology models from different servers; (a) homology modelling by Swiss-Model server; (b) homology modelling by Phyre-2 server; and (c) Homology modelling by ConSurf; homology modelling by RaptorX server. Figure S4. Ramachandran plots of different models; (a) LEP_ Swiss-Model, (b) LEP_Phyre2, (c) LEP_ConSurf, and (d) LEP_RaptorX. (Supplementary Materials)

References

M. W. Schwartz, S. C. Woods, D. Porte Jr., R. J. Seeley, and D. G. Baskin, “Central nervous system control of food intake,” Nature, vol. 404, no. 6778, pp. 661–671, 2000.
View at: Publisher Site | Google Scholar
H.-K. Park and R. S. Ahima, “Physiology of leptin: energy homeostasis, neuroendocrine function and metabolism,” Metabolism, vol. 64, no. 1, pp. 24–34, 2015.
View at: Publisher Site | Google Scholar
S. A. Robertson, G. M. Leinninger, and M. G. Myers Jr., “Molecular and neural mediators of leptin action,” Physiology & Behavior, vol. 94, no. 5, pp. 637–642, 2008.
View at: Publisher Site | Google Scholar
V. M. Rodríguez, C. Picó, M. P. Portillo, M. T. Macarulla, and A. Palou, “Dietary fat source regulates ob gene expression in white adipose tissue of rats under hyperphagic feeding,” British Journal of Nutrition, vol. 87, no. 5, pp. 427–434, 2002.
View at: Publisher Site | Google Scholar
L. Yang and G. S. McKnight, “Hypothalamic PKA regulates leptin sensitivity and adiposity,” Nature Communications, vol. 6, no. 1, p. 8237, 2015.
View at: Publisher Site | Google Scholar
M. G. Myers Jr., R. L. Leibel, R. J. Seeley, and M. W. Schwartz, “Obesity and leptin resistance: distinguishing cause from effect,” Trends in Endocrinology & Metabolism, vol. 21, no. 11, pp. 643–651, 2010.
View at: Publisher Site | Google Scholar
Y. Zhou and L. Rui, “Leptin signaling and leptin resistance,” Frontiers of Medicine, vol. 7, no. 2, pp. 207–222, 2013.
View at: Publisher Site | Google Scholar
C. Ding and S. Jin, “High-throughput methods for SNP genotyping,” Single Nucleotide Polymorphisms, Springer, pp. 245–254, 2009.
View at: Publisher Site | Google Scholar
H. Wohlrab, “The human mitochondrial transport/carrier protein family. Nonsynonymous single nucleotide polymorphisms (nsSNPs) and mutations that lead to human diseases,” Biochimica et Biophysica Acta (BBA)-Bioenergetics, vol. 1757, no. 9-10, pp. 1263–1270, 2006.
View at: Publisher Site | Google Scholar
L. V. Singh, S. Jayakumar, A. Sharma et al., “Comparative screening of single nucleotide polymorphisms in β-casein and κ-casein gene in different livestock breeds of India,” Meta Gene, vol. 4, pp. 85–91, 2015.
View at: Publisher Site | Google Scholar
K. Mahrous, M. Hassanane, H. Shafey, M. A. Mordy, and H. Rushdi, “Association between single nucleotide polymorphism in ovine Calpain gene and growth performance in three Egyptian sheep breeds,” Journal of Genetic Engineering and Biotechnology, vol. 14, no. 2, pp. 233–240, 2016.
View at: Publisher Site | Google Scholar
K. R. Rosenbloom, J. Armstrong, G. P. Barber et al., “The UCSC genome browser database: 2015 update,” Nucleic Acids Research, vol. 43, no. D1, pp. D670–D681, 2015.
View at: Google Scholar
D. R. Zerbino, P. Achuthan, W. Akanni et al., “Ensembl 2018,” Nucleic Acids Research, vol. 46, no. D1, pp. D754–D761, 2018.
View at: Publisher Site | Google Scholar
N.-L. Sim, P. Kumar, J. Hu, S. Henikoff, G. Schneider, and P. C. Ng, “SIFT web server: predicting effects of amino acid substitutions on proteins,” Nucleic Acids Research, vol. 40, no. W1, pp. W452–W457, 2012.
View at: Publisher Site | Google Scholar
P. C. Ng and S. Henikoff, “SIFT: Predicting amino acid changes that affect protein function,” Nucleic Acids Research, vol. 31, no. 13, pp. 3812–3814, 2003.
View at: Publisher Site | Google Scholar
J. Bendl, J. Stourac, O. Salanda et al., “PredictSNP: robust and accurate consensus classifier for prediction of disease-related mutations,” PLoS Computational Biology, vol. 10, no. 1, p. e1003440, 2014.
View at: Publisher Site | Google Scholar
Y. Choi, G. E. Sims, S. Murphy, J. R. Miller, and A. P. Chan, “Predicting the functional effect of amino acid substitutions and indels,” PloS One, vol. 7, no. 10, p. e46688, 2012.
View at: Publisher Site | Google Scholar
J. Turkson, “STAT proteins as novel targets for cancer drug discovery,” Expert Opinion on Therapeutic Targets, vol. 8, no. 5, pp. 409–422, 2004.
View at: Publisher Site | Google Scholar
M. Hecht, Y. Bromberg, and B. Rost, “Better prediction of functional effects for sequence variants,” BMC Genomics, vol. 16, no. 8, pp. 1–12, 2015.
View at: Google Scholar
N. Ota, R. Kurahashi, S. Sano, and K. Takano, “The direction of protein evolution is destined by the stability,” Biochimie, vol. 150, pp. 100–109, 2018.
View at: Publisher Site | Google Scholar
N. Tokuriki, F. Stricher, L. Serrano, and D. S. Tawfik, “How protein stability and new functions trade off,” PLoS Computational Biology, vol. 4, no. 2, p. e1000002, 2008.
View at: Publisher Site | Google Scholar
C.-W. Chen, J. Lin, and Y.-W. Chu, “iStable: off-the-shelf predictor integration for predicting protein stability changes,” BMC Bioinformatics, vol. 14, no. S2, 2013.
View at: Publisher Site | Google Scholar
C. L. Worth, R. Preissner, and T. L. Blundell, “SDM—a server for predicting effects of mutations on protein stability and malfunction,” Nucleic Acids Research, vol. 39, suppl, pp. W215–W222, 2011.
View at: Publisher Site | Google Scholar
C. H. Rodrigues, D. E. Pires, and D. B. Ascher, “DynaMut: predicting the impact of mutations on protein conformation, flexibility and stability,” Nucleic Acids Research, vol. 46, no. W1, 2018.
View at: Google Scholar
H. Ashkenazy, S. Abadi, E. Martz et al., “ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules,” Nucleic Acids Research, vol. 44, no. W1, pp. W344–W350, 2016.
View at: Publisher Site | Google Scholar
H. I. Ahmad, N. Ijaz, G. Afzal et al., “Computational Insights into the Structural and Functional Impacts of nsSNPs of bone morphogenetic proteins,” BioMed Research International, vol. 2022, Article ID 4013729, 17 pages, 2022.
View at: Publisher Site | Google Scholar
H. I. Ahmad, G. Afzal, M. N. Iqbal et al., “Positive selection drives the adaptive evolution of mitochondrial antiviral signaling (MAVS) proteins-mediating innate immunity in mammals,” Frontiers in Veterinary Science, vol. 8, 2021.
View at: Publisher Site | Google Scholar
D. W. Buchan, F. Minneci, T. C. Nugent, K. Bryson, and D. T. Jones, “Scalable web services for the PSIPRED protein analysis workbench,” Nucleic Acids Research, vol. 41, no. W1, pp. W349–W357, 2013.
View at: Publisher Site | Google Scholar
D. T. Jones, “Protein secondary structure prediction based on position-specific scoring matrices 1 1Edited by G. Von Heijne,” Journal of Molecular Biology, vol. 292, no. 2, pp. 195–202, 1999.
View at: Publisher Site | Google Scholar
K. Arnold, L. Bordoli, J. Kopp, and T. Schwede, “The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling,” Bioinformatics, vol. 22, no. 2, pp. 195–201, 2006.
View at: Publisher Site | Google Scholar
H. I. Ahmad, G. Afzal, S. Sadia et al., “Structural and evolutionary adaptations of Nei-like DNA glycosylases proteins involved in base excision repair of oxidative DNA damage in vertebrates,” Oxidative Medicine and Cellular Longevity, vol. 2022, 20 pages, 2022.
View at: Publisher Site | Google Scholar
M. Biasini, S. Bienert, A. Waterhouse et al., “SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information,” Nucleic Acids Research, vol. 42, no. W1, pp. W252–W258, 2014.
View at: Publisher Site | Google Scholar
L. A. Kelley, S. Mezulis, C. M. Yates, M. N. Wass, and M. J. Sternberg, “The Phyre2 web portal for protein modeling, prediction and analysis,” Nature Protocols, vol. 10, no. 6, pp. 845–858, 2015.
View at: Publisher Site | Google Scholar
H. I. Ahmad, A. R. Asif, M. J. Ahmad et al., “Adaptive evolution of peptidoglycan recognition protein family regulates the innate signaling against microbial pathogens in vertebrates,” Microbial Pathogenesis, vol. 147, p. 104361, 2020.
View at: Publisher Site | Google Scholar
M. Källberg, H. Wang, S. Wang et al., “Template-based protein structure modeling using the RaptorX web server,” Nature Protocols, vol. 7, no. 8, pp. 1511–1522, 2012.
View at: Publisher Site | Google Scholar
S. C. Lovell, I. W. Davis, W. B. Arendall III et al., “Structure validation by Cα geometry: ϕ, ψ and Cβ deviation,” Proteins: Structure, Function, and Bioinformatics, vol. 50, no. 3, pp. 437–450, 2003.
View at: Publisher Site | Google Scholar
G. S. H. Prince and T. Dhar, A simulation analysis and screening of deleterious non-synonymous single nucleotide polymorphisms (SNPs) in human CDKN1A gene, bioRxiv, 2018.
View at: Publisher Site
T. Geary, E. McFadin, M. MacNeil et al., “Leptin as a predictor of carcass composition in beef cattle,” Journal of Animal Science, vol. 81, no. 1, pp. 1–8, 2003.
View at: Publisher Site | Google Scholar
H. Bouafi, S. Bencheikh, M. A. L. Krami et al., “Prediction and structural comparison of deleterious coding nonsynonymous single nucleotide polymorphisms (nsSNPs) in human LEP gene associated with obesity,” Biomed Research International, vol. 2019, 10 pages, 2019.
View at: Publisher Site | Google Scholar
R. Nogueiras, M. H. Tschöp, and J. M. Zigman, “Central nervous system regulation of energy metabolism: ghrelin versus leptin,” Annals of the New York Academy of Sciences, vol. 1126, no. 1, pp. 14–19, 2008.
View at: Publisher Site | Google Scholar
F. C. Buchanan, C. J. Fitzsimmons, A. G. Van Kessel, T. D. Thue, D. C. Winkelman-Sim, and S. M. Schmutz, “Association of a missense mutation in the bovine leptin gene with carcass fat content and leptin mRNA levels,” Genetics Selection Evolution, vol. 34, no. 1, p. 105, 2002.
View at: Publisher Site | Google Scholar
T. Yamada, S. I. Kawakami, and N. Nakanishi, “The relationship between plasma leptin concentrations and the distribution of body fat in crossbred steers,” Animal Science Journal, vol. 74, no. 2, pp. 95–100, 2003.
View at: Publisher Site | Google Scholar
B. Kemp, Utilization of the cytosine to thymine missense mutation in the bovine leptin gene, RAK Genetic Consulting Ltd, 2003.
A. Hashemi, K. Mardani, M. Farhadian, I. Ashrafi, and M. Ranjbari, “Allelic polymorphism of Makoei sheep leptin gene identified by polymerase chain reaction and single strand conformation polymorphism,” African Journal of Biotechnology, vol. 10, no. 77, pp. 17903–17906, 2011.
View at: Google Scholar
H. Zhou, J. G. Hickford, and H. Gong, “Identification of allelic polymorphism in the ovine leptin gene,” Molecular Biotechnology, vol. 41, no. 1, pp. 22–25, 2009.
View at: Publisher Site | Google Scholar
M. Marie, P. Findlay, L. Thomas, and C. L. Adam, “Daily patterns of plasma leptin in sheep: effects of photoperiod and food intake,” Journal of Endocrinology, vol. 170, no. 1, pp. 277–286, 2001.
View at: Publisher Site | Google Scholar
A. Hajihosseinlo, A. Hashemi, and S. Sadeghi, “Association between polymorphism in exon 3 of leptin gene and growth traits in the Makooei sheep of Iran,” Livestock Research for Rural Development, vol. 24, no. 9, pp. 543–546, 2012.
View at: Google Scholar
M. Shojaei, M. Mohammad Abadi, M. Asadi Fozi, O. Dayani, A. Khezri, and M. Akhondi, “Association of growth trait and Leptin gene polymorphism in Kermani sheep,” Journal of Cell and Molecular Research, vol. 2, no. 2, pp. 67–73, 2011.
View at: Google Scholar
R. Barzehkar, A. Salehi, and F. Mahjoubi, “Polymorphisms of the ovine leptin gene and its association with growth and carcass traits in three Iranian sheep breeds,” Iranian Journal of Biotechnology, vol. 7, no. 4, pp. 241–246, 2009.
View at: Google Scholar
D. Boucher, M. Palin, F. Castonguay, C. Gariépy, and F. Pothier, “Detection of polymorphisms in the ovine leptin (LEP) gene: association of a single nucleotide polymorphism with muscle growth and meat quality traits,” Canadian Journal of Animal Science, vol. 86, no. 1, pp. 31–35, 2006.
View at: Google Scholar
A. Vignal, D. Milan, M. SanCristobal, and A. Eggen, “A review on SNP and other types of molecular markers and their use in animal genetics,” Genetics Selection Evolution, vol. 34, no. 3, pp. 275–305, 2002.
View at: Publisher Site | Google Scholar
M. Solayman, M. A. Saleh, S. Paul, M. I. Khalil, and S. H. Gan, “In silico analysis of nonsynonymous single nucleotide polymorphisms of the human adiponectin receptor 2 (ADIPOR2) gene,” Computational Biology and Chemistry, vol. 68, pp. 175–185, 2017.
View at: Publisher Site | Google Scholar
J. B. Patel and J. B. Chauhan, “Computational analysis of non-synonymous single nucleotide polymorphism in the bovine cattle kappa-casein (CSN3) gene,” Meta Gene, vol. 15, pp. 1–9, 2018.
View at: Publisher Site | Google Scholar
C. B. Anfinsen, “Principles that govern the folding of protein chains,” Science, vol. 181, no. 4096, pp. 223–230, 1973.
View at: Publisher Site | Google Scholar
H. I. Ahmad, J. Zhou, M. J. Ahmad et al., “Adaptive selection in the evolution of programmed cell death-1 and its ligands in vertebrates,” Aging (Albany NY), vol. 12, no. 4, pp. 3516–3557, 2020.
View at: Publisher Site | Google Scholar
C. L. Araya, D. M. Fowler, W. Chen, I. Muniez, J. W. Kelly, and S. Fields, “A fundamental protein property, thermodynamic stability, revealed solely from large-scale measurements of protein function,” Proceedings of the National Academy of Sciences, vol. 109, no. 42, pp. 16858–16863, 2012.
View at: Publisher Site | Google Scholar
J. D. Bloom, S. T. Labthavikul, C. R. Otey, and F. H. Arnold, “Protein stability promotes evolvability,” Proceedings of the National Academy of Sciences, vol. 103, no. 15, pp. 5869–5874, 2006.
View at: Publisher Site | Google Scholar
J. D. Bloom, A. Raval, and C. O. Wilke, “Thermodynamics of neutral protein evolution,” Genetics, vol. 175, no. 1, pp. 255–266, 2007.
View at: Publisher Site | Google Scholar
S. Bershtein, K. Goldin, and D. S. Tawfik, “Intense neutral drifts yield robust and evolvable consensus proteins,” Journal of Molecular Biology, vol. 379, no. 5, pp. 1029–1044, 2008.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Shishay Girmay et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

339

Downloads

431

Citations