Repurposing DrugBank compounds as potential Plasmodium falciparum class 1a aminoacyl tRNA synthetase multi-stage pan-inhibitors with a specific focus on mitomycin.

Plasmodium falciparum aminoacyl tRNA synthetases (PfaaRSs) are potent antimalarial targets essential for proteome fidelity and overall parasite survival in every stage of the parasite's life cycle. So far, some of these proteins have been singly targeted yielding inhibitor compounds that have been limited by incidences of resistance which can be overcome via pan-inhibition strategies. Hence, herein, for the first time, we report the identification and in vitro antiplasmodial validation of Mitomycin (MMC) as a probable pan-inhibitor of class 1a (arginyl(A)-, cysteinyl(C), isoleucyl(I)-, leucyl(L), methionyl(M), and valyl(V)-) PfaaRSs which hypothetically may underlie its previously reported activity on the ribosomal RNA to inhibit protein translation and biosynthesis. We combined multiple in silico structure-based discovery strategies that first helped identify functional and druggable sites that were preferentially targeted by the compound in each of the plasmodial proteins: Ins1-Ins2 domain in Pf-ARS; anticodon binding domain in Pf-CRS; CP1-editing domain in Pf-IRS and Pf-MRS; C-terminal domain in Pf-LRS; and CP-core region in Pf-VRS. Molecular dynamics studies further revealed that MMC allosterically induced changes in the global structures of each protein. Likewise, prominent structural perturbations were caused by the compound across the functional domains of the proteins. More so, MMC induced systematic alterations in the binding of the catalytic nucleotide and amino acid substrates which culminated in the loss of key interactions with key active site residues and ultimate reduction in the nucleotide-binding affinities across all proteins, as deduced from the binding energy calculations. These altogether confirmed that MMC uniformly disrupted the structure of the target proteins and essential substrates. Further, MMC demonstrated IC50 < 5 µM against the Dd2 and 3D7 strains of parasite making it a good starting point for malarial drug development. We believe that findings from our study will be important in the current search for highly effective multi-stage antimalarial drugs.

Investigation of Multi-Subunit Mycobacterium tuberculosis DNA-Directed RNA Polymerase and Its Rifampicin Resistant Mutants.

Emerging Mycobacterium tuberculosis (Mtb) resistant strains have continued to limit the efficacies of existing antitubercular therapies. More specifically, mutations in the RNA replicative machinery of Mtb, RNA polymerase (RNAP), have been widely linked to rifampicin (RIF) resistance, which has led to therapeutic failures in many clinical cases. Moreover, elusive details on the underlying mechanisms of RIF-resistance caused by Mtb-RNAP mutations have hampered the development of new and efficient drugs that are able to overcome this challenge. Therefore, in this study we attempt to resolve the molecular and structural events associated with RIF-resistance in nine clinically reported missense Mtb RNAP mutations. Our study, for the first time, investigated the multi-subunit Mtb RNAP complex and findings revealed that the mutations commonly disrupted structural-dynamical attributes that may be essential for the protein's catalytic functions, particularly at the ßfork loop 2, ß'zinc-binding domain, the ß' trigger loop and ß'jaw, which in line with previous experimental reports, are essential for RNAP processivity. Complementarily, the mutations considerably perturbed the RIF-BP, which led to alterations in the active orientation of RIF needed to obstruct RNA extension. Consequentially, essential interactions with RIF were lost due to the mutation-induced repositioning with corresponding reductions in the binding affinity of the drug observed in majority of the mutants. We believe these findings will significantly aid future efforts in the discovery of new treatment options with the potential to overcome antitubercular resistance.

Structure-based discovery and in vitro validation of inhibitors of chloride intracellular channel 4 protein.

The use of computer-aided methods have continued to propel accelerated drug discovery across various disease models, interestingly allowing the specific inhibition of pathogenic targets. Chloride Intracellular Channel Protein 4 (CLIC4) is a novel class of intracellular ion channel highly implicated in tumor and vascular biology. It regulates cell proliferation, apoptosis and angiogenesis; and is involved in multiple pathologic signaling pathways. Absence of specific inhibitors however impedes its advancement to translational research. Here, we integrate structural bioinformatics and experimental research approaches for the discovery and validation of small-molecule inhibitors of CLIC4. High-affinity allosteric binders were identified from a library of 1615 Food and Drug Administration (FDA)-approved drugs via a high-performance computing-powered blind-docking approach, resulting in the selection of amphotericin B and rapamycin. NMR assays confirmed the binding and conformational disruptive effects of both drugs while they also reversed stress-induced membrane translocation of CLIC4 and inhibited endothelial cell migration. Structural and dynamics simulation studies further revealed that the inhibitory mechanisms of these compounds were hinged on the allosteric modulation of the catalytic glutathione (GSH)-like site loop and the extended catalytic ß loop which may elicit interference with the catalytic activities of CLIC4. Structure-based insights from this study provide the basis for the selective targeting of CLIC4 to treat the associated pathologies.

