Kinetic and Structural Insights into the Mechanism of Binding of Sulfonamides to Human Carbonic Anhydrase by Computational and Experimental Studies

Computational approaches for the identification of novel metal-binding pharmacophores: advances and challenges

Metalloenzymes are important therapeutic targets for a variety of human diseases. Computational approaches have recently emerged as effective tools to understand metal-ligand interactions and expand the structural diversity of both metalloenzyme inhibitors (MIs) and metal-binding pharmacophores (MBPs). In this review, we highlight key advances in currently available fine-tuning modeling methods and data-driven cheminformatic approaches. We also discuss major challenges to the recognition of novel MBPs and MIs. The evidence provided herein could expedite future computational efforts to guide metalloenzyme-based drug discovery.

Multiple Binding Modes of Inhibitors to Human Carbonic Anhydrases: An Update on the Design of Isoform-Specific Modulators of Activity

Human carbonic anhydrases (hCAs) are widespread zinc enzymes that catalyze the hydration of CO2 to bicarbonate and a proton. Currently, 15 isoforms have been identified, of which only 12 are catalytically active. Given their involvement in numerous physiological and pathological processes, hCAs are recognized therapeutic targets for the development of inhibitors with biomedical applications. However, despite massive development efforts, very few of the presently available hCA inhibitors show selectivity for a specific isoform. X-ray crystallography is a very useful tool for the rational drug design of enzyme inhibitors. In 2012 we published in Chemical Reviews a highly cited review on hCA family (Alterio, V. et al. Chem Rev. 2012, 112, 4421-4468), analyzing about 300 crystallographic structures of hCA/inhibitor complexes and describing the different CA inhibition mechanisms existing up to that date. However, in the period 2012-2023, almost 700 new hCA/inhibitor complex structures have been deposited in the PDB and a large number of new inhibitor classes have been discovered. Based on these considerations, the aim of this Review is to give a comprehensive update of the structural aspects of hCA/inhibitor interactions covering the period 2012-2023 and to recapitulate how this information can be used for the rational design of more selective versions of such inhibitors.

Design and Synthesis of 3-(Phenylsulfonamido)benzamide Derivatives as Potent Carbonic Anhydrase IX Inhibitors: Biological Evaluations and Molecular Modeling Studies

INTRODUCTION

Carbonic anhydrase IX (CAIX) is known to be overexpressed in various tumors and plays a significant role in tumor development and progression.

METHODS

A series of 3-(benzylsulfonamido)benzamides derivatives was synthesized and tested for their CAIX inhibitory activities. The two most active compounds were subjected to cytotoxicity testing against a panel of 60 cancer cell lines.

RESULTS

Many of the synthesized compounds successfully inhibited CAIX activities, exhibiting IC50 values in the low nanomolar range. The most potent CAIX inhibitor was compound 14, with an IC50 of 140 nM. Structure-activity relationship analysis of the synthesized compounds supported with molecular docking revealed strong coordination of sulfonamide moiety with the catalytic Zn2+ metal, hydrophobic interactions of the benzylsulfonamido ring with a hydrophobic pocket, and π- stacking interactions of the aryl ring with an aromatic surface. The two most active analogues (10 and 14) were further tested for their antiproliferative activities in the NCI-60 human tumor cell lines. Notably, compound 14 demonstrated potent growth inhibitory effects against several cancer cell lines.

CONCLUSION

The synthesized analogues represent a novel scaffold for the treatment of different types of cancer by targeting CAIX.

Effect of Water Networks On Ligand Binding: Computational Predictions vs Experiments

Rational drug design focuses on the explanation and prediction of complex formation between therapeutic targets and small-molecule ligands. As a third and often overlooked interacting partner, water molecules play a critical role in the thermodynamics of protein-ligand binding, impacting both the entropy and enthalpy components of the binding free energy and by extension, on-target affinity and bioactivity. The community has realized the importance of binding site waters, as evidenced by the number of computational tools to predict the structure and thermodynamics of their networks. However, quantitative experimental characterization of relevant protein-ligand-water systems, and consequently the validation of these modeling methods, remains challenging. Here, we investigated the impact of solvent exchange from light (H2O) to heavy water (D2O) to provide complete thermodynamic profiling of these ternary systems. Utilizing the solvent isotope effects, we gain a deeper understanding of the energetic contributions of various components. Specifically, we conducted isothermal titration calorimetry experiments on trypsin with a series of p-substituted benzamidines, as well as carbonic anhydrase II (CAII) with a series of aromatic sulfonamides. Significant differences in binding enthalpies found between light vs heavy water indicate a substantial role of the binding site water network in protein-ligand binding. Next, we challenged two conceptually distinct modeling methods, the grid-based WaterFLAP and the molecular dynamics-based MobyWat, by predicting and scoring relevant water networks. The predicted water positions accurately reproduce those in available high-resolution X-ray and neutron diffraction structures of the relevant protein-ligand complexes. Estimated energetic contributions of the identified water networks were corroborated by the experimental thermodynamics data. Besides providing a direct validation for the predictive power of these methods, our findings confirmed the importance of considering binding site water networks in computational ligand design.

