Unusual binding interactions in PDZ domain crystal structures help explain binding mechanisms

PDZ domains most commonly bind the C‐terminus of their protein targets. Typically the C‐terminal four residues of the protein target are considered as the binding motif, particularly the C‐terminal residue (P0) and third‐last residue (P‐2) that form the major contacts with the PDZ domain's “binding groove”. We solved crystal structures of seven human PDZ domains, including five of the seven PDLIM family members. The structures of GRASP, PDLIM2, PDLIM5, and PDLIM7 show a binding mode with only the C‐terminal P0 residue bound in the binding groove. Importantly, in some cases, the P‐2 residue formed interactions outside of the binding groove, providing insight into the influence of residues remote from the binding groove on selectivity. In the GRASP structure, we observed both canonical and noncanonical binding in the two molecules present in the asymmetric unit making a direct comparison of these binding modes possible. In addition, structures of the PDZ domains from PDLIM1 and PDLIM4 also presented here allow comparison with canonical binding for the PDLIM PDZ domain family. Although influenced by crystal packing arrangements, the structures nevertheless show that changes in the positions of PDZ domain side‐chains and the αB helix allow noncanonical binding interactions. These interactions may be indicative of intermediate states between unbound and fully bound PDZ domain and target protein. The noncanonical “perpendicular” binding observed potentially represents the general form of a kinetic intermediate. Comparison with canonical binding suggests that the rearrangement during binding involves both the PDZ domain and its ligand.     

