Fragment-based drug discovery and its application to challenging drug targets

The "Doorstop Pocket" In Thioredoxin Reductases─An Unexpected Druggable Regulator of the Catalytic Machinery

Pyridine nucleotide-disulfide oxidoreductases are underexplored as drug targets, and thioredoxin reductases (TrxRs) stand out as compelling pharmacological targets. Selective TrxR inhibition is challenging primarily due to the reliance on covalent inhibition strategies. Recent studies identified a regulatory and druggable pocket in Schistosoma mansoni thioredoxin glutathione reductase (TGR), a TrxR-like enzyme, and an established drug target for schistosomiasis. This site is termed the "doorstop pocket" because compounds that bind there impede the movement of an aromatic side-chain necessary for the entry and exit of NADPH and NADP+ during enzymatic turnover. This discovery spearheaded the development of new TGR inhibitors with efficacies surpassing those of current schistosomiasis treatment. Targeting the "doorstop pocket" is a promising strategy, as the pocket is present in all members of the pyridine nucleotide-disulfide oxidoreductase family, opening new avenues for exploring therapeutic approaches in diseases where the importance of these enzymes is established, including cancer and inflammatory and infectious diseases.

Crystallographic fragment-binding studies of the Mycobacterium tuberculosis trifunctional enzyme suggest binding pockets for the tails of the acyl-CoA substrates at its active sites and a potential substrate-channeling path between them

