Correlation between binding and dynamics at SH2 domain interfaces

Measurement of 1Hα transverse relaxation rates in proteins: application to solvent PREs

It has recently been demonstrated that accurate near surface electrostatic potentials can be calculated for proteins from solvent paramagnetic relaxation enhancements (PREs) of amide protons measured using spin labels of similar structures but different charges (Yu et al. in Proc Natl Acad Sci 118(25):e2104020118, 2021). Here we develop methodology for extending such measurements to intrinsically disordered proteins at neutral pH where amide spectra are of very poor quality. Under these conditions it is shown that accurate PRE values can be measured using the haCONHA experiment that has been modified for recording ¹H^α transverse relaxation rates. The optimal pulse scheme includes a spin-lock relaxation element for suppression of homonuclear scalar coupled evolution for all ¹H^α protons, except those derived from Ser and Thr residues, and minimizes the radiation damping field from water magnetization that would otherwise increase measured relaxation rates. The robustness of the experiment is verified by developing a second approach using a band selective adiabatic decoupling scheme for suppression of scalar coupling modulations during ¹H^α relaxation and showing that the measured PRE values from the two methods are in excellent agreement. The near surface electrostatic potential of a 103-residue construct comprising the C-terminal intrinsically disordered region of the RNA-binding protein CAPRIN1 is obtained at pH 5.5 using both ¹H^N and ¹H^α-based relaxation rates, and at pH 7.4 where only ¹H^α rates can be quantified, with very good agreement between potentials obtained under all experimental conditions.

Targeting Oncogenic Src Homology 2 Domain-Containing Phosphatase 2 (SHP2) by Inhibiting Its Protein-Protein Interactions

We developed a new class of inhibitors of protein-protein interactions of the SHP2 phosphatase, which is pivotal in cell signaling and represents a central target in the therapy of cancer and rare diseases. Currently available SHP2 inhibitors target the catalytic site or an allosteric pocket but lack specificity or are ineffective for disease-associated SHP2 mutants. Considering that pathogenic lesions cause signaling hyperactivation due to increased levels of SHP2 association with cognate proteins, we developed peptide-based molecules with nanomolar affinity for the N-terminal Src homology domain of SHP2, good selectivity, stability to degradation, and an affinity for pathogenic variants of SHP2 that is 2-20 times higher than for the wild-type protein. The best peptide reverted the effects of a pathogenic variant (D61G) in zebrafish embryos. Our results provide a novel route for SHP2-targeted therapies and a tool for investigating the role of protein-protein interactions in the function of SHP2.

Correlated Response of Protein Side-Chain Fluctuations and Conformational Entropy to Ligand Binding

The heterogeneous fast side-chain dynamics of proteins plays crucial roles in molecular recognition and binding. Site-specific NMR experiments quantify these motions by measuring the model-free order parameter (O_axis²) on a scale of 0 (most flexible) to 1 (least flexible) for each methyl-containing residue of proteins. Here, we have examined ligand-induced variations in the fast side-chain dynamics and conformational entropy of calmodulin (CaM) using five different CaM-peptide complexes. O_axis² of CaM in the ligand-free (O_axis,U²) and ligand-bound (O_axis,B²) states are calculated from molecular dynamics trajectories and conformational energy surfaces obtained using the adaptive biasing force (ABF) method. ΔO_axis² = O_axis,B² - O_axis,U² follows a Gaussian-like unimodal distribution whose second moment is a potential indicator of the binding affinity of these complexes. The probability for the binding-induced O_axis,U² → O_axis,B² transition decreases with increasing magnitude of ΔO_axis², indicating that large flexibility changes are improbable for side chains of CaM after ligand binding. A linear correlation established between ΔO_axis² and the conformational entropy change of the protein makes possible the determination of the conformational entropy of binding of protein-ligand complexes. The results not only underscore the functional importance of fast side-chain fluctuations but also highlight key motional and thermodynamic correlates of protein-ligand binding.

Backbone assignments and conformational dynamics in the S. typhimurium tryptophan synthase α-subunit from solution-state NMR

Backbone assignments for the isolated α-subunit of Salmonella typhimurium tryptophan synthase (TS) are reported based on triple resonance solution-state NMR experiments on a uniformly 2H,13C,15N-labeled sample. From the backbone chemical shifts, secondary structure and random coil index order parameters (RCI-S2) are predicted. Titration with the 3-indole-D-glycerol 3'-phosphate analog, N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), leads to chemical shift perturbations indicative of conformational changes from which an estimate of the dissociation constant is obtained. Comparisons of the backbone chemical-shifts, RCI-S2 values, and site-specific relaxation times with and without F9 reveal allosteric changes including modulation in secondary structures and loop rigidity induced upon ligand binding. A comparison is made to the X-ray crystal structure of the α-subunit in the full TS αββα bi-enzyme complex and to two new X-ray crystal structures of the isolated TS α-subunit reported in this work.

Specificity and regulation of phosphotyrosine signaling through SH2 domains

Phosphotyrosine (pY) signaling is instrumental to numerous cellular processes. pY recognition occurs through specialized protein modules, among which the Src-homology 2 (SH2) domain is the most common. SH2 domains are small protein modules with an invariant fold, and are present in more than a hundred proteins with different function. Here we ask the question of how such a structurally conserved, small protein domain can recognize distinct phosphopeptides with the breath of binding affinity, specificity and kinetic parameters necessary for proper control of pY-dependent signaling and rapid cellular response. We review the current knowledge on structure, thermodynamics and kinetics of SH2-phosphopeptide complexes and conclude that selective phosphopeptide recognition is governed by both structure and dynamics of the SH2 domain, as well as by the kinetics of the binding events. Further studies on the thermodynamic and kinetic properties of SH2-phosphopeptide complexes, beyond their structure, are required to understand signaling regulation.

Dynamically Coupled Residues within the SH2 Domain of FYN Are Key to Unlocking Its Activity

