Structure-based thermodynamic design of peptide ligands: Application to peptide inhibitors of the aspartic protease endothiapepsin

Preparation and application of an economical and environmentally friendly hydrate inhibitor in gas field development

The prevention and control of natural gas hydrates is an important link in ensuring winter production. Traditional thermodynamic inhibitors, like methanol, are commonly utilized due to their low unit costs and pricing, but they come with considerable safety issues when used on-site due to their high toxicity, flammability, and explosive potential. A cost-effective and eco-friendly hydrate inhibitor was created by combining light polyol amine with other ingredients to solve this problem. At a concentration of 30%, the product has a flash point greater than 80°C and a solidification point of -45°C. With success rates of 99% and 100%, respectively, it was used for winter casing pre-injection anti-freezing operations and balancing tank defoamer anti-freezing operations. Experiments have demonstrated the effectiveness of this inhibitor in preventing the formation of natural gas hydrates. In wintertime on-site anti-freezing activities, the projected cost can be substituted for methanol, as it is essentially equivalent to methanol.

The regulatory subunit of PKA-I remains partially structured and undergoes β-aggregation upon thermal denaturation

Background

The regulatory subunit (R) of cAMP-dependent protein kinase (PKA) is a modular flexible protein that responds with large conformational changes to the binding of the effector cAMP. Considering its highly dynamic nature, the protein is rather stable. We studied the thermal denaturation of full-length RIα and a truncated RIα(92-381) that contains the tandem cyclic nucleotide binding (CNB) domains A and B.

Methodology/Principal Findings

As revealed by circular dichroism (CD) and differential scanning calorimetry, both RIα proteins contain significant residual structure in the heat-denatured state. As evidenced by CD, the predominantly α-helical spectrum at 25°C with double negative peaks at 209 and 222 nm changes to a spectrum with a single negative peak at 212–216 nm, characteristic of β-structure. A similar α→β transition occurs at higher temperature in the presence of cAMP. Thioflavin T fluorescence and atomic force microscopy studies support the notion that the structural transition is associated with cross-β-intermolecular aggregation and formation of non-fibrillar oligomers.

Conclusions/Significance

Thermal denaturation of RIα leads to partial loss of native packing with exposure of aggregation-prone motifs, such as the B' helices in the phosphate-binding cassettes of both CNB domains. The topology of the β-sandwiches in these domains favors inter-molecular β-aggregation, which is suppressed in the ligand-bound states of RIα under physiological conditions. Moreover, our results reveal that the CNB domains persist as structural cores through heat-denaturation.

Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

The increasing call for an overall picture of the interactions between the components of a biological system that give rise to the observed function is often summarized by the expression systems biology. Both the interpretative and predictive capabilities of holistic models of biochemical systems, however, depend to a large extent on the level of physico-chemical knowledge of the individual molecular interactions making up the network. This review is focused on the structure-based quantitative characterization of protein-protein interactions, ubiquitous in any biochemical pathway. Recently developed, fast and effective computational methods are reviewed, which allow the assessment of kinetic and thermodynamic features of the association-dissociation processes of protein complexes, both in water soluble and membrane environments. The performance and the accuracy of fast and semi-empirical structure-based methods have reached comparable levels with respect to the classical and more elegant molecular simulations. Nevertheless, the broad accessibility and lower computational cost provide the former methods with the advantageous possibility to perform systems-level analyses including extensive in silico mutagenesis screenings and large-scale structural predictions of multiprotein complexes.

Structural and Functional Properties of Isocitrate Dehydrogenase from the Psychrophilic Bacterium Desulfotalea psychrophila Reveal a Cold-active Enzyme with an Unusual High Thermal Stability

Isocitrate dehydrogenase (IDH) has been studied extensively due to its central role in the Krebs cycle, catalyzing the oxidative NAD(P)(+)-dependent decarboxylation of isocitrate to alpha-ketoglutarate and CO(2). Here, we present the first crystal structure of IDH from a psychrophilic bacterium, Desulfotalea psychrophila (DpIDH). The structural information is combined with a detailed biochemical characterization and a comparative study with IDHs from the mesophilic bacterium Desulfitobacterium hafniense (DhIDH), porcine (PcIDH), human cytosolic (HcIDH) and the hyperthermophilic Thermotoga maritima (TmIDH). DpIDH was found to have a higher melting temperature (T(m)=66.9 degrees C) than its mesophilic homologues and a suboptimal catalytic efficiency at low temperatures. The thermodynamic activation parameters indicated a disordered active site, as seen also for the drastic increase in K(m) for isocitrate at elevated temperatures. A methionine cluster situated at the dimeric interface between the two active sites and a cluster of destabilizing charged amino acids in a region close to the active site might explain the poor isocitrate affinity. On the other hand, DpIDH was optimized for interacting with NADP(+) and the crystal structure revealed unique interactions with the cofactor. The highly acidic surface, destabilizing charged residues, fewer ion pairs and reduced size of ionic networks in DpIDH suggest a flexible global structure. However, strategic placement of ionic interactions stabilizing the N and C termini, and additional ionic interactions in the clasp domain as well as two enlarged aromatic clusters might counteract the destabilizing interactions and promote the increased thermal stability. The structure analysis of DpIDH illustrates how psychrophilic enzymes can adjust their flexibility in dynamic regions during their catalytic cycle without compromising the global stability of the protein.

