An Integrated Computational Framework to Assess the Mutational Landscape of α-L-Iduronidase IDUA Gene

Influence of Cataract Causing Mutations on αA-Crystallin: A Computational Approach

The αA-crystallin protein plays a vital role in maintaining the refractive index and transparency of the eye lens. Significant clinical studies have emerged as the αA-crystallin is prone to aggregation, resulting in the formation of cataracts with varied etiologies due to mutations. This work aims to comprehend the structural and functional role of cataract-causing mutations in αA-crystallin, particularly at N-Terminal and α-Crystallin Domains, using in-silico approaches including molecular dynamics simulation. About 19 mutants of αA-crystallin along with native structure were simulated for 100 ns and the post-simulations analyses reveal pronounced dynamics of αA-crystallin due to the enhanced structure flexibility as its native compactness was lost and is witnessed mainly by the mutants R12L, R21L, R21Q, R54L, R65Q, R116C and R116H. It is observed that αA-crystallin discloses the NTD motions as the dominant one and the same was endorsed by the linear variation between Rg and the center-of-mass of αA-crystallin. Interestingly, such enhanced dynamics of αA-crystallin mutants associated with the structure flexibility is internally modulated by the dynamic exchange of secondary structure elements β-sheets and coils (R2 = 0.619) during simulation. Besides, the observed pronounced dynamics of dimer interface region (β3-L6-β4 segment) of ACD along with CTD dynamics also gains importance. Particularly, the highly dynamic mutants are also characterized by enhanced non-covalent and hydrophobic interactions which renders detrimental effects towards its stability, and favours possible protein unfolding mechanisms. Overall, this study highlights the mutation-mediated structural distortions in αA-crystallin and demands the need for further potential development of inhibitors against cataract formation.

Elucidation of escitalopram oxalate and related antidepressants as putative inhibitors of PTP4A3/PRL-3 protein in hepatocellular carcinoma: A multi-computational investigation

Hepatocellular carcinoma (HCC) persists to be one of the most devastating and deadliest malignancies globally. Recent research into the molecular signaling networks entailed in many malignancies has given some prominent insights that can be leveraged to create molecular therapeutics for combating HCC. Therefore, in the current communication, an in-silico drug repurposing approach has been employed to target the function of PTP4A3/PRL-3 protein in HCC using antidepressants: Fluoxetine hydrochloride, Citalopram, Amitriptyline, Imipramine, and Escitalopram oxalate as the desired ligands. The density function theory (DFT) and chemical absorption, distribution, metabolism, excretion, and toxicity (ADMET) parameters for the chosen ligands were evaluated to comprehend the pharmacokinetics, drug-likeness properties, and bioreactivity of the ligands. The precise interaction mechanism was explored using computational methods such as molecular docking and molecular dynamics (MD) simulation studies to assess the inhibitory effect and the stability of the interactions against the protein of interest. Escitalopram oxalate exhibited a comparatively significant docking score (-7.4 kcal/mol) compared to the control JMS-053 (-6.8 kcal/mol) against the PRL-3 protein. The 2D interaction plots exhibited an array of hydrophobic and hydrogen bond interactions. The findings of the ADMET forecast confirmed that it adheres to Lipinski's rule of five with no violations, and DFT analysis revealed a HOMO-LUMO energy gap of -0.26778 ev, demonstrating better reactivity than the control molecule. The docked complexes were subjected to MD studies (100 ns) showing stable interactions. Considering all the findings, it can be concluded that Escitalopram oxalate and related therapeutics can act as potential pharmacological candidates for targeting the activity of PTP4A3/PRL-3 in HCC.

In silico prediction of mutation sites for anthranilate synthase from Serratia marcesens to deregulate tryptophan feedback inhibition

Allosteric feedback inhibition of the committed step in amino acid biosynthetic pathways is a major concern for production of amino acids at industrial scale. Anthranilate synthase (AS) catalyzes the first reaction of tryptophan biosynthetic pathway found in microorganisms and is feedback inhibited by its own product i.e. tryptophan. Here, we identified new mutant sites in AS using computational mutagenesis approach. MD simulations (20 ns) followed by MMPBSA and per residue decomposition energy analysis identified seven amino acid residues with best binding affinity for tryptophan. All 19 mutant structures were generated for each identified amino acid residue followed by simulation to evaluate effect of mutation on protein stability. Later, molecular docking studies were employed to generate mutant-tryptophan complex and structures with binding energies (kcal/mol) much higher than wild-type AS were selected. Finally, two mutants i.e., S37W and S37H were identified on the basis of positive binding scores and loss of tryptophan binding inside pocket. Further, MD simulations run for 200 ns were performed over these mutant-tryptophan complexes followed by RMSD, RMSF, radius of gyration , solvent accessible surface area , intra-protein hydrogen bond numbers, principal component analysis, free energy landscape (FEL) and secondary structure analysis to rationale effect of mutations on stability of protein. Cross correlation analysis of mutant site amino acids (S37W) with key residues of catalytic site (G325, T326, H395 and G482) was done to evaluate the effect of mutations on catalytic site conformation. Current computational mutagenesis approach predicted two mutants S37W and S37H with proposed deregulated feedback inhibition by tryptophan and retained catalytic activity.Communicated by Ramaswamy H. Sarma.

