Protein-protein binding is often associated with changes in protonation state

pK(a) values of ionizable residues have been calculated using the PROPKA method and structures of 75 protein-protein complexes and their corresponding free forms. These pK(a) values were used to compute changes in protonation state of individual residues, net changes in protonation state of the complex relative to the uncomplexed proteins, and the correction to a binding energy calculated assuming standard protonation states at pH 7. For each complex, two different structures for the uncomplexed form of the proteins were used: the X-ray structures determined for the proteins in the absence of the other protein and the individual protein structures taken from the structure of the complex (referred to as unbound and bound structures, respectively). In 28 and 77% of the cases considered here, protein-protein binding is accompanied by a complete (>95%) or significant (>50%) change in protonation state of at least one residue using unbound structures. Furthermore, in 36 and 61% of the cases, protein-protein binding is accompanied by a complete or significant net change in protonation state of the complex relative to the separated monomers. Using bound structures, the corresponding values are 12, 51, 20, and 48%. Comparison to experimental data suggest that using unbound and bound structures lead to over- and underestimation of binding-induced protonation state changes, respectively. Thus, we conclude that protein-protein binding is often associated with changes in protonation state of amino acid residues and with changes in the net protonation state of the proteins. The pH-dependent correction to the binding energy contributes at least one order of magnitude to the binding constant in 45 and 23%, using unbound and bound structures, respectively.     

A fast and accurate computational approach to protein ionization

We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.     

Electrostatic effects in unfolded staphylococcal nuclease

Structure-based calculations of pKa values and electrostatic free energies of proteins assume that electrostatic effects in the unfolded state are negligible. In light of experimental evidence showing that this assumption is invalid for many proteins, and with increasing awareness that the unfolded state is more structured and compact than previously thought, a detailed examination of electrostatic effects in unfolded proteins is warranted. Here we address this issue with structure-based calculations of electrostatic interactions in unfolded staphylococcal nuclease. The approach involves the generation of ensembles of structures representing the unfolded state, and calculation of Coulomb energies to Boltzmann weight the unfolded state ensembles. Four different structural models of the unfolded state were tested. Experimental proton binding data measured with a variant of nuclease that is unfolded under native conditions were used to establish the validity of the calculations. These calculations suggest that weak Coulomb interactions are an unavoidable property of unfolded proteins. At neutral pH, the interactions are too weak to organize the unfolded state; however, at extreme pH values, where the protein has a significant net charge, the combined action of a large number of weak repulsive interactions can lead to the expansion of the unfolded state. The calculated pKa values of ionizable groups in the unfolded state are similar but not identical to the values in small peptides in water. These studies suggest that the accuracy of structure-based calculations of electrostatic contributions to stability cannot be improved unless electrostatic effects in the unfolded state are calculated explicitly.     

Protonation states and pH titration in the photocycle of photoactive yellow protein

Photoactive yellow protein (PYP) undergoes a light-driven cycle of color and protonation states that is part of a mechanism of bacterial phototaxis. This article concerns functionally important protonation states of PYP and the interactions that stabilize them, and changes in the protonation state during the photocycle. In particular, the chromophore pK(a) is known to be shifted down so that the chromophore is negatively charged in the ground state (dark state) even though it is buried in the protein, while nearby Glu46 has an unusually high pK(a). The photocycle involves changes of one or both of these protonation states. Calculations of pK(a) values and protonation states using a semi-macroscopic electrostatic model are presented for the wild-type and three mutants, in both the ground state and the bleached (I(2)) intermediate state. Calculations allowing multiple H-bonding arrangements around the chromophore also have been carried out. In addition, ground-state pK(a) values of the chromophore have been measured by UV-visible spectroscopy for the wild-type and the same three mutants. Because of the unusual protonation states and strong electrostatic interactions, PYP represents a severe test of the ability of theoretical models to yield correct calculations of electrostatic interactions in proteins. Good agreement between experiment and theory can be obtained for the ground state provided the protein interior is assumed to have a relatively low dielectric constant, but only partial agreement between theory and experiment is obtained for the bleached state. We also present a reinterpretation of previously published data on the pH-dependence of the recovery of the ground state from the bleached state. The new analysis implies a pK(a) value of 6.37 for Glu46 in the bleached state, which is consistent with other available experimental data, including data that only became available after this analysis. The new analysis suggests that signal transduction is modulated by the titration properties of the bleached state, which are in turn determined by electrostatic interactions. Overall, the results of this study provide a quantitative picture of the interactions responsible for the unusual protonation states of the chromophore and Glu46, and of protonation changes upon bleaching.     

