A Case of Fabry Disease With Lacrimal Gland Involvement.

Fabry disease is an X-linked lysosomal storage disease resulting from an error in the glycosphingolipid metabolic pathway, which leads to accumulation of globotriaosylceramide in lysosomes of the skin, kidneys, heart, brain, and other organs. There are no existing reports of histologically proven lacrimal gland involvement in Fabry disease. The authors report the case of a 26-year-old male with Fabry disease who presented with bilateral upper eyelid dermatochalasis, steatoblepharon, and prolapsed lacrimal glands. The patient underwent surgical repair of the upper eyelids and biopsy of the lacrimal glands. The pathologic assessment demonstrated lamellated intracytoplasmic inclusions characteristic of Fabry disease. The prevalence of globotriaosylceramide lacrimal gland deposition in Fabry disease and the effect on lacrimal gland morphology and function have yet to be determined.

Inhibitory mechanism of an antifungal drug, caspofungin against amyloid ß peptide aggregation: Repurposing via neuroinformatics and an experimental approach.

The multifactorial neurological condition called Alzheimer's disease (AD) primarily affects elderly individuals. Despite the calamitous consequences of AD, curative strategies for a regimen to apply remain inadequate as several factors contribute to AD etiology. Drug repurposing is an advance strategy prior to drug discovery as various effective drugs perform through alteration of multiple targets, and the present "poly-pharmacology" can be a curative approach to complex disorders. AD's multifactorial behavior actively encourages the hypothesis for a drug design approach focused on drug repurposing. In this study, we discovered that an antifungal drug, Caspofungin (CAS) is a potent Aß aggregation inhibitor that displays significantly reduced toxicity associated with AD. Drug reprofiling and REMD simulations demonstrated that CAS interacts with the ß-sheet section, known as Aß amyloid fibrils hotspot. CAS leads to destabilization of ß-sheet and, conclusively, in its devaluation. Later, in vitro experiments were acquired in which the fibrillar volume was reduced for CAS-treated Aß peptide. For the first time ever, this study has determined an antifungal agent as the Aß amyloid aggregation's potent inhibitor. Several efficient sequence-reliant potent inhibitors can be developed in future against the amyloid aggregation for different amyloid peptide by the processing and conformational optimization of CAS.

Predicting phosphorylation sites using machine learning by integrating the sequence, structure, and functional information of proteins.

BACKGROUND: Post-translational modification (PTM) is a biological process that alters proteins and is therefore involved in the regulation of various cellular activities and pathogenesis. Protein phosphorylation is an essential process and one of the most-studied PTMs: it occurs when a phosphate group is added to serine (Ser, S), threonine (Thr, T), or tyrosine (Tyr, Y) residue. Dysregulation of protein phosphorylation can lead to various diseases-most commonly neurological disorders, Alzheimer's disease, and Parkinson's disease-thus necessitating the prediction of S/T/Y residues that can be phosphorylated in an uncharacterized amino acid sequence. Despite a surplus of sequencing data, current experimental methods of PTM prediction are time-consuming, costly, and error-prone, so a number of computational methods have been proposed to replace them. However, phosphorylation prediction remains limited, owing to substrate specificity, performance, and the diversity of its features. METHODS: In the present study we propose machine-learning-based predictors that use the physicochemical, sequence, structural, and functional information of proteins to classify S/T/Y phosphorylation sites. Rigorous feature selection, the minimum redundancy/maximum relevance approach, and the symmetrical uncertainty method were employed to extract the most informative features to train the models. RESULTS: The RF and SVM models generated using diverse feature types in the present study were highly accurate as is evident from good values for different statistical measures. Moreover, independent test sets and benchmark validations indicated that the proposed method clearly outperformed the existing methods, demonstrating its ability to accurately predict protein phosphorylation. CONCLUSIONS: The results obtained in the present work indicate that the proposed computational methodology can be effectively used for predicting putative phosphorylation sites further facilitating discovery of various biological processes mechanisms.

