Electrostatic enhancement of diffusion-controlled protein-protein association: comparison of theory and experiment on barnase and barstar

Protein flexibility, not disorder, is intrinsic to molecular recognition

An 'intrinsically disordered protein' (IDP) is assumed to be unfolded in the cell and perform its biological function in that state. We contend that most intrinsically disordered proteins are in fact proteins waiting for a partner (PWPs), parts of a multi-component complex that do not fold correctly in the absence of other components. Flexibility, not disorder, is an intrinsic property of proteins, exemplified by X-ray structures of many enzymes and protein-protein complexes. Disorder is often observed with purified proteins in vitro and sometimes also in crystals, where it is difficult to distinguish from flexibility. In the crowded environment of the cell, disorder is not compatible with the known mechanisms of protein-protein recognition, and, foremost, with its specificity. The self-assembly of multi-component complexes may, nevertheless, involve the specific recognition of nascent polypeptide chains that are incompletely folded, but then disorder is transient, and it must remain under the control of molecular chaperones and of the quality control apparatus that obviates the toxic effects it can have on the cell.

Kinetics of ligand-receptor interaction reveals an induced-fit mode of binding in a cyclic nucleotide-activated protein

Many receptors and ion channels are activated by ligands. One key question concerns the binding mechanism. Does the ligand induce conformational changes in the protein via the induced-fit mechanism? Or does the protein preexist as an ensemble of conformers and the ligand selects the most complementary one, via the conformational selection mechanism? Here, we study ligand binding of a tetrameric cyclic nucleotide-gated channel from Mesorhizobium loti and of its monomeric binding domain (CNBD) using rapid mixing, mutagenesis, and structure-based computational biology. Association rate constants of ∼10(7) M(-1) s(-1) are compatible with diffusion-limited binding. Ligand binding to the full-length CNG channel and the isolated CNBD differ, revealing allosteric control of the CNBD by the effector domain. Finally, mutagenesis of allosteric residues affects only the dissociation rate constant, suggesting that binding follows the induced-fit mechanism. This study illustrates the strength of combining mutational, kinetic, and computational approaches to unravel important mechanistic features of ligand binding.

Use of Molecular Dynamics for the Refinement of an Electrostatic Model for the In Silico Design of a Polymer Antidote for the Anticoagulant Fondaparinux

Molecular dynamics (MD) simulations results are herein incorporated into an electrostatic model used to determine the structure of an effective polymer-based antidote to the anticoagulant fondaparinux. In silico data for the polymer or its cationic binding groups has not, up to now, been available, and experimental data on the structure of the polymer-fondaparinux complex is extremely limited. Consequently, the task of optimizing the polymer structure is a daunting challenge. MD simulations provided a means to gain microscopic information on the interactions of the binding groups and fondaparinux that would have otherwise been inaccessible. This was used to refine the electrostatic model and improve the quantitative model predictions of binding affinity. Once refined, the model provided guidelines to improve electrostatic forces between candidate polymers and fondaparinux in order to increase association rate constants.

Impact of protein/protein interactions on global intermolecular translocation rates of the transcription factors Sox2 and Oct1 between DNA cognate sites analyzed by z-exchange NMR spectroscopy

Oct1 and Sox2 synergistically regulate developmental genes by binding to adjacent sites within promoters. We have investigated the kinetics of global intermolecular translocation of Sox2 and Oct1 between cognate sites located on different DNA molecules by z-exchange NMR spectroscopy. In the Hoxb1 promoter, the Sox2 and Oct1 sites are immediately adjacent to one another, and the intermolecular translocation rates are too slow to be measured by z-exchange spectroscopy. By introducing a 3-bp insertion between the Sox2 and Oct1 sites to mimic the spacing in the FGF4 enhancer, the interprotein contact surface is reduced, and the translocation rates are increased. Interaction between Sox2 and the POU-specific domain (POU(S)) of Oct1 does not affect the translocation mechanism but modulates the rates. Translocation involves only jumping (dissociation and reassociation) for Sox2, but both jumping and direct intersegment transfer (no dissociation into free solution) for Oct1. The dissociation (k(off) ∼1.5 s(-1)) and association (k(on) ∼5.1 × 10(9) m(-1)s(-1)) rate constants for Sox2 are reduced 4-fold and increased 5-fold, respectively, in the presence of Oct1. k(off) (∼3.5 s(-1)) for Oct1 is unaffected by Sox2, whereas k(on) (∼1.3 × 10(9) m(-1)s(-1)) is increased ∼13-fold. The direct intermolecular translocation rate (k(inter) ∼1.8 × 10(4) m(-1)s(-1)) for the POU(S) domain of Oct1 is reduced 2-fold by Sox2, whereas that for the POU homeodomain (POU(HD)) of Oct1 (k(inter) ∼ 1.7 × 10(4) m(-1)s(-1)) remains unaltered, consistent with the absence of contacts between Sox2 and POU(HD). The data suggest a model for the sequence of binding events involved in synergistic gene regulation by Sox2 and Oct1.

Cooperation between catalytic and DNA binding domains enhances thermostability and supports DNA synthesis at higher temperatures by thermostable DNA polymerases

