Dissection of binding between a phosphorylated tyrosine hydroxylase peptide and 14-3-3zeta: A complex story elucidated by NMR

14-3-3 Family of Proteins: Biological Implications, Molecular Interactions, and Potential Intervention in Cancer, Virus and Neurodegeneration Disorders

The 14-3-3 family of proteins are highly conserved acidic eukaryotic proteins (25-32 kDa) abundantly present in the body. Through numerous binding partners, the 14-3-3 is responsible for many essential cellular pathways, such as cell cycle regulation and gene transcription control. Hence, its dysregulation has been linked to the onset of critical illnesses such as cancers, neurodegenerative diseases and viral infections. Interestingly, explorative studies have revealed an inverse correlation of 14-3-3 protein in cancer and neurodegenerative diseases, and the direct manipulation of 14-3-3 by virus to enhance infection capacity has dramatically extended its significance. Of these, COVID-19 has been linked to the 14-3-3 proteins by the interference of the SARS-CoV-2 nucleocapsid (N) protein during virion assembly. Given its predisposition towards multiple essential host signalling pathways, it is vital to understand the holistic interactions between the 14-3-3 protein to unravel its potential therapeutic unit in the future. As such, the general structure and properties of the 14-3-3 family of proteins, as well as their known biological functions and implications in cancer, neurodegeneration, and viruses, were covered in this review. Furthermore, the potential therapeutic target of 14-3-3 proteins in the associated diseases was discussed.

Hiding in plain sight: Complex interaction patterns between Tau and 14-3-3ζ protein variants

Tau protein is an intrinsically disordered protein that plays a key role in Alzheimer's disease (AD). In brains of AD patients, Tau occurs abnormally phosphorylated and aggregated in neurofibrillary tangles (NFTs). Together with Tau, 14-3-3 proteins - abundant cytosolic dimeric proteins - were found colocalized in the NFTs. However, so far, the molecular mechanism of the process leading to pathological changes in Tau structure as well as the direct involvement of 14-3-3 proteins are not well understood. Here, we aimed to reveal the effects of phosphorylation by protein kinase A (PKA) on Tau structural preferences and provide better insight into the interaction between Tau and 14-3-3 proteins. We also addressed the impact of monomerization-inducing phosphorylation of 14-3-3 at S58 on the binding to Tau protein. Using multidimensional nuclear magnetic resonance spectroscopy (NMR), chemical cross-linking analyzed by mass spectrometry (MS) and PAGE, we unveiled differences in their binding affinity, stoichiometry, and interfaces with single-residue resolution. We revealed that the interaction between 14-3-3 and Tau proteins is mediated not only via the 14-3-3 amphipathic binding grooves, but also via less specific interactions with 14-3-3 protein surface and, in the case of monomeric 14-3-3, also partially via the exposed dimeric interface. In addition, the hyperphosphorylation of Tau changes its affinity to 14-3-3 proteins. In conclusion, we propose quite complex interaction mode between the Tau and 14-3-3 proteins.

Characterizing the protein-protein interaction between MDM2 and 14-3-3σ; proof of concept for small molecule stabilization

Mouse Double Minute 2 (MDM2) is a key negative regulator of the tumor suppressor protein p53. MDM2 overexpression occurs in many types of cancer and results in the suppression of WT p53. The 14-3-3 family of adaptor proteins are known to bind MDM2 and the 14-3-3σ isoform controls MDM2 cellular localization and stability to inhibit its activity. Therefore, small molecule stabilization of the 14-3-3σ/MDM2 protein-protein interaction (PPI) is a potential therapeutic strategy for the treatment of cancer. Here, we provide a detailed biophysical and structural characterization of the phosphorylation-dependent interaction between 14-3-3σ and peptides that mimic the 14-3-3 binding motifs within MDM2. The data show that di-phosphorylation of MDM2 at S166 and S186 is essential for high affinity 14-3-3 binding and that the binary complex formed involves one MDM2 di-phosphorylated peptide bound to a dimer of 14-3-3σ. However, the two phosphorylation sites do not simultaneously interact so as to bridge the 14-3-3 dimer in a 'multivalent' fashion. Instead, the two phosphorylated MDM2 motifs 'rock' between the two binding grooves of the dimer, which is unusual in the context of 14-3-3 proteins. In addition, we show that the 14-3-3σ-MDM2 interaction is amenable to small molecule stabilization. The natural product fusicoccin A forms a ternary complex with a 14-3-3σ dimer and an MDM2 di-phosphorylated peptide resulting in the stabilization of the 14-3-3σ/MDM2 PPI. This work serves as a proof-of-concept of the drugability of the 14-3-3/MDM2 PPI and paves the way toward the development of more selective and efficacious small molecule stabilizers.

