In silico analysis and high-risk pathogenic phenotype predictions of non-synonymous single nucleotide polymorphisms in human Crystallin beta A4 gene associated with congenital cataract

Molecular Dynamics Simulation of Kir6.2 Variants Reveals Potential Association with Diabetes Mellitus

Diabetes mellitus (DM) represents a problem for the healthcare system worldwide. DM has very serious complications such as blindness, kidney failure, and cardiovascular disease. In addition to the very bad socioeconomic impacts, it influences patients and their families and communities. The global costs of DM and its complications are huge and expected to rise by the year 2030. DM is caused by genetic and environmental risk factors. Genetic testing will aid in early diagnosis and identification of susceptible individuals or populations using ATP-sensitive potassium (KATP) channels present in different tissues such as the pancreas, myocardium, myocytes, and nervous tissues. The channels respond to different concentrations of blood sugar, stimulation by hormones, or ischemic conditions. In pancreatic cells, they regulate the secretion of insulin and glucagon. Mutations in the KCNJ11 gene that encodes the Kir6.2 protein (a major constituent of KATP channels) were reported to be associated with Type 2 DM, neonatal diabetes mellitus (NDM), and maturity-onset diabetes of the young (MODY). Kir6.2 harbors binding sites for ATP and phosphatidylinositol 4,5-diphosphate (PIP2). The ATP inhibits the KATP channel, while the (PIP2) activates it. A Kir6.2 mutation at tyrosine330 (Y330) was demonstrated to reduce ATP inhibition and predisposes to NDM. In this study, we examined the effect of mutations on the Kir6.2 structure using bioinformatics tools and molecular dynamic simulations (SIFT, PolyPhen, SNAP2, PANTHER, PhD&SNP, SNP&Go, I-Mutant, MuPro, MutPred, ConSurf, HOPE, and GROMACS). Our results indicated that M199R, R201H, R206H, and Y330H mutations influence Kir6.2 structure and function and therefore may cause DM. We conclude that MD simulations are useful techniques to predict the effects of mutations on protein structure. In addition, the M199R, R201H, R206H, and Y330H variant in the Kir6.2 protein may be associated with DM. These results require further verification in protein-protein interactions, Kir6.2 function, and case-control studies.

Identification and In-Silico study of non-synonymous functional SNPs in the human SCN9A gene

Single nucleotide polymorphisms are the most common form of DNA alterations at the level of a single nucleotide in the genomic sequence. Genome-wide association studies (GWAS) were carried to identify potential risk genes or genomic regions by screening for SNPs associated with disease. Recent studies have shown that SCN9A comprises the NaV1.7 subunit, Na+ channels have a gene encoding of 1988 amino acids arranged into 4 domains, all with 6 transmembrane regions, and are mainly found in dorsal root ganglion (DRG) neurons and sympathetic ganglion neurons. Multiple forms of acute hypersensitivity conditions, such as primary erythermalgia, congenital analgesia, and paroxysmal pain syndrome have been linked to polymorphisms in the SCN9A gene. Under this study, we utilized a variety of computational tools to explore out nsSNPs that are potentially damaging to heath by modifying the structure or activity of the SCN9A protein. Over 14 potentially damaging and disease-causing nsSNPs (E1889D, L1802P, F1782V, D1778N, C1370Y, V1311M, Y1248H, F1237L, M936V, I929T, V877E, D743Y, C710W, D623H) were identified by a variety of algorithms, including SNPnexus, SNAP-2, PANTHER, PhD-SNP, SNP & GO, I-Mutant, and ConSurf. Homology modeling, structure validation, and protein-ligand interactions also were performed to confirm 5 notable substitutions (L1802P, F1782V, D1778N, V1311M, and M936V). Such nsSNPs may become the center of further studies into a variety of disorders brought by SCN9A dysfunction. Using in-silico strategies for assessing SCN9A genetic variations will aid in organizing large-scale investigations and developing targeted therapeutics for disorders linked to these variations.

Revisiting of Properties and Modified Polyethylenimine-Based Cancer Gene Delivery Systems

