Multi-cavity molecular descriptor interconnections: Enhanced protocol for prediction of serum albumin drug binding.

The role of human serum albumin (HSA) in the transport of molecules predicates its involvement in the determination of drug distribution and metabolism. Optimization of ADME properties are analogous to HSA binding thus this is imperative to the drug discovery process. Currently, various in silico predictive tools exist to complement the drug discovery process, however, the prediction of possible ligand-binding sites on HSA has posed several challenges. Herein, we present a strong and deeper-than-surface case for the prediction of HSA-ligand binding sites using multi-cavity molecular descriptors by exploiting all experimentally available and crystallized HSA-bound drugs. Unlike previously proposed models found in literature, we established an in-depth correlation between the physicochemical properties of available crystallized HSA-bound drugs and different HSA binding site characteristics to precisely predict the binding sites of investigational molecules. Molecular descriptors such as the number of hydrogen bond donors (nHD), number of heteroatoms (nHet), topological polar surface area (TPSA), molecular weight (MW), and distribution coefficient (LogD) were correlated against HSA binding site characteristics, including hydrophobicity, hydrophilicity, enclosure, exposure, contact, site volume, and donor/acceptor ratio. Molecular descriptors nHD, TPSA, LogD, nHet, and MW were found to possess the most inherent capacities providing baseline information for the prediction of serum albumin binding site. We believe that these associations may form the bedrock for establishing a solid correlation between the physicochemical properties and Albumin binding site architecture. Information presented in this report would serve as critical in provisions of rational drug designing as well as drug delivery, bioavailability, and pharmacokinetics.

Apigenin attenuates depressive-like behavior via modulating monoamine oxidase A enzyme activity in chronically stressed mice.

Chronic stress is a risk factor for depression and is characterized by elevated levels of brain monoamine oxidase A (MAOA). Mounting evidence has shown that MAOA is a biochemical link between stress and depression. Apigenin (API), a natural flavonoid, as demonstrated in vitro inhibitory effect on MAOA, is suggestive of antidepressant-like activity. However, the in vivo inhibitory effect of API on MAOA and how it affects depression still remain unclear. Here, we report the probable mechanisms of action of API in chronic unpredictable mild stress (CUMS)-induced depression in mice. Treatment with API reversed anhedonia, and reduced anxiety and immobility time in behavioral studies. API reduced brain corticosterone and malondialdehyde (MDA) levels but increased brain levels of glutathione and superoxide dismutase. Furthermore, interleukin-6 and tumor necrosis factor-α were attenuated by API. It also restored cell loss and inhibited the activity of MAOA in the hippocampal brain regions and prefrontal cortex. Comparative binding affinity of API for MAOA (-7.7 kcal/mol) through molecular docking studies was greater than that of reference compound, clorgyline (-6.8 kcal/mol). Favorable hydrophobic interactions important to API binding at MAOA binding cavity was revealed to include conventional hydrogen bond (Cys323 and Tyr444), π-Sulfur (Cys323), π-π Stacked (Tyr407), π-π T-shaped (Phe208), π-lone pair and π-alkyl (Ile335, Ile180) interactions. These results suggest that API is a potent, selective, reversible inhibitor of MAOA with capability of attenuating CUMS-induced depression via inhibiting MAOA enzyme activity and altering other pathomechanisms.

Talazoparib Dual-targeting on Poly (ADP-ribose) Polymerase-1 and -16 Enzymes Offers a Promising Therapeutic Strategy in Small Cell Lung Cancer Therapy: Insight from Biophysical Computations.

In recent times, inhibition of poly (ADP-ribose) polymerase (PARP) enzymes by pharmacological drugs has attracted much attention as an anticancer therapy. As reported, PARP-16 has been discovered as a novel anticancer target for small cell lung cancer, and that the inhibition of both PARP-16 and PARP-1 by talazoparib can increase the overall effectiveness of talazoparib in the SCLC treatment. In this study, we employed computational approaches to investigate the differential inhibitory potency of Talazoparib on PARP-1 and PARP-16. Talazoparib has excellent PARP-1 and PARP-16 binding activities, as revealed by the ΔGbind (total binding energy). Pp16-tpb had binding energy of -34.85 kcal/mol, while pp1-tpb had a binding energy of -26.36 kcal/mol. The binding activity of Talazoparib on both PARP-1 and PARP-16 was significantly influenced by van der Waal and electrostatic interactions. Correspondingly, according to the findings of this study, binding residues with total binding energy greater than 1.00 kcal/mol contributed considerably to the Talazoparib's binding activities on PARP-1 and PARP-16. We believe the findings of this study will pave the way for developing dual targeting of PARP enzymes as a strategy for small-cell lung cancer treatment.

Simulation Models for Prediction of Bioavailability of Medicinal Drugs-the Interface Between Experiment and Computation.

