Dynamic Profiling and Binding Affinity Prediction of NBTI Antibacterials against DNA Gyrase Enzyme by Multidimensional Machine Learning and Molecular Dynamics Simulations.

Bacterial type II topoisomerases are well-characterized and clinically important targets for antibacterial chemotherapy. Novel bacterial topoisomerase inhibitors (NBTIs) are a newly disclosed class of antibacterials. Prediction of their binding affinity to these enzymes would be beneficial for de novo design/optimization of new NBTIs. Utilizing in vitro NBTI experimental data, we constructed two comprehensive multidimensional DNA gyrase surrogate models for Staphylococcus aureus (q2 = 0.791) and Escherichia coli (q2 = 0.806). Both models accurately predicted the IC50s of 26 NBTIs from our recent studies. To investigate the NBTI's dynamic profile and binding to both targets, 10 selected NBTIs underwent molecular dynamics (MD) simulations. The analysis of MD production trajectories confirmed key hydrogen-bonding and hydrophobic contacts that NBTIs establish in both enzymes. Moreover, the binding free energies of selected NBTIs were computed by the linear interaction energy (LIE) method employing an in-house derived set of fitting parameters (α = 0.16, ß = 0.029, γ = 0.0, and intercept = -1.72), which are successfully applicable to DNA gyrase of Gram-positive/Gram-negative pathogens. Both methods offer accurate predictions of the binding free energies of NBTIs against S. aureus and E. coli DNA gyrase. We are confident that this integrated modeling approach could be valuable in the de novo design and optimization of efficient NBTIs for combating resistant bacterial pathogens.

The overview of development of novel bacterial topoisomerase inhibitors effective against multidrug-resistant bacteria in an academic environment: From early hits to in vivo active antibacterials.

Antimicrobial resistance caused by the excessive and inappropriate use of antibacterial drugs is a global health concern. Currently, we are walking a fine line between the fact that most bacterial infections can still be cured with the antibiotics known so far, and the emergence of infections with bacteria resistant to several drugs at the same time, against which we no longer have an effective drug. Therefore, new antibacterial drugs are urgently needed to curb the hard-to-treat infections. Our group has developed new antibacterials from the class of novel bacterial topoisomerase inhibitors (NBTIs) that exhibit broad-spectrum antibacterial activity. This article reviews our efforts in developing highly potent NBTIs over the past decade. Following the discovery of an initial hit with potent enzyme inhibitory and broad-spectrum antibacterial activity, an extensive hit-to-lead campaign was conducted with the goal of optimizing physicochemical properties, reducing hERG inhibition, and maintaining antibacterial activity against both Gram-positive and Gram-negative bacteria, with a focus on methicillin-resistant Staphylococcus aureus (MRSA). This optimization strategy resulted in an amide-containing, focused NBTI library with compounds exhibiting potent antibacterial activity against Gram-positive bacteria, reduced hERG inhibition, no cardiotoxicity in in vivo zebrafish model, and favorable in vivo efficacy in a neutropenic murine thigh infection model for MRSA infections.

Exploring Alternative Pathways to Target Bacterial Type II Topoisomerases Using NBTI Antibacterials: Beyond Halogen-Bonding Interactions.

Novel bacterial topoisomerase inhibitors (NBTIs) are a new class of antibacterial agents that target bacterial type II topoisomerases (DNA gyrase and topoisomerase IV). Our recently disclosed crystal structure of an NBTI ligand in complex with DNA gyrase and DNA revealed that the halogen atom in the para position of the phenyl right hand side (RHS) moiety is able to establish strong symmetrical bifurcated halogen bonds with the enzyme; these are responsible for the excellent enzyme inhibitory potency and antibacterial activity of these NBTIs. To further assess the possibility of any alternative interactions (e.g., hydrogen-bonding and/or hydrophobic interactions), we introduced various non-halogen groups at the p-position of the phenyl RHS moiety. Considering the hydrophobic nature of amino acid residues delineating the NBTI's binding pocket in bacterial topoisomerases, we demonstrated that designed NBTIs cannot establish any hydrogen-bonding interactions with the enzyme; hydrophobic interactions are feasible in all respects, while halogen-bonding interactions are apparently the most preferred.

Amide containing NBTI antibacterials with reduced hERG inhibition, retained antimicrobial activity against gram-positive bacteria and in vivo efficacy.

