Chasing protons: how isothermal titration calorimetry, mutagenesis, and pKa calculations trace the locus of charge in ligand binding to a tRNA-binding enzyme

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Smart thrombosis inhibitors without bleeding side effects via charge tunable ligand design

Current treatments to prevent thrombosis, namely anticoagulants and platelets antagonists, remain complicated by the persistent risk of bleeding. Improved therapeutic strategies that diminish this risk would have a huge clinical impact. Antithrombotic agents that neutralize and inhibit polyphosphate (polyP) can be a powerful approach towards such a goal. Here, we report a design concept towards polyP inhibition, termed macromolecular polyanion inhibitors (MPI), with high binding affinity and specificity. Lead antithrombotic candidates are identified through a library screening of molecules which possess low charge density at physiological pH but which increase their charge upon binding to polyP, providing a smart way to enhance their activity and selectivity. The lead MPI candidates demonstrates antithrombotic activity in mouse models of thrombosis, does not give rise to bleeding, and is well tolerated in mice even at very high doses. The developed inhibitor is anticipated to open avenues in thrombosis prevention without bleeding risk, a challenge not addressed by current therapies.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Modeling the Structure of Human tRNA-Guanine Transglycosylase in Complex with 7-Methylguanine and Revealing the Factors that Determine the Enzyme Interaction with Inhibitors

tRNA-guanine transglycosylase, an enzyme catalyzing replacement of guanine with queuine in human tRNA and participating in the translation mechanism, is involved in the development of cancer. However, information on the small-molecule inhibitors that can suppress activity of this enzyme is very limited. Molecular dynamics simulations were used to determine the amino acid residues that provide efficient binding of inhibitors in the active site of tRNA-guanine transglycosylase. It was demonstrated using 7-methylguanine molecule as a probe that the ability of the inhibitor to adopt a charged state in the environment of hydrogen bond acceptors Asp105 and Asp159 plays a key role in complex formation. Formation of the hydrogen bonds and hydrophobic contacts with Gln202, Gly229, Phe109, and Met259 residues are also important. It has been predicted that introduction of the substituents would have a different effect on the ability to inhibit tRNA-guanine transglycosylase, as well as the DNA repair protein poly(ADP-ribose) polymerase 1, which can contribute to the development of more efficient and selective compounds.

Fragment Binding to Kinase Hinge: If Charge Distribution and Local pKa Shifts Mislead Popular Bioisosterism Concepts

Medicinal‐chemistry optimization follows strategies replacing functional groups and attaching larger substituents at a promising lead scaffold. Well‐established bioisosterism rules are considered, however, it is difficult to estimate whether the introduced modifications really match the required properties at a binding site. The electron density distribution and pK_a values are modulated influencing protonation states and bioavailability. Considering the adjacent H‐bond donor/acceptor pattern of the hinge binding motif in a kinase, we studied by crystallography a set of fragments to map the required interaction pattern. Unexpectedly, benzoic acid and benzamidine, decorated with the correct substituents, are totally bioisosteric just as carboxamide and phenolic OH. A mono‐dentate pyridine nitrogen out‐performs bi‐dentate functionalities. The importance of correctly designing pK_a values of attached functional groups by additional substituents at the parent scaffold is rendered prominent.

Fragment Screening Hit Draws Attention to a Novel Transient Pocket Adjacent to the Recognition Site of the tRNA-Modifying Enzyme TGT

Fragment-based lead discovery was applied to tRNA-guanine transglycosylase, an enzyme modifying post-transcriptionally tRNAs in Shigella, the causative agent of shigellosis. TGT inhibition prevents translation of Shigella's virulence factor VirF, hence reducing pathogenicity. One discovered fragment opens a transient subpocket in the preQ₁-recognition site by pushing back an aspartate residue. This step is associated with reorganization of further amino acids structurally transforming a loop adjacent to the recognition site by duplicating the volume of the preQ₁-recognition pocket. We synthesized 6-carboxamido-, 6-hydrazido-, and 4-guanidino-benzimidazoles to target the opened pocket, including a dihydro-imidazoquinazoline with a propyn-1-yl exit vector pointing into the transient pocket and displacing a conserved water network. MD simulations and hydration-site analysis suggest water displacement to contribute favorably to ligand binding. A cysteine residue, exclusively present in bacterial TGTs, serves as gatekeeper of the transient subpocket. It becomes accessible upon pocket opening for selective covalent attachment of electrophilic ligands in eubacterial TGTs.

Broad-scale analysis of thermodynamic signatures in medicinal chemistry: are enthalpy-favored binders the better development option?

