- Research article
- Open Access
Mechanism of PhosphoThreonine/Serine Recognition and Specificity for Modular Domains from All-atom Molecular Dynamics
© Huang and Chang; licensee BioMed Central Ltd. 2011
Received: 24 March 2011
Accepted: 25 May 2011
Published: 25 May 2011
Phosphopeptide-binding domains mediate many vital cellular processes such as signal transduction and protein recognition. We studied three well-known domains important for signal transduction: BRCT repeats, WW domain and forkhead-associated (FHA) domain. The first two recognize both phosphothreonine (pThr) and phosphoserine (pSer) residues, but FHA has high specificity for pThr residues. Here we used molecular dynamics (MD) simulations to reveal how FHA exclusively chooses pThr and how BRCT and WW recognize both pThr/pSer. The work also investigated the energies and thermodynamic information of intermolecular interactions.
Simulations carried out included wide-type and mutated systems. Through analysis of MD simulations, we found that the conserved His residue defines dual loops feature of the FHA domain, which creates a small cavity reserved for only the methyl group of pThr. These well-organized loop interactions directly response to the pThr binding selectivity, while single loop (the 2nd phosphobinding site of FHA) or in combination with α-helix (BRCT repeats) or β-sheet (WW domain) fail to differentiate pThr/pSer.
Understanding the domain pre-organizations constructed by conserved residues and the driving force of domain-phosphopeptide recognition provides structural insight into pThr specific binding, which also helps in engineering proteins and designing peptide inhibitors.
Protein phosphorylation is widely exploited in DNA damage repair, signal transduction, cell growth and cell cycle regulation; the cascades of downstream signals can be triggered by grabbing a certain phosphoprotein [1–6]. Elucidating the characteristics of phosphopeptide recognition is fundamental to study cellular functions . The phosphoproteins are usually classified into two families, phosphotyrosine (pTyr)-containing and phosphoserine (pSer)/phosphothreonine (pThr)-containing sequences which are phosphorylated and dephosphorylated by different categories of kinases (e.g., pThr/pSer kinase and pThr kinase) and phosphatases . Recent studies discovered a few modular domains that particularly recognize pThr/pSer- or pThr-containing sequences, such as the breast-cancer-associated protein BRCA1 C-terminal (BRCT) repeats, WW domain and forkhead-associated (FHA) domain . Among them, the FHA domain differentiates pThr-containing peptides from pSer-containing peptides, although the difference is only one methyl group [10–15]. Because Ser/Thr kinase phosphorylates both residues, the FHA domain can efficiently reduce potential interaction sites by specifically binding to pThr-containing regions. Although experimental structures and recent studies have shown important interactions involve in the binding of the methyl group of pThr [14, 16], the detailed mechanisms of the phosphoresidue recognition of different domains and how the FHA domain can reserve non-polar interactions for a small non-polar methyl group are not fully understood.
The FHA domain is associated with proteins of diverse functions in different organisms. For example, the Rad53-FHA1 domain interacts with phosphorylated Rad9 in response to DNA damage, and the Dun1-FHA domain interacts with SCD1 of Rad53, which leads to activation of Dun1in response to DNA damage response . Instead of binding to a single pThr, the formation of the Dun1-FHA:SCD1 complex requires two phosphoresidues, which suggests a potential mechanism whereby sequential signaling events could be triggered through the recognition of multiple phosphoresidue-binding sites. Similar to Dun1-FHA, Ki67-FHA domain also recognizes dual phosphorylated residues at the same time , and the sequence identities of Dun1-FHA and the Ki67-FHA to Rad53-FHA1 are both 34%.
BRCT repeats in BRCA1 are considered to be related to breast cancer [24–27]. The structure of the tandem BRCA1-BRCT repeats bind to phosphorylated protein that contains pSer or pThr, although binding to pSer is preferred . Several structural studies have revealed a conserved structure for the repeats, mainly composed of α helixes, β sheets and loops that link secondary structures. The phosphoresidue-recognized site is located between β1 and α2 (Figure 1 (b)) [29–31].
The WW domain in Pin1 is essential for mitotic progression . The domain has only about 40 residues and is one of the smallest pThr/pSer binding domains . It specifically binds to pThr-Pro- or pSer-Pro-containing motifs with slightly higher affinity for pThr-Pro-containing peptides . For example, in the Pin1-WW domain, the aromatic rings of Tyr23 and Trp34 define a steric clamp to confer a Pro adjacent to pSer/pThr . The WW domain folds into three anti-parallel β stands and Arg21 and Ser22 residues in the loops between β1 and β2 are the phosphate group recognition sites (Figure 1 (c)) [3, 4, 34].
Because FHA domains are pThr specific modular domains, this study focused on how FHA domains display selective for pThr/pSer residues and comparison with BRCT and WW domains. We study the dynamic and conformational changes of the free domain and the complexes of Rad53-FHA1:Rad9, Dun1-FHA:SCD1, Ki67-FHA:hNIFK, BRCA1-BRCT:BRCH1 and Pin1-WW:CTD systems. We also computed the interaction energy between pThr/pSer and the domains to disclose the driving force of pThr/pSer binding. We propose a model for pThr specificity and potential applications.
Peptide sequences of domain-phosphopeptide complexes
Molecular dynamics simulation
We performed molecular dynamics (MD) simulations using the Amber10 and NAMD2.6 simulation packages with the ff99sb amber force field [36–38]. Because phosphoresidues are not defined in the ff99sb force field, we used the pThr and pSer force field developed by Homeyer et al. . All simulations of wild-type sequences were initialed from six experimental complexes. In this work, we studied six un-mutated complexes and eight mutated structures. The protonation states were checked by the MCCE program . All complexes were solvated in a rectangular box of 12 Å TIP3P water with the tleap program in Amber10 . The placement of counter-ions of Na+ was based on the Columbic potential to keep the whole system neutral. Particle Mesh Ewald was used to consider the long-range electrostatic interactions . Following 10,000 and 20,000 steps of minimization of the water and system, respectively, the systems were gradually heated for each complex from 250 K for 20-ps, 275 K for 20-ps and 300 K for 200-ps. To initial the mutated structures, after equilibrium from 300 K, we manually added or deleted the methyl group and changed the residue name accordingly, then used the Amber program to build the prmtop files for mutants. A quick 100-steps minimization was applied to the mutants, then we preformed 20-ps equilibrium at 300 K. All MD simulations for each wild-type and mutated complex was performed in 1 ns by five different random number seeds to generate different initial velocity. The resulting trajectories were collected every 1 ps and the time step was 2-fs. The NPT ensemble was applied, and periodic boundary conditions were used throughout the MD simulations. A temperature of 300 K was maintained with use of a Langevin thermostat, with a damping constant of 2 ps-1, and the hybrid Nose-Hoover Langevin piston method was used to control the pressure at 1 atm. The SHAKE procedure was used to constrain hydrogen atoms during MD simulations .
Binding energy calculation by MM-PBSA methods and entropy calculation
representing the changes in valance (v) energy (bond, angle, dihedral, and improper dihedral energies), van der Waal (vdw) interactions, Coulombic (Coul) interactions, and solvation free energy ΔWPB/ΔWGB and ΔWnp. We note that the binding energy computed here includes the solvation free energy which considers water entropy, and the valance energies cancel out in Eq. 2 due to the single trajectory approach. The solvation free energy can be further decomposed into the polar term, ΔWPB/ΔWGB, and non-polar cavity term, ΔWnp term [44, 51]. Here, we demonstrate two methods, Poisson-Boltzmann (PB) and Generalized-Born (GB), to estimate the polar solvation term [45, 46]. PB was calculated by solving the PB equation in the PBSA model of Amber11. The dielectric constants of the interior and exterior protein were set to 1 and 80 respectively. GB (igb = 1) was used in the sander program of Amber11 package. The non-polar solvation term was calculated by the solvent-accessible surface area (SASA) model. All energy involved a 40 Å cutoff value for non-bonded interactions.
where X denotes each dihedral angle, such as phi, psi, omega and sidechain.
Results and discussion
Ser and Thr have very similar sidechains, so the Ser/Thr protein kinases phosphorylate the OH group of either residue without discriminating between them. Similarly, most pSer/pThr binding modules, such as BRCT and WW domains, can specifically bind to short pSer/pThr-containing motifs equally well. Not unsurprisingly, one more methyl group of Thr has few effects on the overall binding and molecular recognition. However, uniquely, most FHA domains recognize only pThr residues in target proteins, and substitution of pSer for pThr in model peptides severely weakens binding. The pThr specific recognition may function as a filter to further select the protein partner.
One of our aims was to investigate how FHA domains can distinguish the tiny difference between Ser and Thr in their binding partner. The work helps gain deeper understanding in molecular recognition and provides valuable insight into strategies of protein engineering. Although we have simulated the entire FHA, WW and BRCT systems, we do not attempt to draw conclusions on properties of the whole system. Instead, we focus on the phosphoresidue binding region, in particular pThr or pSer. Multiple 1 ns MD simulations were performed for each system with different initial velocity so that the simulations evolve independently of each other to reduce potential bias caused by the initial conformation. Because the short phosphopeptides are highly flexible, carrying out short simulations allows the peptide to stay in similar local energy well during different runs [54, 55]. The root mean square deviation (RMSD) of selected systems is shown in Figure S1 of Additional file 1.
To obtain an estimate of the differences in interaction energies of the domain with its pSer- or pThr-containing peptides, we post-processed our MD trajectories and computed the domain-peptide interaction energies, including the potential energy and solvation energy. The method is usually called the MM-PBSA method. We computed intermolecular interactions between the whole domain and the entire peptide, termed the "global binding energy calculation". In our global energy calculations, the solvation energy term includes a PB term, WPB, for electrostatic solvation free energy, and a cavity/surface area term, Wnp, for nonpolar solvation free energy. Because we are particularly interested in studying the local region that contributes significantly to pThr or pSer recognition, we also selected residues within 5 to 7 Å of the phosphate group (Table S3 in Additional file 1) to calculate interaction energies between the residues chosen. We called the calculations between these selected residues "local interaction energy calculation". Because the calculations involved only residues near the binding sites, the solvation energies computed by the PB or GB model are similar. As a result, we used the GB model in our local energy calculation to speed up the calculation.
Revealing the specificity of the main pThr-binding site in FHA domains
Global MM-PBSA energy calculations
Local MM-PBSA energy calculations
Of note, although experimental structures demonstrate a pocket to accept the pThr methyl group, the static conformation cannot ascertain that pSer fails to form equally good interactions with the nearby residues, because the protein is dynamic and may fill the space by slightly changing the protein conformations. Nevertheless, our MD simulations show that the cavity is highly suited to pThr, and small changes in this particular residue can diminish the domain-peptide interactions.
Configuration entropy changes
The second phosphoresidue-binding site of Dun1-FHA and Ki67-FHA
Some FHA domains also show the second phosphoresidue-binding site, and knowing whether the second site can discriminate pThr and pSer is of interest. We therefore studied two diphosphoresidue-recognized FHA domains, Dun1-FHA and Ki67-FHA. Both domains have one pThr binding to the main pThr-binding site, but they also have one more phosphoresidue, pSer or pThr, in the peptide sequences. One main difference between the main pThr and the second phosphoresidue-binding site is that the main pThr-binding site consists of two loops that form a well-defined pocket, whereas the second binding site is located in areas with a single loop. To understand whether the difference contributes to residue specificity, we mutated the phosphoresidue available in the experiment, pThr of Dun1-FHA and pSer of Ki67-FHA to pSer and pThr, respectively.
Overall, the global binding energy calculations show that for the second phosphoresidue-binding site, the mutations worsen binding affinities (see Table 2), but the changes are smaller than the values for the main pThr-binding site. However, the local interaction energy calculations do not show the same trend, and the mutation of Dun1-FHA is preferable. Therefore, the calculations do not directly support that the domain strongly prefers either pThr or pSer in the second phosphoresidue-binding site. The local interaction energy calculations suggest that pSer can have good interactions with the domain, which are contributed mainly from the electrostatic attractions, and by losing the methyl group, the van der Waal interactions are weakened, but not significantly. As illustrated in Figure 2 (B2) and (B3), the second phosphoresidue-binding site in Dun1-FHA uses two Arg residues, Arg62 and Arg64, to recognize pThr or pSer and form multiple H-bonds with the phosphate group. Therefore, the electrostatic attractions are the major driving forces in pThr/pSer binding is not surprising . In addition, both Arg residues are located in one single loop, which is a flexible region of FHA domains, so the protein is freely adjustable to adopt both pThr and pSer. Although the methyl group of pThr forms non-polar attractions with the alkane branch of Arg62 shown in Figure 2 (B2), the binding site does not hold a small pocket when pThr is substituted by pSer (Figure 2 (B3)), because the space is filled by the nearby FHA domain sidechains.
The second phosphoresidue-binding site is located in a single loop, β3-β4, and without spatial constraint, the second site allows the FHA domain to rearrange sidechains to optimize both pThr and pSer binding. We therefore examined changes in local flexibility when a different phosphoresidue stays in this binding site, and we focused on dihedral angles of the phosphate group of pThr and pSer. Both residues do not expressly reveal dynamic motions in either complex, but the most populated angles modeled from our MD simulations shift more than 20° (Figure 3 (C1) (C2)). Although the dihedral angle has only one rotameric states in both cases, pThr has smaller vibration range and the configuration entropy is 0.5-1.0 kcal/mol smaller than pSer, presumably due to a bulkier methyl group. The entropy changes between pThr and pSer is less pronounced in the second phosphoresidue-binding site than in the main one. Moreover, the motions of backbone dihedrals remain the same, which indicates the negligible influence of the mutation.
In conclusion, the main pThr-binding site has a unique feature to recognize pThr, and a special pocket built by linking two loops with the conserved His is reserved for the methyl group of pThr, which plays a crucial role in distinguishing between pThr and pSer. However, the second phosphoresidue binding site is positioned in a single loop near the N-terminus, which uses two Arg residues to recognize a phosphate group but lacks a well structured binding cavity to identify only pThr or pSer. Our simulations show that the protein sidechain of the binding site changes when pThr or pSer binds to the domain. Therefore, a single loop used to provide a phosphoresidue-binding site cannot discriminate pThr/pSer but can bind to both residues. Other domains, such as the WW domain, also use a similar strategy, as discussed in the next section. Of note, although not within the scope of this paper, the promiscuous domain has preferences for selected sequences, and the peptide sequences also play an important role in the phosphopeptide-binding site.
BRCT repeats and WW domain recognize both pThr and pSer
The BRCT and WW domains are distinct pSer/pThr binding domains. Not all BRCT and WW domains function as phosphopeptide binding modules, but both have a subset that binds to phosphopeptides. Both domains can bind to particular sequences that contain pSer or pThr, but in general, tandem BRCT domains bind stronger to pSer than pThr and WW domains have a preference for pThr preceding a Pro [28, 33]. Note that pSer or pThr must be followed by Pro, for pSer/pThr-Pro sequences for binding to WW domains. In contrast to FHA domains, which bind exclusively to pThr-containing peptides, BRCT and WW domains do not recognize solely pThr- or pSer-containing peptides. Although proteins that treat Ser and Thr as similar residues may be common, knowing how both domains have a specific or non-specific pThr/pSer recognition is of interest.
Both global binding energy and local interaction energy calculations suggest that the tandem BRCT domains prefer the pSer- than pThr-containing sequence, although the preference is not strong. The local interaction energies shown in Table 3 suggest that pSer can form a more favorable electrostatic attraction, ~0.5 kcal/mol more negative than that contributed by pThr binding, but the difference is relatively small. In addition, our simulations show that most of the time, the methyl group of pThr does not directly interact with the domain, as demonstrated by a representative complex conformation in Figure 2 (C2). As a result, we see a negligible difference in van der Waal interactions in the local interaction energy calculations (see Table 3), and the electrostatic attractions are the main driving forces to recruit phosphopeptides binding to tandem BRCT domains (Table S2 in Additional file 1). The phospho recognition is through forming interactions with Lys1702 in an α helix and Ser1655 and Gly1656 near the loop, where no small cavity is reserved for the methyl group of pThr. Although the methyl group is not directly involved in binding, the local arrangement of the phosphate group is changed, but the overall flexibility of the phosphate group remains similar (Figure 2 (C)).
Although the WW domains are able to recognize both pSer- and pThr-containing peptides, global binding energy calculations suggest that the domain favors pThr because of the more preferable van der Waals attractions. The trend is in agreement with experimental results . However, the local interaction energy calculations show that the favorable van der Waals attractions are mostly compensated by weaker electrostatic interactions. As shown in Figure 2 (D2) and (D3), two conserved aromatic residues of the domain, Tyr and Trp, create a cavity, but no sidechains of phosphopeptides could nicely fit into the cavity during our simulations. Interestingly, the conserved Pro residue adjacent to the phosphoresidue is clamped by Tyr and Trp, which stays in the cavity and further restricts nearby phosphopeptide conformations. The confined region formed by rings of Tyr, Trp and Pro is conserved regardless of the presence of pSer or pThr (see Figure 2 (D2) and (D3)), which also explains the crucial roles of Pro. Because of the bulky ring conformations, an empty space is observed during the course of the simulations. The empty space can be partially filled by the methyl group of pThr, thus resulting in more favorable van der Waals interactions and a less flexible sidechain while pThr is binding. However, the Pro residue but not the methyl group of pThr primarily occupies the cavity in phosphopeptide recognition. Therefore, the domains do not show significant discrimination between pSer and pThr.
Comparisons between FHA domains and tandem BRCT repeats and WW domain
FHA domains use the conserved His to bridge two loops, β4-β5 and β6-β7, to construct the main pThr binding, which have a phosphate group binding site and preserve a small pocket nicely fit by the methyl group of pThr (see Figure 4(C)). Without the methyl group, neither pSer binding nor rearranging sidechains of FHA near the methyl binding site can effectively fill the pocket, which results in unsuccessful binding. The second binding site of FHA domains makes use of a single loop to recognize pThr/pSer, and WW and BRCT domains combine a single loop and a nearby α helix (BRCT domains) or β sheet (WW domains) to bind to pThr/pSer. The structures of these phosphoresidue binding sites allow the protein sidechains to be adjustable to both pThr and pSer residues. Notably, although loops are usually considered flexible regions of a protein, the dual loops in the main pThr-binding site of FHA domains show an interaction network between the loops to form a pre-organized binding cavity for pThr (data not shown). Besides the unique features of using dual loops specifically for pThr binding, all other phosphoresidue binding sites share common characteristics that include a binding site comprised of positive-charged residues to form Coulombic attractions with the phosphate group and geometry complementary in the binding surface.
Of note, the binding affinities of phosphopeptide binding to these domains are in general weak, in the micromolar range (see Table 1); therefore, weakened attractions by a few kcal/mol can completely diminish the phosphopeptides binding. Therefore, although substitution of pThr by pSer mainly reduces van der Waals attractions in the main pThr binding site of FHA domains, the pSer-containing peptide cannot form the domain-peptide complex. Different from other phosphopeptide binding sites, sidechain rearrangements cannot bring other attraction forces to compensate for the loss of the van der Waals interactions because of the rigid structure formed by His and the dual loops. In addition to energy calculations, our local entropy calculations suggest that binding pThr to the main pThr-binding site of FHA can reduce the mobility significantly, which indicates stronger attraction and more geometry confinement. However, the entropy changes between pThr and pSer binding to other domains show smaller differences, which suggests that the system retains a similar dynamic behavior that may help balance energy loss by gaining other attractions.
Modular domains are common regulators in important biological processes. This work studied three important domains for DNA damage responses, FHA, BRCT and WW domains, all with a phosphopeptide binding site to relay the damage signal and trigger further repair. The specific peptide can be recognized by three different proteins: a protein kinase to phosphorylate Ser/Thr, a modular domain that binds to the phosphopeptide for a downstream process, and a phosphatase to dephosphorylate the phosphoresidue. Kinases involved in the DNA damage response, ATM and ATR, can phosphorylate both Ser and Thr of a substrate. Similar to kinases, phosphatases work for both pThr and pSer , and most phosphodomains can also bind to both pThr- and pSer-containing peptides. FHA domains have evolved a simple but remarkable mechanism to specifically recognize pThr to further select particular partners after kinase phosphorylation. For example, a binding partner of FHA domains, the SCD protein family, contains rich Thr-glutamine (TQ) and Ser-glutamine (SQ) repeat motifs . Although the kinase phosphorylates both Thr and Ser, Rad53-FHA1 can bind only to regions that have pThr. How binding to the particular pThr region triggers further responses is unclear, but the pThr-only recognition may play a role in regulation. The mechanism contributing to pThr binding brings insights into how modular domains differentiate pThr and pSer or recognize both residues. The information aids in the design and discovery of phosphopeptides to access the cellular function of the domain-containing proteins. The rigid dual loops centered on the conserved His in the main pThr binding site of FHA domains may be applied to protein engineering that may need to recognize small functional groups.
In this study, we performed dynamic-guided process for FHA, BRCT and WW domain-peptide structures. The components of detailed interaction energies were calculated by MM-PBSA/MM-GBSA method. The main pThr-binding cavity is identical in four different FHA complexes. Our results reveal FHA domain uses the conserved His residue to define a dual loops structure which shows strong favor for pThr because of the geometry of methyl group embedded in deep binding pocket nicely. The dynamics simulations, energy and entropy calculations indicate that the phosphoresidue binding site of FHA is highly suited to pThr, and small changes of pThr to pSer can diminish the domain-peptide interactions due to the pre-organized binding cavity. On the other hand, BRCT repeats and WW domain utilize the combination a single loop with α-helix or β-sheet which allows effectively sidechain rearrangement to accept both pSer and pThr. The results highlight broader implications in recognition pathway of kinase/phosphotase and also help to engineer proteins and design peptide inhibitors.
We thank Dr. Ming-Daw Tsai for valuable suggestions and discussion. This research was supported in part by start-up funds from the University of California, Riverside and the National Science Foundation (MCB-0919586) through TeraGrid resources provided by the National Center for Supercomputing Application under grant number (TG-MCB080039N) and ShaRCS, the University of California Shared Research Computing Services Cluster, which is technically supported by multiple UC information technology divisions and managed by the University of California, Office of the President.
- Westheimer FH: WHY NATURE CHOSE PHOSPHATES. Science. 1987, 235 (4793): 1173-1178. 10.1126/science.2434996.View ArticleGoogle Scholar
- Yaffe MB, Cantley LC: Signal transduction - Grabbing phosphoproteins. Nature. 1999, 402 (6757): 30-31. 10.1038/46925.View ArticleGoogle Scholar
- Yaffe MB, Smerdon SJ: Phosphoserine/threonine binding domains: You can't pSERious?. Structure. 2001, 9 (3): R33-R38. 10.1016/S0969-2126(01)00580-9.View ArticleGoogle Scholar
- Yaffe MB, Smerdon SJ: The use of in vitro peptide-library screens in the analysis of phosphoserine/threonine-binding domain structure and function. Annual Review of Biophysics and Biomolecular Structure. 2004, 33: 225-244. 10.1146/annurev.biophys.33.110502.133346.View ArticleGoogle Scholar
- Seet BT, Dikic I, Zhou MM, Pawson T: Reading protein modifications with interaction domains. 2006, 7:Google Scholar
- Jorgensen C, Linding R: Directional and quantitative phosphorylation networks. Briefings in Functional Genomics & Proteomics. 2008, 7 (1): 17-26. 10.1093/bfgp/eln001.View ArticleGoogle Scholar
- Narayanan A, Jacobson MP: Computational studies of protein regulation by post-translational phosphorylation. Current Opinion in Structural Biology. 2009, 19 (2): 156-163. 10.1016/j.sbi.2009.02.007.View ArticleGoogle Scholar
- Virshup DM, Shenolikar S: From Promiscuity to Precision: Protein Phosphatases Get a Makeover. Molecular Cell. 2009, 33 (5): 537-545. 10.1016/j.molcel.2009.02.015.View ArticleGoogle Scholar
- Diella F, Haslam N, Chica C, Budd A, Michael S, Brown NP, Trave G, Gibson TJ: Understanding eukaryotic linear motifs and their role in cell signaling and regulation. Frontiers in Bioscience. 2008, 13: 6580-6603.View ArticleGoogle Scholar
- Liao H, Byeon IJL, Tsai MD: Structure and function of a new phosphopeptide-binding domain containing the FHA2 of Rad53. 1999, 294:Google Scholar
- Liao H, Yuan CH, Su MI, Yongkiettrakul S, Qin DY, Li HY, Byeon IJL, Pei DH, Tsai MD: Structure of the FHA1 domain of yeast Rad53 and identification of binding sites for both FHA1 and its target protein Rad9. Journal of Molecular Biology. 2000, 304 (5): 941-951. 10.1006/jmbi.2000.4291.View ArticleGoogle Scholar
- Durocher D, Taylor IA, Sarbassova D, Haire LF, Westcott SL, Jackson SP, Smerdon SJ, Yaffe MB: The molecular basis of FHA Domain: Phosphopeptide binding specificity and implications for phospho-dependent signaling mechanisms. 2000, 6:Google Scholar
- Liang XY, Van Doren SR: Mechanistic insights into phosphoprotein-binding FHA domains. 2008, 41:Google Scholar
- Mahajan A, Yuan CH, Lee H, Chen ESW, Wu PY, Tsai MD: Structure and Function of the Phosphothreonine-Specific FHA Domain. Science Signaling. 2009, 2 (58):Google Scholar
- Pennell S, Westcott S, Ortiz-Lombardia M, Patel D, Li JJ, Nott TJ, Mohammed D, Buxton RS, Yaffe MB, Verma C, Smerdon SJ: Structural and Functional Analysis of Phosphothreonine-Dependent FHA Domain Interactions. Structure. 2010, 18 (12): 1587-1595. 10.1016/j.str.2010.09.014.View ArticleGoogle Scholar
- Yuan CH, Yongkiettrakul S, Byeon IJL, Zhou SZ, Tsai MD: Solution structures of two FHA1-phosphothreonine peptide complexes provide insight into the structural basis of the ligand specificity of FHA1 from yeast Rad53. 2001, 314:Google Scholar
- Lee SJ, Schwartz MF, Duong JK, Stern DF: Rad53 phosphorylation site clusters are important for Rad53 regulation and signaling. Molecular and Cellular Biology. 2003, 23 (17): 6300-6314. 10.1128/MCB.23.17.6300-6314.2003.View ArticleGoogle Scholar
- Byeon IJL, Li HY, Song HY, Gronenborn AM, Tsai MD: Sequential phosphorylation and multisite interactions characterize specific target recognition by the FHA domain of Ki67. Nature Structural & Molecular Biology. 2005, 12 (11): 987-993.View ArticleGoogle Scholar
- Hofmann K, Bucher P: The Fha Domain - a Putative Nuclear Signaling Domain Found in Protein-Kinases and Transcription Factors. 1995, 20:Google Scholar
- Byeon IJL, Yongkiettrakul S, Tsai MD: Solution structure of the yeast Rad53 FHA2 complexed with a phosphothreonine peptide pTXXL: Comparison with the structures of FHA2-pYXL and FHA1-pTXXD complexes. 2001, 314:Google Scholar
- Lee GI, Ding ZF, Walker JC, Van Doren SR: NMR structure of the forkhead-associated domain from the Arabidopsis receptor kinase-associated protein phosphatase. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (20): 11261-11266. 10.1073/pnas.2031918100.View ArticleGoogle Scholar
- Ali AAE, Jukes RM, Pearl LH, Oliver AW: Specific recognition of a multiply phosphorylated motif in the DNA repair scaffold XRCC1 by the FHA domain of human PNK. 2009, 37:Google Scholar
- Nott TJ, Kelly G, Stach L, Li JJ, Westcott S, Patel D, Hunt DM, Howell S, Buxton RS, O'Hare HM, Smerdon SJ: An Intramolecular Switch Regulates Phosphoindependent FHA Domain Interactions in Mycobacterium tuberculosis. Science Signaling. 2009, 2 (63):Google Scholar
- Callebaut I, Mornon JP: From BRCA1 to RAP1: A widespread BRCT module closely associated with DNA repair. 1997, 400:Google Scholar
- Glover JNM, Williams RS, Lee MS: Interactions between BRCT repeats and phosphoproteins: tangled up in two. 2004, 29:Google Scholar
- Gough CA, Gojobori T, Imanishi T: Cancer-related mutations in BRCA1-BRCT cause long-range structural changes in protein-protein binding sites: A molecular dynamics study. 2007, 66:Google Scholar
- Watts FZ, Brissett NC: Linking up and interacting with BRCT domains. 2010, 9:Google Scholar
- Manke IA, Lowery DM, Nguyen A, Yaffe MB: BRCT repeats as phosphopeptide-binding modules involved in protein targeting. Science. 2003, 302 (5645): 636-639. 10.1126/science.1088877.View ArticleGoogle Scholar
- Williams RS, Lee MS, Hau DD, Glover JNM: Structural basis of phosphopeptide recognition by the BRCT domain of BRCA1. Nature Structural & Molecular Biology. 2004, 11 (6): 519-525. 10.1038/nsmb776.View ArticleGoogle Scholar
- Clapperton JA, Manke IA, Lowery DM, Ho T, Haire LF, Yaffe MB, Smerdon SJ: Structure and mechanism of BRCA1 BRCT domain recognition of phosphorylated BACH1 with implications for cancer. Nature Structural & Molecular Biology. 2004, 11 (6): 512-518. 10.1038/nsmb775.View ArticleGoogle Scholar
- Shiozaki EN, Gu LC, Yan N, Shi YG: Structure of the BRCT repeats of BRCA1 bound to a BACH1 phosphopeptide: Implications for signaling. Molecular Cell. 2004, 14 (3): 405-412. 10.1016/S1097-2765(04)00238-2.View ArticleGoogle Scholar
- Lu PJ, Zhou XZ, Liou YC, Noel JP, Lu KP: Critical role of WW domain phosphorylation in regulating phosphoserine binding activity and Pin1 function. Journal of Biological Chemistry. 2002, 277 (4): 2381-2384. 10.1074/jbc.C100228200.View ArticleGoogle Scholar
- Otte L, Wiedemann U, Schlegel B, Pires JR, Beyermann M, Schmieder P, Krause G, Volkmer-Engert R, Schneider-Mergener J, Oschkinat H: WW domain sequence activity relationships identified using ligand recognition propensities of 42 WW domains. Protein Science. 2003, 12 (3): 491-500. 10.1110/ps.0233203.View ArticleGoogle Scholar
- Verdecia MA, Bowman ME, Lu KP, Hunter T, Noel JP: Structural basis for phosphoserine-proline recognition by group IVWW domains. Nature Structural Biology. 2000, 7 (8): 639-643. 10.1038/77929.View ArticleGoogle Scholar
- Lee H, Yuan CH, Hammet A, Mahajan A, Chen ESW, Wu MR, Su MI, Heierhorst J, Tsai MD: Diphosphothreonine-specific interaction between an SQ/TQ cluster and an FHA domain in the Rad53-Dun1 kinase cascade. 2008, 30:Google Scholar
- Okur A, Strockbine B, Hornak V, Simmerling C: Using PC clusters to evaluate the transferability of molecular mechanics force fields for proteins. Journal of Computational Chemistry. 2003, 24 (1): 21-31. 10.1002/jcc.10184.View ArticleGoogle Scholar
- Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C: Comparison of multiple amber force fields and development of improved protein backbone parameters. Proteins-Structure Function and Bioinformatics. 2006, 65 (3): 712-725. 10.1002/prot.21123.View ArticleGoogle Scholar
- Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kale L, Schulten K: Scalable molecular dynamics with NAMD. Journal of Computational Chemistry. 2005, 26 (16): 1781-1802. 10.1002/jcc.20289.View ArticleGoogle Scholar
- Homeyer N, Horn AHC, Lanig H, Sticht H: AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine. Journal of Molecular Modeling. 2006, 12 (3): 281-289. 10.1007/s00894-005-0028-4.View ArticleGoogle Scholar
- Song YF, Gunner MR: Using Multiconformation Continuum Electrostatics to Compare Chloride Binding Motifs in alpha-Amylase, Human Serum Albumin, and Omp32. Journal of Molecular Biology. 2009, 387 (4): 840-856. 10.1016/j.jmb.2009.01.038.View ArticleGoogle Scholar
- Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: COMPARISON OF SIMPLE POTENTIAL FUNCTIONS FOR SIMULATING LIQUID WATER. Journal of Chemical Physics. 1983, 79 (2): 926-935. 10.1063/1.445869.View ArticleGoogle Scholar
- Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG: A SMOOTH PARTICLE MESH EWALD METHOD. Journal of Chemical Physics. 1995, 103 (19): 8577-8593. 10.1063/1.470117.View ArticleGoogle Scholar
- Ryckaert JP, Ciccotti G, Berendsen HJC: NUMERICAL-INTEGRATION OF CARTESIAN EQUATIONS OF MOTION OF A SYSTEM WITH CONSTRAINTS - MOLECULAR-DYNAMICS OF N-ALKANES. Journal of Computational Physics. 1977, 23 (3): 327-341. 10.1016/0021-9991(77)90098-5.View ArticleGoogle Scholar
- Sitkoff D, Sharp KA, Honig B: ACCURATE CALCULATION OF HYDRATION FREE-ENERGIES USING MACROSCOPIC SOLVENT MODELS. Journal of Physical Chemistry. 1994, 98 (7): 1978-1988. 10.1021/j100058a043.View ArticleGoogle Scholar
- Tsui V, Case DA: Theory and applications of the generalized Born solvation model in macromolecular Simulations. Biopolymers. 2000, 56 (4): 275-291. 10.1002/1097-0282(2000)56:4<275::AID-BIP10024>3.0.CO;2-E.View ArticleGoogle Scholar
- Bashford D, Case DA: Generalized born models of macromolecular solvation effects. Annual Review of Physical Chemistry. 2000, 51: 129-152. 10.1146/annurev.physchem.51.1.129.View ArticleGoogle Scholar
- Meirovitch H: Recent developments in methodologies for calculating the entropy and free energy of biological systems by computer simulation. Current Opinion in Structural Biology. 2007, 17 (2): 181-186. 10.1016/j.sbi.2007.03.016.View ArticleGoogle Scholar
- Chang CEA, McLaughlin WA, Baron R, Wang W, McCammon JA: Entropic contributions and the influence of the hydrophobic environment in promiscuous protein-protein association. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (21): 7456-7461. 10.1073/pnas.0800452105.View ArticleGoogle Scholar
- Wang J, Cai Q, Li ZL, Zhao HK, Luo R: Achieving energy conservation in Poisson-Boltzmann molecular dynamics: Accuracy and precision with finite-difference algorithms. Chemical Physics Letters. 2009, 468 (4-6): 112-118. 10.1016/j.cplett.2008.12.049.View ArticleGoogle Scholar
- Ye X, Wang J, Luo R: A Revised Density Function for Molecular Surface Calculation in Continuum Solvent Models. Journal of Chemical Theory and Computation. 2010, 6 (4): 1157-1169. 10.1021/ct900318u.View ArticleGoogle Scholar
- Tan C, Tan YH, Luo R: Implicit nonpolar solvent models. Journal of Physical Chemistry B. 2007, 111: 12263-12274. 10.1021/jp073399n.View ArticleGoogle Scholar
- Hnizdo V, Tan J, Killian BJ, Gilson MK: Efficient calculation of configurational entropy from molecular simulations by combining the mutual-information expansion and nearest-neighbor methods. Journal of Computational Chemistry. 2008, 29 (10): 1605-1614. 10.1002/jcc.20919.View ArticleGoogle Scholar
- Ai R, Fatmi MQ, Chang C-eA: T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles. Journal of Computer-Aided Molecular Design. 2010, 24 (10): 819-827. 10.1007/s10822-010-9376-y.View ArticleGoogle Scholar
- Joshi M, Ebalunode JO, Briggs JM: Computational insights into the interaction of the anthrax lethal factor with the N-terminal region of its substrates. Proteins-Structure Function and Bioinformatics. 2009, 75 (2): 323-335. 10.1002/prot.22241.View ArticleGoogle Scholar
- Hou T, Wang J, Li Y, Wang W: Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations. Journal of Chemical Information and Modeling. 2011, 51: 14-View ArticleGoogle Scholar
- Mandell DJ, Chorny I, Groban ES, Wong SE, Levine E, Rapp CS, Jacobson MP: Strengths of hydrogen bonds involving phosphorylated amino acid side chains. Journal of the American Chemical Society. 2007, 129: 820-827. 10.1021/ja063019w.View ArticleGoogle Scholar
- Khanna KK, Jackson SP: DNA double-strand breaks: signaling, repair and the cancer connection. 2001, 27:Google Scholar
- Traven A, Heierhorst J: SQ/TQ cluster domains: concentrated ATM/ATR kinase phosphorylation site regions in DNA-damage-response proteins. Bioessays. 2005, 27 (4): 397-407. 10.1002/bies.20204.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.