Drug repurposing screen to identify inhibitors of the RNA polymerase (nsp12) and helicase (nsp13) from SARS-CoV-2 replication and transcription complex.

Coronaviruses contain one of the largest genomes among the RNA viruses, coding for 14-16 non-structural proteins (nsp) that are involved in proteolytic processing, genome replication and transcription, and four structural proteins that build the core of the mature virion. Due to conservation across coronaviruses, nsps form a group of promising drug targets as their inhibition directly affects viral replication and, therefore, progression of infection. A minimal but fully functional replication and transcription complex was shown to be formed by one RNA-dependent RNA polymerase (nsp12), one nsp7, two nsp8 accessory subunits, and two helicase (nsp13) enzymes. Our approach involved, targeting nsp12 and nsp13 to allow multiple starting point to interfere with virus infection progression. Here we report a combined in-vitro repurposing screening approach, identifying new and confirming reported SARS-CoV-2 nsp12 and nsp13 inhibitors.

"DompeKeys": a set of novel substructure-based descriptors for efficient chemical space mapping, development and structural interpretation of machine learning models, and indexing of large databases.

The conversion of chemical structures into computer-readable descriptors, able to capture key structural aspects, is of pivotal importance in the field of cheminformatics and computer-aided drug design. Molecular fingerprints represent a widely employed class of descriptors; however, their generation process is time-consuming for large databases and is susceptible to bias. Therefore, descriptors able to accurately detect predefined structural fragments and devoid of lengthy generation procedures would be highly desirable. To meet additional needs, such descriptors should also be interpretable by medicinal chemists, and suitable for indexing databases with trillions of compounds. To this end, we developed-as integral part of EXSCALATE, Dompé's end-to-end drug discovery platform-the DompeKeys (DK), a new substructure-based descriptor set, which encodes the chemical features that characterize compounds of pharmaceutical interest. DK represent an exhaustive collection of curated SMARTS strings, defining chemical features at different levels of complexity, from specific functional groups and structural patterns to simpler pharmacophoric points, corresponding to a network of hierarchically interconnected substructures. Because of their extended and hierarchical structure, DK can be used, with good performance, in different kinds of applications. In particular, we demonstrate how they are very well suited for effective mapping of chemical space, as well as substructure search and virtual screening. Notably, the incorporation of DK yields highly performing machine learning models for the prediction of both compounds' activity and metabolic reaction occurrence. The protocol to generate the DK is freely available at https://dompekeys.exscalate.eu and is fully integrated with the Molecular Anatomy protocol for the generation and analysis of hierarchically interconnected molecular scaffolds and frameworks, thus providing a comprehensive and flexible tool for drug design applications.

Diketo acid inhibitors of nsp13 of SARS-CoV-2 block viral replication.

For RNA viruses, RNA helicases have long been recognized to play critical roles during virus replication cycles, facilitating proper folding and replication of viral RNAs, therefore representing an ideal target for drug discovery. SARS-CoV-2 helicase, the non-structural protein 13 (nsp13) is a highly conserved protein among all known coronaviruses, and, at the moment, is one of the most explored viral targets to identify new possible antiviral agents. In the present study, we present six diketo acids (DKAs) as nsp13 inhibitors able to block both SARS-CoV-2 nsp13 enzymatic functions. Among them four compounds were able to inhibit viral replication in the low micromolar range, being active also on other human coronaviruses such as HCoV229E and MERS CoV. The experimental investigation of the binding mode revealed ATP-non-competitive kinetics of inhibition, not affected by substrate-displacement effect, suggesting an allosteric binding mode that was further supported by molecular modelling calculations predicting the binding into an allosteric conserved site located in the RecA2 domain.

MEDIATE - Molecular DockIng at homE: Turning collaborative simulations into therapeutic solutions.

