Dissecting cooperative and additive binding energetics in the affinity maturation pathway of a protein-protein interface

When two proteins associate they form a molecular interface that is a structural and energetic mosaic. Within such interfaces, individual amino acid residues contribute distinct binding energies to the complex. In combination, these energies are not necessarily additive, and significant positive or negative cooperative effects often exist. The basis of reliable algorithms to predict the specificities and energies of protein-protein interactions depends critically on a quantitative understanding of this cooperativity. We have used a model protein-protein system defined by an affinity maturation pathway, comprising variants of a T cell receptor Vbeta domain that exhibit an overall affinity range of approximately 1500-fold for binding to the superantigen staphylococcal enterotoxin C3, in order to dissect the cooperative and additive energetic contributions of residues within an interface. This molecular interaction has been well characterized previously both structurally, by x-ray crystallographic analysis, and energetically, by scanning alanine mutagenesis. Through analysis of group and individual maturation and reversion mutations using surface plasmon resonance spectroscopy, we have identified energetically important interfacial residues, determined their cooperative and additive energetic properties, and elucidated the kinetic and thermodynamic bases for molecular evolution in this system. The summation of the binding free energy changes associated with the individual mutations that define this affinity maturation pathway is greater than that of the fully matured variant, even though the affinity gap between the end point variants is relatively large. Two mutations in particular, both located in the complementarity determining region 2 loop of the Vbeta domain, exhibit negative cooperativity.     

Chemical basis of glycine riboswitch cooperativity

The glycine binding riboswitch forms a unique tandem aptamer structure that binds glycine cooperatively. We employed nucleotide analog interference mapping (NAIM) and mutagenesis to explore the chemical basis of glycine riboswitch cooperativity. Based on the interference pattern, at least two sites appear to facilitate cooperative tertiary interactions, namely, the minor groove of the P1 helix from aptamer 1 and the major groove of the P3a helix from both aptamers. Mutation of these residues altered both the cooperativity and binding affinity of the riboswitch. The data support a model in which the P1 helix of the first aptamer participates in a tertiary interaction important for cooperativity, while nucleotides in the P1 helix of the second aptamer interface with the expression platform. These data have direct analogy to well-characterized mutations in hemoglobin, which provides a framework for considering cooperativity in this RNA-based system.     

Identification of stabilizing point mutations through mutagenesis of destabilized protein libraries

Although there have been recent transformative advances in the area of protein structure prediction, prediction of point mutations that improve protein stability remains challenging. It is possible to construct and screen large mutant libraries for improved activity or ligand binding. However, reliable screens for mutants that improve protein stability do not yet exist, especially for proteins that are well folded and relatively stable. Here, we demonstrate that incorporation of a single, specific, destabilizing mutation termed parent inactivating mutation into each member of a single-site saturation mutagenesis library, followed by screening for suppressors, allows for robust and accurate identification of stabilizing mutations. We carried out fluorescence-activated cell sorting of such a yeast surface display, saturation suppressor library of the bacterial toxin CcdB, followed by deep sequencing of sorted populations. We found that multiple stabilizing mutations could be identified after a single round of sorting. In addition, multiple libraries with different parent inactivating mutations could be pooled and simultaneously screened to further enhance the accuracy of identification of stabilizing mutations. Finally, we show that individual stabilizing mutations could be combined to result in a multi-mutant that demonstrated an increase in thermal melting temperature of about 20 °C, and that displayed enhanced tolerance to high temperature exposure. We conclude that as this method is robust and employs small library sizes, it can be readily extended to other display and screening formats to rapidly isolate stabilized protein mutants.     

Dual regulatory functions of the thin filament revealed by replacement of the troponin I inhibitory peptide with a linker

