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.

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.

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.

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.

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.

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.

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.

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.

Design, synthesis, and evaluation of "dual-site"-binding diarylpyrimidines targeting both NNIBP and the NNRTI adjacent site of the HIV-1 reverse transcriptase.

Inspired by our previous efforts to improve the drug-resistance profiles of HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs), a novel series of "dual-site" binding diarylpyrimidine (DAPY) derivatives targeting both the NNRTI adjacent site and NNRTIs binding pocket (NNIBP) were designed, synthesized, and evaluated for their anti-HIV potency in TZM-bl and MT-4 cells. Eight compounds exhibited moderate to excellent potencies in inhibiting wild-type (WT) HIV-1 replication with EC50 values ranging from 2.45 nM to 5.36 nM, and 14c (EC50 = 2.45 nM) proved to be the most promising inhibitor. Of note, 14c exhibited potent activity against the single mutant strain E138K (EC50 = 10.6 nM), being comparable with ETR (EC50 = 9.80 nM) and 3.5-fold more potent than that of compound 7 (EC50 = 37.3 nM). Moreover, 14c acted as a classical NNRTI with high affinity for WT HIV-1 RT (IC50 = 0.0589 µM). The detailed structure-activity relationships (SARs) of the representative compounds were also determined, and further supported by molecular dynamics simulation. Overall, we envision that the "dual-site"-binding NNRTIs have significant prospects and pave the way for the next round of rational design of potent anti-HIV-1 agents.

Exploiting the tolerant region I of the non-nucleoside reverse transcriptase inhibitor (NNRTI) binding pocket. Part 2: Discovery of diarylpyrimidine derivatives as potent HIV-1 NNRTIs with high Fsp3 values and favorable drug-like properties.

To yield potent HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) with favorable drug-like properties, a series of novel diarylpyrimidine derivatives targeting the tolerant region I of the NNRTI binding pocket were designed, synthesized and biologically evaluated. The most active inhibitor 10c exhibited outstanding antiviral activity against most of the viral panel, being about 2-fold (wild-type, EC50 = 0.0021 µM), 1.7-fold (K103N, EC50 = 0.0019 µM), and slightly more potent (E138K, EC50 = 0.0075 µM) than the NNRTI drug etravirine (ETR). Additionally, 10c was endowed with relatively low cytotoxicity (CC50 = 18.52 µM). More importantly, 10c possessed improved drug-like properties compared to those of ETR with an increased Fsp3 (Fraction of sp3 carbon atoms) value. Furthermore, the molecular dynamics simulation and molecular docking studies were implemented to reveal the binding mode of 10c in the binding pocket. Taken together, 10c is a promising lead compound that is worth further investigation.

Revealing the role of fluorine pharmacophore in chalcone scaffold for shifting the MAO-B selectivity: investigation of a detailed molecular dynamics and quantum chemical study.

The development of highly selective monoamine oxidase-B (MAO-B) inhibitors has great therapeutic benefit in treatment of various neurodegenerative disorders. Recent study documented that shifting of fluorine atom from para to ortho position on the phenyl B ring of heteroaryl chalcones shown a remarkable shift in the selectivity and potency between MAO-A and MAO-B isoforms. Despite the large plethora of the design of new selective MAO-B inhibitors, the current paper illustrates the role and orientation of fluorine atom with remarkable MAO-B selectivity of three compounds (O23, O24 and O25), which differ from all other substituents encountered in the chalcone scaffolds is recently reported by our group. Conformational analyses of differential inhibitory effects of O23, O24 and O25 on MAO-A and MAO-B, differential analyses of complementary interactions at MAO-A/-B active sites and differential analysis of affinity binding and per-residue energy contributions are calculated by molecular dynamics study. Density functional theory based electronic structure calculations were employed with special emphasis to electrostatic potential and frontier molecular orbitals. Results of the current study can be used for lead modification and a new insight for the development of novel fluorinated chalcones for the treatment of various neurodegenerative disorders. Communicated by Ramaswamy H. Sarma.

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.

Investigating the Mechanistic Inhibitory Discrepancies of Novel Halogen and Alkyl Di-Substituted Oxadiazole-Based Dibenzo-Azepine-Dione Derivatives on Poly (ADP-Ribose) Polymerase-1.

Numerous studies have established the involvement of Poly (ADP-ribose) Polymerase-1 (PARP-1) in cancer presenting it as an important therapeutic target over recent years. Although homology among the PARP protein family makes selective targeting difficult, two compounds [d11 (0.939 µM) and d21 (0.047 µM)] with disparate inhibitory potencies against PARP-1 were recently identified. In this study, free energy calculations and molecular simulations were used to decipher underlying mechanisms of differential PARP-1 inhibition exhibited by the two compounds. The thermodynamics calculation revealed that compound d21 had a relatively higher ΔGbind than d11. High involvement of van der Waal and electrostatic effects potentiated the affinity of d21 at PARP-1 active site. More so, incorporated methyl moiety in d11 accounted for steric hindrance which, in turn, prevented complementary interactions of key site residues such as TYR889, MET890, TYR896, TYR907. Conformational studies also revealed that d21 is more stabilized for interactions in the active site compared to d11. We believe that findings from this study would provide an important avenue for the development of selective PARP-1 inhibitors.

Therapeutic Path to Double Knockout: Investigating the Selective Dual-Inhibitory Mechanisms of Adenosine Receptors A1 and A2 by a Novel Methoxy-Substituted Benzofuran Derivative in the Treatment of Parkinson's Disease.

The dual inhibition of adenosine receptors A1 (A1 AR) and A2 (A2A AR) has been considered as an efficient strategy in the treatment of Parkinson's disease (PD). This led to the recent development of a series of methoxy-substituted benzofuran derivatives among which compound 3j exhibited dual-inhibitory potencies in the micromolar range. Therefore, in this study, we seek to resolve the mechanisms by which this novel compound elicits its selective dual targeting against A1 AR and A2A AR. Unique to the binding of 3j in both proteins, from our findings, is the ring-ring interaction elicited by A1Phe275 (→ A2Phe170) with the benzofuran ring of the compound. As observed, this π-stacking interaction contributes notably to the stability of 3j at the active sites of A1 and A2A AR. Besides, conserved active site residues in the proteins such as A1Ala170 (→ A2Ala65), A1Ile173 (→ A2Ile68), A1Val191 (→ A2Val86), A1Leu192 (→ A2Leu87), A1Ala195 (→ A2Ala90), A1Met284 (→ A2Met179), A1Tyr375 (→ A2Tyr369), A1Ile378 (→ A2Ile372), and A1His382 (→ A2His376) were commonly involved with other ring substituents which further complement the dual binding and stability of 3j. This reflects a similar interaction mechanism that involved aromatic (π) interactions. Consequentially, vdW energies contributed immensely to the dual binding of the compound, which culminated in high ΔGbinds that were homogenous in both proteins. Furthermore, 3j commonly disrupted the stable and compact conformation of A1 and A2A AR, coupled with their active sites where Cα deviations were relatively high. Ligand mobility analysis also revealed that both compounds exhibited a similar motion pattern at the active site of the proteins relative to their optimal dual binding. We believe that findings from this study with significantly aid the structure-based design of highly selective dual-inhibitors of A1 and A2A AR.