We have previously introduced a general kinetic approach for comparative study of processivity, thermostability, and resistance to inhibitors of DNA polymerases [Pavlov, A. R., et al. (2002) Proc. Natl. Acad. Sci. U.S.A.99, 13510-13515]. The proposed method was successfully applied to characterize hybrid DNA polymerases created by fusing catalytic DNA polymerase domains with various sequence-nonspecific DNA binding domains. Here we use the developed kinetic analysis to assess basic parameters of DNA elongation by DNA polymerases and to further study the interdomain interactions in both previously constructed and new chimeric DNA polymerases. We show that connecting helix-hairpin-helix (HhH) domains to catalytic polymerase domains can increase thermostability, not only of DNA polymerases from extremely thermophilic species but also of the enzyme from a faculatative thermophilic bacterium Bacillus stearothermophilus. We also demonstrate that addition of Topo V HhH domains extends efficient DNA synthesis by chimerical polymerases up to 105 °C by maintaining processivity of DNA synthesis at high temperatures. We found that reversible high-temperature structural transitions in DNA polymerases decrease the rates of binding of these enzymes to the templates. Furthermore, activation energies and pre-exponential factors of the Arrhenius equation suggest that the mechanism of electrostatic enhancement of diffusion-controlled association plays a minor role in binding of templates to DNA polymerases.

Brownian dynamics study of the association between the 70S ribosome and elongation factor G

Protein synthesis on the ribosome involves a number of external protein factors that bind at its functional sites. One key factor is the elongation factor G (EF-G) that facilitates the translocation of transfer RNAs between their binding sites, as well as advancement of the messenger RNA by one codon. The details of the EF-G/ribosome diffusional encounter and EF-G association pathway still remain unanswered. Here, we applied Brownian dynamics methodology to study bimolecular association in the bacterial EF-G/70S ribosome system. We estimated the EF-G association rate constants at 150 and 300 mM monovalent ionic strengths and obtained reasonable agreement with kinetic experiments. We have also elucidated the details of EF-G/ribosome association paths and found that positioning of the L11 protein of the large ribosomal subunit is likely crucial for EF-G entry to its binding site.

Rapid search for specific sites on DNA through conformational switch of nonspecifically bound proteins

We develop a theory for the rapid search of specific sites on DNA, via a mechanism in which a nonspecifically-bound protein can switch between two conformations. In the "inactive" conformation, the bound protein has favorable, nonspecific interactions with the DNA, but cannot be recognized by the target site. In the "active" conformation, the protein is recognized by the target site but has a very rugged energy surface elsewhere on the DNA. The rate constant for protein binding to the specific site is calculated by an approach in which the protein, after reaching the DNA surface via 3D diffusion, searches for the target site via 1D diffusion while being allowed to escape to the bulk solution. Mindful of the pitfalls of many previous approximate treatments, we validate our approach against a rigorous solution of the problem when the protein has a fixed conformation. In the 1D diffusion toward the specific site, a conformationally switchable protein predominantly adopts the inactive conformation due to the favorable nonspecific interactions with the DNA, thus maximizing the 1D diffusion constant and minimizing the chance of escape to the bulk solution. Once at the target site, a transition to the active conformation allows the protein to be captured. This induced-switch mechanism provides robust speedup of protein-DNA binding rates, and appears to be adopted by many transcription factors and DNA-modifying enzymes.

Direct visualization reveals dynamics of a transient intermediate during protein assembly

Interactions between proteins underlie numerous biological functions. Theoretical work suggests that protein interactions initiate with formation of transient intermediates that subsequently relax to specific, stable complexes. However, the nature and roles of these transient intermediates have remained elusive. Here, we characterized the global structure, dynamics, and stability of a transient, on-pathway intermediate during complex assembly between the Signal Recognition Particle (SRP) and its receptor. We show that this intermediate has overlapping but distinct interaction interfaces from that of the final complex, and it is stabilized by long-range electrostatic interactions. A wide distribution of conformations is explored by the intermediate; this distribution becomes more restricted in the final complex and is further regulated by the cargo of SRP. These results suggest a funnel-shaped energy landscape for protein interactions, and they provide a framework for understanding the role of transient intermediates in protein assembly and biological regulation.

Cortisone-reductase deficiency associated with heterozygous mutations in 11beta-hydroxysteroid dehydrogenase type 1

In peripheral target tissues, levels of active glucocorticoid hormones are controlled by 11β-hydroxysteroid dehydrogenase type 1 (11β-HSD1), a dimeric enzyme that catalyzes the reduction of cortisone to cortisol within the endoplasmic reticulum. Loss of this activity results in a disorder termed cortisone reductase deficiency (CRD), typified by increased cortisol clearance and androgen excess. To date, only mutations in H6PD, which encodes an enzyme supplying cofactor for the reaction, have been identified as the cause of disease. Here we examined the HSD11B1 gene in two cases presenting with biochemical features indicative of a milder form of CRD in whom the H6PD gene was normal. Novel heterozygous mutations (R137C or K187N) were found in the coding sequence of HSD11B1. The R137C mutation disrupts salt bridges at the subunit interface of the 11β-HSD1 dimer, whereas K187N affects a key active site residue. On expression of the mutants in bacterial and mammalian cells, activity was either abolished (K187N) or greatly reduced (R137C). Expression of either mutant in a bacterial system greatly reduced the yield of soluble protein, suggesting that both mutations interfere with subunit folding or dimer assembly. Simultaneous expression of mutant and WT 11β-HSD1 in bacterial or mammalian cells, to simulate the heterozygous condition, indicated a marked suppressive effect of the mutants on both the yield and activity of 11β-HSD1 dimers. Thus, these heterozygous mutations in the HSD11B1 gene have a dominant negative effect on the formation of functional dimers and explain the genetic cause of CRD in these patients.

A biosensor study indicating that entropy, electrostatics, and receptor glycosylation drive the binding interaction between interleukin-7 and its receptor