Experiments suggest that simulations may overestimate electrostatic contributions to the mechanical stability of a fibronectin type III domain

Steered molecular dynamics simulations have previously been used to investigate the mechanical properties of the extracellular matrix protein fibronectin. The simulations suggest that the mechanical stability of the tenth type III domain from fibronectin (FNfn10) is largely determined by a number of critical hydrogen bonds in the peripheral strands. Interestingly, the simulations predict that lowering the pH from 7 to approximately 4.7 will increase the mechanical stability of FNfn10 significantly (by approximately 33 %) due to the protonation of a few key acidic residues in the A and B strands. To test this simulation prediction, we used single-molecule atomic force microscopy (AFM) to investigate the mechanical stability of FNfn10 at neutral pH and at lower pH where these key residues have been shown to be protonated. Our AFM experimental results show no difference in the mechanical stability of FNfn10 at these different pH values. These results suggest that some simulations may overestimate the role played by electrostatic interactions in determining the mechanical stability of proteins.     

Solid phase proteomics: dramatic reinforcement of very weak protein-protein interactions

Very weak protein-protein interactions may play a critical role in cell physiology but they are not easily detectable in "in vitro" experiments. To detect these weak interactions, we have developed a strategy that included: (a) design of a rapid and very effective crosslinking of protein-protein complexes with poly-functional reagents; (b) selective adsorption of very large proteins on lowly activated ionic exchangers, based on the need of a multipoint physical adsorption to incorporate the proteins into the matrix; (c) purification by selective adsorption of protein-protein complexes formed by strong protein-protein interactions, via selective adsorption of the complexes on lowly activated ionic exchangers via multi-protein physical adsorption and leaving the non-associated proteins in the solution; (d) reinforcement of very weak protein-protein interactions by selective adsorption of the complex on lowly activated ionic exchange supports via a synergetic cooperation of the weak protein-protein interaction plus the interactions of both proteins with the support enabling the almost full shifting of the equilibrium towards the association position; (e) control of the aggregation state of proteins like BSA, formed by weak protein-protein interactions. In this last case, it seems that the interaction of the protein molecules placed on the borders of the aggregate with the groups on the support partially stabilizes the whole aggregate, although, some molecules of the aggregate cannot interact with the support. The size of the aggregates may be defined by controlling the concentration of ionised groups on the support: the less activated the supports are, the bigger the complexes. In this way, solid-phase proteomics could be a very interesting tool to detect weak protein-protein interactions.     

Weak protein interactions and pH- and temperature-dependent aggregation of human Fc1

The Fc (fragment crystallizable) is a common structural region in immunoglobulin gamma (IgG) proteins, IgG-based multi-specific platforms, and Fc-fusion platform technologies. Changes in conformational stability, protein-protein interactions, and aggregation of NS0-produced human Fc1 were quantified experimentally as a function of pH (4 to 6) and temperature (30 to 77°C), using a combination of differential scanning calorimetry, laser light scattering, size-exclusion chromatography, and capillary electrophoresis. The Fc1 was O-glycosylated at position 3 (threonine), and confirmed to correspond to the intact IgG1 by comparison with Fc1 produced by cleavage of the parent IgG1. Changing the pH caused large effects for thermal unfolding transitions, but it caused surprisingly smaller effects for electrostatic protein-protein interactions. The aggregation behavior was qualitatively similar across different solution conditions, with soluble dimers and larger oligomers formed in most cases. Aggregation rates spanned approximately 5 orders of magnitude and could be divided into 2 regimes: (i) Arrhenius, unfolding-limited aggregation at temperatures near or above the midpoint-unfolding temperature of the C_H2 domain; (ii) a non-Arrhenius regime at lower temperatures, presumably as a result of the temperature dependence of the unfolding enthalpy for the C_H2 domain. The non-Arrhenius regime was most pronounced for lower temperatures. Together with the weak protein-protein repulsions, these highlight challenges that are expected for maintaining long-term stability of biotechnology products that are based on human Fc constructs.     

Weak protein interactions and pH- and temperature-dependent aggregation of human Fc1