Effects of Phosphorylation on the Activity, Inhibition and Stability of Carbonic Anhydrases

Carbonic anhydrases (CAs) are a metalloenzyme family that have important roles in cellular processes including pH homeostasis and have been implicated in multiple pathological conditions. Small molecule inhibitors have been developed to target carbonic anhydrases, but the effects of post-translational modifications (PTMs) on the activity and inhibition profiles of these enzymes remain unclear. Here, we investigate the effects of phosphorylation, the most prevalent carbonic anhydrase PTM, on the activities and drug-binding affinities of human CAI and CAII, two heavily modified active isozymes. Using serine to glutamic acid (S > E) mutations to mimic the effect of phosphorylation, we demonstrate that phosphomimics at a single site can significantly increase or decrease the catalytic efficiencies of CAs, depending on both the position of the modification and the CA isoform. We also show that the S > E mutation at Ser50 of hCAII decreases the binding affinities of hCAII with well-characterized sulphonamide inhibitors including by over 800-fold for acetazolamide. Our findings suggest that CA phosphorylation may serve as a regulatory mechanism for enzymatic activity, and affect the binding affinity and specificity of small, drug and drug-like molecules. This work should motivate future studies examining the PTM-modification forms of CAs and their distributions, which should provide insights into CA physiopathological functions and facilitate the development of 'modform-specific' carbonic anhydrase inhibitors.

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.

In the footsteps of an inhibitor unbinding from the active site of human carbonic anhydrase II

The crystal structure of human carbonic anhydrase (HCA) II bound to an inhibitor molecule, 6-hydroxy-2-thioxocoumarin (FC5), shows FC5 to be located in a hydrophobic pocket at the active site. The present work employs classical molecular dynamics (MD) simulation to follow the FC5 molecule for 1 μs as it unbinds from its binding location, adopts the path of substrate/product diffusion (path 1) to leave the active site at around 75 ns. It is then found to undergo repeated binding and unbinding at different locations on the surface of the enzyme in water. Several transient excursions through different regions of the enzyme are also observed prior to its exit from the active site. These transient paths are combined with functionally relevant cavities/channels to enlist five additional pathways (path 2-6). Pathways 1-6 are subsequently explored using steered MD and umbrella sampling simulations. A free energy barrier of 0.969 kcal mol^-1 is encountered along path 1, while barriers in the range of 0.57-2.84 kcal mol^-1 are obtained along paths 2, 4 and 5. We also analyze in detail the interaction between FC5 and the enzyme along each path as the former leaves the active site of HCA II. Our results indicate path 1 to be the major exit pathway for FC5, although competing contributions may also come from the paths 2, 4 and 5.Communicated by Ramaswamy H. Sarma.

Reconsidering anion inhibitors in the general context of drug design studies of modulators of activity of the classical enzyme carbonic anhydrase

Inorganic anions inhibit the metalloenzyme carbonic anhydrase (CA, EC 4.2.1.1) generally by coordinating to the active site metal ion. Cyanate was reported as a non-coordinating CA inhibitor but those erroneous results were subsequently corrected by another group. We review the anion CA inhibitors (CAIs) in the more general context of drug design studies and the discovery of a large number of inhibitor classes and inhibition mechanisms, including zinc binders (sulphonamides and isosteres, dithiocabamates and isosteres, thiols, selenols, benzoxaboroles, ninhydrins, etc.); inhibitors anchoring to the zinc-coordinated water molecule (phenols, polyamines, sulfocoumarins, thioxocoumarins, catechols); CAIs occluding the entrance to the active site (coumarins and derivatives, lacosamide), as well as compounds that bind outside the active site. All these new chemotypes integrated with a general procedure for obtaining isoform-selective compounds (the tail approach) has resulted, through the guidance of rigorous X-ray crystallography experiments, in the development of highly selective CAIs for all human CA isoforms with many pharmacological applications.

Design and synthesis of benzenesulfonamide-linked imidazo[2,1-b][1,3,4]thiadiazole derivatives as carbonic anhydrase I and II inhibitors

A novel series of imidazothiadiazole-linked benzenesulfonamide derivatives (5a-t) was synthesized and subjected for screening against the four physiologically and pharmacologically relevant human carbonic anhydrase (hCA) isoforms: hCA I, II, VA, and IX. The compounds selectively inhibited hCA I and II over hCA VA and IX. Furthermore, among the two cytosolic isoforms, hCA II was more effectively inhibited as compared with hCA I. The most active compounds were 5o with K i  = 0.246 µM and 5p with K i  = 0.376 µM against hCA II, whereas compound 5f showed good inhibition against both hCA I and II with K i  = 0.493 and 0.4 µM, respectively. This class of underexplored sulfonamides may be used to design isoform-selective CA inhibitors targeting enzymes of medicinal chemistry interest.

