Intrinsic Thermodynamics and Structure Correlation of Benzenesulfonamides with a Pyrimidine Moiety Binding to Carbonic Anhydrases I,II, VII,XII, and XIII

The early stage of drug discovery is often based on selecting the highest affinity lead compound. To this end the structural and energetic characterization of the binding reaction is important. The binding energetics can be resolved into enthalpic and entropic contributions to the binding Gibbs free energy. Most compound binding reactions are coupled to the absorption or release of protons by the protein or the compound. A distinction between the observed and intrinsic parameters of the binding energetics requires the dissection of the protonation/deprotonation processes. Since only the intrinsic parameters can be correlated with molecular structural perturbations associated with complex formation, it is these parameters that are required for rational drug design. Carbonic anhydrase (CA) isoforms are important therapeutic targets to treat a range of disorders including glaucoma, obesity, epilepsy, and cancer. For effective treatment isoform-specific inhibitors are needed. In this work we investigated the binding and protonation energetics of sixteen [(2-pyrimidinylthio)acetyl]benzenesulfonamide CA inhibitors using isothermal titration calorimetry and fluorescent thermal shift assay. The compounds were built by combining four sulfonamide headgroups with four tailgroups yielding 16 compounds. Their intrinsic binding thermodynamics showed the limitations of the functional group energetic additivity approach used in fragment-based drug design, especially at the level of enthalpies and entropies of binding. Combined with high resolution crystal structural data correlations were drawn between the chemical functional groups on selected inhibitors and intrinsic thermodynamic parameters of CA-inhibitor complex formation.     

Group additivity schemes for the calculation of the partial molar heat capacities and volumes of unfolded proteins in aqueous solution

A critical review is given of the present state of group additivity schemes for the calculation of partial molar volumes and heat capacities of unfolded proteins. The comparison between the experimental values and the predictions based on the different models shows clearly that only the peptide-based additivity scheme represents properly both the absolute values and the temperature dependence of these thermodynamic quantities.     

A new strategy for improved secondary screening and lead optimization using high-resolution SPR characterization of compound-target interactions

Biophysical label-free assays such as those based on SPR are essential tools in generating high-quality data on affinity, kinetic, mechanistic and thermodynamic aspects of interactions between target proteins and potential drug candidates. Here we show examples of the integration of SPR with bioinformatic approaches and mutation studies in the early drug discovery process. We call this combination 'structure-based biophysical analysis'. Binding sites are identified on target proteins using information that is either extracted from three-dimensional structural analysis (X-ray crystallography or NMR), or derived from a pharmacore model based on known binders. The binding site information is used for in silico screening of a large substance library (e.g. available chemical directory), providing virtual hits. The three-dimensional structure is also used for the design of mutants where the binding site has been impaired. The wild-type target and the impaired mutant are then immobilized on different spots of the sensor chip and the interactions of compounds with the wild-type and mutant are compared in order to identify selective binders for the binding site of the target protein. This method can be used as a cost-effective alternative to high-throughput screening methods in cases when detailed binding site information is available. Here, we present three examples of how this technique can be applied to provide invaluable data during different phases of the drug discovery process.     

Additivity in protein-DNA interactions: how good an approximation is it?

Man and Stormo and Bulyk et al. recently presented their results on the study of the DNA binding affinity of proteins. In both of these studies the main conclusion is that the additivity assumption, usually applied in methods to search for binding sites, is not true. In the first study, the analysis of binding affinity data from the Mnt repressor protein bound to all possible DNA (sub)targets at positions 16 and 17 of the binding site, showed that those positions are not independent. In the second study, the authors analysed DNA binding affinity data of the wild-type mouse EGR1 protein and four variants differing on the middle finger. The binding affinity of these proteins was measured to all 64 possible trinucleotide (sub)targets of the middle finger using microarray technology. The analysis of the measurements also showed interdependence among the positions in the DNA target. In the present report, we review the data of both studies and we re- analyse them using various statistical methods, including a comparison with a multiple regression approach. We conclude that despite the fact that the additivity assumption does not fit the data perfectly, in most cases it provides a very good approximation of the true nature of the specific protein–DNA interactions. Therefore, additive models can be very useful for the discovery and prediction of binding sites in genomic DNA.     

