Structural bases of lectin-carbohydrate affinities: comparison with protein-folding energetics.

We have made a comparative structure based analysis of the thermodynamics of lectin-carbohydrate (L-C) binding and protein folding. Examination of the total change in accessible surface area in those processes revealed a much larger decrease in free energy per unit of area buried in the case of L-C associations. According to our analysis, this larger stabilization of L-C interactions arises from a more favorable enthalpy of burying a unit of polar surface area, and from higher proportions of polar areas. Hydrogen bonds present at 14 L-C interfaces were identified, and their overall characteristics were compared to those reported before for hydrogen bonds in protein structures. Three major factors might explain why polar-polar interactions are stronger in L-C binding than in protein folding: (1) higher surface density of hydrogen bonds; (2) better hydrogen-bonding geometry; (3) larger proportion of hydrogen bonds involving charged groups. Theoretically, the binding entropy can be partitioned into three main contributions: entropy changes due to surface desolvation, entropy losses arising from freezing rotatable bonds, and entropic effects that result from restricting translation and overall rotation motions. These contributions were estimated from structural information and added up to give calculated binding entropies. Good correlation between experimental and calculated values was observed when solvation effects were treated according to a parametrization developed by other authors from protein folding studies. Finally, our structural parametrization gave calculated free energies that deviate from experimental values by 1.1 kcal/mol on the average; this amounts to an uncertainty of one order of magnitude in the binding constant.     

Insufficient hydrogen-bond desolvation and prion-related disease.

A structuring and eventual exclusion of water surrounding backbone hydrogen bonds takes place during protein folding as hydrophobic residues cluster around such bonds. Taken as an average over all hydrogen bonds, the extent of desolvation is nearly a constant of motion, as revealed by re-examination of the longest all-atom trajectory with explicit solvent [Y. Duan & P. A. Kollman (1998) Science 282, 740]. Furthermore, this extent of desolvation is preserved across native soluble proteins, except for cellular prion proteins. Thus, a physico-chemical picture of prion-related disease emerges. The epitope for protein-X binding, the region undergoing vast conformational change and the trigger and locker for this change are inferred from the location of under-desolvated hydrogen bonds in the cellular prion protein.     

Predicting protein interaction sites: binding hot-spots in protein-protein and protein-ligand interfaces

MOTIVATION: Protein assemblies are currently poorly represented in structural databases and their structural elucidation is a key goal in biology. Here we analyse clefts in protein surfaces, likely to correspond to binding 'hot-spots', and rank them according to sequence conservation and simple measures of physical properties including hydrophobicity, desolvation, electrostatic and van der Waals potentials, to predict which are involved in binding in the native complex. RESULTS: The resulting differences between predicting binding-sites at protein-protein and protein-ligand interfaces are striking. There is a high level of prediction accuracy (< or =93%) for protein-ligand interactions, based on the following attributes: van der Waals potential, electrostatic potential, desolvation and surface conservation. Generally, the prediction accuracy for protein-protein interactions is lower, with the exception of enzymes. Our results show that the ease of cleft desolvation is strongly predictive of interfaces and strongly maintained across all classes of protein-binding interface.     

Determinants of ligand binding to cAMP-dependent protein kinase

Protein kinases are essential for the regulation of cellular growth and metabolism. Since their dysfunction leads to debilitating diseases, they represent key targets for pharmaceutical research. The rational design of kinase inhibitors requires an understanding of the determinants of ligand binding to these proteins. In the present study, a theoretical model based on continuum electrostatics and a surface-area-dependent nonpolar term is used to calculate binding affinities of balanol derivatives, H-series inhibitors, and ATP analogues toward the catalytic subunit of cAMP-dependent protein kinase (cAPK or protein kinase A). The calculations reproduce most of the experimental trends and provide insight into the driving forces responsible for binding. Nonpolar interactions are found to govern protein-ligand affinity. Hydrogen bonds represent a negligible contribution, because hydrogen bond formation in the complex requires the desolvation of the interacting partners. However, the binding affinity is decreased if hydrogen-bonding groups of the ligand remain unsatisfied in the complex. The disposition of hydrogen-bonding groups in the ligand is therefore crucial for binding specificity. These observations should be valuable guides in the design of potent and specific kinase inhibitors.     