The Fc (fragment crystallizable) is a common structural region in immunoglobulin gamma (IgG) proteins, IgG-based multi-specific platforms, and Fc-fusion platform technologies. Changes in conformational stability, protein-protein interactions, and aggregation of NS0-produced human Fc1 were quantified experimentally as a function of pH (4 to 6) and temperature (30 to 77°C), using a combination of differential scanning calorimetry, laser light scattering, size-exclusion chromatography, and capillary electrophoresis. The Fc1 was O-glycosylated at position 3 (threonine), and confirmed to correspond to the intact IgG1 by comparison with Fc1 produced by cleavage of the parent IgG1. Changing the pH caused large effects for thermal unfolding transitions, but it caused surprisingly smaller effects for electrostatic protein-protein interactions. The aggregation behavior was qualitatively similar across different solution conditions, with soluble dimers and larger oligomers formed in most cases. Aggregation rates spanned approximately 5 orders of magnitude and could be divided into 2 regimes: (i) Arrhenius, unfolding-limited aggregation at temperatures near or above the midpoint-unfolding temperature of the C_H2 domain; (ii) a non-Arrhenius regime at lower temperatures, presumably as a result of the temperature dependence of the unfolding enthalpy for the C_H2 domain. The non-Arrhenius regime was most pronounced for lower temperatures. Together with the weak protein-protein repulsions, these highlight challenges that are expected for maintaining long-term stability of biotechnology products that are based on human Fc constructs.     

Prediction of charge-induced molecular alignment: residual dipolar couplings at pH 3 and alignment in surfactant liquid crystalline phases

Recently we reported that the alignment tensor of a biological macromolecule, which was dissolved in a dilute suspension of highly negatively charged filamentous phage at close to neutral pH, can be predicted from the molecule’s 3D charge distribution and shape (Zweckstetter et al. 2004). Here it is demonstrated that this approach is also applicable to alignment of proteins in liquid crystalline phases formed by filamentous phage at low pH. Residual dipolar couplings (RDCs) predicted by our simple electrostatic model for the B1 domain of protein G in fd phage at pH 3 fit very well with the experimental values. The sign of charge–shape predicted one-bond ¹H–¹⁵N dipolar couplings for the B1 domain of protein G (GB1) was inverted at pH 3 compared to neutral pH, in agreement with experimental observations. Our predictions indicate that this is a feature specific for GB1. In addition, it is shown that RDCs induced in the protein ubiquitin by the presence of a positively charged surfactant system comprising cetylpyridinium bromide/hexanol/sodium bromide can be predicted accurately by a simple electrostatic alignment model. This shows that steric and electrostatic interactions dominate weak alignment of biomolecules for a wide range of pH values both in filamentous phage and in surfactant liquid crystalline phases.     

Coarse-grained strategy for modeling protein stability in concentrated solutions

We present a coarse-grained approach for modeling the thermodynamic stability of single-domain globular proteins in concentrated aqueous solutions. Our treatment derives effective protein-protein interactions from basic structural and energetic characteristics of the native and denatured states. These characteristics, along with the intrinsic (i.e., infinite dilution) thermodynamics of folding, are calculated from elementary sequence information using a heteropolymer collapse theory. We integrate this information into Reactive Canonical Monte Carlo simulations to investigate the connections between protein sequence hydrophobicity, protein-protein interactions, protein concentration, and the thermodynamic stability of the native state. The model predicts that sequence hydrophobicity can affect how protein concentration impacts native-state stability in solution. In particular, low hydrophobicity proteins are primarily stabilized by increases in protein concentration, whereas high hydrophobicity proteins exhibit richer nonmonotonic behavior. These trends appear qualitatively consistent with the available experimental data. Although factors such as pH, salt concentration, and protein charge are also important for protein stability, our analysis suggests that some of the nontrivial experimental trends may be driven by a competition between destabilizing hydrophobic protein-protein attractions and entropic crowding effects.     

Calbindin D(9k): a protein optimized for calcium binding at neutral pH.

The binding of calcium ions by EF-hand proteins depends strongly on the electrostatic interactions between Ca(2+) ions and negatively charged residues of these proteins. We have investigated the pH dependence of the binding of Ca(2+) ions by calbindin D(9k). This protein offers a unique possibility for interpretation of such data since the pK(a) values of all ionizable groups are known. The binding is independent of pH between 7 and 9, where maximum calcium affinity is observed. An abrupt decrease in the binding affinity is observed at pH values below 7. This decrease is due to protonation of acidic groups, leading to modification of protein charges. The pH dependence of the product of the two macroscopic Ca(2+)-binding constants can be formally described by the involvement of two acidic groups with pK(a) = 6.6. Monte Carlo calculations show that the reduction of Ca(2+) binding is strictly determined by variable electrostatic interactions due to pH-dependent changes not only in the binding sites, but also of the overall charge of the protein.     