A new era of medical technology in cancer treatment is a directly specific modification of gene expression in tumor cells by nucleic acid delivery. Currently, the main challenge to achieving this goal is to find a non-toxic, safe, and effective strategy for gene transfer to cancer cells. Synthetic composites based on cationic polymers have historically been favored in bioengineering due to their ability to mimic bimolecular structures. Among them, polyethylenimines (PEIs) with superior properties such as a wide range of molecular weight and a flexible structure may propel the development of functional combinations in the biomedical and biomaterial fields. Here, in this review, we will focus on the recent progressions in the formulation optimization of PEI-based polyplex in gene delivery to treat cancer. Also, the effect of PEI's intrinsic characteristics such as structure, molecular weight, and positive charges which influence the gene delivery efficiency will be discussed.

Computational analysis of the functional and structural impact of the most deleterious missense mutations in the human Protein C

Protein C (PC) is a vitamin K-dependent factor that plays a crucial role in controlling anticoagulant processes and acts as a cytoprotective agent to promote cell survival. Several mutations in human PC are associated with decreased protein production or altered protein structure, resulting in PC deficiency. In this study, we conducted a comprehensive analysis of nonsynonymous single nucleotide polymorphisms in human PC to prioritize and confirm the most high-risk mutations predicted to cause disease. Of the 340 missense mutations obtained from the NCBI database, only 26 were classified as high-risk mutations using various bioinformatic tools. Among these, we identified that 12 mutations reduced the stability of protein, and thereby had the greatest potential to disturb protein structure and function. Molecular dynamics simulations revealed moderate alterations in the structural stability, flexibility, and secondary structural organization of the serine protease domain of human PC for five missense mutations (L305R, W342C, G403R, V420E, and W444C) when compared to the native structure that could maybe influence its interaction with other molecules. Protein-protein interaction analyses demonstrated that the occurrence of these five mutations can affect the regular interaction between PC and activated factor V. Therefore, our findings assume that these mutants can be used in the identification and development of therapeutics for diseases associated with PC dysfunction, although assessment the effect of these mutations need to be proofed in in-vitro.

Identification of potential mutational hotspots in serratiopeptidase to address its poor pH tolerance issue

Serratiopeptidase is the multifunctionality metalloendopeptidase extensively employed in biopharmaceutical and industrial biotechnology. Despite its poor pH tolerance, serratiopeptidase must withstand the highly acidic environment of the gastrointestinal tract to be used as a potent anti-inflammatory and analgesic medication. In earlier studies, post-translational deamination related mutations showed alteration in the net charge of protein's surface. Therefore, the current study aimed to enhance the acid resistance of serratiopeptidase via implementing computational interventions to screen out the most stable mutational hotspot. The methodology used in this study is as follows: (a) Higher accessibility to surface (b) 4 Å away from active site region to avoid interference with its proteolytic activity, and (c) By converting non-conserved amide residues to acidic residues. A docking study has been conducted to establish the substrate specificity and binding affinity to native and mutant proteins. The docking outcomes were then validated using molecular dynamic simulations to clarify each mutant's molecular stability and conformation while preserving their activity. The results showed that N412D is the best-screened mutant with negative electrostatic potential that can alter the overall charge on the protein's surface with increased H+ ions. Alteration in overall charge leads to protein surface more acidic that causes a common ion effect in stomach pH and act as a buffer which could stabilize the serratiopeptidase amid extreme pH.Communicated by Ramaswamy H. Sarma.

Predicting the effects of rare genetic variants on oncogenic signaling pathways: A computational analysis of HRAS protein function

