NMR characterization and ligand binding site of the stem-loop 2 motif from the Delta variant of SARS-CoV-2.

The stem-loop 2 motif (s2m) in SARS-CoV-2 (SCoV-2) is located in the 3'-UTR. Although s2m has been reported to display characteristics of a mobile genomic element that might lead to an evolutionary advantage, its function has remained unknown. The secondary structure of the original SCoV-2 RNA sequence (Wuhan-Hu-1) was determined by NMR in late 2020, delineating the base-pairing pattern and revealing substantial differences in secondary structure compared to SARS-CoV-1 (SCoV-1). The existence of a single G29742-A29756 mismatch in the upper stem of s2m leads to its destabilization and impedes a complete NMR analysis. With Delta, a variant of concern has evolved with one mutation compared to the original sequence that replaces G29742 by U29742. We show here that this mutation results in a more defined structure at ambient temperature accompanied by a rise in melting temperature. Consequently, we were able to identify >90% of the relevant NMR resonances using a combination of selective RNA labeling and filtered 2D NOESY as well as 4D NMR experiments. We present a comprehensive NMR analysis of the secondary structure, (sub)nanosecond dynamics, and ribose conformation of s2m Delta based on heteronuclear 13C NOE and T 1 measurements and ribose carbon chemical shift-derived canonical coordinates. We further show that the G29742U mutation in Delta has no influence on the druggability of s2m compared to the Wuhan-Hu-1 sequence. With the assignment at hand, we identify the flexible regions of s2m as the primary site for small molecule binding.

Exploring Binding Pockets in the Conformational States of the SARS-CoV-2 Spike Trimers for the Screening of Allosteric Inhibitors Using Molecular Simulations and Ensemble-Based Ligand Docking.

Understanding mechanisms of allosteric regulation remains elusive for the SARS-CoV-2 spike protein, despite the increasing interest and effort in discovering allosteric inhibitors of the viral activity and interactions with the host receptor ACE2. The challenges of discovering allosteric modulators of the SARS-CoV-2 spike proteins are associated with the diversity of cryptic allosteric sites and complex molecular mechanisms that can be employed by allosteric ligands, including the alteration of the conformational equilibrium of spike protein and preferential stabilization of specific functional states. In the current study, we combine conformational dynamics analysis of distinct forms of the full-length spike protein trimers and machine-learning-based binding pocket detection with the ensemble-based ligand docking and binding free energy analysis to characterize the potential allosteric binding sites and determine structural and energetic determinants of allosteric inhibition for a series of experimentally validated allosteric molecules. The results demonstrate a good agreement between computational and experimental binding affinities, providing support to the predicted binding modes and suggesting key interactions formed by the allosteric ligands to elicit the experimentally observed inhibition. We establish structural and energetic determinants of allosteric binding for the experimentally known allosteric molecules, indicating a potential mechanism of allosteric modulation by targeting the hinges of the inter-protomer movements and blocking conformational changes between the closed and open spike trimer forms. The results of this study demonstrate that combining ensemble-based ligand docking with conformational states of spike protein and rigorous binding energy analysis enables robust characterization of the ligand binding modes, the identification of allosteric binding hotspots, and the prediction of binding affinities for validated allosteric modulators, which is consistent with the experimental data. This study suggested that the conformational adaptability of the protein allosteric sites and the diversity of ligand bound conformations are both in play to enable efficient targeting of allosteric binding sites and interfere with the conformational changes.

Multivalent calix[4]arene-based mannosylated dendrons as new FimH ligands and inhibitors.

We report the synthesis and biological characterization of a novel class of multivalent glycoconjugates as hit compounds for the design of new antiadhesive therapies against urogenital tract infections (UTIs) caused by uropathogenic E. coli strains (UPEC). The first step of UTIs is the molecular recognition of high mannose N-glycan expressed on the surface of urothelial cells by the bacterial lectin FimH, allowing the pathogen adhesion required for mammalian cell invasion. The inhibition of FimH-mediated interactions is thus a validated strategy for the treatment of UTIs. To this purpose, we designed and synthesized d-mannose multivalent dendrons supported on a calixarene core introducing a significant structural change from a previously described family of dendrimers bearing the same dendrons units on a flexible pentaerythritol scaffold core. The new molecular architecture increased the inhibitory potency against FimH-mediated adhesion processes by about 16 times, as assessed by yeast agglutination assay. Moreover, the direct molecular interaction of the new compounds with FimH protein was assessed by on-cell NMR experiments acquired in the presence of UPEC cells.