Characterization of PDZ domain‐peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein–protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 217 peptides using a method called MIECSVM, which energetically characterizes the domain‐peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross‐validation and leave‐one‐domain‐out test showed that the MIEC‐SVM using all 44 PDZ‐peptide residue pairs at the interaction interface outperformed the sequence‐based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (～50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC‐SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ‐peptide recognition is established. The success of the MIEC‐SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein‐ peptide interactions and understanding protein recognition from a structural and energetic viewpoint.     

Characterization of PDZ domain-peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein-protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 2387 peptides using a method called MIEC-SVM, which energetically characterizes the domain-peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross-validation and leave-one-domain-out test showed that the MIEC-SVM using all 44 PDZ-peptide residue pairs at the interaction interface outperformed the sequence-based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (~50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC-SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ-peptide recognition is established. The success of the MIEC-SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein-peptide interactions and understanding protein recognition from a structural and energetic viewpoint.     

Causality,transfer entropy,and allosteric communication landscapes in proteins with harmonic interactions

A fast and approximate method of generating allosteric communication landscapes in proteins is presented by using Schreiber's entropy transfer concept in combination with the Gaussian Network Model of proteins. Predictions of the model and the allosteric communication landscapes generated show that information transfer in proteins does not necessarily take place along a single path, but an ensemble of pathways is possible. The model emphasizes that knowledge of entropy only is not sufficient for determining allosteric communication and additional information based on time delayed correlations should be introduced, which leads to the presence of causality in proteins. The model provides a simple tool for mapping entropy sink‐source relations into pairs of residues. By this approach, residues that should be manipulated to control protein activity may be determined. This should be of great importance for allosteric drug design and for understanding the effects of mutations on function. The model is applied to determine allosteric communication in three proteins, Ubiquitin, Pyruvate Kinase, and the PDZ domain. Predictions are in agreement with molecular dynamics simulations and experimental evidence. Proteins 2017; 85:1056–1064. © 2017 Wiley Periodicals, Inc.     

On evolutionary conservation of thermodynamic coupling in proteins

The inherent complexity of thermodynamic coupling in proteins presents a major challenge in understanding and engineering protein function. Recent work has argued that the study of proteins can be simplified by the use of correlated mutations in the evolutionary record to locate a small subset of thermodynamically coupled residues that participate in functionally important, evolutionarily conserved energetic pathways. To test this hypothesis, we examined the predictions of correlated mutation algorithms for a number of proteins for which coupling between residues has been determined by analysis of double mutant cycles. We find that correlated mutation algorithms can find residue pairs that are physically close and that physically close residue pairs tend to be thermodynamically coupled. We find little evidence, however, for the hypothesis that thermodynamic coupling is limited to the subset of evolutionarily constrained residue positions.     

Molecular determinants for the complex binding specificity of the PDZ domain in PICK1

PICK1 (protein interacting with C kinase 1) contains a single PDZ domain known to mediate interaction with the C termini of several receptors, transporters, ion channels, and kinases. In contrast to most PDZ domains, the PICK1 PDZ domain interacts with binding sequences classifiable as type I (terminating in (S/T)XPhi; X, any residue) as well as type II (PhiXPhi; Phi, any hydrophobic residue). To enable direct assessment of the affinity of the PICK1 PDZ domain for its binding partners we developed a purification scheme for PICK1 and a novel quantitative binding assay based on fluorescence polarization. Our results showed that the PICK1 PDZ domain binds the type II sequence presented by the human dopamine transporter (-WLKV) with an almost 15-fold and >100-fold higher affinity than the type I sequences presented by protein kinase Calpha (-QSAV) and the beta(2)-adrenergic receptor (-DSLL), respectively. Mutational analysis of Lys(83) in the alphaB1 position of the PDZ domain suggested that this residue mimics the function of hydrophobic residues present in this position in regular type II PDZ domains. The PICK1 PDZ domain was moreover found to prefer small hydrophobic residues in the C-terminal P(0) position of the ligand. Molecular modeling predicted a rank order of (Val > Ile > Leu) that was verified experimentally with up to a approximately 16-fold difference in binding affinity between a valine and a leucine in P(0). The results define the structural basis for the unusual binding pattern of the PICK1 PDZ domain by substantiating the critical role of the alphaB1 position (Lys(83)) and of discrete side chain differences in position P(0) of the ligands.     

Severity of osteogenesis imperfecta and structure of a collagen-like peptide modeling a lethal mutation site

We show that there are correlations between the severities of osteogenesis imperfecta (OI) phenotypes and changes in the residues near the mutation site. Our results show the correlations between the severity of various forms of the inherited disease OI and alteration of residues near the site of OI causing mutations. Among our many observed correlations are particularly striking ones between the presence of nearby proline residues and lethal mutations, and the presence of nearby alanines residues and nonlethal mutations. We investigated the possibility that these correlations have a structural basis using molecular dynamics simulations of collagen-like molecules designed to mimic the site of a lethal OI mutation in collagen type I. Our significant finding is that interchain hydrogen bonding is greatly affected by variations in residue type. We found that the strength of hydrogen bond networks between backbone atoms on different chains depends on the local residue sequence and is weaker in proline-rich regions of the molecule. We also found that an alanine at a site near an OI mutation causes less structural disruption than a proline, and that residue side chains also form interchain hydrogen bonds with frequencies that are dependent on residue type. For example, arginine side chains form strong hydrogen bonds with the backbone of the subsequent peptide chain, while lysine and glutamine less frequently form similar hydrogen bonds. This decrease in the observed hydrogen bond frequency correlates with a decrease in the experimentally determined thermal stability. We contrasted general structural properties of model collagen peptides with and without the mutation to examine the effect of the single-point mutation on the surrounding residues.     

Proteolytic activation of Chlamydia trachomatis HTRA is mediated by PDZ1 domain interactions with protease domain loops L3 and LC and beta strand β5

Chlamydia trachomatis is a bacterial pathogen responsible for one of the most prevalent sexually transmitted infections worldwide. Its unique development cycle has limited our understanding of its pathogenic mechanisms. However, CtHtrA has recently been identified as a potential C. trachomatis virulence factor. CtHtrA is a tightly regulated quality control protein with a monomeric structural unit comprised of a chymotrypsin-like protease domain and two PDZ domains. Activation of proteolytic activity relies on the C-terminus of the substrate allosterically binding to the PDZ1 domain, which triggers subsequent conformational change and oligomerization of the protein into 24-mers enabling proteolysis. This activation is mediated by a cascade of precise structural arrangements, but the specific CtHtrA residues and structural elements required to facilitate activation are unknown. Using in vitro analysis guided by homology modeling, we show that the mutation of residues Arg362 and Arg224, predicted to disrupt the interaction between the CtHtrA PDZ1 domain and loop L3, and between loop L3 and loop LD, respectively, are critical for the activation of proteolytic activity. We also demonstrate that mutation to residues Arg299 and Lys160, predicted to disrupt PDZ1 domain interactions with protease loop LC and strand β5, are also able to influence proteolysis, implying their involvement in the CtHtrA mechanism of activation. This is the first investigation of protease loop LC and strand β5 with respect to their potential interactions with the PDZ1 domain. Given their high level of conservation in bacterial HtrA, these structural elements may be equally significant in the activation mechanism of DegP and other HtrA family members.     

Change in allosteric network affects binding affinities of PDZ domains: analysis through perturbation response scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation.     

Highly covarying residues have a functional role in antibody constant domains

The ability to generate and design antibodies recognizing specific targets has revolutionized the pharmaceutical industry and medical imaging. Engineering antibody therapeutics in some cases requires modifying their constant domains to enable new and altered interactions. Engineering novel specificities into antibody constant domains has proved challenging due to the complexity of inter‐domain interactions. Covarying networks of residues that tend to cluster on the protein surface and near binding sites have been identified in some proteins. However, the underlying role these networks play in the protein resulting in their conservation remains unclear in most cases. Resolving their role is crucial, because residues in these networks are not viable design targets if their role is to maintain the fold of the protein. Conversely, these networks of residues are ideal candidates for manipulating specificity if they are primarily involved in binding, such as the myriad interdomain interactions maintained within antibodies. Here, we identify networks of evolutionarily‐related residues in C‐class antibody domains by evaluating covariation, a measure of propensity with which residue pairs vary dependently during evolution. We computationally test whether mutation of residues in these networks affects stability of the folded antibody domain, determining their viability as design candidates. We find that members of covarying networks cluster at domain‐domain interfaces, and that mutations to these residues are diverse and frequent during evolution, precluding their importance to domain stability. These results indicate that networks of covarying residues exist in antibody domains for functional reasons unrelated to thermodynamic stability, making them ideal targets for antibody design. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Synergy and anti-cooperativity in allostery: Molecular dynamics study of WT and oncogenic KRAS-RGL1

This study focuses on investigating the effects of an oncogenic mutation (G12V) on the stability and interactions within the KRAS-RGL1 protein complex. The KRAS-RGL1 complex is of particular interest due to its relevance to KRAS-associated cancers and the potential for developing targeted drugs against the KRAS system. The stability of the complex and the allosteric effects of specific residues are examined to understand their roles as modulators of complex stability and function. Using molecular dynamics simulations, we calculate the mutual information, MI, between two neighboring residues at the interface of the KRAS-RGL1 complex, and employ the concept of interaction information, II, to measure the contribution of a third residue to the interaction between interface residue pairs. Negative II indicates synergy, where the presence of the third residue strengthens the interaction, while positive II suggests anti-cooperativity. Our findings reveal that MI serves as a dominant factor in determining the results, with the G12V mutation increasing the MI between interface residues, indicating enhanced correlations due to the formation of a more compact structure in the complex. Interestingly, although II plays a role in understanding three-body interactions and the impact of distant residues, it is not significant enough to outweigh the influence of MI in determining the overall stability of the complex. Nevertheless, II may nonetheless be a relevant factor to consider in future drug design efforts. This study provides valuable insights into the mechanisms of complex stability and function, highlighting the significance of three-body interactions and the impact of distant residues on the binding stability of the complex. Additionally, our findings demonstrate that constraining the fluctuations of a third residue consistently increases the stability of the G12V variant, making it challenging to weaken complex formation of the mutated species through allosteric manipulation. The novel perspective offered by this approach on protein dynamics, function, and allostery has potential implications for understanding and targeting other protein complexes involved in vital cellular processes. The results contribute to our understanding of the effects of oncogenic mutations on protein–protein interactions and provide a foundation for future therapeutic interventions in the context of KRAS-associated cancers and beyond.     

The interaction of Jagged‐1 cytoplasmic tail with afadin PDZ domain is local,folding‐independent,and tuned by phosphorylation

Jagged‐1, one of the five Notch ligands in man, is a membrane‐spanning protein made of a large extracellular region and a 125‐residue cytoplasmic tail bearing a C‐terminal PDZ recognition motif (¹²¹³RMEYIV¹²¹⁸). Binding of Jagged‐1 intracellular region to the PDZ domain of afadin, a protein located at cell–cell adherens junctions, couples Notch signaling with the adhesion system and the cytoskeleton. Using NMR chemical shift perturbation and surface plasmon resonance, we studied the interaction between the PDZ domain of afadin (AF6_PDZ) and a series of polypeptides comprising the PDZ‐binding motif. Chemical shift mapping of AF6_PDZ upon binding of ligands of different length (6, 24, and 133 residues) showed that the interaction is strictly local and involves only the binding groove in the PDZ. The recombinant protein corresponding to the entire intracellular region of Jagged‐1, J1_ic, is mainly disordered in solution, and chemical shift mapping of J1_ic in the presence of AF6_PDZ showed that binding is not coupled to folding. Binding studies on a series of 24‐residue peptides phosphorylated at different positions showed that phosphorylation of the tyrosine at position ‐2 of the PDZ‐binding motif decreases its affinity for AF6_PDZ, and may play a role in the modulation of this interaction. Finally, we show that the R1213Q mutation located in the PDZ‐binding motif and associated with extrahepatic biliary atresia increases the affinity for AF6_PDZ, suggesting that this syndrome may arise from an imbalance in the coupling of Notch signaling to the cytoskeleton. Copyright © 2010 John Wiley & Sons, Ltd.     

Redox-regulated affinity of the third PDZ domain in the phosphotyrosine phosphatase PTP-BL for cysteine-containing target peptides

PDZ domains are protein-protein interaction modules that are crucial for the assembly of structural and signalling complexes. They specifically bind to short C-terminal peptides and occasionally to internal sequences that structurally resemble such peptide termini. The binding of PDZ domains is dominated by the residues at the P(0) and P(-2) position within these C-terminal targets, but other residues are also important in determining specificity. In this study, we analysed the binding specificity of the third PDZ domain of protein tyrosine phosphatase BAS-like (PTP-BL) using a C-terminal combinatorial peptide phage library. Binding of PDZ3 to C-termini is preferentially governed by two cysteine residues at the P(-1) and P(-4) position and a valine residue at the P(0) position. Interestingly, we found that this binding is lost upon addition of the reducing agent dithiothrietol, indicating that the interaction is disulfide-bridge-dependent. Site-directed mutagenesis of the single cysteine residue in PDZ3 revealed that this bridge formation does not occur intermolecularly, between peptide and PDZ3 domain, but rather is intramolecular. These data point to a preference of PTP-BL PDZ3 for cyclic C-terminal targets, which may suggest a redox state-sensing role at the cell cortex.     

The Cancer Mutation D83V Induces an α-Helix to β-Strand Conformation Switch in MEF2B

MEF2B is a major target of somatic mutations in non-Hodgkin lymphoma. Most of these mutations are non-synonymous substitutions of surface residues in the MADS-box/MEF2 domain. Among them, D83V is the most frequent mutation found in tumor cells. The link between this hotspot mutation and cancer is not well understood. Here we show that the D83V mutation induces a dramatic α-helix to β-strand switch in the MEF2 domain. Located in an α-helix region rich in β-branched residues, the D83V mutation not only removes the extensive helix stabilization interactions but also introduces an additional β-branched residue that further shifts the conformation equilibrium from α-helix to β-strand. Cross-database analyses of cancer mutations and chameleon sequences revealed a number of well-known cancer targets harboring β-strand favoring mutations in chameleon α-helices, suggesting a commonality of such conformational switch in certain cancers and a new factor to consider when stratifying the rapidly expanding cancer mutation data.     

Canonical and Noncanonical Sites Determine NPT2A Binding Selectivity to NHERF1 PDZ1

Na⁺/H⁺ Exchanger Regulatory Factor-1 (NHERF1) is a scaffolding protein containing 2 PDZ domains that coordinates the assembly and trafficking of transmembrane receptors and ion channels. Most target proteins harboring a C-terminus recognition motif bind more-or-less equivalently to the either PDZ domain, which contain identical core-binding motifs. However some substrates such as the type II sodium-dependent phosphate co-transporter (NPT2A), uniquely bind only one PDZ domain. We sought to define the structural determinants responsible for the specificity of interaction between NHERF1 PDZ domains and NPT2A. By performing all-atom/explicit-solvent molecular dynamics (MD) simulations in combination with biological mutagenesis, fluorescent polarization (FP) binding assays, and isothermal titration calorimetry (ITC), we found that in addition to canonical interactions of residues at 0 and -2 positions, Arg at the -1 position of NPT2A plays a critical role in association with Glu43 and His27 of PDZ1 that are absent in PDZ2. Experimentally introduced mutation in PDZ1 (Glu43Asp and His27Asn) decreased binding to NPT2A. Conversely, introduction of Asp183Glu and Asn167His mutations in PDZ2 promoted the formation of favorable interactions yielding micromolar K _Ds. The results describe novel determinants within both the PDZ domain and outside the canonical PDZ-recognition motif that are responsible for discrimination of NPT2A between two PDZ domains. The results challenge general paradigms for PDZ recognition and suggest new targets for drug development.     

Specific myosin heavy chain mutations suppress troponin I defects in Drosophila muscles

We show that specific mutations in the head of the thick filament molecule myosin heavy chain prevent a degenerative muscle syndrome resulting from the hdp2 mutation in the thin filament protein troponin I. One mutation deletes eight residues from the actin binding loop of myosin, while a second affects a residue at the base of this loop. Two other mutations affect amino acids near the site of nucleotide entry and exit in the motor domain. We document the degree of phenotypic rescue each suppressor permits and show that other point mutations in myosin, as well as null mutations, fail to suppress the hdp2 phenotype. We discuss mechanisms by which the hdp2 phenotypes are suppressed and conclude that the specific residues we identified in myosin are important in regulating thick and thin filament interactions. This in vivo approach to dissecting the contractile cycle defines novel molecular processes that may be difficult to uncover by biochemical and structural analysis. Our study illustrates how expression of genetic defects are dependent upon genetic background, and therefore could have implications for understanding gene interactions in human disease.     

Comparison of the PKCalpha and the PKCepsilon C1b domains: identification of residues critical for PKCepsilon-mediated neurite induction

We showed earlier that over-expression of protein kinase C (PKC) epsilon induces neurite outgrowth. The effect is mediated by a region (PKCepsilonPSC1V3) encompassing the pseudosubstrate, the two C1 domains and part of the V3 region, and is independent of the catalytic activity of the enzyme. In this region, residues immediately N-terminal of the C1b domain are crucial for neurite outgrowth. However, in this study we show that the PKCepsilon C1b domain itself is necessary for neurite induction, since a mutant in which the PKCepsilon C1b domain has been replaced with the C1b domain from PKCalpha, PKCepsilonPSC1a(alphaC1b)V3 lacks neurite-inducing capacity. The molecular basis for the importance of the PKCepsilon C1b domain was investigated by mutation studies of the PKCalpha C1b domain. Point mutations were done in the PKCalpha C1b domain of the PKCepsilonPSC1a(alphaC1b)V3 construct, in which the PKCalpha residues were mutated into the corresponding residues in PKCepsilon. This highlighted residues in the C-terminal part of the primary sequence of the C1b domain, located in the base of the C1b domain, as important for neurite outgrowth. The mutations S48P, D32K and L49N all influenced neurite induction positively. Furthermore, the mutation of L49N alone was sufficient to make PKCepsilonPSC1a(alphaC1b)V3 neuritogenic in phorbol ester-stimulated cells, and mutation of this residue in full-length PKCepsilon into the corresponding residue in PKCalpha, N291L reduced the neurite-inducing effect of PKCepsilon. In conclusion, we have identified residues in the PKCepsilon C1b domain, in particular Asn49, that are essential for neurite outgrowth.     

Probing subtle conformational changes induced by phosphorylation and point mutations in the TIR domains of TLR2 and TLR3

Extensive research performed on Toll‐like receptor (TLR) signaling has identified residues in the Toll/interleukin‐1 receptor (TIR) domains that are essential for its proper functioning. Among these residues, those in BB loop are particularly significant as single amino acid mutations in this region can cause drastic changes in downstream signaling. However, while the effect of these mutations on the function is well studied (like the P681H mutation in TLR2, the A795P mutation in TLR3, and the P714H mutation in TLR4), their influence on the dynamics and inter‐residue networks is not well understood. The effects of local perturbations induced by these mutations could propagate throughout the TIR domain, influencing interactions with other TIR domain‐containing proteins. The identification of these subtle changes in inter‐residue interactions can provide new insights and structural rationale for how single‐point mutations cause drastic changes in TIR–TIR interactions. We employed molecular dynamics simulations and protein structure network (PSN) analyses to investigate the structural transitions with special emphasis on TLR2 and TLR3. Our results reveal that phosphorylation of the Tyr 759 residue in the TIR domain of TLR3 introduces rigidity to its BB loop. Subtle differences in the intra BB loop hydrogen bonding network between TLR3 and TLR2 are also observed. The PSN analyses indicate that the TIR domain is highly connected and pinpoints key differences in the inter‐residue interactions between the wild‐type and mutant TIR domains, suggesting that TIR domain structure is prone to allosteric effects, consistent with the current view of the influence of allostery on TLR signaling.     

Signaling pathways of PDZ2 domain: a molecular dynamics interaction correlation analysis

PDZ domains are found in many signaling proteins. One of their functions is to provide scaffolds for forming membrane-associated protein complexes by binding to the carboxyl termini of their partners. PDZ domains are thought also to play a signal transduction role by propagating the information that binding has occurred to remote sites. In this study, a molecular dynamics (MD) simulation-based approach, referred to as an interaction correlation analysis, is applied to the PDZ2 domain to identify the possible signal transduction pathways. A residue correlation matrix is constructed from the interaction energy correlations between all residue pairs obtained from the MD simulations. Two continuous interaction pathways, starting at the ligand binding pocket, are identified by a hierarchical clustering analysis of the residue correlation matrix. One pathway is mainly localized at the N-terminal side of helix alpha1 and the adjacent C-terminus of loop beta1-beta2. The other pathway is perpendicular to the central beta-sheet and extends toward the side of PDZ2 domain opposite to the ligand binding pocket. The results complement previous studies based on multiple sequence analysis, NMR, and MD simulations. Importantly, they reveal the energetic origin of the long-range coupling. The PDZ2 results, as well as the earlier rhodopsin analysis, show that the interaction correlation analysis is a robust approach for determining pathways of intramolecular signal transduction.     

Critical contacts between the eukaryotic initiation factor 2B (eIF2B) catalytic domain and both eIF2beta and -2gamma mediate guanine nucleotide exchange

Diverse guanine nucleotide exchange factors (GEFs) regulate the activity of GTP binding proteins. One of the most complicated pairs is eukaryotic initiation factor 2B (eIF2B) and eIF2, which function during protein synthesis initiation in eukaryotes. We have mutated conserved surface residues within the eIF2B GEF domain, located at the eIF2Bepsilon C terminus. Extensive genetic and biochemical characterization established how these residues contribute to GEF activity. We find that the universally conserved residue E569 is critical for activity and that even a conservative E569D substitution is lethal in vivo. Several mutations within residues close to E569 have no discernible effect on growth or GCN4 expression, but an alanine substitution at the adjacent L568 is cold sensitive and deregulates GCN4 activity at 15 degrees C. The mutation of W699, found on a separate surface approximately 40 A from E569, is also lethal. Binding studies show that W699 is critical for interaction with eIF2beta, while L568 and E569 are not. In contrast, all three residues are critical for interaction with eIF2gamma. These data show that multiple contacts between eIF2gamma and eIF2Bepsilon mediate nucleotide exchange.