Molecular dynamics, residue network analysis, and cross-correlation matrix to characterize the deleterious missense mutations in GALE causing galactosemia III

Mutations in LIFR rewire the JAK/STAT signaling pathway: A study unveiling mechanistic details of Stüve-Wiedemann syndrome

Stüve-Wiedemann syndrome (SWS), a rare autosomal recessive disorder, characterized by diminutive size, curvature of the elongated bones, bent fingers, episodes of heightened body temperature, respiratory distress or periods of breath-holding, and challenges with feeding, especially causes fatality in infants. SWS is an outcome of potential missense mutations in the leukemia inhibitory factor receptor gene reflected as numerous amino acid mutations at protein level. Employing in silico tools and techniques like mutational screening with Pred_MutHTP, I-Mutant2.0, PANTHER.db, PolyPhen, to classify mutations as deleterious/destabilizing, in conjunction with experimental data analysis, P136A and S279P emerged as 'effect'-causing mutations. Pre-existing knowledge suggests, SWS progression is effectuated conformationally altered and dysfunctional LIFR, unable to bind to LIF and further form the LIF/LIFR/gp130 signalling complex. To gain functional insights into the effect of the said mutations on the wild type protein, an all-atom, explicit, solvent molecular dynamics simulation was performed following docking approaches. Consequently, referring to the RMSD, RMSF, protein dynamic network analysis, energy landscape plots and domain motion analysis, it was revealed that unbound LIFR_WT was more prone to LIF binding as usual whereas the mutants exhibited considerable domain closure to inhibit LIF binding. We conducted binding affinity analysis via MM/GBSA and dissociation constant estimation after LIFR-LIF docking and found the WT_complex to be more stable and compact as a whole when compared to the flexible mutant complexes thus being associated with SWS. Our study offers a route for understanding molecular level implications upon LIFR mutations which opens an avenue for therapeutic interventions.

Molecular Dynamics Simulation Studies of Beta-Glucogallin and Dihydro Dehydro Coniferyl Alcohol from Syzygium cumini for its Antimicrobial Activity on Staphylococcus aureus

With the escalating threat of antimicrobial resistance (AMR), discovering novel therapeutic agents against resistant pathogens like Staphylococcus aureus is crucial. This study explores phytochemicals from Syzygium cumini for their potential efficacy against AMR S. aureus infections, elucidating their mechanisms through in silico methods. We investigated 83 compounds from S. cumini, sourced from PubMed, using rigorous docking analysis against the ATP binding domain AgrC of S. aureus with AMdock with Autodock Vina v1.5.2. Drug-likeness predictions were assessed using SwissADME v2023 and Pass online v2.0. Molecular dynamics (MD) simulations identified promising compounds, focusing on stability and interaction dynamics. Beta-Glucogallin (BEG) and Dihydro Dehydro Coniferyl alcohol (DIH) emerged as significant hits. MD simulations with GROMACS v2020.6 revealed stable BEG and DIH complexes with AgrC, forming six hydrogen bonds with six key amino acids (Arg-303, Asp-338, Glu-342, Glu-384, Lys-389, Gly-396), indicating strong and stable bonding. The binding affinities for DIH and BEG are -73.474 ± 11.104 kJ/mol and -6.319 ± 18.823 kJ/mol with 4BXI, respectively. Our findings highlight BEG and DIH as promising candidates against AMR S. aureus infections, showing favourable binding affinities and stable interactions with AgrC. This study underscores the importance of natural products in combating AMR and demonstrates the utility of computational methodologies in drug discovery. Further experimental validation is warranted to fully exploit these phytochemicals' therapeutic potential.

A new trivalent recombinant protein for type 2 diabetes mellitus with oral delivery potential: design, expression, and experimental validation