Novel in-frame duplication variant characterization in late infantile metachromatic leukodystrophy using whole-exome sequencing and molecular dynamics simulation

Metachromatic leukodystrophy (MLD) is a neurodegenerative lysosomal storage disease caused by a deficiency in the arylsulfatase A (ARSA). ARSA deficiency leads to sulfatide accumulation, which involves progressive demyelination. The profound impact of early diagnosis on MLD treatment options necessitates the development of new or updated analysis tools and approaches. In this study, to identify the genetic etiology in a proband from a consanguineous family with MLD presentation and low ARSA activity, we employed Whole-Exome Sequencing (WES) followed by co-segregation analysis using Sanger sequencing. Also, MD simulation was utilized to study how the variant alters the structural behavior and function of the ARSA protein. GROMACS was applied and the data was analyzed by RMSD, RMSF, Rg, SASA, HB, atomic distance, PCA, and FEL. Variant interpretation was done based on the American College of Medical Genetics and Genomics (ACMG) guidelines. WES results showed a novel homozygous insertion mutation, c.109_126dup (p.Asp37_Gly42dup), in the ARSA gene. This variant is located in the first exon of ARSA, fulfilling the criteria of being categorized as likely pathogenic, according to the ACMG guidelines and it was also found to be co-segregating in the family. The MD simulation analysis revealed this mutation influenced the structure and the stabilization of ARSA and led to the protein function impairment. Here, we report a useful application of WES and MD to identify the causes of a neurometabolic disorder.

An Allosteric Inhibitory Potential of Triterpenes from Combretum racemosum on the Structural and Functional Dynamics of Plasmodium falciparum Lactate Dehydrogenase Binding Landscape

Multidrug resistance is a significant drawback in malaria treatment, and mutations in the active sites of the many critical antimalarial drug targets have remained challenging. Therefore, this has necessitated the global search for new drugs with new mechanisms of action. Plasmodium falciparum lactate dehydrogenase (pfLHD), a glycolytic enzyme, has emerged as a potential target for developing new drugs due to the parasite reliance on glycolysis for energy. Strong substrate-binding is required in pfLDH enzymatic catalysis; however, there is a lack of information on small molecules' inhibitory mechanism bound to the substrate-binding pocket. Therefore, this study investigated a potential allosteric inhibition of pfLDH by targeting the substrate-binding site. The structural and functional behaviour of madecassic acid (MA), the most promising among the six triterpenes bound to pfLDH, were unravelled using molecular dynamic simulations at 300 ns to gain insights into its mechanism of binding and inhibition and chloroquine as a standard drug. The docking studies identified that the substrate site has the preferred position for the compounds even in the absence of a co-factor. The bound ligands showed comparably higher binding affinity at the substrate site than at the co-factor site. Mechanistically, a characteristic loop implicated in the enzyme catalytic activity was identified at the substrate site. This loop accommodates key interacting residues (LYS174, MET175, LEU177 and LYS179) pivotal in the MA binding and inhibitory action. The MA-bound pfLHD average RMSD (1.60 Å) relative to chloroquine-bound pfLHD RMSD (2.00 Å) showed higher stability for the substrate pocket, explaining the higher binding affinity (-33.40 kcal/mol) observed in the energy calculations, indicating that MA exhibited profound inhibitory activity. The significant pfLDH loop conformational changes and the allostery substrate-binding landscape suggested inhibiting the enzyme function, which provides an avenue for designing antimalarial compounds in the future studies of pfLDH protein as a target.

c.1898C>G/p.Ser633Trp Mutation in Alpha-L-Iduronidase: Clinical and Structural Implications

