Universal two-point interaction of mediator KIX with 9aaTAD activation domains.

The nine-amino-acid activation domain (9aaTAD) is defined by a short amino acid pattern including two hydrophobic regions (positions p3-4 and p6-7). The KIX domain of mediator transcription CBP interacts with the 9aaTAD domains of transcription factors MLL, E2A, NF-kB, and p53. In this study, we analyzed the 9aaTADs-KIX interactions by nuclear magnetic resonance. The positions of three KIX helixes α1-α2-α3 are influenced by sterically-associated hydrophobic I611, L628, and I660 residues that are exposed to solvent. The positions of two rigid KIX helixes α1 and α2 generate conditions for structural folding in the flexible KIX-L12-G2 regions localized between them. The three KIX I611, L628, and I660 residues interact with two 9aaTAD hydrophobic residues in positions p3 and p4 and together build a hydrophobic core of five residues (5R). Numerous residues in 9aaTAD position p3 and p4 could provide this interaction. Following binding of the 9aaTAD to KIX, the hydrophobic I611, L628, and I660 residues are no longer exposed to solvent and their position changes inside the hydrophobic core together with position of KIX α1-α2-α3 helixes. The new positions of the KIX helixes α1 and α2 allow the KIX-L12-G2 enhanced formation. The second hydrophobic region of the 9aaTAD (positions p6 and p7) provides strong binding with the KIX-L12-G2 region. Similarly, multiple residues in 9aaTAD position p6 and p7 could provide this interaction. In conclusion, both 9aaTAD regions p3, p4 and p6, p7 provide co-operative and highly universal binding to mediator KIX. The hydrophobic core 5R formation allows new positions of the rigid KIX α-helixes and enables the enhanced formation of the KIX-L12-G2 region. This contributes to free energy and is the key for the KIX-9aaTAD binding. Therefore, the 9aaTAD-KIX interactions do not operate under the rigid key-and-lock mechanism what explains the 9aaTAD natural variability.

Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions.

Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.

Quantitation of Human 14-3-3ζ Dimerization and the Effect of Phosphorylation on Dimer-monomer Equilibria.

14-3-3 proteins are universal regulatory proteins and their function depends on their oligomeric form which may alter between the monomeric, homodimeric and heterodimeric states. The populations of individual oligomeric forms are controlled by Kd values of the dimer-monomer equilibria between the involved isoforms. This complex picture is extended by post-translational modifications, e.g. phosphorylation. In this work, we describe the equilibria between monomers, homo- and heterodimers of the 14-3-3ζ isoform in the unmodified and phosphorylated form. To cover a wide range of dimerization affinities, we combined solution NMR, microscale thermophoresis, native PAGE, and a set of novel fluorescence assays. Using a FRET based assay, we also determined the kinetic parameters of dimerization. We found that phosphorylation of 14-3-3ζ at Ser58 increases its homodimeric Kd value by 6 orders of magnitude. The presented assays allow to efficiently monitor 14-3-3ζ dimerization as a function of external factors, such as temperature, salt concentration, and client protein binding. For instance, we obtained values of both transient and equilibrium thermodynamic constants for the dimerization, and observed a substantial decrease of 14-3-3ζ dimer dissociation rate upon binding to the doubly phosphorylated regulatory domain of tyrosine hydroxylase. In summary, our work provides a conceptual framework to characterise the isoform exchanges of homo- and heterodimers which can significantly deepen our knowledge about the regulatory function of 14-3-3 proteins.

Phosphorylation of the regulatory domain of human tyrosine hydroxylase 1 monitored using non-uniformly sampled NMR.

Human tyrosine hydroxylase 1 (hTH1) activity is regulated by phosphorylation of its regulatory domain (RD-hTH1) and by an interaction with the 14-3-3 protein. The RD-hTH1 is composed of a structured region (66-169) preceded by an intrinsically disordered protein region (IDP, hTH1_65) containing two phosphorylation sites (S19 and S40) which are highly relevant for its increase in activity. The NMR signals of the IDP region in the non-phosphorylated, singly phosphorylated (pS40) and doubly phosphorylated states (pS19_pS40) were assigned by non-uniformly sampled spectra with increased dimensionality (5D). The structural changes induced by phosphorylation were analyzed by means of secondary structure propensities. The phosphorylation kinetics of the S40 and S19 by kinases PKA and PRAK respectively were monitored by non-uniformly sampled time-resolved NMR spectroscopy followed by their quantitative analysis.