The interaction between interleukin-7 (IL-7) and its α-receptor, IL-7Rα, plays fundamental roles in the development, survival, and homeostasis of B- and T-cells. N-Linked glycosylation of human IL-7Rα enhances its binding affinity for human IL-7 300-fold versus that of the nonglycosylated receptor through an allosteric mechanism. The N-glycans of IL-7Rα do not participate directly in the binding interface with IL-7. This biophysical study involves dissection of the properties of binding of IL-7 to both nonglycosylated and glycosylated forms of the IL-7Rα extracellular domain (ECD) as functions of salt, pH, and temperature using surface plasmon resonance (SPR) spectroscopy. Interactions of IL-7 with both IL-7Rα variants display weaker binding affinities with increasing salt concentrations primarily reflected by changes in the first on rates of a two-step reaction pathway. The electrostatic parameter of the IL-7-IL-7Rα interaction is not driven by complementary charge interactions through residues at the binding interface or N-glycan composition of IL-7Rα, but presumably by favorable global charges of the two proteins. van't Hoff analysis indicates both IL-7-IL-7Rα interactions are driven by large favorable entropy changes and smaller unfavorable (nonglycosylated complex) and favorable (glycosylated complex) enthalpy changes. Eyring analysis of the IL-7-IL-7Rα interactions reveals different reaction pathways and barriers for the transition-state thermodynamics with the enthalpy and entropy changes of IL-7 binding to nonglycosylated and glycosylated IL-7Rα. There were no discernible heat capacity changes for the equilibrium or transition-state binding thermodynamics of the IL-7-IL-7Rα interactions. The results suggest that the unbound nonglycosylated IL-7Rα samples an extensive conformational landscape relative to the unbound glycosylated IL-7Rα, potentially explaining the switch from a "conformationally controlled" reaction (k(1) ∼ 10(2) M(-1) s(-1)) for the nonglycosylated interaction to a "diffusion-controlled" reaction (k(1) ∼ 10(6) M(-1) s(-1)) for the glycosylated interaction. Thus, a large favorable entropy change, a global favorable electrostatic component, and glycosylation of the receptor, albeit not at the interface, contribute significantly to the interaction between IL-7 and the IL-7Rα ECD.

Structural model of the pre-pore ring-like structure of Panton-Valentine leukocidin: providing dimensionality to biophysical and mutational data

Panton-Valentine leukocidin (PVL) is a bipartite toxin that plays an important role in the pathogenesis of methicillin-resistant Staphylococcus aureus. Recent clinical data suggest a correlation between PVL and severe cases of S. aureus pneumonia. A clear understanding of the structure and function of PVL is critical to the development of novel, effective treatments. Here, we report an all-atom model of the macromolecular structure of Panton-Valentine leukocidin in its octameric, pre-pore conformation that confirms and extends our understanding of the toxin's mechanism of action.

Effect of the ordered interfacial water layer in protein complex formation: A nonlocal electrostatic approach

Using a nonlocal electrostatic approach that incorporates the short-range structure of the contacting media, we evaluated the electrostatic contribution to the energy of the complex formation of two model proteins. In this study, we have demonstrated that the existence of an ordered interfacial water layer at the protein-solvent interface reduces the charging energy of the proteins in the aqueous solvent, and consequently increases the electrostatic contribution to the protein binding (change in free energy upon the complex formation of two proteins). This is in contrast with the finding of the continuum electrostatic model, which suggests that electrostatic interactions are not strong enough to compensate for the unfavorable desolvation effects.

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.

The impact of viral RNA on the association rates of capsid protein assembly: bacteriophage MS2 as a case study

A large number of single-stranded RNA viruses, which form a major class of all viruses, co-assemble their protein container and their genomic material. The multiple roles of the viral genome in this process are presently only partly understood. Recent experimental results indicate that RNA, in addition to its function as a repository for genetic information, could play important functional roles during the assembly of the viral protein containers. An investigation of the impact of genomic RNA on the association of the protein subunits may therefore provide further insights into the mechanism of virus assembly. We study here the impact of viral RNA on the association rates of the capsid proteins during virus assembly. As a case study, we consider the viral capsid of bacteriophage MS2, which is formed from 60 asymmetric (AB) and 30 symmetric (CC) protein dimers. Using Brownian dynamics simulations, we investigate the effect of the binding of an RNA stem-loop (the translational repressor) on the association rates of the capsid protein dimers. Our analysis shows that translational repressor binding results in self-association of AB dimers being inhibited, whilst association of AB with CC dimers is greatly enhanced. This provides an explanation for experimental results in which an alternating assembly pattern of AB and CC dimer addition to the growing assembly intermediate has been observed to be the dominant mode of assembly. The presence of the RNA hence dramatically decreases the number of dominant assembly pathways and thereby reduces the complexity of the self-assembly process of these viruses.

Assembly dynamics and stability of the pneumococcal epsilon zeta antitoxin toxin (PezAT) system from Streptococcus pneumoniae

The pneumococcal epsilon zeta antitoxin toxin (PezAT) system is a chromosomally encoded, class II toxin antitoxin system from the human pathogen Streptococcus pneumnoniae. Neutralization of the bacteriotoxic protein PezT is carried out by complex formation with its cognate antitoxin PezA. Here we study the stability of the inhibitory complex in vivo and in vitro. We found that toxin release is impeded in Escherichia coli and Bacillus subtilis due to the proteolytic resistance of PezA once bound to PezT. These findings are supported by in vitro experiments demonstrating a strong thermodynamic stabilization of both proteins upon binding. A detailed kinetic analysis of PezAT assembly revealed that these particular features of PezAT are based on a strong, electrostatically guided binding mechanism leading to a stable toxin antitoxin complex with femtomolar affinity. Our data show that PezAT complex formation is distinct to all other conventional toxin antitoxin modules and a controlled mode of toxin release is required for activation.

Rate theories for biologists