Src kinase activity is controlled by various mechanisms involving a coordinated movement of kinase and regulatory domains. Notwithstanding the extensive knowledge related to the backbone dynamics, little is known about the more subtle side-chain dynamics within the regulatory domains and their role in the activation process. Here, we show through experimental methyl dynamic results and predicted changes in side-chain conformational couplings that the SH2 structure of Fyn contains a dynamic network capable of propagating binding information. We reveal that binding the phosphorylated tail of Fyn perturbs a residue cluster near the linker connecting the SH2 and SH3 domains of Fyn, which is known to be relevant in the regulation of the activity of Fyn. Biochemical perturbation experiments validate that those residues are essential for inhibition of Fyn, leading to a gain of function upon mutation. These findings reveal how side-chain dynamics may facilitate the allosteric regulation of the different members of the Src kinase family.

Analysis of Multidomain Protein Dynamics

Proteins with a modular architecture of multiple domains connected by linkers often exhibit diversity in the relative positions of domains, while the domain tertiary structure remains unchanged. The biological function of these modular proteins, or the regulation of their activity, depends on the variation in domain orientation and separation. Accordingly, careful characterization of interdomain motion and correlated fluctuations of multidomain systems is relevant for understanding the functional behavior of modular proteins. Molecular dynamics (MD) simulations provides a powerful approach to study these motions in atomic detail. Nevertheless, the common procedure for analyzing fluctuations from MD simulations after rigid-body alignment fails for multidomain proteins; it greatly overestimates correlated positional fluctuations in the presence of relative domain motion. We show here that expressing the atomic motions of a multidomain protein as a combination of displacement within the domain reference frame and motion of the relative domains correctly separates the internal motions to allow a useful description of correlated fluctuations. We illustrate the methodology of separating the domain fluctuations and local fluctuations by application to the tandem SH2 domains of human Syk protein kinase and by characterizing an effect of phosphorylation on the dynamics. Correlated motions are assessed from a distance covariance rather than the more common vector-coordinate covariance. The approach makes it possible to calculate the proper correlations in fluctuations internal to a domain as well as between domains.

Protein conformational exchange measured by 1H R1ρ relaxation dispersion of methyl groups

Activated dynamics plays a central role in protein function, where transitions between distinct conformations often underlie the switching between active and inactive states. The characteristic time scales of these transitions typically fall in the microsecond to millisecond range, which is amenable to investigations by NMR relaxation dispersion experiments. Processes at the faster end of this range are more challenging to study, because higher RF field strengths are required to achieve refocusing of the exchanging magnetization. Here we describe a rotating-frame relaxation dispersion experiment for (1)H spins in methyl (13)CHD2 groups, which improves the characterization of fast exchange processes. The influence of (1)H-(1)H rotating-frame nuclear Overhauser effects (ROE) is shown to be negligible, based on a comparison of R 1ρ relaxation data acquired with tilt angles of 90° and 35°, in which the ROE is maximal and minimal, respectively, and on samples containing different (1)H densities surrounding the monitored methyl groups. The method was applied to ubiquitin and the apo form of calmodulin. We find that ubiquitin does not exhibit any (1)H relaxation dispersion of its methyl groups at 10 or 25 °C. By contrast, calmodulin shows significant conformational exchange of the methionine methyl groups in its C-terminal domain, as previously demonstrated by (1)H and (13)C CPMG experiments. The present R 1ρ experiment extends the relaxation dispersion profile towards higher refocusing frequencies, which improves the definition of the exchange correlation time, compared to previous results.

The description of protein internal motions aids selection of ligand binding poses by the INPHARMA method

Protein internal motions influence observables of NMR experiments. The effect of internal motions occurring at the sub-nanosecond timescale can be described by NMR order parameters. Here, we report that the use of order parameters derived from Molecular Dynamics (MD) simulations of two holo-structures of Protein Kinase A increase the discrimination power of INPHARMA, an NMR based methodology that selects docked ligand orientations by maximizing the correlation of back-calculated to experimental data. By including internal motion in the back-calculation of the INPHARMA transfer, we obtain a more realistic description of the system, which better represents the experimental data. Furthermore, we propose a set of generic order parameters, derived from MD simulations of globular proteins, which can be used in the back-calculation of INPHARMA NOEs for any protein-ligand complex, thus by-passing the need of obtaining system-specific order parameters for new protein-ligand complexes.

Probing side-chain dynamics in proteins by the measurement of nine deuterium relaxation rates per methyl group

We demonstrate the feasibility of the measurement of up to nine deuterium spin relaxation rates in 13CHD2 and 13CH2D methyl isotopomers of small proteins. In addition to five measurable 2H relaxation rates in a 13CH2D methyl group (Millet, O.; Muhandiram, D. R.; Skrynnikov, N. R.; Kay, L. E. J. Am. Chem. Soc. 2002, 124, 6439-48), the measurement of additional four rates of (nearly) single-exponentially decaying magnetization terms in methyl groups of the 13CHD2 variety is reported. Consistency relationships between 2H spin relaxation rates measured in the two different types of methyl groups are derived and verified experimentally for a subset of methyl-containing side chains in the protein ubiquitin. A detailed comparison of methyl-bearing side-chain dynamics parameters obtained from relaxation measurements in 13CH2D and 13CHD2 methyls of ubiquitin at 10, 27, and 40 °C reveals that transverse 2H relaxation rates in 13CHD2 groups are reliable and accurate reporters of the amplitudes of methyl 3-fold axis motions (S(axis)2) for protein molecules with global molecular tumbling times τ(C) >~9 ns. For smaller molecules, simple correction of transverse 2H relaxation rates in 13CHD2 groups is sufficient for the derivation of robust measures of order. Residue-specific distributions of S(axis)2 are consistent with atomic-detail molecular dynamics (MD) results. Both 13CHD2- and 13CH2D-derived S(axis)2 values are in good overall agreement with those obtained from 1 μs MD simulations at all the three temperatures, although some differences in the site-specific temperature dependence between MD- and 2H-relaxation-derived S(axis)2 values are observed.

Toward a predictive understanding of slow methyl group dynamics in proteins