Striated muscles are relaxed under low Ca(2+) concentration conditions due to actions of the thin filament protein troponin. To investigate this regulatory mechanism, an 11-residue segment of cardiac troponin I previously termed the inhibitory peptide region was studied by mutagenesis. Several mutant troponin complexes were characterized in which specific effects of the inhibitory peptide region were abrogated by replacements of 4-10 residues with Gly-Ala linkers. The mutations greatly impaired two of troponin's actions under low Ca(2+) concentration conditions: inhibition of myosin subfragment 1 (S1)-thin filament MgATPase activity and cooperative suppression of myosin S1-ADP binding to thin filaments with low myosin saturation. Inhibitory peptide replacement diminished but did not abolish the Ca(2+) dependence of the ATPase rate; ATPase rates were at least 2-fold greater when Ca(2+) rather than EGTA was present. This residual regulation was highly cooperative as a function of Ca(2+) concentration, similar to the degree of cooperativity observed with WT troponin present. Other effects of the mutations included 2-fold or less increases in the apparent affinity of the thin filament regulatory Ca(2+) sites, similar decreases in the affinity of troponin for actin-tropomyosin regardless of Ca(2+), and increases in myosin S1-thin filament ATPase rates in the presence of saturating Ca(2+). The overall results indicate that cooperative myosin binding to Ca(2+)-free thin filaments depends upon the inhibitory peptide region but that a cooperatively activating effect of Ca(2+) binding does not. The findings suggest that these two processes are separable and involve different conformational changes in the thin filament.     

Catalytic activation of human glucokinase by substrate binding: residue contacts involved in the binding of D-glucose to the super-open form and conformational transitions

alpha-D-Glucose activates glucokinase (EC 2.7.1.1) on its binding to the active site by inducing a global hysteretic conformational change. Using intrinsic tryptophan fluorescence as a probe on the alpha-D-glucose induced conformational changes in the pancreatic isoform 1 of human glucokinase, key residues involved in the process were identified by site-directed mutagenesis. Single-site W-->F mutations enabled the assignment of the fluorescence enhancement (DeltaF/F(0)) mainly to W99 and W167 in flexible loop structures, but the biphasic time course of DeltaF/F(0) is variably influenced by all tryptophan residues. The human glucokinase-alpha-D-glucose association (K(d) = 4.8 +/- 0.1 mm at 25 degrees C) is driven by a favourable entropy change (DeltaS = 150 +/- 10 J.mol(-1).K(-1)). Although X-ray crystallographic studies have revealed the alpha-d-glucose binding residues in the closed state, the contact residues that make essential contributions to its binding to the super-open conformation remain unidentified. In the present study, we combined functional mutagenesis with structural dynamic analyses to identify residue contacts involved in the initial binding of alpha-d-glucose and conformational transitions. The mutations N204A, D205A or E256A/K in the L-domain resulted in enzyme forms that did not bind alpha-D-glucose at 200 mm and were essentially catalytically inactive. Our data support a molecular dynamic model in which a concerted binding of alpha-D-glucose to N204, N231 and E256 in the super-open conformation induces local torsional stresses at N204/D205 propagating towards a closed conformation, involving structural changes in the highly flexible interdomain connecting region II (R192-N204), helix 5 (V181-R191), helix 6 (D205-Y215) and the C-terminal helix 17 (R447-K460).     

The role of protein dynamics in increasing binding affinity for an engineered protein-protein interaction established by H/D exchange mass spectrometry

It is generally accepted that protein and solvation dynamics play fundamental roles in the mechanisms of protein-protein binding; however, assessing their contribution meaningfully has not been straightforward. Here, hydrogen/deuterium exchange mass spectrometry (H/D-Ex) was employed to assess the role of dynamics for a high-affinity human growth hormone variant (hGHv) and the wild-type growth hormone (wt-hGH) each binding to the extracellular domain of their receptor (hGHbp). Comparative analysis of the transient fluctuations in the bound and unbound states revealed that helix-1 of hGHv undergoes significant transient unfolding in its unbound state, a characteristic that was not found in wt-hGH or apparent in the temperature factor data from the X-ray analysis of the unbound hGHv structure. In addition, upon hormone binding, an overall increase in stability was observed for the beta-sheet structure of hGHbp which included sites distant from the binding interface. On the basis of the stability, binding kinetics, and thermodynamic data presented, the increase in the binding free energy of hGHv is primarily generated by factors that appear to increase the energy of the unbound state relative to the free energy of the bound complex. This implies that an alternate route to engineer new interactions aiming to increase protein-protein association energies may be achieved by introducing certain mutations that destabilize one of the interacting molecules without destabilizing the resulting bound complex. Importantly, although the hGHv molecule is less stable than its wt-hGH counterpart, its resulting active ternary complex with two copies of hGHbp has comparable stability to the wt complex.     