Bivalent peptides as models for multimeric targets of PDZ domains

PDZ domains are among the most common modules in eukaryotic, including human, genomes. They are found exclusively in large, multidomain cytosolic proteins--often with other domains that belong to a variety of families--and are involved in a plethora of physiological and pathophysiological events. PDZ domains mediate protein-protein interactions by binding to solvent-exposed and extended C-terminal short fragments of membrane-associated proteins, such as receptors and ion channels. Most of what is known about the mechanisms of target binding by PDZ domains is inferred from studies that involve isolated recombinant PDZ domains and short synthetic peptides that represent the targets. These binary systems constitute an obvious oversimplification and disregard factors such as noncanonical modes of binding and enhanced affinity due to multimeric interactions mediated by clusters and oligomers of PDZ-domain-containing proteins. We have tested whether the interaction between a dimeric form of PDZ domain that mimics a functional dimeric guanine nucleotide exchange factor, PDZ-RhoGEF (PDZ-containing RhoA-specific guanine nucleotide exchange factor) or LARG (leukemia-associated RhoA specific guanine nucleotide exchange factor), and a bivalent peptide that mimics the dimer of the plexin B receptor, could enhance the interaction between the two moieties. Peptide dimerization was achieved by cross-linking the N-terminal ends of peptides attached to Wang resin with poly(ethylene glycol) spacers (30-45 Angstroms in length). The interaction of dimeric PDZ domains with dimeric peptides resulted in an up to 20-fold increase in affinity compared to the simple binary system. This is consistent with the notion that multimerization of both receptors and PDZ-containing proteins might constitute an important regulatory mechanism.

The nature of the free energy barriers to two-state folding

We introduce a simple procedure to analyze the temperature dependence of the folding and unfolding rates of two-state proteins. We start from the simple transition-state-like rate expression: k = D(eff)exp(-DeltaG(TS)/RT), in which upper and lower bounds for the intra-chain effective diffusion coefficient (D(eff)) are obtained empirically using the timescales of elementary processes in protein folding. From the changes in DeltaG(TS) as a function of temperature, we calculate enthalpies and heat capacities of activation, together with the more elusive entropies of activation. We then estimate the conformational entropy of the transition state by extrapolation to the temperature at which the solvation entropy vanishes by cancellation between polar and apolar terms. This approach is based on the convergence temperatures for the entropy of solvating apolar (approximately 385 K) and polar groups (approximately 335 K), the assumption that the structural properties of the transition state are somewhere in between the unfolded and folded states, and the established relationship between observed heat capacity and solvent accessibility.1 To circumvent the lack of structural information about transition states, we use the empirically determined heat capacities of activation as constraints to identify the extreme values of the transition state conformational entropy that are consistent with experiment. The application of this simple approach to six two-state folding proteins for which there is temperature-dependent data available in the literature provides important clues about protein folding. For these six proteins, we obtain an average equilibrium cost in conformational entropy of -4.3 cal x mol(-1)K(-1)per residue, which is in close agreement to previous empirical and computational estimates of the same quantity. Furthermore, we find that all these proteins have a conformationally diverse transition state, with more than half of the conformational entropy of the unfolded state. In agreement with predictions from theory and computer simulations, the transition state signals the change from a regime dominated by loss in conformational entropy to one driven by the gain in stabilization free energy (i.e., including protein interactions and solvation effects). Moreover, the height of the barrier is determined by how much stabilization free energy is realized at that point, which is related to the relative contribution of local versus non-local interactions. A remarkable observation is that the fraction of conformational entropy per residue that is present in the transition state is very similar for the six proteins in this study. Based on this commonality, we propose that the observed change in thermodynamic regime is connected to a change in the pattern of structure formation: from one driven by formation of pairwise interactions to one dominated by coupling of the networks of interactions involved in forming the protein core. In this framework, the barrier to two-state folding is crossed when the folding protein reaches a "critical native density" that allows expulsion of remaining interstitial water and consolidation of the core. The principle of critical native density should be general for all two-state proteins, but can accommodate different folding mechanisms depending on the particularities of the structure and sequence.