Myocilin-associated Glaucoma: A Historical Perspective and Recent Research Progress.

Glaucoma a debilitating disease, is globally the second most common kind of permanent blindness. Primary open-angle glaucoma (POAG) is its most prevalent form and is often linked with alterations in the myocilin gene (MYOC). MYOC encodes the myocilin protein, which is expressed throughout the body, but primarily in trabecular meshwork (TM) tissue in the eyes. TM is principally involved in regulating intraocular pressure (IOP), and elevated IOP is the main risk factor associated with glaucoma. The myocilin protein's function remains unknown; however, mutations compromise its folding and processing inside TM cells, contributing to the glaucoma phenotype. While glaucoma is a complex disease with various molecules and factors as contributing causes, the role played by myocilin has been the most widely studied. The current review describes the present understanding of myocilin and its association with glaucoma and aims to shift the focus toward developing targeted therapies for treating glaucoma patients with variations in MYOC.

1,3,4-oxadiazole conjugates of capsaicin as potent NorA efflux pump inhibitors of Staphylococcus aureus.

NorA efflux pump pertaining to the major facilitator superfamily (MFS) is known to play a key role in antibiotic and biocide resistance in Staphylococcus aureus (S. aureus). It accounts for the extrusion of antibiotics like fluoroquinolones (e.g. ciprofloxacin). Several compounds including synthetic and natural products have been identified as potential NorA efflux pump inhibitors (EPIs) and found to restore the antibacterial activity of antibiotics. However, none of the reported EPIs have reached to clinical approval probably due to their high toxicity profiles. Considering the NorA efflux pump inhibitory potential of capsaicin, a series of capsaicin-based 1,3,4 oxadiazole conjugates were prepared and evaluated for ciprofloxacin activity potentiating effect. Among the new capsaicinoids tested, 17i displayed a minimum effective concentration (MEC) of 12.5 µg/mL against NorA overexpressing S. aureus strain (SA1199B), whereas capsaicin showed MEC of 50 µg/mL. The kill kinetics curve for the combination showed that ciprofloxacin at a sub-inhibitory concentration (0.25 × MIC) was equipotent in effect, to its MIC. 17i has significantly decreased the ethidium bromide efflux confirming NorA inhibition as the mode of action. Mutation prevention concentration of the ciprofloxacin was reduced in combination with 17i.In silico studies revealed the binding efficiency and binding affinity of 17i with NorA. This compound may serve as a template for the further drug discovery.

Structure-based drug repurposing to inhibit the DNA gyrase of Mycobacterium tuberculosis.

Drug repurposing is an alternative avenue for identifying new drugs to treat tuberculosis (TB). Despite the broad-range of anti-tubercular drugs, the emergence of multi-drug-resistant and extensively drug-resistant strains of Mycobacterium tuberculosis (Mtb) H37Rv, as well as the significant death toll globally, necessitates the development of new and effective drugs to treat TB. In this study, we have employed a drug repurposing approach to address this drug resistance problem by screening the drugbank database to identify novel inhibitors of the Mtb target enzyme, DNA gyrase. The compounds were screened against the ATPase domain of the gyrase B subunit (MtbGyrB47), and the docking results showed that echinacoside, doxorubicin, epirubicin, and idarubicin possess high binding affinities against MtbGyrB47. Comprehensive assessment using fluorescence spectroscopy, surface plasmon resonance spectroscopy (SPR), and circular dichroism (CD) titration studies revealed echinacoside as a potent binder of MtbGyrB47. Furthermore, ATPase, and DNA supercoiling assays exhibited an IC50 values of 2.1-4.7 µM for echinacoside, doxorubicin, epirubicin, and idarubicin. Among these compounds, the least MIC90 of 6.3 and 12 µM were observed for epirubicin and echinacoside, respectively, against Mtb. Our findings indicate that echinacoside and epirubicin targets mycobacterial DNA gyrase, inhibit its catalytic cycle, and retard mycobacterium growth. Further, these compounds exhibit potential scaffolds for optimizing novel anti-mycobacterial agents that can act on drug-resistant strains.