Deciphering the molecular mechanisms of selective non-covalency demonstrated differentially by 9-Allylnaphtho[1,8-ef]isoindole-7,8,10(9H)-trione (C11) against fibroblast growth factor receptors 1-4.

The specific inhibition of aberrant Fibroblast Growth Factor Receptors (FGFRs) has been identified as a feasible strategy to therapeutically ameliorate their respective carcinogenic involvements. High homology among these proteins has however limited efforts towards the discovery of selective small-molecule compounds due to undesirable effects elicited by pan-FGFR inhibitors. A recent study showed the selective activity of a new compound C11 which was >52 times more potent against FGFR1 than FGFR2 and FGFR3, and 4 times than FGFR4. This C11 selective non-covalency was investigated in this study using computational methods since it has remained unresolved. Structural findings revealed that C11 enhanced structural perturbations in FGFR1 with less prominent effects in other FGFRs. High deviations also characterized the C11-bound active pocket of FGFR1 with notable fluctuations across the constituent P-loop, αC helix, hinge region, catalytic, and activation loops. These induced motions were essential for optimal C11 motion an d positioning of its phenalenone ring and prop-2-en-l-yl moiety at the FGFR1 active pocket to interact stably and strongly with A564FGFR1, L484FGFR1, Y563FGFR1, and E562FGFR1 which as well had high energy contributions. C11 exhibited highly unstable binding in F GFRs2-3 with a more steady interaction with FGFR4. Free binding energy (ΔGbind) analyses further estimated the highest interaction energy for C11-FGFR1 with favorable desolvation energy that indicated a deep hydrophobic pocket binding for C11 in FGFR1 compared to other FGFRs. We believe rational insights from this study will contribute to the structure-based design of highly specific FGFR1 inhibitors.Communicated by Ramaswamy H. Sarma.

From a Computational Perspective: Elucidating the Neurotherapeutic and Inhibitory properties of LRRK2 Kinase Domain by a benzothiazole-based compound.

BACKGROUND: Parkinson's disease (PD) is one of the most prominent neurodegenerative diseases hence the continual search for viable and effective treatment options. The pathogeny of PD is driven by many key proteins among which is the recently identified Leucine-rich repeated kinase 2 (LRRK2). Going forward, the onus lies on identifying small-molecule inhibitors that can halt its pathogenic involvement, and, importantly, possess the capacity to cross the blood-brain barrier (BBB). Although several compounds have been identified over the past decade for their potencies, a major limitation remains the inability of the majority to cross the blood-brain barrier (BBB). A novel series of benzothiazole-based compounds with varying LRRK2 inhibitory activities were recently synthesized, with one compound 14 (CPD14) that notably inhibited LRRK2 and promoted neuronal progenitor proliferation. METHODS: Here, we implemented molecular modelling and computational simulation methods to characterize CPD14 inhibitory mechanisms and dynamics against LRRK2. More so, we employed pharmacokinetic parameters to evaluate the biological activity and CNS-suitability of CPD14. RESULTS: Molecular dynamics evaluation revealed that CPD14 elicited disruptive effects on the secondary structure of LRRK2, including its catalytic kinase domain. Interaction analyses at the binding site further revealed crucial residues for the affinity binding and stability of CPD14, further supported by a highly favorable binding energy (ΔG). Pharmacokinetic predictions revealed the CNS-suitability of CPD14 based on its adherence to Lipinski's rule of 5 for neurogenic compounds. Also, CPD14 exhibited inhibitory tendencies against transcription proteins such as signal transducer and activation transcription (STAT) protein and STAT3; complementary mechanisms that could account for its in vitro potency. CONCLUSION: These findings, taken together, will aid the pharmacological and pharmacokinetic optimization of novel LRRK2 inhibitors for the treatment of PD.

Piecing the Fragments Together: Dynamical Insights into the Enhancement of BRD4-BD1 (BET Protein) Druggability in Cancer Chemotherapy Using Novel 8-methyl-pyrrolo[1,2-a]pyrazin-1(2H)-one Derivatives.