Novel bacterial topoisomerase inhibitors (NBTIs) are new promising antimicrobials for the treatment of multidrug-resistant bacterial infections. In recent years, many new NBTIs have been discovered, however most of them struggle with the same issue - the balance between antibacterial activity and hERG-related toxicity. We started a new campaign by optimizing the previous series of NBTIs, followed by the design and synthesis of a new, amide-containing focused NBTI library to reduce hERG inhibition and maintain antibacterial activity against Gram-positive bacteria, including methicillin-resistant Staphylococcus aureus (MRSA). This optimization strategy yielded the lead compound 12 that exhibits potent antibacterial activity against Gram-positive bacteria, reduced hERG inhibition, no cardiotoxicity in zebrafish model, and a favorable in vivo efficacy in a neutropenic murine thigh infection model of MRSA infection.

Replica-exchange optimization of antibody fragments.

In the framework of the rational design of macromolecules capable of binding to a specific target for biosensing applications, we here further develop an evolutionary protocol designed to optimize the binding affinity of protein binders. In particular we focus on the optimization of the binding portion of small antibody fragments known as nanobodies (or VHH) and choose the hen egg white lysozyme (HEWL) as our target. By implementing a replica exchange scheme for this optimization, we show that an initial hit is not needed and similar solutions can be found by either optimizing an already known anti-HEWL VHH or a randomly selected binder (here a VHH selective towards another macromolecule). While we believe that exhaustive searches of the mutation space are most appropriate when only few key residues have to be optimized, in case a lead binder is not available the proposed evolutionary algorithm should be instead the method of choice.

Diminishing hERG inhibitory activity of aminopiperidine-naphthyridine linked NBTI antibacterials by structural and physicochemical optimizations.

Novel bacterial topoisomerase inhibitors (NBTIs) are an important new class of antibacterials targeting bacterial type II topoisomerases (DNA gyrase and topoisomerase IV). Notwithstanding their potent antibacterial activity, they suffer from a detrimental class-related hERG blockage. In this study, we designed and synthesized an optimized library of NBTIs comprising different linker moieties that exhibit reduced hERG inhibition and retain inhibitory potencies on DNA gyrase and topoisomerase IV of Staphylococcus aureus and Escherichia coli, respectively, as well as potent antibacterial activities. Substitution of the linker's tertiary amine with polar groups outcome in diminished hERG inhibition. Compound 17 expresses nanomolar enzyme inhibitory potency and antibacterial activity against both Gram-positive and Gram-negative bacteria as well as reduced hERG inhibition relative to our previously published NBTI analogs. Here, we point to some important NBTI's structural features that influence their hERG inhibitory activity.

The Structural Features of Novel Bacterial Topoisomerase Inhibitors That Define Their Activity on Topoisomerase IV.

The continued emergence of bacterial resistance has created an urgent need for new and effective antibacterial agents. Bacterial type II topoisomerases, such as DNA gyrase and topoisomerase IV (topoIV), are well-validated targets for antibacterial chemotherapy. The novel bacterial topoisomerase inhibitors (NBTIs) represent one of the new promising classes of antibacterial agents. They can inhibit both of these bacterial targets; however, their potencies differ on the targets among species, making topoIV probably a primary target of NBTIs in Gram-negative bacteria. Therefore, it is important to gain an insight into the NBTIs key structural features that govern the topoIV inhibition. However, in Gram-positive bacteria, topoIV is also a significant target for achieving dual-targeting, which in turn contributes to avoiding bacterial resistance caused by single-target mutations. In this perspective, we address the structure-activity relationship guidelines for NBTIs that target the topoIV enzyme in Gram-positive and Gram-negative bacteria.

Computational Evolution Protocol for Peptide Design.

Computational peptide design is useful for therapeutics, diagnostics, and vaccine development. To select the most promising peptide candidates, the key is describing accurately the peptide-target interactions at the molecular level. We here review a computational peptide design protocol whose key feature is the use of all-atom explicit solvent molecular dynamics for describing the different peptide-target complexes explored during the optimization. We describe the milestones behind the development of this protocol, which is now implemented in an open-source code called PARCE. We provide a basic tutorial to run the code for an antibody fragment design example. Finally, we describe three additional applications of the method to design peptides for different targets, illustrating the broad scope of the proposed approach.

Structurally Optimized Potent Dual-Targeting NBTI Antibacterials with an Enhanced Bifurcated Halogen-Bonding Propensity.