Cooperativity during the formation of peptide/MHC class II complexes

To generate an effective immune response, class II major histocompatibility complex molecules (MHCII) must present a diverse array of peptide ligands for recognition by T lymphocytes. Peptide/MHCII complexes are stabilized by hydrophobic anchoring of peptide side chains to pockets in the MHCII protein and the formation of hydrogen bonds to the peptide backbone. Many current models of peptide/MHCII association assume an additive and independent contribution of the interactions between major MHCII pockets and corresponding side chains in the peptide. However, significant conformational rearrangements occur in both the peptide and MHCII during binding. Therefore, we hypothesize that peptide binding to MHCII could be viewed as a folding process in which both molecules cooperate to produce the final conformation. To directly test this hypothesis, we adapt a serial mutagenesis strategy to study cooperativity in the interaction of the human MHCII HLA-DR1 and a peptide derived from influenza hemagglutinin. Substitutions in either the peptide or HLA-DR1 that are predicted to interfere with hydrogen bond formation show cooperative effects on complex stability and affinity. Substitution of a peptide side chain that provides a hydrophobic contact also contributes to the cooperative effect, suggesting a role for all energetic sources in the folding process. We propose that cooperativity throughout the peptide-binding groove reflects the folding of segments of the MHCII molecule into helices around the peptide with a concomitant folding of the peptide into a polyproline helix. The implications of cooperativity for peptide/MHCII structure and epitope selection are discussed.     

The role of gamma-carboxyglutamyl residues in the positive cooperative binding of Ca2+ to blood coagulation factor X

1. The calcium binding properties of factor X and its analogous decarboxyprotein have been compared with the aid of flow rate dialysis and ultraviolet difference spectroscopy. 2. Factor X binds approx. 20 mol of calcium per mol of protein. The first four sites exhibit positive cooperativity. 3. Changes in the ultraviolet difference spectrum when Ca2+ is bound suggest a conformational change. 4. In decarboxyfactor X low affinity of Ca2+ and no ligand-induced conformational change was observed. It is concluded that the presence of gamma-carboxyglutamate residues is a prerequisite for positive cooperative Ca2+ binding.     

Identification of novel peptide antagonists for GPIIb/IIIa from a conformationally constrained phage peptide library.

Methods have recently been developed to present vast libraries of random peptides on the surface of filamentous phage. To introduce a degree of conformational constraint into random peptides, a library of hexapeptides flanked by cysteine residues (capable of forming cyclic disulfides) was constructed. This library was screened using the platelet glycoprotein, IIb/IIIa, which mediates the aggregation of platelets through binding of fibrinogen. A variety of peptides containing the sequence Arg-Gly-Asp or Lys-Gly-Asp were discovered and synthesized. The cyclic, disulfide-bonded forms of the peptides bound IIb/IIIa with dissociation constants in the nanomolar range, while reduced forms or an analogue in which Ser replaced the Cys residues bound considerably less tightly. These results demonstrate the feasibility for introducing conformational constraints into random peptide libraries and also demonstrates the potential for using phage peptide libraries to discover pharmacologically active lead compounds.     

The tryptophan residues of aspartate transcarbamylase: site-directed mutagenesis and time-resolved fluorescence spectroscopy.