Glucagon-like peptide-1 (GLP-1) receptor agonists are increasingly used in clinical practice for the management of type 2 diabetes mellitus. However, the extremely short half-life of GLP-1 and the need for subcutaneous administration limit its clinical application. Thus, half-life extension and alternative delivery methods are highly desired. DARPin domains with high affinity for human serum albumin (HSA) have been selected for the half-life extension of therapeutic peptides and proteins. In the present study, novel trivalent fusion proteins as long-acting GLP-1 receptor agonists with potential for oral delivery were computationally engineered by incorporating a protease-resistant modified GLP-1, an anti-human serum albumin DARPin, and an approved cell-penetrating peptide (Penetratin, Tat, and Polyarginine) linked either by rigid or flexible linkers. Theoretical studies and molecular dynamics simulation results suggested that mGLP1-DARPin-Pen has acceptable quality and stability. Moreover, the potential affinity of the selected fusion proteins for GLP-1 receptor and human serum albumin was explored by molecular docking. The recombinant construct was cloned into the pET28a vector and expressed in Escherichia coli. SDS-PAGE analysis of the purified fusion protein matched its molecular size and was confirmed by western blot analysis. The results demonstrated that the engineered fusion protein could bind HSA with high affinity. Importantly, insulin secretion assays using a mouse pancreatic β-cell line (β-TC6) revealed that the engineered trivalent fusion protein retained the ability to stimulate cellular insulin secretion. Immunofluorescence microscopy analysis indicated the CPP-dependent cellular uptake of mGLP1-DARPin-Pen. These findings demonstrated that mGLP1-DARPin-Pen is a highly potent oral drug candidate that could be particularly useful in the treatment of type 2 diabetes mellitus.Communicated by Ramaswamy H. Sarma.

Structural effect of the H992D/H418D mutation of angiotensin-converting enzyme in the Indian population: implications for health and disease

The Non synonymous SNPs (nsSNPs) of the renin-angiotensin-system (RAS) pathway, unique to the Indian population were investigated in view of its importance as an endocrine system. nsSNPs of the RAS pathway genes were mined from the IndiGenome database. Damaging nsSNPs were predicted using SIFT, PredictSNP, SNP and GO, Snap2 and Protein Variation Effect Analyzer. Loss of function was predicted based on protein stability change using I mutant, PremPS and CONSURF. The structural impact of the nsSNPs was predicted using HOPE and Missense3d followed by modeling, refinement, and energy minimization. Molecular Dynamics studies were carried out using Gromacsv2021.1. 23 Indian nsSNPs of the RAS pathway genes were selected for structural analysis and 8 were predicted to be damaging. Further sequence analysis showed that HEMGH zinc binding motif changes to HEMGD in somatic ACE-C domain (sACE-C) H992D and Testis ACE (tACE) H418D resulted in loss of zinc coordination, which is essential for enzymatic activity in this metalloprotease. There was a loss of internal interactions around the zinc coordination residues in the protein structural network. This was also confirmed by Principal Component Analysis, Free Energy Landscape and residue contact maps. Both mutations lead to broadening of the AngI binding cavity. The H992D mutation in sACE-C is likely to be favorable for cardiovascular health, but may lead to renal abnormalities with secondary impact on the heart. H418D in tACE is potentially associated with male infertility.Communicated by Ramaswamy H. Sarma.

Exploring the natural products chemical space to abrogate the F3L-dsRNA interface of monkeypox virus to enhance the immune responses using molecular screening and free energy calculations