Specific keratinase derived designer peptides potently inhibit Aß aggregation resulting in reduced neuronal toxicity and apoptosis.

Compelling evidence implicates self-assembly of amyloid-ß (Aß1-42) peptides into soluble oligomers and fibrils as a major underlying event in Alzheimer's disease (AD) pathogenesis. Herein, we employed amyloid-degrading keratinase (kerA) enzyme as a key Aß1-42-binding scaffold to identify five keratinase-guided peptides (KgPs) capable of interacting with and altering amyloidogenic conversion of Aß1-42 The KgPs showed micromolar affinities with Aß1-42 and abolished its sigmoidal amyloidogenic transition, resulting in abrogation of fibrillogenesis. Comprehensive assessment using dynamic light scattering (DLS), atomic force microscopy (AFM) and Fourier-transform infrared (FTIR) spectroscopy showed that KgPs induced the formation of off-pathway oligomers comparatively larger than the native Aß1-42 oligomers but with a significantly reduced cross-ß signature. These off-pathway oligomers exhibited low immunoreactivity against oligomer-specific (A11) and fibril-specific (OC) antibodies and rescued neuronal cells from Aß1-42 oligomer toxicity as well as neuronal apoptosis. Structural analysis using molecular docking and molecular dynamics (MD) simulations showed two preferred KgP binding sites (Lys16-Phe20 and Leu28-Val39) on the NMR ensembles of monomeric and fibrillar Aß1-42, indicating an interruption of crucial hydrophobic and aromatic interactions. Overall, our results demonstrate a new approach for designing potential anti-amyloid molecules that could pave way for developing effective therapeutics against AD and other amyloid diseases.

Pyrazinamide drug resistance in RpsA mutant (∆438A) of Mycobacterium tuberculosis: Dynamics of essential motions and free-energy landscape analysis.

Pyrazinamide is an essential first-line antitubercular drug which plays pivotal role in tuberculosis treatment. It is a prodrug that requires amide hydrolysis by mycobacterial pyrazinamidase enzyme for conversion into pyrazinoic acid (POA). POA is known to target ribosomal protein S1 (RpsA), aspartate decarboxylase (PanD), and some other mycobacterial proteins. Spontaneous chromosomal mutations in RpsA have been reported for phenotypic resistance against pyrazinamide. We have constructed and validated 3D models of the native and Δ438A mutant form of RpsA protein. RpsA protein variants were then docked to POA and long range molecular dynamics simulations were carried out. Per residue binding free-energy calculations, free-energy landscape analysis, and essential dynamics analysis were performed to outline the mechanism underlying the high-level PZA resistance conferred by the most frequently occurring deletion mutant of RpsA. Our study revealed the conformational modulation of POA binding site due to the disruptive collective modes of motions and increased conformational flexibility in the mutant than the native form. Residue wise MMPBSA decomposition and protein-drug interaction pattern revealed the difference of energetically favorable binding site in the wild-type (WT) protein in comparison with the mutant. Analysis of size and shape of minimal energy landscape area delineated higher stability of the WT complex than the mutant form. Our study provides mechanistic insights into pyrazinamide resistance in Δ438A RpsA mutant, and the results arising out of this study will pave way for design of novel and effective inhibitors targeting the resistant strains of Mycobacterium tuberculosis.

Structural basis for drug resistance mechanisms against anaplastic lymphoma kinase.