Thermodynamic profiles of ligand binding, particularly enthalpically favored binding signatures, have been suggested as a criterion to support the decision-making process around which compounds to select for further optimization in drug development. The concept was enthusiastically taken up, but turned out to be too superficial, either because many aspects determining thermodynamic profiles are insufficiently appreciated or because it is difficult to compare such data on a global scale. The impact of water, changes in protonation states, along with buffer dependencies and incompatible measurement conditions that are far from standard conditions hamper such broad-scale comparisons. However, thermodynamic signatures can make us aware of the impact of these aspects and provide important hints for improving our understanding of the binding process and defining criteria for drug optimization.

Overview of the SAMPL6 host-guest binding affinity prediction challenge

Accurately predicting the binding affinities of small organic molecules to biological macromolecules can greatly accelerate drug discovery by reducing the number of compounds that must be synthesized to realize desired potency and selectivity goals. Unfortunately, the process of assessing the accuracy of current computational approaches to affinity prediction against binding data to biological macromolecules is frustrated by several challenges, such as slow conformational dynamics, multiple titratable groups, and the lack of high-quality blinded datasets. Over the last several SAMPL blind challenge exercises, host-guest systems have emerged as a practical and effective way to circumvent these challenges in assessing the predictive performance of current-generation quantitative modeling tools, while still providing systems capable of possessing tight binding affinities. Here, we present an overview of the SAMPL6 host-guest binding affinity prediction challenge, which featured three supramolecular hosts: octa-acid (OA), the closely related tetra-endo-methyl-octa-acid (TEMOA), and cucurbit[8]uril (CB8), along with 21 small organic guest molecules. A total of 119 entries were received from ten participating groups employing a variety of methods that spanned from electronic structure and movable type calculations in implicit solvent to alchemical and potential of mean force strategies using empirical force fields with explicit solvent models. While empirical models tended to obtain better performance than first-principle methods, it was not possible to identify a single approach that consistently provided superior results across all host-guest systems and statistical metrics. Moreover, the accuracy of the methodologies generally displayed a substantial dependence on the system considered, emphasizing the need for host diversity in blind evaluations. Several entries exploited previous experimental measurements of similar host-guest systems in an effort to improve their physical-based predictions via some manner of rudimentary machine learning; while this strategy succeeded in reducing systematic errors, it did not correspond to an improvement in statistical correlation. Comparison to previous rounds of the host-guest binding free energy challenge highlights an overall improvement in the correlation obtained by the affinity predictions for OA and TEMOA systems, but a surprising lack of improvement regarding root mean square error over the past several challenge rounds. The data suggests that further refinement of force field parameters, as well as improved treatment of chemical effects (e.g., buffer salt conditions, protonation states), may be required to further enhance predictive accuracy.

Sugar Acetonides are a Superior Motif for Addressing the Large, Solvent-Exposed Ribose-33 Pocket of tRNA-Guanine Transglycosylase

The intestinal disease shigellosis caused by Shigella bacteria affects over 120 million people annually. There is an urgent demand for new drugs as resistance against common antibiotics emerges. Bacterial tRNA-guanine transglycosylase (TGT) is a druggable target and controls the pathogenicity of Shigella flexneri. We report the synthesis of sugar-functionalized lin-benzoguanines addressing the ribose-33 pocket of TGT from Zymomonas mobilis. Ligand binding was analyzed by isothermal titration calorimetry and X-ray crystallography. Pocket occupancy was optimized by variation of size and protective groups of the sugars. The participation of a polycyclic water-cluster in the recognition of the sugar moiety was revealed. Acetonide-protected ribo- and psicofuranosyl derivatives are highly potent, benefiting from structural rigidity, good solubility, and metabolic stability. We conclude that sugar acetonides have a significant but not yet broadly recognized value in drug development.

Soaking suggests "alternative facts": Only co-crystallization discloses major ligand-induced interface rearrangements of a homodimeric tRNA-binding protein indicating a novel mode-of-inhibition