The aromatic amino acid hydroxylases: Structures, catalysis, and regulation of phenylalanine hydroxylase, tyrosine hydroxylase, and tryptophan hydroxylase

The aromatic amino acid hydroxylases phenylalanine hydroxylase, tyrosine hydroxylase, and tryptophan hydroxylase are non-heme iron enzymes that catalyze key physiological reactions. This review discusses the present understanding of the common catalytic mechanism of these enzymes and recent advances in understanding the relationship between their structures and their regulation.

Improving IDP theoretical chemical shift accuracy and efficiency through a combined MD/ADMA/DFT and machine learning approach

This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.

Recent advances in structural studies of 14-3-3 protein complexes

Being phosphopeptide-binding hubs, 14-3-3 proteins coordinate multiple cellular processes in eukaryotes, including the regulation of apoptosis, cell cycle, ion channels trafficking, transcription, signal transduction, and hormone biosynthesis. Forming constitutive α-helical dimers, 14-3-3 proteins predominantly recognize specifically phosphorylated Ser/Thr sites within their partners; this generally stabilizes phosphotarget conformation and affects its activity, intracellular distribution, dephosphorylation, degradation and interactions with other proteins. Not surprisingly, 14-3-3 complexes are involved in the development of a range of diseases and are considered promising drug targets. The wide interactome of 14-3-3 proteins encompasses hundreds of different phosphoproteins, for many of which the interaction is well-documented in vitro and in vivo but lack the structural data that would help better understand underlying regulatory mechanisms and develop new drugs. Despite obtaining structural information on 14-3-3 complexes is still lagging behind the research of 14-3-3 interactions on a proteome-wide scale, recent works provided some advances, including methodological improvements and accumulation of new interesting structural data, that are discussed in this review.

Phosphorylated and Phosphomimicking Variants May Differ—A Case Study of 14-3-3 Protein

Protein phosphorylation is a critical mechanism that biology uses to govern cellular processes. To study the impact of phosphorylation on protein properties, a fully and specifically phosphorylated sample is required although not always achievable. Commonly, this issue is overcome by installing phosphomimicking mutations at the desired site of phosphorylation. 14-3-3 proteins are regulatory protein hubs that interact with hundreds of phosphorylated proteins and modulate their structure and activity. 14-3-3 protein function relies on its dimeric nature, which is controlled by Ser58 phosphorylation. However, incomplete Ser58 phosphorylation has obstructed the detailed study of its effect so far. In the present study, we describe the full and specific phosphorylation of 14-3-3ζ protein at Ser58 and we compare its characteristics with phosphomimicking mutants that have been used in the past (S58E/D). Our results show that in case of the 14-3-3 proteins, phosphomimicking mutations are not a sufficient replacement for phosphorylation. At physiological concentrations of 14-3-3ζ protein, the dimer-monomer equilibrium of phosphorylated protein is much more shifted towards monomers than that of the phosphomimicking mutants. The oligomeric state also influences protein properties such as thermodynamic stability and hydrophobicity. Moreover, phosphorylation changes the localization of 14-3-3ζ in HeLa and U251 human cancer cells. In summary, our study highlights that phosphomimicking mutations may not faithfully represent the effects of phosphorylation on the protein structure and function and that their use should be justified by comparing to the genuinely phosphorylated counterpart.

Free energy calculations on the stability of the 14-3-3ζ protein

Mutations of cysteine are often introduced to e.g. avoid formation of non-physiological inter-molecular disulfide bridges in in-vitro experiments, or to maintain specificity in labeling experiments. Alanine or serine is typically preferred, which usually do not alter the overall protein stability, when the original cysteine was surface exposed. However, selecting the optimal mutation for cysteines in the hydrophobic core of the protein is more challenging. In this work, the stability of selected Cys mutants of 14-3-3ζ was predicted by free-energy calculations and the obtained data were compared with experimentally determined stabilities. Both the computational predictions as well as the experimental validation point at a significant destabilization of mutants C94A and C94S. This destabilization could be attributed to the formation of hydrophobic cavities and a polar solvation of a hydrophilic side chain. A L12E, M78K double mutant was further studied in terms of its reduced dimerization propensity. In contrast to naïve expectations, this double mutant did not lead to the formation of strong salt bridges, which was rationalized in terms of a preferred solvation of the ionic species. Again, experiments agreed with the calculations by confirming the monomerization of the double mutants. Overall, the simulation data is in good agreement with experiments and offers additional insight into the stability and dimerization of this important family of regulatory proteins.