Drug resistance to anaplastic lymphoma kinase (ALK) inhibitors (crizotinib and ceritinib) is caused by mutation in the region encoding kinase domain of ALK. Compounds with potential ability to inhibit all strains of ALK are a solution to tackle the problem of drug resistance. In this study, we delineated positions of residues possessing the ability to make ALK drug resistant upon mutation by assessing them using five parameters (conservation index, binding-site root-mean-square deviation, protein structure stability, change in ATP, and drug-binding affinity). Four residual positions (Leu 1122, Thr 1151, Phe 1245, and Gly 1269) were ascertained. This study will be beneficial for designing drugs with better proficiency against ALK and the issues of drug resistance. This study can be taken as a pipeline for investigating drug-resistant mutations in other diseases as well.

Computational models for the prediction of adverse cardiovascular drug reactions.

BACKGROUND: Predicting adverse drug reactions (ADRs) has become very important owing to the huge global health burden and failure of drugs. This indicates a need for prior prediction of probable ADRs in preclinical stages which can improve drug failures and reduce the time and cost of development thus providing efficient and safer therapeutic options for patients. Though several approaches have been put forward for in silico ADR prediction, there is still room for improvement. METHODS: In the present work, we have used machine learning based approach for cardiovascular (CV) ADRs prediction by integrating different features of drugs, biological (drug transporters, targets and enzymes), chemical (substructure fingerprints) and phenotypic (therapeutic indications and other identified ADRs), and their two and three level combinations. To recognize quality and important features, we used minimum redundancy maximum relevance approach while synthetic minority over-sampling technique balancing method was used to introduce a balance in the training sets. RESULTS: This is a rigorous and comprehensive study which involved the generation of a total of 504 computational models for 36 CV ADRs using two state-of-the-art machine-learning algorithms: random forest and sequential minimization optimization. All the models had an accuracy of around 90% and the biological and chemical features models were more informative as compared to the models generated using chemical features. CONCLUSIONS: The results obtained demonstrated that the predictive models generated in the present study were highly accurate, and the phenotypic information of the drugs played the most important role in drug ADRs prediction. Furthermore, the results also showed that using the proposed method, different drugs properties can be combined to build computational predictive models which can effectively predict potential ADRs during early stages of drug development.

Structural basis for isoniazid resistance in KatG double mutants of Mycobacterium tuberculosis.

The failure of drugs for effective treatment against infectious diseases can be attributed to resistant forms of causative agents. The evasive nature of Mycobacterium tuberculosis is partly associated to its physical features, such as having a thick cell wall and incorporation of beneficial mutations leading to drug resistance. The pro drug Isoniazid (INH) interacts with an enzyme catalase peroxidase to get converted into its active form and upon activation stops the cell wall synthesis thus killing the Mycobacterium. The most common mutation i.e. S315T leads to high degree of drug resistance by virtue of its position in the active site. Here, we have characterized the prominent attributes of two double mutant isolates S315 T/D194G and S315T/M624V which are multi drug resistant and extremely drug resistant, respectively. Protein models were generated using the crystal structure which were then subjected to energy minimization and long term molecular dynamics simulations. Further, computational analysis showed decreasing ability of INH binding to the mutants in order of: Native > S315T/D194G > S315T/M624V. Also, a trend was observed that as the docking score and binding area decreased, there was a significant increase in the distortion of the 3D geometry of the mutants as observed by PCA analysis.

Molecular dynamics simulations revealed structural differences among WRKY domain-DNA interaction in barley (Hordeum vulgare).