Some of the rate theories that are most useful for modeling biological processes are reviewed. By delving into some of the details and subtleties in the development of the theories, the review will hopefully help the reader gain a more than superficial perspective. Examples are presented to illustrate how rate theories can be used to generate insight at the microscopic level into biomolecular behaviors. An attempt is made to clear up a number of misconceptions in the literature regarding popular rate theories, including the appearance of Planck's constant in the transition-state theory and the Smoluchowski result as an upper limit for protein-protein and protein-DNA association rate constants. Future work in combining the implementation of rate theories through computer simulations with experimental probes of rate processes, and in modeling effects of intracellular environments so that theories can be used for generating rate constants for systems biology studies is particularly exciting.

Steric and electrostatic complementarity in the assembly of two-dimensional virus arrays

A highly ordered assembly of biological molecules provides a powerful means to study the organizational principles of objects at the nanoscale. Two-dimensional cowpea mosaic virus arrays were assembled in an ordered manner on mica using osmotic depletion effects and a drop-and-dry method. The packing of the virus array was controlled systematically from rhombic packing to hexagonal packing by modulating the concentrations of poly(ethylene glycol) surfactant in the virus solutions. The orientation and packing symmetry of the virus arrays were found to be tuned by the concentrations of surfactants in the sample solutions. A phenomenological model for the present system is proposed to explain the assembly array morphology under the influence of the surfactant. Steric and electrostatic complementarity of neighboring virus capsids is found to be the key factors in controlling the symmetry of packing.

Mechanistic details of a protein-protein association pathway revealed by paramagnetic relaxation enhancement titration measurements

Protein-protein association generally proceeds via the intermediary of a transient, lowly populated, encounter complex ensemble. The mechanism whereby the interacting molecules in this ensemble locate their final stereospecific structure is poorly understood. Further, a fundamental question is whether the encounter complex ensemble is an effectively homogeneous population of nonspecific complexes or whether it comprises a set of distinct structural and thermodynamic states. Here we use intermolecular paramagnetic relaxation enhancement (PRE), a technique that is exquisitely sensitive to lowly populated states in the fast exchange regime, to characterize the mechanistic details of the transient encounter complex interactions between the N-terminal domain of Enzyme I (EIN) and the histidine-containing phosphocarrier protein (HPr), two major bacterial signaling proteins. Experiments were conducted at an ionic strength of 150 mM NaCl to eliminate any spurious nonspecific associations not relevant under physiological conditions. By monitoring the dependence of the intermolecular transverse PRE (Gamma(2)) rates measured on (15)N-labeled EIN on the concentration of paramagnetically labeled HPr, two distinct types of encounter complex configurations along the association pathway are identified and dissected. The first class, which is in equilibrium with and sterically occluded by the specific complex, probably involves rigid body rotations and small translations near or at the active site. In contrast, the second class of encounter complex configurations can coexist with the specific complex to form a ternary complex ensemble, which may help EIN compete with other HPr binding partners in vivo by increasing the effective local concentration of HPr even when the active site of EIN is occupied.

pH-dependent association of proteins. The test case of monoclonal antibody HyHEL-5 and its antigen hen egg white lysozyme

We describe a method for determining diffusion-controlled rate constants for protein-protein association that explicitly includes the solution pH. The method combines the transient-complex theory for computing electrostatically enhanced association rates with an approach based on a rigorous thermodynamic cycle and partition functions for energy levels characterizing protonation states of associating proteins and their complexes. To test our method, we determine the pH-dependent kinetics of association of the HyHEL-5 antibody with its antigen hen egg white lysozyme. It was shown experimentally that their association rate constant depends on pH, increasing linearly in the pH range 6-8 and saturating or even exhibiting a flat maximum in the pH range 8-10. The presented methodology leads to a qualitative agreement with the experimental data. Our approach allows one to study diffusion-controlled protein-protein association under different pH conditions by taking into account the ensembles of protonation states rather than just the most probable protonation state of each protein.

The effect of different force applications on the protein-protein complex Barnase-Barstar

Steered molecular dynamics simulations are a tool to examine the energy landscape of protein-protein complexes by applying external forces. Here, we analyze the influence of the velocity and geometry of the probing forces on a protein complex using this tool. With steered molecular dynamics, we probe the stability of the protein-protein complex Barnase-Barstar. The individual proteins are mechanically labile. The Barnase-Barstar binding site is more stable than the folds of the individual proteins. By using different force protocols, we observe a variety of responses of the system to the applied tension.

Improved binding of raf to Ras.GDP is correlated with biological activity

The GTP-binding protein Ras plays a central role in the regulation of various cellular processes, acting as a molecular switch that triggers signaling cascades. Only Ras bound to GTP is able to interact strongly with effector proteins like Raf kinase, phosphatidylinositol 3-kinase, and RalGDS, whereas in the GDP-bound state, the stability of the complex is strongly decreased, and signaling is interrupted. To determine whether this process is only controlled by the stability of the complex, we used computer-aided protein design to improve the interaction between Ras and effector. We challenged the Ras.Raf complex in this study because Raf among all effectors shows the highest Ras affinity and the fastest association kinetics. The proposed mutations were characterized as to their changes in dynamics and binding strength. We demonstrate that Ras-Raf interaction can only be improved at the cost of a loss in specificity of Ras.GTP versus Ras.GDP. As shown by NMR spectroscopy, the Raf mutation A85K leads to a shift of Ras switch I in the GTP-bound as well as in the GDP-bound state, thereby increasing the complex stability. In a luciferase-based reporter gene assay, Raf A85K is associated with higher signaling activity, which appears to be a mere matter of Ras-Raf affinity.

Knitting and untying the protein network: modulation of protein ensembles as a therapeutic strategy