Aspartate transcarbamylase (EC 2.1.3.2) contains two tryptophan residues in position 209 and 284 of the catalytic chains (c) and no such chromophore in the regulatory chains (r). Thus, as a dodecamer [(c3)2(r2)3] the native enzyme molecule contains 12 tryptophan residues. The present study of the regulatory conformational changes in this enzyme is based on the fluorescence properties of these intrinsic probes. Site-directed mutagenesis was used in order to differentiate the respective contributions of the two tryptophans to the fluorescence properties of the enzyme and to identify the mobility of their environment in the course of the different regulatory processes. Each of these tryptophan residues gives two independent fluorescence decays, suggesting that the catalytic subunit exists in two slightly different conformational states. The binding of the substrate analog N-phosphonacetyl-L-aspartate promotes the same fluorescence signal whether or not the catalytic subunits are associated with the regulatory subunits, suggesting that the substrate-induced conformational change of the catalytic subunit is the essential trigger for the quaternary structure transition involved in cooperativity. The binding of the substrate analog affects mostly the environment of tryptophan 284, while the binding of the activator ATP affects mostly the environment of tryptophan 209, confirming that this activator acts through a mechanism different from that involved in homotropic cooperativity.     

Interaction between the peripheral site residues of human butyrylcholinesterase, D70 and Y332, in binding and hydrolysis of substrates

Human butyrylcholinesterase displays substrate activation with positively charged butyrylthiocholine (BTC) as the substrate. Peripheral anionic site (PAS) residues D70 and Y332 appear to be involved in the initial binding of charged substrates and in activation control. To determine the contribution of PAS residues to binding and hydrolysis of quaternary substrates and activation control, the single mutants D70G/Y and Y332F/A/D and the double mutants Y332A/D70G and Y332D/D70Y were studied. Steady-state hydrolysis of the charged substrates, BTC and succinyldithiocholine, and the neutral ester o-nitrophenyl butyrate was measured. In addition, inhibition of wild-type and mutant enzymes by tetramethylammonium was investigated, at low concentrations of BTC. Single and double mutants of D70 and Y332 showed little or no substrate activation, suggesting that both residues were important for activation control. The effects of double mutations on D70 and Y332 were complex. Double-mutant cycle analysis provided evidence for interaction between these residues. The category of interaction (either synergistic, additive, partially additive or antagonistic) was found to depend on the nature of the substrate and on measured binding or kinetic parameters. This complexity reflects both the cross-talk between residues involved in the sequential formation of productive Michaelian complexes and the effect of peripheral site residues on catalysis. It is concluded that double mutations on the PAS induce a conformational change in the active site gorge of butyrylcholinesterase that can alter both substrate binding and enzyme acylation.     

Elimination of the C-cap in ubiquitin - structure, dynamics and thermodynamic consequences

Single amino acid substitutions rarely produce substantial changes in protein structure. Here we show that substitution of the C-cap residue in the alpha-helix of ubiquitin with proline (34P variant) leads to dramatic structural changes. The resulting conformational perturbation extends over the last two turns of the alpha-helix and leads to enhanced flexibility for residues 27-37. Thermodynamic analysis of this ubiquitin variant using differential scanning calorimetry reveals that the thermal unfolding transition remains highly cooperative, exhibiting two-state behavior. Similarities with the wild type in the thermodynamic parameters (heat capacity change upon unfolding and m-value) of unfolding monitored by DSC and chemical denaturation suggests that the 34P variant has comparable buried surface area. The hydrophobic core of 34P variant is not packed as well as that of the wild type protein as manifested by a lower enthalpy of unfolding. The increased mobility of the polypeptide chain of this ubiquitin variant allows the transient opening of the hydrophobic core as evidenced by ANS binding. Taken together, these results suggest exceptional robustness of cooperativity in protein structures.     

Enthalpy-entropy compensation and cooperativity as thermodynamic epiphenomena of structural flexibility in ligand-receptor interactions