INTRODUCTION: Collaborative computing has attracted great interest in the possibility of joining the efforts of researchers worldwide. Its relevance has further increased during the pandemic crisis since it allows for the strengthening of scientific collaborations while avoiding physical interactions. Thus, the E4C consortium presents the MEDIATE initiative which invited researchers to contribute via their virtual screening simulations that will be combined with AI-based consensus approaches to provide robust and method-independent predictions. The best compounds will be tested, and the biological results will be shared with the scientific community. AREAS COVERED: In this paper, the MEDIATE initiative is described. This shares compounds' libraries and protein structures prepared to perform standardized virtual screenings. Preliminary analyses are also reported which provide encouraging results emphasizing the MEDIATE initiative's capacity to identify active compounds. EXPERT OPINION: Structure-based virtual screening is well-suited for collaborative projects provided that the participating researchers work on the same input file. Until now, such a strategy was rarely pursued and most initiatives in the field were organized as challenges. The MEDIATE platform is focused on SARS-CoV-2 targets but can be seen as a prototype which can be utilized to perform collaborative virtual screening campaigns in any therapeutic field by sharing the appropriate input files.

DHFR Inhibitors Display a Pleiotropic Anti-Viral Activity against SARS-CoV-2: Insights into the Mechanisms of Action.

During the COVID-19 pandemic, drug repurposing represented an effective strategy to obtain quick answers to medical emergencies. Based on previous data on methotrexate (MTX), we evaluated the anti-viral activity of several DHFR inhibitors in two cell lines. We observed that this class of compounds showed a significant influence on the virus-induced cytopathic effect (CPE) partly attributed to the intrinsic anti-metabolic activity of these drugs, but also to a specific anti-viral function. To elucidate the molecular mechanisms, we took advantage of our EXSCALATE platform for in-silico molecular modelling and further validated the influence of these inhibitors on nsp13 and viral entry. Interestingly, pralatrexate and trimetrexate showed superior effects in counteracting the viral infection compared to other DHFR inhibitors. Our results indicate that their higher activity is due to their polypharmacological and pleiotropic profile. These compounds can thus potentially give a clinical advantage in the management of SARS-CoV-2 infection in patients already treated with this class of drugs.

Target Prediction by Multiple Virtual Screenings: Analyzing the SARS-CoV-2 Phenotypic Screening by the Docking Simulations Submitted to the MEDIATE Initiative.

Phenotypic screenings are usually combined with deconvolution techniques to characterize the mechanism of action for the retrieved hits. These studies can be supported by various computational analyses, although docking simulations are rarely employed. The present study aims to assess if multiple docking calculations can prove successful in target prediction. In detail, the docking simulations submitted to the MEDIATE initiative are utilized to predict the viral targets involved in the hits retrieved by a recently published cytopathic screening. Multiple docking results are combined by the EFO approach to develop target-specific consensus models. The combination of multiple docking simulations enhances the performances of the developed consensus models (average increases in EF1% value of 40% and 25% when combining three and two docking runs, respectively). These models are able to propose reliable targets for about half of the retrieved hits (31 out of 59). Thus, the study emphasizes that docking simulations might be effective in target identification and provide a convincing validation for the collaborative strategies that inspire the MEDIATE initiative. Disappointingly, cross-target and cross-program correlations suggest that common scoring functions are not specific enough for the simulated target.

Computational Insights into the Sequence-Activity Relationships of the NGF(1-14) Peptide by Molecular Dynamics Simulations.

The Nerve Growth Factor (NGF) belongs to the neurothrophins protein family involved in the survival of neurons in the nervous system. The interaction of NGF with its high-affinity receptor TrkA mediates different cellular pathways related to Alzheimer's disease, pain, ocular dysfunction, and cancer. Therefore, targeting NGF-TrkA interaction represents a valuable strategy for the development of new therapeutic agents. In recent years, experimental studies have revealed that peptides belonging to the N-terminal domain of NGF are able to partly mimic the biological activity of the whole protein paving the way towards the development of small peptides that can selectively target specific signaling pathways. Hence, understanding the molecular basis of the interaction between the N-terminal segment of NGF and TrkA is fundamental for the rational design of new peptides mimicking the NGF N-terminal domain. In this study, molecular dynamics simulation, binding free energy calculations and per-residue energy decomposition analysis were combined in order to explore the molecular recognition pattern between the experimentally active NGF(1-14) peptide and TrkA. The results highlighted the importance of His4, Arg9 and Glu11 as crucial residues for the stabilization of NGF(1-14)-TrkA interaction, thus suggesting useful insights for the structure-based design of new therapeutic peptides able to modulate NGF-TrkA interaction.