Proteins constitute the working machinery and structural support of all organisms. In performing a given function, they must adopt highly specific structures that can change with their level of activity, often through the direct or indirect action of other proteins. Indeed, proteins typically function within an ensemble, rather than individually. Hence, they must be sufficiently flexible to interact with each other and execute diverse tasks. The discovery that errors within these groups can ultimately cause disease has led to a paradigm shift in drug discovery, from an emphasis on single protein targets to a holistic approach whereby entire ensembles are targeted.

Fruitful and futile encounters along the association reaction between proteins

The association reaction between pairs of proteins proceeds through an encounter complex that develops into the final complex. Here, we combined Brownian dynamics simulations with experimental studies to analyze the structures of the encounter complexes along the association reaction between TEM1-beta-lactamase and its inhibitor, beta-lactamase-inhibitor protein. The encounter complex can be considered as an ensemble of short-lived low free-energy states that are stabilized primarily by electrostatic forces and desolvation. For the wild-type, the simulation showed two main encounter regions located outside the physical binding site. One of these regions was located near the experimentally determined transition state. To validate whether these encounters are fruitful or futile, we examined three groups of mutations that altered the encounter. The first group consisted of mutations that increased the experimental rate of association through electrostatic optimization. This resulted in an increase in the size of the encounter region located near the experimentally determined transition state, as well as a decrease in the energy of this region and an increase in the number of successful trajectories (i.e., encounters that develop into complex). A second group of mutations was specifically designed to either increase or decrease the size and energy of the second encounter complex, but either way it did not affect k(on). A third group of mutations consisted of residues that increased k(on) without significantly affecting the encounter complexes. These results indicate that the size and energy of the encounter regions are only two of several parameters that lead to fruitful association, and that electrostatic optimization is a major driving force in fast association.

Binding mechanism of an SH3 domain studied by NMR and ITC

Complexes between Src-homology 3 domains and proline-rich target peptides can have lifetimes on the order of milliseconds, making them too short-lived for kinetic characterization by conventional methods. Nuclear magnetic resonance (NMR) dynamics experiments are ideally suited to study such rapid binding equilibria, and additionally provide information on partly bound intermediate states. We used NMR together with isothermal titration calorimetry (ITC) to characterize the interaction of the SH3 domain from the Fyn tyrosine kinase with a 12-residue peptide at temperatures between 10 and 50 degrees C. NMR data at all temperatures are consistent with an effectively two-state binding reaction, such that any intermediates are either very weakly populated or exchange extremely rapidly with the free or bound forms. Dissociation rate constants, determined by CPMG and ZZ-exchange NMR experiments, range from k(off)(10 degrees C) = 4.5 s(-1) to k(off)(50 degrees C) = 331 s(-1). ITC data at all temperatures follow a simple two-state interaction model. Binding is favored enthalpically, with a dissociation enthalpy, DeltaH(D)(30 degrees C) = 15.4 kcal mol(-1), and disfavored entropically, with a dissociation entropy, DeltaS(D)(30 degrees C) = 20.0 cal mol(-1) K(-1). The free protein and peptide have significantly higher heat capacity than the bound complex, DeltaC(p) = 352 cal mol(-1) K(-1), which is consistent with the largely hydrophobic character of the binding interface. An Eyring plot of k(off) values gives an activation enthalpy of dissociation, DeltaH(D)(double dagger)(30 degrees C) = 19.3 kcal mol(-1) and exhibits slight curvature consistent with the ITC-derived value of DeltaC(p). The curvature suggests that nonpolar residues of the hydrophobic interface are solvated in the transition state for dissociation. Association rate constants were calculated as k(on) = k(off)/K(D), and range from k(on)(10 degrees C) = 1.03 x 10(8) M(-1) s(-1) to k(on)(50 degrees C) = 2.0 x 10(8) M(-1) s(-1), with an apparent activation enthalpy, DeltaH(A)(double dagger) = 3.4 kcal mol(-1). Both the magnitudes and temperature dependence of k(on) values are consistent with a diffusion-limited association mechanism. The combination of NMR and ITC data sheds light on how the Fyn tyrosine kinase is activated by binding to proline-rich targets, and represents a powerful approach for characterizing transient protein/ligand interactions.

Dissection of the high rate constant for the binding of a ribotoxin to the ribosome

Restrictocin belongs to a family of site-specific ribonucleases that kill cells by inactivating the ribosome. The restrictocin-ribosome binding rate constant was observed to exceed 10(10) M(-1) s(-1). We have developed a transient-complex theory to model the binding rates of protein-protein and protein-RNA complexes. The theory predicts the rate constant as k(a) = k(a0) exp(-DeltaG(el)*/k(B)T), where k(a0) is the basal rate constant for reaching the transient complex, located at the outer boundary of the bound state, by random diffusion, and DeltaG(el)* is the average electrostatic interaction free energy of the transient complex. Here, we applied the transient-complex theory to dissect the high restrictocin-ribosome binding rate constant. We found that the binding rate of restrictocin to the isolated sarcin/ricin loop is electrostatically enhanced by approximately 300-fold, similar to results found in other protein-protein and protein-RNA complexes. The ribosome provides an additional 10,000-fold rate enhancement because of two synergistic mechanisms afforded by the distal regions of the ribosome. First, they provide additional electrostatic attraction with restrictocin. Second, they reposition the transient complex into a region where local electrostatic interactions of restrictocin with the sarcin/ricin loop are particularly favorable. Our calculations rationalize a host of experimental observations and identify a strategy for designing proteins that bind their targets with high speed.

Principal determinants leading to transition state formation of a protein-protein complex, orientation trumps side-chain interactions