Thermodynamic characterization of the binding of tetrahydropterins to phenylalanine hydroxylase

Phenylalanine hydroxylase (PAH) is the key enzyme in the catabolism of L-Phe. The natural cofactor of PAH, 6R-tetrahydrobiopterin (BH4), negatively regulates the enzyme activity in addition to being an essential cosubstrate for catalysis. The analogue 6-methyltetrahydropterin (6M-PH4) is effective in catalysis but does not regulate PAH. Here, the thermodynamics of binding of BH4 and 6M-PH4 to human PAH have been studied by isothermal titration calorimetry. At neutral pH and 25 degrees C, BH4 binds to PAH with higher affinity (Kd = 0.75 +/- 0.18 microM) than 6M-PH4 (Kd = 16.5 +/- 2.7 microM). While BH4 binding is a strongly exothermic process (DeltaH = -11.8 +/- 0.4 kcal/mol) accompanied by an entropic penalty (-TDeltaS = 3.4 +/- 0.4 kcal/mol), 6M-PH4 binding is both enthalpically (DeltaH = -3.3 +/- 0.3 kcal/mol) and entropically (-TDeltaS = -3.2 kcal/mol) driven. No significant changes in binding affinity were observed in the 5-35 degrees C temperature range for both pterins at neutral pH, but the enthalpic contribution increased with temperature rendering a heat capacity change (DeltaCp) of -357 +/- 26 cal/mol/K for BH4 and -63 +/- 12 cal/mol/K for 6M-PH4. Protons do not seem to be taken up or released upon pterin binding. Structure-based energetics calculations applied on the molecular dynamics simulated structures of the complexes suggest that in the case of BH4 binding, the conformational rearrangement of the N-terminal tail of PAH contribute with favorable enthalpic and unfavorable entropic contributions to the intrinsic thermodynamic parameters of binding. The entropic penalty is most probably associated to the reduction of conformational flexibility at the protein level and disappears for the L-Phe activated enzyme. The calculated energetic parameters aid to elucidate the molecular mechanism for cofactor recognition and the regulation of PAH by the dihydroxypropyl side chain of BH4.

Energetics of the cooperative assembly of cell division protein FtsZ and the nucleotide hydrolysis switch

FtsZ is the first protein recruited to the bacterial division site, where it forms the cytokinetic Z ring. We have determined the functional energetics of FtsZ assembly, employing FtsZ from the thermophilic Archaea Methanococcus jannaschii bound to GTP, GMPCPP, GDP, or GMPCP, under different solution conditions. FtsZ oligomerizes in a magnesium-insensitive manner. FtsZ cooperatively assembles with magnesium and GTP or GMPCPP into large polymers, following a nucleated condensation polymerization mechanism, under nucleotide hydrolyzing and non-hydrolyzing conditions. The effect of temperature on the critical concentration indicates polymer elongation with an apparent heat capacity change of -800 +/- 100 cal mol-1 K-1 and positive enthalpy and entropy changes, compatible with axial hydrophobic contacts of each FtsZ in the polymer, and predicts optimal polymer stability near 75 degrees C. Assembly entails the binding of one medium affinity magnesium ion and the uptake of one proton per FtsZ. Interestingly, GDP- or GMPCP-liganded FtsZ cooperatively form helically curved polymers, with an elongation only 1-2 kcal mol-1 more unfavorable than the straight polymers formed with nucleotide triphosphate, suggesting a physiological requirement for FtsZ polymerization inhibitors. This GTP hydrolysis switch should provide the basic properties for FtsZ polymer disassembly and its functional dynamics.

The integration of genomic and structural information in the development of high affinity plasmepsin inhibitors

The plasmepsins are key enzymes in the life cycle of the Plasmodium parasites responsible for malaria. Since plasmepsin inhibition leads to parasite death, these enzymes have been acknowledged to be important targets for the development of new antimalarial drugs. The development of effective plasmepsin inhibitors, however, is compounded by their genomic diversity which gives rise not to a unique target for drug development but to a family of closely related targets. Successful drugs will have to inhibit not one but several related enzymes with high affinity. Structure-based drug design against heterogeneous targets requires a departure from the classic 'lock-and-key' paradigm that leads to the development of conformationally constrained molecules aimed at a single target. Drug molecules designed along those principles are usually rigid and unable to adapt to target variations arising from naturally occurring genetic polymorphisms or drug-induced resistant mutations. Heterogeneous targets need adaptive drug molecules, characterised by the presence of flexible elements at specific locations that sustain a viable binding affinity against existing or expected polymorphisms. Adaptive ligands have characteristic thermodynamic signatures that distinguish them from their rigid counterparts. This realisation has led to the development of rigorous thermodynamic design guidelines that take advantage of correlations between the structure of lead compounds and the enthalpic and entropic components of the binding affinity. In this paper, we discuss the application of the thermodynamic approach to the development of high affinity (K(i) - pM) plasmepsin inhibitors. In particular, a family of allophenylnorstatine-based compounds is evaluated for their potential to inhibit a wide spectrum of plasmepsins.