The pH-induced release of iron from transferrin investigated with a continuum electrostatic model.

A reduction in pH induces the release of iron from transferrin in a process that involves a conformational change in the protein from a closed to an open form. Experimental evidence suggests that there must be changes in the protonation states of certain, as yet not clearly identified, residues in the protein accompanying this conformational change. Such changes in protonation states of residues and the consequent changes in electrostatic interactions are assumed to play a large part in the mechanism of release of iron from transferrin. Using the x-ray crystal structures of human ferri- and apo-lactoferrin, we calculated the pKa values of the titratable residues in both the closed (iron-loaded) and open (iron-free) conformations with a continuum electrostatic model. With the knowledge of a residue's pKa value, its most probable protonation state at any specified pH may be determined. The preliminary results presented here are in good agreement with the experimental observation that the binding of ferric iron and the synergistic anion bicarbonate/carbonate results in the release of approximately three H+ ions. It is suggested that the release of these three H+ ions may be accounted for, in most part, by the deprotonation of the bicarbonate and residues Tyr-92, Lys-243, Lys-282, and Lys-285 together with the protonation of residues Asp-217 and Lys-277.     

Screening of protein-ligand interactions by affinity chromatography

This paper examines affinity chromatography (AC) as an alternative tool for the determination of protein-ligand interactions for the particular case in which the ligand is the same protein. The methodology is less labor-intensive and more sample-efficient than traditional methods used to measure the second virial coefficient (B(22)), a parameter commonly used to evaluate protein-protein interactions. The chromatographic capacity factor (k') was studied for lysozyme and equine serum albumin for a wide range of experimental solution conditions such as crystallizing agent concentration, protein concentration and pH. Parallel experiments using AC to determine k' and static light scattering (SLS) to determine B(22) showed that the two parameters were highly correlated. Two different column volumes ( approximately 1 and approximately 0.1 mL) were tested and gave essentially the same values for k', showing the feasibility of miniaturization.     

Electrostatic contributions to residue-specific protonation equilibria and proton binding capacitance for a small protein

Charge-charge interactions in proteins are important in a host of biological processes. Here we use 13C NMR chemical shift data for individual aspartate and glutamate side chain carboxylate groups to accurately detect site-specific protonation equilibria in a variant of the B1 domain of protein G (PGB1-QDD). Carbon chemical shifts are dominated by changes in the electron distribution within the side chain and therefore excellent reporters of the charge state of individual groups, and the data are of high precision. We demonstrate that it is possible to detect local charge interactions within this small protein domain that stretch and skew the chemical shift titration curves away from "ideal" behavior and introduce a framework for the analysis of such convoluted data to study local charge-charge interactions and electrostatic coupling. It is found that, due to changes in electrostatic potential, the proton binding affinity, Ka, of each carboxyl group changes throughout the titration process and results in a linearly pH dependent pKa value. This result could be readily explained by calculations of direct charge-charge interactions based on Coulomb's law. In addition, the slope of pKa versus pH was dependent on screening by salt, and this dependence allowed the selective study of charge-charge interactions. For PGB1-QDD, it was established that mainly differences in self-energy, and not direct charge-charge interactions, are responsible for shifted pKa values within the protein environment.     

Energetics of protein homodimerization: effects of water sequestering on the formation of beta-lactoglobulin dimer

Transient protein-protein interactions are functionally relevant as a control mechanism in a variety of biological processes. Analysis of the 3D structure of protein-protein complexes indicates that water molecules trapped at the interface are very common; however, their role in the stability and specificity of protein homodimer interactions has been not addressed yet. To provide new insights into the energetic bases that govern the formation of highly hydrated interfaces, the dissociation process of bovine beta lg variant A at a neutral pH was characterized here thermodynamically by conducting dilution experiments with an isothermal titration calorimeter. Association was enthalpically driven throughout the temperature range spanned. DeltaH and deltaC(p) were significantly more negative than estimates based on surface area changes, suggesting the occurrence of effects additional to the dehydration of the contact surfaces between subunits. Near-UV CD spectra proved to be independent of protein concentration, indicating a rigid body-like association. Furthermore, the process proved not to be coupled to significant changes in the protonation state of ionizable groups or counterion exchange. In contrast, both osmotic stress experiments and a computational analysis of the dimer's 3D structure indicated that a large number of water molecules are incorporated into the interface upon association. Numerical estimates considering the contributions of interface area desolvation and water immobilization accounted satisfactorily for the experimental deltaC(p). Thus, our study highlights the importance of explicitly considering the effects of water sequestering to perform a proper quantitative analysis of the formation of homodimers with highly hydrated interfaces.     

Automated computational framework for the analysis of electrostatic similarities of proteins

Charge plays an important role in protein-protein interactions. In the case of excessively charged proteins, their electrostatic potentials contribute to the processes of recognition and binding with other proteins or ligands. We present an automated computational framework for determining the contribution of each charged amino acid to the electrostatic properties of proteins, at atomic resolution level. This framework involves computational alanine scans, calculation of Poisson-Boltzmann electrostatic potentials, calculation of electrostatic similarity distances (ESDs), hierarchical clustering analysis of ESDs, calculation of solvation free energies of association, and visualization of the spatial distributions of electrostatic potentials. The framework is useful to classify families of mutants with similar electrostatic properties and to compare them with the parent proteins in the complex. The alanine scan mutants introduce perturbations in the local electrostatic properties of the proteins and aim in delineating the contribution of each mutated amino acid in the spatial distribution of electrostatic potential, and in biological function when electrostatics is a dominant contributing factor in protein-protein interactions. The framework can be used to design new proteins with tailored electrostatic properties, such as immune system regulators, inhibitors, and vaccines, and in guiding experimental studies. We present an example for the interaction of the immune system protein C3d (the d-fragment of complement protein C3) with its receptor CR2, and we discuss our data in view of a binding site controversy.     

Positive aspects of negative design: simultaneous selection of specificity and interaction stability

The energetic determinants that drive specific protein-protein interactions are not entirely understood. We describe simultaneous in vivo selection of specific and stable interactions using homologous peptides which compete with protein libraries for an interaction with a target molecule. Library members binding to their target, and promoting cell growth, must outcompete competitor interactions with the target (i.e., competition) and evade binding to the competitors (i.e., negative design). We term this a competitive and negative design initiative (CANDI). We combined CANDI with a protein-fragment complementation assay (PCA) and observed major specificity improvements, by driving selection of winning library members that bind their target with maximum efficacy, ensuring that otherwise energetically accessible alternatives are inaccessible. CANDI-PCA has been used with libraries targeted at coiled coil regions of oncogenic AP-1 components cJun and cFos. We demonstrate that comparable hydrophobic and electrostatic contributions in desired species are compromised in nondesired species when CANDI is executed, demonstrating that both core and electrostatic residues are required to direct specific interactions. Major energetic differences (>or=5.6 kcal/mol) are observed between desired and nondesired interaction stabilities for a CANDI-PCA derived peptide relative to a conventional PCA derived helix, with significantly more stability (3.2 kcal/mol) than the wild-type cJun-cFos complex. As a negative control, a library lacking a residue repertoire able to generate a specific and stable helix was tested. Negative protein design has broad implications in generating specific and therapeutically relevant peptide-based drugs, proteins able to act with minimal cross-talk to homologues or analogues, and in nanobiotechnological design.     

Empirical parametrization of pK values for carboxylic acids in proteins using a genetic algorithm

Considerable effort has been devoted to the development of theoretical electrostatic methods to predict the pK values of ionizable residues in proteins. However, predictions appear often to be still at the qualitative or semi-quantitative level. We believe that, with the increasing number experimentally available pK values for proteins of known structure, an alternative approach becomes feasible: the empirical parametrization of the experimental protein pK database. Of course, in the long term, this empirical approach is no substitute for rigorous electrostatic analysis but, in the short term, it may prove to have useful predictive power and it may help to pinpoint the main structural determinants of pK values in proteins. Here we demonstrate the feasibility of the parametrization approach by fitting (using a genetic algorithm as fitting tool) the database for carboxylic acid pK values in proteins on the basis of an empirical equation that takes into account the two following kinds of effects: (1) long-range charge-charge interactions; (2) interactions of the given carboxylic acid group with its environment in the protein, which are described in terms of contributions from the different kind of atoms present in the protein (atomic contributions).