Amid the ongoing monkeypox outbreak, there is an urgent need for the rapid development of effective therapeutic interventions capable of countering the immune evasion mechanisms employed by the monkeypox virus (MPXV). The evasion strategy involves the binding of the F3L protein to dsRNA, resulting in diminished interferon (IFN) production. Consequently, our current research focuses on utilizing virtual drug screening techniques to target the RNA binding domain of the F3L protein. Out of the 954 compounds within the South African natural compound database, only four demonstrated notable docking scores: -6.55, -6.47, -6.37, and -6.35 kcal/mol. The dissociation constant (KD) analysis revealed a stronger binding affinity of the top hits 1-4 (-5.34, -5.32, -5.29, and -5.36 kcal/mol) with the F3L in the MPXV. All-atom simulations of the top-ranked hits 1 to 4 consistently exhibited stable dynamics, suggesting their potential to interact effectively with interface residues. This was further substantiated through analyses of parameters such as radius of gyration (Rg), Root Mean Square Fluctuation, and hydrogen bonding. Cumulative assessments of binding free energy confirmed the top-performing candidates among all the compounds, with values of -35.90, -52.74, -28.17, and -32.11 kcal/mol for top hits 1-4, respectively. These results indicate that compounds top hit 1-4 could hold significant promise for advancing innovative drug therapies, suggesting their suitability for both in vivo and in vitro experiments.

Discovery of novel HDAC8 inhibitors from natural compounds by in silico high throughput screening

A class I histone deacetylase HDAC8 is associated with several diseases, including cancer, intellectual impairment and parasite infection. Most of the HDAC inhibitors that have so far been found to inhibit HDAC8 limit their efficacy in the clinic by producing toxicities. It is therefore very desirable to develop specific HDAC8 inhibitors. The emergence of HDAC inhibitors derived from natural sources has become quite popular. In recent decades, it has been shown that naturally occurring HDAC inhibitors have strong anticancer properties. A total of 0.2 million natural compounds were screened against HDAC8 from the Universal Natural Product Database (UNPD). Molecular docking was performed for these natural compounds and the top six hits were obtained. In addition, molecular dynamics (MD) simulations were used to evaluate the structural stability and binding affinity of the inhibitors, which showed that the protein-ligand complexes remained stable throughout the 100 ns simulation. MM-PBSA method demonstrated that the selected compounds have high affinity towards HDAC8. We infer from our findings that Hit-1 (-29.35 kcal mol-1), Hit-2 (-29.15 kcal mol-1) and Hit-6 (-30.28 kcal mol-1) have better binding affinity and adhesion to ADMET (absorption, distribution, metabolism, excretion and toxicity) characteristics against HDAC8. To compare our discussions and result in an effective way. We performed molecular docking, MD and MM-PBSA analysis for the FDA-approved drug romidepsin. The above results show that our hits show better binding affinity than the compound romidepsin (-12.03 ± 4.66 kcal mol-1). The important hotspot residues Asp29, Ile34, Trp141, Phe152, Asp267, Met274 and Tyr306 have significantly contributed to the protein-ligand interaction. These findings suggest that in vitro testing and additional optimization may lead to the development of HDAC8 inhibitors.Communicated by Ramaswamy H. Sarma.

Computational refinement identifies functional destructive single nucleotide polymorphisms associated with human retinoid X receptor gene

Alterations in the nuclear retinoid X receptor (RXRs) signalling have been implicated in neurodegenerative disorders such as Alzheimer's disease, Parkinson's disease, stroke, multiple sclerosis and glaucoma. Single nucleotide polymorphisms (SNPs) are the main cause underlying single nucleic acid variations which in turn determine heterogeneity within various populations. These genetic polymorphisms have been suggested to associate with various degenerative disorders in population-wide analysis. This bioinformatics study was designed to investigate, search, retrieve and identify deleterious SNPs which may affect the structure and function of various RXR isoforms through a computational and molecular modelling approach. Amongst the 1,813 retrieved SNPs several were found to be deleterious with rs140464195_G139R, rs368400425_R358W and rs368586400_L383F RXRα mutant variants being the most detrimental ones causing changes in the interatomic interactions and decreasing the flexibility of the mutant proteins. Molecular genetics analysis identified seven missense mutations in RXRα/β/γ isoforms. Two novel mutations SNP IDs (rs1588299621 and rs1057519958) were identified in RXRα isoform. We used several in silico prediction tools such as SIFT, PolyPhen, I-Mutant, Protein Variation Effect Analyzer (PROVEAN), PANTHER, SNP&Go, PhD-SNP and SNPeffect to predict pathogenicity and protein stability associated with RXR mutations. The structural assessment by DynaMut tool revealed that hydrogen bonds were affected along with hydrophobic and carbonyl interactions resulting in reduced flexibility at the mutated residue positions but ultimately stabilizing the molecule as a whole. Summarizing, analysis of the missense mutations in RXR isoforms showed a mix of conclusive and inconclusive genotype-phenotype correlations suggesting the use of sophisticated computational analysis tools for studying RXR variants.Communicated by Ramaswamy H. Sarma.