Structural parameterization of the binding enthalpy of small ligands

A major goal in ligand and drug design is the optimization of the binding affinity of selected lead molecules. However, the binding affinity is defined by the free energy of binding, which, in turn, is determined by the enthalpy and entropy changes. Because the binding enthalpy is the term that predominantly reflects the strength of the interactions of the ligand with its target relative to those with the solvent, it is desirable to develop ways of predicting enthalpy changes from structural considerations. The application of structure/enthalpy correlations derived from protein stability data has yielded inconsistent results when applied to small ligands of pharmaceutical interest (MW < 800). Here we present a first attempt at an empirical parameterization of the binding enthalpy for small ligands in terms of structural information. We find that at least three terms need to be considered: (1) the intrinsic enthalpy change that reflects the nature of the interactions between ligand, target, and solvent; (2) the enthalpy associated with any possible conformational change in the protein or ligand upon binding; and, (3) the enthalpy associated with protonation/deprotonation events, if present. As in the case of protein stability, the intrinsic binding enthalpy scales with changes in solvent accessible surface areas. However, an accurate estimation of the intrinsic binding enthalpy requires explicit consideration of long-lived water molecules at the binding interface. The best statistical structure/enthalpy correlation is obtained when buried water molecules within 5-7 A of the ligand are included in the calculations. For all seven protein systems considered (HIV-1 protease, dihydrodipicolinate reductase, Rnase T1, streptavidin, pp60c-Src SH2 domain, Hsp90 molecular chaperone, and bovine beta-trypsin) the binding enthalpy of 25 small molecular weight peptide and nonpeptide ligands can be accounted for with a standard error of +/-0.3 kcal x mol(-1).

L-phenylalanine binding and domain organization in human phenylalanine hydroxylase: a differential scanning calorimetry study

Human phenylalanine hydroxylase (hPAH) is a tetrameric enzyme that catalyzes the hydroxylation of L-phenylalanine (L-Phe) to L-tyrosine; a dysfunction of this enzyme causes phenylketonuria. Each subunit in hPAH contains an N-terminal regulatory domain (Ser2-Ser110), a catalytic domain (Asp112-Arg410), and an oligomerization domain (Ser411-Lys452) including dimerization and tetramerization motifs. Two partially overlapping transitions are seen in differential scanning calorimetry (DSC) thermograms for wild-type hPAH in 0.1 M Na-Hepes buffer, 0.1 M NaCl, pH 7.0. Although these transitions are irreversible, studies on their scan-rate dependence support that the equilibrium thermodynamics analysis is permissible in this case. Comparison with the DSC thermograms for truncated forms of the enzyme, studies on the protein and L-Phe concentration effects on the transitions, and structure-energetic calculations based on a modeled structure support that the thermal denaturation of hPAH occurs in three stages: (i) unfolding of the four regulatory domains, which is responsible for the low-temperature calorimetric transition; (ii) unfolding of two (out of the four) catalytic domains, which is responsible for the high-temperature transition; and (iii) irreversible protein denaturation, which is likely responsible for the observed exothermic distortion in the high-temperature side of the high-temperature transition. Stages 1 and 2 do not appear to be two-state processes. We present an approach to the analysis of ligand effects on DSC transition temperatures, which is based on the general binding polynomial formalism and is not restricted to two-state transitions. Application of this approach to the L-Phe effect on the DSC thermograms for hPAH suggests that (i) there are no binding sites for L-Phe in the regulatory domains; therefore, contrary to the common belief, the activation of PAH by L-Phe seems to be the result of its homotropic cooperative binding to the active sites. (ii) The regulatory domain appears to be involved in cooperativity through its interactions with the catalytic and oligomerization domains; thus, upon regulatory domain unfolding, the cooperativity in the binding of L-Phe to the catalytic domains seems to be lost and the value of the L-Phe concentration corresponding to half-saturation is increased. Overall, our results contribute to the understanding of the conformational stability and the substrate-induced cooperative activation of this important enzyme.

Identification and characterization of allophenylnorstatine-based inhibitors of plasmepsin II, an antimalarial target

