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

Microsecond Timescale Conformational Dynamics of a Small-Molecule Ligand within the Active Site of a Protein

The possible internal dynamics of non-isotope-labeled small-molecule ligands inside a target protein is inherently difficult to capture. Whereas high crystallographic temperature factors can denote either static disorder or motion, even moieties with very low B-factors can be subject to vivid motion between symmetry-related sites. Here we report the experimental identification of internal μs timescale dynamics of a high-affinity, natural-abundance ligand tightly bound to the enzyme human carbonic anhydrase II (hCAII) even within a crystalline lattice. The rotamer jumps of the ligand's benzene group manifest themselves both, in solution and fast magic-angle spinning solid-state NMR 1 H R1ρ relaxation dispersion, for which we obtain further mechanistic insights from molecular-dynamics (MD) simulations. The experimental confirmation of rotameric jumps in bound ligands within proteins in solution or the crystalline state may improve understanding of host-guest interactions in biology and supra-molecular chemistry and may facilitate medicinal chemistry for future drug campaigns.

Do Go Chasing Waterfalls: Enoyl Reductase (FabI) in Complex with Inhibitors Stabilizes the Tetrameric Structure and Opens Water Channels

The enzyme enoyl-ACP reductase (FabI) is the limiting step of the membrane's fatty acid biosynthesis in bacteria and a druggable target for novel antibacterial agents. The FabI active form is a homotetramer, which displays the highest affinity to inhibitors. Herein, molecular dynamics studies were carried out using the structure of FabI in complex with known inhibitors to investigate their effects on tetramerization. Our results suggest that multimerization is essential for the integrity of the catalytic site and that inhibitor binding enables the multimerization by stabilizing the substrate binding loop (SBL, L:195-200) coupled with changes in the H4/5 (QR interface). We also observed that AFN-1252 (naphtpyridinone derivative) promotes unique conformational changes affecting monomer-monomer interfaces. These changes are induced by AFN-1252 interaction with key residues in the binding sites (Ala95, Tyr146, and Tyr156). In addition, the analysis of water trajectories indicated that AFN-1252 complexes allow more water molecules to enter the binding site than triclosan and MUT056399 complexes. FabI-AFN-1252 simulations show accumulation of water molecules near the Tyr146/147 pocket, which can become a hotspot to the design of novel FabI inhibitors.

The feasibility of determining kinetic constants from isothermal titration calorimetry data

Isothermal titration calorimetry (ITC) has long been established as an excellent means to determine the thermodynamic parameters of biomolecular interactions. More recently, efforts have focused on exploiting the power/time trace (the "thermogram") resulting from ITC experiments to glean kinetic association and dissociation rates for these interactions. The ability to do so rests on the ability of algorithms to simulate with high accuracy the output of the calorimeter. Thus, several critical factors must be taken into account: the injection protocol, the kinetics of the interaction, accurate discovery of the instrumental response to heat signals, and the addition of unrelated signals. All of these aspects of extracting kinetic constants from thermograms have been considered and addressed in the current work. To validate the resultant methods, we performed several ITC experiments, titrating small-molecule inhibitors into solutions of bovine carbonic anhydrase II or titrating lysozyme into solutions of anti-lysozyme nanobodies. We found that our methods could arrive at kinetic constants that were close to the known values for these interactions taken from other methods. Finally, the effort to improve ITC kinetic characterizations uncovered a set of best practices for both the calorimetric experiment and the subsequent analyses (termed "kinetically optimized ITC" or "KO-ITC") that is detailed in this work.

Computational Prediction of the Binding Pose of Metal-Binding Pharmacophores

Computational modeling of inhibitors for metalloenzymes in virtual drug development campaigns has proven challenging. To overcome this limitation, a technique for predicting the binding pose of metal-binding pharmacophores (MBPs) is presented. Using a combination of density functional theory (DFT) calculations and docking using a genetic algorithm, inhibitor binding was evaluated in silico and compared with inhibitor-enzyme cocrystal structures. The predicted binding poses were found to be consistent with the cocrystal structures. The computational strategy presented represents a useful tool for predicting metalloenzyme-MBP interactions.

Decisive role of water and protein dynamics in residence time of p38α MAP kinase inhibitors

Target residence time plays a crucial role in the pharmacological activity of small molecule inhibitors. Little is known, however, about the underlying causes of inhibitor residence time at the molecular level, which complicates drug optimization processes. Here, we employ all-atom molecular dynamics simulations (~400 μs in total) to gain insight into the binding modes of two structurally similar p38α MAPK inhibitors (type I and type I½) with short and long residence times that otherwise show nearly identical inhibitory activities in the low nanomolar IC₅₀ range. Our results highlight the importance of protein conformational stability and solvent exposure, buried surface area of the ligand and binding site resolvation energy for residence time. These findings are further confirmed by simulations with a structurally diverse short residence time inhibitor SB203580. In summary, our data provide guidance in compound design when aiming for inhibitors with improved target residence time.

The molecular determinants of the residence time of a small molecule inhibitor at its target protein are not well understood. Here, Pantsar et al. show that the target protein’s conformational stability and solvent exposure are key factors governing the target residence time of kinase inhibitors.

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.

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.

Applications of isothermal titration calorimetry in pure and applied research from 2016 to 2020

The last 5 years have seen a series of advances in the application of isothermal titration microcalorimetry (ITC) and interpretation of ITC data. ITC has played an invaluable role in understanding multiprotein complex formation including proteolysis-targeting chimeras (PROTACS), and mitochondrial autophagy receptor Nix interaction with LC3 and GABARAP. It has also helped elucidate complex allosteric communication in protein complexes like trp RNA-binding attenuation protein (TRAP) complex. Advances in kinetics analysis have enabled the calculation of kinetic rate constants from pre-existing ITC data sets. Diverse strategies have also been developed to study enzyme kinetics and enzyme-inhibitor interactions. ITC has also been applied to study small molecule solvent and solute interactions involved in extraction, separation, and purification applications including liquid-liquid separation and extractive distillation. Diverse applications of ITC have been developed from the analysis of protein instability at different temperatures, determination of enzyme kinetics in suspensions of living cells to the adsorption of uremic toxins from aqueous streams.

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.

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.