Mucopolysaccharidosis type I is a rare autosomal recessive genetic disease caused by deficient activity of α-L-iduronidase. As a consequence of low or absent activity of this enzyme, glycosaminoglycans accumulate in the lysosomal compartments of multiple cell types throughout the body. Mucopolysaccharidosis type I has been classified into 3 clinical subtypes, ranging from a severe Hurler form to the more attenuated Hurler-Scheie and Scheie phenotypes. Over 200 gene variants causing the various forms of mucopolysaccharidosis type I have been reported. DNA isolated from dried blood spot was used to sequencing of all exons of the IDUA gene from a patient with a clinical phenotype of severe mucopolysaccharidosis type I syndrome. Enzyme activity of α-L-iduronidase was quantified by fluorimetric assay. Additionally, a molecular dynamics simulation approach was used to determine the effect of the Ser633Trp mutation on the structure and dynamics of the α-L-iduronidase. The DNA sequencing analysis and enzymatic activity shows a c.1898C>G mutation associated a patient with a homozygous state and α-L-iduronidase activity of 0.24 μmol/L/h, respectively. The molecular dynamics simulation analysis shows that the p.Ser633Trp mutation on the α-L-iduronidase affect significant the temporal and spatial properties of the different structural loops, the N-glycan attached to Asn372 and amino acid residues around the catalytic site of this enzyme. Low enzymatic activity observed for p.Ser633Trp variant of the α-L-iduronidase seems to lead to severe mucopolysaccharidosis type I phenotype, possibly associated with a perturbation of the structural dynamics in regions of the enzyme close to the active site.

Weak spots inhibition in the Mycobacterium tuberculosis antigen 85C target for antitubercular drug design through selective irreversible covalent inhibitor-SER124

Mycobacterium tuberculosis (Mtb) encoded secreted antigen 85 enzymes (Ag85A/Ag85B/Ag85C) play that critical roles in the virulence, survival and drug-resistant TB of the pathogen. Ag85 proteins are potential antitubercular drug targets because they are essential in the catalytic synthesis of trehalose moieties and mycolic acid attachment to the Mtb cell wall. Recently, experimental protocols led to the discovery of a selective covalent Ag85 inhibitor, β-isomer monocyclic enolphosphorus Cycliphostin (CyC_8β) compound, which targets the Ag85 serine 124 to exhibit a promising therapeutic activity. For the first time, our study unravelled the structural features among Mtb Ag85C homologs and motions and dynamics of Ag85C when the CyC_8β bound covalently and in open model conformations to the protein using bioinformatics tools and integrated Molecular dynamics simulations. Comparative Ag85C sequence analysis revealed conserved regions; 70% active site, 90% Adeniyi loop L1 and 50% loop L2, which acts as a switch between open and closed conformations. The average C-α atoms RMSD (2.05 Å) and RMSF (0.9 Å) revealed instability and high induced flexibility in the CyC_8β covalent-bound compared to the apo and open model systems, which displayed more stability and lower fluctuations. DSSP showed structural transitions of α-helices to bend and loops to 3₁₀-helices in the bound systems. SASA of CyC_8β covalent bound showed active site hydrophobic residues exposure to huge solvent. Therefore, these findings present the potential opportunity hotspots in Ag85C protein that would aid the structure-based design of novel chemical entities capable of resulting in potent antitubercular drugs.Communicated by Ramaswamy H. Sarma.

Insight into the structural and functional analysis of the impact of missense mutation on cytochrome P450 oxidoreductase

Cytochrome P450 oxidoreductase (POR) is a steroidogenic and drug-metabolizing enzyme. It helps in the NADPH dependent transfer of electrons to cytochrome P450 (CYP) enzymes for their biological activity. In this study, we employed integrative computational approaches to decipher the impact of proline to leucine missense mutation at position 384 (P384L) in the connecting/hinge domain region which is essential for the catalytic activity of POR. Analysis of protein stability using DUET, MUpro, CUPSAT, I-Mutant2.0, iStable and SAAFEC servers predicted that mutation might alter the structural stability of POR. The significant conformational changes induced by the mutation to the POR structure were analyzed by long-range molecular dynamics simulation. The results revealed that missense mutation decreased the conformational stability of POR as compared to wild type (WT). The PCA based FEL analysis described the mutant-specific conformational alterations and dominant motions essential for the biological activity of POR. The LIGPLOT interaction analysis showed the different binding architecture of FMN, FAD, and NADPH as a result of mutation. The increased number of hydrogen bonds in the FEL conformation of WT proved the strong binding of cofactors in the binding pocket as compared to the mutant. The porcupine plot analysis associated with cross-correlation analysis depicted the high-intensity flexible motion exhibited by functionally important FAD and NADPH binding domain regions. The computational findings unravel the impact of mutation at the structural level, which could be helpful in understanding the molecular mechanism of drug metabolism.