Cytopathic SARS-CoV-2 screening on VERO-E6 cells in a large-scale repurposing effort.

Worldwide, there are intensive efforts to identify repurposed drugs as potential therapies against SARS-CoV-2 infection and the associated COVID-19 disease. To date, the anti-inflammatory drug dexamethasone and (to a lesser extent) the RNA-polymerase inhibitor remdesivir have been shown to be effective in reducing mortality and patient time to recovery, respectively, in patients. Here, we report the results of a phenotypic screening campaign within an EU-funded project (H2020-EXSCALATE4COV) aimed at extending the repertoire of anti-COVID therapeutics through repurposing of available compounds and highlighting compounds with new mechanisms of action against viral infection. We screened 8702 molecules from different repurposing libraries, to reveal 110 compounds with an anti-cytopathic IC50 < 20 µM. From this group, 18 with a safety index greater than 2 are also marketed drugs, making them suitable for further study as potential therapies against COVID-19. Our result supports the idea that a systematic approach to repurposing is a valid strategy to accelerate the necessary drug discovery process.

Extensive Sampling of Molecular Dynamics Simulations to Identify Reliable Protein Structures for Optimized Virtual Screening Studies: The Case of the hTRPM8 Channel.

(1) Background: Virtual screening campaigns require target structures in which the pockets are properly arranged for binding. Without these, MD simulations can be used to relax the available target structures, optimizing the fine architecture of their binding sites. Among the generated frames, the best structures can be selected based on available experimental data. Without experimental templates, the MD trajectories can be filtered by energy-based criteria or sampled by systematic analyses. (2) Methods: A blind and methodical analysis was performed on the already reported MD run of the hTRPM8 tetrameric structures; a total of 50 frames underwent docking simulations by using a set of 1000 ligands including 20 known hTRPM8 modulators. Docking runs were performed by LiGen program and involved the frames as they are and after optimization by SCRWL4.0. For each frame, all four monomers were considered. Predictive models were developed by the EFO algorithm based on the sole primary LiGen scores. (3) Results: On average, the MD simulation progressively enhances the performance of the extracted frames, and the optimized structures perform better than the non-optimized frames (EF1% mean: 21.38 vs. 23.29). There is an overall correlation between performances and volumes of the explored pockets and the combination of the best performing frames allows to develop highly performing consensus models (EF1% = 49.83). (4) Conclusions: The systematic sampling of the entire MD run provides performances roughly comparable with those previously reached by using rationally selected frames. The proposed strategy appears to be helpful when the lack of experimental data does not allow an easy selection of the optimal structures for docking simulations. Overall, the reported docking results confirm the relevance of simulating all the monomers of an oligomer structure and emphasize the efficacy of the SCRWL4.0 method to optimize the protein structures for docking calculations.

Paclitaxel binds and activates C5aR1: A new potential therapeutic target for the prevention of chemotherapy-induced peripheral neuropathy and hypersensitivity reactions.

Chemotherapy-induced peripheral neuropathy (CIPN) and hypersensitivity reactions (HSRs) are among the most frequent and impairing side effects of the antineoplastic agent paclitaxel. Here, we demonstrated that paclitaxel can bind and activate complement component 5a receptor 1 (C5aR1) and that this binding is crucial in the etiology of paclitaxel-induced CIPN and anaphylaxis. Starting from our previous data demonstrating the role of interleukin (IL)-8 in paclitaxel-induced neuronal toxicity, we searched for proteins that activate IL-8 expression and, by using the Exscalate platform for molecular docking simulations, we predicted the high affinity of C5aR1 with paclitaxel. By in vitro studies, we confirmed the specific and competitive nature of the C5aR1-paclitaxel binding and found that it triggers intracellularly the NFkB/P38 pathway and c-Fos. In F11 neuronal cells and rat dorsal root ganglia, C5aR1 inhibition protected from paclitaxel-induced neuropathological effects, while in paclitaxel-treated mice, the absence (knock-out mice) or the inhibition of C5aR1 significantly ameliorated CIPN symptoms-in terms of cold and mechanical allodynia-and reduced the chronic pathological state in the paw. Finally, we found that C5aR1 inhibition can counteract paclitaxel-induced anaphylactic cytokine release in macrophages in vitro, as well as the onset of HSRs in mice. Altogether these data identified C5aR1 as a key mediator and a new potential pharmacological target for the prevention and treatment of CIPN and HSRs induced by paclitaxel.

