MM-GB/SA rescoring of docking poses

Histone lactylation inhibits RARγ expression in macrophages to promote colorectal tumorigenesis through activation of TRAF6-IL-6-STAT3 signaling

Macrophages are phenotypically and functionally diverse in the tumor microenvironment (TME). However, how to remodel macrophages with a protumor phenotype and how to manipulate them for therapeutic purposes remain to be explored. Here, we show that in the TME, RARγ is downregulated in macrophages, and its expression correlates with poor prognosis in patients with colorectal cancer (CRC). In macrophages, RARγ interacts with tumor necrosis factor receptor-associated factor 6 (TRAF6), which prevents TRAF6 oligomerization and autoubiquitination, leading to inhibition of nuclear factor κB signaling. However, tumor-derived lactate fuels H3K18 lactylation to prohibit RARγ gene transcription in macrophages, consequently enhancing interleukin-6 (IL-6) levels in the TME and endowing macrophages with tumor-promoting functions via activation of signal transducer and activator of transcription 3 (STAT3) signaling in CRC cells. We identified that nordihydroguaiaretic acid (NDGA) exerts effective antitumor action by directly binding to RARγ to inhibit TRAF6-IL-6-STAT3 signaling. This study unravels lactate-driven macrophage function remodeling by inhibition of RARγ expression and highlights NDGA as a candidate compound for treating CRC.

Thiadiazolidinone (TDZD) Analogs Inhibit Aggregation-Mediated Pathology in Diverse Neurodegeneration Models, and Extend C. elegans Life- and Healthspan

Chronic, low-grade inflammation has been implicated in aging and age-dependent conditions, including Alzheimer's disease, cardiomyopathy, and cancer. One of the age-associated processes underlying chronic inflammation is protein aggregation, which is implicated in neuroinflammation and a broad spectrum of neurodegenerative diseases such as Alzheimer's, Huntington's, and Parkinson's diseases. We screened a panel of bioactive thiadiazolidinones (TDZDs) from our in-house library for rescue of protein aggregation in human-cell and C. elegans models of neurodegeneration. Among the tested TDZD analogs, PNR886 and PNR962 were most effective, significantly reducing both the number and intensity of Alzheimer-like tau and amyloid aggregates in human cell-culture models of pathogenic aggregation. A C. elegans strain expressing human Aβ1-42 in muscle, leading to AD-like amyloidopathy, developed fewer and smaller aggregates after PNR886 or PNR962 treatment. Moreover, age-progressive paralysis was reduced 90% by PNR886 and 75% by PNR962, and "healthspan" (the median duration of spontaneous motility) was extended 29% and 62%, respectively. These TDZD analogs also extended wild-type C. elegans lifespan by 15-30% (p < 0.001), placing them among the most effective life-extension drugs. Because the lead drug in this family, TDZD-8, inhibits GSK3β, we used molecular-dynamic tools to assess whether these analogs may also target GSK3β. In silico modeling predicted that PNR886 or PNR962 would bind to the same allosteric pocket of inactive GSK3β as TDZD-8, employing the same pharmacophore but attaching with greater avidity. PNR886 and PNR962 are thus compelling candidate drugs for treatment of tau- and amyloid-associated neurodegenerative diseases such as AD, potentially also reducing all-cause mortality.

Combined strategies in structure-based virtual screening

The identification and optimization of lead compounds are inalienable components in drug design and discovery pipelines. As a powerful computational approach for the identification of hits with novel structural scaffolds, structure-based virtual screening (SBVS) has exhibited a remarkably increasing influence in the early stages of drug discovery. During the past decade, a variety of techniques and algorithms have been proposed and tested with different purposes in the scope of SBVS. Although SBVS has been a common and proven technology, it still shows some challenges and problems that are needed to be addressed, where the negative influence regardless of protein flexibility and the inaccurate prediction of binding affinity are the two major challenges. Here, focusing on these difficulties, we summarize a series of combined strategies or workflows developed by our group and others. Furthermore, several representative successful applications from recent publications are also discussed to demonstrate the effectiveness of the combined SBVS strategies in drug discovery campaigns.

Computational elucidation of the binding mechanisms of curcumin analogues as bacterial RecA inhibitors

Antimicrobial resistance (AMR) presents as a serious threat to global public health, which urgently demands action to develop alternative antimicrobial strategies with minimized selective pressure. The bacterial SOS response regulator RecA has emerged as a promising target in the exploration of new classes of antibiotic adjuvants, as RecA has been implicated in bacterial mutagenesis and thus AMR development through its critical roles in error-prone DNA repair. The natural product curcumin has been reported to be an effective RecA inhibitor in several Gram-negative bacteria, but details on the underlying mechanisms are wanting. In order to bridge the gap in how curcumin operates as a RecA inhibitor, we used computational approaches to model interactions between RecA protein and curcumin analogues. We first identified potential binding sites on E. coli RecA protein and classified them into four major binding pockets based on biological literature and computational findings from multiple in silico calculations. In docking analysis, curcumin-thalidomide hybrids were predicted to be superior binders of RecA compared with bis-(arylmethylidene)acetone curcumin analogues, which was further confirmed by MMGBSA calculations. Overall, this work provides mechanistic insights into bacterial RecA protein as a target for curcumin-like compounds and offers a theoretical basis for rational design and development of future antibiotic adjuvants.