NMR and X-ray analysis of structural additivity in metal binding site-swapped hybrids of rubredoxin

Background  

Chimeric hybrids derived from the rubredoxins of Pyrococcus furiosus (Pf) and Clostridium pasteurianum (Cp) provide a robust system for the characterization of protein conformational stability and dynamics in a differential mode. Interchange of the seven nonconserved residues of the metal binding site between the Pf and Cp rubredoxins yields a complementary pair of hybrids, for which the sum of the thermodynamic stabilities is equal to the sum for the parental proteins. Furthermore, the increase in amide hydrogen exchange rates for the hyperthermophile-derived metal binding site hybrid is faithfully mirrored by a corresponding decrease for the complementary hybrid that is derived from the less thermostable rubredoxin, indicating a degree of additivity in the conformational fluctuations that underlie these exchange reactions.     

Enthalpy screen of drug candidates

The enthalpic and entropic contributions to the binding affinity of drug candidates have been acknowledged to be important determinants of the quality of a drug molecule. These quantities, usually summarized in the thermodynamic signature, provide a rapid assessment of the forces that drive the binding of a ligand. Having access to the thermodynamic signature in the early stages of the drug discovery process will provide critical information towards the selection of the best drug candidates for development. In this paper, the Enthalpy Screen technique is presented. The enthalpy screen allows fast and accurate determination of the binding enthalpy for hundreds of ligands. As such, it appears to be ideally suited to aid in the ranking of the hundreds of hits that are usually identified after standard high throughput screening.     

Stability of macromolecular complexes

Macromolecular interactions are crucial in numerous biologic processes, yet few general principles are available that establish firm expectations for the strength of these interactions or the expected contribution of specific forces. The simplest principle would be a monotonic increase in interactions as the size of the interface grows. The exact relationship might be linear or nonlinear depending on the nature of the forces involved. Simple "linear-free energy" relationships based on atomic properties have been well documented, for example, additivity for the interaction of small molecules with solvent, and, recently, have been explored for ligand-receptor interactions. Horton and Lewis propose such additivity based on buried surface area for protein-protein complexes. We investigated macromolecular interactions and found that the highest-affinity complexes do not fulfill this simple expectation. Instead, binding free energies of the tightest macromolecular complexes are roughly constant, independent of interface size, with the notable exception of DNA duplexes. By comparing these results to an earlier study of protein-ligand interactions we find that: (1) The maximum affinity is approximately 1.5 kcal/mol per nonhydrogen atom or 120 cal/mol A(2) of buried surface area, comparable to results of our earlier work; (2) the lack of an increase in affinity with interface size is likely due to nonthermodynamic factors, such as functional and evolutionary constraints rather than some fundamental physical limitation. The implication of these results have some importance for molecular design because they suggest that: (1) The stability of any given complex can be increased significantly if desired; (2) small molecule inhibitors of macromolecular interactions are feasible; and (3) different functional classes of protein-protein complexes exhibit differences in maximal stability, perhaps in response to differing evolutionary pressures. These results are consistent with the widespread observation that proteins have not evolved to maximize thermodynamic stability, but are only marginally stable.     

Higher throughput calorimetry: opportunities, approaches and challenges

Higher throughput thermodynamic measurements can provide value in structure-based drug discovery during fragment screening, hit validation, and lead optimization. Enthalpy can be used to detect and characterize ligand binding, and changes that affect the interaction of protein and ligand can sometimes be detected more readily from changes in the enthalpy of binding than from the corresponding free-energy changes or from protein-ligand structures. Newer, higher throughput calorimeters are being incorporated into the drug discovery process. Improvements in titration calorimeters come from extensions of a mature technology and face limitations in scaling. Conversely, array calorimetry, an emerging technology, shows promise for substantial improvements in throughput and material utilization, but improved sensitivity is needed.     