BACKGROUND: Fragment-based drug discovery in recent times has been explored in the design of highly potent therapeutics. METHODS: In this study, we explored the inhibitory dynamics of Compound 38 (Cpd38), a newly synthesized Bromodomain-containing protein 4 bromodomain 1 (BRD4-BD1) protein inhibitor derived from the synthetic coupling of Fragment 47 (Fgt47) into ABBV-075 scaffold. Using dynamic simulation methods, we unraveled the augmentative effects of chemical fragmentation on improved BRD4- BD1 inhibition. RESULTS: Findings from this study revealed that although Fgt47 exhibited a considerable ΔGbind, its incorporation into the difluoro-phenoxy pyridine scaffold (Cpd38) notably enhanced the binding affinity. Time-based analyses of interaction dynamics further revealed that the bulkiness of Cpd38 favored its interaction at the BRD4-BD1 active site relative to the fragment. Strikingly, compared to Fgt47, Cpd38 demonstrated high mobility, which could have enabled it to bind optimally and complementarily with key residues of the active site such as Ile146, Asn140, Cys136, Tyr98, Leu94, Val87, Phe83, and Trp81. DISCUSSION: On the contrary, the majority of these interactions were gradually lost in Fgt47, which could further indicate the essence of coupling it with the difluoro-phenoxy pyridine scaffold. Furthermore, Cpd38 had a more altering effect on BRD4-BDI relative to Fgt47, which could also be a result of its higher inhibitory activity. CONCLUSION: Conclusively, the design of highly potent therapeutics could be facilitated by the incorporation of pharmacologically active small molecule fragments into the scaffold of existing drugs.

Probing Protein-protein Interactions and Druggable Site Identification: Mechanistic Binding Events Between Ubiquitin and Zinc Finger with UFM1-specific Peptidase Domain Protein (ZUFSP).

BACKGROUND: Deubiquitinating enzymes (DUBs) protein family have been implicated in some deregulated pathways involved in carcinogeneses, such as cell cycle, gene expression, and DNA damage response (DDR). Zinc finger with UFM1-specific peptidase domain protein (ZUFSP) is one of the recently discovered members of the DUBs. OBJECTIVES: To identify and cross-validate the ZUFSP binding site using the bioinformatic tools, including SiteMap&Metapocket, respectively. To understand the molecular basis of complementary ZUFSP-Ub interaction and associated structural events using MD Simulation. METHODS: In this study, four binding pockets were predicted, characterized, and cross-validated based on physiochemical features such as site score, druggability score, site volume, and site size. Also, a molecular dynamics simulation technique was employed to determine the impact of ubiquitin-binding on ZUFSP. RESULTS: Site 1 with a site score 1.065, Size 102, D scores 1.00, and size volume 261 was predicted to be the most druggable site. Structural studies revealed that upon ubiquitin-binding, the motional movement of ZUFSP was reduced when compared to the unbound ZUFSP. Also, the ZUFSP helical arm (ZHA) domain orient in such a way that it moves closer to the Ub; this orientation enables the formation of a UBD which is very peculiar to ZUFSP. CONCLUSION: The impact of ubiquitin on ZUFSP movement and the characterization of its predicted druggable site can be targeted in the development of therapeutics.

Unveiling the mechanistic roles of chlorine substituted phthalazinone-based compounds containing chlorophenyl moiety towards the differential inhibition of poly (ADP-ribose) polymerase-1 in the treatment of lung cancer.

PARP-1 has become an attractive target in cancer treatment owing to its significant role in breast and ovarian cancers. The design of highly selective and effective poly (ADP ribose) polymerase-1 inhibitors has significant therapeutic advantages and has remained the core of several PARP-1-based drug discovery research. The pharmacophoric relevance of a chlorine substituent in a recent study led to the design of compounds 11c (meta-chlorophenyl) and 11d (para-chlorophenyl). In this study, we resolved the mechanistic effects of the changes in chlorine positional orientation, which underlie the inhibitory potencies and selectivity exhibited disparately by 11c and 11d. Compared to 11d, among other multiple higher-affinity complementary interactions with key site residues, the meta-Cl positioning in 11c facilitated its optimal motion and orientation towards conserved residues Arg878 and Asp766 with consistent pi-cation and pi-anion interactions, respectively, thereby favoring the stability of the ligand towards PARP-1. These could account for the higher inhibitory potency exhibited by 11c relative to 11d against PARP-1. The thermodynamics calculation revealed that 11c had a relatively higher total binding energy (ΔGbind) than 11d. We also observed that 11d displayed high deviations, compared to 11c, indicative of its unstable binding orientation. Furthermore, we reported in this study that the high involvement of electrostatic and van der Waal effects potentiated the binding affinity and strength of 11c (ΔEvdW = -50.58 and ΔEele = -27.20) relative to 11d (ΔEvdW = -49.46 and ΔEele = -19.96) at PARP-1 binding pocket. We believe the findings in this current study would provide valuable insights into the design of selective PARP-1 inhibitors containing chlorine substituent for cancer treatment, including lung cancer.Communicated by Ramaswamy H. Sarma.