Association of Multiple Phosphorylated Proteins with the 14-3-3 Regulatory Hubs: Problems and Perspectives

14-3-3 proteins are well-known universal regulators binding a vast number of partners by recognizing their phosphorylated motifs, typically located within the intrinsically disordered regions. The abundance of such phosphomotifs ensures the involvement of 14-3-3 proteins in sophisticated protein-protein interaction networks that govern vital cellular processes. Thousands of 14-3-3 partners have been either experimentally identified or predicted, but the spatiotemporal hierarchy of the processes based on 14-3-3 interactions is not clearly understood. This is exacerbated by the lack of available structural information on full regulatory complexes involving 14-3-3, which resist high-resolution structural studies due to the presence of intrinsically disordered regions. Although deducing three-dimensional structures is of particular urgency, structural advances are lagging behind the rate at which novel 14-3-3 partners are discovered. Here I attempted to critically review the current state of the field and in particular to dissect the unknowns, focusing on questions that could help in moving the frontiers forward.

Exploring the binding pathways of the 14-3-3ζ protein: Structural and free-energy profiles revealed by Hamiltonian replica exchange molecular dynamics with distancefield distance restraints

The 14-3-3 protein family performs regulatory functions in eukaryotic organisms by binding to a large number of phosphorylated protein partners. Whilst the binding mode of the phosphopeptides within the primary 14-3-3 binding site is well established based on the crystal structures of their complexes, little is known about the binding process itself. We present a computational study of the process by which phosphopeptides bind to the 14-3-3ζ protein. Applying a novel scheme combining Hamiltonian replica exchange molecular dynamics and distancefield restraints allowed us to map and compare the most likely phosphopeptide-binding pathways to the 14-3-3ζ protein. The most important structural changes to the protein and peptides involved in the binding process were identified. In order to bind phosphopeptides to the primary interaction site, the 14-3-3ζ adopted a newly found wide-opened conformation. Based on our findings we additionally propose a secondary interaction site on the inner surface of the 14-3-3ζ dimer, and a direct interference on the binding process by the flexible C-terminal tail. A minimalistic model was designed to allow for the efficient calculation of absolute binding affinities. Binding affinities calculated from the potential of mean force along the binding pathway are in line with the available experimental estimates for two of the studied systems.

Quantitative mapping of microtubule-associated protein 2c (MAP2c) phosphorylation and regulatory protein 14-3-3ζ-binding sites reveals key differences between MAP2c and its homolog Tau

Microtubule-associated protein 2c (MAP2c) is involved in neuronal development and is less characterized than its homolog Tau, which has various roles in neurodegeneration. Using NMR methods providing single-residue resolution and quantitative comparison, we investigated molecular interactions important for the regulatory roles of MAP2c in microtubule dynamics. We found that MAP2c and Tau significantly differ in the position and kinetics of sites that are phosphorylated by cAMP-dependent protein kinase (PKA), even in highly homologous regions. We determined the binding sites of unphosphorylated and phosphorylated MAP2c responsible for interactions with the regulatory protein 14-3-3ζ. Differences in phosphorylation and in charge distribution between MAP2c and Tau suggested that both MAP2c and Tau respond to the same signal (phosphorylation by PKA) but have different downstream effects, indicating a signaling branch point for controlling microtubule stability. Although the interactions of phosphorylated Tau with 14-3-3ζ are supposed to be a major factor in microtubule destabilization, the binding of 14-3-3ζ to MAP2c enhanced by PKA-mediated phosphorylation is likely to influence microtubule-MAP2c binding much less, in agreement with the results of our tubulin co-sedimentation measurements. The specific location of the major MAP2c phosphorylation site in a region homologous to the muscarinic receptor-binding site of Tau suggests that MAP2c also may regulate processes other than microtubule dynamics.

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.

Structural insights into the regulation of aromatic amino acid hydroxylation

The aromatic amino acid hydroxylases phenylalanine hydroxylase, tyrosine hydroxylase, and tryptophan hydroxylase are homotetramers, with each subunit containing a homologous catalytic domain and a divergent regulatory domain. The solution structure of the regulatory domain of tyrosine hydroxylase establishes that it contains a core ACT domain similar to that in phenylalanine hydroxylase. The isolated regulatory domain of tyrosine hydroxylase forms a stable dimer, while that of phenylalanine hydroxylase undergoes a monomer-dimer equilibrium, with phenylalanine stabilizing the dimer. These solution properties are consistent with the regulatory mechanisms of the two enzymes, in that phenylalanine hydroxylase is activated by phenylalanine binding to an allosteric site, while tyrosine hydroxylase is regulated by binding of catecholamines in the active site.