BACKGROUND: The WRKY transcription factors are a class of DNA-binding proteins involved in diverse plant processes play critical roles in response to abiotic and biotic stresses. Genome-wide divergence analysis of WRKY gene family in Hordeum vulgare provided a framework for molecular evolution and functional roles. So far, the crystal structure of WRKY from barley has not been resolved; moreover, knowledge of the three-dimensional structure of WRKY domain is pre-requisites for exploring the protein-DNA recognition mechanisms. Homology modelling based approach was used to generate structures for WRKY DNA binding domain (DBD) and its variants using AtWRKY1 as a template. Finally, the stability and conformational changes of the generated model in unbound and bound form was examined through atomistic molecular dynamics (MD) simulations for 100 ns time period. RESULTS: In this study, we investigated the comparative binding pattern of WRKY domain and its variants with W-box cis-regulatory element using molecular docking and dynamics (MD) simulations assays. The atomic insight into WRKY domain exhibited significant variation in the intermolecular hydrogen bonding pattern, leading to the structural anomalies in the variant type and differences in the DNA-binding specificities. Based on the MD analysis, residual contribution and interaction contour, wild-type WRKY (HvWRKY46) were found to interact with DNA through highly conserved heptapeptide in the pre- and post-MD simulated complexes, whereas heptapeptide interaction with DNA was missing in variants (I and II) in post-MD complexes. Consequently, through principal component analysis, wild-type WRKY was also found to be more stable by obscuring a reduced conformational space than the variant I (HvWRKY34). Lastly, high binding free energy for wild-type and variant II allowed us to conclude that wild-type WRKY-DNA complex was more stable relative to variants I. CONCLUSIONS: The results of our study revealed complete dynamic and structural information about WRKY domain-DNA interactions. However, no structure base information reported to date for WRKY variants and their mechanism of interaction with DNA. Our findings highlighted the importance of selecting a sequence to generate newer transgenic plants that would be increasingly tolerance to stress conditions.

Dynamics of Dof domain-DNA interaction in wheat: Insights from atomistic simulations and free energy landscape.

DNA-binding one zinc finger protein (Dof) is a plant-specific transcription factor involved in numerous biological processes. In the current study, the plausible mechanism underlying Dof domain-DNA interaction in wheat was investigated using extensive molecular dynamics (MD) simulations analysis. We depicted that one key residue Lys29, possessing the ability to disturb the interaction between Dof domain-DNA upon substitution to Arg29. Frequent conformational changes were observed in Lys29Arg (K29R)-DNA complex during the entire MD simulation period, which significantly altered the interactions, thereby indicating the importance of Lys29 in complex stability. Principal component analysis and free energy landscape results also suggested strong affinity between wild-type Dof domain and DNA due to restricted atomic movement. Our study not only substantiates the structural and mechanistic insights of Dof transcription factor but also provides new avenues toward employment of these key amino acid residues in genetic engineering for development of abiotic stress tolerance crop plant.

Simulation Based Investigation of Deleterious nsSNPs in ATXN2 Gene and Its Structural Consequence Toward Spinocerebellar Ataxia.

Spinocerebellar degeneration, termed as ataxia is a neurological disorder of central nervous system, characterized by limb in-coordination and a progressive gait. The patient also demonstrates specific symptoms of muscle weakness, slurring of speech, and decreased vibration senses. Expansion of polyglutamine trinucleotide (CAG) within ATXN2 gene with 35 or more repeats, results in spinocerebellar ataxia type-2. Protein ataxin-2 coded by ATXN2 gene has been reported to have a crucial role in translation of the genetic information through sequestering the histone acetyl transferases (HAT) resulting in a state of hypo-acetylation. In the present study, we have evaluated the outcome for 122 non synonymous single nucleotide polymorphisms (nsSNPs) reported within ATXN2 gene through computational tools such as SIFT, PolyPhen 2.0, PANTHER, I-mutant 2.0, Phd-SNP, Pmut, MutPred. The apo and mutant (L305V and Q339L) form of structures for the ataxin-2 protein were modeled for gaining insights toward 3D spatial arrangement. Further, molecular dynamics simulations and structural analysis were performed to observe the brunt of disease associated nsSNPs toward the strength and secondary properties of ataxin-2 protein structure. Our results showed that, L305V is a highly deleterious and disease causing point substitution. Analysis based on RMSD, RMSF, Rg, SASA, number of hydrogen bonds (NH bonds), covariance matrix trace, projection analysis for eigen vector demonstrated a significant instability and conformation along with rise in mutant flexibility values in comparison to the apo form of ataxin-2 protein. The study provides a blue print of computational methodologies to examine the ataxin-blend SNPs. J. Cell. Biochem. 119: 499-510, 2018. © 2017 Wiley Periodicals, Inc.