Identification of potential inhibitors against pathogenic missense mutations of PMM2 using a structure-based virtual screening approach

The autosomal recessive phosphomannomutase 2-congenital disorder of glycosylation (PMM2-CDG) is characterized by defective functioning of the PMM2 enzyme, which is necessary for the conversion of mannose-6-phosphate into mannose-1-phosphate. Here, a computational pipeline was drawn to identify the most significant mutations, and further, we used a virtual screening approach to identify a new lead compound to treat the identified significant mutations. We searched for missense mutation data related to PMM2-CDG in HGMD®, UniProt, and ClinVar. Our search yielded a total of 103 mutations, of which 91 are missense mutations. The D65Y, I132N, I132T, and F183S mutations were classified as deleterious, destabilizing, and altering the biophysical properties using the PredictSNP, iStable, and Align GVGD in silico prediction tools. Additionally, we applied a multistep protocol to screen for an alternative lead compound to the existing CID2876053 (1-(3-chlorophenyl)-3,3-bis(pyridine-2-yl)urea) with affinity to these identified significant mutants. Two compounds, CHEMBL1491007 (6-chloro-4-phenyl-3-(4-pyridin-2-ylpiperazin-1-yl)-1H-quinolin-2-one) and CHEMBL3653029 (5-chloro-4-[6-[(3-fluorophenyl)methylamino]pyridin-2-yl]-N-(piperidin-4-ylmethyl)pyridin-2-amine), exhibited the highest binding affinity with the selected mutants and were chosen for further analysis. Through molecular docking, molecular dynamics simulation, and MMPBSA analysis, we found that the known compound, i.e. CID2876053, has stronger interaction with the D65Y mutant. The newly identified lead compound CHEMBL1491007 showed stronger interaction with the I132N and I132T mutants, whereas the most deleterious mutant, F183S, showed stronger interaction with CHEMBL3653029. This study is expected to aid in the field of precision medicine, and further to in vivo and in vitro analysis of these lead compounds might shed light on the treatment of PMM2-CDG. Communicated by Ramaswamy H. Sarma.

Thompson loop: opportunities for antitubercular drug design by targeting the weak spot in demethylmenaquinone methyltransferase protein

Graphical superimposed snapshots of the Thompson novel loop (yellow) of menG protein: apo (A) and bound (B) systems. The loop switches between open and closed conformations; critical for therapeutic activity.

Computational modelling approaches as a potential platform to understand the molecular genetics association between Parkinson's and Gaucher diseases

Gaucher's disease (GD) is a genetic disorder in which glucocerebroside accumulates in cells and specific organs. It is broadly classified into type I, type II and type III. Patients with GD are at high risk of Parkinson's disease (PD), and the clinical and pathological presentation of GD patients with PD is almost identical to idiopathic PD. Several experimental models like cell culture, animal models, and transgenic mice models were used to understand the molecular mechanism behind GD and PD association; however, such mechanism remains unclear. In this context, based on literature reports, we identified the most common mutations K198T, E326K, T369M, N370S, V394L, D409H, L444P, and R496H, in the Glucosylceramidase (GBA) protein that are known to cause GD1, and represent a risk of developing PD. However, to date, no computational analyses have designed to elucidate the potential functional role of GD mutations with increased risk of PD. The present computational pipeline allows us to understand the structural and functional significance of these GBA mutations with PD. Based on the published data, the most common and severe mutations were E326K, N370S, and L444P, which further selected for our computational analysis. PredictSNP and iStable servers predicted L444P mutant to be the most deleterious and responsible for the protein destabilization, followed by the N370S mutation. Further, we used the structural analysis and molecular dynamics approach to compare the most frequent deleterious mutations (N370S and L444P) with the mild mutation E326K. The structural analysis demonstrated that the location of E326K and N370S in the alpha helix region of the protein whereas the mutant L444P was in the starting region of the beta sheet, which might explain the predicted pathogenicity level and destabilization effect of the L444P mutant. Finally, Molecular Dynamics (MD) at 50 ns showed the highest deviation and fluctuation pattern in the L444P mutant compared to the two mutants E326K and N370S and the native protein. This was consistent with more loss of intramolecular hydrogen bonds and less compaction of the radius of gyration in the L444P mutant. The proposed study is anticipated to serve as a potential platform to understand the mechanism of the association between GD and PD, and might facilitate the process of drug discovery against both GD and PD.