Characterization of raloxifene as a potential pharmacological agent against SARS-CoV-2 and its variants.

The new coronavirus SARS-CoV-2 is the causative agent of the COVID-19 pandemic, which so far has caused over 6 million deaths in 2 years, despite new vaccines and antiviral medications. Drug repurposing, an approach for the potential application of existing pharmaceutical products to new therapeutic indications, could be an effective strategy to obtain quick answers to medical emergencies. Following a virtual screening campaign on the most relevant viral proteins, we identified the drug raloxifene, a known Selective Estrogen Receptor Modulator (SERM), as a new potential agent to treat mild-to-moderate COVID-19 patients. In this paper we report a comprehensive pharmacological characterization of raloxifene in relevant in vitro models of COVID-19, specifically in Vero E6 and Calu-3 cell lines infected with SARS-CoV-2. A large panel of the most common SARS-CoV-2 variants isolated in Europe, United Kingdom, Brazil, South Africa and India was tested to demonstrate the drug's ability in contrasting the viral cytopathic effect (CPE). Literature data support a beneficial effect by raloxifene against the viral infection due to its ability to interact with viral proteins and activate protective estrogen receptor-mediated mechanisms in the host cells. Mechanistic studies here reported confirm the significant affinity of raloxifene for the Spike protein, as predicted by in silico studies, and show that the drug treatment does not directly affect Spike/ACE2 interaction or viral internalization in infected cell lines. Interestingly, raloxifene can counteract Spike-mediated ADAM17 activation in human pulmonary cells, thus providing new insights on its mechanism of action. A clinical study in mild to moderate COVID-19 patients (NCT05172050) has been recently completed. Our contribution to evaluate raloxifene results on SARS-CoV-2 variants, and the interpretation of the mechanisms of action will be key elements to better understand the trial results, and to design new clinical studies aiming to evaluate the potential development of raloxifene in this indication.

Natural Compounds Inhibit SARS-CoV-2 nsp13 Unwinding and ATPase Enzyme Activities.

SARS-CoV-2 infection is still spreading worldwide, and new antiviral therapies are an urgent need to complement the approved vaccine preparations. SARS-CoV-2 nps13 helicase is a validated drug target participating in the viral replication complex and possessing two associated activities: RNA unwinding and 5'-triphosphatase. In the search of SARS-CoV-2 direct antiviral agents, we established biochemical assays for both SARS-CoV-2 nps13-associated enzyme activities and screened both in silico and in vitro a small in-house library of natural compounds. Myricetin, quercetin, kaempferol, and flavanone were found to inhibit the SARS-CoV-2 nps13 unwinding activity at nanomolar concentrations, while licoflavone C was shown to block both SARS-CoV-2 nps13 activities at micromolar concentrations. Mode of action studies showed that all compounds are nsp13 noncompetitive inhibitors versus ATP, while computational studies suggested that they can bind both nucleotide and 5'-RNA nsp13 binding sites, with licoflavone C showing a unique pattern of interaction with nsp13 amino acid residues. Overall, we report for the first time natural flavonoids as selective inhibitors of SARS-CoV-2 nps13 helicase with low micromolar activity.

"Molecular Anatomy": a new multi-dimensional hierarchical scaffold analysis tool.

The scaffold representation is widely employed to classify bioactive compounds on the basis of common core structures or correlate compound classes with specific biological activities. In this paper, we present a novel approach called "Molecular Anatomy" as a flexible and unbiased molecular scaffold-based metrics to cluster large set of compounds. We introduce a set of nine molecular representations at different abstraction levels, combined with fragmentation rules, to define a multi-dimensional network of hierarchically interconnected molecular frameworks. We demonstrate that the introduction of a flexible scaffold definition and multiple pruning rules is an effective method to identify relevant chemical moieties. This approach allows to cluster together active molecules belonging to different molecular classes, capturing most of the structure activity information, in particular when libraries containing a huge number of singletons are analyzed. We also propose a procedure to derive a network visualization that allows a full graphical representation of compounds dataset, permitting an efficient navigation in the scaffold's space and significantly contributing to perform high quality SAR analysis. The protocol is freely available as a web interface at https://ma.exscalate.eu .