Plasmepsin II is a key enzyme in the life cycle of the Plasmodium parasites responsible for malaria, a disease that afflicts more than 300 million individuals annually. Since plasmepsin II inhibition leads to starvation of the parasite, it has been acknowledged as an important target for the development of new antimalarials. In this paper, we identify and characterize high-affinity inhibitors of plasmepsin II based upon the allophenylnorstatine scaffold. The best compound, KNI-727, inhibits plasmepsin II with a K(i) of 70 nM and a 22-fold selectivity with respect to the highly homologous human enzyme cathepsin D. KNI-727 binds to plasmepsin II in a process favored both enthalpically and entropically. At 25 degrees C, the binding enthalpy (DeltaH) is -4.4 kcal/mol and the entropic contribution (-TDeltaS) to the Gibbs energy is -5.56 kcal/mol. Structural stability measurements of plasmepsin II were also utilized to characterize inhibitor binding. High-sensitivity differential scanning calorimetry experiments performed in the absence of inhibitors indicate that, at pH 4.0, plasmepsin II undergoes thermal denaturation at 63.3 degrees C. The structural stability of the enzyme increases with inhibitor concentration in a manner for which the binding energetics of the inhibitor can quantitatively account. The effectiveness of the best compounds in killing the malaria parasite was validated by performing cytotoxicity assays in red blood cells infected with Plasmodium falciparum. EC50s ranging between 6 and 10 microM (3-6 microg/mL) were obtained. These experiments demonstrate the viability of the allophenylnorstatine scaffold in the design of powerful and selective plasmepsin inhibitors.

Isothermal titration calorimetry in drug discovery

Isothermal titration calorimetry (ITC) follows the heat change when a test compound binds to a target protein. It allows precise measurement of affinity. The method is direct, making interpretation facile, because there is no requirement for competing molecules. Titration in the presence of other ligands rapidly provides information on the mechanism of action of the test compound, identifying the intermolecular complexes that are relevant for structure-based design. Calorimetry allows measurement of stoichiometry and so evaluation of the proportion of the sample that is functional. ITC can characterize protein fragments and catalytically inactive mutant enzymes. It is the only technique which directly measures the enthalpy of binding (delta H degree). Interpretation of delta H degree and its temperature dependence (delta Cp) is usually qualitative, not quantitative. This is because of complicated contributions from linked equilibria and a single change in structure giving modification of several physicochemical properties. Measured delta H degree values allow characterization of proton movement linked to the association of protein and ligand, giving information on the ionization of groups involved in binding. Biochemical systems characteristically exhibit enthalpy-entropy compensation where increased bonding is offset by an entropic penalty, reducing the magnitude of change in affinity. This also causes a lack of correlation between the free energy of binding (delta G degree) and delta H degree. When characterizing structure-activity relationships (SAR), most groups involved in binding can be detected as contributing to delta H degree, but not to affinity. Large enthalpy changes may reflect a modified binding mode, or protein conformation changes. Thus, delta H degree values may highlight a potential discontinuity in SAR, so that experimental structural data are likely to be particularly valuable in molecular design.

The binding energetics of first- and second-generation HIV-1 protease inhibitors: implications for drug design

KNI-764 is a powerful HIV-1 protease inhibitor with a reported low susceptibility to the effects of protease mutations commonly associated with drug resistance. In this paper the binding thermodynamics of KNI-764 to the wild-type and drug-resistant mutant V82F/I84V are presented and the results compared to those obtained with existing clinical inhibitors. KNI-764 binds to the wild-type HIV-1 protease with very high affinity (3.1 x 10(10) M(-1) or 32 pM) in a process strongly favored by both enthalpic and entropic contributions to the Gibbs energy of binding (Delta G = -RTlnK(a)). When compared to existing clinical inhibitors, the binding affinity of KNI-764 is about 100 fold higher than that of indinavir, saquinavir, and nelfinavir, but comparable to that of ritonavir. Unlike the existing clinical inhibitors, which bind to the protease with unfavorable or only slightly favorable enthalpy changes, the binding of KNI-764 is strongly exothermic (-7.6 kcal/mol). The resistant mutation V82F/I84V lowers the binding affinity of KNI-764 26-fold, which can be accounted almost entirely by a less favorable binding enthalpy to the mutant. Since KNI-764 binds to the wild type with extremely high affinity, even after a 26-fold decrease, it still binds to the resistant mutant with an affinity comparable to that of other inhibitors against the wild type. These results indicate that the effectiveness of this inhibitor against the resistant mutant is related to two factors: extremely high affinity against the wild type achieved by combining favorable enthalpic and entropic interactions, and a mild effect of the protease mutation due to the presence of flexible structural elements at critical locations in the inhibitor molecule. The conclusions derived from the HIV-1 protease provide important thermodynamic guidelines that can be implemented in general drug design strategies.