Prediction of Protein-Protein Interaction Sites Using Electrostatic Desolvation Profiles

Protein-protein complex formation involves removal of water from the interface region. Surface regions with a small free energy penalty for water removal or desolvation may correspond to preferred interaction sites. A method to calculate the electrostatic free energy of placing a neutral low-dielectric probe at various protein surface positions has been designed and applied to characterize putative interaction sites. Based on solutions of the finite-difference Poisson equation, this method also includes long-range electrostatic contributions and the protein solvent boundary shape in contrast to accessible-surface-area-based solvation energies. Calculations on a large set of proteins indicate that in many cases (>90%), the known binding site overlaps with one of the six regions of lowest electrostatic desolvation penalty (overlap with the lowest desolvation region for 48% of proteins). Since the onset of electrostatic desolvation occurs even before direct protein-protein contact formation, it may help guide proteins toward the binding region in the final stage of complex formation. It is interesting that the probe desolvation properties associated with residue types were found to depend to some degree on whether the residue was outside of or part of a binding site. The probe desolvation penalty was on average smaller if the residue was part of a binding site compared to other surface locations. Applications to several antigen-antibody complexes demonstrated that the approach might be useful not only to predict protein interaction sites in general but to map potential antigenic epitopes on protein surfaces.     

New insights into the molecular basis of lectin-carbohydrate interactions: A calorimetric and structural study of the association of hevein to oligomers of N-acetylglucosamine

Isothermal titration calorimetry was used to characterize thermodynamically the association of hevein, a lectin from the rubber tree latex, with the dimer and trimer of N-acetylglucosamine (GlcNAc). Considering the changes in polar and apolar accessible surface areas due to complex formation, we found that the experimental binding heat capacities can be reproduced adequately by means of parameters used in protein-unfolding studies. The same conclusion applies to the association of the lectin concanavalin A with methyl-α-mannopyranoside. When reduced by the polar area change, binding enthalpy values show a minimal dispersion around 100°C. These findings resemble the convergence observed in protein-folding events; however, the average of reduced enthalpies for lectin-carbohydrate associations is largely higher than that for the folding of proteins. Analysis of hydrogen bonds present at lectin-carbohydrate interfaces revealed geometries closer to ideal values than those observed in protein structures. Thus, the formation of more energetic hydrogen bonds might well explain the high association enthalpies of lectin-carbohydrate systems. We also have calculated the energy associated with the desolvation of the contact zones in the binding molecules and from it the binding enthalpy in vacuum. This latter resulted 20% larger than the interaction energy derived from the use of potential energy functions. Proteins 29:467–477, 1997. © 1997 Wiley-Liss, Inc.     

Synchronized conformational fluctuations and binding site desolvation during molecular recognition

Binding site desolvation is a poorly understood prerequisite to ligand binding. Although structural fluctuations may be expected to have an important role, little is known about which fluctuations are important or the mechanism by which they promote desolvation. This investigation examines whether and how specific structural fluctuations contribute to desolvation of the ligand binding site in glycopeptide antibiotics. Backbone peptide group rotations in vancomycin, known to occur by experimental observation, were examined in this work with a two-dimensional adaptive umbrella sampling molecular dynamics simulation technique. Results indicate that energetic barriers to rotation are relatively small for two of the peptide groups intimately involved in ligand recognition. When they occur, these rotations strip water molecules away from key hydrogen bond donors and simultaneously cause significant distortions in the macrocyclic rings of the antibiotic that force water into and out of the binding site. Both events are intricately synchronized on the molecular level and have consequences that are clearly necessary to prepare the binding site for receiving a ligand. These results suggest that previously reported observations concerning structural dynamics and binding kinetics in these compounds are mechanistically linked, and they illustrate a heretofore unrecognized degree of preorganization, complexity, and synchronization that may be involved in specific molecular recognition. They also suggest that strategies for increasing antibiotic affinity through covalent dimerization may be counterproductive.     