Ligand binding is a thermodynamically cooperative process in many biochemical systems characterized by the conformational flexibility of the reactants. However, the contribution of conformational entropy to cooperativity of ligation needs to be elucidated. Here, we perform kinetic and thermodynamic analyses on a panel of cycle-mutated peptides, derived from influenza H3 HA(306-319), interacting with wild type and a mutant HLA-DR. We observe that, within a certain range of peptide affinity, this system shows isothermal entropy-enthalpy compensation (iEEC). The incremental increases in conformational entropy measured as disruptive mutations are added in the ligand or receptor are more than sufficient in magnitude to account for the experimentally observed lack of free-energy decrease cooperativity. Beyond this affinity range, compensation is not observed, and therefore, the ability of the residual interactions to form a stable complex decreases in an exponential fashion. Taken together, our results indicate that cooperativity and iEEC constitute the thermodynamic epiphenomena of the structural fluctuation that accompanies ligand-receptor complex formation in flexible systems. Therefore, ligand binding affinity prediction needs to consider how each source of binding energy contributes synergistically to the folding and kinetic stability of the complex in a process based on the trade-off between structural tightening and restraint of conformational mobility.     

Small angle X-ray studies reveal that Aspergillus niger glucoamylase has a defined extended conformation and can form dimers in solution

The industrially important glucoamylase 1 is an exo-acting glycosidase with substrate preference for alpha-1,4 and alpha-1,6 linkages at non-reducing ends of starch. It consists of a starch binding and a catalytic domain interspersed by a highly glycosylated polypeptide linker. The linker function is poorly understood and structurally undescribed, and data regarding domain organization and intramolecular functional cooperativity are conflicting or non-comprehensive. Here, we report a combined small angle x-ray scattering and calorimetry study of Aspergillus niger glucoamylase 1, glucoamylase 2, which lacks a starch binding domain, and an engineered low-glycosylated variant of glucoamylase 1 with a short linker. Low resolution solution structures show that the linker adopts a compact structure rendering a well defined extended overall conformation to glucoamylase. We demonstrate that binding of a short heterobidentate inhibitor simultaneously directed toward the catalytic and starch binding domains causes dimerization of glucoamylase and not, as suggested previously, an intramolecular conformational rearrangement mediated by linker flexibility. Our results suggest that glucoamylase functions via transient dimer formation during hydrolysis of insoluble substrates and address the question of the cooperative effect of starch binding and hydrolysis.     

Improving the stability of an antibody variable fragment by a combination of knowledge-based approaches: validation and mechanisms

Numerous approaches have been described to obtain variable fragments of antibodies (Fv or scFv) that are sufficiently stable for their applications. Here, we combined several knowledge-based methods to increase the stability of pre-existing scFvs by design. Firstly, the consensus sequence approach was used in a non-stringent way to predict a large basic set of potentially stabilizing mutations. These mutations were then prioritized by other methods of design, mainly the formation of additional hydrogen bonds, an increase in the hydrophilicity of solvent exposed residues, and previously described mutations in other antibodies. We validated this combined method with antibody mAbD1.3, directed against lysozyme. Fourteen potentially stabilizing mutations were designed and introduced into scFvD1.3 by site-directed mutagenesis, either individually or in combinations. We characterized the effects of the mutations on the thermodynamic stability of scFvD1.3 by experiments of unfolding with urea, monitored by spectrofluorometry, and tested the additivity of their effects by double-mutant cycles. We also quantified the individual contributions of the resistance to denaturation ([urea](1/2)) and cooperativity of unfolding (m) to the variations of stability and the energy of coupling between mutations by a novel approach. Most mutations (75%) were stabilizing and none was destabilizing. The progressive recombination of the mutations into the same molecule of scFvD1.3 showed that their effects were mostly additive or synergistic, provided a large overall increase in protein stability (9.1 kcal/mol), and resulted in a highly stable scFvD1.3 derivative. The mechanisms of the mutations and of their combinations involved variations in the resistance to denaturation, cooperativity of unfolding, and likely residual structures of the denatured state, which was constrained by two disulfide bonds. This combined method should be applicable to any recombinant antibody fragment, through a single step of mutagenesis.     

In vitro evolution of a hyperstable Gbeta1 variant