A Blueprint for High Affinity SARS-CoV-2 Mpro Inhibitors from Activity-Based Compound Library Screening Guided by Analysis of Protein Dynamics.

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ∼30 000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ∼200 virtual screenings of compound libraries on selected protein structures, we redefine the protein's druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.

Combining Different Docking Engines and Consensus Strategies to Design and Validate Optimized Virtual Screening Protocols for the SARS-CoV-2 3CL Protease.

The 3CL-Protease appears to be a very promising medicinal target to develop anti-SARS-CoV-2 agents. The availability of resolved structures allows structure-based computational approaches to be carried out even though the lack of known inhibitors prevents a proper validation of the performed simulations. The innovative idea of the study is to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns. Docking simulations using four different programs (Fred, Glide, LiGen, and PLANTS) were performed investigating the role of both multiple binding modes (by binding space) and multiple isomers/states (by developing the corresponding isomeric space). The computed docking scores were used to develop consensus models, which allow an in-depth comparison of the resulting performances. On average, the reached performances revealed the different sensitivity to isomeric differences and multiple binding modes between the four docking engines. In detail, Glide and LiGen are the tools that best benefit from isomeric and binding space, respectively, while Fred is the most insensitive program. The obtained results emphasize the fruitful role of combining various docking tools to optimize the predictive performances. Taken together, the performed simulations allowed the rational development of highly performing virtual screening workflows, which could be further optimized by considering different 3CL-Pro structures and, more importantly, by including true SARS-CoV-2 3CL-Pro inhibitors (as learning set) when available.

Identification of Inhibitors of SARS-CoV-2 3CL-Pro Enzymatic Activity Using a Small Molecule in Vitro Repurposing Screen.

Compound repurposing is an important strategy for the identification of effective treatment options against SARS-CoV-2 infection and COVID-19 disease. In this regard, SARS-CoV-2 main protease (3CL-Pro), also termed M-Pro, is an attractive drug target as it plays a central role in viral replication by processing the viral polyproteins pp1a and pp1ab at multiple distinct cleavage sites. We here report the results of a repurposing program involving 8.7 K compounds containing marketed drugs, clinical and preclinical candidates, and small molecules regarded as safe in humans. We confirmed previously reported inhibitors of 3CL-Pro and have identified 62 additional compounds with IC50 values below 1 µM and profiled their selectivity toward chymotrypsin and 3CL-Pro from the Middle East respiratory syndrome virus. A subset of eight inhibitors showed anticytopathic effect in a Vero-E6 cell line, and the compounds thioguanosine and MG-132 were analyzed for their predicted binding characteristics to SARS-CoV-2 3CL-Pro. The X-ray crystal structure of the complex of myricetin and SARS-Cov-2 3CL-Pro was solved at a resolution of 1.77 Å, showing that myricetin is covalently bound to the catalytic Cys145 and therefore inhibiting its enzymatic activity.

Binding Mode Exploration of B1 Receptor Antagonists' by the Use of Molecular Dynamics and Docking Simulation-How Different Target Engagement Can Determine Different Biological Effects.

The kinin B1 receptor plays a critical role in the chronic phase of pain and inflammation. The development of B1 antagonists peaked in recent years but almost all promising molecules failed in clinical trials. Little is known about these molecules' mechanisms of action and additional information will be necessary to exploit the potential of the B1 receptor. With the aim of contributing to the available knowledge of the pharmacology of B1 receptors, we designed and characterized a novel class of allosteric non-peptidic inhibitors with peculiar binding characteristics. Here, we report the binding mode analysis and pharmacological characterization of a new allosteric B1 antagonist, DFL20656. We analyzed the binding of DFL20656 by single point mutagenesis and radioligand binding assays and we further characterized its pharmacology in terms of IC50, B1 receptor internalization and in vivo activity in comparison with different known B1 antagonists. We highlighted how different binding modes of DFL20656 and a Merck compound (compound 14) within the same molecular pocket can affect the biological and pharmacological properties of B1 inhibitors. DFL20656, by its peculiar binding mode, involving tight interactions with N114, efficiently induced B1 receptor internalization and evoked a long-lasting effect in an in vivo model of neuropathic pain. The pharmacological characterization of different B1 antagonists highlighted the effects of their binding modes on activity, receptor occupancy and internalization. Our results suggest that part of the failure of most B1 inhibitors could be ascribed to a lack of knowledge about target function and engagement.