Role of pncA gene mutations W68R and W68G in pyrazinamide resistance.

Mycobacterium tuberculosis (Mtb) resistance toward anti-tuberculosis drugs is a widespread problem. Pyrazinamide (PZA) is a first line antitubercular drug that kills semi-dormant bacilli when converted into its activated form, that is, pyrazinoic acid (POA) by Pyrazinamidase (PZase) enzyme coded by pncA gene. In this study, we conducted several analyses on native and mutant structures (W68R, W68G) of PZase before and after docking with the PZA drug to explore the molecular mechanism behind PZA resistance caused due to pncA mutations. Structural changes caused by mutations were studied with respect to their effects on functionality of protein. Docking was performed to analyze the protein-drug binding and comparative analysis was done to observe how the mutations affect drug binding affinity and binding site on protein. Native PZase protein was observed to have the maximum binding affinity in terms of docking score as well as shape complementarity in comparison to the mutant forms. Molecular dynamics simulation analyses showed that mutation in the 68th residue of protein results in a structural change at its active site which further affects the biological function of protein, that is, conversion of PZA to POA. Mutations in the protein thereby led to PZA resistance in the bacterium due to the inefficient binding.

Computational Screening and Exploration of Disease-Associated Genes in Alzheimer's Disease.

Alzheimer's is a neurodegenerative disease affecting large populations worldwide characterized mainly by progressive loss of memory along with various other symptoms. The foremost cause of the disease is still unclear, however various mechanisms have been proposed to cause the disease that include amyloid hypothesis, tau hypothesis, and cholinergic hypothesis in addition to genetic factors. Various genes have been known to be involved which are APOE, PSEN1, PSEN2, and APP among others. In the present study, we have used computational methods to examine the pathogenic effects of non-synonymous single nucleotide polymorphisms (SNPs) associated with ABCA7, CR1, MS4A6A, CD2AP, PSEN1, PSEN2, and APP genes. The SNPs were obtained from dbSNP database followed by identification of deleterious SNPs and prediction of their functional impact. Prediction of disease-associated mutations was performed and the impact of the mutations on the stability of the protein was carried out. To study the structural significance of the computationally prioritized mutations on the proteins, molecular dynamics simulation studies were carried out. On analysis, the SNPs with IDs rs76282929 ABCA7; CR1 rs55962594; MS4A6A rs601172; CD2AP rs61747098; PSEN1 rs63750231, rs63750265, rs63750526, rs63750577, rs63750687, rs63750815, rs63750900, rs63751037, rs63751163, rs63751399; PSEN2 rs63749851; and APP rs63749964, rs63750066, rs63750734, and rs63751039 were predicted to be deleterious and disease-associated having significant structural impact on the proteins. The current study proposes a precise computational methodology for the identification of disease-associated SNPs. J. Cell. Biochem. 118: 1471-1479, 2017. © 2016 Wiley Periodicals, Inc.

Mechanistic Principles Behind Molecular Mechanism of Rifampicin Resistance in Mutant RNA Polymerase Beta Subunit of Mycobacterium tuberculosis.

Evolution of drug-resistant Mycobacterium strains threatens the TB treatment and control programs globally. Rifampicin (RIF) is an important first line antitubercular drug. Resistance to Rifampicin is caused mainly by mutations in its target RNA polymerase beta subunit protein (RpoB). RpoB contains a Rifampicin resistance determining region (RRDR) and has several potent sites for mutations. In this study, we have investigated mutations of a single site (H451) to eight different amino acids, involved in RIF resistance. Long-term molecular dynamics simulations were performed on wild type (WT) and mutant protein structures and various structural analysis were carried out to elucidate the dynamic behavior of WT and mutant forms. Essential dynamics uncovered the difference in conformational flexibility and collective modes of motions between WT and mutants. MMPBSA calculations and interaction pattern analysis revealed the binding site relocation in some mutants. This study presents an exhaustive analysis of RIF binding to the WT and mutant RpoB and clearly highlights structural mechanism for differences in stable binding of Rifampicin with WT than the mutant targets. J. Cell. Biochem. 118: 4594-4606, 2017. © 2017 Wiley Periodicals, Inc.