Inhibition of Neutral Sphingomyelinase 2 by Novel Small Molecule Inhibitors Results in Decreased Release of Extracellular Vesicles by Vascular Smooth Muscle Cells and Attenuated Calcification.

Vascular calcification (VC) is an important contributor and prognostic factor in the pathogenesis of cardiovascular diseases. VC is an active process mediated by the release of extracellular vesicles by vascular smooth muscle cells (VSMCs), and the enzyme neutral sphingomyelinase 2 (nSMase2 or SMPD3) plays a key role. Upon activation, the enzyme catalyzes the hydrolysis of sphingomyelin, thereby generating ceramide and phosphocholine. This conversion mediates the release of exosomes, a type of extracellular vesicles (EVs), which ultimately forms the nidus for VC. nSMase2 therefore represents a drug target, the inhibition of which is thought to prevent or halt VC progression. In search of novel druglike small molecule inhibitors of nSMase2, we have used virtual ligand screening to identify potential ligands. From an in-silico collection of 48,6844 small druglike molecules, we selected 996 compounds after application of an in-house multi-step procedure combining different filtering and docking procedures. Selected compounds were functionally tested in vitro; from this, we identified 52 individual hit molecules that inhibited nSMase2 activity by more than 20% at a concentration of 150 µM. Further analysis showed that five compounds presented with IC50s lower than 2 µM. Of these, compounds ID 5728450 and ID 4011505 decreased human primary VSMC EV release and calcification in vitro. The hit molecules identified here represent new classes of nSMase2 inhibitors that may be developed into lead molecules for the therapeutic or prophylactic treatment of VC.

Identification of antiparasitic drug targets using a multi-omics workflow in the acanthocephalan model.

BACKGROUND: With the expansion of animal production, parasitic helminths are gaining increasing economic importance. However, application of several established deworming agents can harm treated hosts and environment due to their low specificity. Furthermore, the number of parasite strains showing resistance is growing, while hardly any new anthelminthics are being developed. Here, we present a bioinformatics workflow designed to reduce the time and cost in the development of new strategies against parasites. The workflow includes quantitative transcriptomics and proteomics, 3D structure modeling, binding site prediction, and virtual ligand screening. Its use is demonstrated for Acanthocephala (thorny-headed worms) which are an emerging pest in fish aquaculture. We included three acanthocephalans (Pomphorhynchus laevis, Neoechinorhynchus agilis, Neoechinorhynchus buttnerae) from four fish species (common barbel, European eel, thinlip mullet, tambaqui). RESULTS: The workflow led to eleven highly specific candidate targets in acanthocephalans. The candidate targets showed constant and elevated transcript abundances across definitive and accidental hosts, suggestive of constitutive expression and functional importance. Hence, the impairment of the corresponding proteins should enable specific and effective killing of acanthocephalans. Candidate targets were also highly abundant in the acanthocephalan body wall, through which these gutless parasites take up nutrients. Thus, the candidate targets are likely to be accessible to compounds that are orally administered to fish. Virtual ligand screening led to ten compounds, of which five appeared to be especially promising according to ADMET, GHS, and RO5 criteria: tadalafil, pranazepide, piketoprofen, heliomycin, and the nematicide derquantel. CONCLUSIONS: The combination of genomics, transcriptomics, and proteomics led to a broadly applicable procedure for the cost- and time-saving identification of candidate target proteins in parasites. The ligands predicted to bind can now be further evaluated for their suitability in the control of acanthocephalans. The workflow has been deposited at the Galaxy workflow server under the URL tinyurl.com/yx72rda7 .

Structural Characterization of KOR Inactive and Active States for 3D Pharmacology and Drug Discovery.