Thiazide and other Cl-benzenesulfonamide-bearing clinical drug affinities for human carbonic anhydrases

Twelve carbonic anhydrase (CA) isoforms catalyze carbon dioxide hydration to bicarbonate and acid protons and are responsible for many biological functions in human body. Despite their vital functions, they are also responsible for, or implicated in, numerous ailments and diseases such as glaucoma, high altitude sickness, and cancer. Because CA isoforms are highly homologous, clinical drugs designed to inhibit enzymatic activity of a particular isoform, can also bind to others with similar affinity causing toxic side effects. In this study, the affinities of twelve CA isoforms have been determined for nineteen clinically used drugs used to treat hypertension related diseases, i.e. thiazides, indapamide, and metolazone. Their affinities were determined using a fluorescent thermal shift assay. Stopped flow assay and isothermal titration calorimetry were also employed on a subset of compounds and proteins to confirm inhibition of CA enzymatic activity and verify the quantitative agreement between different assays. The findings of this study showed that pharmaceuticals could bind to human CA isoforms with variable affinities and inhibit their catalytic activity, even though the drug was intended to interact with a different (non-CA) protein target. Relatively minor structural changes of the compounds may cause significant changes in affinity and selectivity for a particular CA isoform.

Structural Characterization of Carbonic Anhydrase-Arylsulfonamide Complexes Using Ultraviolet Photodissociation Mass Spectrometry

Numerous mass spectrometry-based strategies ranging from hydrogen-deuterium exchange to ion mobility to native mass spectrometry have been developed to advance biophysical and structural characterization of protein conformations and determination of protein-ligand interactions. In this study, we focus on the use of ultraviolet photodissociation (UVPD) to examine the structure of human carbonic anhydrase II (hCAII) and its interactions with arylsulfonamide inhibitors. Carbonic anhydrase, which catalyzes the conversion of carbon dioxide to bicarbonate, has been the target of countless thermodynamic and kinetic studies owing to its well-characterized active site, binding cavity, and mechanism of inhibition by hundreds of ligands. Here, we showcase the application of UVPD for evaluating structural changes of hCAII upon ligand binding on the basis of variations in fragmentation of hCAII versus hCAII-arylsulfonamide complexes, particularly focusing on the hydrophobic pocket. To extend the coverage in the midregion of the protein sequence, a supercharging agent was added to the solutions to increase the charge states of the complexes. The three arylsulfonamides examined in this study largely shift the fragmentation patterns in similar ways, despite their differences in binding affinities.

Modulating the Efficacy of Carbonic Anhydrase Inhibitors through Fluorine Substitution

The insertion of fluorine atoms and/or fluoroalkyl groups can lead to many beneficial effects in biologically active molecules, such as enhanced metabolic stability, bioavailability, lipophilicity, and membrane permeability, as well as a strengthening of protein-ligand binding interactions. However, this "magic effect" of fluorine atom(s) insertion can often be meaningless. Taking advantage of the wide range of data coming from the quest for carbonic anhydrase (CA) fluorinated inhibitors, this Minireview attempts to give "general guidelines" on how to wisely insert fluorine atom(s) within an inhibitor moiety to precisely enhance or disrupt ligand-protein interactions, depending on the target location of the fluorine substitution in the ligand. Multiple approaches such as ITC, kinetic and inhibition studies, X-ray crystallography, and NMR spectroscopy are useful in dissecting single binding contributions to the overall observed effect. The exploitation of innovative directions made in the field of protein and ligand-based fluorine NMR screening is also discussed to avoid misconduct and finely tune the exploitation of selective fluorine atom insertion in the future.

The Glitazone Class of Drugs as Carbonic Anhydrase Inhibitors-A Spin-Off Discovery from Fragment Screening

The approved drugs that target carbonic anhydrases (CA, EC 4.2.1.1), a family of zinc metalloenzymes, comprise almost exclusively of primary sulfonamides (R-SO2NH2) as the zinc binding chemotype. New clinical applications for CA inhibitors, particularly for hard-to-treat cancers, has driven a growing interest in the development of novel CA inhibitors. We recently discovered that the thiazolidinedione heterocycle, where the ring nitrogen carries no substituent, is a new zinc binding group and an alternate CA inhibitor chemotype. This heterocycle is curiously also a substructure of the glitazone class of drugs used in the treatment options for type 2 diabetes. Herein, we investigate and characterise three glitazone drugs (troglitazone 11, rosiglitazone 12 and pioglitazone 13) for binding to CA using native mass spectrometry, protein X-ray crystallography and hydrogen-deuterium exchange (HDX) mass spectrometry, followed by CA enzyme inhibition studies. The glitazone drugs all displayed appreciable binding to and inhibition of CA isozymes. Given that thiazolidinediones are not credited as a zinc binding group nor known as CA inhibitors, our findings indicate that CA may be an off-target of these compounds when used clinically. Furthermore, thiazolidinediones may represent a new opportunity for the development of novel CA inhibitors as future drugs.