We designed and synthesized an optimized library of novel bacterial topoisomerase inhibitors with p-halogenated phenyl right-hand side fragments and significantly enhanced and balanced dual-targeted DNA gyrase and topoisomerase IV activities of Staphylococcus aureus and Escherichia coli. By increasing the electron-withdrawing properties of the p-halogenated phenyl right-hand side fragment and maintaining a similar lipophilicity and size, an increased potency was achieved, indicating that the antibacterial activities of this series of novel bacterial topoisomerase inhibitors against all target enzymes are determined by halogen-bonding rather than van der Waals interactions. They show nanomolar enzyme inhibitory and whole-cell antibacterial activities against S. aureus and methicillin-resistant S. aureus (MRSA) strains. However, due to the relatively high substrate specificity for the bacterial efflux pumps, they tend to be less potent against E. coli and other Gram-negative pathogens.

A Fine-Tuned Lipophilicity/Hydrophilicity Ratio Governs Antibacterial Potency and Selectivity of Bifurcated Halogen Bond-Forming NBTIs.

Herein, we report the design of a focused library of novel bacterial topoisomerase inhibitors (NBTIs) based on innovative mainly monocyclic right-hand side fragments active against DNA gyrase and Topo IV. They exhibit a very potent and wide range of antibacterial activity, even against some of the most concerning hard-to-treat pathogens for which new antibacterials are urgently needed, as reported by the WHO and CDC. NBTIs enzyme activity and whole cell potency seems to depend on the fine-tuned lipophilicity/hydrophilicity ratio that governs the permeability of those compounds through the bacterial membranes. Lipophilicity of NBTIs is apparently optimal for passing through the membrane of Gram-positive bacteria, but the higher, although not excessive lipophilicity and suitable hydrophilicity seems to determine the passage through Gram-negative bacterial membranes. However, due to the considerable hERG inhibition, which is still at least two orders of magnitude away from MICs, continued optimization is required to realize their full potential.

Potent DNA gyrase inhibitors bind asymmetrically to their target using symmetrical bifurcated halogen bonds.

Novel bacterial type II topoisomerase inhibitors (NBTIs) stabilize single-strand DNA cleavage breaks by DNA gyrase but their exact mechanism of action has remained hypothetical until now. We have designed a small library of NBTIs with an improved DNA gyrase-binding moiety resulting in low nanomolar inhibition and very potent antibacterial activity. They stabilize single-stranded cleavage complexes and, importantly, we have obtained the crystal structure where an NBTI binds gyrase-DNA in a single conformation lacking apparent static disorder. This directly proves the previously postulated NBTI mechanism of action and shows that they stabilize single-strand cleavage through asymmetric intercalation with a shift of the scissile phosphate. This crystal stucture shows that the chlorine forms a halogen bond with the backbone carbonyls of the two symmetry-related Ala68 residues. To the best of our knowledge, such a so-called symmetrical bifurcated halogen bond has not been identified in a biological system until now.

Structure-Function Relationships in the Human P-Glycoprotein (ABCB1): Insights from Molecular Dynamics Simulations.

P-Glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette superfamily of transporters, and it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping compounds out of cells. P-gp contributes to a reduction in toxicity, and has broad substrate specificity. It is involved in the failure of many cancer and antiviral chemotherapies due to the phenomenon of multidrug resistance (MDR), in which the membrane transporter removes chemotherapeutic drugs from target cells. Understanding the details of the ligand-P-gp interaction is therefore critical for the development of drugs that can overcome the MDR phenomenon, for the early identification of P-gp substrates that will help us to obtain a more effective prediction of toxicity, and for the subsequent outdesign of substrate properties if needed. In this work, a series of molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the effects of binding different compounds on the conformational dynamics of P-gp. The results revealed significant differences in the behaviour of P-gp in the presence of active and non-active compounds within the binding pocket, as different patterns of movement were identified that could be correlated with conformational changes leading to the activation of the translocation mechanism. The predicted ligand-P-gp interactions are in good agreement with the available experimental data, as well as the estimation of the binding-free energies of the studied complexes, demonstrating the validity of the results derived from the MD simulations.

Homology Modeling of the Human P-glycoprotein (ABCB1) and Insights into Ligand Binding through Molecular Docking Studies.