The Mycobacterium tuberculosis trifunctional enzyme (MtTFE) is an α2β2 tetrameric enzyme in which the α-chain harbors the 2E-enoyl-CoA hydratase (ECH) and 3S-hydroxyacyl-CoA dehydrogenase (HAD) active sites, and the β-chain provides the 3-ketoacyl-CoA thiolase (KAT) active site. Linear, medium-chain and long-chain 2E-enoyl-CoA molecules are the preferred substrates of MtTFE. Previous crystallographic binding and modeling studies identified binding sites for the acyl-CoA substrates at the three active sites, as well as the NAD binding pocket at the HAD active site. These studies also identified three additional CoA binding sites on the surface of MtTFE that are different from the active sites. It has been proposed that one of these additional sites could be of functional relevance for the substrate channeling (by surface crawling) of reaction intermediates between the three active sites. Here, 226 fragments were screened in a crystallographic fragment-binding study of MtTFE crystals, resulting in the structures of 16 MtTFE-fragment complexes. Analysis of the 121 fragment-binding events shows that the ECH active site is the `binding hotspot' for the tested fragments, with 41 binding events. The mode of binding of the fragments bound at the active sites provides additional insight into how the long-chain acyl moiety of the substrates can be accommodated at their proposed binding pockets. In addition, the 20 fragment-binding events between the active sites identify potential transient binding sites of reaction intermediates relevant to the possible channeling of substrates between these active sites. These results provide a basis for further studies to understand the functional relevance of the latter binding sites and to identify substrates for which channeling is crucial.

Using DNA-encoded libraries of fragments for hit discovery of challenging therapeutic targets

INTRODUCTION

The effectiveness of Fragment-based drug design (FBDD) for targeting challenging therapeutic targets has been hindered by two factors: the small library size and the complexity of the fragment-to-hit optimization process. The DNA-encoded library (DEL) technology offers a compelling and robust high-throughput selection approach to potentially address these limitations.

AREA COVERED

In this review, the authors propose the viewpoint that the DEL technology matches perfectly with the concept of FBDD to facilitate hit discovery. They begin by analyzing the technical limitations of FBDD from a medicinal chemistry perspective and explain why DEL may offer potential solutions to these limitations. Subsequently, they elaborate in detail on how the integration of DEL with FBDD works. In addition, they present case studies involving both de novo hit discovery and full ligand discovery, especially for challenging therapeutic targets harboring broad drug-target interfaces.

EXPERT OPINION

The future of DEL-based fragment discovery may be promoted by both technical advances and application scopes. From the technical aspect, expanding the chemical diversity of DEL will be essential to achieve success in fragment-based drug discovery. From the application scope side, DEL-based fragment discovery holds promise for tackling a series of challenging targets.

Streamlining Computational Fragment-Based Drug Discovery through Evolutionary Optimization Informed by Ligand-Based Virtual Prescreening

Recent advances in computational methods provide the promise of dramatically accelerating drug discovery. While mathematical modeling and machine learning have become vital in predicting drug-target interactions and properties, there is untapped potential in computational drug discovery due to the vast and complex chemical space. This paper builds on our recently published computational fragment-based drug discovery (FBDD) method called fragment databases from screened ligand drug discovery (FDSL-DD). FDSL-DD uses in silico screening to identify ligands from a vast library, fragmenting them while attaching specific attributes based on predicted binding affinity and interaction with the target subdomain. In this paper, we further propose a two-stage optimization method that utilizes the information from prescreening to optimize computational ligand synthesis. We hypothesize that using prescreening information for optimization shrinks the search space and focuses on promising regions, thereby improving the optimization for candidate ligands. The first optimization stage assembles these fragments into larger compounds using genetic algorithms, followed by a second stage of iterative refinement to produce compounds with enhanced bioactivity. To demonstrate broad applicability, the methodology is demonstrated on three diverse protein targets found in human solid cancers, bacterial antimicrobial resistance, and the SARS-CoV-2 virus. Combined, the proposed FDSL-DD and a two-stage optimization approach yield high-affinity ligand candidates more efficiently than other state-of-the-art computational FBDD methods. We further show that a multiobjective optimization method accounting for drug-likeness can still produce potential candidate ligands with a high binding affinity. Overall, the results demonstrate that integrating detailed chemical information with a constrained search framework can markedly optimize the initial drug discovery process, offering a more precise and efficient route to developing new therapeutics.

PDBe CCDUtils: an RDKit-based toolkit for handling and analysing small molecules in the Protein Data Bank

While the Protein Data Bank (PDB) contains a wealth of structural information on ligands bound to macromolecules, their analysis can be challenging due to the large amount and diversity of data. Here, we present PDBe CCDUtils, a versatile toolkit for processing and analysing small molecules from the PDB in PDBx/mmCIF format. PDBe CCDUtils provides streamlined access to all the metadata for small molecules in the PDB and offers a set of convenient methods to compute various properties using RDKit, such as 2D depictions, 3D conformers, physicochemical properties, scaffolds, common fragments, and cross-references to small molecule databases using UniChem. The toolkit also provides methods for identifying all the covalently attached chemical components in a macromolecular structure and calculating similarity among small molecules. By providing a broad range of functionality, PDBe CCDUtils caters to the needs of researchers in cheminformatics, structural biology, bioinformatics and computational chemistry.

The current role and evolution of X-ray crystallography in drug discovery and development

INTRODUCTION

Macromolecular X-ray crystallography and cryo-EM are currently the primary techniques used to determine the three-dimensional structures of proteins, nucleic acids, and viruses. Structural information has been critical to drug discovery and structural bioinformatics. The integration of artificial intelligence (AI) into X-ray crystallography has shown great promise in automating and accelerating the analysis of complex structural data, further improving the efficiency and accuracy of structure determination.

AREAS COVERED

This review explores the relationship between X-ray crystallography and other modern structural determination methods. It examines the integration of data acquired from diverse biochemical and biophysical techniques with those derived from structural biology. Additionally, the paper offers insights into the influence of AI on X-ray crystallography, emphasizing how integrating AI with experimental approaches can revolutionize our comprehension of biological processes and interactions.

EXPERT OPINION

Investing in science is crucially emphasized due to its significant role in drug discovery and advancements in healthcare. X-ray crystallography remains an essential source of structural biology data for drug discovery. Recent advances in biochemical, spectroscopic, and bioinformatic methods, along with the integration of AI techniques, hold the potential to revolutionize drug discovery when effectively combined with robust data management practices.

Structurally modified Cyclovirobuxine-D Buxus alkaloids as effective analgesic agents through Cav3.2 T-Type calcium channel inhibition

Cyclovirobuxine-D (CVB-D) is a Buxus alkaloid and a major active constituent in the Chinese medicinal herb Buxus microphylls. Traditionally, the natural alkaloid cyclovirobuxine-D has a long history of use as a traditional Chinese medicine for cardiovascular diseases as well as to treat a wide variety of medical conditions. As we found that CVB-D inhibited T-type calcium channels, we designed and synthesized a variety of fragments and analogues and evaluated them for the first time as new Cav3.2 inhibitors. Compounds 2-7 exhibited potency against Cav 3.2 channels, and two of them were more active than their parent molecules. As a result of the in vivo experiments, both compounds 3 and 4 showed significantly reduced writhes in the acetic acid-induced writhing test. Studies of molecular modeling have identified possible mechanism(s) of Cav3.2 binding. Moreover, the relationship between structure and activity was studied in a preliminary manner. Our results indicated that compounds 3 and 4 could play an important role in the discovery and development of novel analgesics.

Emerging approaches to CDK inhibitor development, a structural perspective

Aberrant activity of the cyclin-dependent kinase family is frequently noted in a number of diseases identifying them as potential targets for drug development. However, current CDK inhibitors lack specificity owing to the high sequence and structural conservation of the ATP binding cleft across family members, highlighting the necessity of finding novel modes of CDK inhibition. The wealth of structural information regarding CDK assemblies and inhibitor complexes derived from X-ray crystallographic studies has been recently complemented through the use of cryo-electron microscopy. These recent advances have provided insights into the functional roles and regulatory mechanisms of CDKs and their interaction partners. This review explores the conformational malleability of the CDK subunit, the importance of SLiM recognition sites in CDK complexes, the progress made in chemically induced CDK degradation and how these studies can contribute to CDK inhibitor design. Additionally, fragment-based drug discovery can be utilised to identify small molecules that bind to allosteric sites on the CDK surface employing interactions which mimic those of native protein-protein interactions. These recent structural advances in CDK inhibitor mechanisms and in chemical probes which do not occupy the orthosteric ATP binding site can provide important insights for targeted CDK therapies.

Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here, we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF (UniProtKB Q9X0C6). We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 μM. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve of the receiver operating characteristics of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.

Exploring species-level infant gut bacterial biodiversity by meta-analysis and formulation of an optimized cultivation medium

In vitro gut cultivation models provide host-uncoupled, fast, and cost-efficient solutions to investigate the effects of intrinsic and extrinsic factors impacting on both composition and functionality of the intestinal microbial ecosystem. However, to ensure the maintenance and survival of gut microbial players and preserve their functions, these systems require close monitoring of several variables, including oxygen concentration, pH, and temperature, as well as the use of a culture medium satisfying the microbial nutritional requirements. In this context, in order to identify the macro- and micro-nutrients necessary for in vitro cultivation of the infant gut microbiota, a meta-analysis based on 1669 publicly available shotgun metagenomic samples corresponding to fecal samples of healthy, full-term infants aged from a few days to three years was performed to define the predominant species characterizing the “infant-like” gut microbial ecosystem. A subsequent comparison of growth performances was made using infant fecal samples that contained the most abundant bacterial taxa of the infant gut microbiota, when cultivated on 18 different culture media. This growth analysis was performed by means of flow cytometry-based bacterial cell enumeration and shallow shotgun sequencing, which allowed the formulation of an optimized growth medium, i.e., Infant Gut Super Medium (IGSM), which maintains and sustains the infant gut microbial biodiversity under in vitro growth conditions. Furthermore, this formulation was used to evaluate the in vitro effect of two drugs commonly used in pediatrics, i.e., acetaminophen and simethicone, on the taxonomic composition of the infant gut microbiota.

Structural Analysis and Development of Notum Fragment Screening Hits

The Wnt signaling suppressor Notum is a promising target for osteoporosis, Alzheimer's disease, and colorectal cancers. To develop novel Notum inhibitors, we used an X-ray crystallographic fragment screen with the Diamond-SGC Poised Library (DSPL) and identified 59 fragment hits from the analysis of 768 data sets. Fifty-eight of the hits were found bound at the enzyme catalytic pocket with potencies ranging from 0.5 to >1000 μM. Analysis of the fragments' diverse binding modes, enzymatic inhibitory activities, and chemical properties led to the selection of six hits for optimization, and five of these resulted in improved Notum inhibitory potencies. One hit, 1-phenyl-1,2,3-triazole 7, and its related cluster members, have shown promising lead-like properties. These became the focus of our fragment development activities, resulting in compound 7d with IC50 0.0067 μM. The large number of Notum fragment structures and their initial optimization provided an important basis for further Notum inhibitor development.

Fast fragment- and compound-screening pipeline at the Swiss Light Source

Over the last two decades, fragment-based drug discovery (FBDD) has emerged as an effective and efficient method to identify new chemical scaffolds for the development of lead compounds. X-ray crystallography can be used in FBDD as a tool to validate and develop fragments identified as binders by other methods. However, it is also often used with great success as a primary screening technique. In recent years, technological advances at macromolecular crystallography beamlines in terms of instrumentation, beam intensity and robotics have enabled the development of dedicated platforms at synchrotron sources for FBDD using X-ray crystallography. Here, the development of the Fast Fragment and Compound Screening (FFCS) platform, an integrated next-generation pipeline for crystal soaking, handling and data collection which allows crystallography-based screening of protein crystals against hundreds of fragments and compounds, at the Swiss Light Source is reported.

Dual-target inhibitors of poly (ADP-ribose) polymerase-1 for cancer therapy: Advances, challenges, and opportunities

PARP1 plays a crucial role in DNA damage repair, making it an essential target for cancer therapy. PARP1 inhibitors are widely used to treat BRCA-deficient malignancies, and six PARP inhibitors have been approved for clinical use. However, excluding the great clinical success of PARP inhibitors, the concomitant toxicity, drug resistance, and limited scope of application restrict their clinical efficacy. To find solutions to these problems, dual-target inhibitors have shown great potential. In recent years, several studies have linked PAPR1 to other primary cancer targets. Many dual-target inhibitors have been developed using structural fusion, linkage, or library construction methods, overcoming the defects of many single-target inhibitors of PARP1 and achieving great success in clinical cancer therapy. This review summarizes the advance of dual-target PARP1 inhibitors in recent years, focusing on their structural optimization process, structure-activity relationships (SARs), and in vitro or in vivo analysis results.

A Lead-Based Fragment Library Screening of the Glycosyltransferase WaaG from Escherichia coli

Glucosyl transferase I (WaaG) in E. coli catalyzes the transfer of an α-d-glucosyl group to the inner core of the lipopolysaccharide (LPS) and plays an important role in the biogenesis of the outer membrane. If its activity could be inhibited, the integrity of the outer membrane would be compromised and the bacterium would be susceptible to antibiotics that are normally prevented from entering the cell. Herein, three libraries of molecules (A, B and C) were docked in the binding pocket of WaaG, utilizing the docking binding affinity as a filter to select fragment-based compounds for further investigations. From the results of the docking procedure, a selection of compounds was investigated by molecular dynamics (MD) simulations to obtain binding free energy (BFE) and K_D values for ligands as an evaluation for the binding to WaaG. Derivatives of 1,3-thiazoles (A7 and A4) from library A and 1,3,4-thiadiazole (B33) from library B displayed a promising profile of BFE, with K_D < mM, viz., 0.11, 0.62 and 0.04 mM, respectively. Further root-mean-square-deviation (RMSD), electrostatic/van der Waals contribution to the binding and H-bond interactions displayed a favorable profile for ligands A4 and B33. Mannose and/or heptose-containing disaccharides C1–C4, representing sub-structures of the inner core of the LPS, were also investigated by MD simulations, and compound C4²⁻ showed a calculated K_D = 0.4 µM. In the presence of UDP-Glc²⁻, the best-docked pose of disaccharide C4²⁻ is proximate to the glucose-binding site of WaaG. A study of the variation in angle and distance was performed on the different portions of WaaG (N-, the C- domains and the hinge region). The Spearman correlation coefficient between the two variables was close to unity, where both variables increase in the same way, suggesting a conformational rearrangement of the protein during the MD simulation, revealing molecular motions of the enzyme that may be part of the catalytic cycle. Selected compounds were also analyzed by Saturation Transfer Difference (STD) NMR experiments. STD effects were notable for the 1,3-thiazole derivatives A4, A8 and A15 with the apo form of the protein as well as in the presence of UDP for A4.

Identification of a small compound that specifically inhibits Zika virus in vitro and in vivo by targeting the NS2B-NS3 protease

Zika virus (ZIKV) has rapid become a global threat, but no ZIKV-specific vaccines or drugs are currently available. In this study, inhibitors of ZIKV NS2B-NS3 protease were screened from a library containing 4,452 compound fragments. One of the compounds, 6-bromo-1,2-naphthalenedione, exhibited high specific inhibition against ZIKV NS2B-NS3 protease, but had no inhibitory effects against other viral proteases. A microscale thermophoresis (MST) assay confirmed that the compound bound to ZIKV NS2B-NS3 protein with a binding constant (Kd) of 12.26 μM. Indirect immunofluorescence assays, Western blots, and plaque assays indicated that the compound inhibited virus replication in cells. Virus titer was reduced by more than 75% when the compound was present at 1 μM. A time-of-addition assay showed that inhibition occurred at the virus replication stage, but not at the adsorption or invasion stages. The half cytotoxicity concentration (CC50) of the compound on HeLa, Vero, and BHK-21 cells were 445.44 μM, 123.87 μM, and 123.64 μM, respectively. In vivo tests using infected AG129 mice demonstrated that treatment with the compound reduced mortality by up to 60%. Mice treated with the compound showed a reduction in histopathological lesions in brain, testis, and ovary. Viral RNA, IL-1β, and IL-6 mRNA levels decreased significantly in these tissues. In summary, this study has identified a small compound with high and specific inhibitory effects on ZIKV. The compound can be used as a therapeutic agent and is also an ideal starting point for drug optimization.

Fragment based drug design

Rational approaches towards drug development have emerged as one of the most promising ways among the tedious conventional procedures with the aim of redefining the drug discovery process. The need of current medical system is demanding a much precise, faster and reliable approaches in parallel to faster growing technology for development of drugs with more intrinsic action and fewer side effects. Systematic and well-defined diagnostic studies have revealed the specific causes of disease and related targets for drug development. Designing a drug as per the specific target, studying it in-silico prior to its development has been proved as an added benefit to the studies. Many approaches like structure based drug design, fragment based drug design and ligand based drug design are been in practice for the drug discovery and development with the similar fundamental theory. Fragment based drug design utilizes a library of fragments designed specifically for the concerned target and these fragments are studied further before screening with virtual methods as well as with biophysical methods. The process follows a well-defined pathway which moulds a fragment into a perfect drug candidate. In this chapter we have tried to cover all the basic aspects of fragment based drug design and related technologies.

Peptide Inhibitors of Kv1.5: An Option for the Treatment of Atrial Fibrillation

The human voltage gated potassium channel Kv1.5 that conducts the IKur current is a key determinant of the atrial action potential. Its mutations have been linked to hereditary forms of atrial fibrillation (AF), and the channel is an attractive target for the management of AF. The development of IKur blockers to treat AF resulted in small molecule Kv1.5 inhibitors. The selectivity of the blocker for the target channel plays an important role in the potential therapeutic application of the drug candidate: the higher the selectivity, the lower the risk of side effects. In this respect, small molecule inhibitors of Kv1.5 are compromised due to their limited selectivity. A wide range of peptide toxins from venomous animals are targeting ion channels, including mammalian channels. These peptides usually have a much larger interacting surface with the ion channel compared to small molecule inhibitors and thus, generally confer higher selectivity to the peptide blockers. We found two peptides in the literature, which inhibited IKur: Ts6 and Osu1. Their affinity and selectivity for Kv1.5 can be improved by rational drug design in which their amino acid sequences could be modified in a targeted way guided by in silico docking experiments.

Fragment-Based Ligand Discovery Applied to the Mycolic Acid Methyltransferase Hma (MmaA4) from Mycobacterium tuberculosis: A Crystallographic and Molecular Modelling Study

The mycolic acid biosynthetic pathway represents a promising source of pharmacological targets in the fight against tuberculosis. In Mycobacterium tuberculosis, mycolic acids are subject to specific chemical modifications introduced by a set of eight S-adenosylmethionine dependent methyltransferases. Among these, Hma (MmaA4) is responsible for the introduction of oxygenated modifications. Crystallographic screening of a library of fragments allowed the identification of seven ligands of Hma. Two mutually exclusive binding modes were identified, depending on the conformation of residues 147–154. These residues are disordered in apo-Hma but fold upon binding of the S-adenosylmethionine (SAM) cofactor as well as of analogues, resulting in the formation of the short η1-helix. One of the observed conformations would be incompatible with the presence of the cofactor, suggesting that allosteric inhibitors could be designed against Hma. Chimeric compounds were designed by fusing some of the bound fragments, and the relative binding affinities of initial fragments and evolved compounds were investigated using molecular dynamics simulation and generalised Born and Poisson–Boltzmann calculations coupled to the surface area continuum solvation method. Molecular dynamics simulations were also performed on apo-Hma to assess the structural plasticity of the unliganded protein. Our results indicate a significant improvement in the binding properties of the designed compounds, suggesting that they could be further optimised to inhibit Hma activity.

Fragment- and structure-based drug discovery for developing therapeutic agents targeting the DNA Damage Response

Cancer will directly affect the lives of over one-third of the population. The DNA Damage Response (DDR) is an intricate system involving damage recognition, cell cycle regulation, DNA repair, and ultimately cell fate determination, playing a central role in cancer etiology and therapy. Two primary therapeutic approaches involving DDR targeting include: combinatorial treatments employing anticancer genotoxic agents; and synthetic lethality, exploiting a sporadic DDR defect as a mechanism for cancer-specific therapy. Whereas, many DDR proteins have proven "undruggable", Fragment- and Structure-Based Drug Discovery (FBDD, SBDD) have advanced therapeutic agent identification and development. FBDD has led to 4 (with ∼50 more drugs under preclinical and clinical development), while SBDD is estimated to have contributed to the development of >200, FDA-approved medicines. Protein X-ray crystallography-based fragment library screening, especially for elusive or "undruggable" targets, allows for simultaneous generation of hits plus details of protein-ligand interactions and binding sites (orthosteric or allosteric) that inform chemical tractability, downstream biology, and intellectual property. Using a novel high-throughput crystallography-based fragment library screening platform, we screened five diverse proteins, yielding hit rates of ∼2-8% and crystal structures from ∼1.8 to 3.2 Å. We consider current FBDD/SBDD methods and some exemplary results of efforts to design inhibitors against the DDR nucleases meiotic recombination 11 (MRE11, a.k.a., MRE11A), apurinic/apyrimidinic endonuclease 1 (APE1, a.k.a., APEX1), and flap endonuclease 1 (FEN1).

Probing the Surface of a Parasite Drug Target Thioredoxin Glutathione Reductase Using Small Molecule Fragments

Fragment screening is a powerful drug discovery approach particularly useful for enzymes difficult to inhibit selectively, such as the thiol/selenol-dependent thioredoxin reductases (TrxRs), which are essential and druggable in several infectious diseases. Several known inhibitors are reactive electrophiles targeting the selenocysteine-containing C-terminus and thus often suffering from off-target reactivity in vivo. The lack of structural information on the interaction modalities of the C-terminus-targeting inhibitors, due to the high mobility of this domain and the lack of alternative druggable sites, prevents the development of selective inhibitors for TrxRs. In this work, fragments selected from actives identified in a large screen carried out against Thioredoxin Glutathione Reductase from Schistosoma mansoni (SmTGR) were probed by X-ray crystallography. SmTGR is one of the most promising drug targets for schistosomiasis, a devastating, neglected disease. Utilizing a multicrystal method to analyze electron density maps, structural analysis, and functional studies, three binding sites were characterized in SmTGR: two sites are close to or partially superposable with the NADPH binding site, while the third one is found between two symmetry related SmTGR subunits of the crystal lattice. Surprisingly, one compound bound to this latter site stabilizes, through allosteric effects mediated by the so-called guiding bar residues, the crucial redox active C-terminus of SmTGR, making it finally visible at high resolution. These results further promote fragments as small molecule probes for investigating functional aspects of the target protein, exemplified by the allosteric effect on the C-terminus, and providing fundamental chemical information exploitable in drug discovery.

CASTELO: clustered atom subtypes aided lead optimization-a combined machine learning and molecular modeling method

Background

Drug discovery is a multi-stage process that comprises two costly major steps: pre-clinical research and clinical trials. Among its stages, lead optimization easily consumes more than half of the pre-clinical budget. We propose a combined machine learning and molecular modeling approach that partially automates lead optimization workflow in silico, providing suggestions for modification hot spots.

Results

The initial data collection is achieved with physics-based molecular dynamics simulation. Contact matrices are calculated as the preliminary features extracted from the simulations. To take advantage of the temporal information from the simulations, we enhanced contact matrices data with temporal dynamism representation, which are then modeled with unsupervised convolutional variational autoencoder (CVAE). Finally, conventional and CVAE-based clustering methods are compared with metrics to rank the submolecular structures and propose potential candidates for lead optimization.

Conclusion

With no need for extensive structure-activity data, our method provides new hints for drug modification hotspots which can be used to improve drug potency and reduce the lead optimization time. It can potentially become a valuable tool for medicinal chemists.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04214-4.

Identification of new building blocks by fragment screening for discovering GyrB inhibitors

DNA gyrase is an essential DNA topoisomerase that exists only in bacteria. Since novobiocin was withdrawn from the market, new scaffolds and new mechanistic GyrB inhibitors are urgently needed. In this study, we employed fragment screening and X-ray crystallography to identify new building blocks, as well as their binding mechanisms, to support the discovery of new GyrB inhibitors. In total, 84 of the 618 chemical fragments were shown to either thermally stabilize the ATPase domain of Escherichia coli GyrB or inhibit the ATPase activity of E. coli gyrase. Among them, the IC₅₀ values of fragments 10 and 23 were determined to be 605.3 μM and 446.2 μM, respectively. Cocrystal structures of the GyrB ATPase domain with twelve fragment hits were successfully determined at a high resolution. All twelve fragments were deeply inserted in the pocket and formed H-bonds with Asp73 and Thr165, and six fragments formed an additional H-bond with the backbone oxygen of Val71. Fragment screening further highlighted the capability of Asp73, Thr165 and Val71 to bind chemicals and provided diverse building blocks for the design of GyrB inhibitors.

Divergent On-DNA Transformations from DNA-Linked Piperidones

A group of highly efficient and divergent transformations for constructing multiple DNA-linked chemotypes based on a piperidone core were successfully developed. We reported the first procedure for the synthesis of a DNA-conjugated piperidine intermediate under basic conditions. Subsequently, this substructure was subjected to additional reactions to generate several privileged scaffolds, including 4-aminopiperidine, fused [1,2,4]triazolo[1,5-a]pyrimidine, and a quinoline derivative. These transformations paved the way for constructing focused scaffold-based DNA-encoded libraries with druglike properties.

Organic Antifungal Drugs and Targets of Their Action

In recent decades, there has been a significant increase in the number of fungal diseases. This is due to a wide spectrum of action, immunosuppressants and other group drugs. In terms of frequency, rapid spread and globality, fungal infections are approaching acute respiratory infections. Antimycotics are medicinal substances endorsed with fungicidal or fungistatic properties. For the treatment of fungal diseases, several groups of compounds are used that differ in their origin (natural or synthetic), molecular targets and mechanism of action, antifungal effect (fungicidal or fungistatic), indications for use (local or systemic infections), and methods of administration (parenteral, oral, outdoor). Several efforts have been made by various medicinal chemists around the world for the development of antifungal drugs with high efficacy with the least toxicity and maximum selectivity in the area of antifungal chemotherapy. The pharmacokinetic properties of the new antimycotics are also important: the ability to penetrate biological barriers, be absorbed and distributed in tissues and organs, get accumulated in tissues affected by micromycetes, undergo drug metabolism in the intestinal microflora and human organs, and in the kinetics of excretion from the body. There are several ways to search for new effective antimycotics: - Obtaining new derivatives of the already used classes of antimycotics with improved activity properties. - Screening of new chemical classes of synthetic antimycotic compounds. - Screening of natural compounds. - Identification of new unique molecular targets in the fungal cell. - Development of new compositions and dosage forms with effective delivery vehicles. The methods of informatics, bioinformatics, genomics and proteomics were extensively investigated for the development of new antimycotics. These techniques were employed in finding and identification of new molecular proteins in a fungal cell; in the determination of the selectivity of drugprotein interactions, evaluation of drug-drug interactions and synergism of drugs; determination of the structure-activity relationship (SAR) studies; determination of the molecular design of the most active, selective and safer drugs for the humans, animals and plants. In medical applications, the methods of information analysis and pharmacogenomics allow taking into account the individual phenotype of the patient, the level of expression of the targets of antifungal drugs when choosing antifungal agents and their dosage. This review article incorporates some of the most significant studies covering the basic structures and approaches for the synthesis of antifungal drugs and the directions for their further development.

MCSS-Based Predictions of Binding Mode and Selectivity of Nucleotide Ligands

Computational fragment-based approaches are widely used in drug design and discovery. One of their limitations is the lack of performance of docking methods, mainly the scoring functions. With the emergence of fragment-based approaches for single-stranded RNA ligands, we analyze the performance in docking and screening powers of an MCSS-based approach. The performance is evaluated on a benchmark of protein-nucleotide complexes where the four RNA residues are used as fragments. The screening power can be considered the major limiting factor for the fragment-based modeling or design of sequence-selective oligonucleotides. We show that the MCSS sampling is efficient even for such large and flexible fragments. Hybrid solvent models based on some partial explicit representations improve both the docking and screening powers. Clustering of the n best-ranked poses can also contribute to a lesser extent to better performance. A detailed analysis of molecular features suggests various ways to optimize the performance further.

Design, synthesis, and biological evaluation of quinazolin-4(3H)-one derivatives co-targeting poly(ADP-ribose) polymerase-1 and bromodomain containing protein 4 for breast cancer therapy

This study was aimed to design the first dual-target small-molecule inhibitor co-targeting poly (ADP-ribose) polymerase-1 (PARP1) and bromodomain containing protein 4 (BRD4), which had important cross relation in the global network of breast cancer, reflecting the synthetic lethal effect. A series of new BRD4 and PARP1 dual-target inhibitors were discovered and synthesized by fragment-based combinatorial screening and activity assays that together led to the chemical optimization. Among these compounds, 19d was selected and exhibited micromole enzymatic potencies against BRD4 and PARP1, respectively. Compound 19d was further shown to efficiently modulate the expression of BRD4 and PARP1. Subsequently, compound 19d was found to induce breast cancer cell apoptosis and stimulate cell cycle arrest at G1 phase. Following pharmacokinetic studies, compound 19d showed its antitumor activity in breast cancer susceptibility gene 1/2 (BRCA1/2) wild-type MDA-MB-468 and MCF-7 xenograft models without apparent toxicity and loss of body weight. These results together demonstrated that a highly potent dual-targeted inhibitor was successfully synthesized and indicated that co-targeting of BRD4 and PARP1 based on the concept of synthetic lethality would be a promising therapeutic strategy for breast cancer.

LEADS-FRAG: A Benchmark Data Set for Assessment of Fragment Docking Performance

Fragment-based drug design is a popular approach in drug discovery, which makes use of computational methods such as molecular docking. To assess fragment placement performance of molecular docking programs, we constructed LEADS-FRAG, a benchmark data set containing 93 high-quality protein-fragment complexes that were selected from the Protein Data Bank using a rational and unbiased process. The data set contains fully prepared protein and fragment structures and is publicly available. Moreover, we used LEADS-FRAG for evaluating the small-molecule docking programs AutoDock, AutoDock Vina, FlexX, and GOLD for their fragment docking performance. GOLD in combination with the scoring function ChemPLP and AutoDock Vina performed best and generated near-native conformations (root mean square deviation <1.5 Å) for more than 50% of the data set considering the top-ranked docking pose. Taking into account all docking poses, the tested programs generated near-native conformations for up to 86% of the fragments in LEADS-FRAG. By rescoring all docking poses with the GOLD scoring functions and the Protein-Ligand Informatics force field, the number of near-native conformations increased up to 40% with respect to the top-rescored poses. Our results show that conventional small-molecule docking programs achieve a satisfactory fragment docking performance when utilizing rescoring.

Fragment-Based Stabilizers of Protein-Protein Interactions through Imine-Based Tethering

Small‐molecule stabilization of protein–protein interactions (PPIs) is a promising concept in drug discovery, however the question how to identify or design chemical starting points in a “bottom‐up” approach is largely unanswered. We report a novel concept for identifying initial chemical matter for PPI stabilization based on imine‐forming fragments. The imine bond offers a covalent anchor for site‐directed fragment targeting, whereas its transient nature enables efficient analysis of structure–activity relationships. This bond enables fragment identification and optimisation using protein crystallography. We report novel fragments that bind specifically to a lysine at the PPI interface of the p65‐subunit‐derived peptide of NF‐κB with the adapter protein 14‐3‐3. Those fragments that subsequently establish contacts with the p65‐derived peptide, rather than with 14‐3‐3, efficiently stabilize the 14‐3‐3/p65 complex and offer novel starting points for molecular glues.

Fragment-based Differential Targeting of PPI Stabilizer Interfaces

Stabilization of protein-protein interactions (PPIs) holds great potential for therapeutic agents, as illustrated by the successful drugs rapamycin and lenalidomide. However, how such interface-binding molecules can be created in a rational, bottom-up manner is a largely unanswered question. We report here how a fragment-based approach can be used to identify chemical starting points for the development of small-molecule stabilizers that differentiate between two different PPI interfaces of the adapter protein 14-3-3. The fragments discriminately bind to the interface of 14-3-3 with the recognition motif of either the tumor suppressor protein p53 or the oncogenic transcription factor TAZ. This X-ray crystallography driven study shows that the rim of the interface of individual 14-3-3 complexes can be targeted in a differential manner with fragments that represent promising starting points for the development of specific 14-3-3 PPI stabilizers.

MDA-9/Syntenin (SDCBP): Novel gene and therapeutic target for cancer metastasis

The primary cause of cancer-related death from solid tumors is metastasis. While unraveling the mechanisms of this complicated process continues, our ability to effectively target and treat it to decrease patient morbidity and mortality remains disappointing. Early detection of metastatic lesions and approaches to treat metastases (both pharmacological and genetic) are of prime importance to obstruct this process clinically. Metastasis is complex involving both genetic and epigenetic changes in the constantly evolving tumor cell. Moreover, many discrete steps have been identified in metastatic spread, including invasion, intravasation, angiogenesis, attachment at a distant site (secondary seeding), extravasation and micrometastasis and tumor dormancy development. Here, we provide an overview of the metastatic process and highlight a unique pro-metastatic gene, melanoma differentiation associated gene-9/Syntenin (MDA-9/Syntenin) also called syndecan binding protein (SDCBP), which is a major contributor to the majority of independent metastatic events. MDA-9 expression is elevated in a wide range of carcinomas and other cancers, including melanoma, glioblastoma multiforme and neuroblastoma, suggesting that it may provide an appropriate target to intervene in metastasis. Pre-clinical studies confirm that inhibiting MDA-9 either genetically or pharmacologically profoundly suppresses metastasis. An additional benefit to blocking MDA-9 in metastatic cells is sensitization of these cells to a second therapeutic agent, which converts anti-invasion effects to tumor cytocidal effects. Continued mechanistic and therapeutic insights hold promise to advance development of truly effective therapies for metastasis in the future.

Promising 2,6,9-Trisubstituted Purine Derivatives for Anticancer Compounds: Synthesis, 3D-QSAR, and Preliminary Biological Assays

We designed, synthesized, and evaluated novel 2,6,9-trisubstituted purine derivatives for their prospective role as antitumor compounds. Using simple and efficient methodologies, 31 compounds were obtained. We tested these compounds in vitro to draw conclusions about their cell toxicity on seven cancer cells lines and one non-neoplastic cell line. Structural requirements for antitumor activity on two different cancer cell lines were analyzed with SAR and 3D-QSAR. The 3D-QSAR models showed that steric properties could better explain the cytotoxicity of compounds than electronic properties (70% and 30% of contribution, respectively). From this analysis, we concluded that an arylpiperazinyl system connected at position 6 of the purine ring is beneficial for cytotoxic activity, while the use of bulky systems at position C-2 of the purine is not favorable. Compound 7h was found to be an effective potential agent when compared with a currently marketed drug, cisplatin, in four out of the seven cancer cell lines tested. Compound 7h showed the highest potency, unprecedented selectivity, and complied with all the Lipinski rules. Finally, it was demonstrated that 7h induced apoptosis and caused cell cycle arrest at the S-phase on HL-60 cells. Our study suggests that substitution in the purine core by arylpiperidine moiety is essential to obtain derivatives with potential anticancer activity.

High-throughput structures of protein-ligand complexes at room temperature using serial femtosecond crystallography

High-throughput X-ray crystal structures of protein-ligand complexes are critical to pharmaceutical drug development. However, cryocooling of crystals and X-ray radiation damage may distort the observed ligand binding. Serial femtosecond crystallography (SFX) using X-ray free-electron lasers (XFELs) can produce radiation-damage-free room-temperature structures. Ligand-binding studies using SFX have received only modest attention, partly owing to limited beamtime availability and the large quantity of sample that is required per structure determination. Here, a high-throughput approach to determine room-temperature damage-free structures with excellent sample and time efficiency is demonstrated, allowing complexes to be characterized rapidly and without prohibitive sample requirements. This yields high-quality difference density maps allowing unambiguous ligand placement. Crucially, it is demonstrated that ligands similar in size or smaller than those used in fragment-based drug design may be clearly identified in data sets obtained from <1000 diffraction images. This efficiency in both sample and XFEL beamtime opens the door to true high-throughput screening of protein-ligand complexes using SFX.

BIM-46174 fragments as potential ligands of G proteins

The 5,6,7,8-tetrahydroimidazo[1,2-a]pyrazine derivative BIM-46174 has received attention as Gα_q inhibitor. We conducted structural reductions to monocyclic and bicyclic substructures to explore the chemical space of BIM fragments and to gain insights into the pharmacophore of BIM-type Gα_q inhibitors. Two piperazin-2-one-containing fragments and a small library of bicyclic lactams featuring fused pyrazine and diazepine rings were synthesized and evaluated. The results of a second messenger-based cellular assay indicate that the entire BIM structure is required for efficient Gα_q inhibition.

Pharmacological targeting of α-synuclein and TPPP/p25 in Parkinson's disease: challenges and opportunities in a Nutshell

With the aging of population, neurological disorders, and especially disorders involving defects in protein conformation (also known as proteopathies) pose a serious socio-economic problem. So far there is no effective treatment for most proteopathies, including Parkinson's disease (PD). The mechanism underlying PD pathogenesis is largely unknown, and the hallmark proteins, α-synuclein (SYN) and tubulin polymerization promoting protein (TPPP/p25) are challenging drug targets. These proteins are intrinsically disordered with high conformational plasticity, and have diverse physiological and pathological functions. In the healthy brain, SYN and TPPP/p25 occur in neurons and oligodendrocytes, respectively; however, in PD and multiple system atrophy, they are co-enriched and co-localized in both cell types, thereby marking pathogenesis. Although large inclusions appear at a late disease stage, small, soluble assemblies of SYN promoted by TPPP/p25 are pathogenic. In the light of these issues, we established a new innovative strategy for the validation of a specific drug target based upon the identification of contact surfaces of the pathological SYN-TPPP/p25 complex that may lead to the development of peptidomimetic foldamers suitable for pharmaceutical intervention.

Cofactor Analogues as Active Site Probes in Lysine Acetyltransferases

Lysine acetyltransferases (KATs, also termed histone acetyltransferases, HATs) catalyze the acetylation of substrate lysine residues by employing the cofactor acetyl-coenzyme A (AcCoA), thereby providing a dynamic control mechanism of protein function. Because of their major involvement in cell development and homeostasis, small-molecule modulators of KAT activity are urgently needed to assess their therapeutic potential and for probing their underlying biology. Recent advances in the field suggest that targeting the cofactor binding site represents a promising strategy for identifying potent and selective ligands. Here, we present the synthesis of two functional cofactor-based chemical probes and their usage as mechanistic tools in a broadly applicable assay platform. A fluorescence polarization (FP)-based binding assay was combined with biolayer interferometry competition analysis and a FP competition activity immunoassay to enable easy, reliable, and profound evaluation of ligands that target the KAT cofactor binding site.

Structural Insights on Fragment Binding Mode Conservation

Aiming at a deep understanding of fragment binding to ligandable targets, we performed a large scale analysis of the Protein Data Bank. Binding modes of 1832 drug-like ligands and 1079 fragments to 235 proteins were compared. We observed that the binding modes of fragments and their drug-like superstructures binding to the same protein are mostly conserved, thereby providing experimental evidence for the preservation of fragment binding modes during molecular growing. Furthermore, small chemical changes in the fragment are tolerated without alteration of the fragment binding mode. The exceptions to this observation generally involve conformational variability of the molecules. Our data analysis also suggests that, provided enough fragments have been crystallized within a protein, good interaction coverage of the binding pocket is achieved. Last, we extended our study to 126 crystallization additives and discuss in which cases they provide information relevant to structure-based drug design.

Structure-based drug design: aiming for a perfect fit

Knowledge of the three-dimensional structure of therapeutically relevant targets has informed drug discovery since the first protein structures were determined using X-ray crystallography in the 1950s and 1960s. In this editorial we provide a brief overview of the powerful impact of structure-based drug design (SBDD), which has its roots in computational and structural biology, with major contributions from both academia and industry. We describe advances in the application of SBDD for integral membrane protein targets that have traditionally proved very challenging. We emphasize the major progress made in fragment-based approaches for which success has been exemplified by over 30 clinical drug candidates and importantly three FDA-approved drugs in oncology. We summarize the articles in this issue that provide an excellent snapshot of the current state of the field of SBDD and fragment-based drug design and which offer key insights into exciting new developments, such as the X-ray free-electron laser technology, cryo-electron microscopy, open science approaches and targeted protein degradation. We stress the value of SBDD in the design of high-quality chemical tools that are used to interrogate biology and disease pathology, and to inform target validation. We emphasize the need to maintain the scientific rigour that has been traditionally associated with structural biology and extend this to other methods used in drug discovery. This is particularly important because the quality and robustness of any form of contributory data determines its usefulness in accelerating drug design, and therefore ultimately in providing patient benefit.