Fragment-based discovery of 6-azaindazoles as inhibitors of bacterial DNA ligase

Fast free energy estimates from λ-dynamics with bias-updated Gibbs sampling

Relative binding free energy calculations have become an integral computational tool for lead optimization in structure-based drug design. Classical alchemical methods, including free energy perturbation or thermodynamic integration, compute relative free energy differences by transforming one molecule into another. However, these methods have high operational costs due to the need to perform many pairwise perturbations independently. To reduce costs and accelerate molecular design workflows, we present a method called λ-dynamics with bias-updated Gibbs sampling. This method uses dynamic biases to continuously sample between multiple ligand analogues collectively within a single simulation. We show that many relative binding free energies can be determined quickly with this approach without compromising accuracy. For five benchmark systems, agreement to experiment is high, with root mean square errors near or below 1.0 kcal mol-1. Free energy results are consistent with other computational approaches and within statistical noise of both methods (0.4 kcal mol-1 or less). Notably, large efficiency gains over thermodynamic integration of 18-66-fold for small perturbations and 100-200-fold for whole aromatic ring substitutions are observed. The rapid determination of relative binding free energies will enable larger chemical spaces to be more readily explored and structure-based drug design to be accelerated.

Electrostatic Complementarity in Structure-Based Drug Design

Optimization of electrostatic complementarity is an important strategy in structure-based drug discovery for improving the affinity of molecules against a specific protein target. In this Miniperspective we identify examples where deliberate optimization of protein-ligand electrostatic complementarity or intramolecular electrostatic interactions gave improvements in target affinity (up to 250-fold), physicochemical properties, in vitro properties, and off-target selectivity. We also look retrospectively at a series of factor Xa inhibitors that show an almost 8000-fold range in potency that can be correlated with the calculated electrostatic potential (ESP) surfaces. Recent developments using a graph-convolutional deep neural network to rapidly generate high quality ESP surfaces have the potential to make this useful tool more accessible for a wider audience within the field of medicinal chemistry.

Synthesis and Antiproliferative Activity of 2,4,6,7-Tetrasubstituted-2H-pyrazolo[4,3-c]pyridines

A library of 2,4,6,7-tetrasubstituted-2H-pyrazolo[4,3-c]pyridines was prepared from easily accessible 1-phenyl-3-(2-phenylethynyl)-1H-pyrazole-4-carbaldehyde via an iodine-mediated electrophilic cyclization of intermediate 4-(azidomethyl)-1-phenyl-3-(phenylethynyl)-1H-pyrazoles to 7-iodo-2,6-diphenyl-2H-pyrazolo[4,3-c]pyridines followed by Suzuki cross-couplings with various boronic acids and alkylation reactions. The compounds were evaluated for their antiproliferative activity against K562, MV4-11, and MCF-7 cancer cell lines. The most potent compounds displayed low micromolar GI50 values. 4-(2,6-Diphenyl-2H-pyrazolo[4,3-c]pyridin-7-yl)phenol proved to be the most active, induced poly(ADP-ribose) polymerase 1 (PARP-1) cleavage, activated the initiator enzyme of apoptotic cascade caspase 9, induced a fragmentation of microtubule-associated protein 1-light chain 3 (LC3), and reduced the expression levels of proliferating cell nuclear antigen (PCNA). The obtained results suggest a complex action of 4-(2,6-diphenyl-2H-pyrazolo[4,3-c]pyridin-7-yl)phenol that combines antiproliferative effects with the induction of cell death. Moreover, investigations of the fluorescence properties of the final compounds revealed 7-(4-methoxyphenyl)-2,6-diphenyl-2H-pyrazolo[4,3-c]pyridine as the most potent pH indicator that enables both fluorescence intensity-based and ratiometric pH sensing.

Compound Uptake into E. coli Can Be Facilitated by N-Alkyl Guanidiniums and Pyridiniums

Multidrug-resistant Gram-negative bacterial infections are on the rise, and with no FDA approvals for new classes of broad-spectrum antibiotics in over 50 years, these infections constitute a major threat to human health. A significant challenge is the inability of most compounds to accumulate in Gram-negative bacteria. Recently developed predictive guidelines show that appending a primary amine to an appropriately shaped compound can enhance Gram-negative accumulation. Here, we report that other positively charged nitrogen functional groups, namely, N-alkyl guanidiniums and pyridiniums, can also facilitate compound uptake into Gram-negative bacteria. The accumulation of a set of 60 nonantibiotic compounds, consisting of 20 primary amines and their corresponding guanidiniums and pyridiniums, was assessed in Escherichia coli. We also installed these alternate functional groups onto antibiotic scaffolds and assessed their accumulation and antibacterial activity in Gram-negative bacteria. The results suggest that other positively-charged, nitrogen-containing functional groups should be considered when designing antibiotics with Gram-negative activity.