The development of the most recent generation of molecular mechanics force fields promises an increasingly predictive understanding of the protein dynamics-function relationship. Based on extensive validation against various types of experimental data, the AMBER force field ff99SB was benchmarked in recent years as a favorable force field for protein simulations. Recent improvements of the side chain and backbone potentials, made by different groups, led to the ff99SB-ILDN and ff99SBnmr1 force fields, respectively. The combination of these potentials into a unified force field, termed ff99SBnmr1-ILDN, was used in this study to perform a microsecond time scale molecular dynamics simulation of free ubiquitin in explicit solvent for validation against an extensive set of experimental NMR methyl group residual dipolar couplings. Our results show a high level of consistency between the experimental data and the values predicted from the molecular dynamics trajectory reflecting a systematically improved performance of ff99SBnmr1-ILDN over the original ff99SB force field. Moreover, the unconstrained ff99SBnmr1-ILDN MD ensemble achieves a similar level of agreement as the recently introduced EROS ensemble, which was constructed based on a large body of NMR data as constraints, including the methyl residual dipolar couplings. This suggests that ff99SBnmr1-ILDN provides a high-quality representation of the motions of methyl-bearing protein side chains, which are sensitive probes of protein-protein and protein-ligand interactions.

Protein-ligand interactions: thermodynamic effects associated with increasing nonpolar surface area

Thermodynamic parameters were determined for complex formation between the Grb2 SH2 domain and Ac-pTyr-Xaa-Asn derived tripeptides in which the Xaa residue is an α,α-cycloaliphatic amino acid that varies in ring size from three- to seven-membered. Although the six- and seven-membered ring analogs are approximately equipotent, binding affinities of those having three- to six-membered rings increase incrementally with ring size because increasingly more favorable binding enthalpies dominate increasingly less favorable binding entropies, a finding consistent with an enthalpy-driven hydrophobic effect. Crystallographic analysis reveals that the only significant differences in structures of the complexes are in the number of van der Waals contacts between the domain and the methylene groups in the Xaa residues. There is a positive correlation between buried nonpolar surface area and binding free energy and enthalpy, but not with ΔC(p). Displacing a water molecule from a protein-ligand interface is not necessarily reflected in a favorable change in binding entropy. These findings highlight some of the fallibilities associated with commonly held views of relationships of structure and energetics in protein-ligand interactions and have significant implications for ligand design.

Mathematical treatment of adiabatic fast passage pulses for the computation of nuclear spin relaxation rates in proteins with conformational exchange

Although originally designed for broadband inversion and decoupling in NMR spectroscopy, recent methodological developments have introduced adiabatic fast passage (AFP) pulses into the field of protein dynamics. AFP pulses employ a frequency sweep, and have not only superior inversion properties with respect to offset effects, but they are also easily implemented into a pulse sequence. As magnetization is dragged from the +z to the -z direction, Larmor precession is impeded since magnetization becomes spin-locked, which is a potentially useful feature for the investigation of microsecond to millisecond dynamics. A major drawback of these pulses as theoretical prediction is concerned, however, results from their time-dependent offset: simulations of spin density matrices under the influence of a time-dependent Hamiltonian with non-commuting elements are costly in terms of computational time, rendering data analysis impracticable. In this paper we suggest several ways to reduce the computational time without compromising accuracy with respect to effects such as cross-correlated relaxation and modulation of the chemical shift.

Protein dynamics and allostery: an NMR view

Allostery, the process by which distant sites within a protein system are energetically coupled, is an efficient and ubiquitous mechanism for activity regulation. A purely mechanical view of allostery invoking only structural changes has developed over the decades as the classical view of the phenomenon. However, a fast growing list of examples illustrate the intimate link between internal motions over a wide range of time scales and function in protein-ligand interactions. Proteins respond to perturbations by redistributing their motions and they use fluctuating conformational states for binding and conformational entropy as a carrier of allosteric energy to modulate association with ligands. In several cases allosteric interactions proceed with minimal or no structural changes. We discuss emerging paradigms for the central role of protein dynamics in allostery.

Solution structure and dynamics of the I214V mutant of the rabbit prion protein

BACKGROUND

The conformational conversion of the host-derived cellular prion protein (PrP(C)) into the disease-associated scrapie isoform (PrP(Sc)) is responsible for the pathogenesis of transmissible spongiform encephalopathies (TSEs). Various single-point mutations in PrP(C)s could cause structural changes and thereby distinctly influence the conformational conversion. Elucidation of the differences between the wild-type rabbit PrP(C) (RaPrP(C)) and various mutants would be of great help to understand the ability of RaPrP(C) to be resistant to TSE agents.

METHODOLOGY/PRINCIPAL FINDINGS

We determined the solution structure of the I214V mutant of RaPrP(C)(91-228) and detected the backbone dynamics of its structured C-terminal domain (121-228). The I214V mutant displays a visible shift of surface charge distribution that may have a potential effect on the binding specificity and affinity with other chaperones. The number of hydrogen bonds declines dramatically. Urea-induced transition experiments reveal an obvious decrease in the conformational stability. Furthermore, the NMR dynamics analysis discloses a significant increase in the backbone flexibility on the pico- to nanosecond time scale, indicative of lower energy barrier for structural rearrangement.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that both the surface charge distribution and the intrinsic backbone flexibility greatly contribute to species barriers for the transmission of TSEs, and thereby provide valuable hints for understanding the inability of the conformational conversion for RaPrP(C).

NMR backbone dynamics of VEK-30 bound to the human plasminogen kringle 2 domain

To gain insights into the mechanisms for the tight and highly specific interaction of the kringle 2 domain of human plasminogen (K2(Pg)) with a 30-residue internal peptide (VEK-30) from a group A streptococcal M-like protein, the dynamic properties of free and bound K2(Pg) and VEK-30 were investigated using backbone amide (15)N-NMR relaxation measurements. Dynamic parameters, namely the generalized order parameter, S(2), the local correlation time, tau(e), and the conformational exchange contribution, R(ex), were obtained for this complex by Lipari-Szabo model-free analysis. The results show that VEK-30 displays distinctly different dynamic behavior as a consequence of binding to K2(Pg), manifest by decreased backbone flexibility, particularly at the binding region of the peptide. In contrast, the backbone dynamics parameters of K2(Pg) displayed similar patterns in the free and bound forms, but, nonetheless, showed interesting differences. Based on our previous structure-function studies of this interaction, we also made comparisons of the VEK-30/K2(Pg) dynamics results from different kringle modules complexed with small lysine analogs. The differences in dynamics observed for kringles with different ligands provide what we believe to be new insights into the interactions responsible for protein-ligand recognition and a better understanding of the differences in binding affinity and binding specificity of kringle domains with various ligands.