Protein-drug interactions: characterization of inhibitor binding in complexes of DHFR with trimethoprim and related derivatives

Structural and thermodynamic interactions for the binding of trimethoprim and related congeners to the binary complex of dihydrofolate reductase (from chicken) and NADPH are explored using free energy simulation methods. Good agreement between structures from experimental X-ray refinement and molecular dynamics simulations is found for the complexes. Agreement with thermodynamic measurements is found as well. Our thermodynamic calculations suggest that entropic contributions and desolvation thermodynamics can play a crucial role in overall binding, and that extreme care must be taken in the use of simple model building to rationalize or predict protein-drug binding.     

Exploring additivity effects of double mutations on the binding affinity of protein‐protein complexes

Additivity in binding affinity of protein‐protein complexes refers to the change in free energy of binding (ΔΔG_bind) for double (or multiple) mutations which is approximately equal to the sum of their corresponding single mutation ΔΔG_bind values. In this study, we have explored the additivity effect of double mutants, which shows a linear relationship between the binding affinity of double and sum of single mutants with a correlation of 0.90. However, the comparison of ΔΔG_bind values showed a mean absolute deviation of 0.86 kcal/mol, and 25.6% of the double mutants show a deviation of more than 1 kcal/mol, which are identified as non‐additive. The additivity effects have been analyzed based on the influence of structural features such as accessible surface area, long range order, binding propensity change, surrounding hydrophobicity, flexibility, atomic contacts between the mutations and distance between the 2 mutations. We found that non‐additive mutations tend to be closer to each other and have more contacts. We have also used machine learning methods to discriminate additive and non‐additive mutations using structure‐based features, which showed the accuracies in the range of 0.77–0.92 for protein‐protein complexes belonging to different functions. Further, we have compared the additivity effects of protein stability along with binding affinity and explored the similarities and differences between them. The results obtained in this study provide insights into the effects of various structural features on binding affinity of double mutants, and will aid the development of accurate methods to predict the binding affinity of double mutants.     

Ligand binding affinity determined by temperature-dependent circular dichroism: cyclin-dependent kinase 2 inhibitors

To support drug discovery efforts for cyclin-dependent kinase 2 (CDK2), a moderate-throughput binding assay that can rank order or estimate the affinity of lead inhibitors has been developed. The method referred to as temperature-dependent circular dichroism (TdCD) uses the classical temperature-dependent unfolding of proteins by circular dichroism (CD) to measure the degree of protein unfolding in the absence and presence of potential inhibitors. The midpoint of unfolding is the Tm value. Rank ordering the affinity and predictions of the dissociation constant of compounds is obtained by measuring the increase in Tm for different protein-inhibitor complexes. This is the first time an extensive characterization of the TdCD method has been described for characterizing lead inhibitors in a drug discovery mode. The method has several favorable properties. Using the new six-cell Peltier temperature controller for the Jasco 810 spectropolarimeter, one can determine the affinity of 12-18 compounds per day. The method also requires only 20-40 microg protein per sample and can be used to estimate the affinity of compounds with dissociation constants of picomolar to micromolar. An important property of the method for lead discovery is that dissociation constants of approximately 5 microM can be estimated from a single experiment using a low concentration of compound such as 20 microM, which is generally low enough for most small molecules to be soluble for testing. In addition, the method does not require labeling the compound or protein. Although other methods such as isothermal titration calorimetry (ITC) can provide a full thermodynamic characterization of binding, ITC requires 1-2 mg protein per sample, cannot readily determine binding constants below nanomolar values, is most versatile with soluble compounds, and has a throughput of two to three experiments per day. The ITC method is not usually used in a high-throughput drug discovery mode; however, using the thermodynamic information from several ITC experiments can make the TdCD method very robust in determining reliable binding constants. Using the kinase inhibitors BMS-250595, purvalanol B, AG-12275, flavopiridol, and several other compounds, it is demonstrated that one can obtain excellent comparisons between the Kd values of binding to CDK2 obtained by TdCD and ITC.     

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.     