The binding transition state (TS) is the rate-limiting step for transient molecular interactions. This important step in the molecular recognition process, however, is largely understood only at a qualitative level. To establish a more quantitative picture of the TS structure, we exploit a set of biophysical techniques that have provided major insights in protein folding applications. As a model system representing the large class of "weakly charged" protein-protein interactions, we examine the binding of a variety of human growth hormone (hGH) variants to the human growth hormone receptor (hGHR) and the human prolactin receptor (hPRLR). hGH variants were chosen to probe different features of the TS structure, based on their highly reengineered interfaces. Both Eyring and urea (m value) analyses suggest that the majority of binding surface burial occurs after TS. A comprehensive phi analysis showed that individual hGH interface residues do not contribute energetically to the stability of the TS, but there is a TS "hot spot" in the receptor. Zinc dependence studies that take advantage of an endogenous tetracoordinated interfacial metal binding demonstrate that surfaces of the molecules have attained a high orientational complementarity by the time the TS is reached. The model that best fits these data are that a "knobs-into-holes" process precisely aligns the two molecular interfaces in forming the TS structure. Surprisingly, most of the thermodynamic character of the binding reaction is focused in the fine-tuning process occurring after TS.

Coulomb energies of protein-protein complexes with monopole-free charge distributions

From an analysis of Coulomb energy distributions of a large set of protein-protein complexes we show that the positive tail in the energy distribution disappears when the monopole-monopole term, the only energy term independent of inter-subunit orientations, is removed. This indicates that unfavorable Coulomb energies associated with subunit orientations are excluded in protein-protein complexes. The overall result remained unchanged when solvent effects were included. Our results have important bearing on the restriction of subunit orientations in protein-protein complexes and complement a recent work [K. Brock, K. Talley, K. Coley, P. Kundrotas, E. Alexov, Optimization of electrostatic interactions in protein-protein complexes, Biophys. J. 93 (2007) 3340-3352.] which showed that Coulomb energy of interaction in protein-protein complexes is sequence-optimized.

Association rate constants of ras-effector interactions are evolutionarily conserved

Evolutionary conservation of protein interaction properties has been shown to be a valuable indication for functional importance. Here we use homology interface modeling of 10 Ras-effector complexes by selecting ortholog proteins from 12 organisms representing the major eukaryotic branches, except plants. We find that with increasing divergence time the sequence similarity decreases with respect to the human protein, but the affinities and association rate constants are conserved as predicted by the protein design algorithm, FoldX. In parallel we have done computer simulations on a minimal network based on Ras-effector interactions, and our results indicate that in the absence of negative feedback, changes in kinetics that result in similar binding constants have strong consequences on network behavior. This, together with the previous results, suggests an important biological role, not only for equilibrium binding constants but also for kinetics in signaling processes involving Ras-effector interactions. Our findings are important to take into consideration in system biology approaches and simulations of biological networks.

Cellular signal transductions processes are based on protein interactions. Proteins can either associate transiently with each other or form stable complexes, and the strength of the interaction is described by the affinity (the affinity is the ratio between the rate of dissociation and association). Protein complexes with similar affinities can bind and dissociate with different rates, and these rates describe the kinetic properties of protein binding. These kinetic rates are important for signaling; however, to what extent individual changes in such rate constants are biologically important or whether the affinity is more crucial might be different in different signaling processes. In this study we analyze whether association rates are conserved during evolution, because evolutionary conservation of protein biochemical properties is usually a valuable indication of its importance. We analyzed the binding of Ras proteins to effector domains, which are central proteins in many signal transduction pathways, in different organisms. On the basis of homology modeling and energy calculations we find that association rates are conserved, although the sequence similarity decreases compared to the human protein. Our finding should encourage further analysis of the importance of kinetics for cellular signal transduction.

Visualizing lowly-populated regions of the free energy landscape of macromolecular complexes by paramagnetic relaxation enhancement

Many biological macromolecular interactions proceed via lowly-populated, highly transient species that arise from rare excursions between the minimum free energy configuration and other local minima of the free energy landscape. Little is known about the structural properties of such lowly-occupied states since they are difficult to trap and hence inaccessible to conventional structural and biophysical techniques. Yet these states play a crucial role in a variety of dynamical processes including molecular recognition and binding, allostery, induced-fit and self-assembly. Here we highlight recent progress in paramagnetic nuclear magnetic resonance to detect, visualize and characterize lowly-populated transient species at equilibrium. The underlying principle involves the application of paramagnetic relaxation enhancement (PRE) in the fast exchange regime. Under these conditions the footprint of the minor species can be observed in the PRE profiles measured for the major species, providing distances between the paramagnetic label and protons of interest are shorter in the minor species than the major one. Ensemble simulated annealing refinement directly against the PRE data permits one to obtain structural data on the minor species. We have used the PRE (a) to detect and characterize the stochastic target search process whereby a sequence-specific transcription factor (the Hox-D9 homeodomain) binds to non-cognate DNA sites as a means of enhancing the rate of specific association via intramolecular sliding and intermolecular translocation; (b) to directly visualize the distribution of non-specific transient encounter complexes involved in the formation of stereospecific protein-protein complexes; (c) to detect and visualize ultra-weak self-association of a protein, a process that is relevant to early nucleation events involved in the formation of higher order structures; and (d) to determine the structure of a minor species for a multidomain protein (maltose binding protein) where large interdomain motions are associated with ligand binding, thereby shedding direct light on the fundamental question of allostery versus induced fit in this system. The PRE offers unique opportunities to directly probe and explore in structural terms lowly-populated regions of the free energy landscape and promises to yield fundamental new insights into biophysical processes.