Loops govern SH2 domain specificity by controlling access to binding pockets

Cellular functions require specific protein-protein interactions that are often mediated by modular domains that use binding pockets to engage particular sequence motifs in their partners. Yet, how different members of a domain family select for distinct sequence motifs is not fully understood. The human genome encodes 120 Src homology 2 (SH2) domains (in 110 proteins), which mediate protein-protein interactions by binding to proteins with diverse phosphotyrosine (pTyr)-containing sequences. The structure of the SH2 domain of BRDG1 bound to a peptide revealed a binding pocket that was blocked by a loop residue in most other SH2 domains. Analysis of 63 SH2 domain structures suggested that the SH2 domains contain three binding pockets, which exhibit selectivity for the three positions after the pTyr in a peptide, and that SH2 domain loops defined the accessibility and shape of these pockets. Despite sequence variability in the loops, we identified conserved structural features in the loops of SH2 domains responsible for controlling access to these surface pockets. We engineered new loops in an SH2 domain that altered specificity as predicted. Thus, selective blockage of binding subsites or pockets by surface loops provides a molecular basis by which the diverse modes of ligand recognition by the SH2 domain may have evolved and provides a framework for engineering SH2 domains and designing SH2-specific inhibitors.

Measurement of signs of chemical shift differences between ground and excited protein states: a comparison between H(S/M)QC and R1rho methods

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR spectroscopy has emerged as a powerful tool for quantifying the kinetics and thermodynamics of millisecond exchange processes between a major, populated ground state and one or more minor, low populated and often invisible 'excited' conformers. Analysis of CPMG data-sets also provides the magnitudes of the chemical shift difference(s) between exchanging states (|Deltavarpi|), that inform on the structural properties of the excited state(s). The sign of Deltavarpi is, however, not available from CPMG data. Here we present one-dimensional NMR experiments for measuring the signs of (1)H(N) and (13)C(alpha) Deltavarpi values using weak off-resonance R (1rho ) relaxation measurements, extending the spin-lock approach beyond previous applications focusing on the signs of (15)N and (1)H(alpha) shift differences. The accuracy of the method is established by using an exchanging system where the invisible, excited state can be converted to the visible, ground state by altering conditions so that the signs of Deltavarpi values obtained from the spin-lock approach can be validated with those measured directly. Further, the spin-lock experiments are compared with the established H(S/M)QC approach for measuring the signs of chemical shift differences. For the Abp1p and Fyn SH3 domains considered here it is found that while H(S/M)QC measurements provide signs for more residues than the spin-lock data, the two different methodologies are complementary, so that combining both approaches frequently produces signs for more residues than when the H(S/M)QC method is used alone.

Dynamic activation of an allosteric regulatory protein

Allosteric regulation is used as a very efficient mechanism to control protein activity in most biological processes, including signal transduction, metabolism, catalysis and gene regulation. Allosteric proteins can exist in several conformational states with distinct binding or enzymatic activity. Effectors are considered to function in a purely structural manner by selectively stabilizing a specific conformational state, thereby regulating protein activity. Here we show that allosteric proteins can be regulated predominantly by changes in their structural dynamics. We have used NMR spectroscopy and isothermal titration calorimetry to characterize cyclic AMP (cAMP) binding to the catabolite activator protein (CAP), a transcriptional activator that has been a prototype for understanding effector-mediated allosteric control of protein activity. cAMP switches CAP from the 'off' state (inactive), which binds DNA weakly and non-specifically, to the 'on' state (active), which binds DNA strongly and specifically. In contrast, cAMP binding to a single CAP mutant, CAP-S62F, fails to elicit the active conformation; yet, cAMP binding to CAP-S62F strongly activates the protein for DNA binding. NMR and thermodynamic analyses show that despite the fact that CAP-S62F-cAMP(2) adopts the inactive conformation, its strong binding to DNA is driven by a large conformational entropy originating in enhanced protein motions induced by DNA binding. The results provide strong evidence that changes in protein motions may activate allosteric proteins that are otherwise structurally inactive.

Ligand-induced changes in the structure and dynamics of Escherichia coli peptide deformylase

Peptide deformylase (PDF) is an enzyme that is responsible for removing the formyl group from nascently synthesized polypeptides in bacteria, attracting much attention as a potential target for novel antibacterial agents. Efforts to develop potent inhibitors of the enzyme have progressed on the basis of classical medicinal chemistry, combinatorial chemistry, and structural approaches, yet the validity of PDF as an antibacterial target hangs, in part, on the ability of inhibitors to selectively target this enzyme in favor of structurally related metallohydrolases. We have used (15)N NMR spectroscopy and isothermal titration calorimetry to investigate the high-affinity interaction of EcPDF with actinonin, a naturally occurring potent EcPDF inhibitor. Backbone amide chemical shifts, residual dipolar couplings, hydrogen-deuterium exchange, and (15)N relaxation reveal structural and dynamic effects of ligand binding in the immediate vicinity of the ligand-binding site as well as at remote sites. A comparison of the crystal structures of free and actinonin-bound EcPDF with the solution data suggests that most of the consequences of the ligand binding to the protein are lost or obscured during crystallization. The results of these studies improve our understanding of the thermodynamic global minimum and have important implications for structure-based drug design.

Structure, dynamics, and binding thermodynamics of the v-Src SH2 domain: implications for drug design