The structure of the human kappa opioid receptor (KOR) in complex with the long-acting antagonist JDTic was solved crystallographically in 2012 and, along with structures of other opioid receptors, revolutionized our understanding of opioid system function and pharmacology. More recently, active state KOR structure was also determined, giving important insights into activation mechanisms of the receptor. In this review, we will discuss how the understanding of atomistic structures of KOR established a key platform for deciphering details of subtype and functional selectivity of KOR-targeting ligands and for discovery of new chemical probes with potentially beneficial pharmacological profiles.

A simple and sensitive detection of the binding ligands by using the receptor aggregation and NMR spectroscopy: a test case of the maltose binding protein.

Protein-ligand interaction is one of the highlights of molecular recognition. The most popular application of this type of interaction is drug development which requires a high throughput screening of a ligand that binds to the target protein. Our goal was to find a binding ligand with a simple detection, and once this type of ligand was found, other methods could then be used to measure the detailed kinetic or thermodynamic parameters. We started with the idea that the ligand NMR signal would disappear if it was bound to the non-tumbling mass. In order to create the non-tumbling mass, we tried the aggregates of a target protein, which was fused to the elastin-like polypeptide. We chose the maltose binding proteinas a test case, and we tried it with several sugars, which included maltose, glucose, sucrose, lactose, galactose, maltotriose, and ß-cyclodextrin. The maltose signal in the H-1 NMR spectrum disappeared completely as hoped around the protein to ligand ratio of 1:3 at 298 K where the proteins aggregated. The protein signals also disappeared upon aggregation except for the fast-moving part, which resulted in a cleaner background than the monomeric form. Since we only needed to look for a disappearing signal amongst those from the mixture, it should be useful in high throughput screening. Other types of sugars except for the maltotriose and ß-cyclodextrin, which are siblings of the maltose, did not seem to bind at all. We believe that our system would be especially more effective when dealing with a smaller target protein, so both the protein and the bound ligand would lose their signals only when the aggregates formed. We hope that our proposed method would contribute to accelerating the development of the potent drug candidates by simultaneously identifying several binders directly from a mixture.

On-cell saturation transfer difference NMR for the identification of FimH ligands and inhibitors.

We describe the development of an on-cell NMR method for the rapid screening of FimH ligands and the structural identification of ligand binding epitopes. FimH is a mannose-binding bacterial adhesin expressed at the apical end of type 1 pili of uropathogenic bacterial strains and responsible for their d-mannose sensitive adhesion to host mammalian epithelial cells. Because of these properties, FimH is a key virulence factor and an attractive therapeutic target for urinary tract infection. We prepared synthetic d-mannose decorated dendrimers, we tested their ability to prevent the FimH-mediated yeast agglutination, and thus we used the compounds showing the best inhibitory activity as models of FimH multivalent ligands to set up our NMR methodology. Our experimental protocol, based on on-cell STD NMR techniques, is a suitable tool for the screening and the epitope mapping of FimH ligands aimed at the development of new antiadhesive and diagnostic tools against urinary tract infection pathogens. Notably, the study is carried out in a physiological environment, i.e. at the surface of living pathogen cells expressing FimH.

Insights into Dynamic Mechanism of Ligand Binding to Peroxisome Proliferator-Activated Receptor γ toward Potential Pharmacological Applications.

Peroxisome proliferator-activated receptor γ (PPARγ) is a member of the nuclear receptor superfamily, which regulates the transcription of a variety of genes involved in lipid and glucose metabolism, inflammation, and cell proliferation. These functions correlate with the onset of type-2 diabetes, obesity, and immune disorders, which makes PPARγ a promising target for drug development. The majority of PPARγ functions are regulated by binding of small molecule ligands, which cause conformational changes of PPARγ followed by coregulator recruitment. The ligand-binding domain (LBD) of PPARγ contains a large Y-shaped cavity that can be occupied by various classes of compounds such as full agonists, partial agonists, natural lipids, and in some cases, a combination of multiple molecules. Several crystal structure studies have revealed the binding modes of these compounds in the LBD and insight into the resulting conformational changes. Notably, the apo form of the PPARγ LBD contains a highly mobile region that can be stabilized by ligand binding. Furthermore, recent biophysical investigations have shed light on the dynamic mechanism of how ligands induce conformational changes in PPARγ and result in functional output. This information may be useful for the design of new and repurposed structures of ligands that serve a different function from original compounds and more potent pharmacological effects with less undesirable clinical outcomes. This review provides an overview of the peculiar characteristics of the PPARγ LBD by examining a series of structural studies focused on the dynamic mechanism of binding and the potential applications of strategies for ligand screening and chemical labeling.