Genetically Encoded Fragment-Based Discovery from Phage-Displayed Macrocyclic Libraries with Genetically Encoded Unnatural Pharmacophores

Genetically encoded macrocyclic peptide libraries with unnatural pharmacophores are valuable sources for the discovery of ligands for many targets of interest. Traditionally, generation of such libraries employs "early stage" incorporation of unnatural building blocks into the chemically or translationally produced macrocycles. Here, we describe a divergent late-stage approach to such libraries starting from readily available starting material: genetically encoded libraries of peptides. A diketone linchpin 1,5-dichloropentane-2,4-dione converts peptide libraries displayed on phage to 1,3-diketone bearing macrocyclic peptides (DKMP): shelf-stable precursors for Knorr pyrazole synthesis. Ligation of diverse hydrazine derivatives onto DKMP libraries displayed on phage that carries silent DNA-barcodes yields macrocyclic libraries in which the amino acid sequence and the pharmacophore are encoded by DNA. Selection of this library against carbonic anhydrase enriched macrocycles with benzenesulfonamide pharmacophore and nanomolar Kd. The methodology described in this manuscript can graft diverse pharmacophores into many existing genetically encoded phage libraries and significantly increase the value of such libraries in molecular discoveries.

Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors

Equilibrium binding constants (K_b) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured K_b values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured K_b (and equivalent K_d) values and the corresponding free energy change ΔG_b are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric K_b distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of K_b, K_d and ΔG_b are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of K_b, K_d, and ΔG_b values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.

Benzyl alcohol inhibits carbonic anhydrases by anchoring to the zinc coordinated water molecule

Up to date alcohols have been scarcely investigated as carbonic anhydrase (CA) inhibitors. To get more insights into the CA inhibition properties of this class of molecules, in this paper, by means of inhibition assays and X-ray crystallographic studies we report a detailed characterization of the CA inhibition properties and the binding mode to human CA II of benzyl alcohol. Results show that, although possessing a very simple scaffold, this molecule acts as a micromolar CA II inhibitor, which anchors to the enzyme active site by means of an H-bond interaction with the zinc bound solvent molecule. Taken together our results clearly indicate primary alcohols as a class of CA inhibitors that deserve to be more investigated.

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 pKa 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 pKa values of attached functional groups by additional substituents at the parent scaffold is rendered prominent.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

Genetic design of co-expressed Mesorhizobium loti carbonic anhydrase and chaperone GroELS to enhancing carbon dioxide sequestration

Mesorhizobium loti carbonic anhydrase (MlCA), an intrinsically high catalytic enzyme, has been employed for carbon dioxide capture and sequestration. However, recombinant expression of MlCA in Escherichia coli often forms inclusion bodies. Hence, protein partners such as fusion-tags and molecular chaperones are involved in regarding reduce the harshness of protein folding. TrxA-tag and GroELS have been chosen to co-express with MlCA in E. coli under an inducible T7 promoter or a constitutive J23100 promoter to compare productivity and activity. The results possessed that coupling protein partners effectively increased soluble MlCA up to 2.9-folds under T7 promoter, thus enhancing the CA activity by 120% and achieving a 5.2-folds turnover rate. Besides, it has also shifted the optimum temperature from 40 °C to 50 °C, promoted stability in the broad pH range (4.5 to 9.5) and the presence of various metal ions. Based on the in vitro assay and isothermal titration calorimetry (ITC) analysis, GroELS enhancing CA activity was due to change the intrinsic thermodynamic properties of the enzyme from endothermic to exothermic reaction (i.e., ∆H = 89.8 to -121.8 kJ/mol). Therefore, the collaboration of TrxA-MlCA with GroELS successfully augmented CO₂ biomineralization.

Simultaneous determination of thermodynamic and kinetic data by isothermal titration calorimetry

BACKGROUND

Thermodynamic and binding kinetic data increasingly support and guide the drug optimization process.

METHODS

Because ITC thermograms contain binding thermodynamic and kinetic information, an efficient protocol for the simultaneous extraction of thermodynamic and kinetic data for 1:1 protein ligand reactions from AFFINImeter kinITC in one single experiment are presented.

RESULTS