Group additivity calculations of the thermodynamic properties of unfolded proteins in aqueous solution: a critical comparison of peptide-based and HKF models

A recent paper in this journal [Amend and Helgeson, Biophys. Chem. 84 (2000) 105] presented a new group additivity model to calculate various thermodynamic properties of unfolded proteins in aqueous solution. The parameters given for the revised Helgeson-Kirkham-Flowers (HKF) equations of state for all the constituent groups of unfolded proteins can be used, in principle, to calculate the partial molar heat capacity, C(o)p.2, and volume, V2(0), at infinite dilution of any polypeptide. Calculations of the values of C(o)p.2 and V2(0) for several polypeptides have been carried out to test the predictive utility of the HKF group additivity model. The results obtained are in very poor agreement with experimental data, and also with results calculated using a peptide-based group additivity model. A critical assessment of these two additivity models is presented.     

Calculation of the standard molal thermodynamic properties of aqueous biomolecules at elevated temperatures and pressures II. Unfolded proteins

Equations of state for completely unfolded proteins have been generated from group additivity algorithms and the revised Helgeson-Kirkham-Flowers (HKF) equations of state to compute the standard molal thermodynamic properties of these molecules at elevated temperatures and pressures. The requisite equations of state parameters were computed from those of groups retrieved by regression of experimental calorimetric and densimetric data reported in the literature. This approach permits calculation of the standard molal thermodynamic properties as a function of temperature and pressure for any completely unfolded protein for which the amino acid sequence is known. Calculations of this kind have been carried out for 11 thermophilic proteins. The thermodynamic properties reported below can be combined with those for protein unfolding to compute the corresponding properties of completely folded (i.e. native) proteins.     

Determinants of Affinity and Proteolytic Stability in Interactions of Kunitz Family Protease Inhibitors with Mesotrypsin

An important functional property of protein protease inhibitors is their stability to proteolysis. Mesotrypsin is a human trypsin that has been implicated in the proteolytic inactivation of several protein protease inhibitors. We have found that bovine pancreatic trypsin inhibitor (BPTI), a Kunitz protease inhibitor, inhibits mesotrypsin very weakly and is slowly proteolyzed, whereas, despite close sequence and structural homology, the Kunitz protease inhibitor domain of the amyloid precursor protein (APPI) binds to mesotrypsin 100 times more tightly and is cleaved 300 times more rapidly. To define features responsible for these differences, we have assessed the binding and cleavage by mesotrypsin of APPI and BPTI reciprocally mutated at two nonidentical residues that make direct contact with the enzyme. We find that Arg at P₁ (versus Lys) favors both tighter binding and more rapid cleavage, whereas Met (versus Arg) at P′₂ favors tighter binding but has minimal effect on cleavage. Surprisingly, we find that the APPI scaffold greatly enhances proteolytic cleavage rates, independently of the binding loop. We draw thermodynamic additivity cycles analyzing the interdependence of P₁ and P′₂ substitutions and scaffold differences, finding multiple instances in which the contributions of these features are nonadditive. We also report the crystal structure of the mesotrypsin·APPI complex, in which we find that the binding loop of APPI displays evidence of increased mobility compared with BPTI. Our data suggest that the enhanced vulnerability of APPI to mesotrypsin cleavage may derive from sequence differences in the scaffold that propagate increased flexibility and mobility to the binding loop.     

Collagen model peptides: Sequence dependence of triple-helix stability