For the efficient pathogenesis of Shigella, the causative agent of bacillary dysentery, full functionality of tRNA-guanine transglycosylase (TGT) is mandatory. TGT performs post-transcriptional modifications of tRNAs in the anticodon loop taking impact on virulence development. This suggests TGT as a putative target for selective anti-shigellosis drug therapy. Since bacterial TGT is only functional as homodimer, its activity can be inhibited either by blocking its active site or by preventing dimerization. Recently, we discovered that in some crystal structures obtained by soaking the full conformational adaptation most likely induced in solution upon ligand binding is not displayed. Thus, soaked structures may be misleading and suggest irrelevant binding modes. Accordingly, we re-investigated these complexes by co-crystallization. The obtained structures revealed large conformational rearrangements not visible in the soaked complexes. They result from spatial perturbations in the ribose-34/phosphate-35 recognition pocket and, consequently, an extended loop-helix motif required to prevent access of water molecules into the dimer interface loses its geometric integrity. Thermodynamic profiles of ligand binding in solution indicate favorable entropic contributions to complex formation when large conformational adaptations in the dimer interface are involved. Native MS titration experiments reveal the extent to which the homodimer is destabilized in the presence of each inhibitor. Unexpectedly, one ligand causes a complete rearrangement of subunit packing within the homodimer, never observed in any other TGT crystal structure before. Likely, this novel twisted dimer is catalytically inactive and, therefore, suggests that stabilizing this non-productive subunit arrangement may be used as a further strategy for TGT inhibition.

Molecular imaging of fibrosis using a novel collagen-binding peptide labelled with 99mTc on SPECT/CT

Fibrosis, closely related to chronic various diseases, is a pathological process characterised by the accumulation of collagen (largely collagen type I). Non-invasive methods are necessary for the diagnosis and follow-up of fibrosis. This study aimed to develop a collagen-targeted probe for the molecular imaging of fibrosis. We identified CPKESCNLFVLKD (CBP1495) as an original collagen-binding peptide using isothermal titration calorimetry and enzyme-linked immunosorbent assay. CBP1495 effectively bound to collagen type I (K _d = 861 nM) and (GPO)₉ (K _d = 633 nM), a collagen mimetic peptide. Western blot and histochemistry validated CBP1495 targeting collagen in vitro and ex vivo. (Gly-(D)-Ala-Gly-Gly) was introduced to CBP1495 for coupling ^99mTc. Labelling efficiency of ^99mTc-CBP1495 was 95.06 ± 1.08 %. The physico-chemical properties, tracer kinetics and biodistribution of ^99mTc-CBP1495 were carried out, and showed that the peptide stably chelated ^99mTc in vitro and in vivo. SPECT/CT imaging with ^99mTc-CBP1495 was performed in rat fibrosis models, and revealed that ^99mTc-CBP1495 significantly accumulated in fibrotic lungs or livers of rats. Finally, ^99mTc-CBP1495 uptake and hydroxyproline (Hyp), a specific amino acid of collagen, were quantitatively analysed. The results demonstrated that ^99mTc-CBP1495 uptake was positvely correlated with Hyp content in lungs (P < 0.0001, r ² = 0.8266) or livers (P < 0.0001, r ² = 0.7581). Therefore, CBP1495 is a novel collagen-binding peptide, and ^99mTc-labelled CBP1495 may be a promising radiotracer for the molecular imaging of fibrosis.

Applications of isothermal titration calorimetry - the research and technical developments from 2011 to 2015

Isothermal titration calorimetry is a widely used biophysical technique for studying the formation or dissociation of molecular complexes. Over the last 5 years, much work has been published on the interpretation of isothermal titration calorimetry (ITC) data for single binding and multiple binding sites. As over 80% of ITC papers are on macromolecules of biological origin, this interpretation is challenging. Some researchers have attempted to link the thermodynamics constants to events at the molecular level. This review highlights work carried out using binding sites characterized using x-ray crystallography techniques that allow speculation about individual bond formation and the displacement of individual water molecules during ligand binding and link these events to the thermodynamic constants for binding. The review also considers research conducted with synthetic binding partners where specific binding events like anion-π and π-π interactions were studied. The revival of assays that enable both thermodynamic and kinetic information to be collected from ITC data is highlighted. Lastly, published criticism of ITC research from a physical chemistry perspective is appraised and practical advice provided for researchers unfamiliar with thermodynamics and its interpretation. Copyright © 2016 John Wiley & Sons, Ltd.

An Immucillin-Based Transition-State-Analogous Inhibitor of tRNA-Guanine Transglycosylase (TGT)

Shigellosis is one of the most severe diarrheal diseases worldwide without any efficient treatment so far. The enzyme tRNA-guanine transglycosylase (TGT) has been identified as a promising target for small-molecule drug design. Herein, we report a transition-state analogue, a small, immucillin-derived inhibitor, as a new lead structure with a novel mode of action. The complex inhibitor synthesis was accomplished in 18 steps with an overall yield of 3 %. A co-crystal structure of the inhibitor bound to Z. mobilis TGT confirmed the predicted conformation of the immucillin derivative in the enzyme active site.

Application of isothermal titration calorimetry as a tool to study natural product interactions

The study of molecular interactions of natural products by isothermal titration calorimetry (ITC) is a potent tool to get new insights of the underpinning driving forces.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.