Conformational changes and flexibility in T-cell receptor recognition of peptide-MHC complexes

A necessary feature of the immune system, TCR (T-cell receptor) cross-reactivity has been implicated in numerous autoimmune pathologies and is an underlying cause of transplant rejection. Early studies of the interactions of alphabeta TCRs (T-cell receptors) with their peptide-MHC ligands suggested that conformational plasticity in the TCR CDR (complementarity determining region) loops is a dominant contributor to T-cell cross-reactivity. Since these initial studies, the database of TCRs whose structures have been solved both bound and free is now large enough to permit general conclusions to be drawn about the extent of TCR plasticity and the types and locations of motion that occur. In the present paper, we review the conformational differences between free and bound TCRs, quantifying the structural changes that occur and discussing their possible roles in specificity and cross-reactivity. We show that, rather than undergoing major structural alterations or 'folding' upon binding, the majority of TCR CDR loops shift by relatively small amounts. The structural changes that do occur are dominated by hinge-bending motions, with loop remodelling usually occurring near loop apexes. As predicted from previous studies, the largest changes are in the hypervariable CDR3alpha and CDR3beta loops, although in some cases the germline-encoded CDR1alpha and CDR2alpha loops shift in magnitudes that approximate those of the CDR3 loops. Intriguingly, the smallest shifts are in the germline-encoded loops of the beta-chain, consistent with recent suggestions that the TCR beta domain may drive ligand recognition.

Thermodynamics of T-cell receptor-peptide/MHC interactions: progress and opportunities

alphabeta T-cell receptors (TCRs) recognize peptide antigens presented by class I or class II major histocompatibility complex molecules (pMHC). Here we review the use of thermodynamic measurements in the study of TCR-pMHC interactions, with attention to the diversity in binding thermodynamics and how this is related to the variation in TCR-pMHC interfaces. We show that there is no enthalpic or entropic signature for TCR binding; rather, enthalpy and entropy changes vary in a compensatory manner that reflects a narrow free energy window for the interactions that have been characterized. Binding enthalpy and entropy changes do not correlate with structural features such as buried surface area or the number of hydrogen bonds within TCR-pMHC interfaces, possibly reflecting the myriad of contributors to binding thermodynamics, but likely also reflecting a reliance on van't Hoff over calorimetric measurements and the unaccounted influence of equilibria linked to binding. TCR-pMHC binding heat capacity changes likewise vary considerably. In some cases, the heat capacity changes are consistent with conformational differences between bound and free receptors, but there is little data indicating these conformational differences represent the need to organize disordered CDR loops. In this regard, we discuss how thermodynamics may provide additional insight into conformational changes occurring upon TCR binding. Finally, we highlight opportunities for the further use of thermodynamic measurements in the study of TCR-pMHC interactions, not only for understanding TCR binding in general, but also for understanding specifics of individual interactions and the engineering of TCRs with desired molecular recognition properties.

Replica exchange simulations of transient encounter complexes in protein-protein association

Recent paramagnetic relaxation enhancement (PRE) studies on several weakly interacting protein complexes have unequivocally demonstrated the existence of transient encounter complexes. Here, we present a computational method to study protein-protein binding by creating equilibrium ensembles that include both specific and nonspecific protein complexes. In a joint analysis of simulation and experiment we explore the physical nature and underlying physicochemical characteristics of encounter complexes involving three protein-protein interactions of the bacterial phosphotransferase system. Replica exchange Monte Carlo simulations using a coarse-grained energy function recover the structures of the specific complexes and produce binding affinities in good agreement with experiment. Together with the specific complex, a relatively small number of distinct nonspecific complexes largely accounts for the measured PRE data. The combined relative population of the latter is less than approximately 10%. The binding interfaces of the specific and nonspecific complexes differ primarily in size but exhibit similar amino acid compositions. We find that the overall funnel-shaped energy landscape of complex formation is dominated by the specific complex, a small number of structured nonspecific complexes, and a diffuse cloud of loosely bound complexes connecting the specific and nonspecific binding sites with each other and the unbound state. Nonspecific complexes may not only accelerate the binding kinetics by enhancing the rate of success of random diffusional encounters but also play a role in protein function as alternative binding modes.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

Crystal structural analysis of protein-protein interactions drastically destabilized by a single mutation

The complex of barnase (bn) and barstar (bs), which has been widely studied as a model for quantitative analysis of protein-protein interactions, is significantly destabilized by a single mutation, namely, bs Asp39 --> Ala, which corresponds to a change of 7.7 kcal x mol(-1) in the free energy of binding. However, there has been no structural information available to explain such a drastic destabilization. In the present study, we determined the structure of the mutant complex at 1.58 A resolution by X-ray crystallography. The complex was similar to the wild-type complex in terms of overall and interface structures; however, the hydrogen bond network mediated by water molecules at the interface was significantly different. Several water molecules filled the cavity created by the mutation and consequently caused rearrangement of the hydrated water molecules at the interface. The water molecules were redistributed into a channel-like structure that penetrated into the complex. Furthermore, molecular dynamics simulations showed that the mutation increased the mobility of water molecules at the interface. Since such a drastic change in hydration was not observed in other mutant complexes of bn and bs, the significant destabilization of the interaction may be due to this channel-like structure of hydrated water molecules.

Binding-induced folding of a natively unstructured transcription factor

