Exploring the ring potential of 2,4-diaminopyrimidine derivatives towards the identification of novel caspase-1 inhibitors in Alzheimer's disease therapy

Molecular Modeling Identification of Key Secondary Metabolites from Xylopia aethiopica as Promising Therapeutics Targeting Essential Measles Viral Proteins

This study computationally screened three key compounds (vanillin (VAN), oxophoebine (OPB), and dihydrochalcone (DHC)) derived from Xylopia aethiopica (Guinea pepper), a medicinal plant with known antiviral activity, against key druggable measles virus (MV) proteins (fusion protein (FUP), haemagglutinin protein (HMG), and phosphoprotein (PSP)). Each molecular species was subjected to a 100 ns molecular dynamics (MD) simulation following docking, and a range of postdynamic parameters including free binding energy and pharmacokinetic properties were determined. The docking scores of the resulting OPB-FUP (-5.4 kcal/mol), OPB-HMG (-8.1 kcal/mol), and OPB-PSP (-8.0 kcal/mol) complexes were consistent with their respective binding energy values (-25.37, -28.74, and -40.68 kcal/mol), and higher than that of the reference standard, ribavirin (RBV) in each case. Furthermore, all the investigated compounds were thermodynamically compact and stable, especially HMG of MV, and this observation could be attributed to the resulting intermolecular interactions in each system. Overall, OPB may possess inhibitory properties against MV glycoproteins (FUP and HMG) and PSP that play important roles in the replication of MV and measles pathogenesis. While OPB could serve as a scaffold for the development of novel MV fusion and entry inhibitors, further in vitro and in vivo evaluation is highly recommended.

Bridging the Gap in Malaria Parasite Resistance, Current Interventions, and the Way Forward from in Silico Perspective: A Review

The past decade has seen most antimalarial drugs lose their clinical potency stemming from parasite resistance. Despite immense efforts by researchers to mitigate this global scourge, a breakthrough is yet to be achieved, as most current malaria chemotherapies suffer the same fate. Though the etiology of parasite resistance is not well understood, the parasite's complex life has been implicated. A drug-combination therapy with artemisinin as the central drug, artemisinin-based combination therapy (ACT), is currently the preferred malaria chemotherapy in most endemic zones. The emerging concern of parasite resistance to artemisinin, however, has compromised this treatment paradigm. Membrane-bound Ca²⁺-transporting ATPase and endocytosis pathway protein, Kelch13, among others, are identified as drivers in plasmodium parasite resistance to artemisinin. To mitigate parasite resistance to current chemotherapy, computer-aided drug design (CADD) techniques have been employed in the discovery of novel drug targets and the development of small molecule inhibitors to provide an intriguing alternative for malaria treatment. The evolution of plasmepsins, a class of aspartyl acid proteases, has gained tremendous attention in drug discovery, especially the non-food vacuole. They are expressed at multi-stage of the parasite's life cycle and involve in hepatocytes' egress, invasion, and dissemination of the parasite within the human host, further highlighting their essentiality. In silico exploration of non-food vacuole plasmepsin, PMIX and PMX unearthed the dual enzymatic inhibitory mechanism of the WM382 and 49c, novel plasmepsin inhibitors presently spearheading the search for potent antimalarial. These inhibitors impose structural compactness on the protease, distorting the characteristic twist motion. Pharmacophore modeling and structure activity of these compounds led to the generation of hits with better affinity and inhibitory prowess towards PMIX and PMX. Despite these headways, the major obstacle in targeting PM is the structural homogeneity among its members and to human Cathepsin D. The incorporation of CADD techniques described in the study at early stages of drug discovery could help in selective inhibition to augment malaria chemotherapy.

Network Pharmacology- and Molecular Dynamics Simulation-Based Bioprospection of Aspalathus linearis for Type-2 Diabetes Care

The medicinal herb Aspalathus linearis (rooibos) is globally recognized in type-2 diabetes mellitus (T2DM) treatment due to its known and distinctive compounds. This work utilized network pharmacology (NP) coupled with molecular dynamics simulation in gaining new insight into the anti-diabetic molecular mechanism of action of rooibos teas. It looked at the interactions between rooibos constituents with various relevant protein receptors and signaling routes associated with T2DM progression. The initial analysis revealed 197 intersecting gene targets and 13 bioactive rooibos constituents linked to T2DM. The interactions between proteins and compounds to the target matrix were generated with the Cystoscope platform and STRING database. These analyses revealed intersecting nodes active in T2DM and hypoxia-inducible factor 1 (HIF-1) as an integral receptors target. In addition, KEGG analysis identified 11 other pathways besides the hub HIF-1 signaling route which may also be targeted in T2DM progression. In final molecular docking and dynamics simulation analysis, a significant binding affinity was confirmed for key compound-protein matrices. As such, the identified rooibos moieties could serve as putative drug candidates for T2DM control and therapy. This study shows rooibos constituents' interaction with T2DM-linked signaling pathways and target receptors and proposes vitexin, esculin and isovitexin as well as apigenin and kaempferol as respective pharmacologically active rooibos compounds for the modulation of EGFR and IGF1R in the HIF-1 signaling pathway to maintain normal homeostasis and function of the pancreas and pancreatic β-cells in diabetics.