Computation of electrostatic complements to proteins: a case of charge stabilized binding.

Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.     

Conformational selection through electrostatics: Free energy simulations of GTP and GDP binding to archaeal initiation factor 2

Archaeal Initiation Factor 2 is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and their dependence on the conformational state, molecular dynamics free energy simulations were performed. The ON state specificity was predicted to be weak, with a GTP/GDP binding free energy difference of -1 kcal/mol, favoring GTP. The OFF state specificity is larger, 4 kcal/mol, favoring GDP. The overall effects result from a competition among many interactions in several complexes. To interpret them, we use a simpler, dielectric continuum model. Several effects are robust with respect to the model details. Both nucleotides have a net negative charge, so that removing them from solvent into the binding pocket carries a desolvation penalty, which is large for the ON state, and strongly disfavors GTP binding compared to GDP. Short-range interactions between the additional GTP phosphate group and ionized sidechains in the binding pocket offset most, but not all of the desolvation penalty; more distant groups also contribute significantly, and the switch 1 loop only slightly. The desolvation penalty is lower for the more open, wetter OFF state, and the GTP/GDP difference much smaller. Short-range interactions in the binding pocket and with more distant groups again make a significant contribution. Overall, the simulations help explain how conformational selection is achieved with a single phosphate group.     

Using physical potentials and learned models to distinguish native binding interfaces from de novo designed interfaces that do not bind

Protein–protein interactions are a fundamental aspect of many biological processes. The advent of recombinant protein and computational techniques has allowed for the rational design of proteins with novel binding capabilities. It is therefore desirable to predict which designed proteins are capable of binding in vitro. To this end, we have developed a learned classification model that combines energetic and non‐energetic features. Our feature set is adapted from specialized potentials for aromatic interactions, hydrogen bonds, electrostatics, shape, and desolvation. A binding model built on these features was initially developed for CAPRI Round 21, achieving top results in the independent assessment. Here, we present a more thoroughly trained and validated model, and compare various support‐vector machine kernels. The Gaussian kernel model classified both high‐resolution complexes and designed nonbinders with 79–86% accuracy on independent test data. We also observe that multiple physical potentials for dielectric‐dependent electrostatics and hydrogen bonding contribute to the enhanced predictive accuracy, suggesting that their combined information is much greater than that of any single energetics model. We also study the change in predictive performance as the model features or training data are varied, observing unusual patterns of prediction in designed interfaces as compared with other data types. Proteins 2013; 81:1919–1930. © 2013 Wiley Periodicals, Inc.     

Algorithms for computational solvent mapping of proteins

Computational mapping methods place molecular probes (small molecules or functional groups) on a protein surface to identify the most favorable binding positions by calculating an interaction potential. We have developed a novel computational mapping program called CS-Map (computational solvent mapping of proteins), which differs from earlier mapping methods in three respects: (i) it initially moves the ligands on the protein surface toward regions with favorable electrostatics and desolvation, (ii) the final scoring potential accounts for desolvation, and (iii) the docked ligand positions are clustered, and the clusters are ranked on the basis of their average free energies. To understand the relative importance of these factors, we developed alternative algorithms that use the DOCK and GRAMM programs for the initial search. Because of the availability of experimental solvent mapping data, lysozyme and thermolysin are considered as test proteins. Both DOCK and GRAMM speed up the initial search, and the combined algorithms yield acceptable mapping results. However, the DOCK-based approaches place the consensus site farther from its experimentally determined position than CS-Map, primarily because of the lack of a solvation term in the initial search. The GRAMM-based program also finds the correct consensus site for thermolysin. We conclude that good sampling is the most important requirement for successful mapping, but accounting for desolvation and clustering of ligand positions also help to reduce the number of false positives.     

Wrapping Effects within a Proposed Function-Rescue Strategy for the Y220C Oncogenic Mutation of Protein p53