An integrated investigation of structural and pathway alteration caused by PIK3CA and TP53 mutations identified in cfDNA of metastatic breast cancer

In peripheral blood, cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA), which indicates molecular abnormalities in metastatic breast tumor tissue. The sequencing of cfDNA of Metastatic Breast Cancer (MBC) patients allows assessment of therapy response and noninvasive treatment. In the proposed study, clinically significant alterations in PIK3CA and TP53 genes associated with MBC resulting in a missense substitution of His1047Arg and Arg282Trp from an next-generation sequencing-based multi-gene panel were reported in a cfDNA of a patient with MBC. To investigate the impact of the reported mutation, we used molecular docking, molecular dynamics simulation, network analysis, and pathway analysis. Molecular Docking analysis determined the distinct binding pattern revealing H1047R-ATP complex has a higher number of Hydrogen bonds (H-bonds) and binding affinity with a slight difference compared to the PIK3CA-ATP complex. Following, molecular dynamics simulation for 200 ns, of which H1047R-ATP complex resulted in the instability of PIK3CA. Similarly, for TP53 mutant R282W, the zinc-free state (apo) and zinc-bounded (holo) complexes were investigated for conformational change between apo and holo complexes, of which the holo complex mutant R282W was unstable. To validate the conformational change of PIK3CA and TP53, 80% mutation of H1047R in the kinase domain of p110α expressed ubiquitously in PIK3CA protein that alters PI3K pathway, while R282W mutation in DNA binding helix (H2) region of P53 protein inhibits the transcription factor in P53 pathway causing MBC. According to our findings, the extrinsic (hypoxia, oxidative stress, and acidosis); intrinsic factors (MYC amplification) in PIK3CA and TP53 mutations will provide potential insights for developing novel therapeutic methods for MBC therapy.

Structure-based virtual screening to identify potential lipase inhibitors to reduce lipid storage in Wolman disorder

Wolman disorder (WD) was first described in Iranian-Jewish (IJ) children, and it is caused by a deficiency of the lysosomal acid lipase (LAL). Newborns with WD are healthy and active at birth but soon develop severe malnutrition symptoms and often die before 1 year. In particular, spleens, livers, bone marrows, intestines, adrenal glands, and lymph nodes accumulate harmful amounts of lipids. G87V mutation in LIPA is responsible for Wolman disorder. Some reports suggest that δ-tocopherol can reduce lipid accumulation in cholesterol storage disorders. Hence, we used δ-tocopherol for the virtual screening process in this study. Initially, the lead compounds were docked with native and G87V mutant LIPA. Subsequently, the ADME and toxicity parameters for screened compounds were determined to ensure the safety profiles. Finally, the molecular dynamics simulations result indicated that dl-alpha-Tocopherol-13C3, a molecule obtained from the PubChem database, is identified as a potential and stable lead molecule that could be effective against the G87V mutant form of LIPA.

Novel variants in GALE cause syndromic macrothrombocytopenia by disrupting glycosylation and thrombopoiesis