Transcription factors are central components of the intracellular regulatory networks that control gene expression. An increasingly recognized phenomenon among human transcription factors is the formation of structure upon target binding. Here, we study the folding and binding of the pKID domain of CREB to the KIX domain of the co-activator CBP. Our simulations of a topology-based Gō-type model predict a coupled folding and binding mechanism, and the existence of partially bound intermediates. From transition-path and Phi-value analyses, we find that the binding transition state resembles the unstructured state in solution, implying that CREB becomes structured only after committing to binding. A change of structure following binding is reminiscent of an induced-fit mechanism and contrasts with models in which binding occurs to pre-structured conformations that exist in the unbound state at equilibrium. Interestingly, increasing the amount of structure in the unbound pKID reduces the rate of binding, suggesting a "fly-casting"-like process. We find that the inclusion of attractive non-native interactions results in the formation of non-specific encounter complexes that enhance the on-rate of binding, but do not significantly change the binding mechanism. Our study helps explain how being unstructured can confer an advantage in protein target recognition. The simulations are in general agreement with the results of a recently reported nuclear magnetic resonance study, and aid in the interpretation of the experimental binding kinetics.

Electrostatic rate enhancement and transient complex of protein-protein association

The association of two proteins is bounded by the rate at which they, via diffusion, find each other while in appropriate relative orientations. Orientational constraints restrict this rate to approximately 10(5)-10(6) M(-1) s(-1). Proteins with higher association rates generally have complementary electrostatic surfaces; proteins with lower association rates generally are slowed down by conformational changes upon complex formation. Previous studies (Zhou, Biophys J 1997;73:2441-2445) have shown that electrostatic enhancement of the diffusion-limited association rate can be accurately modeled by $k_{\bf D}$ = $k_{D}0\ {exp} ( - \langle U_{el} \rangle;{\star}/k_{B} T),$ where k(D) and k(D0) are the rates in the presence and absence of electrostatic interactions, respectively, U(el) is the average electrostatic interaction energy in a "transient-complex" ensemble, and k(B)T is the thermal energy. The transient-complex ensemble separates the bound state from the unbound state. Predictions of the transient-complex theory on four protein complexes were found to agree well with the experiment when the electrostatic interaction energy was calculated with the linearized Poisson-Boltzmann (PB) equation (Alsallaq and Zhou, Structure 2007;15:215-224). Here we show that the agreement is further improved when the nonlinear PB equation is used. These predictions are obtained with the dielectric boundary defined as the protein van der Waals surface. When the dielectric boundary is instead specified as the molecular surface, electrostatic interactions in the transient complex become repulsive and are thus predicted to retard association. Together these results demonstrate that the transient-complex theory is predictive of electrostatic rate enhancement and can help parameterize PB calculations.

Prediction of salt and mutational effects on the association rate of U1A protein and U1 small nuclear RNA stem/loop II

We have developed a computational approach for predicting protein-protein association rates (Alsallaq and Zhou, Structure 2007, 15, 215). Here we expand the range of applicability of this approach to protein-RNA binding and report the first results for protein-RNA binding rates predicted from atomistic modeling. The system studied is the U1A protein and stem/loop II of the U1 small nuclear RNA. Experimentally it was observed that the binding rate is significantly reduced by increasing salt concentration while the dissociation changes little with salt concentration, and charges distant from the binding site make marginal contribution to the binding rate. These observations are rationalized. Moreover, predicted effects of salt and charge mutations are found to be in quantitative agreement with experimental results.

Evaluation of the influence of the internal aqueous solvent structure on electrostatic interactions at the protein-solvent interface by nonlocal continuum electrostatic approach

The dielectric properties of the polar solvent on the protein-solvent interface at small intercharge distances are still poorly explored. To deconvolute this problem and to evaluate the pair-wise electrostatic interaction (PEI) energies of the point charges located at the protein-solvent interface we used a nonlocal (NL) electrostatic approach along with a static NL dielectric response function of water. The influence of the aqueous solvent microstructure (determined by a strong nonelectrostatic correlation effect between water dipoles within the orientational Debye polarization mode) on electrostatic interactions at the interface was studied in our work. It was shown that the PEI energies can be significantly higher than the energies evaluated by the classical (local) consideration, treating water molecules as belonging to the bulk solvent with a high dielectric constant. Our analysis points to the existence of a rather extended, effective low-dielectric interfacial water shell on the protein surface. The main dielectric properties of this shell (effective thickness together with distance- and orientation-dependent dielectric permittivity function) were evaluated. The dramatic role of this shell was demonstrated when estimating the protein association rate constants.

Predicted effects of missense mutations on native-state stability account for phenotypic outcome in phenylketonuria, a paradigm of misfolding diseases

Phenylketonuria (PKU) is a genetic disease caused by mutations in human phenylalanine hydroxylase (PAH). Most missense mutations result in misfolding of PAH, increased protein turnover, and a loss of enzymatic function. We studied the prediction of the energetic impact on PAH native-state stability of 318 PKU-associated missense mutations, using the protein-design algorithm FoldX. For the 80 mutations for which expression analyses have been performed in eukaryote systems, in most cases we found substantial overall correlations between the mutational energetic impact and both in vitro residual activities and patient metabolic phenotype. This finding confirmed that the decrease in protein stability is the main molecular pathogenic mechanism in PKU and the determinant for phenotypic outcome. Metabolic phenotypes have been shown to be better predicted than in vitro residual activities, probably because of greater stringency in the phenotyping process. Finally, all the remaining 238 PKU missense mutations compiled at the PAH locus knowledgebase (PAHdb) were analyzed, and their phenotypic outcomes were predicted on the basis of the energetic impact provided by FoldX. Residues in exons 7-9 and in interdomain regions within the subunit appear to play an important structural role and constitute hotspots for destabilization. FoldX analysis will be useful for predicting the phenotype associated with rare or new mutations detected in patients with PKU. However, additional factors must be considered that may contribute to the patient phenotype, such as possible effects on catalysis and interindividual differences in physiological and metabolic processes.