Machine Learning Assisted Approach for Finding Novel High Activity Agonists of Human Ectopic Olfactory Receptors.

Olfactory receptors (ORs) constitute the largest superfamily of G protein-coupled receptors (GPCRs). ORs are involved in sensing odorants as well as in other ectopic roles in non-nasal tissues. Matching of an enormous number of the olfactory stimulation repertoire to its counterpart OR through machine learning (ML) will enable understanding of olfactory system, receptor characterization, and exploitation of their therapeutic potential. In the current study, we have selected two broadly tuned ectopic human OR proteins, OR1A1 and OR2W1, for expanding their known chemical space by using molecular descriptors. We present a scheme for selecting the optimal features required to train an ML-based model, based on which we selected the random forest (RF) as the best performer. High activity agonist prediction involved screening five databases comprising ~23 M compounds, using the trained RF classifier. To evaluate the effectiveness of the machine learning based virtual screening and check receptor binding site compatibility, we used docking of the top target ligands to carefully develop receptor model structures. Finally, experimental validation of selected compounds with significant docking scores through in vitro assays revealed two high activity novel agonists for OR1A1 and one for OR2W1.

A comprehensive in vitro fluorescence anisotropy assay system for screening ligands of the jasmonate COI1-JAZ co-receptor in plants.

The phytohormone (+)-7-iso-jasmonoyl-l-isoleucine regulates many developmental and stress responses in plants and induces protein-protein interactions between COI1, the F-box component of E3 ubiquitin ligase, and jasmonate ZIM domain (JAZ) repressors. These interactions cause JAZ degradation and activate jasmonate (JA), leading to plant defense responses, growth inhibition, and senescence. Thirteen JAZ subtypes are encoded in the Arabidopsis thaliana genome, but a detailed understanding of the physiological functions of these JAZ subtypes remains unclear, partially because of the genetic redundancy of JAZ genes. One strategy to elucidate the complex JA signaling pathways is to develop a reliable and comprehensive binding assay system of the ligands with all combinations of the co-receptors. Herein, we report the development of a fluorescence anisotropy-based in vitro binding assay system to screen for the ligands of the COI1-JAZ co-receptors. Our assay enabled the first quantitative analysis of the affinity values and JAZ-subtype selectivity of various endogenous JA derivatives, such as coronatine, jasmonic acid, and 12-hydroxyjasmonoyl-l-isoleucine. Because of its high signal-to-noise ratio and convenient mix-and-read assay system, our screening approach can be used in plate reader-based assays of both agonists and antagonists of COI1-JAZ co-receptors.

Structure based discovery of novel hexokinase 2 inhibitors.

Hexokinase 2 (HK2) is over-expressed in most of human cancers and has been proved to be a promising target for cancer therapy. In this study, based on the structure of HK2, we screened over 6 millions of compounds to obtain the lead. A total of 26 (E)-N'-(2,3,4-trihydroxybenzylidene) arylhydrazide derivatives were then designed, synthesized, and evaluated for their HK2 enzyme activity and IC50 values against two cancer cell lines. Most of the 26 target compounds showed excellently in vitro activity. Among them, compound 3j showed the strongest inhibitory effects on HK2 enzyme activity with an IC50 of 0.53 ± 0.13 µM and exhibited the most potent growth inhibition against SW480 cells with an IC50 of 7.13 ± 1.12 µM, which deserves further studies.

Discovery of New Immune Checkpoints: Family Grows Up.

