Addressing docking pose selection with structure-based deep learning: Recent advances, challenges and opportunities.

Molecular docking is a widely used technique in drug discovery to predict the binding mode of a given ligand to its target. However, the identification of the near-native binding pose in docking experiments still represents a challenging task as the scoring functions currently employed by docking programs are parametrized to predict the binding affinity, and, therefore, they often fail to correctly identify the ligand native binding conformation. Selecting the correct binding mode is crucial to obtaining meaningful results and to conveniently optimizing new hit compounds. Deep learning (DL) algorithms have been an area of a growing interest in this sense for their capability to extract the relevant information directly from the protein-ligand structure. Our review aims to present the recent advances regarding the development of DL-based pose selection approaches, discussing limitations and possible future directions. Moreover, a comparison between the performances of some classical scoring functions and DL-based methods concerning their ability to select the correct binding mode is reported. In this regard, two novel DL-based pose selectors developed by us are presented.

Discovery of a New Class of 1-(4-Sulfamoylbenzoyl)piperidine-4-carboxamides as Human Carbonic Anhydrase Inhibitors.

A series of 1-(4-sulfamoylbenzoyl)piperidine-4-carboxamides deriving from substituted piperazines/benzylamines was designed, synthesized, and tested on human carbonic anhydrase (hCA). The inhibitory activity of the new sulfonamides was analyzed using acetazolamide (AAZ) as a standard inhibitor against hCA I, II, IX, and XII. Several sulfonamides showed both inhibitory activity at low nanomolar concentrations and selectivity against the cytosolic hCA II isoform, and the same trend was observed on the tumor-associated hCA IX and XII. The benzenesulfonamido carboxamides 11 and 15 were the most potent of the piperazino- and benzylamino-based series, respectively. Docking and molecular dynamics studies related the high selectivity of compound 11 toward the tumor-associated hCA isoforms to its capability to participate in favorable interactions within hCA IX and hCA XII active sites, whereas no such interactions were detected within both hCA I and hCA II isoforms.

Syntheses, reactivity, and biological applications of coumarins.

This comprehensive review, covering 2021-2023, explores the multifaceted chemical and pharmacological potential of coumarins, emphasizing their significance as versatile natural derivatives in medicinal chemistry. The synthesis and functionalization of coumarins have advanced with innovative strategies. This enabled the incorporation of diverse functional fragments or the construction of supplementary cyclic architectures, thereby the biological and physico-chemical properties of the compounds obtained were enhanced. The unique chemical structure of coumarine facilitates binding to various targets through hydrophobic interactions pi-stacking, hydrogen bonding, and dipole-dipole interactions. Therefore, this important scaffold exhibits promising applications in uncountable fields of medicinal chemistry (e.g., neurodegenerative diseases, cancer, inflammation).

Discovery of a septin-4 covalent binder with antimetastatic activity in a mouse model of melanoma.

Cancer spreading through metastatic processes is one of the major causes of tumour-related mortality. Metastasis is a complex phenomenon which involves multiple pathways ranging from cell metabolic alterations to changes in the biophysical phenotype of cells and tissues. In the search for new effective anti-metastatic agents, we modulated the chemical structure of the lead compound AA6, in order to find the structural determinants of activity, and to identify the cellular target responsible of the downstream anti-metastatic effects observed. New compounds synthesized were able to inhibit in vitro B16-F10 melanoma cell invasiveness, and one selected compound, CM365, showed in vivo anti-metastatic effects in a lung metastasis mouse model of melanoma. Septin-4 was identified as the most likely molecular target responsible for these effects. This study showed that CM365 is a promising molecule for metastasis prevention, remarkably effective alone or co-administered with drugs normally used in cancer therapy, such as paclitaxel.

Human carnosinases: A brief history, medicinal relevance, and in silico analyses.

Carnosine, an endogenous dipeptide, has been found to have a plethora of medicinal properties, such as antioxidant, antiageing, and chelating effects, but with one downside: a short half-life. Carnosinases and two hydrolytic enzymes, which remain enigmatic, are responsible for these features. Hence, here we emphasize why research is valuable for better understanding crucial concepts like ageing, neurodegradation, and cancerogenesis, given that inhibition of carnosinases might significantly prolong carnosine bioavailability and allow its further use in medicine. Herein, we explore the literature regarding carnosinases and present a short in silico analysis aimed at elucidating the possible recognition pattern between CN1 and its ligands.

