Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure–thermodynamics and structure–kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure–energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure–thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design

Switching the Inhibitor‐Enzyme Recognition Profile via Chimeric Carbonic Anhydrase XII

A key part of the optimization of small molecules in pharmaceutical inhibitor development is to vary the molecular design to enhance complementarity of chemical features of the compound with the positioning of amino acids in the active site of a target enzyme. Typically this involves iterations of synthesis, to modify the compound, and biophysical assay, to assess the outcomes. Selective targeting of the anti-cancer carbonic anhydrase isoform XII (CA XII), this process is challenging because the overall fold is very similar across the twelve CA isoforms. To enhance drug development for CA XII we used a reverse engineering approach where mutation of the key six amino acids in the active site of human CA XII into the CA II isoform was performed to provide a protein chimera (chCA XII) which is amenable to structure-based compound optimization. Through determination of structural detail and affinity measurement of the interaction with over 60 compounds we observed that the compounds that bound CA XII more strongly than CA II, switched their preference and bound more strongly to the engineered chimera, chCA XII, based on CA II, but containing the 6 key amino acids from CA XII, behaved as CA XII in its compound recognition profile. The structures of the compounds in the chimeric active site also resembled those determined for complexes with CA XII, hence validating this protein engineering approach in the development of new inhibitors

PLBD: protein–ligand binding database of thermodynamic and kinetic intrinsic parameters

We introduce a protein–ligand binding database (PLBD) that presents thermodynamic and kinetic data of reversible protein interactions with small molecule compounds. The manually curated binding data are linked to protein–ligand crystal structures, enabling structure–thermodynamics correlations to be determined. The database contains over 5500 binding datasets of 556 sulfonamide compound interactions with the 12 catalytically active human carbonic anhydrase isozymes defined by fluorescent thermal shift assay, isothermal titration calorimetry, inhibition of enzymatic activity and surface plasmon resonance. In the PLBD, the intrinsic thermodynamic parameters of interactions are provided, which account for the binding-linked protonation reactions. In addition to the protein–ligand binding affinities, the database provides calorimetrically measured binding enthalpies, providing additional mechanistic understanding. The PLBD can be applied to investigations of protein–ligand recognition and could be integrated into small molecule drug design

Predation on an Upper Trophic Marine Predator, the Steller Sea Lion: Evaluating High Juvenile Mortality in a Density Dependent Conceptual Framework

The endangered western stock of the Steller sea lion (Eumetopias jubatus) – the largest of the eared seals – has declined by 80% from population levels encountered four decades ago. Current overall trends from the Gulf of Alaska to the Aleutian Islands appear neutral with strong regional heterogeneities. A published inferential model has been used to hypothesize a continuous decline in natality and depressed juvenile survival during the height of the decline in the mid-late 1980's, followed by the recent recovery of juvenile survival to pre-decline rates. However, these hypotheses have not been tested by direct means, and causes underlying past and present population trajectories remain unresolved and controversial. We determined post-weaning juvenile survival and causes of mortality using data received post-mortem via satellite from telemetry transmitters implanted into 36 juvenile Steller sea lions from 2005 through 2011. Data show high post-weaning mortality by predation in the eastern Gulf of Alaska region. To evaluate the impact of such high levels of predation, we developed a conceptual framework to integrate density dependent with density independent effects on vital rates and population trajectories. Our data and model do not support the hypothesized recent recovery of juvenile survival rates and reduced natality. Instead, our data demonstrate continued low juvenile survival in the Prince William Sound and Kenai Fjords region of the Gulf of Alaska. Our results on contemporary predation rates combined with the density dependent conceptual framework suggest predation on juvenile sea lions as the largest impediment to recovery of the species in the eastern Gulf of Alaska region. The framework also highlights the necessity for demographic models based on age-structured census data to incorporate the differential impact of predation on multiple vital rates