The ABCB1 transporter also known as P-glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette super-family of transporters; it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping the compounds out of cells. P-gp contributes to a decrease of toxicity and possesses broad substrate specificity. It is involved in the failure of numerous anticancer and antiviral chemotherapies due to the multidrug resistance (MDR) phenomenon, where it removes the chemotherapeutics out of the targeted cells. Understanding the details of the ligand-P-gp interaction is therefore crucial for the development of drugs that might overcome the MRD phenomenon and for obtaining a more effective prediction of the toxicity of certain compounds. In this work, an in silico modeling was performed using homology modeling and molecular docking methods with the aim of better understanding the ligand-P-gp interactions. Based on different mouse P-gp structural templates from the PDB repository, a 3D model of the human P-gp (hP-gp) was constructed by means of protein homology modeling. The homology model was then used to perform molecular docking calculations on a set of thirteen compounds, including some well-known compounds that interact with P-gp as substrates, inhibitors, or both. The sum of ranking differences (SRD) was employed for the comparison of the different scoring functions used in the docking calculations. A consensus-ranking scheme was employed for the selection of the top-ranked pose for each docked ligand. The docking results showed that a high number of π interactions, mainly π-sigma, π-alkyl, and π-π type of interactions, together with the simultaneous presence of hydrogen bond interactions contribute to the stability of the ligand-protein complex in the binding site. It was also observed that some interacting residues in hP-gp are the same when compared to those observed in a co-crystallized ligand (PBDE-100) with mouse P-gp (PDB ID: 4XWK). Our in silico approach is consistent with available experimental results regarding P-gp efflux transport assay; therefore it could be useful in the prediction of the role of new compounds in systemic toxicity.

Two Decades of Successful SAR-Grounded Stories of the Novel Bacterial Topoisomerase Inhibitors (NBTIs).

The emergence of bacterial resistance against life-saving medicines has forced the scientific community and pharmaceutical industry to take actions in the quest for novel antibacterials. These should not only overcome the existing bacterial resistance but also provide at least interim effective protection against emerging bacterial infections. Research into DNA gyrase and topoisomerase IV inhibitors has become a particular focus, with the description of a new class of bacterial topoisomerase type II inhibitors known as "novel bacterial topoisomerase inhibitors", NBTIs. Elucidation of the key structural modifications incorporated into these inhibitors and the impact these can have on their general physicochemical properties are detailed in this review. This defines novel bacterial topoisomerase inhibitors with promising antibacterial activities and potencies, which thus represent one potential example of the future "drugs for bad bugs", as identified by the World Health Organization.

Multiclass Classifier for P-Glycoprotein Substrates, Inhibitors, and Non-Active Compounds.

P-glycoprotein (P-gp) is a transmembrane protein that actively transports a wide variety of chemically diverse compounds out of the cell. It is highly associated with the ADMET (absorption, distribution, metabolism, excretion and toxicity) properties of drugs/drug candidates and contributes to decreasing toxicity by eliminating compounds from cells, thereby preventing intracellular accumulation. Therefore, in the drug discovery and toxicological assessment process it is advisable to pay attention to whether a compound under development could be transported by P-gp or not. In this study, an in silico multiclass classification model capable of predicting the probability of a compound to interact with P-gp was developed using a counter-propagation artificial neural network (CP ANN) based on a set of 2D molecular descriptors, as well as an extensive dataset of 2512 compounds (1178 P-gp inhibitors, 477 P-gp substrates and 857 P-gp non-active compounds). The model provided a good classification performance, producing non error rate (NER) values of 0.93 for the training set and 0.85 for the test set, while the average precision (AvPr) was 0.93 for the training set and 0.87 for the test set. An external validation set of 385 compounds was used to challenge the model's performance. On the external validation set the NER and AvPr values were 0.70 for both indices. We believe that this in silico classifier could be effectively used as a reliable virtual screening tool for identifying potential P-gp ligands.

Cyclohexyl amide-based novel bacterial topoisomerase inhibitors with prospective GyrA-binding fragments.

Aim: Novel bacterial topoisomerase inhibitors (NBTIs) are a promising class of bacterial topoisomerase II inhibitors that are gaining more and more importance mainly because of their excellent antibacterial activity, as well as their lack of cross-resistance to quinolones. Results: Described here is the synthesis and biological evaluation of a tiny series of new virtually assembled NBTIs containing synthetically feasible right-hand side fragments capable of binding the GyrA subunit of the bacterial DNA gyrase-DNA complex. Conclusion: NBTI variants with incorporated 1-phenylpyrazole right-hand side moiety show suitable antibacterial activity against Gram-positive Staphylococcus aureus, with confirmed selectivity over the human topoisomerase IIα enzyme.