Molecular Interactions of Pyrazine-Based Compounds to Proteins

Pyrazine-based compounds are of great importance in medicinal chemistry. Due to their heteroaromatic nature, they uniquely combine properties of heteroatoms (polar interactions) with the properties of aromatic moieties (nonpolar interactions). This review summarizes results of a systematic analysis of RCSB PDB database focused on important binding interactions of pyrazine-based ligands cocrystallized in protein targets. The most frequent interaction of pyrazine was hydrogen bond to pyrazine nitrogen atom as an acceptor, followed by weak hydrogen bond with pyrazine hydrogen as donor. We also identified intramolecular hydrogen bonds within pyrazine ligands, π-interactions, coordination to metal ions, and few halogen bonds in chloropyrazines. In many cases the binding mode of the pyrazine fragment was complex, involving a combination of several interactions. We conclude that pyrazine as a molecular fragment should not be perceived as a simple aromatic isostere but rather as a readily interacting moiety of drug-like molecules with high potential for interactions to proteins.

Perspective on Antibacterial Lead Identification Challenges and the Role of Hypothesis-Driven Strategies

For the past three decades, the pharmaceutical industry has undertaken many diverse approaches to discover novel antibiotics, with limited success. We have witnessed and personally experienced many mistakes, hurdles, and dead ends that have derailed projects and discouraged scientists and business leaders. Of the many factors that affect the outcomes of screening campaigns, a lack of understanding of the properties that drive efflux and permeability requirements across species has been a major barrier for advancing hits to leads. Hits that possess bacterial spectrum have seldom also possessed druglike properties required for developability and safety. Persistence in solving these two key barriers is necessary for the reinvestment into discovering antibacterial agents. This perspective narrates our experience in antibacterial discovery-our lessons learned about antibacterial challenges as well as best practices for screening strategies. One of the tenets that guides us is that drug discovery is a hypothesis-driven science. Application of this principle, at all steps in the antibacterial discovery process, should improve decision making and possibly the odds of what has become, in recent decades, an increasingly challenging endeavor with dwindling success rates.

FragPELE: Dynamic Ligand Growing within a Binding Site. A Novel Tool for Hit-To-Lead Drug Design

The early stages of drug discovery rely on hit-to-lead programs, where initial hits undergo partial optimization to improve binding affinities for their biological target. This is an expensive and time-consuming process, requiring multiple iterations of trial and error designs, an ideal scenario for applying computer simulation. However, most state-of-the-art modeling techniques fail to provide a fast and reliable answer to the Induced-Fit protein-ligand problem. To aid in this matter, we present FragPELE, a new tool for in silico hit-to-lead drug design, capable of growing a fragment from a bound core while exploring the protein-ligand conformational space. We tested the ability of FragPELE to predict crystallographic data, even in cases where cryptic sub-pockets open because of the presence of particular R-groups. Additionally, we evaluated the potential of the software on growing and scoring five congeneric series from the 2015 FEP+ dataset, comparing them to FEP+, SP and Induced-Fit Glide, and MMGBSA simulations. Results show that FragPELE could be useful not only for finding new cavities and novel binding modes in cases where standard docking tools cannot but also to rank ligand activities in a reasonable amount of time and with acceptable precision.

Antimicrobial Activity of Substituted Benzopentathiepin-6-amines

A number of substituted benzopentathiepin-6-amines and their analogues without a polysulfur ring were synthesized and evaluated in vitro for antimicrobial activity against a panel of reference bacterial and fungal strains. Trifluoroacetamide 14 demonstrated high antibacterial activity against Staphylococcus aureus (MRSA strain) with a MIC of 4 μg/mL, which was four-fold higher than the activity of a reference drug amoxicillin. This compound was also most active against the Candida albicans fungus (MIC of 1 μg ml^-1), whereas amide 17 containing a morpholine substituent was most active against the Cryptococcus neoformans fungus (MIC of 2 μg ml^-1). These compounds have no hemolytic activity and are low cytotoxic. Replacement of the pentathiepine ring with 1,3-dithiolan-2-one or 1,3-dithiolane moieties leads to loss of antimicrobial activity. Based on the QSAR analysis and molecular docking data, bacterial DNA ligase might be one of the targets for the antibacterial activity of substituted benzopentathiepin-6-amines against S. aureus.

The Macromolecular Machines that Duplicate the Escherichia coli Chromosome as Targets for Drug Discovery