Atomistic insights into the selective therapeutic activity of 6-(2,4-difluorophenoxy)-5-((ethylmethyl)pyridine-3-yl)-8-methylpyrrolo[1,2-a]pyrazin-1(2H)-one towards bromodomain-containing proteins.

Cross-target effect has been one of the major mechanisms of drug toxicity, this has necessitated the design of inhibitors that are specifically tailored to target particular biomolecules. 6-(2,4-difluorophenoxy)-5-((ethylmethyl)pyridine-3-yl)-8-methylpyrrolo[1,2-a] pyrazin-1(2H)-one (Cpd38) is an inhibitor possessing high inhibition rate and tailored specificity towards bromodomain-containing protein 4 (BRD4). In this research, we used an array of computational techniques to provide insight at the atomistic level the specific targeting of BRD4 by Cpd38 relative to the binding of Cpd38 with E1A binding protein P300 (EP300); another bromodomain-containing protein (BCP). Comparatively, binding of Cpd38 improved the conformational stability and compactness of BRD4 protein when compared to the Cpd38 bound EP300. Also, Cpd38 induced a conformational change in the active site of BRD4 that facilitated a complementary pose between Cpd38 and BRD4 suitable for effective atomistic interactions. Expectedly, thermodynamic calculations revealed that the Cpd38-BRD4 system had higher binding energy (-36.11 Kcal/mol) than the Cpd38-EP300 system with a free binding energy of -15.86 Kcal/mol. Noteworthy is the opposing role Trp81 (acting as hydrogen bond acceptor) and Pro1074 (acting as hydrogen bond donor) found on the WPF and LPF loops respectively play in maintaining Cpd38 stability. Furthermore, the hydrogen bond acceptor/donator ratio was approximately 4:1 in Cpd38-BRD4 system compared with 2:1 in Cpd38-EP300 system. Taken together, atomistic insights and structural perspectives detailed in this report supplements the experimental report supporting the improved selectivity of Cpd38 for BRD4 ahead of other BCPs while providing leeway for the future design of BET selective agents with better pharmacological profile.

Integrative immunoinformatics paradigm for predicting potential B-cell and T-cell epitopes as viable candidates for subunit vaccine design against COVID-19 virulence.

BACKGROUND: The increase in global mortality rates from SARS-COV2 (COVID-19) infection has been alarming thereby necessitating the continual search for viable therapeutic interventions. Due to minimal microbial components, subunit (peptide-based) vaccines have demonstrated improved efficacies in stimulating immunogenic responses by host B- and T-cells. METHODS: Integrative immunoinformatics algorithms were used to determine linear and discontinuous B-cell epitopes from the S-glycoprotein sequence. End-point selection of the most potential B-cell epitope was based on highly essential physicochemical attributes. NetCTL-I and NetMHC-II algorithms were used to predict probable MHC-I and II T-cell epitopes for globally frequent HLA-A∗O2:01, HLA-B∗35:01, HLA-B∗51:01 and HLA-DRB1∗15:02 molecules. Highly probable T-cell epitopes were selected based on their high propensities for C-terminal cleavage, transport protein (TAP) processing and MHC-I/II binding. RESULTS: Preferential epitope binding sites were further identified on the HLA molecules using a blind peptide-docking method. Phylogenetic analysis revealed close relativity between SARS-CoV-2 and SARS-CoV S-protein. LALHRSYLTPGDSSSGWTAGAA242→263 was the most probable B-cell epitope with optimal physicochemical attributes. MHC-I antigenic presentation pathway was highly favourable for YLQPRTFLL269-277 (HLA-A∗02:01), LPPAYTNSF24-32 (HLA-B∗35:01) and IPTNFTISV714-721 (HLA-B∗51:01). Also, LTDEMIAQYTSALLA865-881 exhibited the highest binding affinity to HLA-DR B1∗15:01 with core interactions mediated by IAQYTSALL870-878. COVID-19 YLQPRTFLL269-277 was preferentially bound to a previously undefined site on HLA-A∗02:01 suggestive of a novel site for MHC-I-mediated T-cell stimulation. CONCLUSION: This study implemented combinatorial immunoinformatics methods to model B- and T-cell epitopes with high potentials to trigger immunogenic responses to the S protein of SARS-CoV-2.

