Introduction:The COVID-19 pandemic resulted in disastrous human and economic costs, mainly due to the initial lack of specific treatments. Complementary to immunotherapies, drug repurposing is possibly the best option to arrive at COVID-19 treatments in the short term. Areas covered: Repurposing prospects undergoing clinical trials or with some level of evidence emerging from clinical studies are overviewed. The authors discuss some possible intellectual property and commercial barriers to drug repurposing, and strategies to facilitate equitable access to incoming therapeutic solutions, highlighting the importance of collaborative drug discovery models. Based on a critical analysis of the available literature about in silico screens against SARS-CoV-2 main protease, the authors illustrate how frequently overconfident conclusions are being drawn in COVID-19-related literature. Expert opinion: Most of the current clinical trials on potential COVID-19 treatments are, in fact, drug repurposing examples. In October 2020, the FDA approved a repurposed antiviral, remdesivir, as the first treatment for COVID-19. Considering the high expectations invested in approaching therapeutic solutions, the scientific community must be careful not to raise unrealistic expectations. Today more than ever, the conclusions drawn in scientific reports have to be fully supported by the level of evidence, avoiding any sort of unfounded speculation.

The scientific community is working against the clock to arrive at therapeutic interventions to treat patients with COVID-19. Among the strategies for drug discovery, virtual screening approaches have the capacity to search potential hits within millions of chemical structures in days, with the appropriate computing infrastructure. In this article, we first analyzed the published research targeting the inhibition of the main protease (Mpro), one of the most studied targets of SARS-CoV-2, by docking-based methods. An alarming finding was the lack of an adequate validation of the docking protocols (i.e., pose prediction and virtual screening accuracy) before applying them in virtual screening campaigns. The performance of the docking protocols was tested at some level in 57.7% of the 168 investigations analyzed. However, we found only three examples of a complete retrospective analysis of the scoring functions to quantify the virtual screening accuracy of the methods. Moreover, only two publications reported some experimental evaluation of the proposed hits until preparing this manuscript. All of these findings led us to carry out a retrospective performance validation of three different docking protocols, through the analysis of their pose prediction and screening accuracy. Surprisingly, we found that even though all tested docking protocols have a good pose prediction, their screening accuracy is quite limited as they fail to correctly rank a test set of compounds. These results highlight the importance of conducting an adequate validation of the docking protocols before carrying out virtual screening campaigns, and to experimentally confirm the predictions made by the models before drawing bold conclusions. Finally, successful structure-based drug discovery investigations published during the redaction of this manuscript allow us to propose the inclusion of target flexibility and consensus scoring as alternatives to improve the accuracy of the methods.

The clustering of small molecules implies the organization of a group of chemical structures into smaller subgroups with similar features. Clustering has important applications to sample chemical datasets or libraries in a representative manner (e.g., to choose, from a virtual screening hit list, a chemically diverse subset of compounds to be submitted to experimental confirmation, or to split datasets into representative training and validation sets when implementing machine learning models). Most strategies for clustering molecules are based on molecular fingerprints and hierarchical clustering algorithms. Here, two open-source in-house methodologies for clustering of small molecules are presented: iterative Random subspace Principal Component Analysis clustering (iRaPCA), an iterative approach based on feature bagging, dimensionality reduction, and K-means optimization; and Silhouette Optimized Molecular Clustering (SOMoC), which combines molecular fingerprints with the Uniform Manifold Approximation and Projection (UMAP) and Gaussian Mixture Model algorithm (GMM). In a benchmarking exercise, the performance of both clustering methods has been examined across 29 datasets containing between 100 and 5000 small molecules, comparing these results with those given by two other well-known clustering methods, Ward and Butina. iRaPCA and SOMoC consistently showed the best performance across these 29 datasets, both in terms of within-cluster and between-cluster distances. Both iRaPCA and SOMoC have been implemented as free Web Apps and standalone applications, to allow their use to a wide audience within the scientific community.

The transient receptor potential vanilloid 1 (TRPV1) receptor is a nonselective cation channel, known to be involved in the regulation of many important physiological and pathological processes. In the last few years, it has been proposed as a promising target to develop novel anticonvulsant compounds. However, thermoregulatory effects associated with the channel inhibition have hampered the path for TRPV1 antagonists to become marketed drugs. In this regard, we conducted a structure-based virtual screening campaign to find potential TRPV1 modulators among approved drugs, which are known to be safe and thermally neutral. To this end, different docking models were developed and validated by assessing their pose and score prediction powers. Novobiocin, montelukast, and cinnarizine were selected from the screening as promising candidates for experimental testing and all of them exhibited nanomolar inhibitory activity. Moreover, the in vivo profiles showed promising results in at least one of the three models of seizures tested.

Trypanosoma cruzi carbonic anhydrase (TcCA) has recently emerged as an interesting target for the design of new compounds to treat Chagas disease. In this study we report the results of a structure-based virtual screening campaign to identify novel and selective TcCA inhibitors. The combination of properly validated computational methodologies such as comparative modelling, molecular dynamics and docking simulations allowed us to find high potency hits, with KI values in the nanomolar range. The compounds also showed trypanocidal effects against T. cruzi epimastigotes and trypomastigotes. All the candidates are selective for inhibiting TcCA over the human isoform CA II, which is encouraging in terms of possible therapeutic safety and efficacy.

No abstract

The human voltage-gated proton channel (hHv1) is a highly selective ion channel codified by the HVCN1 gene. It plays a fundamental role in several physiological processes such as innate and adaptive immunity, insulin secretion, and sperm capacitation. Moreover, in humans, a higher hHv1 expression/function has been reported in several types of cancer cells. Here we report a multitemplate homology model of the hHv1 channel, built and refined as a dimer in Rosetta. The model was then subjected to extensive Gaussian accelerated molecular dynamics (GaMD) for enhanced conformational sampling, and representative snapshots were extracted by clustering analysis. Combining different structure-and sequence-based methodologies, we predicted a putative ATP-binding site located on the intracellular portion of the channel. Furthermore, GaMD simulations of the ATP-bound dimeric hHv1 model showed that ATP interacts with a cluster of positively charged residues from the cytoplasmic N and C terminal segments. According to the in silico predictions, we found that 3 mM intracellular ATP significantly increases the H + current mediated by the hHv1 channel expressed in HEK293 cells and measured by the patch-clamp technique in an inside-out configuration (2.86 ± 0.63 fold over control at +40 mV). When ATP was added on the extracellular side, it was not able to activate the channel supporting the idea that the ATP-binding site resides in the intracellular face of the hHV1 channel. In a physiological and pathophysiological context, this ATP-mediated modulation could integrate the cell metabolic state with the H + efflux, especially in cells where hHv1 channels are relevant for pH regulation, such as pancreatic β-cells, immune cells, and cancer cells.

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.