Bespoke library docking for 5-HT2A receptor agonists with antidepressant activity.

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

Experimental and Computational Study of Aryl-thiosemicarbazones Inhibiting Cruzain Reveals Reversible Inhibition and a Stepwise Mechanism.

Trypanosoma cruzi is a parasite that infects about 6-7 million people worldwide, mostly in Latin America, causing Chagas disease. Cruzain, the main cysteine protease of T. cruzi, is a validated target for developing drug candidates for Chagas disease. Thiosemicarbazones are one of the most relevant warheads used in covalent inhibitors targeting cruzain. Despite its relevance, the mechanism of inhibition of cruzain by thiosemicarbazones is unknown. Here, we combined experiments and simulations to unveil the covalent inhibition mechanism of cruzain by a thiosemicarbazone-based inhibitor (compound 1). Additionally, we studied a semicarbazone (compound 2), which is structurally similar to compound 1 but does not inhibit cruzain. Assays confirmed the reversibility of inhibition by compound 1 and suggested a two-step mechanism of inhibition. The Ki was estimated to be 36.3 µM and Ki* to be 11.5 µM, suggesting the pre-covalent complex to be relevant for inhibition. Molecular dynamics simulations of compounds 1 and 2 with cruzain were used to propose putative binding modes for the ligands. One-dimensional (1D) quantum mechanics/molecular mechanics (QM/MM) potential of mean force (PMF) and gas-phase energies showed that the attack of Cys25-S- on the C═S or C═O bond yields a more stable intermediate than the attack on the C═N bond of the thiosemicarbazone/semicarbazone. Two-dimensional (2D) QM/MM PMF revealed a putative reaction mechanism for compound 1, involving the proton transfer to the ligand, followed by the Cys25-S- attack at C═S. The ΔG and energy barrier were estimated to be -1.4 and 11.7 kcal/mol, respectively. Overall, our results shed light on the inhibition mechanism of cruzain by thiosemicarbazones.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations.

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC50 = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Synthesis and biological evaluation of potential inhibitors of the cysteine proteases cruzain and rhodesain designed by molecular simplification.

Analogues of 8-chloro-N-(3-morpholinopropyl)-5H-pyrimido[5,4-b]indol-4-amine 1, a known cruzain inhibitor, were synthesized using a molecular simplification strategy. Five series of analogues were obtained: indole, pyrimidine, quinoline, aniline and pyrrole derivatives. The activity of the compounds was evaluated against the enzymes cruzain and rhodesain as well as against Trypanosoma cruzi amastigote and trypomastigote forms. The 4-aminoquinoline derivatives showed promising activity against both enzymes, with IC50 values ranging from 15 to 125µM. These derivatives were selective inhibitors for the parasitic proteases, being unable to inhibit mammalian cathepsins B and S. The most active compound against cruzain (compound 5a; IC50=15µM) is considerably more synthetically accessible than 1, while retaining its ligand efficiency. As observed for the original lead, compound 5a was shown to be a competitive enzyme inhibitor. In addition, it was also active against T. cruzi (IC50=67.7µM). Interestingly, the pyrimidine derivative 4b, although inactive in enzymatic assays, was highly active against T. cruzi (IC50=3.1µM) with remarkable selectivity index (SI=128) compared to uninfected fibroblasts. Both 5a and 4b exhibit drug-like physicochemical properties and are predicted to have a favorable ADME profile, therefore having great potential as candidates for lead optimization in the search for new drugs to treat Chagas disease.

Screening the Pathogen Box to Discover and Characterize New Cruzain and TbrCatL Inhibitors.

Chagas disease and Human African Trypanosomiasis, caused by Trypanosoma cruzi and T. brucei, respectively, pose relevant health challenges throughout the world, placing 65 to 70 million people at risk each. Given the limited efficacy and severe side effects associated with current chemotherapy, new drugs are urgently needed for both diseases. Here, we report the screening of the Pathogen Box collection against cruzain and TbrCatL, validated targets for Chagas disease and Human African Trypanosomiasis, respectively. Enzymatic assays were applied to screen 400 compounds, validate hits, determine IC50 values and, when possible, mechanisms of inhibition. In this case, 12 initial hits were obtained and ten were prioritized for follow-up. IC50 values were obtained for six of them (hit rate = 1.5%) and ranged from 0.46 ± 0.03 to 27 ± 3 µM. MMV687246 was found to be a mixed inhibitor of cruzain (Ki = 57 ± 6 µM) while MMV688179 was found to be a competitive inhibitor of cruzain with a nanomolar potency (Ki = 165 ± 63 nM). A putative binding mode for MMV688179 was obtained by docking. The six hits discovered against cruzain and TbrCatL are of great interest for further optimization by the medicinal chemistry community.

Synthesis, Design, and Structure-Activity Relationship of a Benzenesulfonylpiperazine Series against Trypanosoma cruzi.

Chagas disease is a neglected tropical disease, endemic in Latin America and caused by the protozoan parasite Trypanosoma cruzi. Available treatments show low cure efficacy during the chronic phase of the disease and cause a series of side effects, reinforcing the need to develop new drugs against Chagas disease. In this work, we describe the optimization of a trypanocidal hit compound recently reported in phenotypic high-throughput screening studies against Trypanosoma cruzi. A hit-to-lead process was initiated and a structure-activity relationship against Trypanosoma cruzi was obtained after the synthesis and biological evaluation of 22 new benzenesulfonylpiperazine derivatives. From this structure-activity relationship study, we identified three compounds with a promising predicted ADMET profile and potency comparable to the reference drug benznidazole, which are candidates for further development towards therapies for Chagas disease.

Fragment Binding to the Nsp3 Macrodomain of SARS-CoV-2 Identified Through Crystallographic Screening and Computational Docking.

The SARS-CoV-2 macrodomain (Mac1) within the non-structural protein 3 (Nsp3) counteracts host-mediated antiviral ADP-ribosylation signalling. This enzyme is a promising antiviral target because catalytic mutations render viruses non-pathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of diverse fragment libraries resulted in 214 unique macrodomain-binding fragments, out of 2,683 screened. An additional 60 molecules were selected from docking over 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several crystallographic and docking fragment hits were validated for solution binding using three biophysical techniques (DSF, HTRF, ITC). Overall, the 234 fragment structures presented explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.