Improved Antibacterial Activity of 1,3,4-Oxadiazole-Based Compounds That Restrict Staphylococcus aureus Growth Independent of LtaS Function.

The lipoteichoic acid (LTA) biosynthesis pathway has emerged as a promising antimicrobial therapeutic target. Previous studies identified the 1,3,4 oxadiazole compound 1771 as an LTA inhibitor with activity against Gram-positive pathogens. We have succeeded in making six 1771 derivatives and, through subsequent hit validation, identified the incorporation of a pentafluorosulfanyl substituent as central in enhancing activity. Our newly described derivative, compound 13, showed a 16- to 32-fold increase in activity compared to 1771 when tested against a cohort of multidrug-resistant Staphylococcus aureus strains while simultaneously exhibiting an improved toxicity profile against mammalian cells. Molecular techniques were employed in which the assumed target, lipoteichoic acid synthase (LtaS), was both deleted and overexpressed. Neither deletion nor overexpression of LtaS altered 1771 or compound 13 susceptibility; however, overexpression of LtaS increased the MIC of Congo red, a previously identified LtaS inhibitor. These data were further supported by comparing the docking poses of 1771 and derivatives in the LtaS active site, which indicated the possibility of an additional target(s). Finally, we show that both 1771 and compound 13 have activity that is independent of LtaS, extending to cover Gram-negative species if the outer membrane is first permeabilized, challenging the classification that these compounds are strict LtaS inhibitors.

Heterocyclic Compounds as Synthetic Tyrosinase Inhibitors: Recent Advances.

Tyrosinase is a copper-containing enzyme which is widely distributed in nature (e.g., bacteria, mammals, fungi) and involved in two consecutive steps of melanin biosynthesis. In humans, an excessive production of melanin can determine hyperpigmentation disorders as well as neurodegenerative processes in Parkinson's disease. The development of molecules able to inhibit the high activity of the enzyme remain a current topic in medicinal chemistry, because the inhibitors reported so far present several side effects. Heterocycle-bearing molecules are largely diffuse in this sense. Due to their importance as biologically active compounds, we decided to report a comprehensive review of synthetic tyrosinase inhibitors possessing heterocyclic moieties reported within the last five years. For the reader's convenience, we classified them as inhibitors of mushroom tyrosinase (Agaricus bisporus) and human tyrosinase.

The VEGA web service: multipurpose online tools for molecular modelling and docking analyses.

The paper presents the VEGA Online web service, which includes a set of freely available tools deriving from the development of the VEGA suite of programs. In detail, the paper is focused on two tools: the VEGA Web Edition (WE) and the Score tool. The former is a versatile file format converter including relevant features for 2D/3D conversion, for surface mapping and for editing/preparing input files. The Score application allows rescoring docking poses and in particular includes the MLP Interactions Scores (MLPInS) for describing hydrophobic interactions. To the best of our knowledge, this web service is the only available resource by which one can calculate both the virtual log P of a given input molecule according to the MLP approach plus the corresponding MLP surface.

Ensemble of structure and ligand-based classification models for hERG liability profiling.

Drug-induced cardiotoxicity represents one of the most critical safety concerns in the early stages of drug development. The blockade of the human ether-à-go-go-related potassium channel (hERG) is the most frequent cause of cardiotoxicity, as it is associated to long QT syndrome which can lead to fatal arrhythmias. Therefore, assessing hERG liability of new drugs candidates is crucial to avoid undesired cardiotoxic effects. In this scenario, computational approaches have emerged as useful tools for the development of predictive models able to identify potential hERG blockers. In the last years, several efforts have been addressed to generate ligand-based (LB) models due to the lack of experimental structural information about hERG channel. However, these methods rely on the structural features of the molecules used to generate the model and often fail in correctly predicting new chemical scaffolds. Recently, the 3D structure of hERG channel has been experimentally solved enabling the use of structure-based (SB) strategies which may overcome the limitations of the LB approaches. In this study, we compared the performances achieved by both LB and SB classifiers for hERG-related cardiotoxicity developed by using Random Forest algorithm and employing a training set containing 12789 hERG binders. The SB models were trained on a set of scoring functions computed by docking and rescoring calculations, while the LB classifiers were built on a set of physicochemical descriptors and fingerprints. Furthermore, models combining the LB and SB features were developed as well. All the generated models were internally validated by ten-fold cross-validation on the TS and further verified on an external test set. The former revealed that the best performance was achieved by the LB model, while the model combining the LB and the SB attributes displayed the best results when applied on the external test set highlighting the usefulness of the integration of LB and SB features in correctly predicting unseen molecules. Overall, our predictive models showed satisfactory performances providing new useful tools to filter out potential cardiotoxic drug candidates in the early phase of drug discovery.