The first generation of immune checkpoint inhibitors (ICIs) including anti-CTLA-4 and anti-PD-1/anti-PD-L1 has achieved profound and great success. Till 2019 Q1, there are nine ICIs landing the oncology market: Ipilimumab (anti-CTLA-4, Bristol-Myers Squibb), Nivolumab (anti-PD-1, Bristol-Myers Squibb), Pembrolizumab (anti-PD-1, Merck), Atezolizumab (anti-PD-L1, Roche/Genentech), Durvalumab (anti-PD-L1, Astra Zeneca), Tremelimumab (anti-CTLA-4, Astra Zeneca), Cemiplimab (anti-PD-1, Sanofi/Regeneron), Toripalimab (anti-PD-1, Junshi), and Sintilimab (anti-PD-1, Innovent), which have covered the majority of hematologic and solid malignancies' indication. Beyond the considerable benefits for the patients, frustrated boundary still exists: limited response rate in monotherapy in late-stage population, poor effectiveness in neoplasms with immune desert and immune excluded types, and immune-related toxicities, some are life-threatened and with higher incidence in I-O combination regiment. Moreover, clinicians observed some cases switching to progression after achieving partial or complete response, indicating treatment failure or drug resistance. So people begin looking for the next generation of immune checkpoint members.

Protein X-ray Crystallography and Drug Discovery.

With the advent of structural biology in the drug discovery process, medicinal chemists gained the opportunity to use detailed structural information in order to progress screening hits into leads or drug candidates. X-ray crystallography has proven to be an invaluable tool in this respect, as it is able to provide exquisitely comprehensive structural information about the interaction of a ligand with a pharmacological target. As fragment-based drug discovery emerged in the recent years, X-ray crystallography has also become a powerful screening technology, able to provide structural information on complexes involving low-molecular weight compounds, despite weak binding affinities. Given the low numbers of compounds needed in a fragment library, compared to the hundreds of thousand usually present in drug-like compound libraries, it now becomes feasible to screen a whole fragment library using X-ray crystallography, providing a wealth of structural details that will fuel the fragment to drug process. Here, we review theoretical and practical aspects as well as the pros and cons of using X-ray crystallography in the drug discovery process.

Chemical space of Escherichia coli dihydrofolate reductase inhibitors: New approaches for discovering novel drugs for old bugs.

Escherichia coli Dihydrofolate reductase is an important enzyme that is essential for the survival of the Gram-negative microorganism. Inhibitors designed against this enzyme have demonstrated application as antibiotics. However, either because of poor bioavailability of the small-molecules resulting from their inability to cross the double membrane in Gram-negative bacteria or because the microorganism develops resistance to the antibiotics by mutating the DHFR target, discovery of new antibiotics against the enzyme is mandatory to overcome drug-resistance. This review summarizes the field of DHFR inhibition with special focus on recent efforts to effectively interface computational and experimental efforts to discover novel classes of inhibitors that target allosteric and active-sites in drug-resistant variants of EcDHFR.

Identification of a pyrimidinetrione derivative as the potent DprE1 inhibitor by structure-based virtual ligand screening.

Despite the increasing need of new antituberculosis drugs, the number of agents approved for the market has fallen to an all-time low. In response to the emerging drug resistance followed, structurally unique chemical entities will be highlighted. decaprenylphosphoryl-ß-d-ribose oxidase (DprE1) participating in the biosynthesis of mycobacterium cell wall is a highly vulnerable and validated antituberculosis target. On the basis of it, a systematic strategy was applied to identify a high-quality lead compound (compound 50) that inhibits the essential enzyme DprE1, thus blocking the synthesis of the mycobacterial cell wall to kill M. tuberculosis in vitro and in vivo. Correspondingly, the rational design and synthetic strategy for compound 50 was reported. Notably, the compound 50 has been confirmed to be no toxicity. Altogether, our data suggest the compound 50 targeting DprE1 is a promising candidate for the tuberculosis (TB) therapy.

Study of adenylyl cyclase-GαS interactions and identification of novel AC ligands.

