Structure of human TRPM8 channel.

TRPM8 is a non-selective cation channel permeable to both monovalent and divalent cations that is activated by multiple factors, such as temperature, voltage, pressure, and changes in osmolality. It is a therapeutic target for anticancer drug development, and its modulators can be utilized for several pathological conditions. Here, we present a cryo-electron microscopy structure of a human TRPM8 channel in the closed state that was solved at 2.7 Å resolution. Our structure comprises the most complete model of the N-terminal pre-melastatin homology region. We also visualized several lipids that are bound by the protein and modeled how the human channel interacts with icilin. Analyses of pore helices in available TRPM structures showed that all these structures can be grouped into different closed, desensitized and open state conformations based on the register of the pore helix S6 which positions particular amino acid residues at the channel constriction.

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.

Conformational response to ligand binding of TMPRSS2, a protease involved in SARS-CoV-2 infection: Insights through computational modeling.

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS-CoV-2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors.

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.

The SARS-CoV-2 spike protein binds and modulates estrogen receptors.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein binds angiotensin-converting enzyme 2 as its primary infection mechanism. Interactions between S and endogenous proteins occur after infection but are not well understood. We profiled binding of S against >9000 human proteins and found an interaction between S and human estrogen receptor α (ERα). Using bioinformatics, supercomputing, and experimental assays, we identified a highly conserved and functional nuclear receptor coregulator (NRC) LXD-like motif on the S2 subunit. In cultured cells, S DNA transfection increased ERα cytoplasmic accumulation, and S treatment induced ER-dependent biological effects. Non-invasive imaging in SARS-CoV-2-infected hamsters localized lung pathology with increased ERα lung levels. Postmortem lung experiments from infected hamsters and humans confirmed an increase in cytoplasmic ERα and its colocalization with S in alveolar macrophages. These findings describe the discovery of a S-ERα interaction, imply a role for S as an NRC, and advance knowledge of SARS-CoV-2 biology and coronavirus disease 2019 pathology.

The Spike Mutants Website: A Worldwide Used Resource against SARS-CoV-2.

A large number of SARS-CoV-2 mutations in a short period of time has driven scientific research related to vaccines, new drugs, and antibodies to combat the new variants of the virus. Herein, we present a web portal containing the structural information, the tridimensional coordinates, and the molecular dynamics trajectories of the SARS-CoV-2 spike protein and its main variants. The Spike Mutants website can serve as a rapid online tool for investigating the impact of novel mutations on virus fitness. Taking into account the high variability of SARS-CoV-2, this application can help the scientific community when prioritizing molecules for experimental assays, thus, accelerating the identification of promising drug candidates for COVID-19 treatment. Below we describe the main features of the platform and illustrate the possible applications for speeding up the drug discovery process and hypothesize new effective strategies to overcome the recurrent mutations in SARS-CoV-2 genome.

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.

Structural Evolution of Delta (B.1.617.2) and Omicron (BA.1) Spike Glycoproteins.

The vast amount of epidemiologic and genomic data that were gathered as a global response to the COVID-19 pandemic that was caused by SARS-CoV-2 offer a unique opportunity to shed light on the structural evolution of coronaviruses and in particular on the spike (S) glycoprotein, which mediates virus entry into the host cell by binding to the human ACE2 receptor. Herein, we carry out an investigation into the dynamic properties of the S glycoprotein, focusing on the much more transmissible Delta and Omicron variants. Notwithstanding the great number of mutations that have accumulated, particularly in the Omicron S glycoprotein, our data clearly showed the conservation of some structural and dynamic elements, such as the global motion of the receptor binding domain (RBD). However, our studies also revealed structural and dynamic alterations that were concentrated in the aa 627-635 region, on a small region of the receptor binding motif (aa 483-485), and the so-called "fusion-peptide proximal region". In particular, these last two S regions are known to be involved in the human receptor ACE2 recognition and membrane fusion. Our structural evidence, therefore, is likely involved in the observed different transmissibility of these S mutants. Finally, we highlighted the role of glycans in the increased RBD flexibility of the monomer in the up conformation of Omicron.

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.