Interaction of Aggregated Cationic Porphyrins with Human Serum Albumin.

The interaction of an equilibrium mixture of monomeric and aggregated cationic trans-5,15-bis(N-methylpyridinium-4-yl)-10,15-bis-diphenylporphine (t-H2Pagg) chloride salt with human serum albumin (HSA) has been investigated through UV/Vis absorption, fluorescence emission, circular dichroism and resonant light scattering techniques. The spectroscopic evidence reveals that both the monomeric t-H2Pagg and its aggregates bind instantaneously to HSA, leading to the formation of a tight adduct in which the porphyrin is encapsulated within the protein scaffold (S430) and to clusters of aggregated porphyrins in electrostatic interaction with the charged biomolecules. These latter species eventually interconvert into the final S430 species following pseudo-first-order kinetics. Molecular docking simulations have been performed to get some insights into the nature of the final adduct. Analogously to hemin bound to HSA, the obtained model supports favorable interactions of the porphyrin in the same 1B subdomain of the protein. Hydrophobic and van der Waals energy terms are the main contributions to the calculated ΔGbind value of -117.24 kcal/mol.

An integrated metabolomic and proteomic approach for the identification of covalent inhibitors of the main protease (Mpro) of SARS-COV-2 from crude natural extracts.

Mpro represents one of the most promising drug targets for SARS-Cov-2, as it plays a crucial role in the maturation of viral polyproteins into functional proteins. HTS methods are currently used to screen Mpro inhibitors, and rely on searching chemical databases and compound libraries, meaning that they only consider previously structurally clarified and isolated molecules. A great advancement in the hit identification strategy would be to set-up an approach aimed at exploring un-deconvoluted mixtures of compounds such as plant extracts. Hence, the aim of the present study is to set-up an analytical platform able to fish-out bioactive molecules from complex natural matrices even where there is no knowledge on the constituents. The proposed approach begins with a metabolomic step aimed at annotating the MW of the matrix constituents. A further metabolomic step is based on identifying those natural electrophilic compounds able to form a Michael adduct with thiols, a peculiar chemical feature of many Mpro inhibitors that covalently bind the catalytic Cys145 in the active site, thus stabilizing the complex. A final step consists of incubating recombinant Mpro with natural extracts and identifying compounds adducted to the residues within the Mpro active site by bottom-up proteomic analysis (nano-LC-HRMS). Data analysis is based on two complementary strategies: (i) a targeted search applied by setting the adducted moieties identified as Michael acceptors of Cys as variable modifications; (ii) an untargeted approach aimed at identifying the whole range of adducted peptides containing Cys145 on the basis of the characteristic b and y fragment ions independent of the adduct. The method was set-up and then successfully tested to fish-out bioactive compounds from the crude extract of Scutellaria baicalensis, a Chinese plant containing the catechol-like flavonoid baicalin and its corresponding aglycone baicalein which are well-established inhibitors of Mpro. Molecular dynamics (MD) simulations were carried out in order to explore the binding mode of baicalin and baicalein, within the SARS-CoV-2 Mpro active site, allowing a better understanding of the role of the nucleophilic residues (i.e. His41, Cys145, His163 and His164) in the protein-ligand recognition process.

Synthesis, molecular docking and antibacterial activity of an oxadiazole-based lipoteichoic acid inhibitor and its metabolites.

Amongst drug resistant Gram-positive bacteria, Staphylococcus aureus is a pathogen of great concern as it is the leading cause of life-threatening nosocomial and community acquired infections which are often associated with implanted medical devices. The biosynthesis of lipotheicoic acid (LTA) by S. aureus has been recognized as a promising antibacterial target, owing its critical role in the growth and survival of Gram-positive bacteria. Here we report for the first time the chemical synthesis and characterisation of an oxadiazole based compound (1771), previously described as an inhibitor of LTA biosynthesis by targeting Lta synthase enzyme (LtaS). To investigate its controversial mode of action, we also performed molecular docking studies, which indicated that 1771 behaves as a competitive inhibitor against LtaS. We also synthesised and evaluated the antimicrobial activity of 1771 metabolites which we have identified from its decomposition in mouse serum, proving that the biological activity was caused by intact 1771.