Adenylyl cyclases (ACs) are membrane bound enzymes that catalyze the production of cAMP from ATP in response to the activation by G-protein Gαs. Different isoforms of ACs are ubiquitously expressed in different tissues involved in regulatory mechanisms in response to specific stimulants. There are 9 AC isoforms present in humans, with AC5 and AC6 proposed to play a vital role in cardiac functions. The activity of AC6 is sensitive to nitric oxide, such that nitrosylation of the protein might regulate its function. However, the information on structural determinants of nitrosylation in ACs and how they interact with Gαs is limited. Here we used homology modeling to build a molecular model of human AC6 bound to Gαs. Based on this 3D model, we predict the nitrosylation amenable cysteines, and identify potential novel ligands of AC6 using virtual ligand screening. Our model suggests Cys1004 in AC6 (subunit C2) and Cys174 in Gαs present at the AC-Gαs interface as the possible residues that might undergo reversible nitrosylation. Docking analysis predicted novel ligands of AC6 that include forskolin-based compounds and its derivatives. Further work involving site-directed mutagenesis of the predicted residues will allow manipulation of AC activity using novel ligands, and crucial insights on the role of nitrosylation of these proteins in pathophysiological conditions.

Affinity capillary electrophoresis for identification of active drug candidates in myotonic dystrophy type 1.

Myotonic dystrophy type 1 (DM1) is an autosomal dominantly inherited degenerative disease with a slow progression. At the present, there is no commercially available treatment, but sustained effort is currently undertaken for the development of a promising lead compound. In the present paper we report the development of a fast, versatile, and cost-effective affinity capillary electrophoresis (ACE) method for the screening and identification of potential drug candidates targeting pathological ARN probes relevant for DM1. The affinity studies were conducted in physiologically relevant conditions using 50 mM HEPES buffer (pH 7.4) in a fused silica capillary dynamically coated with poly(ethylene oxide), by testing a library of potential ligands against (CUG)50 RNA as target probe with a total run time of 4-5 h/ligand. For the most promising ligands, their affinity parameters were assessed and some results formerly reported on the affinity of pentamidine (PTMD) and neomycin against CUG repeats were confirmed. To the best of the authors' knowledge, the estimated binding stoichiometry for some of the tested compounds (i.e., ~ 121:1 for PTMD against the tested RNA probe) is reported for the first time. Additionally, the potential of a novel pentamidine like compound, namely 1,2-ethane bis-1-amino-4-benzamidine (EBAB) with much lower in vivo toxicity than its parent compound has also been confirmed studying its effect on a live cell model by fluorescence microscopy. Further tests, such as the evaluation of the rescue in the mis-splicing of the involved genes, can be performed to corroborate the potential therapeutic value of EBAB in DM1 treatment. Graphical abstract ᅟ.

Steroids from Ganoderma sinense as new natural inhibitors of cancer-associated mutant IDH1.

Isocitrate dehydrogenase (IDH) is one of the key enzymes in the tricarboxylic acid cycle, and IDH mutations have been associated with many cancers, including glioblastoma, sarcoma, acute myeloid leukemia, etc. Three natural steroids 1-3 from Ganoderma sinense, a unique and rare edible-medicinal fungi in China, were found as potential IDH1 inhibitors by virtual ligand screening method. Among the three compounds, 3 showed the highest binding affinity to IDH1 with significant calculated binding free energy. Enzymatic kinetics demonstrated that 3 inhibited mutant enzyme in a noncompetitive manner. The half effective concentration of 3 for reducing the concentration of D-2HG in HT1080 cells was 35.97 µM. The levels of histone H3K9me3 methylation in HT1080 cells were reduced by treating with 3. Furthermore, knockdown of mutant IDH1 in HT1080 cells decreased the anti-proliferative sensitivity to 3. In short, our findings highlight that compound 3 may have clinical potential in tumor therapies as an effective inhibitor of mutant IDH1.

[In silico Search for Tubulin Polymerization Inhibitors].

Cytostatic colchicine is widely used in the treatment of Familial Mediterranean fever, but it has several side effects. For finding new, more effective drugs with higher affinity and diminishside effects we carried out virtual screening of potential inhibitors of the main target of colchicine, the polymerization of tubulin by evaluating affinity 25745 compounds, structurally related to the colchicine. We have identified 11 commercially available compounds with higher affinity to tubulin. Compounds with highest binding scores include trimethoxybenzene and its derivatives; these compounds bind to the same site in similar orientation. Information provided can form the basis for design of new cytostatics.