Antiproliferative Activity of (-)-Isopulegol-based 1,3-Oxazine, 1,3-Thiazine and 2,4-Diaminopyrimidine Derivatives

A series of novel heterocyclic structures, namely 1,3‐oxazines, 1,3‐thiazines and 2,4‐diaminopyrimidines, were designed and synthesised. The bioassay tests demonstrated that, among these analogues, 2,4‐diaminopyridine derivatives showed significant antiproliferative activity against different human cancer cell lines (A2780, SiHa, HeLa, MCF‐7 and MDA‐MB‐231). Pyrimidines substituted with N²‐(p‐trifluoromethyl)aniline, in particular, displayed a potent inhibitory effect on the growth of cancer cells. Structure–activity relationships were also studied from the aspects of stereochemistry on the aminodiol moiety as well as exploring the effects of substituents on the pyrimidine scaffold.

Multi-Condition QSAR Model for the Virtual Design of Chemicals with Dual Pan-Antiviral and Anti-Cytokine Storm Profiles

Respiratory viruses are infectious agents, which can cause pandemics. Although nowadays the danger associated with respiratory viruses continues to be evidenced by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the virus responsible for the current COVID-19 pandemic, other viruses such as SARS-CoV-1, the influenza A and B viruses (IAV and IBV, respectively), and the respiratory syncytial virus (RSV) can lead to globally spread viral diseases. Also, from a biological point of view, most of these viruses can cause an organ-damaging hyperinflammatory response known as the cytokine storm (CS). Computational approaches constitute an essential component of modern drug development campaigns, and therefore, they have the potential to accelerate the discovery of chemicals able to simultaneously inhibit multiple molecular and nonmolecular targets. We report here the first multicondition model based on quantitative structure-activity relationships and an artificial neural network (mtc-QSAR-ANN) for the virtual design and prediction of molecules with dual pan-antiviral and anti-CS profiles. Our mtc-QSAR-ANN model exhibited an accuracy higher than 80%. By interpreting the different descriptors present in the mtc-QSAR-ANN model, we could retrieve several molecular fragments whose assembly led to new molecules with drug-like properties and predicted pan-antiviral and anti-CS activities.

Caspase-1 and Gasdermin D Afford the Optimal Targets with Distinct Switching Strategies in NLRP1b Inflammasome-Induced Cell Death

Inflammasomes are essential complexes of innate immune system, which form the first line of host defense against pathogens. Mounting evidence accumulates that inflammasome signaling is highly correlated with coronavirus disease 2019 (COVID-19). However, there remains a significant gap in our understanding of the regulatory mechanism of inflammasome signaling. Combining mathematical modeling with experimental analysis of NLRP1b inflammasome signaling, we found that only the expression levels of caspase-1 and GSDMD have the potential to individually switch cell death modes. Reduction of caspase-1 or GSDMD switches cell death from pyroptosis to apoptosis. Caspase-1 and GSDMD present different thresholds and exert distinct pathway choices in switching death modes. Pyroptosis switches to apoptosis with an extremely low threshold level of caspase-1, but with a high threshold of GSDMD. Caspase-1-impaired cells employ ASC-caspase-8-dependent pathway for apoptosis, while GSDMD-impaired cells primarily utilize caspase-1-dependent pathway. Additionally, caspase-1 and GSDMD can severally ignite the cooccurrence of pyroptosis and apoptosis. Landscape topography unravels that the cooccurrence is dramatically different in caspase-1- and GSDMD-impaired cells. Besides pyroptosis state and apoptosis state, a potential new “coexisting” state in single cells is proposed when GSDMD acts as the driving force of the landscape. The “seesaw model” is therefore proposed, which can well describe the death states that are controlled by caspase-1 or GSDMD in single cells. Our study sheds new light on NLRP1b inflammasome signaling and uncovers the switching mechanisms among various death modes, providing potential clues to guide the development of more rational control strategies for diseases.

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.