The SARS-CoV-2 spike protein binds and modulates estrogen receptors.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein binds angiotensin-converting enzyme 2 (ACE2) at the cell surface, which constitutes the primary mechanism driving SARS-CoV-2 infection. Molecular interactions between the transduced S and endogenous proteins likely occur post-infection, but such interactions are not well understood. We used an unbiased primary screen to profile the binding of full-length S against >9,000 human proteins and found significant S-host protein interactions, including one between S and human estrogen receptor alpha (ERα). After confirming this interaction in a secondary assay, we used bioinformatics, supercomputing, and experimental assays to identify a highly conserved and functional nuclear receptor coregulator (NRC) LXD-like motif on the S2 subunit and an S-ERα binding mode. In cultured cells, S DNA transfection increased ERα cytoplasmic accumulation, and S treatment induced ER-dependent biological effects and ACE2 expression. Noninvasive multimodal PET/CT imaging in SARS-CoV-2-infected hamsters using [ 18 F]fluoroestradiol (FES) localized lung pathology with increased ERα lung levels. Postmortem experiments in lung tissues from SARS-CoV-2-infected hamsters and humans confirmed an increase in cytoplasmic ERα expression and its colocalization with S protein in alveolar macrophages. These findings describe the discovery and characterization of a novel S-ERα interaction, imply a role for S as an NRC, and are poised to advance knowledge of SARS-CoV-2 biology, COVID-19 pathology, and mechanisms of sex differences in the pathology of infectious disease.

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.

Repurposing the estrogen receptor modulator raloxifene to treat SARS-CoV-2 infection.

The ongoing coronavirus disease 2019 (COVID-19) pandemic caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) necessitates strategies to identify prophylactic and therapeutic drug candidates to enter rapid clinical development. This is particularly true, given the uncertainty about the endurance of the immune memory induced by both previous infections or vaccines, and given the fact that the eradication of SARS-CoV-2 might be challenging to reach, given the attack rate of the virus, which would require unusually high protection by a vaccine. Here, we show how raloxifene, a selective estrogen receptor modulator with anti-inflammatory and antiviral properties, emerges as an attractive candidate entering clinical trials to test its efficacy in early-stage treatment COVID-19 patients.

SCoV2-MD: a database for the dynamics of the SARS-CoV-2 proteome and variant impact predictions.

SCoV2-MD (www.scov2-md.org) is a new online resource that systematically organizes atomistic simulations of the SARS-CoV-2 proteome. The database includes simulations produced by leading groups using molecular dynamics (MD) methods to investigate the structure-dynamics-function relationships of viral proteins. SCoV2-MD cross-references the molecular data with the pandemic evolution by tracking all available variants sequenced during the pandemic and deposited in the GISAID resource. SCoV2-MD enables the interactive analysis of the deposited trajectories through a web interface, which enables users to search by viral protein, isolate, phylogenetic attributes, or specific point mutation. Each mutation can then be analyzed interactively combining static (e.g. a variety of amino acid substitution penalties) and dynamic (time-dependent data derived from the dynamics of the local geometry) scores. Dynamic scores can be computed on the basis of nine non-covalent interaction types, including steric properties, solvent accessibility, hydrogen bonding, and other types of chemical interactions. Where available, experimental data such as antibody escape and change in binding affinities from deep mutational scanning experiments are also made available. All metrics can be combined to build predefined or custom scores to interrogate the impact of evolving variants on protein structure and function.

"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 .

Altered Local Interactions and Long-Range Communications in UK Variant (B.1.1.7) Spike Glycoprotein.

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the "UK variant", in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type; this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic "hotspots" and epitopes targeted by neutralizing antibodies.