Calculation of entropy changes in biological processes: folding, binding, and oligomerization

Changes in configurational entropy represent one of the major contributions to the thermodynamics of folding, binding, and oligomerization. Methods have been developed to estimate changes in the entropy of the backbone and side chains, and for the loss of translational entropy. These methods have been used in combination with empirical methods that provide estimates of the changes in entropy of solvation as well as estimates of the changes of enthalpy. The results of such calculations are in excellent agreement with experimentally observed values.

Docking molecules by families to increase the diversity of hits in database screens: computational strategy and experimental evaluation

Molecular docking programs screen chemical databases for novel ligands that fit protein binding sites. When one compound fits the site well, close analogs typically do the same. Therefore, many of the compounds that are found in such screens resemble one another. This reduces the variety and novelty of the compounds suggested. In an attempt to increase the diversity of docking hit lists, the Available Chemicals Directory was grouped into families of related structures. All members of every family were docked and scored, but only the best scoring molecule of a high-ranking family was allowed in the hit list. The identity and scores of the other members of these families were recorded as annotations to the best family member, but they were not independently ranked. This family-based docking method was compared with molecule-by-molecule docking in screens against the structures of thymidylate synthase, dihydrofolate reductase (DHFR), and the cavity site of the mutant T4 lysozyme Leu99 --> Ala (L99A). In each case, the diversity of the hit list increased, and more families of known ligands were found. To investigate whether the newly identified hits were sensible, we tested representative examples experimentally for binding to L99A and DHFR. Of the six compounds tested against L99A, five bound to the internal cavity. Of the seven compounds tested against DHFR, six inhibited the enzyme with apparent K(i) values between 0.26 and 100 microM. The segregation of potential ligands into families of related molecules is a simple technique to increase the diversity of candidates suggested by database screens. The general approach should be applicable to most docking methods. Proteins 2001;42:279-293.

Comparison of binding energies of SrcSH2-phosphotyrosyl peptides with structure-based prediction using surface area based empirical parameterization

The prediction of binding energies from the three-dimensional (3D) structure of a protein-ligand complex is an important goal of biophysics and structural biology. Here, we critically assess the use of empirical, solvent-accessible surface area-based calculations for the prediction of the binding of Src-SH2 domain with a series of tyrosyl phosphopeptides based on the high-affinity ligand from the hamster middle T antigen (hmT), where the residue in the pY+ 3 position has been changed. Two other peptides based on the C-terminal regulatory site of the Src protein and the platelet-derived growth factor receptor (PDGFR) are also investigated. Here, we take into account the effects of proton linkage on binding, and test five different surface area-based models that include different treatments for the contributions to conformational change and protein solvation. These differences relate to the treatment of conformational flexibility in the peptide ligand and the inclusion of proximal ordered solvent molecules in the surface area calculations. This allowed the calculation of a range of thermodynamic state functions (deltaCp, deltaS, deltaH, and deltaG) directly from structure. Comparison with the experimentally derived data shows little agreement for the interaction of SrcSH2 domain and the range of tyrosyl phosphopeptides. Furthermore, the adoption of the different models to treat conformational change and solvation has a dramatic effect on the calculated thermodynamic functions, making the predicted binding energies highly model dependent. While empirical, solvent-accessible surface area based calculations are becoming widely adopted to interpret thermodynamic data, this study highlights potential problems with application and interpretation of this type of approach. There is undoubtedly some agreement between predicted and experimentally determined thermodynamic parameters: however, the tolerance of this approach is not sufficient to make it ubiquitously applicable.

Side-chain conformational entropy in protein unfolded states

The largest force disfavoring the folding of a protein is the loss of conformational entropy. A large contribution to this entropy loss is due to the side-chains, which are restricted, although not immobilized, in the folded protein. In order to accurately estimate the loss of side-chain conformational entropy that occurs upon folding it is necessary to have accurate estimates of the amount of entropy possessed by side-chains in the ensemble of unfolded states. A new scale of side-chain conformational entropies is presented here. This scale was derived from Monte Carlo computer simulations of small peptide models. It is demonstrated that the entropies are independent of host peptide length. This new scale has the advantage over previous scales of being more precise with low standard errors. Better estimates are obtained for long (e.g., Arg and Lys) and rare (e.g., Trp and Met) side-chains. Excellent agreement with previous side-chain entropy scales is achieved, indicating that further advancements in accuracy are likely to be small at best. Strikingly, longer side-chains are found to possess a smaller fraction of the theoretical maximum entropy available than short side-chains. This indicates that rotations about torsions after chi(2) are significantly affected by side-chain interactions with the polypeptide backbone. This finding invalidates previous assumptions about side-chain-backbone interactions. Proteins 2000;40:443-450.