Glycosylation is recognized as a key process for proper megakaryopoiesis and platelet formation. The enzyme uridine diphosphate (UDP)-galactose-4-epimerase, encoded by GALE, is involved in galactose metabolism and protein glycosylation. Here, we studied 3 patients from 2 unrelated families who showed lifelong severe thrombocytopenia, bleeding diathesis, mental retardation, mitral valve prolapse, and jaundice. Whole-exome sequencing revealed 4 variants that affect GALE, 3 of those previously unreported (Pedigree A, p.Lys78ValfsX32 and p.Thr150Met; Pedigree B, p.Val128Met; and p.Leu223Pro). Platelet phenotype analysis showed giant and/or grey platelets, impaired platelet aggregation, and severely reduced alpha and dense granule secretion. Enzymatic activity of the UDP-galactose-4-epimerase enzyme was severely decreased in all patients. Immunoblotting of platelet lysates revealed reduced GALE protein levels, a significant decrease in N-acetyl-lactosamine (LacNAc), showing a hypoglycosylation pattern, reduced surface expression of gylcoprotein Ibα-IX-V (GPIbα-IX-V) complex and mature β1 integrin, and increased apoptosis. In vitro studies performed with patients-derived megakaryocytes showed normal ploidy and maturation but decreased proplatelet formation because of the impaired glycosylation of the GPIbα and β1 integrin, and reduced externalization to megakaryocyte and platelet membranes. Altered distribution of filamin A and actin and delocalization of the von Willebrand factor were also shown. Overall, this study expands our knowledge of GALE-related thrombocytopenia and emphasizes the critical role of GALE in the physiological glycosylation of key proteins involved in platelet production and function.

Controlling cell proliferation by targeting cyclin-dependent kinase 6 using drug repurposing approach

Cyclin-dependent kinase 6 (CDK6) is an essential kinase in cell cycle progression, which is a viable target for inhibitors in various malignancies, including breast cancer. This study aimed to virtually screen efficient compounds as new leads in treating breast cancer using a drug repurposing approach. Apoptosis regulatory compounds were taken from the seleckchem database. Molecular docking experiments were carried out in the presence of abemaciclib, a routinely used FDA drug. Compared to conventional drugs, the two compounds demonstrated a higher binding affinity for CDK6. Compounds (N-benzyl-6-[(4-hydroxyphenyl)methyl]-8-(naphthalen-1-ylmethyl)-4,7-dioxo-3,6,9,9a-tetrahydro-2H-pyrazino[1,2-a]pyrimidine-1-carboxamide) and (1'-[4-[1-(4-fluorophenyl)indol-3-yl]butyl]spiro[1H-2-benzofuran-3,4'-piperidine]) were discovered to have an inhibitory effect against CDK6 at -8.49 and -6.78kcal/mol, respectively, compared to -8.09kcal/mol of the control molecule, the interacting residues of these two new compounds were found to fall within the binding site of the CDK6 molecule. Both compounds exhibited equal ADME features compared with abemaciclib and would be well distributed and metabolized by the body with an appropriate druglikeness range. Lastly, molecular dynamics was initiated for 200ns for the selected potent inhibitors and abemaciclib as complexed with CDK6. The RMSD, RMSF, Rg, H-Bond interactions, SASA, PCA, FEL, and MM/PBSA analysis were performed for the complexes to assess the stability, fluctuations, radius of gyration, hydrogen bond interaction, solvent accessibility, essential dynamics, free energy landscape, and MM/PBSA. The selected two compounds are small molecules in the appropriate druglikeness range. The results observed in molecular docking and molecular dynamics simulations were most promising for two compounds, suggesting their potent inhibitory effect against CDK6. We propose that these candidate compounds can undergo in vitro validation and in vivo testing for their further use against cancer.

Unraveling the Structural Changes in the DNA-Binding Region of Tumor Protein p53 (TP53) upon Hotspot Mutation p53 Arg248 by Comparative Computational Approach