SH2 domains provide fundamental recognition sites in tyrosine kinase-mediated signaling pathways which, when aberrant, give rise to disease states such as cancer, diabetes, and immune deficiency. Designing specific inhibitors that target the SH2 domain-binding site, however, have presented a major challenge. Despite well over a decade of intensive research, clinically useful SH2 domain inhibitors have yet to become available. A better understanding of the structural, dynamic, and thermodynamic contributions to ligand binding of individual SH2 domains will provide some insight as to whether inhibitor development is possible. We report the first high resolution solution structure of the apo-v-Src SH2 domain. This is accompanied by the analysis of backbone dynamics and pK(a) values within the apo- and peptide-bound states. Our results indicate that the phosphotyrosine (pY) pocket is tightly structured and hence not adaptable to exogenous ligands. On the other hand, the pocket which accommodates residues proximal and C-terminal of the pY (pY + 3) or so-called specificity determining region, is a large dynamic-binding surface. This appears to allow a high level of promiscuity in binding. Binding of a series of synthetic, phosphotyrosyl, peptidomimetic compounds designed to explore interactions in the pY + 3 pocket further demonstrates the ability of the SH2 domain to accommodate diverse ligands. The thermodynamic parameters of these interactions show dramatic enthalpy/entropy compensation. These data suggest that the v-Src SH2 domain does not have a highly specific secondary-binding site, which clearly presents a major hurdle to design selective inhibitors.

Furanose dynamics in the HhaI methyltransferase target DNA studied by solution and solid-state NMR relaxation

Both solid-state and solution NMR relaxation measurements are routinely used to quantify the internal dynamics of biomolecules, but in very few cases have these two techniques been applied to the same system, and even fewer attempts have been made so far to describe the results obtained through these two methods through a common theoretical framework. We have previously collected both solution 13C and solid-state 2H relaxation measurements for multiple nuclei within the furanose rings of several nucleotides of the DNA sequence recognized by HhaI methyltransferase. The data demonstrated that the furanose rings within the GCGC recognition sequence are very flexible, with the furanose rings of the cytidine, which is the methylation target, experiencing the most extensive motions. To interpret these experimental results quantitatively, we have developed a dynamic model of furanose rings based on the analysis of solid-state 2H line shapes. The motions are modeled by treating bond reorientations as Brownian excursions within a restoring potential. By applying this model, we are able to reproduce the rates of 2H spin-lattice relaxation in the solid and 13C spin-lattice relaxation in solution using comparable restoring force constants and internal diffusion coefficients. As expected, the 13C relaxation rates in solution are less sensitive to motions that are slower than overall molecular tumbling than to the details of global molecular reorientation, but are somewhat more sensitive to motions in the immediate region of the Larmor frequency. Thus, we conclude that the local internal motions of this DNA oligomer in solution and in the hydrated solid state are virtually the same, and we validate an approach to the conjoint analysis of solution and solid-state NMR relaxation and line shapes data, with wide applicability to many biophysical problems.

Mechanistic insights into phosphoprotein-binding FHA domains

[Structure: see text]. FHA domains are protein modules that switch signals in diverse biological pathways by monitoring the phosphorylation of threonine residues of target proteins. As part of the effort to gain insight into cellular avoidance of cancer, FHA domains involved in the cellular response to DNA damage have been especially well-characterized. The complete protein where the FHA domain resides and the interaction partners determine the nature of the signaling. Thus, a key biochemical question is how do FHA domains pick out their partners from among thousands of alternatives in the cell? This Account discusses the structure, affinity, and specificity of FHA domains and the formation of their functional structure. Although FHA domains share sequence identity at only five loop residues, they all fold into a beta-sandwich of two beta-sheets. The conserved arginine and serine of the recognition loops recognize the phosphorylation of the threonine targeted. Side chains emanating from loops that join beta-strand 4 with 5, 6 with 7, or 10 with 11 make specific contacts with amino acids of the ligand that tailor sequence preferences. Many FHA domains choose a partner in extended conformation, somewhat according to the residue three after the phosphothreonine in sequence (pT + 3 position). One group of FHA domains chooses a short carboxylate-containing side chain at pT + 3. Another group chooses a long, branched aliphatic side chain. A third group prefers other hydrophobic or uncharged polar side chains at pT + 3. However, another FHA domain instead chooses on the basis of pT - 2, pT - 3, and pT + 1 positions. An FHA domain from a marker of human cancer instead chooses a much longer protein fragment that adds a beta-strand to its beta-sheet and that presents hydrophobic residues from a novel helix to the usual recognition surface. This novel recognition site and more remote sites for the binding of other types of protein partners were predicted for the entire family of FHA domains by a bioinformatics approach. The phosphopeptide-dependent dynamics of an FHA domain, SH2 domain, and PTB domain suggest a common theme: rigid, preformed binding surfaces support van der Waals contacts that provide favorable binding enthalpy. Despite the lack of pronounced conformational changes in FHA domains linked to binding events, more subtle adjustments may be possible. In the one FHA domain tested, phosphothreonine peptide binding is accompanied by increased flexibility just outside the binding site and increased rigidity across the beta-sandwich. The folding of the same FHA domain progresses through near-native intermediates that stabilize the recognition loops in the center of the phosphoprotein-binding surface; this may promote rigidity in the interface and affinity for targets phosphorylated on threonine.

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.

Quantifying Lipari-Szabo modelfree parameters from 13CO NMR relaxation experiments

It is proposed to obtain effective Lipari-Szabo order parameters and local correlation times for relaxation vectors of protein (13)CO nuclei by carrying out a (13)CO-R(1) auto relaxation experiment, a transverse (13)CO CSA/13CO-13Calpha CSA/dipolar cross correlation and a transverse (13)CO CSA/(13)CO-(15)N CSA/dipolar cross correlation experiment. Given the global rotational correlation time from (15)N relaxation experiments, a new program COMFORD (CO-Modelfree Fitting Of Relaxation Data) is presented to fit the (13)CO data to an effective order parameter S2CO, an effective local correlation time and the orientation of the CSA tensor with respect to the molecular frame. It is shown that the effective S2CO is least sensitive to rotational fluctuations about an imaginary Calpha-Calpha axis and most sensitive to rotational fluctuations about an imaginary axis parallel to the NH bond direction. As such, the Calpha-Calpha information is fully complementary to the (15)N relaxation order parameter, which is least sensitive to fluctuations about the NH axis and most sensitive to fluctuations about the Calpha-Calpha axis. The new paradigm is applied on data of Ca(2+) saturated Calmodulin, and on available literature data for Ubiquitin. Our data indicate that the S2CO order parameters rapport on slower, and sometimes different, motions than the (15)N relaxation order parameters. The CO local correlation times correlate well with the calmodulin's secondary structure.