The HRAS gene plays a crucial role in regulating essential cellular processes for life, and this gene's misregulation is linked to the development of various types of cancers. Nonsynonymous single nucleotide polymorphisms (nsSNPs) within the coding region of HRAS can cause detrimental mutations that disrupt wild-type protein function. In the current investigation, we have employed in-silico methodologies to anticipate the consequences of infrequent genetic variations on the functional properties of the HRAS protein. We have discovered a total of 50 nsSNPs, of which 23 were located in the exon region of the HRAS gene and denoting that they were expected to cause harm or be deleterious. Out of these 23, 10 nsSNPs ([G60V], [G60D], [R123P], [D38H], [I46T], [G115R], [R123G], [P11OL], [A59L], and [G13R]) were identified as having the most delterious effect based on results of SIFT analysis and PolyPhen2 scores ranging from 0.53 to 69. The DDG values -3.21 kcal/mol to 0.87 kcal/mol represent the free energy change associated with protein stability upon mutation. Interestingly, we identified that the three mutations (Y4C, T58I, and Y12E) were found to improve the structural stability of the protein. We performed molecular dynamics (MD) simulations to investigate the structural and dynamic effects of HRAS mutations. Our results showed that the stable model of HRAS had a significantly lower energy value of -18756 kj/mol compared to the initial model of -108915 kj/mol. The RMSD value for the wild-type complex was 4.40 Å, and the binding energies for the G60V, G60D, and D38H mutants were -107.09 kcal/mol, -109.42 kcal/mol, and -107.18 kcal/mol, respectively as compared to wild-type HRAS protein had -105.85 kcal/mol. The result of our investigation presents convincing corroboration for the potential functional significance of nsSNPs in augmenting HRAS expression and adding to the activation of malignant oncogenic signalling pathways.

ADGRL3 genomic variation implicated in neurogenesis and ADHD links functional effects to the incretin polypeptide GIP

Attention deficit/hyperactivity disorder (ADHD) is the most common childhood neurodevelopmental disorder. Single nucleotide polymorphisms (SNPs) in the Adhesion G Protein-Coupled Receptor L3 (ADGRL3) gene are associated with increased susceptibility to developing ADHD worldwide. However, the effect of ADGRL3 non-synonymous SNPs (nsSNPs) on the ADGRL3 protein function is vastly unknown. Using several bioinformatics tools to evaluate the impact of mutations, we found that nsSNPs rs35106420, rs61747658, and rs734644, previously reported to be associated and in linkage with ADHD in disparate populations from the world over, are predicted as pathogenic variants. Docking analysis of rs35106420, harbored in the ADGLR3-hormone receptor domain (HRM, a common extracellular domain of the secretin-like GPCRs family), showed that HRM interacts with the Glucose-dependent insulinotropic polypeptide (GIP), part of the incretin hormones family. GIP has been linked to the pathogenesis of diabetes mellitus, and our analyses suggest a potential link to ADHD. Overall, the comprehensive application of bioinformatics tools showed that functional mutations in the ADGLR3 gene disrupt the standard and wild ADGRL3 structure, most likely affecting its metabolic regulation. Further in vitro experiments are granted to evaluate these in silico predictions of the ADGRL3-GIP interaction and dissect the complexity underlying the development of ADHD.

Association of single nucleotide polymorphism variations in CRYAA and CRYAB genes with congenital cataract in Pakistani population

Background

The purpose of present study was to analyze the association of single nucleotide polymorphism (SNPs) variant in CRYAA and CRYAB genes with Congenital Cataract.

Method

Total 196 blood samples of children were collected, out of which 102 samples were congenital cataract (case group) and 94 samples were normal individuals (control group). Genomic DNA was extracted by using optimized inorganic method. Tetra primers for SNPs were designed and TETRA-ARMs assay was performed on both groups. Genotypic, allelic frequency and haplotype analyses were obtained by using SNPstats software.

Results

The coordination of genotypic and allelic frequencies of CRYAA and CRYAB genes variants and the association between case and control groups showed increased risk of congenital cataract in children who contained rs13053109 G > C variant of CRYAA in all models (all P > 0.05). This depicts the evident difference between the frequencies of case and control groups. The haplotype analysis of SNPs rs3761382, rs7278468 and rs13051039 of CRYAA gene showed weak linkage disequilibrium between the 3 SNPs (r2 < 0.8). The haplotype CTC indicated the high risk of congenital cataract in infants based of its p value (OR = 1.60 95% CI = 0.11-22.64, P > 0.05).

Conclusion

The variation in CRYAA gene can be the risk factor for congenital cataract in infants.

In silico study of missense variants of FANCA, FANCC and FANCG genes reveals high risk deleterious alleles predisposing to Fanconi anemia pathogenesis