Soluble proteins must protect their structural integrity from water attack by wrapping interactions which imply the clustering of nonpolar residues around the backbone hydrogen bonds. Thus, poorly wrapped hydrogen bonds constitute defects which have been identified as promoters of protein associations since they favor the removal of hydrating molecules. More specifically, a recent study of our group has shown that wrapping interactions allow the successful identification of protein binding hot spots. Additionally, we have also shown that drugs disruptive of protein-protein interfaces tend to mimic the wrapping behavior of the protein they replace. Within this context, in this work we study wrapping three body interactions related to the oncogenic Y220C mutation of the tumor suppressor protein p53. Our computational results rationalize the oncogenic nature of the Y220C mutation, explain the binding of a drug-like molecule already designed to restore the function of p53 and provide clues to help improve this function-rescue strategy and to apply in other drug design or re-engineering techniques.     

Desolvation is a likely origin of robust enthalpic barriers to protein folding

Experimental data from global analyses of temperature (T) and denaturant dependence of the folding rates of small proteins led to an intrinsic enthalpic folding barrier hypothesis: to a good approximation, the T-dependence of folding rate under constant native stability conditions is Arrhenius. Furthermore, for a given protein, the slope of isostability folding rate versus 1/T is essentially independent of native stability. This hypothesis implies a simple relationship between chevron and Eyring plots of folding that is easily discernible when both sets of rates are expressed as functions of native stability. Using experimental data in the literature, we verify the predicted chevron-Eyring relationship for 14 proteins and determine their intrinsic enthalpic folding barriers, which vary approximately from 15 kcal/mol to 40 kcal/mol for different proteins. These enthalpic barriers do not appear to correlate with folding rates, but they exhibit correlation with equilibrium unfolding enthalpy at room temperature. Intrinsic enthalpic barriers with similarly high magnitudes apply as well to at least two cases of peptide-peptide and peptide-protein association, suggesting that these barriers are a hallmark of certain general and fundamental kinetic processes during folding and binding. Using a class of explicit-chain C(alpha) protein models with constant elementary enthalpic desolvation barriers between C(alpha) positions, we show that small microscopic pairwise desolvation barriers, which are a direct consequence of the particulate nature of water, can act cooperatively to give rise to a significant overall enthalpic barrier to folding. This theoretical finding provides a physical rationalization for the high intrinsic enthalpic barriers in protein folding energetics. Ramifications of entropy-enthalpy compensation in hydrophobic association for the height of enthalpic desolvation barrier are discussed.     

Interpreting trends in the binding of cyclic ureas to HIV-1 protease

The design of new HIV protease inhibitors requires an improved understanding of the physical basis of inhibitor/protein binding. Here, the binding affinities of seven aliphatic cyclic ureas to HIV-1 protease are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed U-shaped trend of binding free energy as a function of aliphatic chain length. Interestingly, the decrease in affinity for the longest chains is attributable primarily to the energy cost of partly desolvating charged aspartic and arginine groups at the mouths of the active site. Even aliphatic chains too short to contact these charged groups directly are subject to considerable desolvation penalties. We are not aware of other systems where binding affinity trends have been attributed to long-ranged electrostatic desolvation of ionized groups. A generalized Born/surface area solvation model yields a much smaller change in desolvation energy with chain length and, therefore, does not reproduce the experimental binding affinity trends. This result suggests that the generalized Born model should be used with caution for complex, partly desolvated systems like protein binding sites. We also find that changing the assumed protonation state of the active site aspartyl dyad significantly affects the computed binding affinity trends. The protonation state of the aspartyl dyad in the presence of cyclic ureas is discussed in light of the observation that the monoprotonated state reproduces the experimental results best.     

Free energy landscapes of encounter complexes in protein-protein association

