Correlation analysis for protein evolutionary family based on amino acid position mutations and application in PDZ domain

Allostery Frustrates the Experimentalist

Proteins interact with other proteins, with nucleic acids, lipids, carbohydrates and various small molecules in the living cell. These interactions have been quantified and structurally characterized in numerous studies such that we today have a comprehensive picture of protein structure and function. However, proteins are dynamic and even folded proteins are likely more heterogeneous than they appear in most descriptions. One property of proteins that relies on dynamics and heterogeneity is allostery, the ability of a protein to change structure and function upon ligand binding to an allosteric site. Over the last decades the concept of allostery was broadened to embrace all types of long-range interactions across a protein including purely entropic changes without a conformational change in single protein domains. But with this re-definition came a problem: How do we measure allostery? In this opinion, we discuss some caveats arising from the quantitative description of single-domain allostery from an experimental perspective and how the limitations cannot be separated from the definition of allostery per se. Furthermore, we attempt to tie together allostery with the concept of frustration in an effort to investigate the links between these two complex, and yet general, properties of proteins. We arrive at the conclusion that the sensitivity to perturbation of allosteric networks in single protein domains is too large for the networks to be of significant biological relevance.

Investigating the allosteric response of the PICK1 PDZ domain to different ligands with all-atom simulations

The PDZ family is comprised of small modular domains that play critical roles in the allosteric modulation of many cellular signaling processes by binding to the C-terminal tail of different proteins. As dominant modular proteins that interact with a diverse set of peptides, it is of particular interest to explore how different binding partners induce different allosteric effects on the same PDZ domain. Because the PICK1 PDZ domain can bind different types of ligands, it is an ideal test case to answer this question and explore the network of interactions that give rise to dynamic allostery. Here, we use all-atom molecular dynamics simulations to explore dynamic allostery in the PICK1 PDZ domain by modeling two PICK1 PDZ systems: PICK1 PDZ-DAT and PICK1 PDZ-GluR2. Our results suggest that ligand binding to the PICK1 PDZ domain induces dynamic allostery at the αA helix that is similar to what has been observed in other PDZ domains. We found that the PICK1 PDZ-ligand distance is directly correlated with both dynamic changes of the αA helix and the distance between the αA helix and βB strand. Furthermore, our work identifies a hydrophobic core between DAT/GluR2 and I35 as a key interaction in inducing such dynamic allostery. Finally, the unique interaction patterns between different binding partners and the PICK1 PDZ domain can induce unique dynamic changes to the PICK1 PDZ domain. We suspect that unique allosteric coupling patterns with different ligands may play a critical role in how PICK1 performs its biological functions in various signaling networks.

Three Binding Conformations of BIO124 in the Pocket of the PICK1 PDZ Domain

The PDZ family has drawn attention as possible drug targets because of the domains’ wide ranges of function and highly conserved binding pockets. The PICK1 PDZ domain has been proposed as a possible drug target because the interactions between the PICK1 PDZ domain and the GluA2 subunit of the AMPA receptor have been shown to progress neurodegenerative diseases. BIO124 has been identified as a sub µM inhibitor of the PICK1–GluA2 interaction. Here, we use all-atom molecular dynamics simulations to reveal the atomic-level interaction pattern between the PICK1 PDZ domain and BIO124. Our simulations reveal three unique binding conformations of BIO124 in the PICK1 PDZ binding pocket, referred to here as state 0, state 1, and state 2. Each conformation is defined by a unique hydrogen bonding network and a unique pattern of hydrophobic interactions between BIO124 and the PICK1 PDZ domain. Interestingly, each conformation of BIO124 results in different dynamic changes to the PICK1 PDZ domain. Unlike states 1 and 2, state 0 induces dynamic coupling between BIO124 and the αA helix. Notably, this dynamic coupling with the αA helix is similar to what has been observed in other PDZ–ligand complexes. Our analysis indicates that the interactions formed between BIO124 and I35 may be the key to inducing dynamic coupling with the αA helix. Lastly, we suspect that the conformational shifts observed in our simulations may affect the stability and thus the overall effectiveness of BIO124. We propose that a physically larger inhibitor may be necessary to ensure sufficient interactions that permit stable binding between a drug and the PICK1 PDZ domain.

Allosterism in the PDZ Family

Dynamic allosterism allows the propagation of signal throughout a protein. The PDZ (PSD-95/Dlg1/ZO-1) family has been named as a classic example of dynamic allostery in small modular domains. While the PDZ family consists of more than 200 domains, previous efforts have primarily focused on a few well-studied PDZ domains, including PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ. Taken together, experimental and computational studies have identified regions of these domains that are dynamically coupled to ligand binding. These regions include the αA helix, the αB lower-loop, and the αC helix. In this review, we summarize the specific residues on the αA helix, the αB lower-loop, and the αC helix of PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ that have been identified as participants in dynamic allostery by either experimental or computational approaches. This review can serve as an index for researchers to look back on the previously identified allostery in the PDZ family. Interestingly, our summary of previous work reveals clear consistencies between the domains. While the PDZ family has a low sequence identity, we show that some of the most consistently identified allosteric residues within PTP-BL PDZ2 and PSD-95 PDZ3 domains are evolutionarily conserved. These residues include A46/A347, V61/V362, and L66/L367 on PTP-BL PDZ2 and PSD-95 PDZ3, respectively. Finally, we expose a need for future work to explore dynamic allostery within (1) PDZ domains with multiple binding partners and (2) multidomain constructs containing a PDZ domain.