Cell death modes are specified by the crosstalk dynamics within pyroptotic and apoptotic signaling

The crosstalk between pyroptosis and apoptosis pathways plays crucial roles in homeostasis, cancer, and other pathologies. However, its molecular regulatory mechanisms for cell death decision-making remain to be elucidated. Based on the recent experimental studies, we developed a core regulatory network model of the crosstalk between pyroptosis and apoptosis pathways. Sensitivity analysis and bifurcation analysis were performed to assess the death mode switching of the network. Both the approaches determined that only the level of caspase-1 or gasdermin D (GSDMD) has the potential to individually change death modes. The decrease of caspase-1 or GSDMD switches cell death from pyroptosis to apoptosis. Seven biochemical reactions among the 21 reactions in total that are essential for determining cell death modes are identified by using sensitivity analysis. While with bifurcation analysis of state transitions, nine reactions are suggested to be able to efficiently switch death modes. Monostability, bistability, and tristability are observed under different conditions. We found that only the reaction that caspase-1 activation induced by stimuli can trigger tristability. Six and two of the nine reactions are identified to be able to induce bistability and monostability, respectively. Moreover, the concurrence of pyroptosis and apoptosis is observed not only within proper bistable ranges, but also within tristable ranges, implying two potentially distinct regulatory mechanisms. Taken together, this work sheds new light on the crosstalk between pyroptosis and apoptosis and uncovers the regulatory mechanisms of various stable state transitions, which play important roles for the development of potential control strategies for disease prevention and treatment.

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, when 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, 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 as 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.

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.

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.

Recent studies have shown that inhibition of the hSIRT2 enzyme provides favorable effects in neurodegenerative diseases such as Alzheimer's disease.

Transcription-translation error: In-silico investigation of the structural and functional impact of deleterious single nucleotide polymorphisms in GULP1 gene

Nonsynonymous single nucleotide polymorphisms (nsSNPs) are one of the most common forms of mutations known to disrupt the product of translation thereby altering the protein structure-function relationship. GULP1 (PTB domain-containing engulfment adaptor protein 1) is an evolutionarily conserved adaptor protein that has been associated with glycated hemoglobin (HbA1c) in Genome-Wide Association Studies (GWAS). In order to understand the role of GULP1 in the etiology of diabetes, it is important to study some functional nsSNPs present within the GULP1 protein. We, therefore, used a SNPinformatics approach to retrieve, classify, and determine the stability effect of some nsSNPs. Y27C, G142D, A144T, and Y149C were jointly predicted by the pathogenic-classifying tools to be disease-causing, however, only G142D, A144T, and Y149C had their structural architecture perturbed as predicted by I-MUTANT and MuPro. Interestingly, G142D and Y149C occur at positions 142 and 149 of GULP1 which coincidentally are found within the binding site of GULP1. Protein-Protein interaction analysis also revealed that GULP1 interacted with 10 proteins such as Cell division cycle 5-like protein (CDC5L), ADP-ribosylation factor 6 (ARF6), Arf-GAP with coiled-coil (ACAP1), and Multiple epidermal growth factor-like domains protein 10 (MEGF10), etc. Taken together, rs1357922096, rs1264999716, and rs128246649 could be used as genetic biomarkers for the diagnosis of diabetes. However, being a computational study, these nsSNPs require experimental validation to explore their metabolic involvement in the pathogenesis of diseases.

Computational and drug target analysis of functional single nucleotide polymorphisms associated with Haemoglobin Subunit Beta (HBB) gene

There is overwhelming evidence implicating Haemoglobin Subunit Beta (HBB) protein in the onset of beta thalassaemia. In this study for the first time, we used a combined SNP informatics and computer algorithms such as Neural network, Bayesian network, and Support Vector Machine to identify deleterious non-synonymous Single Nucleotide Polymorphisms (nsSNPs) present in the HBB gene. Our findings highlight three major mutation points (R31G, W38S, and Q128P) within the HBB gene sequence that have significant statistical and computational associations with the onset of beta thalassaemia. The dynamic simulation study revealed that R31G, W38S, and Q128P elicited high structural perturbation and instability, however, the wild type protein was considerably stable. Ten compounds with therapeutic potential against HBB were also predicted by structure-based virtual screening. Interestingly, the instability caused by the mutations was reversed upon binding to a ligand. This study has been able to predict potential deleterious mutants that can be further explored in the understanding of the pathological basis of beta thalassaemia and the design of tailored inhibitors.