Prospecting the therapeutic edge of a novel compound (B12) over berberine in the selective targeting of Retinoid X Receptor in colon cancer.

The Retinoid X Receptor (RXR) is an attractive target in the treatment of colon cancer. Different therapeutic binders with high potency have been used to specifically target RXR. Among these compounds is a novel analogue of berberine, B12. We provided structural and molecular insights into the therapeutic activity properties of B12 relative to its parent compound, berberine, using force field estimations and thermodynamic calculations. Upon binding of B12 to RXR, the high instability elicited by RXR was markedly reduced; similar observation was seen in the berberine-bound RXR. However, our analysis revealed that B12 could have a more stabilizing effect on RXR when compared to berberine. Interestingly, the mechanistic behaviour of B12 in the active site of RXR opposed its impact on RXR protein. This disparity could be due to the bond formation and breaking elicited between B12/berberine and the active site residues. We observed that B12 and berberine could induce a disparate conformational change in regions Gly250-Asp258 located on the His-RXRα/LBD domain. Comparatively, the high agonistic and activation potential reported for B12 compared to berberine might be due to its superior binding affinity as evidenced in the thermodynamic estimations. The total affinity for B12 (-25.76 kcal/mol) was contributed by electrostatic interactions from Glu243 and Glu239. Also, Arg371, which plays a crucial role in the activity of RXR, formed a strong hydrogen bond with B12; however, a weak interaction was elicited between Arg371 and berberine. Taken together, our study has shown the RXRα activating potential of B12, and findings from this study could provide a framework in the future design of RXRα binders specifically tailored in the selective treatment of colon cancer.