The oral drug bioavailability (BA) problems have remained inevitable over the years, impairing drug efficacy and indirectly leading to eventual human morbidity and mortality. However, some conventional lab-based methods improve drug absorption leading to enhanced BA, and the recent experimental techniques are up-and-coming. Nevertheless, some have inherent drawbacks in improving the efficacy of poorly insoluble and low impermeable drugs. Drug BA and strategies to overcome these challenges were briefly highlighted. This review has significantly unravelled the different computational models for studying and predicting drug bioavailability. Several computational approaches provide mechanistic insights into the oral drug delivery system simulation of descriptors like solubility, permeability, transport protein-ligand interactions, and molecular structures. The in silico techniques have long been known still are just being applied to unravel drug bioavailability issues. Many publications have reported novel applications of the computational models towards achieving improved drug BA, including predicting gastrointestinal tract (GIT) drug absorption properties and passive intestinal membrane permeability, thus maximizing time and resources. Also, the classical molecular simulation models for free solvation energies of soluble-related processes such as solubilization, dissolutions, supersaturation, and precipitation have been used in virtual screening studies. A few of the tools are GastroPlusTM that supports biowaiver for drugs, mainly BCS class III and predicts drug compounds' absorption and pharmacokinetic process; SimCyp® simulator for mechanistic modelling and simulation of drug formulation processes; pharmacodynamics analysis on non-linear mixed-effects modelling; and mathematical models, predicting absorption potential/maximum absorption dose. This review provides in silico-experiment annexation in the drug bioavailability enhancement, possible insights that lead to critical opinion on the applications and reliability of the various in silico models as a growing tool for drug development and discovery, thus accelerating drug development processes.

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.

Unravelling the Mechanistic Role of Quinazolinone Pharmacophore in the Inhibitory Activity of Bis-quinazolinone Derivative on Tankyrase-1 in the Treatment of Colorectal Cancer (CRC) and Non-small Cell Lung Cancer (NSCLC): A Computational Approach.

In recent years, tankyrase inhibition has gained a great focus as an anti-cancer strategy due to their modulatory effect on WNT/ß-catenin pathway implicated in many malignancies, including colorectal cancer (CRC) and non-small cell lung cancer (NSCLC). Based on the structural homology in the catalytic domain of PARP enzymes, bis-quinazolinone 5 (Cpd 5) was designed to be a potent selective tankyrase inhibitor. In this study, we employed molecular dynamics simulations and binding energy analysis to decipher the underlying mechanism of TNK-1 inhibition by Cpd 5 in comparison with a known selective tankyrase, IWR-1. The Cpd 5 had a relatively higher ΔGbind than IWR-1 from the thermodynamics analysis, revealing the better inhibitory activity of Cpd 5 compared to IWR-1. High involvement of solvation energy (ΔGsol) and the van der Waals energy (ΔEvdW) potentiated the affinity of Cpd 5 at TNK-1 active site. Interestingly, the keto group and the N3 atom of the quinazolinone nucleus of Cpd 5, occupying the NAM subsite, was able to form H-bond with Gly1185, thereby favoring the better stability and higher inhibitory efficacy of Cpd 5 relative to IWR-1. Our analysis proved that the firm binding of Cpd 5 was achieved by the quinazolinone groups via the hydrophobic interactions with the side chains of key site residues at the two subsite regions: His1201, Phe1188, Ala1191, and Ile1192 at the AD subsite and Tyr1224, Tyr1213, and Ala1215 at the NAM subsite. Thus, Cpd 5 is dominantly bound through π-π stacked interactions and other hydrophobic interactions. We believe that findings from this study would provide an important rationale towards the structure-based design of improved selective tankyrase inhibitors in cancer therapy.

Mechanistic Insights into the Selective Dual BET and PLK1 Inhibitory Activity of a Novel Benzamide Compound in Castration-Resistant Prostrate Cancer.

Though multifactorial, BET and PLK1 proteins have been found to be key players in the oncogenic process leading to castration-resistant prostate cancer through regulation of AR and MYC-mediated transcription. Hence, dual inhibition of these proteins appears to be an auspicious approach for CRPC therapy. WNY0824 has been reported to exhibit nanomolar range inhibition as well as significant anti-proliferative activity on AR-positive CRPC cells in vitro. However, structural, and mechanistic events associated with its dual inhibitory and anti-proliferative mechanisms remain unclear. Utilizing integrative computer-assisted atomistic techniques, analyses revealed that the dual-inhibitory activity of WNY0824 against BRD4 and PLK1 proteins is mediated by conserved residues present in the binding cavities of both proteins which are shown to elicit various strong intermolecular interactions and thus favour binding affinity. Also, binding orientation of the ligand at the protein binding cavities allowed for important hydrophobic interactions which resulted in high binding free energy of -42.50 kcal/mol and -51.64 kcal/mol towards BRD4 and PLK1, respectively. While van der Waals interactions are very important to ligand binding in BRD4-WNY complex, electrostatic interactions are pertinent to PLK1-WNY complex. Intriguingly, WNY0824 triggered conformational alterations in both proteins through increased structural instability, decreased structural compactness and mitigation in exposure of residues to solvent surface area. Consequently, critical interactions peculiar to the oncogenic activities of BRD4 and PLK1 were inhibited, a phenomenon that results in an antagonism of CRPC progression. The mechanistic insights presented in this report would further assist in the structure-based design of improved inhibitors useful in CRPC therapy.

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.

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.

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.