Computational Studies of SARS-CoV-2 3CLpro: Insights from MD Simulations.

Given the enormous social and health impact of the pandemic triggered by severe acute respiratory syndrome 2 (SARS-CoV-2), the scientific community made a huge effort to provide an immediate response to the challenges posed by Coronavirus disease 2019 (COVID-19). One of the most important proteins of the virus is an enzyme, called 3CLpro or main protease, already identified as an important pharmacological target also in SARS and Middle East respiratory syndrome virus (MERS) viruses. This protein triggers the production of a whole series of enzymes necessary for the virus to carry out its replicating and infectious activities. Therefore, it is crucial to gain a deeper understanding of 3CLpro structure and function in order to effectively target this enzyme. All-atoms molecular dynamics (MD) simulations were performed to examine the different conformational behaviors of the monomeric and dimeric form of SARS-CoV-2 3CLpro apo structure, as revealed by microsecond time scale MD simulations. Our results also shed light on the conformational dynamics of the loop regions at the entry of the catalytic site. Studying, at atomic level, the characteristics of the active site and obtaining information on how the protein can interact with its substrates will allow the design of molecules able to block the enzymatic function crucial for the virus.

A Comprehensive Mapping of the Druggable Cavities within the SARS-CoV-2 Therapeutically Relevant Proteins by Combining Pocket and Docking Searches as Implemented in Pockets 2.0.

(1) Background: Virtual screening studies on the therapeutically relevant proteins of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) require a detailed characterization of their druggable binding sites, and, more generally, a convenient pocket mapping represents a key step for structure-based in silico studies; (2) Methods: Along with a careful literature search on SARS-CoV-2 protein targets, the study presents a novel strategy for pocket mapping based on the combination of pocket (as performed by the well-known FPocket tool) and docking searches (as performed by PLANTS or AutoDock/Vina engines); such an approach is implemented by the Pockets 2.0 plug-in for the VEGA ZZ suite of programs; (3) Results: The literature analysis allowed the identification of 16 promising binding cavities within the SARS-CoV-2 proteins and the here proposed approach was able to recognize them showing performances clearly better than those reached by the sole pocket detection; and (4) Conclusions: Even though the presented strategy should require more extended validations, this proved successful in precisely characterizing a set of SARS-CoV-2 druggable binding pockets including both orthosteric and allosteric sites, which are clearly amenable for virtual screening campaigns and drug repurposing studies. All results generated by the study and the Pockets 2.0 plug-in are available for download.

Combining Molecular Dynamics and Docking Simulations to Develop Targeted Protocols for Performing Optimized Virtual Screening Campaigns on The hTRPM8 Channel.

BACKGROUND: There is an increasing interest in TRPM8 ligands of medicinal interest, the rational design of which can be nowadays supported by structure-based in silico studies based on the recently resolved TRPM8 structures. Methods: The study involves the generation of a reliable hTRPM8 homology model, the reliability of which was assessed by a 1.0 µs MD simulation which was also used to generate multiple receptor conformations for the following structure-based virtual screening (VS) campaigns; docking simulations utilized different programs and involved all monomers of the selected frames; the so computed docking scores were combined by consensus approaches based on the EFO algorithm. Results: The obtained models revealed very satisfactory performances; LiGen™ provided the best results among the tested docking programs; the combination of docking results from the four monomers elicited a markedly beneficial effect on the computed consensus models. Conclusions: The generated hTRPM8 model appears to be amenable for successful structure-based VS studies; cross-talk modulating effects between interacting monomers on the binding sites can be accounted for by combining docking simulations as performed on all the monomers; this strategy can have general applicability for docking simulations involving quaternary protein structures with multiple identical binding pockets.