DNA replication is an essential process. Although the fundamental strategies to duplicate chromosomes are similar in all free-living organisms, the enzymes of the three domains of life that perform similar functions in DNA replication differ in amino acid sequence and their three-dimensional structures. Moreover, the respective proteins generally utilize different enzymatic mechanisms. Hence, the replication proteins that are highly conserved among bacterial species are attractive targets to develop novel antibiotics as the compounds are unlikely to demonstrate off-target effects. For those proteins that differ among bacteria, compounds that are species-specific may be found. Escherichia coli has been developed as a model system to study DNA replication, serving as a benchmark for comparison. This review summarizes the functions of individual E. coli proteins, and the compounds that inhibit them.

Targeting DNA Replication and Repair for the Development of Novel Therapeutics against Tuberculosis

Mycobacterium tuberculosis is the etiological agent of tuberculosis (TB), an infectious disease which results in approximately 10 million incident cases and 1.4 million deaths globally each year, making it the leading cause of mortality from infection. An effective frontline combination chemotherapy exists for TB; however, this regimen requires the administration of four drugs in a 2 month long intensive phase followed by a continuation phase of a further 4 months with two of the original drugs, and is only effective for the treatment of drug-sensitive TB. The emergence and global spread of multidrug-resistant (MDR) as well as extensively drug-resistant (XDR) strains of M. tuberculosis, and the complications posed by co-infection with the human immunodeficiency virus (HIV) and other co-morbidities such as diabetes, have prompted urgent efforts to develop shorter regimens comprising new compounds with novel mechanisms of action. This demands that researchers re-visit cellular pathways and functions that are essential to M. tuberculosis survival and replication in the host but which are inadequately represented amongst the targets of current anti-mycobacterial agents. Here, we consider the DNA replication and repair machinery as a source of new targets for anti-TB drug development. Like most bacteria, M. tuberculosis encodes a complex array of proteins which ensure faithful and accurate replication and repair of the chromosomal DNA. Many of these are essential; so, too, are enzymes in the ancillary pathways of nucleotide biosynthesis, salvage, and re-cycling, suggesting the potential to inhibit replication and repair functions at multiple stages. To this end, we provide an update on the state of chemotherapeutic inhibition of DNA synthesis and related pathways in M. tuberculosis. Given the established links between genotoxicity and mutagenesis, we also consider the potential implications of targeting DNA metabolic pathways implicated in the development of drug resistance in M. tuberculosis, an organism which is unusual in relying exclusively on de novo mutations and chromosomal rearrangements for evolution, including the acquisition of drug resistance. In that context, we conclude by discussing the feasibility of targeting mutagenic pathways in an ancillary, “anti-evolution” strategy aimed at protecting existing and future TB drugs.

Biochemical and Structural Characterisation of DNA Ligases from Bacteria and Archaea

DNA ligases are enzymes that seal breaks in the backbones of DNA, leading to them being essential for the survival of all organisms. DNA ligases have been studied from many different types of cells and organisms and shown to have diverse sizes and sequences, with well conserved specific sequences that are required for enzymatic activity. A significant number of DNA ligases have been isolated or prepared in recombinant forms and, here, we review their biochemical and structural characterisation. All DNA ligases contain an essential lysine that transfers an adenylate group from a co-factor to the 5'-phosphate of the DNA end that will ultimately be joined to the 3'-hydroxyl of the neighbouring DNA strand. The essential DNA ligases in bacteria use nicotinamide adenine dinucleotide ( β -NAD+) as their co-factor whereas those that are essential in other cells use adenosine-5'-triphosphate (ATP) as their co-factor. This observation suggests that the essential bacterial enzyme could be targeted by novel antibiotics and the complex molecular structure of β -NAD+ affords multiple opportunities for chemical modification. Several recent studies have synthesised novel derivatives and their biological activity against a range of DNA ligases has been evaluated as inhibitors for drug discovery and/or non-natural substrates for biochemical applications. Here, we review the recent advances that herald new opportunities to alter the biochemical activities of these important enzymes. The recent development of modified derivatives of nucleotides highlights that the continued combination of structural, biochemical and biophysical techniques will be useful in targeting these essential cellular enzymes.

The Search for Herbal Antibiotics: An In-Silico Investigation of Antibacterial Phytochemicals

Recently, the emergence and spread of pathogenic bacterial resistance to many antibiotics (multidrug-resistant strains) have been increasing throughout the world. This phenomenon is of great concern and there is a need to find alternative chemotherapeutic agents to combat these antibiotic-resistant microorganisms. Higher plants may serve as a resource for new antimicrobials to replace or augment current therapeutic options. In this work, we have carried out a molecular docking study of a total of 561 antibacterial phytochemicals listed in the Dictionary of Natural Products, including 77 alkaloids (17 indole alkaloids, 27 isoquinoline alkaloids, 4 steroidal alkaloids, and 28 miscellaneous alkaloids), 99 terpenoids (5 monoterpenoids, 31 sesquiterpenoids, 52 diterpenoids, and 11 triterpenoids), 309 polyphenolics (87 flavonoids, 25 chalcones, 41 isoflavonoids, 5 neoflavonoids, 12 pterocarpans, 10 chromones, 7 condensed tannins, 11 coumarins, 30 stilbenoids, 2 lignans, 5 phenylpropanoids, 13 xanthones, 5 hydrolyzable tannins, and 56 miscellaneous phenolics), 30 quinones, and 46 miscellaneous phytochemicals, with six bacterial protein targets (peptide deformylase, DNA gyrase/topoisomerase IV, UDP-galactose mutase, protein tyrosine phosphatase, cytochrome P450 CYP121, and NAD⁺-dependent DNA ligase). In addition, 35 known inhibitors were docked with their respective targets for comparison purposes. Prenylated polyphenolics showed the best docking profiles, while terpenoids had the poorest. The most susceptible protein targets were peptide deformylases and NAD⁺-dependent DNA ligases.