Ensemble-Based Analysis of the Dynamic Allostery in the PSD-95 PDZ3 Domain in Relation to the General Variability of PDZ Structures

PDZ domains are abundant interaction hubs found in a number of different proteins and they exhibit characteristic differences in their structure and ligand specificity. Their internal dynamics have been proposed to contribute to their biological activity via changes in conformational entropy upon ligand binding and allosteric modulation. Here we investigate dynamic structural ensembles of PDZ3 of the postsynaptic protein PSD-95, calculated based on previously published backbone and side-chain S2 order parameters. We show that there are distinct but interdependent structural rearrangements in PDZ3 upon ligand binding and the presence of the intramolecular allosteric modulator helix α3. We have also compared these rearrangements in PDZ1-2 of PSD-95 and the conformational diversity of an extended set of PDZ domains available in the PDB database. We conclude that although the opening-closing rearrangement, occurring upon ligand binding, is likely a general feature for all PDZ domains, the conformer redistribution upon ligand binding along this mode is domain-dependent. Our findings suggest that the structural and functional diversity of PDZ domains is accompanied by a diversity of internal motional modes and their interdependence.

Seeking allosteric networks in PDZ domains

Ever since Ranganathan and coworkers subjected the covariation of amino acid residues in the postsynaptic density-95/Discs large/Zonula occludens 1 (PDZ) domain family to a statistical correlation analysis, PDZ domains have represented a paradigmatic family to explore single domain protein allostery. Nevertheless, several theoretical and experimental studies in the past two decades have contributed contradicting results with regard to structural localization of the allosteric networks, or even questioned their actual existence in PDZ domains. In this review, we first describe theoretical and experimental approaches that were used to probe the energetic network(s) in PDZ domains. We then compare the proposed networks for two well-studied PDZ domains namely the third PDZ domain from PSD-95 and the second PDZ domain from PTP-BL. Our analysis highlights the contradiction between the different methods and calls for additional work to better understand these allosteric phenomena.

Bi-factor analysis based on noise-reduction (BIFANR): a new algorithm for detecting coevolving amino acid sites in proteins

Previous statistical analyses have shown that amino acid sites in a protein evolve in a correlated way instead of independently. Even though located distantly in the linear sequence, the coevolved amino acids could be spatially adjacent in the tertiary structure, and constitute specific protein sectors. Moreover, these protein sectors are independent of one another in structure, function, and even evolution. Thus, systematic studies on protein sectors inside a protein will contribute to the clarification of protein function. In this paper, we propose a new algorithm BIFANR (Bi-factor Analysis Based on Noise-reduction) for detecting protein sectors in amino acid sequences. After applying BIFANR on S1A family and PDZ family, we carried out internal correlation test, statistical independence test, evolutionary rate analysis, evolutionary independence analysis, and function analysis to assess the prediction. The results showed that the amino acids in certain predicted protein sector are closely correlated in structure, function, and evolution, while protein sectors are nearly statistically independent. The results also indicated that the protein sectors have distinct evolutionary directions. In addition, compared with other algorithms, BIFANR has higher accuracy and robustness under the influence of noise sites.

Structural position correlation analysis (SPCA) for protein family

Background

The proteins in a family, which perform the similar biological functions, may have very different amino acid composition, but they must share the similar 3D structures, and keep a stable central region. In the conservative structure region similar biological functions are performed by two or three catalytic residues with the collaboration of several functional residues at key positions. Communication signals are conducted in a position network, adjusting the biological functions in the protein family.

Methodology

A computational approach, namely structural position correlation analysis (SPCA), is developed to analyze the correlation relationship between structural segments (or positions). The basic hypothesis of SPCA is that in a protein family the structural conservation is more important than the sequence conservation, and the local structural changes may contain information of biology functional evolution. A standard protein P⁽⁰⁾ is defined in a protein family, which consists of the most-frequent amino acids and takes the average structure of the protein family. The foundational variables of SPCA is the structural position displacements between the standard protein P⁽⁰⁾ and individual proteins P_i of the family. The structural positions are organized as segments, which are the stable units in structural displacements of the protein family. The biological function differences of protein members are determined by the position structural displacements of individual protein P_i to the standard protein P⁽⁰⁾. Correlation analysis is used to analyze the communication network among segments.

Conclusions

The structural position correlation analysis (SPCA) is able to find the correlation relationship among the structural segments (or positions) in a protein family, which cannot be detected by the amino acid sequence and frequency-based methods. The functional communication network among the structural segments (or positions) in protein family, revealed by SPCA approach, well illustrate the distantly allosteric interactions, and contains valuable information for protein engineering study.