The effort to apply this protocol requires the same time as for the standard protocol but increases the precision of both thermodynamic and kinetic data.

CONCLUSIONS

The protocol enables reliable extraction of both thermodynamic and kinetic data for 1:1 protein-ligand binding reactions with improved precision compared to the 'standard protocol'.

GENERAL SIGNIFICANCE

Thermodynamic and kinetic data are recorded under exactly the same conditions in solution without any labeling or immobilization from a protein sample that is not 100% active and would otherwise render the extraction of kinetic parameters impossible.

A structure-kinetic relationship study using matched molecular pair analysis

The lifetime of a binary drug–target complex is increasingly acknowledged as an important parameter for drug efficacy and safety. With a better understanding of binding kinetics and better knowledge about kinetic parameter optimization, intentionally induced prolongation of the drug–target residence time through structural changes of the ligand could become feasible. In this study we assembled datasets from 21 publications and the K4DD (Kinetic for Drug Discovery) database to conduct large scale data analysis. This resulted in 3812 small molecules annotated to 78 different targets from five protein classes (GPCRs: 273, kinases: 3238, other enzymes: 240, HSPs: 160, ion channels: 45). Performing matched molecular pair (MMP) analysis to further investigate the structure–kinetic relationship (SKR) in this data collection allowed us to identify a fundamental contribution of a ligand's polarity to its association rate, and in selected cases, also to its dissociation rate. However, we furthermore observed that the destabilization of the transition state introduced by increased polarity is often accompanied by simultaneous destabilization of the ground state resulting in an unaffected or even worsened residence time. Supported by a set of case studies, we provide concepts on how to alter ligands in ways to trigger on-rates, off-rates, or both.

A large-scale study employing matched molecular pair (MMP) analysis to uncover the contribution of a compound's polarity to its association and dissociation rates.

Optimization of Acetazolamide-Based Scaffold as Potent Inhibitors of Vancomycin-Resistant Enterococcus

Vancomycin-resistant enterococci (VRE) are the second leading cause of hospital-acquired infections (HAIs) attributed to a drug-resistant bacterium in the United States, and resistance to the frontline treatments is well documented. To combat VRE, we have repurposed the FDA-approved carbonic anhydrase drug acetazolamide to design potent antienterococcal agents. Through structure-activity relationship optimization we have arrived at two leads possessing improved potency against clinical VRE strains from MIC = 2 μg/mL (acetazolamide) to MIC = 0.007 μg/mL (22) and 1 μg/mL (26). Physicochemical properties were modified to design leads that have either high oral bioavailability to treat systemic infections or low intestinal permeability to treat VRE infections in the gastrointestinal tract. Our data suggest the intracellular targets for the molecules are putative α-carbonic and γ-carbonic anhydrases, and homology modeling and molecular dynamics simulations were performed. Together, this study presents potential anti-VRE therapeutic options to provide alternatives for problematic VRE infections.

A Proof-of-Concept Fragment Screening of a Hit-Validated 96-Compounds Library against Human Carbonic Anhydrase II

Fragment screening is a powerful tool to identify and characterize binding pockets in proteins. We herein present the results of a proof-of-concept screening campaign of a versatile 96-entry fragment library from our laboratory against the drug target and model protein human carbonic anhydrase II. The screening revealed a novel chemotype for carbonic anhydrase inhibition, as well as less common non-covalent interaction types and unexpected covalent linkages. Lastly, different runs of the PanDDA tool reveal a practical hint for its application.

The Influence of Varying Fluorination Patterns on the Thermodynamics and Kinetics of Benzenesulfonamide Binding to Human Carbonic Anhydrase II

The fluorination of lead-like compounds is a common tool in medicinal chemistry to alter molecular properties in various ways and with different goals. We herein present a detailed study of the binding of fluorinated benzenesulfonamides to human Carbonic Anhydrase II by complementing macromolecular X-ray crystallographic observations with thermodynamic and kinetic data collected with the novel method of kinITC. Our findings comprise so far unknown alternative binding modes in the crystalline state for some of the investigated compounds as well as complex thermodynamic and kinetic structure-activity relationships. They suggest that fluorination of the benzenesulfonamide core is especially advantageous in one position with respect to the kinetic signatures of binding and that a higher degree of fluorination does not necessarily provide for a higher affinity or more favorable kinetic binding profiles. Lastly, we propose a relationship between the kinetics of binding and ligand acidity based on a small set of compounds with similar substitution patterns.

Conformational Changes in Alkyl Chains Determine the Thermodynamic and Kinetic Binding Profiles of Carbonic Anhydrase Inhibitors