We report the computer generation of a high-density map of the thermodynamic properties of the diffusion-accessible encounter conformations of four receptor-ligand protein pairs, and use it to study the electrostatic and desolvation components of the free energy of association. Encounter complex conformations are generated by sampling the translational/rotational space of the ligand around the receptor, both at 5-A and zero surface-to-surface separations. We find that partial desolvation is always an important effect, and it becomes dominant for complexes in which one of the reactants is neutral or weakly charged. The interaction provides a slowly varying attractive force over a small but significant region of the molecular surface. In complexes with no strong charge complementarity this region surrounds the binding site, and the orientation of the ligand in the encounter conformation with the lowest desolvation free energy is similar to the one observed in the fully formed complex. Complexes with strong opposite charges exhibit two types of behavior. In the first group, represented by barnase/barstar, electrostatics exerts strong orientational steering toward the binding site, and desolvation provides some added adhesion within the local region of low electrostatic energy. In the second group, represented by the complex of kallikrein and pancreatic trypsin inhibitor, the overall stability results from the rather nonspecific electrostatic attraction, whereas the affinity toward the binding region is determined by desolvation interactions.     

Thermodynamic analysis of binding between mouse major urinary protein-I and the pheromone 2-sec-butyl-4,5-dihydrothiazole

The mouse pheromone 2-sec-butyl-4,5-dihydrothiazole (SBT) binds to an occluded, nonpolar cavity in the mouse major urinary protein-I (MUP-I). The thermodynamics of this interaction have been characterized using isothermal titration calorimetry (ITC). MUP-I-SBT binding is accompanied by a large favorable enthalpy change (DeltaH = -11.2 kcal/mol at 25 degrees C), an unfavorable entropy change (-TDeltaS = 2.8 kcal/mol at 25 degrees C), and a negative heat capacity change [DeltaC(p)() = -165 cal/(mol K)]. Thermodynamic analysis of binding between MUP-I and several 2-alkyl-4,5-dihydrothiazole ligands indicated that the alkyl chain contributes more favorably to the enthalpy and less favorably to the entropy of binding than would be expected on the basis of the hydrophobic desolvation of short-chain alcohols. However, solvent transfer experiments indicated that desolvation of SBT is accompanied by a net unfavorable change in enthalpy (DeltaH = +1.0 kcal/mol) and favorable change in entropy (-TDeltaS = -1.8 kcal/mol). These results are discussed in terms of the possible physical origins of the binding thermodynamics, including (1) hydrophobic desolvation of both the protein and the ligand, (2) formation of a buried water-mediated hydrogen bond network between the protein and ligand, (3) formation of strong van der Waals interactions, and (4) changes in the structure, dynamics, and/or hydration of the protein upon binding.     

Oncogenic mutations and packing defects in protein structure

Oncogenic mutations in expressed proteins are of primary interest to understand tumor formation but their structural consequences bearing on protein function are not clearly understood. In this contribution I report on two illustrative examples, p21ras and p57, revealing that such mutations have an effect on specific structural deficiencies in the packing of the protein structure, i. e., on backbone hydrogen bonds insufficiently shielded from water attack. These structural deficiencies in the wild type are typically "corrected intermolecularly" by protein complexation or protein-ligand association. However, in the oncogenic mutants, these binding signals are partially or completely suppressed: the mutated residues properly wrap or desolvate the hydrogen bonds intramolecularly. Thus, the interactivity of the proteins becomes impaired: their binding affinity decreases sharply, as there is no thermodynamic benefit from removing water surrounding properly desolvated hydrogen bonds. The results, specialized for p21ras and p53, reveal how oncogenic mutations determine a hindrance to GAP-induced hydrolysis (p21) and decrease binding affinity for DNA (p53). Furthermore, the oncogenic potential of mutations in residues not directly engaged in the interface electrostatics is assessed. The results suggest that a high sensitivity of structural defects to genetic accident might be a necessary condition to establish the existence of a proto-oncogene, an angle that merits a systematic study.     

Oncogenic Mutations and Packing Defects in Protein Structure

Abstract