An in-vitro selection strategy was used to obtain strongly stabilized variants of the beta1 domain of protein G (Gbeta1). In a two-step approach, first candidate positions with a high potential for stabilization were identified in Gbeta1 libraries that were created by error-prone PCR, and then, after randomization of these positions by saturation mutagenesis, strongly stabilized variants were selected. For both steps the in-vitro selection method Proside was employed. Proside links the stability of a protein with the infectivity of a filamentous phage. Ultimately, residues from the two best selected variants were combined in a single Gbeta1 molecule. This variant with the four mutations E15V, T16L, T18I, and N37L showed an increase of 35.1 degrees C in the transition midpoint and of 28.5 kJ mol(-1) (at 70 degrees C) in the Gibbs free energy of stabilization. It was considerably more stable than the best variant from a previous Proside selection, in which positions were randomized that had originally been identified by computational design. Only a single substitution (T18I) was found in both selections. The best variants from the present selection showed a higher cooperativity of thermal unfolding, as indicated by an increase in the enthalpy of unfolding by about 60 kJ mol(-1). This increase is apparently correlated with the presence of Leu residues that were selected at the positions 16 and 37.     

Structure‐based design of combinatorial mutagenesis libraries

The development of protein variants with improved properties (thermostability, binding affinity, catalytic activity, etc.) has greatly benefited from the application of high‐throughput screens evaluating large, diverse combinatorial libraries. At the same time, since only a very limited portion of sequence space can be experimentally constructed and tested, an attractive possibility is to use computational protein design to focus libraries on a productive portion of the space. We present a general‐purpose method, called “Structure‐based Optimization of Combinatorial Mutagenesis ” (SOCoM ), which can optimize arbitrarily large combinatorial mutagenesis libraries directly based on structural energies of their constituents. SOCoM chooses both positions and substitutions, employing a combinatorial optimization framework based on library‐averaged energy potentials in order to avoid explicitly modeling every variant in every possible library. In case study applications to green fluorescent protein, β‐lactamase, and lipase A, SOCoM optimizes relatively small, focused libraries whose variants achieve energies comparable to or better than previous library design efforts, as well as larger libraries (previously not designable by structure‐based methods) whose variants cover greater diversity while still maintaining substantially better energies than would be achieved by representative random library approaches. By allowing the creation of large‐scale combinatorial libraries based on structural calculations, SOCoM promises to increase the scope of applicability of computational protein design and improve the hit rate of discovering beneficial variants. While designs presented here focus on variant stability (predicted by total energy), SOCoM can readily incorporate other structure‐based assessments, such as the energy gap between alternative conformational or bound states.     

Activation of phenylalanine hydroxylase: effect of substitutions at Arg68 and Cys237

Phenylalanine hydroxylase (PAH) is a multidomain tetrameric enzyme that displays positive cooperative substrate binding. This cooperative response is believed to be of physiological significance as a mechanism that controls L-Phe homeostasis in blood. The substrate induces an activating conformational change in the enzyme affecting the secondary, tertiary, and quaternary structures. Chemical modification and substitution with a negatively charged residue of Cys237 in human PAH (hPAH) also result in activation of the enzyme. As seen in the modeled structure of full-length hPAH, Cys237 is located in the catalytic domain close to residues in the oligomerization and regulatory domains of an adjacent subunit in the dimer, notably to Arg68. This residue is located in a prominent loop (68-75), which also has contacts with the dimerization motif from the same subunit. To investigate further the involvement of Cys237 and Arg68 in the activation of the enzyme, we have prepared mutants of hPAH at these positions, with substitutions of different charge and size. The mutations C237D, R68A, and C237A cause an increase of the basal activity and affinity for L-Phe, while the mutation C237R results in reduced affinity for the substrate and elimination of the positive cooperativity. The conformational changes induced by the mutations were studied by far-UV circular dichroism, fluorescence spectroscopy, and molecular dynamics simulations. All together, our results indicate that the activating mutations induce a series of conformational changes including both the displacement of the inhibitory N-terminal sequence (residues 19-33) that covers the active site and the domain movements around the hinge region Arg111-Thr117, in addition to the rearrangement of the loop 68-75. The same conformational changes appear to be involved in the activation of PAH induced by L-Phe.