The vital tissue homeostasis regulator p53 forms a tetramer when it binds to DNA and regulates the genes that mediate essential biological processes such as cell-cycle arrest, senescence, DNA repair, and apoptosis. Missense mutations in the core DNA-binding domain (109-292) simultaneously cause the loss of p53 tumor suppressor function and accumulation of the mutant p53 proteins that are carcinogenic. The most common p53 hotspot mutation at codon 248 in the DNA-binding region, where arginine (R) is substituted by tryptophan (W), glycine (G), leucine (L), proline (P), and glutamine (Q), is reported in various cancers. However, it is unclear how the p53 Arg248 mutation with distinct amino acid substitution affects the structure, function, and DNA binding affinity. Here, we characterized the pathogenicity and protein stability of p53 hotspot mutations at codon 248 using computational tools PredictSNP, Align GVGD, HOPE, ConSurf, and iStable. We found R248W, R248G, and R248P mutations highly deleterious and destabilizing. Further, we subjected all five R248 mutant-p53-DNA and wt-p53-DNA complexes to molecular dynamics simulation to investigate the structural stability and DNA binding affinity. From the MD simulation analysis, we observed increased RMSD, RMSF, and Rg values and decreased protein-DNA intermolecular hydrogen bonds in the R248-p53-DNA than the wt-p53-DNA complexes. Likewise, due to high SASA values, we observed the shrinkage of proteins in R248W, R248G, and R248P mutant-p53-DNA complexes. Compared to other mutant p53-DNA complexes, the R248W, R248G, and R248P mutant-p53-DNA complexes showed more structural alteration. MM-PBSA analysis showed decreased binding energies with DNA in all five R248-p53-DNA mutants than the wt-p53-DNA complexes. Henceforth, we conclude that the amino acid substitution of Arginine with the other five amino acids at codon 248 reduces the p53 protein's affinity for DNA and may disrupt cell division, resulting in a gain of p53 function. The proposed study influences the development of rationally designed molecular-targeted treatments that improve p53-based therapeutic outcomes in cancer.

Selectivity mechanism of muscarinic acetylcholine receptor antagonism through in silico investigation

Structures of muscarinic acetylcholine receptors illustrate the strikingly high degree of homology of the residues among isoforms, thus leading to difficulty in achieving subtype selectivity when targeting these receptors and causing undesired side effects when treating the corresponding diseases. Considering the urgent need for more selective and potency therapies, this study is aimed at revealing the selectivity mechanism against M4/5 via in silico strategies, revealing crucial molecular interactions such as hydrogen bonds and pi-cation interaction formed between the key residues TYR416, ASN417, and TRP435 of M4, respectively, hydrophobic pocket formed by the key residues, especially CYS484 of M5. Besides, the water around TYR416^M4 and ASN459^M5, which can be replaced by substituent groups which can form the hydrogen bond interaction network by simulated bridging water and the water around ASP112^M4, whose replacement maybe not contribute to the increase in binding affinity of the compound, may affect the inhibitory selectivity among M4/5 in the aspect of the solvent. Moreover, from the point of inhibitors, compounds with a positively ionizable group could selectively bind to M4 receptors, while hydrophobic molecules may bind preferably to M5. We believe that the current study would provide a basis for the design of subsequent M4/5 selective antagonists.

Deciphering the effect of mutations in MMAA protein causing methylmalonic acidemia-A computational approach

Methylmalonic acidemia (MMA) is a rare genetic disorder affecting multiple body systems. We aimed to investigate the pathogenic mutations in MMAA that are associated with isolated methylmalonic acidemia to identify the structural behavior of MMAA upon mutation. The algorithms such as PredictSNP, iStable, ConSurf, and Align GVGD were employed to analyze the consequence of the mutations. Molecular docking was carried out for the native MMAA, L89P, G274D, and R359G to interpret its interactions with the GDP substrate. The docked complexes were simulated for 200ns aiding GROMACS in apprehending the behavior of MMAA upon mutation and GDP binding. After simulation, cα disruptions were observed using the RMSF plot, which indicated that several regions of mutant MMAAs have highly fluctuated. The gyration and H-bond plots were used to understand the compactness and intermolecular interaction with the GDP molecule. The MDS analysis showed that the mutations L89P, G274D, and R359G are highly unstable even after GDP binding, with the least compactness, fewer H-bonds, and larger conformational cα motions. Our study provided structural and dynamic insights into MMAA protein, which further helps to characterize these mutants and provide potential treatment strategies for MMA patients.