Double Mutants in DNA Gyrase Lead to Ofloxacin Resistance in Mycobacterium tuberculosis.

Fluoroquinolones are among the most important classes of highly effective antibacterial drugs, exhibiting wide range of activity to cure infectious diseases. Ofloxacin is second generation fluoroquinolone approved by FDA for the treatment of tuberculosis by selectively inhibiting DNA gyrase. However, the emergence of drug resistance owing to mutations in DNA gyrase poses intimidating challenge for the effective therapy of this drug. The double mutants GyrAA90V GyrBD500N and GyrAA90V GyrBT539N are reported to be implicated in conferring higher levels of OFX resistance. The present study was designed to unravel the molecular principles behind development of resistance by the bug against fluoroquinolones. Our results highlighted that polar interactions play critical role in the development of drug resistance and highlight the significant correlation between the free energy calculations predicted by MM-PBSA and stability of the ligand-bound complexes. Modifications at the OFX binding pocket due to amino acid substitution leads to fewer hydrogen bonds in mutants DNA gyrase-OFX complex, which determined the low susceptibility of the ligand in inhibiting the mutant protein. This study provides a structural rationale to the mutation-based resistance to ofloxacin and will pave way for development potent fluoroquinolone-based resistant-defiant drugs. J. Cell. Biochem. 118: 2950-2957, 2017. © 2017 Wiley Periodicals, Inc.

Natural polyphenolic inhibitors against the antiapoptotic BCL-2.

The apoptotic mechanism is regulated by the BCL-2 family of proteins, such as BCL-2 or Bcl-xL, which block apoptosis while Bad, Bak, Bax, Bid, Bim or Hrk induce apoptosis. The overexpression of BCL-2 was found to be related to the progression of cancer and also providing resistance towards chemotherapeutic treatments. In the present study, we found that all polyphenols (apigenin, fisetin, galangin and luteolin) bind to the hydrophobic groove of BCL-2 and the interaction is stable throughout MD simulation run. Luteolin was found to bind with highest negative binding energy and thus, claimed highest potency towards BCL-2 inhibition followed by fisetin. The hydrophobic interactions were found to be critical for stable complex formation as revealed by the vdW energy and ligplot analysis. Finally, on the basis of data obtained during the study, it can be concluded that these polyphenols have the potential to be used as lead molecules for BCL-2 inhibition.

Molecular dynamics-based identification of novel natural mortalin-p53 abrogators as anticancer agents.

INTRODUCTION: Cancer is one of the leading causes of mortality worldwide that requires attention in terms of extensive study and research. Eradication of mortalin-p53 interaction that leads to the inhibition of transcriptional activation or blocking of p53 from functioning as a suppressor and induction of nuclear translocation of p53 can prove to be one of the useful approaches for cancer management. RESULTS: In this study, we used structure-based approach to target the p53-binding domain of mortalin in order to prevent mortalin-p53 complex formation. We screened compounds from ZINC database against the modeled mortalin protein using Glide virtual screening. The top two compounds, DTOM (ZINC 28639308) and TTOM (ZINC 38143676) with Glide score of -12.27 and -12.16, respectively, were identified with the potential to abrogate mortalin-p53 interaction. Finally, molecular dynamics simulations were used to analyze the dynamic stability of the ligand-bound complex and it was observed that residues Tyr196, Asn198, Val264 and Thr267 were involved in intermolecular interactions in both the simulated ligand-bound complexes, and thus, these residues may have a paramount role in stabilizing the binding of the ligands with the protein. CONCLUSION: These detailed insights can further facilitate the development of potent inhibitors against mortalin-p53 complex.