Evaluation of energetic and dynamic coupling networks in a PDZ domain protein

A number of computational and experimental studies have identified intramolecular communication "pathways" or "networks" important for transmitting allostery. Here, we have used mutagenesis and NMR relaxation methods to investigate the scope and nature of the communication networks found in the second post-synaptic density-95/discs large/zonula occludens-1 (PDZ) domain of the human protein tyrosine phosphatase 1E protein (hPTP1E) (PDZ2). It was found that most mutations do not have a significant energetic contribution to peptide ligand binding. Three mutants that showed significant changes in binding also displayed context-dependent dynamic effects. Both a mutation at a partially exposed site (H71Y) and a buried core position (I35V) had a limited response in side-chain (2)H-based dynamics when compared to wild-type PDZ2. In contrast, a change at a second core position (I20F) that had previously been shown to be part of an energetic and dynamic network, resulted in extensive changes in side-chain dynamics. This response is reminiscent to that seen previously upon peptide ligand binding. These results shed light on the nature of the PDZ2 dynamic network and suggest that position 20 in PDZ2 acts as a "hub" that is energetically and dynamically critical for transmitting changes in dynamics throughout the PDZ domain.

Multiscale modeling of macromolecular conformational changes combining concepts from rigidity and elastic network theory

The development of a two-step approach for multiscale modeling of macromolecular conformational changes is based on recent developments in rigidity and elastic network theory. In the first step, static properties of the macromolecule are determined by decomposing the molecule into rigid clusters by using the graph-theoretical approach FIRST and an all-atom representation of the protein. In this way, rigid clusters are not limited to consist of residues adjacent in sequence or secondary structure elements as in previous studies. Furthermore, flexible links between rigid clusters are identified and can be modeled as such subsequently. In the second step, dynamical properties of the molecule are revealed by the rotations-translations of blocks approach (RTB) using an elastic network model representation of the coarse-grained protein. In this step, only rigid body motions are allowed for rigid clusters, whereas links between them are treated as fully flexible. The approach was tested on a data set of 10 proteins that showed conformational changes on ligand binding. For efficiency, coarse-graining the protein results in a remarkable reduction of memory requirements and computational times by factors of 9 and 27 on average and up to 25 and 125, respectively. For accuracy, directions and magnitudes of motions predicted by our approach agree well with experimentally determined ones, despite embracing in extreme cases >50% of the protein into one rigid cluster. In fact, the results of our method are in general comparable with when no or a uniform coarse-graining is applied; and the results are superior if the movement is dominated by loop or fragment motions. This finding indicates that explicitly distinguishing between flexible and rigid regions is advantageous when using a simplified protein representation in the second step. Finally, motions of atoms in rigid clusters are also well predicted by our approach, which points to the need to consider mobile protein regions in addition to flexible ones when modeling correlated motions.

Relation between native ensembles and experimental structures of proteins

Different experimental structures of the same protein or of proteins with high sequence similarity contain many small variations. Here we construct ensembles of "high-sequence similarity Protein Data Bank" (HSP) structures and consider the extent to which such ensembles represent the structural heterogeneity of the native state in solution. We find that different NMR measurements probing structure and dynamics of given proteins in solution, including order parameters, scalar couplings, and residual dipolar couplings, are remarkably well reproduced by their respective high-sequence similarity Protein Data Bank ensembles; moreover, we show that the effects of uncertainties in structure determination are insufficient to explain the results. These results highlight the importance of accounting for native-state protein dynamics in making comparisons with ensemble-averaged experimental data and suggest that even a modest number of structures of a protein determined under different conditions, or with small variations in sequence, capture a representative subset of the true native-state ensemble.

NMR studies of Escherichia coli acyl carrier protein: dynamic and structural differences of the apo- and holo-forms

Two indicators of conformational variability of Escherichia coli acyl carrier protein (ACP) have been investigated, namely backbone dynamics and chemical shift variations of ACP. Hydrophobic interactions between the 4'-PP prosthetic group and the hydrophobic pocket enclosed by the amphipathic helices resulted in chemical shift perturbations in the residues near the prosthetic group binding sites and contact sites in the hydrophobic pockets upon conversion from apo- to holo-forms. At pH 7.9, destabilization of ACP due to negative charge repulsions and the deprotonated state of His 75 resulted in observed chemical shift changes in the C-terminal region. Model-free analysis showed that the alpha(1)alpha(2) loop region near the prosthetic group binding site in ACP shows the greatest flexibility (lowest S(2) values) and this result may suggest these flexibilities are required for structural rearrangements when the acyl chain binds to the prosthetic group of ACP. Flexibility of ACP shown in this study is essential for its ability to interact with functionally different enzyme partners specifically and weakly in the rapid delivery of acyl chain from one partner to another.

Structural basis for the requirement of two phosphotyrosine residues in signaling mediated by Syk tyrosine kinase

The protein-tyrosine kinase Syk couples immune recognition receptors to multiple signal transduction pathways, including the mobilization of calcium and the activation of NFAT. The ability of Syk to regulate signaling is influenced by its phosphorylation on tyrosine residues within the linker B region. The phosphorylation of both Y342 and Y346 is necessary for optimal signaling from the B cell receptor for antigen. The SH2 domains of multiple signaling proteins share the ability to bind this doubly phosphorylated site. The NMR structure of the C-terminal SH2 domain of PLCgamma (PLCC) bound to a doubly phosphorylated Syk peptide reveals a novel mode of phosphotyrosine recognition. PLCC undergoes extensive conformational changes upon binding to form a second phosphotyrosine-binding pocket in which pY346 is largely desolvated and stabilized through electrostatic interactions. The formation of the second binding pocket is distinct from other modes of phosphotyrosine recognition in SH2-protein association. The dependence of signaling on simultaneous phosphorylation of these two tyrosine residues offers a new mechanism to fine-tune the cellular response to external stimulation.