Thermodynamics and kinetics of protein-ligand binding are both important aspects for the design of novel drug molecules. Presently, thermodynamic data are collected with isothermal titration calorimetry, while kinetic data are mostly derived from surface plasmon resonance. The new method of kinITC provides both thermodynamic and kinetic data from calorimetric titration measurements. The present study demonstrates the convenient collection of calorimetric data suitable for both thermodynamic and kinetic analysis for two series of congeneric ligands of human carbonic anhydrase II and correlates these findings with structural data obtained by macromolecular crystallography to shed light on the importance of shape complementarity for thermodynamics and kinetics governing a protein-ligand binding event. The study shows how minute chemical alterations change preferred ligand conformation and can be used to manipulate thermodynamic and kinetic signatures of binding. They give rise to the observation that analogous n-alkyl and n-alkyloxy derivatives of identical chain length swap their binding kinetic properties at unchanged binding affinity.

Characterizing the Interactions of Dietary Condensed Tannins with Bile Salts

This study aimed to reveal the mechanisms underlying the interaction between condensed tannins (CTs) and bile salts. The interaction mechanism was analyzed by transmission electron microscopy, exposure to various physicochemical conditions, electrophoresis, fluorescence spectroscopy, isothermal titration calorimetry, and molecular modeling. A new complex was formed from CTs and bile salts. The complex showed a negative enthalpy change and a positive entropy change, demonstrating that the main thermodynamic driving force was both entropy and enthalpy and indicating that binding occurred through hydrogen bonds and hydrophobic interactions. The analysis of the effects of CTs on the stability and digestion properties of bile salt emulsions indicated that CTs were able to inhibit lipid digestion to an extent. Our findings may provide evidence that foods rich in CTs offer health benefits by aggregating with bile salts and reducing the absorption of fat.

Application of ITC-Based Characterization of Thermodynamic and Kinetic Association of Ligands With Proteins in Drug Design

A comprehensive characterization of the thermodynamic and kinetic profiling of ligands binding to a given target protein is crucial for the hit selection as well as the hit-to-lead-to-drug evolution. Isothermal titration calorimetry (ITC), widely known as an invaluable tool to measure the thermodynamic data, has recently found its way to determine the binding kinetics too. The extensive application of ITC in measurement of both thermodynamic and kinetic data manifests unique roles of ITC in drug discovery and development. This mini-review concentrates on elaborating how to gain the thermodynamic and kinetic data using ITC, highlighting the importance of these data in lead discovery and optimization, and intends to provide an overview of the technical and conceptual advances that offer unprecedented access to protein–ligand recognition by ITC measurement.

KinITC-One Method Supports both Thermodynamic and Kinetic SARs as Exemplified on FimH Antagonists

Affinity data, such as dissociation constants (K_D ) or inhibitory concentrations (IC₅₀ ), are widely used in drug discovery. However, these parameters describe an equilibrium state, which is often not established in vivo due to pharmacokinetic effects and they are therefore not necessarily sufficient for evaluating drug efficacy. More accurate indicators for pharmacological activity are the kinetics of binding processes, as they shed light on the rate of formation of protein-ligand complexes and their half-life. Nonetheless, although highly desirable for medicinal chemistry programs, studies on structure-kinetic relationships (SKR) are still rare. With the recently introduced analytical tool kinITC this situation may change, since not only thermodynamic but also kinetic information of the binding process can be deduced from isothermal titration calorimetry (ITC) experiments. Using kinITC, ITC data of 29 mannosides binding to the bacterial adhesin FimH were re-analyzed to make their binding kinetics accessible. To validate these kinetic data, surface plasmon resonance (SPR) experiments were conducted. The kinetic analysis by kinITC revealed that the nanomolar affinities of the FimH antagonists arise from both (i) an optimized interaction between protein and ligand in the bound state (reduced off-rate constant k_off ) and (ii) a stabilization of the transition state or a destabilization of the unbound state (increased on-rate constant k_on ). Based on congeneric ligand modifications and structural input from co-crystal structures, a strong relationship between the formed hydrogen-bond network and k_off could be concluded, whereas electrostatic interactions and conformational restrictions upon binding were found to have mainly an impact on k_on .

Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery

Simple comparative correlation analyses and quantitative structure-kinetics relationship (QSKR) models highlight the interplay of kinetic rates and binding affinity as an essential feature in drug design and discovery. The choice of the molecular series, and their structural variations, used in QSKR modeling is fundamental to understanding the mechanistic implications of ligand and/or drug-target binding and/or unbinding processes. Here, we discuss the implications of linear correlations between kinetic rates and binding affinity constants and the relevance of the computational approaches to QSKR modeling.

An overview of recent molecular dynamics applications as medicinal chemistry tools for the undruggable site challenge