Oncogenic mutations in expressed proteins are of primary interest to understand tumor formation but their structural consequences bearing on protein function are not clearly understood. In this contribution I report on two illustrative examples, p21ras and p57, revealing that such mutations have an effect on specific structural deficiencies in the packing of the protein structure, i. e., on backbone hydrogen bonds insufficiently shielded from water attack. These structural deficiencies in the wild type are typically “corrected intermolecu- larly” by protein complexation or protein-ligand association. However, in the oncogenic mutants, these binding signals are partially or completely suppressed: the mutated residues properly wrap or desolvate the hydrogen bonds intramolecularly. Thus, the interactivity of the proteins becomes impaired: their binding affinity decreases sharply, as there is no thermodynamic benefit from removing water surrounding properly desolvated hydrogen bonds. The results, specialized for p21ras and p53, reveal how oncogenic mutations determine a hindrance to GAP-induced hydrolysis (p21) and decrease binding affinity for DNA (p53). Furthermore, the oncogenic potential of mutations in residues not directly engaged in the interface electrostatics is assessed. The results suggest that a high sensitivity of structural defects to genetic accident might be a necessary condition to establish the existence of a proto-oncogene, an angle that merits a systematic study.     

Probing molecular docking in a charged model binding site

A model binding site was used to investigate charge-charge interactions in molecular docking. This simple site, a small (180A(3)) engineered cavity in cyctochrome c peroxidase (CCP), is negatively charged and completely buried from solvent, allowing us to explore the balance between electrostatic energy and ligand desolvation energy in a system where many of the common approximations in docking do not apply. A database with about 5300 molecules was docked into this cavity. Retrospective testing with known ligands and decoys showed that overall the balance between electrostatic interaction and desolvation energy was captured. More interesting were prospective docking scre"ens that looked for novel ligands, especially those that might reveal problems with the docking and energy methods. Based on screens of the 5300 compound database, both high-scoring and low-scoring molecules were acquired and tested for binding. Out of 16 new, high-scoring compounds tested, 15 were observed to bind. All of these were small heterocyclic cations. Binding constants were measured for a few of these, they ranged between 20microM and 60microM. Crystal structures were determined for ten of these ligands in complex with the protein. The observed ligand geometry corresponded closely to that predicted by docking. Several low-scoring alkyl amino cations were also tested and found to bind. The low docking score of these molecules owed to the relatively high charge density of the charged amino group and the corresponding high desolvation penalty. When the complex structures of those ligands were determined, a bound water molecule was observed interacting with the amino group and a backbone carbonyl group of the cavity. This water molecule mitigates the desolvation penalty and improves the interaction energy relative to that of the "naked" site used in the docking screen. Finally, six low-scoring neutral molecules were also tested, with a view to looking for false negative predictions. Whereas most of these did not bind, two did (phenol and 3-fluorocatechol). Crystal structures for these two ligands in complex with the cavity site suggest reasons for their binding. That these neutral molecules do, in fact bind, contradicts previous results in this site and, along with the alkyl amines, provides instructive false negatives that help identify weaknesses in our scoring functions. Several improvements of these are considered.     

Non-covalent binding of azo compound to peptide chain: interactions of biebrich scarlet and naphthochrome green with four model proteins

We studied the non-specific interactions of two azo compounds: biebrich scarlet (BS) and naphthochrome green (NG), with four model proteins: bovine serum albumin, ovalbumin, poly-l-lysine and hemoglobin by UV-VIS spectrometry, fluorophotometry and circular dichroism melting technique. The optimal acidities of NG and BS for binding to proteins correspond to the physiological pHs of skin and gastro tissues. The saturation binding numbers of BS and NG on peptide chains were determined and the effects of electrolytes and temperature were investigated. These interactions were fitted by the Temkin absorption model and their thermodynamic parameters were calculated. The different bindings of BS and NG to proteins were compared from their molecular structures. We inferred that an ion-pair electrostatic interaction first fixes azo compounds to basic amino acid residues and subsequent binding involves the collective action of other non-covalent bonds: hydrogen bond, van der Waals force, and hydrophobic interaction. This combination of bonds caused a change of secondary conformation of protein from β-sheet to helix and the possible process was illustrated. The potential protein toxicity resulting from such a non-specific binding was analyzed. Besides, the interaction of BS with peptide chains was applied to protein assay.