Effects of protein-pheromone complexation on correlated chemical shift modulations

Major urinary protein (MUP) is a pheromone-carrying protein of the lipocalin family. Previous studies by isothermal titration calorimetry (ITC) show that the affinity of MUP for the pheromone 2-methoxy-3-isobutylpyrazine (IBMP) is mainly driven by enthalpy, with a small unfavourable entropic contribution. Entropic terms can be attributed in part to changes in internal motions of the protein upon binding. Slow internal motions can lead to correlated or anti-correlated modulations of the isotropic chemical shifts of carbonyl C' and amide N nuclei. Correlated chemical shift modulations (CSM/CSM) in MUP have been determined by measuring differences of the transverse relaxation rates of zero- and double-quantum coherences ZQC{C'N} and DQC{C'N}, and by accounting for the effects of correlated fluctuations of dipole-dipole couplings (DD/DD) and chemical shift anisotropies (CSA/CSA). The latter can be predicted from tensor parameters of C' and N nuclei that have been determined in earlier work. The effects of complexation on slow time-scale protein dynamics can be determined by comparing the temperature dependence of the relaxation rates of APO-MUP (i.e., without ligand) and HOLO-MUP (i.e., with IBMP as a ligand).

PhosphoThr peptide binding globally rigidifies much of the FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

A net increase in the backbone rigidity of the kinase-interacting FHA domain (KI-FHA) from the Arabidopsis receptor kinase-associated protein phosphatase (KAPP) accompanies the binding of a phosphoThr peptide from its CLV1 receptor-like kinase partner, according to (15)N NMR relaxation at 11.7 and 14.1 T. All of the loops of free KI-FHA display evidence of nanosecond-scale motions. Many of these same residues have residual dipolar couplings that deviate from structural predictions. Binding of the CLV1 pT868 peptide seems to reduce nanosecond-scale fluctuations of all loops, including half of the residues of recognition loops. Residues important for affinity are found to be rigid, i.e., conserved residues and residues of the subsite for the key pT+3 peptide position. This behavior parallels SH2 and PTB domain recognition of pTyr peptides. PhosphoThr peptide binding increases KI-FHA backbone rigidity (S(2)) of three recognition loops, a loop nearby, seven strands from the beta-sandwich, and a distal loop. Compensating the trend of increased rigidity, binding enhances fast mobility at a few sites in four loops on the periphery of the recognition surface and in two loops on the far side of the beta-sandwich. Line broadening evidence of microsecond- to millisecond-scale fluctuations occurs across the six-stranded beta-sheet and nearby edges of the beta-sandwich; this forms a network connected by packing of interior side chains and H-bonding. A patch of the slowly fluctuating residues coincides with the site of segment-swapped dimerization in crystals of the FHA domain of human Chfr. Phosphopeptide binding introduces microsecond- to millisecond-scale fluctuations to more residues of the long 8/9 recognition loop of KI-FHA. The rigidity of this FHA domain appears to couple as a whole to pThr peptide binding.

What Contributions to Protein Side-chain Dynamics are Probed by NMR Experiments? A Molecular Dynamics Simulation Analysis

Molecular dynamics simulations of the structurally homologous proteins TNfn3 and FNfn10 have been used to investigate the contributions to side-chain dynamics measured by NMR relaxation experiments. The results reproduce the variation in core side-chain dynamics observed by NMR and highlight the relevance of anharmonic motion and transitions between local minima for explaining NMR side-chain order parameters. A method is described for calculating converged order parameters by use of replica exchange molecular dynamics in conjunction with an implicit solvent model. These simulations allow the influence of various factors, such as the flexibility of side-chains and their free volume, on the mobility to be tested by perturbing the system. Deletion mutations are found to have the largest effect on the more densely packed FNfn10. Some counterintuitive effects are seen, such as an increase in order parameters close to deletion mutation sites, but these can be rationalized in terms of direct interactions with the modified side-chains. A statistical analysis of published order parameters supports the conclusions drawn from the simulations.

NMR studies of protein structure and dynamics

Recent advances in solution NMR spectroscopy have significantly extended the spectrum of problems that can now be addressed with this technology. In particular, studies of proteins with molecular weights on the order of 100 kDa are now possible at a level of detail that was previously reserved for much smaller systems. An example of the sort of information that is now accessible is provided in a study of malate synthase G, a 723 residue enzyme that has been a focal point of research efforts in my laboratory. Details of the labeling schemes that have been employed and optimal experiments for extraction of structural and dynamics information on this protein are described. NMR studies of protein dynamics, in principle, give insight into the relation between motion and function. A description of deuterium-based spin relaxation methods for the investigation of side chain dynamics is provided. Examples where millisecond (ms) time scale dynamics play an important role and where relaxation dispersion NMR spectroscopy has been particularly informative, including applications involving the membrane enzyme PagP and mutants of the Fyn SH3 domain that fold on a ms time scale, are presented.

Structural insight into the binding diversity between the Tyr-phosphorylated human ephrinBs and Nck2 SH2 domain

The binding interaction between the Nck2 SH2 domain and the phosphorylated ephrinB initiates a critical pathway for the reverse signaling network mediated by Eph receptor-ephrinB. Previously, the NMR structure and Tyr phosphorylations of the human ephrinB cytoplasmic domain have been studied. To obtain a complete story, it would be of significant interest to determine the structure of the Nck2 SH2 domain that shows a low sequence identity to other SH2 domains with known structures. Here, we report the determination of the solution structure of the human Nck2 SH2 domain and investigate its interactions with three phosphorylated ephrinB fragments by NMR spectroscopy. The results indicate that: 1) although the human Nck2 SH2 domain adopts a core tertiary fold common to all SH2 domains, it owns some unique properties such as a shorter C-terminal helix and unusual electrostatic potential surface. However, the most striking finding is that the C-terminal tail of the human Nck2 SH2 domain adopts a short antiparallel beta-sheet that, to the best of our knowledge, has never been identified in other SH2 domains. The truncation study suggests that one function of the C-terminal tail is to control the folding/solubility of the SH2 domain. 2) In addition to [Tyr(P)304]ephrinB2(301-322) and [Tyr(P)316]ephrinB2(301-322), here we identified [Tyr(P)330]ephrinB2(324-333) also capable of binding to the SH2 domain. The detailed NMR study indicated that the binding mechanisms for the three ephrinB fragments might be different. The binding with [Tyr(P)304]-ephrinB2(301-322) and [Tyr(P)316]ephrinB2(301-322) might be mostly involved in the residues over the N-half of the SH2 domain and provoked a significant increase in the backbone and side chain dynamics of the SH2 domain on the microsecond-millisecond time scale. In contrast, binding with [Tyr(P)330]ephrinB2(324-333) might have most residues over both halves engaged but induced less profound conformational dynamics on the mus-ms time scale.

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.