Identification of potential inhibitors, conformational dynamics, and mechanistic insights into mutant Kirsten rat sarcoma virus (G13D) driven cancers

The mutations at the hotspot region of K-Ras result in the progression of cancer types. Our study aimed to explore the small molecule inhibitors against the G13D mutant K-Ras model with anti-cancerous activity from food and drug administration (FDA)-approved drug compounds. We implemented several computational strategies such as pharmacophore-based virtual screening, molecular docking, absorption, distribution, metabolism and excretion features, and molecular simulation to ensure the identified hit compounds have potential efficacy against G13D K-Ras. We found that the FDA-approved compounds, namely, azelastine, dihydrocodeine, paroxetine, and tramadol, are potential candidates to inhibit the action of G13D mutant K-Ras. All four compounds exhibited similar binding patterns of sotorasib, and a structural binding mechanism with significant hydrophobic contacts. The descriptor features from the QikProp of all four compounds are within allowable limits compared to sotorasib drug. Consequently, a molecular simulation result emphasized that the dihydrocodeine and tramadol exhibited less fluctuation, minimal basin, significant h-bonds, and potent inhibition against G13D K-Ras. As a result, the current research identifies prospective K-Ras inhibitors that could be further improved with biochemical analysis for precision medicine against K-Ras-driven cancers.

Mutations in G6PC2 gene with increased risk for development of type 2 diabetes: Understanding via computational approach

An increase in the fast blood glucose (FBG) levels has been linked to an increased risk of developing a chronic condition, type 2 diabetes (T2D). The mutation in the G6PC2 gene was identified to have a lead role in the modulation of FBG levels. The abnormal regulation of this enzyme influences glucose-stimulated insulin secretion (GSIS), which controls the insulin levels corresponding to the system's glucose level. This study focuses on the mutations at the G6PC2 gene, which cause the variation from normal expression levels and increase the risk of T2D. We examined the non-synonymous single nucleotide polymorphisms (nsSNPs) present in the G6PC2 and subjected them to pathogenicity, stability, residue conservation, and membrane simulation. The individual representation of surrounding amino acids in the mutant (I63T) model showed the loss of hydrophobic interactions compared to the native G6PC2. In addition, the trajectory results from the membrane simulation exhibited reduced stability, and the least compactness was identified for the I63T mutant model. Our study shed light on the structural and conformational changes at the transmembrane region due to the I63T mutation in G6PC2. Additionally, the Gibbs free energy landscape analysis against the two principal components showed structural differences and decreased the conformational stability of the I63T mutant model compared to the native. Like those presented in this study, dynamical simulations may indeed be crucial to comprehending the structural insights of G6PC2 mutations in cardiovascular-associated mortality and T2D.

Computational validation of ABCB1 gene polymorphism and its effect on tacrolimus dose concentration/levels in renal transplant individuals of South India