Among the 22 Fanconi anemia (FA) reported genes, 90% of mutational spectra were found in three genes, namely FANCA (64%), FANCC (12%) and FANCG (8%). Therefore, this study aimed to identify the high-risk deleterious variants in three selected genes (FANCA, FANCC, and FANCG) through various computational approaches. The missense variant datasets retrieved from the UCSC genome browser were analyzed for their pathogenicity, stability, and phylogenetic conservancy. A total of 23 alterations, of which 16 in FANCA, 6 in FANCC and one variant in FANCG, were found to be highly deleterious. The native and mutant structures were generated, which demonstrated a profound impact on the respective proteins. Besides, their pathway analysis predicted many other pathways in addition to the Fanconi anemia pathway, homologous recombination, and mismatch repair pathways. Hence, this is the first comprehensive study that can be useful for understanding the genetic signatures in the development of FA.

Antibacterial kinetics and phylogenetic analysis of Aloe vera plants

Uncontrolled use of antibiotics has resulted in the emergence of resistant bacteria. It has necessitated the evaluation of antibacterial activities and phylo-diversity of Aloe vera (also called Aloe barbadensis) plants as antimicrobial agent in Nigeria. Biotyped enteric bacilli of 251 strains obtained from fecal samples of patients with various gastro-intestinal complications are profiled for antibiogram. Resistant biotypes were assayed for susceptibility to Aloe vera latex and further evaluated for time-kill kinetics and phylo-diversity. More than 30% of enteric bacilli, including Citrobacter freundii, Escherichia coli and Proteus mirabilis were resistant to cotrimoxazole, ciprofloxacin, and tetracycline respectively at MIC >16 µg/ml (p=0.004). Aloe vera latex significantly inhibited 39.5% resistant enteric biotypes with a significant average reduction of the viable count at 1xMIC and 2xMIC to less than 3.0 Log10CFU/mL after 24 hours. Flavonoids, alkaloids, terpenoids and anthraquinine in anti-enteric sap significantly correlated and regressed with antibacterial activity (p<0.05), while two of the antimicrobial Aloe vera plants showed phylogenetic relatedness with other homologous. Anti-bacteria efficacy of some Nigerian Aloe vera latex could provide alternative therapy, while its phylo-diversity and genomic profiling would offer a promising avenue for identification and development of antimicrobial agents as drug candidates for natural antibiotics.

Genetic analysis of 39 erythrocytosis and hereditary hemochromatosis-associated genes in the Slovenian family with idiopathic erythrocytosis

BACKGROUND

Erythrocytosis is a condition with an excessive number of erythrocytes, accompanied by an elevated haemoglobin and/or haematocrit value. Congenital erythrocytosis has a diverse genetic background with several genes involved in erythropoiesis. In clinical practice, nine genes are usually examined, but in approximately 70% of patients, no causative mutation can be identified. In this study, we screened 39 genes, aiming to identify potential disease-driving variants in the family with erythrocytosis of unknown cause.

PATIENTS AND METHODS

Two affected family members with elevated haemoglobin and/or haematocrit and negative for acquired causes and one healthy relative from the same family were selected for molecular-genetic analysis of 24 erythrocytosis and 15 hereditary haemochromatosis-associated genes with targeted NGS. The identified variants were further analysed for pathogenicity using various bioinformatic tools and review of the literature.

RESULTS

Of the 12 identified variants, two heterozygous variants, the missense variant c.471G>C (NM_022051.2) (p.(Gln157His)) in the EGLN1 gene and the intron variant c.2572-13A>G (NM_004972.3) in the JAK2 gene, were classified as low-frequency variants in European population. None of the two variants were present in a healthy family member. Variant c.2572-13A>G has potential impact on splicing by one prediction tool.

CONCLUSION

For the first time, we included 39 genes in the erythrocytosis clinical panel and identified two potential disease-driving variants in the Slovene family studied. Based on the reported functional in vitro studies combined with our bioinformatics analysis, we suggest further functional analysis of variant in the JAK2 gene and evaluation of a cumulative effect of both variants.