Change in protein flexibility upon complex formation: Analysis of Ras-Raf using molecular dynamics and a molecular framework approach

Changes in flexibility upon protein-protein complex formation of H-Ras and the Ras-binding domain of C-Raf1 have been investigated using the molecular framework approach FIRST (Floppy Inclusion and Rigid Substructure Topology) and molecular dynamics simulations (MD) of in total approximately 35 ns length. In a computational time of about one second, FIRST identifies flexible and rigid regions in a single, static three-dimensional molecular framework, whose vertices represent protein atoms and whose edges represent covalent and non-covalent (hydrogen bond and hydrophobic) constraints and fixed bond angles within the protein. The two methods show a very good agreement with respect to the identification of changes in flexibility in both binding partners on a local scale. This implies that flexibility can be successfully predicted by identifying which bonds limit motion within a molecule and how they are coupled. In particular, as identified by MD, the beta-sheet in Raf shows considerably more pronounced orientational correlations in the bound state compared to the unbound state. Similarly, FIRST assigns the beta-sheet to the largest rigid cluster of the complex. Interestingly, FIRST allows us to identify that interactions across the interface (but not conformational changes upon complex formation) result in the observed rigidification. Since regions of the beta-sheet of Raf that do not interact directly with Ras become rigidified, this also demonstrates the long-range aspect to rigidity percolation. Possible implications of the change of flexibility of the Ras-binding domain of Raf on the activation of Raf upon complex formation are discussed. Finally, the sensitivity of FIRST results with respect to the representation of non-covalent interactions used as constraints is probed.

Ligand-dependent dynamics and intramolecular signaling in a PDZ domain

Allosteric communication is a fundamental process that proteins use to propagate signals from one site to functionally important distal sites. Although allostery is usually associated with multimeric proteins and enzymes, "long-range" communication may be a fundamental property of proteins. In some cases, communication occurs with minimal structural change. PDZ (post-synaptic density-95/discs large/zonula occludens-1) domains are small, protein-protein binding modules that can use multiple surfaces for docking diverse molecules. Furthermore, these domains have long-range energetic couplings that link the ligand-binding site to distal regions of the structure. Here, we show that allosteric behavior in a representative member of the PDZ domain family may be directly detected using side-chain methyl dynamics measurements. The changes in side-chain dynamics parameters in the second PDZ domain from the human tyrosine phosphatase 1E (hPTP1E) were determined upon binding a peptide target. Long-range dynamic effects were detected that correspond to previously observed pair-wise energetic couplings. These results provide one of the first experimental examples for the potential role of ps-ns timescale dynamics in propagating long-range signals within a protein, and reinforce the idea that dynamic fluctuations in proteins contribute to allosteric signal transduction.

Sequence, structure and energetic determinants of phosphopeptide selectivity of SH2 domains

Here, we present an approach for the prediction of binding preferences of members of a large protein family for which structural information for a number of family members bound to a substrate is available. The approach involves a number of steps. First, an accurate multiple alignment of sequences of all members of a protein family is constructed on the basis of a multiple structural superposition of family members with known structure. Second, the methods of continuum electrostatics are used to characterize the energetic contribution of each residue in a protein to the binding of its substrate. Residues that make a significant contribution are mapped onto the protein sequence and are used to define a "binding site signature" for the complex being considered. Third, sequences whose structures have not been determined are checked to see if they have binding-site signatures similar to one of the known complexes. Predictions of binding affinity to a given substrate are based on similarities in binding-site signature. An important component of the approach is the introduction of a context-specific substitution matrix suitable for comparison of binding-site residues. The methods are applied to the prediction of phosphopeptide selectivity of SH2 domains. To this end, the energetic roles of all protein residues in 17 different complexes of SH2 domains with their cognate targets are analyzed. The total number of residues that make significant contributions to binding is found to vary from nine to 19 in different complexes. These energetically important residues are found to contribute to binding through a variety of mechanisms, involving both electrostatic and hydrophobic interactions. Binding-site signatures are found to involve residues in different positions in SH2 sequences, some of them as far as 9A away from a bound peptide. Surprisingly, similarities in the signatures of different domains do not correlate with whole-domain sequence identities unless the latter is greater than 50%. An extensive comparison with the optimal binding motifs determined by peptide library experiments, as well as other experimental data indicate that the similarity in binding preferences of different SH2 domains can be deduced on the basis of their binding-site signatures. The analysis provides a rationale for the empirically derived classification of SH2 domains described by Songyang & Cantley, in that proteins in the same group are found to have similar residues at positions important for binding. Confident predictions of binding preference can be made for about 85% of SH2 domain sequences found in SWISSPROT. The approach described in this work is quite general and can, in principle, be used to analyze binding preferences of members of large protein families for which structural information for a number of family members is available. It also offers a strategy for predicting cross-reactivity of compounds designed to bind to a particular target, for example in structure-based drug design.

Specificity in Signal Transduction

Over the last two decades, a new and unifying concept of cellular organization has emerged in which modular protein-protein interactions provide an underlying framework through which signaling pathways are assembled and controlled. In this scheme, posttranslational modifications such as phosphorylation commonly exert their biological effects by regulating molecular interactions, exemplified by the ability of phosphotyrosine sites to bind selectively to SH2 domains. Although these interactions are rather simple in isolation, they can nonetheless be exploited to generate complex cellular systems. Here, I discuss experiments that have led to this view of dynamic cellular behavior and identify some current and future areas of interest in cell signaling.