Ligand-Based Discovery of a Small Molecule as Inhibitor of α-Synuclein Amyloid Formation.

α-Synuclein (α-Syn) aggregates are implicated in Parkinson's disease (PD), so inhibitors of α-Syn aggregation have been intensively explored. It has been demonstrated that small molecules might be able to reduce α-Syn aggregation in fibrils, thus exerting neuroprotective effects in models of PD. To expand our knowledge about the structural requirements for blocking the recognition process into the oligomeric assembly of α-Syn aggregates, we performed a ligand-based virtual screening procedure using two well-known α-Syn aggregation inhibitors, SynuClean-D and ZPD-2, as query compounds. A collection of thirty-four compounds bearing distinct chemical functionalities and mutual chemical features were studied in a Th-T fluorescence test, thus identifying 5-(2,6-dinitro-4-(trifluoromethyl)benzyl)-1-methyl-1H-tetrazole (named MeSC-04) as a potent α-Syn amyloid formation inhibitor that demonstrated similar behavior when compared to SynuClean-D in the thioflavin-T-monitored kinetic assays, with both molecules reducing the number and size of amyloid fibrils, as evidenced by electron microscopy. Molecular modeling studies suggested the binding mode of MeSC-04 through the identification of putative druggable pockets on α-syn fibrils and a subsequent consensus docking methodology. Overall, this work could furnish new insights in the development of α-Syn amyloid inhibitors from synthetic sources.

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.

Discovery of Neuroprotective Agents Based on a 5-(4-Pyridinyl)-1,2,4-triazole Scaffold.

Parkinson's disease (PD) is characterized by the death of dopaminergic neurons. The common histopathological hallmark in PD patients is the formation of intracellular proteinaceous accumulations. The main constituent of these inclusions is alpha-synuclein (α-syn), an intrinsically disordered protein that in pathological conditions creates amyloid aggregates that lead to neurotoxicity and neurodegeneration. The main goal of our study was to optimize our previously identified α-syn aggregation inhibitors of 5-(4-pyridinyl)-1,2,4-triazole chemotype in terms of in vivo efficacy. Our efforts resulted in the identification of ethyl 2-((4-amino-5-(pyridin-4-yl)-4H-1,2,4-triazol-3-yl)thio)acetate (15), which displayed the ability to prevent 1-methyl-4-phenyl-1,2,3,6-tetrahydropiridine-induced bradykinesia as well as to affect the levels of PD markers after the administration of the same neurotoxin. In addition to the in vivo evaluation, for the 5-(4-pyridinyl)-1,2,4-triazole-based compounds, we measured the prevention of the fibrillization process using light scattering and a ThT binding assay; these compounds have been shown to slightly reduce the α-syn aggregation.

Investigation on Hydrazonobenzenesulfonamides as Human Carbonic Anhydrase I, II, IX and XII Inhibitors.

A small series of hydrazonobenzenesulfonamides was designed, synthesized and studied for their human carbonic anhydrase (hCA) inhibitory activity. The synthesized compounds were evaluated against hCA I, II, IX and XII isoforms using acetazolamide (AAZ) as the standard inhibitor. Various hydrazonosulfonamide derivatives showed inhibitory activity at low nanomolar levels with selectivity against the cytosolic hCA II isoform, as well as the transmembrane, tumor-associated enzymes hCA IX and XII. The most potent and selective hydrazones 8, 9, 10, 11, 19 and 24 were docked into isoforms I, II, IX and XII to better understand their activity and selectivity for the different CA isoforms.

Exploration of the 2,3-dihydroisoindole pharmacophore for inhibition of the influenza virus PA endonuclease.