The triple helix is a specialized protein motif, found in all collagens as well as in noncollagenous proteins involved in host defense. Peptides will adopt a triple-helical conformation if the sequence contains its characteristic features of Gly as every third residue and a high content of Pro and Hyp residues. Such model peptides have proved amenable to structural studies by x-ray crystallography and NMR spectroscopy, suitable for thermodynamic and kinetic analysis, and a valuable tool in characterizing the binding activities of the collagen triple helix. A systematic approach to understanding the amino acid sequence dependence of the collagen triple helix has been initiated, based on a set of host-guest peptides of the form, (Gly-Pro-Hyp)(3)-Gly-X-Y-(Gly-Pro-Hyp)(4). Comparison of their thermal stabilities has led to a propensity scale for the X and Y positions, and the additivity of contributions of individual residues is now under investigation. The local and global stability of the collagen triple helix is normally modulated by the residues in the X and Y positions, with every third position occupied by Gly in fibril-forming collagens. However, in collagen diseases, such as osteogenesis imperfecta, a single Gly may be substituted by another residue. Host-guest studies where the Gly is replaced by various amino acids suggest that the identity of the residue in the Gly position affects the degree of destabilization and the clinical severity of the disease.     

Structural basis of neurophysin hormone specificity: Geometry, polarity, and polarizability in aromatic ring interactions

The structural origins of the specificity of the neurophysin hormone-binding site for an aromatic residue in peptide position 2 were explored by analyzing the binding of a series of peptides in the context of the crystal structure of liganded neurophysin. A new modeling method for describing the van der Waals surface of binding sites assisted in the analysis. Particular attention was paid to the unusually large (5 kcal/mol) difference in binding free energy between Phe and Leu in position 2, a value representing more than three times the maximum expected based on hydrophobicity alone, and additionally remarkable since modeling indicated that the Leu side chain was readily accommodated by the binding pocket. Although evidence was obtained of a weak thermodynamic linkage between the binding interactions of the residue 2 side chain and of the peptide alpha-amino group, two factors are considered central. (1) The bound Leu side chain can establish only one-third of the van der Waals contacts available to a Phe side chain. (2) The bound Phe side chain appears to be additionally stabilized relative to Leu by more favorable dipole and induced dipole interactions with nonaromatic polar and sulfur ligands in the binding pocket, as evidenced by examination of its interactions in the pocket, analysis of the detailed energetics of transfer of Phe and Leu side chains from water to other phases, and comparison with thermodynamic and structural data for the binding of residue 1 side chains in this system. While such polar interactions of aromatic rings have been previously observed, the present results suggest their potential for significant thermodynamic contributions to protein structure and ligand recognition.     

Binding dynamics and energetic insight into the molecular forces driving nucleotide binding by guanylate kinase

Plasmodium deoxyguanylate pathways are an attractive area of investigation for future metabolic and drug discovery studies due to their unique substrate specificities. We investigated the energetic contribution to guanylate kinase substrate binding and the forces underlying ligand recognition. In the range from 20 to 35°C, the thermodynamic profiles displayed marked decrease in binding enthalpy, while the free energy of binding showed little changes. GMP produced a large binding heat capacity change of -356 cal mol(-1) K(-1), indicating considerable conformational changes upon ligand binding. Interestingly, the calculated ΔCp was -32 cal mol(-1) K(-1), indicating that the accessible surface area is not the central change in substrate binding, and that other entropic forces, including conformational changes, are more predominant. The thermodynamic signature for GMP is inconsistent with rigid-body association, while dGMP showed more or less rigid-body association. These binding profiles explain the poor catalytic efficiency and low affinity for dGMP compared with GMP. At low temperature, the ligands bind to the receptor site under the effect of hydrophobic forces. Interestingly, by increasing the temperature, the entropic forces gradually vanish and proceed to a nonfavorable contribution, and the interaction occurs mainly through bonding, electrostatic forces, and van der Waals interactions.     

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.