Thermodynamic characterization of the interaction between TRAF2 and tumor necrosis factor receptor peptides by isothermal titration calorimetry

The tumor necrosis factor receptor (TNFR) superfamily can induce diverse biological effects, including cell survival, proliferation, differentiation, and apoptosis. The major signal transducers for TNFRs are the family of TNF receptor associated factors (TRAFs). The direct interaction between TRAFs and the intracellular tails of TNFRs is the first step of this signal relay process. Structural studies have revealed a trimeric nature of TRAF2 and a symmetrical mode of receptor binding, suggesting the involvement of trivalent TNFR2-receptor interaction in the signal transduction. In this study, using isothermal titration calorimetry (ITC), we report thermodynamic characterization of the interaction between TRAF2 and monomeric peptide sequences from TNFR members, including TNFR2, CD40, CD30, Ox40, and 4-1BB, and the Epstein-Barr virus (EBV)-transforming protein, latent infection membrane protein-1 (LMP1). The dissociation constants of the interaction were shown to range between 40 microM and 1.9 mM, which are substantially weaker than most protein-peptide interactions. The interaction is entirely driven by exothermic enthalpy, consistent with the abundance of polar contacts. The enthalpy of the interaction has a significant temperature dependence (DeltaCp = -245 cal/mol small middle dotK). The unfavorable entropy in the interaction and the comparison with structural energetics calculations suggest the involvement of conformational rearrangement in the interaction. The low affinity of TRAF2 to monomeric receptor peptides further supports the importance of avidity contribution in TRAF2 recruitment by these receptors upon ligand-induced trimerization or higher order oligomerization.

X-ray structures of small ligand-FKBP complexes provide an estimate for hydrophobic interaction energies

A new crystal form of native FK506 binding protein (FKBP) has been obtained which has proved useful in ligand binding studies. Three different small molecule ligand complexes and the native enzyme have been determined at higher resolution than 2.0 A. Dissociation constants of the related small molecule ligands vary from 20 mM for dimethylsulphoxide to 200 microM for tetrahydrothiophene 1-oxide. Comparison of the four available crystal structures shows that the protein structures are identical to within experimental error, but there are differences in the water structure in the active site. Analysis of the calculated buried surface areas of these related ligands provides an estimated van der Waals contribution to the binding energy of -0.5 kJ/A(2) for non-polar interactions between ligand and protein.

Tricyclic farnesyl protein transferase inhibitors: crystallographic and calorimetric studies of structure-activity relationships

Crystallographic and thermodynamic studies of farnesyl protein transferase (FPT) complexed with novel tricyclic inhibitors provide insights into the observed SAR for this unique class of nonpeptidic FPT inhibitors. The crystallographic structures reveal a binding pattern conserved across the mono-, di-, and trihalogen series. In the complexes, the tricycle spans the FPT active site cavity and interacts with both protein atoms and the isoprenoid portion of bound farnesyl diphosphate. An amide carbonyl, common to the tricyclic compounds described here, participates in a water-mediated hydrogen bond to the protein backbone. Ten high-resolution crystal structures of inhibitors complexed with FPT are reported. Included are crystallographic data for FPT complexed with SCH 66336, a compound currently undergoing clinical trials as an anticancer agent (SCH 66336, 4-[2-[4-(3,10-dibromo-8-chloro-6,11-dihydro-5H-benzo[5, 6]cyclohepta[1, 2-b]pyridin-11-yl)-1-piperidinyl]-2-oxoethyl]-1-piperidinecarbo xamide ). Thermodynamic binding parameters show favorable enthalpies of complex formation and small net entropic contributions as observed for 4-[2-[4-(3,10-dibromo-8-chloro-6,11-dihydro-11H-benzo[5, 6]cyclohepta[1, 2-b]pyridin-11-ylidene)-1-piperidinyl]-2-oxoethyl]pyridine N-oxide where DeltaH degrees bind = -12.5 kcal/mol and TDeltaS degrees bind = -1.5 kcal/mol.

Ligand solvation in molecular docking