Design, synthesis and anticancer activity of dihydropyrimidinone–semicarbazone hybrids as potential human DNA ligase 1 inhibitors

A series of rationally designed new class of hLig1 inhibitors with potentin vitroanti-cancer properties is presented.

Detection of secondary binding sites in proteins using fragment screening

Proteins need to be tightly regulated as they control biological processes in most normal cellular functions. The precise mechanisms of regulation are rarely completely understood but can involve binding of endogenous ligands and/or partner proteins at specific locations on a protein that can modulate function. Often, these additional secondary binding sites appear separate to the primary binding site, which, for example for an enzyme, may bind a substrate. In previous work, we have uncovered several examples in which secondary binding sites were discovered on proteins using fragment screening approaches. In each case, we were able to establish that the newly identified secondary binding site was biologically relevant as it was able to modulate function by the binding of a small molecule. In this study, we investigate how often secondary binding sites are located on proteins by analyzing 24 protein targets for which we have performed a fragment screen using X-ray crystallography. Our analysis shows that, surprisingly, the majority of proteins contain secondary binding sites based on their ability to bind fragments. Furthermore, sequence analysis of these previously unknown sites indicate high conservation, which suggests that they may have a biological function, perhaps via an allosteric mechanism. Comparing the physicochemical properties of the secondary sites with known primary ligand binding sites also shows broad similarities indicating that many of the secondary sites may be druggable in nature with small molecules that could provide new opportunities to modulate potential therapeutic targets.

Accurate Binding Free Energy Predictions in Fragment Optimization

Predicting protein-ligand binding free energies is a central aim of computational structure-based drug design (SBDD)--improved accuracy in binding free energy predictions could significantly reduce costs and accelerate project timelines in lead discovery and optimization. The recent development and validation of advanced free energy calculation methods represents a major step toward this goal. Accurately predicting the relative binding free energy changes of modifications to ligands is especially valuable in the field of fragment-based drug design, since fragment screens tend to deliver initial hits of low binding affinity that require multiple rounds of synthesis to gain the requisite potency for a project. In this study, we show that a free energy perturbation protocol, FEP+, which was previously validated on drug-like lead compounds, is suitable for the calculation of relative binding strengths of fragment-sized compounds as well. We study several pharmaceutically relevant targets with a total of more than 90 fragments and find that the FEP+ methodology, which uses explicit solvent molecular dynamics and physics-based scoring with no parameters adjusted, can accurately predict relative fragment binding affinities. The calculations afford R(2)-values on average greater than 0.5 compared to experimental data and RMS errors of ca. 1.1 kcal/mol overall, demonstrating significant improvements over the docking and MM-GBSA methods tested in this work and indicating that FEP+ has the requisite predictive power to impact fragment-based affinity optimization projects.

Base-modified NAD and AMP derivatives and their activity against bacterial DNA ligases

We report the chemical synthesis and conformational analysis of a collection of 2-, 6- and 8-substituted derivatives of β-NAD(+) and AMP, and their biochemical evaluation against NAD(+)-dependent DNA ligases from Escherichia coli and Mycobacterium tuberculosis. Bacterial DNA ligases are validated anti-microbial targets, and new strategies for their inhibition are therefore of considerable scientific and practical interest. Our study includes several pairs of β-NAD(+) and AMP derivatives with the same substitution pattern at the adenine base. This has enabled the first direct comparison of co-substrate and inhibitor behaviour against bacterial DNA ligases. Our results suggest that an additional substituent in position 6 or 8 of the adenine base in β-NAD(+) is detrimental for activity as either co-substrate or inhibitor. In contrast, substituents in position 2 are not only tolerated, but appear to give rise to a new mode of inhibition, which targets the conformational changes these DNA ligases undergo during catalysis. Using a molecular modelling approach, we highlight that these findings have important implications for our understanding of ligase mechanism and inhibition, and may provide a promising starting point for the rational design of a new class of inhibitors against NAD(+)-dependent DNA ligases.