Insight into the Therapeutic Potential of a Bicyclic Hydroxypyridone Compound 2-[(2,4-Dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one as COMT Inhibitor in the Treatment of Parkinson's Disease: A Molecular Dynamic Simulation Approach.

Parkinson's disease (PD) is one of the most targeted neurodegenerative diseases in clinical research. Awareness of research is due to its increasing number of affected people worldwide. The pathology of PD has been linked to several key proteins upregulation such as the catechol O-Methyltransferase (COMT). Hence, the synthesis of compounds possessing inhibitory capacity has been the frontline of research in recent years. Several compounds have been synthesized among which is the nitrocatechol. However, major limitations associated with the nitrocatechol scaffold include the inability to possess adequate CNS penetration properties and hepatic toxicity associated with the compounds. However, a series of bicyclic hydroxypyridones compounds were synthesized to evaluate their inhibitory potentials on COMT protein with compound 38 (c38) 2-[(2,4-dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one shown to have a 40 fold increase level coverage in its IC50 over brain exposure when compared to the other synthesized compound. The molecular dynamics method was employed to understand the nature of interaction exhibited by c38. Molecular mechanics of c38 revealed a disruptive effect on the secondary structure of COMT protein. Per residue decomposition analysis revealed similar crucial residues involved in the favorable binding of c38 and tolcapone implicated its increased inhibitory capacity on COMT in preventing PD. Free binding energy (ΔGbind ) of c38 further revealed the inhibitory capacity towards COMT protein in comparison to the FDA approved tolcapone. Ligand mobility analysis of both compounds showed a timewise different mobility pattern across the simulation time frame at the active site pocket of the protein connoting the different inhibitory potency exhibited by c38 and tolcapone. Findings from this study revealed optimization of c38 could facilitate the discovery of new compounds with enhanced inhibitory properties towards COMT in treating PD.

Immunoinformatics prediction of potential B-cell and T-cell epitopes as effective vaccine candidates for eliciting immunogenic responses against Epstein-Barr virus.

BACKGROUND: The ongoing search for viable treatment options to curtail Epstein Barr Virus (EBV) pathogenicity has necessitated a paradigmatic shift towards the design of peptide-based vaccines. Potential B-cell and T-cell epitopes were predicted for nine antigenic EBV proteins that mediate epithelial cell-attachment and spread, capsid self-assembly, DNA replication and processivity. METHODS: Predictive algorithms incorporated in the Immune Epitope Database (IEDB) resources were used to determine potential B-cell epitopes based on their physicochemical attributes. These were combined with a string-kernel method and an antigenicity predictive AlgPred tool to enhance accuracy in the end-point selection of highly potential antigenic EBV B-cell epitopes. NetCTL 1.2 algorithms enabled the prediction of probable T-cell epitopes which were structurally modeled and subjected to blind peptide-protein docking with HLA-A*02:01. All-atom molecular dynamics (MD) simulation and Molecular Mechanics Generalized-Born Surface Area methods were used to investigate interaction dynamics and affinities of predicted T-cell peptide-protein complexes. RESULTS: Computational predictions and sequence overlapping analysis yielded 18 linear (continuous) and discontinuous (conformational) subunit epitopes from the antigenic proteins with characteristic surface accessibility, flexibility and antigenicity, and predictive scores above the threshold value (1) set. A novel site was identified on HLA-A*02:01 with preferential affinity binding for modeled BMRF2, BXLF1 and BGLF4 T-cell epitopes. Interaction dynamics and energies were also computed in addition to crucial residues that mediated complex formation and stability. CONCLUSION: This study implemented an integrative meta-analytical approach to model highly probable B-cell and T-cell epitopes as potential peptide-vaccine candidates for the treatment of EBV-related diseases.

Piece of the puzzle: Remdesivir disassembles the multimeric SARS-CoV-2 RNA-dependent RNA polymerase complex.

The recently emerged SARS-like coronavirus (SARS-CoV-2) has continued to spread rapidly among humans with alarming upsurges in global mortality rates. A major key to tackling this virus is to disrupt its RNA replication process as previously reported for Remdesivir (Rem-P3). In this study, we theorize, using computational simulations, novel mechanisms that may underlie the binding of Rem-P3 to SARS-CoV-2 RdRp-NSPs complex; a multimeric assembly that drives viral RNA replication in human hosts. Findings revealed that while ATP-binding stabilized the replicative tripartite, Rem-P3 disintegrated the RdRp-NSP complex, starting with the detachment of the NSP7-NSP8 heterodimer followed by minimal displacement of the second NSP8 subunit (NSP8II). More so, Rem-P3 interacted with a relatively higher affinity (ΔGbind) while inducing high perturbations across the RdRp-NSP domains. D452, T556, V557, S682, and D760 were identified for their crucial roles in stacking the cyano-adenosine and 3,4-dihydroxyoxolan rings of Rem-P3 while its flexible P3 tail extended towards the palm domain blocking D618 and K798; a residue-pair identified for essential roles in RNA replication. However, ATP folded away from D618 indicative of a more coordinated binding favorable for nucleotide polymerization. We believe findings from this study will significantly contribute to the structure-based design of novel disruptors of the SARS-CoV-2 RNA replicative machinery.

East to West not North-West: Structure-Based Mechanistic Resolution of 8-Hydroxyl Replacement and Resulting Effects on the Activities of Imidazole-Based Heme Oxygenase-1 Inhibitors.

Upregulation of Heme Oxygenase-1 (HO-1) has been widely implicated in cancer growth and chemoresistance. This explains the numerous drug discovery efforts aimed at mitigating its pro-carcinogenic roles till date. In a recent study, two selective azole-based HO-1 inhibitors (Cpd1 and Cpd2) were synthesized, which exhibited differential inhibitory potencies of ~200µm. Interestingly, variations in the affinities of these compounds were determined by their positioning across specific regions of the HO-1 binding domain, pin-pointing a pharmacological interrelationship that remains unresolved. Therefore, in this study, using molecular dynamics simulations and binding free energy calculations, we investigate how dynamical orientations of these compounds influence their binding affinities at the active HO-1 domain. Findings revealed favorable binding for the bromobenzene and imidazole substituents of Cpd1 at the western and eastern regions of the HO-1 active domain. The constituent hydroxyl group was coordinated by residues Asp140 and Arg136 over the simulation period. On the contrary, stable binding of the bromobenzene and imidazole substituents were negated by the optimal orientations of the benzyl substituent, which extended into the northeastern region. These were supported by the displacement of Asp140 and Arg136, crucial for hydrogen bond formation in Cpd1. Also, we observed that Cpd2 exhibited high deviations indicative of an unstable binding relative to Cpd1. This further supports the presumption that Cpd2 was systematically oriented away from the active HO-1 region, a phenomenon that was due to the optimal motions of the benzyl group at the northeastern regions. The highlight of our findings is that the benzyl substituent in Cpd2 elicited negative effects on HO-1, vis a vis, instability, displacement of crucial residues, and low binding energy when compared to Cpd1. Findings pave the way for future drug discovery efforts related to HO-1 inhibition in cancer therapy.

Exploring the structural basis and atomistic binding mechanistic of the selective antagonist blockade at D3 dopamine receptor over D2 dopamine receptor.

More recently, there has been a paradigm shift toward selective drug targeting in the treatment of neurological disorders, including drug addiction, schizophrenia, and Parkinson's disease mediated by the different dopamine receptor subtypes. Antagonists with higher selectivity for D3 dopamine receptor (D3DR) over D2 dopamine receptor (D2DR) have been shown to attenuate drug-seeking behavior and associated side effects compared to non-subtype selective antagonists. However, high conservations among constituent residues of both proteins, particularly at the ligand-binding pockets, remain a challenge to therapeutic drug design. Recent studies have reported the discovery of two small-molecules R-VK4-40 and Y-QA31 which substantially inhibited D3DR with >180-fold selectivity over D2DR. Therefore, in this study, we seek to provide molecular and structural insights into these differential binding mechanistic using meta-analytic computational simulation methods. Findings revealed that R-VK4-40 and Y-QA31 adopted shallow binding modes and were more surface-exposed at D3DR while on the contrary, they exhibited deep hydrophobic pocket binding at D2DR. Also, two non-conserved residues; Tyr361.39 and Ser18245.51 were identified in D3DR, based on their crucial roles and contributions to the selective binding of R-VK4-40 and Y-QA31. Importantly, both antagonists exhibited high affinities in complex with D3DR compared to D2DR, while van der Waals energies contributed majorly to their binding and stability. Structural analyses also revealed the distinct stabilizing effects of both compounds on D3DR secondary architecture relative to D2DR. Therefore, findings herein pinpointed the origin and mechanistic of selectivity of the compounds, which may assist in the rational design of potential small molecules of the D2 -like dopamine family receptor subtype with improved potency and selectivity.

A probable means to an end: exploring P131 pharmacophoric scaffold to identify potential inhibitors of Cryptosporidium parvum inosine monophosphate dehydrogenase.

Compound P131 has been established to inhibit Cryptosporidium parvum's inosine monophosphate dehydrogenase (CpIMPDH). Its inhibitory activity supersedes that of paromomycin, which is extensively used in treating cryptosporidiosis. Through the per-residue energy decomposition approach, crucial moieties of P131 were identified and subsequently adopted to create a pharmacophore model for virtual screening in the ZINC database. This search generated eight ADMET-compliant hits that were examined thoroughly to fit into the active site of CpIMPDH via molecular docking. Three compounds ZINC46542062, ZINC58646829, and ZINC89780094, with favorable docking scores of - 8.3 kcal/mol, - 8.2 kcal/mol, and - 7.5 kcal/mol, were selected. The potential inhibitory mechanism of these compounds was probed using molecular dynamics simulation and Molecular Mechanics Generalized Poisson Boltzmann Surface Area (MM/PBSA) analyses. Results revealed that one of the hits (ZINC46542062) exhibited a lower binding free energy of - 39.52 kcal/mol than P131, which had - 34.6 kcal/mol. Conformational perturbation induced by the binding of the identified hits to CpIMPDH was similar to P131, suggesting a similarity in inhibitory mechanisms. Also, in silico investigation of the properties of the hit compounds implied superior physicochemical properties with regards to their synthetic accessibility, lipophilicity, and number of hydrogen bond donors and acceptors in comparison with P131. ZINC46542062 was identified as a promising hit compound with the highest binding affinity to the target protein and favorable physicochemical and pharmacokinetic properties relative to P131. The identified compounds can serve as a basis for conducting further experimental investigations toward the development of anticryptosporidials, which can overcome the challenges of existing therapeutic options. Graphical abstract P131 and the identified compounds docked in the NAD+ binding site of Cryptosporidium parvum IMPDH.

A Mechanistic Probe into the Dual Inhibition of T. cruzi Glucokinase and Hexokinase in Chagas Disease Treatment - A Stone Killing Two Birds?

Glucokinase (GLK) and Hexokinase (HK) have been characterized as essential targets in Trypanosoma cruzi (Tc)-mediated infection. A recent study reported the propensity of the concomitant inhibition of TcGLK and TcHK by compounds GLK2-003 and GLK2-004, thereby presenting an efficient approach in Chagas disease treatment. We investigated this possibility using atomic and molecular scaling methods. Sequence alignment of TcGLK and TcHK revealed that both proteins shared approximately 33.3 % homology in their glucose/inhibitor binding sites. The total binding free energies of GLK2-003 and GLK2-004 were favorable in both proteins. PRO92 and THR185 were pivotal to the binding and stabilization of the ligands in TcGLK, likewise their conserved counterparts, PRO163 and THR237 in TcHK. Both compounds also induced a similar pattern of perturbations in both TcGLK and TcHK secondary structure. Findings from this study therefore provide insights into the underlying mechanisms of dual inhibition exhibited by the compounds. These results can pave way to discover and optimize novel dual Tc inhibitors with favorable pharmacokinetics properties eventuating in the mitigation of Chagas disease.

Structural alterations in the catalytic core of hSIRT2 enzyme predict therapeutic benefits of Garcinia mangostana derivatives in Alzheimer's disease: molecular dynamics simulation study.

Recent studies have shown that inhibition of the hSIRT2 enzyme provides favorable effects in neurodegenerative diseases such as Alzheimer's disease. Prenylated xanthone phytochemicals including α-mangostin, ß-mangostin and γ-mangostin obtained from Garcinia mangostana, a well-established tropical plant, have been shown experimentally to inhibit sirtuin enzymatic activity. However, the molecular mechanism of this sirtuin inhibition has not been reported. Using comprehensive integrated computational techniques, we provide molecular and timewise dynamical insights into the structural alterations capable of facilitating therapeutically beneficial effects of these phytochemicals at the catalytic core of the hSIRT2 enzyme. Findings revealed the enhanced conformational stability and compactness of the hSIRT2 catalytic core upon binding of γ-mangostin relative to the apoenzyme and better than α-mangostin and ß-mangostin. Although thermodynamic calculations revealed favorable binding of all the phytochemicals to the hSIRT2 enzyme, the presence of only hydroxy functional groups on γ-mangostin facilitated the occurrence of additional hydrogen bonds involving Pro115, Phe119, Asn168 and His187 which are absent in α-mangostin- and ß-mangostin-bound systems. Per-residue energy contributions showed that van der Waals and more importantly electrostatic interactions are involved in catalytic core stability with Phe96, Tyr104 and Phe235 notably contributing π-π stacking, π-π T shaped and π-sigma interactions. Cumulatively, our study revealed the structural alterations leading to inhibition of hSIRT2 catalysis and findings from this study could be significantly important for the future design and development of sirtuin inhibitors in the management of Alzheimer's disease.

Molecular Basis of P131 Cryptosporidial-IMPDH Selectivity-A Structural, Dynamical and Mechanistic Stance.

Cryptosporidiosis accounts for a surge in infant (<5 years) mortality and morbidity. To date, several drug discovery efforts have been put in place to develop effective therapeutic options against the causative parasite. Based on a recent report, P131 spares inosine monophosphate dehydrogenase (IMPDH) in a eukaryotic model (mouse IMPDH (mIMPDH)) while binding selectively to the NAD+ site in Cryptosporidium parvum (CpIMPDH). However, no structural detail exists on the underlining mechanisms of P131-CpIMPDH selective targeting till date. To this effect, we investigate the selective inhibitory dynamics of P131 in CpIMPDH relative to mIMPDH via molecular biocomputation methods. Pairwise sequence alignment revealed prominent variations at the NAD+ binding regions of both proteins that accounted for disparate P131 binding activities. The influence of these variations was further revealed by the MM/PBSA energy estimations coupled with per-residue energy decomposition which monitored the systematic binding of the compound. Furthermore, relative high-affinity interactions occurred at the CpIMPDH NAD+ site which were majorly mediated by SER22, VAL24, PRO26, SER354, GLY357, and TYR358 located on chain D. These residues are unique to the parasite IMPDH form and not in the eukaryotic protein, highlighting variations that account for preferential P131 binding. Molecular insights provided herein corroborate previous experimental reports and further underpin the basis of CpIMPDH inhibitor selectivity. Findings from this study could present attractive prospects toward the design of novel anticryptosporidials with improved selectivity and binding affinity against parasitic targets.