Molecular dynamics (MD) has become increasingly popular due to the development of hardware and software solutions and the improvement in algorithms, which allowed researchers to scale up calculations in order to speed them up. MD simulations are usually used to address protein folding issues or protein-ligand complex stability through energy profile analysis over time. In recent years, the development of new tools able to deeply explore a potential energy surface (PES) has allowed researchers to focus on the dynamic nature of the binding recognition process and binding-induced protein conformational changes. Moreover, modern approaches have been demonstrated to be effective and reliable in calculating some kinetic and thermodynamic parameters behind the host-guest recognition process. Starting from all of these considerations, several efforts have been made in order to integrate MD within the virtual screening process in drug discovery. Knowledge retrieved from MD can, in fact, be exploited as a starting point to build pharmacophores or docking constraints in the early stage of the screening campaign as well as to define key features, in order to unravel hidden binding modes and help the optimisation of the molecular structure of a lead compound. Based on these outcomes, researchers are nowadays using MD as an invaluable tool to discover and target previously considered undruggable binding sites, including protein-protein interactions and allosteric sites on a protein surface. As a matter of fact, the use of MD has been recognised as vital to the discovery of selective protein-protein interaction modulators. The use of a dynamic overview on how the host-guest recognition occurs and of the relative conformational modifications induced allows researchers to optimise small molecules and small peptides capable of tightly interacting within the cleft between two proteins. In this review, we aim to present the most recent applications of MD as an integrated tool to be used in the rational design of small molecules or small peptides able to modulate undruggable targets, such as allosteric sites and protein-protein interactions.

Novel fluorinated carbonic anhydrase IX inhibitors reduce hypoxia-induced acidification and clonogenic survival of cancer cells

Human carbonic anhydrase (CA) IX has emerged as a promising anticancer target and a diagnostic biomarker for solid hypoxic tumors. Novel fluorinated CA IX inhibitors exhibited up to 50 pM affinity towards the recombinant human CA IX, selectivity over other CAs, and direct binding to Zn(II) in the active site of CA IX inducing novel conformational changes as determined by X-ray crystallography. Mass spectrometric gas-analysis confirmed the CA IX-based mechanism of the inhibitors in a CRISPR/Cas9-mediated CA IX knockout in HeLa cells. Hypoxia-induced extracellular acidification was significantly reduced in HeLa, H460, MDA-MB-231, and A549 cells exposed to the compounds, with the IC₅₀ values up to 1.29 nM. A decreased clonogenic survival was observed when hypoxic H460 3D spheroids were incubated with our lead compound. These novel compounds are therefore promising agents for CA IX-specific therapy.

Dynamics and Interactions of a 29 kDa Human Enzyme Studied by Solid-State NMR

Solid-state NMR has been employed for characterization of a broad range of biomacromolecules and supramolecular assemblies. However, because of limitations in sensitivity and resolution, the size of the individual monomeric units has rarely exceeded 15 kDa. As such, enzymes, which are often more complex and comprise long peptide chains, have not been easily accessible, even though manifold desirable information could potentially be provided by solid-state NMR studies. Here, we demonstrate that more than 1200 backbone and side-chain chemical shifts can be reliably assessed from minimal sample quantities for a 29 kDa human enzyme of the carbonic anhydrase family, giving access to its backbone dynamics and intermolecular interactions with a small-molecule inhibitor. The possibility of comprehensive assessment of enzymes in this molecular-weight regime without molecular-tumbling-derived limitations enables the study of residue-specific properties important for their mode of action as well as for pharmacological interference in this and many other enzymes.

"To Be or Not to Be" Protonated: Atomic Details of Human Carbonic Anhydrase-Clinical Drug Complexes by Neutron Crystallography and Simulation

Human carbonic anhydrases (hCAs) play various roles in cells, and have been drug targets for decades. Sequence similarities of hCA isoforms necessitate designing specific inhibitors, which requires detailed structural information for hCA-inhibitor complexes. We present room temperature neutron structures of hCA II in complex with three clinical drugs that provide in-depth analysis of drug binding, including protonation states of the inhibitors, hydration water structure, and direct visualization of hydrogen-bonding networks in the enzyme's active site. All sulfonamide inhibitors studied bind to the Zn metal center in the deprotonated, anionic, form. Other chemical groups of the drugs can remain neutral or be protonated when bound to hCA II. MD simulations have shown that flexible functional groups of the inhibitors may alter their conformations at room temperature and occupy different sub-sites. This study offers insights into the design of specific drugs to target cancer-related hCA isoform IX.

"Seriously Sweet": Acesulfame K Exhibits Selective Inhibition Using Alternative Binding Modes in Carbonic Anhydrase Isoforms