Solvation plays an important role in ligand-protein association and has a strong impact on comparisons of binding energies for dissimilar molecules. When databases of such molecules are screened for complementarity to receptors of known structure, as often occurs in structure-based inhibitor discovery, failure to consider ligand solvation often leads to putative ligands that are too highly charged or too large. To correct for the different charge states and sizes of the ligands, we calculated electrostatic and non-polar solvation free energies for molecules in a widely used molecular database, the Available Chemicals Directory (ACD). A modified Born equation treatment was used to calculate the electrostatic component of ligand solvation. The non-polar component of ligand solvation was calculated based on the surface area of the ligand and parameters derived from the hydration energies of apolar ligands. These solvation energies were subtracted from the ligand-receptor interaction energies. We tested the usefulness of these corrections by screening the ACD for molecules that complemented three proteins of known structure, using a molecular docking program. Correcting for ligand solvation improved the rankings of known ligands and discriminated against molecules with inappropriate charge states and sizes.

Crystal structure of the abl-SH3 domain complexed with a designed high-affinity peptide ligand: implications for SH3-ligand interactions

The Abl-SH3 domain is implicated in negative regulation of the Abl kinase by mediating protein-protein interactions. High-affinity SH3 ligands could compete for these interactions and specifically activate the Abl kinase, providing control and a better understanding of the molecular interactions that underlie diseases where SH3 domains are involved. The p41 peptide (APSYSPPPPP) is a member of a group of peptide ligands designed to bind specifically the Abl-SH3 domain. It binds to Abl-SH3 with a Kd of 1.5 microM, whereas its affinity for the Fyn-SH3 domain is 273 microM. We have determined the crystal structure of the Abl-SH3 domain in complex with the high-affinity peptide ligand p41 at 1.6 A resolution. In the crystal structure, this peptide adopts a polyproline type II helix conformation through residue 5 to 10, and it binds in type I orientation to the Abl-SH3 domain. The tyrosine side-chain in position 4 of the peptide is hydrogen bonded to two residues in the RT-loop of the Abl-SH3 domain. The tight fit of this side-chain into the RT-loop pocket is enhanced by conformational adjustment of the main chain at position 5. The SH3 ligand peptides can be divided into two distinct parts. The N-terminal part binds to the SH3 domain in the region formed by the valley between the nSrc and RT-loops. It determines the specificity for different SH3 domains. The C-terminal part adopts a polyproline type II helix conformation. This binds to a well-conserved hydrophobic surface of the SH3 domain. Analysis of two "half"-peptides, corresponding to these ligand parts, shows that both are essential components for strong binding to the SH3 domains. The crystal structure of the Abl-SH3:p41 complex explains the high affinity and specificity of the p41 peptide towards the Abl-SH3 domain, and reveals principles that will be exploited for future design of small, high-affinity ligands to interfere efficiently with the in vivo regulation of Abl kinase activity.

Structure-based drug design

Most of the techniques used in structure-based drug design have experienced significant improvements in the past few years, resulting in a marked enhancement of the speed and the efficacy of this approach. At the same time, it was thought that the future of drug design lay in strategies involving solely combinatorial chemistry. It is becoming evident, however, that the development of future drugs will use a combination of methods that will contain a major component of structure-based design.

Molecular basis of resistance to HIV-1 protease inhibition: a plausible hypothesis

The binding thermodynamics of the HIV-1 protease inhibitor acetyl pepstatin and the substrate Val-Ser-Gln-Asn-Tyr-Pro-Ile-Val-Gln, corresponding to one of the cleavage sites in the gag, gag-pol polyproteins, have been measured by direct microcalorimetric analysis. The results indicate that the binding of the peptide substrate or peptide inhibitor is entropically driven; i.e., it is characterized by an unfavorable enthalpy and a favorable entropy change, in agreement with a structure-based thermodynamic analysis based upon an empirical parameterization of the energetics. Dissection of the binding enthalpy indicates that the intrinsic interactions are favorable and that the unfavorable enthalpy originates from the energy cost of rearranging the flap region in the protease molecule. In addition, the binding is coupled to a negative heat capacity change. The dominant binding force is the increase in solvent entropy that accompanies the burial of a significant hydrophobic surface. Comparison of the binding energetics obtained for the substrate with that obtained for synthetic nonpeptide inhibitors indicates that the major difference is in the magnitude of the conformational entropy change. In solution, the peptide substrate has a higher flexibility than the synthetic inhibitors and therefore suffers a higher conformational entropy loss upon binding. This higher entropy loss accounts for the lower binding affinity of the substrate. On the other hand, due to its higher flexibility, the peptide substrate is more amenable to adapt to backbone rearrangements or subtle conformational changes induced by mutations in the protease. The synthetic inhibitors are less flexible, and their capacity to adapt is more restricted. The expected result is a more pronounced effect of mutations on the binding affinity of the synthetic inhibitors. On the basis of the thermodynamic differences in the mode of binding of substrate and synthetic inhibitors, it appears that a key factor to understanding resistance is given by the relative balance of the different forces that contribute to the binding free energy and, in particular, the balance between conformational and solvation entropy.