Patients with end-stage renal failure require hemodialysis and peritoneal dialysis; however, kidney transplantation is considered a better treatment option for renal failure patients, improving their quality of life and longevity. Among several potent immunosuppressive agents, tacrolimus (TAC) has shown progressive improvement in the graft survival rates after renal transplantation. Fifty kidney transplant patients undergoing TAC immunosuppressive treatment were included. The human genomic DNA was isolated using the phenol-chloroform extraction procedure. CYP3A5*6, CYP3A5*2, and ABCB1 exon 21 G2677 T/A polymorphisms were genotyped using the polymerase chain reaction-restriction fragment length polymorphism method. Fisher's exact test and Chi-square analysis were performed to analyze the data, where p < 0.05 was considered statistically significant. In addition, we implemented bioinformatics studies on ABCB1 protein to determine the mutation's effect sequentially and structurally. Among the genotyped single nucleotide polymorphisms (SNPs), SNPs of CYP3A5*2 and CYP3A5*6 did not vary in the studied population. The concentration/dose (C/D) ratio of TT genotype of the ABCB1 gene was higher (95% CI: 177.38-269.46) when compared to TA and AA. However, there were no substantial differences between the ABCB1 genotypes and TAC C/D ratio (p = 0.953). The TAC dose mg/kg/day (p = 0.002) and C/D ratio (p = 0.004) exhibited a statistically significant difference. However, no significant difference was found with respect to the ABCB1 gene between the non-toxicity and toxicity groups. Mutation and residue interaction analysis results showed that the S893T mutation destabilizes the ABCB1 protein, thus reducing the protein's flexibility. The present study demonstrated a substantial relationship between the TAC dose and C/D ratio, including the non-toxicity and toxicity groups. However, no possible correlation was observed between the ABCB1 gene polymorphism and renal transplant.

Computational Analysis of Missense Variants in the Human Transmembrane Protease Serine 2 (TMPRSS2) and SARS-CoV-2

The human transmembrane protease serine 2 (TMPRSS2) protein plays an important role in prostate cancer progression. It also facilitates viral entry into target cells by proteolytically cleaving and activating the S protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In the current study, we used different available tools like SIFT, PolyPhen2.0, PROVEAN, SNAP2, PMut, MutPred2, I-Mutant Suite, MUpro, iStable, ConSurf, ModPred, SwissModel, PROCHECK, Verify3D, and TM-align to identify the most deleterious variants and to explore possible effects on the TMPRSS2 stability, structure, and function. The six missense variants tested were evaluated to have deleterious effects on the protein by SIFT, PolyPhen2.0, PROVEAN, SNAP2, and PMut. Additionally, V160M, G181R, R240C, P335L, G432A, and D435Y variants showed a decrease in stability by at least 2 servers; G181R, G432A, and D435Y are highly conserved and identified posttranslational modifications sites (PTMs) for proteolytic cleavage and ADP-ribosylation using ConSurf and ModPred servers. The 3D structure of TMPRSS2 native and mutants was generated using 7 meq as a template from the SwissModeller group, refined by ModRefiner, and validated using the Ramachandran plot. Hence, this paper can be advantageous to understand the association between these missense variants rs12329760, rs781089181, rs762108701, rs1185182900, rs570454392, and rs867186402 and susceptibility to SARS-CoV-2.

Computational investigation to identify potent inhibitors of the GTPase-Kirsten RAt sarcoma virus (K-Ras) mutants G12C and G12D

K-Ras mutations are frequent in various cancer types, and according to recent research, K-Ras possesses four-drug targeting sites. This increased our interest in finding potential small molecule inhibitors with anticancer activity to treat K-Ras-driven cancers. We utilized integrated bioinformatic strategies, such as XP docking, MM-GBSA, cell-line cytotoxicity prediction, ADMET, and molecular simulation, to discover potential inhibitors of G12C and G12D mutants compared to sotorasib, which is a recent FDA-approved inhibitor of G12C. We identified compounds, such as flupentixol, amlodipine, and fluvoxamine, for the G12C mutant and paroxetine, flupentixol, and zuclopenthixol for the G12D mutant with significant inhibitory functions. All five compounds bound to the H95 cryptic groove of mutant K-Ras with high efficiency and, like sotorasib, retained a novel binding mechanism with additional hydrophobic interactions at the molecular level. Furthermore, the simulation studies suggested that the binding of flupentixol and amlodipine to G12C stabilizes switch I and switch II. In contrast, paroxetine and flupentixol to G12D showed a similar trend compared to sotorasib complexes. Thus, despite the very dynamic functionality of K-Ras switches I and II, the binding of shortlisted compounds is highly stable. Therefore, the reported study provides potential drug candidates for K-Ras inhibition that can be further developed with in vitro and in vivo evidence for targeted therapy.