Seasonal influenza A and B viruses represent a global concern. Antiviral drugs are crucial to treat severe influenza in high-risk patients and prevent virus spread in case of a pandemic. The emergence of viruses showing drug resistance, in particular for the recently licensed polymerase inhibitor baloxavir marboxil, drives the need for developing alternative antivirals. The endonuclease activity residing in the N-terminal domain of the polymerase acidic protein (PAN) is crucial for viral RNA synthesis and a validated target for drug design. Its function can be impaired by molecules bearing a metal-binding pharmacophore (MBP) able to coordinate the two divalent metal ions in the active site. In the present work, the 2,3-dihydro-6,7-dihydroxy-1H-isoindol-1-one scaffold is explored for the inhibition of influenza virus PA endonuclease. The structure-activity relationship was analysed by modifying the substituents on the lipophilic moiety linked to the MBP. The new compounds exhibited nanomolar inhibitory activity in a FRET-based enzymatic assay, and a few compounds (15-17, 21) offered inhibition in the micromolar range, in a cell-based influenza virus polymerase assay. When investigated against a panel of PA-mutant forms, compound 17 was shown to retain full activity against the baloxavir-resistant I38T mutant. This was corroborated by docking studies providing insight into the binding mode of this novel class of PA inhibitors.

Evaluation of 4-(4-Fluorobenzyl)piperazin-1-yl]-Based Compounds as Competitive Tyrosinase Inhibitors Endowed with Antimelanogenic Effects.

There is a considerable attention for the development of inhibitors of tyrosinase (TYR) as therapeutic strategy for the treatment of hyperpigmentation disorders in humans. Continuing in our efforts to identify TYR inhibitors, we describe the design, synthesis and pharmacophore exploration of new small molecules structurally characterized by the presence of the 4-fluorobenzylpiperazine moiety as key pharmacophoric feature for the inhibition of TYR from Agaricus bisporus (AbTYR). Our investigations resulted in the discovery of the competitive inhibitor [4-(4-fluorobenzyl)piperazin-1-yl]-(3-chloro-2-nitro-phenyl)methanone 26 (IC50 =0.18 µM) that proved to be ∼100-fold more active than reference compound kojic acid (IC50 =17.76 µM). Notably, compound 26 exerted antimelanogenic effect on B16F10 cells in absence of cytotoxicity. Docking analysis suggested its binding mode into AbTYR and into modelled human TYR.

Exploring Molecular Contacts of MUC1 at CIN85 Binding Interface to Address Future Drug Design Efforts.

The modulation of protein-protein interactions (PPIs) by small molecules represents a valuable strategy for pharmacological intervention in several human diseases. In this context, computer-aided drug discovery techniques offer useful resources to predict the network of interactions governing the recognition process between protein partners, thus furnishing relevant information for the design of novel PPI modulators. In this work, we focused our attention on the MUC1-CIN85 complex as a crucial PPI controlling cancer progression and metastasis. MUC1 is a transmembrane glycoprotein whose extracellular domain contains a variable number of tandem repeats (VNTRs) regions that are highly glycosylated in normal cells and under-glycosylated in cancer. The hypo-glycosylation fosters the exposure of the backbone to new interactions with other proteins, such as CIN85, that alter the intracellular signalling in tumour cells. Herein, different computational approaches were combined to investigate the molecular recognition pattern of MUC1-CIN85 PPI thus unveiling new structural information useful for the design of MUC1-CIN85 PPI inhibitors as potential anti-metastatic agents.

In Silico Strategy for Targeting the mTOR Kinase at Rapamycin Binding Site by Small Molecules.

Computer aided drug-design methods proved to be powerful tools for the identification of new therapeutic agents. We employed a structure-based workflow to identify new inhibitors targeting mTOR kinase at rapamycin binding site. By combining molecular dynamics (MD) simulation and pharmacophore modelling, a simplified structure-based pharmacophore hypothesis was built starting from the FKBP12-rapamycin-FRB ternary complex retrieved from RCSB Protein Data Bank (PDB code 1FAP). Then, the obtained model was used as filter to screen the ZINC biogenic compounds library, containing molecules derived from natural sources or natural-inspired compounds. The resulting hits were clustered according to their similarity; moreover, compounds showing the highest pharmacophore fit-score were chosen from each cluster. The selected molecules were subjected to docking studies to clarify their putative binding mode. The binding free energy of the obtained complexes was calculated by MM/GBSA method and the hits characterized by the lowest ΔGbind values were identified as potential mTOR inhibitors. Furthermore, the stability of the resulting complexes was studied by means of MD simulation which revealed that the selected compounds were able to form a stable ternary complex with FKBP12 and FRB domain, thus underlining their potential ability to inhibit mTOR with a rapamycin-like mechanism.