Human carbonic anhydrase IX (CA IX) is upregulated in neoplastic tissues; as such, it is studied as a drug target for anticancer chemotherapy. Inhibition of CA IX has been shown to be therapeutically favorable in terms of reducing tumor growth. Previously, saccharin, a commonly used artificial sweetener, has been observed to selectively inhibit CA IX over other CA isoforms. In this study, X-ray crystallography showed acesulfame potassium (Ace K) binding directly to the catalytic zinc in CA IX (mimic) and through a bridging water in CA II. This modulation in binding is reflected in the binding constants, with Ace K inhibiting CA IX but not other CA isoforms. Hence, this study establishes the potential of Ace K (an FDA approved food additive) as a lead compound in the design and development of CA IX specific inhibitors.

Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme

Human carbonic anhydrase II (HCA II) is an enzyme that catalyzes the reversible hydration of CO2 into bicarbonate (HCO3-) and a proton (H+) as well as other reactions at an extremely high rate. This enzyme plays fundamental roles in human physiology/pathology, such as controlling the pH level in cells and so on. However, the binding mechanism between apo-HCA II and CO2 or other ligands as well as related conformational changes remains poorly understood, and atomic investigation into it could promote our understanding of related internal physiological/pathological mechanisms. In this study, long-time atomic molecular dynamics simulations as well as the clustering and free-energy analysis were performed to reveal the dynamics of apo-HCA II as well as the mechanism upon ligand binding. Our simulations indicate that the crystallographic B-factors considerably underestimate the loop dynamics: multiple conformations can be adopted by loops 1 and 2, especially for loop 1 because loop 1 is one side of the binding pocket, and its left-to-right movement can compress or extend the binding pocket, leading to one inactive (closed) state, three intermediate (semiopen) states, and one active (open) state; CO2 cannot get into the binding pocket of the inactive state but can get into those of intermediate and active states. The coexistence of multiple conformational states proposes a possible conformational selection model for the binding mechanism between apo-HCA II and CO2 or other ligands, revising our previous view of its functional mechanism of conformational change upon ligand binding and offering valuable structural insights into the workings of HCA II.

Decoupled Associative and Dissociative Processes in Strong yet Highly Dynamic Host-Guest Complexes

Kinetics and thermodynamics in supramolecular systems are intimately linked, yet both are independently important for application in sensing assays and stimuli-responsive switching/self-healing of materials. Host-guest interactions are of particular interest in many water-based materials, sensing, and drug delivery applications. Herein we investigate the binding dynamics of a variety of electron-rich aromatic moieties forming hetero-ternary complexes with the macrocycle cucurbit[8]uril (CB[8]) and an auxiliary guest, dimethyl viologen, with high selectivity and equilibrium binding constants (K_eq up to 10¹⁴ M^-2). Using stopped-flow spectrofluorimetry, association rate constants were observed to approach the diffusion limit and were found to be insensitive to the structure of the guest. Conversely, the dissociation rate constants of the ternary complexes varied dramatically with the guest structure and were correlated with the thermodynamic binding selectivity. Hence differing molecular features were found to contribute to the associative and dissociative processes, mimicking naturally occurring reactions and giving rise to a decoupling of these kinetic parameters. Moreover, we demonstrate the ability to exploit these phenomena and selectively perturb the associative process with external stimuli (e.g., viscosity and pressure). Significantly, these complexes exhibit increased binding equilibria with increasing pressure, with important implications for the application of the CB[8] ternary complex for the formation of hydrogels, as these gels exhibit unprecedented pressure-insensitive rheological properties. A high degree of flexibility therefore exists in the design of host-guest systems with tunable kinetic and thermodynamic parameters for tailor-made applications across a broad range of fields.

Characterizing Drug-Target Residence Time with Metadynamics: How To Achieve Dissociation Rate Efficiently without Losing Accuracy against Time-Consuming Approaches

Drug-target residence time plays a vital role in drug efficacy. However, there is still no effective strategy to predict drug residence time. Here, we propose to use the optimized (or minimized) structures derived from holo-state proteins to calculate drug residence time, which could give a comparable or even better prediction accuracy compared with those calculated utilizing a large number of molecular dynamics (MD) structures based on the Poisson process. Besides, in addition to the Poisson process, one may use fewer samples for predicting residence time due to the reason that, in a large extent, the calculated drug residence time is stable and independent of the number of samples used for the prediction. With remarkably reduced computational load, the proposed strategy may be promising for large-scale drug residence time prediction, such as post-processing in virtual screening (VS) and lead compound optimization.

Mechanistic Explanation of the Weak Carbonic Anhydrase's Esterase Activity

In order to elucidate the elementary mechanism of the promiscuous esterase activity of human carbonic anhydrase (h-CA), we present an accurate theoretical investigation on the hydrolysis of fully-acetylated d-glucose functionalized as sulfamate. This h-CA's inhibitor is of potential relevance in cancer therapy. The study has been performed within the framework of three-layer ONIOM (QM-high:QM'-medium:MM-low) hybrid approach. The computations revealed that the hydrolysis process is not energetically favored, in agreement with the observed weak carbonic anhydrase's esterase activity.