Molecular mechanism of MBX2319 inhibition of Escherichia coli AcrB multidrug efflux pump and comparison with other inhibitors

Efflux pumps of the resistance nodulation division (RND) superfamily, such as AcrB, make a major contribution to multidrug resistance in Gram-negative bacteria. The development of inhibitors of the RND pumps would improve the efficacy of current and next-generation antibiotics. To date, however, only one inhibitor has been cocrystallized with AcrB. Thus, in silico structure-based analysis is essential for elucidating the interaction between other inhibitors and the efflux pumps. In this work, we used computer docking and molecular dynamics simulations to study the interaction between AcrB and the compound MBX2319, a novel pyranopyridine efflux pump inhibitor with potent activity against RND efflux pumps of Enterobacteriaceae species, as well as other known inhibitors (D13-9001, 1-[1-naphthylmethyl]-piperazine, and phenylalanylarginine-β-naphthylamide) and the binding of doxorubicin to the efflux-defective F610A variant of AcrB. We also analyzed the binding of a substrate, minocycline, for comparison. Our results show that MBX2319 binds very tightly to the lower part of the distal pocket in the B protomer of AcrB, strongly interacting with the phenylalanines lining the hydrophobic trap, where the hydrophobic portion of D13-9001 was found to bind by X-ray crystallography. Additionally, MBX2319 binds to AcrB in a manner that is similar to the way in which doxorubicin binds to the F610A variant of AcrB. In contrast, 1-(1-naphthylmethyl)-piperazine and phenylalanylarginine-β-naphthylamide appear to bind to somewhat different areas of the distal pocket in the B protomer of AcrB than does MBX2319. However, all inhibitors (except D13-9001) appear to distort the structure of the distal pocket, impairing the proper binding of substrates.

Comprehending renin inhibitor's binding affinity using structure-based approaches

The performance of several structure-based design (SBD) approaches in predicting the binding affinity of diverse small molecule inhibitors co-crystallized to human renin was assessed to ascertain the modeling tool and method of choice required when dealing with structure-based lead optimization projects. Most of the SBD approaches investigated here were able to provide qualitative guidance, but quantitative accuracy as well as decisive discrimination between [in]actives is still not within reach. Such an outcome suggests that the current methods need improvement to capture the overall physics of the binding phenomenon for consistent applications in a lead optimization setting.

Molecular modeling, quantum polarized ligand docking and structure-based 3D-QSAR analysis of the imidazole series as dual AT(1) and ET(A) receptor antagonists

AIM

Both endothelin ETA receptor antagonists and angiotensin AT1 receptor antagonists lower blood pressure in hypertensive patients. A dual AT1 and ETA receptor antagonist may be more efficacious antihypertensive drug. In this study we identified the mode and mechanism of binding of imidazole series of compounds as dual AT1 and ETA receptor antagonists.

METHODS

Molecular modeling approach combining quantum-polarized ligand docking (QPLD), MM/GBSA free-energy calculation and 3D-QSAR analysis was used to evaluate 24 compounds as dual AT1 and ETA receptor antagonists and to reveal their binding modes and structural basis of the inhibitory activity. Pharmacophore-based virtual screening and docking studies were performed to identify more potent dual antagonists.

RESULTS

3D-QSAR models of the imidazole compounds were developed from the conformer generated by QPLD, and the resulting models showed a good correlation between the predicted and experimental activity. The visualization of the 3D-QSAR model in the context of the compounds under study revealed the details of the structure-activity relationship: substitution of methoxymethyl and cyclooctanone might increase the activity against AT1 receptor, while substitution of cyclohexone and trimethylpyrrolidinone was important for the activity against ETA receptor; addition of a trimethylpyrrolidinone to compound 9 significantly reduced its activity against AT1 receptor but significantly increased its activity against ETA receptor, which was likely due to the larger size and higher intensities of the H-bond donor and acceptor regions in the active site of ETA receptor. Pharmacophore-based virtual screening followed by subsequent Glide SP, XP, QPLD and MM/GBSA calculation identified 5 potential lead compounds that might act as dual AT1 and ETA receptor antagonists.

CONCLUSION

This study may provide some insights into the development of novel potent dual ETA and AT1 receptor antagonists. As a result, five compounds are found to be the best dual antagonists against AT1R and ETA receptors.