Revisiting fish toxicity of active pharmaceutical ingredients: Mechanistic insights from integrated ligand-/structure-based assessments on acetylcholinesterase.

The release of active pharmaceutical ingredients (APIs) into the environment is of great concern for aquatic ecosystem as many of these chemicals are designed to exert biological activity. Hence, their impact on non-target organisms like fish would not be surprising. In this respect, we revisited fish toxicity data of pharmaceuticals to generate linear and non-linear quantitative structure-toxicity relationships (QSTRs). We predicted fish lethality data from the validated QSTR models for 120 APIs with no experimental fish toxicity data. Toxicity of APIs on aquatic organisms is not fully characterized. Therefore, to provide a mechanistic insight for the assessment of API's toxicity to fish, the outcome of the derived QSTR models was integrated with structure-based toxicophore and molecular docking studies, utilizing the biomarker enzyme acetylcholinesterase originating from fish Torpedo californica (TcAChE). Toxicophore virtual screening of 60 chemicals with pT > 0 identified 23 hits as potential TcAChE binders with binding free energies ranging from -6.5 to -12.9 kcal/mol. The TcAChE-ligand interaction analysis revealed a good nesting of all 23 hits within TcAChE binding site through establishing strong lipophilic and hydrogen bonding interactions with the surrounding key amino acid residues. Among the chemicals passing the criteria of our integrated approach, majority of APIs belong noticeably to the Central Nervous System class. The screened chemicals displayed not only comprehensive toxicophore coverage, but also strong binding affinities according to the docking calculations, mainly due to interactions with TcAChE's key amino acid residues Tyr121, Tyr130, Tyr334, Trp84, Phe290, Phe330, Phe331, Ser122, and Ser200. Moreover, we propose here that binding of pharmaceuticals to AChE might have a potential in triggering molecular initiating events for adverse outcome pathways (AOPs), which in turn can play an important role for future screening of APIs lacking fish lethality data.

Insight into structural requirements for selective and/or dual CXCR3 and CXCR4 allosteric modulators.

Based on the previously published pyrazolopyridine-based hit compound for which negative allosteric modulation of both CXCR3 and CXCR4 receptors was disclosed, we designed, synthesized and biologically evaluated a set of novel, not only negative, but also positive allosteric modulators with preserved pyrazolopyridine core. Compound 9e is a dual negative modulator, inhibiting G protein activity of both receptors. For CXCR4 receptor para-substituted aromatic group of compounds distinguishes between negative and positive modulation. Para-methoxy substitution leads to functional antagonism, while para-chloro triggers agonism. Additionally, we discovered that chemotaxis is not completely correlated with G protein pathways. This is the first work in which we have on a series of compounds successfully demonstrated that it is possible to produce selective as well as dual-acting modulators of chemokine receptors, which is very promising for future research in the field of discovery of selective or dual modulators of chemokine receptors.

Structure-based design of novel combinatorially generated NBTIs as potential DNA gyrase inhibitors against various Staphylococcus aureus mutant strains.

Although intercalating agents such as quinolones have had proven therapeutic success as antibacterial agents for more than 40 years, new forms of quinolone-based resistance in bacteria are continually emerging. To alleviate this problem, a new class of antibacterials is urgently needed; recently, novel bacterial topoisomerase inhibitors (NBTIs) have been found to be particularly important. Based on 67 experimentally evaluated NBTIs against wild-type (WT) DNA gyrase originating from Staphylococcus aureus, a predictive QSAR model was initially constructed and validated and was later used for in silico prediction of biological activities for an in house designed compound library of 548 novel drug-like NBTI combinatorial analogs. To evaluate the influence of gyrA alterations on NBTI resistance, various mutant homology models were constructed; meanwhile, their resistance profiles were assessed and validated relative to that of WT enzyme by structure-based virtual screening (VS) of known NBTIs. Surprisingly, the M121K mutant model was recognized as the most selective due to an additional established cation-π interaction between K121-NH3+ (not found in the WT) and the aromatic moiety of the NBTI right-hand site (RHS) fragment; this finding was additionally supported by VS of our combinatorially generated NBTIs. Moreover, we identified several attractive, synthetically feasible RHS building blocks that may enable the development of new NBTIs.