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.

Shaping the future of preclinical development of successful disease-modifying drugs against Alzheimer's disease: a systematic review of tau propagation models

The transcellular propagation of the aberrantly modified protein tau along the functional brain network is a key hallmark of Alzheimer's disease and related tauopathies. Inoculation-based tau propagation models can recapitulate the stereotypical spread of tau and reproduce various types of tau inclusions linked to specific tauopathy, albeit with varying degrees of fidelity. With this systematic review, we underscore the significance of judicious selection and meticulous functional, biochemical, and biophysical characterization of various tau inocula. Furthermore, we highlight the necessity of choosing suitable animal models and inoculation sites, along with the critical need for validation of fibrillary pathology using confirmatory staining, to accurately recapitulate disease-specific inclusions. As a practical guide, we put forth a framework for establishing a benchmark of inoculation-based tau propagation models that holds promise for use in preclinical testing of disease-modifying drugs.

Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region

An increase in phosphorylation of the Tau protein is associated with Alzheimer's disease (AD) progression through unclear molecular mechanisms. In general, phosphorylation modifies the interaction of intrinsically disordered proteins, such as Tau, with other proteins; however, elucidating the structural basis of this regulation mechanism remains challenging. The bridging integrator-1 gene is an AD genetic determinant whose gene product, BIN1, directly interacts with Tau. The proline-rich motif recognized within a Tau(210-240) peptide by the SH3 domain of BIN1 (BIN1 SH3) is defined as ₂₁₆PTPP₂₁₉, and this interaction is modulated by phosphorylation. Phosphorylation of T217 within the Tau(210-240) peptide led to a 6-fold reduction in the affinity, while single phosphorylation at either T212, T231, or S235 had no effect on the interaction. Nonetheless, combined phosphorylation of T231 and S235 led to a 3-fold reduction in the affinity, although these phosphorylations are not within the BIN1 SH3-bound region of the Tau peptide. Using nuclear magnetic resonance (NMR) spectroscopy, these phosphorylations were shown to affect the local secondary structure and dynamics of the Tau(210-240) peptide. Models of the (un)phosphorylated peptides were obtained from molecular dynamics (MD) simulation validated by experimental data and showed compaction of the phosphorylated peptide due to increased salt bridge formation. This dynamic folding might indirectly impact the BIN1 SH3 binding by a decreased accessibility of the binding site. Regulation of the binding might thus not only be due to local electrostatic or steric effects from phosphorylation but also to the modification of the conformational properties of Tau.

Biomolecular Complexation on the "Wrong Side": A Case Study of the Influence of Salts and Sugars on the Interactions between Bovine Serum Albumin and Sodium Polystyrene Sulfonate

In the protein purification, drug delivery, food industry, and biotechnological applications involving protein-polyelectrolyte complexation, proper selection of co-solutes and solution conditions plays a crucial role. The onset of (bio)macromolecular complexation occurs even on the so-called "wrong side" of the protein isoionic point where both the protein and the polyelectrolyte are net like-charged. To gain mechanistic insights into the modulatory role of salts (NaCl, NaBr, and NaI) and sugars (sucrose and sucralose) in protein-polyelectrolyte complexation under such conditions, interaction between bovine serum albumin (BSA) and sodium polystyrene sulfonate (NaPSS) at pH = 8.0 was studied by a combination of isothermal titration calorimetry, fluorescence spectroscopy, circular dichroism, and thermodynamic modeling. The BSA-NaPSS complexation proceeds by two binding processes (first, formation of intrapolymer complexes and then formation of interpolymer complexes), both driven by favorable electrostatic interactions between the negatively charged sulfonic groups (-SO₃^-) of NaPSS and positively charged patches on the BSA surface. Two such positive patches were identified, each responsible for one of the two binding processes. The presence of salts screened both short-range attractive and long-range repulsive electrostatic interactions between both macromolecules, resulting in a nonmonotonic dependence of the binding affinity on the total ionic strength for both binding processes. In addition, distinct anion-specific effects were observed (NaCl < NaBr < NaI). The effect of sugars was less pronounced: sucrose had no effect on the complexation, but its chlorinated analogue, sucralose, promoted it slightly due to the screening of long-range repulsive electrostatic interactions between BSA and NaPSS. Although short-range non-electrostatic interactions are frequently mentioned in the literature in relation to BSA or NaPSS, we found that the main driving force of complexation on the "wrong side" are electrostatic interactions.

A new fibrillization mechanism of β-lactoglobulin in glycine solutions

Even though amyloid aggregates were discovered many years ago the mechanism of their formation is still a mystery. Because of their connection to many of untreatable neurodegenerative diseases the motivation for finding a common aggregation path is high. We report a new high heat induced fibrillization path of a model protein β-lactoglobulin (BLG) when incubated in glycine instead of water at pH 2. By combining atomic force microscopy (AFM), transmission emission microscopy (TEM), dynamic light scattering (DLS) and circular dichroism (CD) we predict that the basic building blocks of fibrils made in glycine are not peptides, but rather spheroid oligomers of different height that form by stacking of ring-like structures. Spheroid oligomers linearly align to form fibrils by opening up and combining. We suspect that glycine acts as an hydrolysation inhibitor which consequently promotes a different fibrillization path. By combining the known data on fibrillization in water with our experimental conclusions we come up with a new fibrillization scheme for BLG. We show that by changing the fibrillization conditions just by small changes in buffer composition can dramatically change the aggregation pathway and the effect of buffer shouldn't be neglected. Fibrils seen in our study are also gaining more and more attention because of their pore-like structure and a possible cytotoxic mechanism by forming pernicious ion-channels. By preparing them in a simple model system as BLG we opened a new way to study their formation.

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.

NMR Studies of Tau Protein in Tauopathies

Tauopathies, including Alzheimer's disease (AD), are the most troublesome of all age-related chronic conditions, as there are no well-established disease-modifying therapies for their prevention and treatment. Spatio-temporal distribution of tau protein pathology correlates with cognitive decline and severity of the disease, therefore, tau protein has become an appealing target for therapy. Current knowledge of the pathological effects and significance of specific species in the tau aggregation pathway is incomplete although more and more structural and mechanistic insights are being gained using biophysical techniques. Here, we review the application of NMR to structural studies of various tau forms that appear in its aggregation process, focusing on results obtained from solid-state NMR. Furthermore, we discuss implications from these studies and their prospective contribution to the development of new tauopathy therapies.

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.

Cysteine Modification by Ebselen Reduces the Stability and Cellular Levels of 14-3-3 Proteins

The 14-3-3 proteins constitute a family of adaptor proteins with many binding partners and biological functions, and they are considered promising drug targets in cancer and neuropsychiatry. By screening 1280 small-molecule drugs using differential scanning fluorimetry (DSF), we found 15 compounds that decreased the thermal stability of 14-3-3ζ Among these compounds, ebselen was identified as a covalent, destabilizing ligand of 14-3-3 isoforms ζ, ε, γ, and η Ebselen bonding decreased 14-3-3ζ binding to its partner Ser19-phosphorylated tyrosine hydroxylase. Characterization of site-directed mutants at cysteine residues in 14-3-3ζ (C25, C94, and C189) by DSF and mass spectroscopy revealed covalent modification by ebselen of all cysteines through a selenylsulfide bond. C25 appeared to be the preferential site of ebselen interaction in vitro, whereas modification of C94 was the main determinant for protein destabilization. At therapeutically relevant concentrations, ebselen and ebselen oxide caused decreased 14-3-3 levels in SH-SY5Y cells, accompanied with an increased degradation, most probably by the ubiquitin-dependent proteasome pathway. Moreover, ebselen-treated zebrafish displayed decreased brain 14-3-3 content, a freezing phenotype, and reduced mobility, resembling the effects of lithium, consistent with its proposed action as a safer lithium-mimetic drug. Ebselen has recently emerged as a promising drug candidate in several medical areas, such as cancer, neuropsychiatric disorders, and infectious diseases, including coronavirus disease 2019. Its pleiotropic actions are attributed to antioxidant effects and formation of selenosulfides with critical cysteine residues in proteins. Our work indicates that a destabilization of 14-3-3 may affect the protein interaction networks of this protein family, contributing to the therapeutic potential of ebselen. SIGNIFICANCE STATEMENT: There is currently great interest in the repurposing of established drugs for new indications and therapeutic targets. This study shows that ebselen, which is a promising drug candidate against cancer, bipolar disorder, and the viral infection coronavirus disease 2019, covalently bonds to cysteine residues in 14-3-3 adaptor proteins, triggering destabilization and increased degradation in cells and intact brain tissue when used in therapeutic concentrations, potentially explaining the behavioral, anti-inflammatory, and antineoplastic effects of this drug.

Polyethylene glycol perturbs the unfolding of CRABP I: A correlation between experimental and theoretical approach

The importance of macromolecules paves the way towards a detailed molecular level investigation as all most all cellular processes occurring at the interior of cells in the form of proteins, enzymes, and other biological molecules are significantly affected because of their crowding. Thus, exploring the role of crowding environment on the denaturation and renaturation kinetics of protein molecules is of great importance. Here, CRABP I (cellular retinoic acid binding protein I) is employed as a model protein along with different molecular weights of Polyethylene glycol (PEG) as molecular crowders. The experimental evaluations are done by accessing the protein secondary structure analysis using circular dichroism (CD) spectroscopy and unfolding kinetics using intrinsic fluorescence of CRABP I at 37 °C to mimic the in vivo crowding environment. The unfolding kinetics results indicated that both PEG 2000 and PEG 4000 act as stabilizers by retarding the unfolding kinetic rates. Both kinetic and stability outcomes presented the importance of crowding environment on stability and kinetics of CRABP I. The molecular dynamics (MD) studies revealed that thirteen PEG 2000 molecules assembled during the 500 ns simulation, which increases the stability and percentage of β-sheet. The experimental findings are well supported by the molecular dynamics simulation results.

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.

Structure and Functions of Microtubule Associated Proteins Tau and MAP2c: Similarities and Differences

The stability and dynamics of cytoskeleton in brain nerve cells are regulated by microtubule associated proteins (MAPs), tau and MAP2. Both proteins are intrinsically disordered and involved in multiple molecular interactions important for normal physiology and pathology of chronic neurodegenerative diseases. Nuclear magnetic resonance and cryo-electron microscopy recently revealed propensities of MAPs to form transient local structures and long-range contacts in the free state, and conformations adopted in complexes with microtubules and filamentous actin, as well as in pathological aggregates. In this paper, we compare the longest, 441-residue brain isoform of tau (tau40), and a 467-residue isoform of MAP2, known as MAP2c. For both molecules, we present transient structural motifs revealed by conformational analysis of experimental data obtained for free soluble forms of the proteins. We show that many of the short sequence motifs that exhibit transient structural features are linked to functional properties, manifested by specific interactions. The transient structural motifs can be therefore classified as molecular recognition elements of tau40 and MAP2c. Their interactions are further regulated by post-translational modifications, in particular phosphorylation. The structure-function analysis also explains differences between biological activities of tau40 and MAP2c.

Functionally specific binding regions of microtubule-associated protein 2c exhibit distinct conformations and dynamics

Microtubule-associated protein 2c (MAP2c) is a 49-kDa intrinsically disordered protein regulating the dynamics of microtubules in developing neurons. MAP2c differs from its sequence homologue Tau in the pattern and kinetics of phosphorylation by cAMP-dependent protein kinase (PKA). Moreover, the mechanisms through which MAP2c interacts with its binding partners and the conformational changes and dynamics associated with these interactions remain unclear. Here, we used NMR relaxation and paramagnetic relaxation enhancement techniques to determine the dynamics and long-range interactions within MAP2c. The relaxation rates revealed large differences in flexibility of individual regions of MAP2c, with the lowest flexibility observed in the known and proposed binding sites. Quantitative conformational analyses of chemical shifts, small-angle X-ray scattering (SAXS), and paramagnetic relaxation enhancement measurements disclosed that MAP2c regions interacting with important protein partners, including Fyn tyrosine kinase, plectin, and PKA, adopt specific conformations. High populations of polyproline II and α-helices were found in Fyn- and plectin-binding sites of MAP2c, respectively. The region binding the regulatory subunit of PKA consists of two helical motifs bridged by a more extended conformation. Of note, although MAP2c and Tau did not differ substantially in their conformations in regions of high sequence identity, we found that they differ significantly in long-range interactions, dynamics, and local conformation motifs in their N-terminal domains. These results highlight that the N-terminal regions of MAP2c provide important specificity to its regulatory roles and indicate a close relationship between MAP2c's biological functions and conformational behavior.

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.

Conformational dynamics are a key factor in signaling mediated by the receiver domain of a sensor histidine kinase from Arabidopsis thaliana

Multistep phosphorelay (MSP) cascades mediate responses to a wide spectrum of stimuli, including plant hormonal signaling, but several aspects of MSP await elucidation. Here, we provide first insight into the key step of MSP-mediated phosphotransfer in a eukaryotic system, the phosphorylation of the receiver domain of the histidine kinase CYTOKININ-INDEPENDENT 1 (CKI1_RD) from Arabidopsis thaliana We observed that the crystal structures of free, Mg²⁺-bound, and beryllofluoridated CKI1_RD (a stable analogue of the labile phosphorylated form) were identical and similar to the active state of receiver domains of bacterial response regulators. However, the three CKI1_RD variants exhibited different conformational dynamics in solution. NMR studies revealed that Mg²⁺ binding and beryllofluoridation alter the conformational equilibrium of the β3-α3 loop close to the phosphorylation site. Mutations that perturbed the conformational behavior of the β3-α3 loop while keeping the active-site aspartate intact resulted in suppression of CKI1 function. Mechanistically, homology modeling indicated that the β3-α3 loop directly interacts with the ATP-binding site of the CKI1 histidine kinase domain. The functional relevance of the conformational dynamics observed in the β3-α3 loop of CKI1_RD was supported by a comparison with another A. thaliana histidine kinase, ETR1. In contrast to the highly dynamic β3-α3 loop of CKI1_RD, the corresponding loop of the ETR1 receiver domain (ETR1_RD) exhibited little conformational exchange and adopted a different orientation in crystals. Biochemical data indicated that ETR1_RD is involved in phosphorylation-independent signaling, implying a direct link between conformational behavior and the ability of eukaryotic receiver domains to participate in MSP.

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.

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

Human tyrosine hydroxylase activity is regulated by phosphorylation of its N-terminus and by an interaction with the modulator 14-3-3 proteins. We investigated the binding of singly or doubly phosphorylated and thiophosphorylated peptides, comprising the first 50 amino acids of human tyrosine hydroxylase, isoform 1 (hTH1), that contain the critical interaction domain, to 14-3-3?, by (31)P NMR. Single phosphorylation at S19 generates a high affinity 14-3-3? binding epitope, whereas singly S40-phosphorylated peptide interacts with 14-3-3? one order-of-magnitude weaker than the S19-phosphorylated peptide. Analysis of the binding data revealed that the 14-3-3? dimer and the S19- and S40-doubly phosphorylated peptide interact in multiple ways, with three major complexes formed: 1), a single peptide bound to a 14-3-3? dimer via the S19 phosphate with the S40 phosphate occupying the other binding site; 2), a single peptide bound to a 14-3-3? dimer via the S19 phosphorous with the S40 free in solution; or 3), a 14-3-3? dimer with two peptides bound via the S19 phosphorous to each binding site. Our system and data provide information as to the possible mechanisms by which 14-3-3 can engage binding partners that possess two phosphorylation sites on flexible tails. Whether these will be realized in any particular interacting pair will naturally depend on the details of each system.

Sequence and structural determinants of human APOBEC3H deaminase and anti-HIV-1 activities

Background

Human APOBEC3H (A3H) belongs to the A3 family of host restriction factors, which are cytidine deaminases that catalyze conversion of deoxycytidine to deoxyuridine in single-stranded DNA. A3 proteins contain either one (A3A, A3C, A3H) or two (A3B, A3D, A3F, A3G) Zn-binding domains. A3H has seven haplotypes (I-VII) that exhibit diverse biological phenotypes and geographical distribution in the human population. Its single Zn-coordinating deaminase domain belongs to a phylogenetic cluster (Z3) that is different from the Z1- and Z2-type domains in other human A3 proteins. A3H HapII, unlike A3A or A3C, has potent activity against HIV-1. Here, we sought to identify the determinants of A3H HapII deaminase and antiviral activities, using site-directed sequence- and structure-guided mutagenesis together with cell-based, biochemical, and HIV-1 infectivity assays.

Results

We have constructed a homology model of A3H HapII, which is similar to the known structures of other A3 proteins. The model revealed a large cluster of basic residues (not present in A3A or A3C) that are likely to be involved in nucleic acid binding. Indeed, RNase A pretreatment of 293T cell lysates expressing A3H was shown to be required for detection of deaminase activity, indicating that interaction with cellular RNAs inhibits A3H catalytic function. Similar observations have been made with A3G. Analysis of A3H deaminase substrate specificity demonstrated that a 5′ T adjacent to the catalytic C is preferred. Changing the putative nucleic acid binding residues identified by the model resulted in reduction or abrogation of enzymatic activity, while substituting Z3-specific residues in A3H to the corresponding residues in other A3 proteins did not affect enzyme function. As shown for A3G and A3F, some A3H mutants were defective in catalysis, but retained antiviral activity against HIV-1vif (−) virions. Furthermore, endogenous reverse transcription assays demonstrated that the E56A catalytic mutant inhibits HIV-1 DNA synthesis, although not as efficiently as wild type.

Conclusions

The molecular and biological activities of A3H are more similar to those of the double-domain A3 proteins than to those of A3A or A3C. Importantly, A3H appears to use both deaminase-dependent and -independent mechanisms to target reverse transcription and restrict HIV-1 replication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-014-0130-8) contains supplementary material, which is available to authorized users.

Deliberate attenuation of chikungunya virus by adaptation to heparan sulfate-dependent infectivity: a model for rational arboviral vaccine design

Mosquito-borne chikungunya virus (CHIKV) is a positive-sense, single-stranded RNA virus from the genus Alphavirus, family Togaviridae, which causes fever, rash and severe persistent polyarthralgia in humans. Since there are currently no FDA licensed vaccines or antiviral therapies for CHIKV, the development of vaccine candidates is of critical importance. Historically, live-attenuated vaccines (LAVs) for protection against arthropod-borne viruses have been created by blind cell culture passage leading to attenuation of disease, while maintaining immunogenicity. Attenuation may occur via multiple mechanisms. However, all examined arbovirus LAVs have in common the acquisition of positively charged amino acid substitutions in cell-surface attachment proteins that render virus infection partially dependent upon heparan sulfate (HS), a ubiquitously expressed sulfated polysaccharide, and appear to attenuate by retarding dissemination of virus particles in vivo. We previously reported that, like other wild-type Old World alphaviruses, CHIKV strain, La Réunion, (CHIKV-LR), does not depend upon HS for infectivity. To deliberately identify CHIKV attachment protein mutations that could be combined with other attenuating processes in a LAV candidate, we passaged CHIKV-LR on evolutionarily divergent cell-types. A panel of single amino acid substitutions was identified in the E2 glycoprotein of passaged virus populations that were predicted to increase electrostatic potential. Each of these substitutions was made in the CHIKV-LR cDNA clone and comparisons of the mutant viruses revealed surface exposure of the mutated residue on the spike and sensitivity to competition with the HS analog, heparin, to be primary correlates of attenuation in vivo. Furthermore, we have identified a mutation at E2 position 79 as a promising candidate for inclusion in a CHIKV LAV.

With the adaptation of chikungunya virus (CHIKV) to transmission by the Aedes albopictus mosquito, a pandemic has occurred resulting in four to six million human infections, and the virus continues to become endemic in new regions, most recently in the Caribbean. CHIKV can cause debilitating polyarthralgia, lasting for weeks to years, and there are currently no licensed vaccines or antiviral therapies available. While an investigational live-attenuated vaccine (LAV) exists, problems with reactogenicity have precluded its licensure. The purpose of the current study was to: i) devise an in vitro passage procedure that reliably generates a panel of CHIKV envelope glycoprotein mutations for screening as vaccine candidates; ii) determine the position of the mutations in the three-dimensional structure of the alphavirus spike complex and their effect on electrostatic potential; iii) determine the attenuation characteristics of each mutation in a murine model of CHIKV musculoskeletal disease; and iv) to identify in vitro assays examining the dependency of infection upon HS that correlate with attenuation and localization in the glycoprotein spike. This approach provides a paradigm for the rational design of future LAVs for CHIKV and other mosquito-borne viruses, by deliberately selecting and combining attenuating processes.

Frequent mutation of receptor protein tyrosine phosphatases provides a mechanism for STAT3 hyperactivation in head and neck cancer

The underpinnings of STAT3 hyperphosphorylation resulting in enhanced signaling and cancer progression are incompletely understood. Loss-of-function mutations of enzymes that dephosphorylate STAT3, such as receptor protein tyrosine phosphatases, which are encoded by the PTPR gene family, represent a plausible mechanism of STAT3 hyperactivation. We analyzed whole exome sequencing (n = 374) and reverse-phase protein array data (n = 212) from head and neck squamous cell carcinomas (HNSCCs). PTPR mutations are most common and are associated with significantly increased phospho-STAT3 expression in HNSCC tumors. Expression of receptor-like protein tyrosine phosphatase T (PTPRT) mutant proteins induces STAT3 phosphorylation and cell survival, consistent with a "driver" phenotype. Computational modeling reveals functional consequences of PTPRT mutations on phospho-tyrosine-substrate interactions. A high mutation rate (30%) of PTPRs was found in HNSCC and 14 other solid tumors, suggesting that PTPR alterations, in particular PTPRT mutations, may define a subset of patients where STAT3 pathway inhibitors hold particular promise as effective therapeutic agents.

Role of water in molecular docking simulations of cytochrome P450 2D6

Active-site water molecules form an important component in biological systems, facilitating promiscuous binding or an increase in specificity and affinity. Taking water molecules into account in computational approaches to drug design or site-of-metabolism predictions is currently far from straightforward. In this study, the effects of including water molecules in molecular docking simulations of the important metabolic enzyme cytochrome P450 2D6 are investigated. The structure and dynamics of water molecules that are present in the active site simultaneously with a selected substrate are described, and based on this description, water molecules are selected to be included in docking experiments into multiple protein conformations. Apart from the parent substrate, 11 similar and 53 dissimilar substrates are included to investigate the transferability of active-site hydration sites between substrates. The role of water molecules appears to be highly dependent on the protein conformation and the substrate.

Efficient free energy calculations for compounds with multiple stable conformations separated by high energy barriers

Compounds with high intramolecular energy barriers represent challenging targets for free energy calculations because of the difficulty to obtain sufficient conformational sampling. Existing approaches are therefore computationally very demanding, thus preventing practical applications for such compounds. We present an enhanced sampling-one step perturbation method (ES-OS) to tackle this problem in a highly efficient way. A single molecular dynamics simulation of a judiciously chosen reference state (using two sets of soft-core interactions) is sufficient to determine conformational distributions of chemically similar compounds and the free energy differences between them. The ES-OS method is applied to a set of five biologically relevant 8-substituted GTP analogs having high energy barriers between the anti and the syn conformations of the base with respect to the ribose part. The reliability of ES-OS is verified by comparing the results to Hamiltonian replica exchange simulations of GTP and 8-Br-GTP and the experimentally determined 3J(C4,H1') coupling constant for GMP in water. Additional simulations in vacuum and octanol allow us to calculate differences in the solvation free energies and in lipophilicities (log P). Free energy contributions from individual conformational regions are also calculated, and their relationship with the overall free energy is derived leading to a set of multiconformational free energy formulas. These relationships are of general applicability and can be used in free energy calculations for a more diverse set of compounds.

A library of fluorescent peptides for exploring the substrate specificities of prolyl isomerases

To fully explore the substrate specificities of prolyl isomerases, we synthesized a library of 20 tetrapeptides that are labeled with a 2-aminobenzoyl (Abz) group at the amino terminus and a p-nitroanilide (pNA) group at the carboxy terminus. In this peptide library of the general formula Abz-Ala-Xaa-Pro-Phe-pNA, the position Xaa before the proline is occupied by all 20 proteinogenic amino acids. A conformational analysis of the peptide by molecular dynamics simulations and by NMR spectroscopy showed that the mutual distance between the Abz and pNA moieties in the peptides depends on the isomeric state of the Xaa-Pro bond. In the cis, but not in the trans form, there are significant chemical shift changes of the Abz and pNA moieties, because their aromatic rings are close to each other. This proximity also leads to a strong quenching of Abz fluorescence, which, in combination with a solvent jump, was used to devise a sensitive assay for prolyl isomerases. Unlike the traditional assay, it is not coupled with peptide proteolysis and thus can be employed for protease-sensitive prolyl isomerases as well. The peptide library was used to provide a complete set of P1-site specificities for prototypic human members of the three prolyl isomerase families, FKBP12, cyclophilin 18, and parvulin 14. In a second application, the substrate specificity of SlyD, a protease-sensitive prolyl isomerase from Escherichia coli, was characterized and compared with that of human FKBP12 as well as with homologues from other bacteria.

Impact of plasticity and flexibility on docking results for cytochrome P450 2D6: a combined approach of molecular dynamics and ligand docking

Cytochrome P450s (CYPs) exhibit a large plasticity and flexibility in the active site allowing for the binding of a large variety of substrates. The impact of plasticity and flexibility on ligand binding is investigated by docking 65 known CYP2D6 substrates to an ensemble of 2500 protein structures. The ensemble was generated by molecular dynamics simulations of CYP2D6 in complex with five representative substrates. The effect of induced fit, the conformation of Phe483, and thermal motion on the accuracy of site of metabolism (SOM) predictions is analyzed. For future predictions, the three most essential CYP2D6 structures were selected which are suitable for different kinds of ligands. We have developed a binary decision tree to decide which protein structure to dock the ligand into, such that each ligand needs to be docked only once, leading to successful SOM prediction in 80% of the substrates.

Virtual screening and prediction of site of metabolism for cytochrome P450 1A2 ligands

With the availability of an increasing number of high resolution 3D structures of human cytochrome P450 enzymes, structure-based modeling tools are more readily used. In this study we explore the possibilities of using docking and scoring experiments on cytochrome P450 1A2. Three different questions have been addressed: 1. Binding orientations and conformations were successfully predicted for various substrates. 2. A virtual screen was performed with satisfying enrichment rates. 3. A classification of individual compounds into active and inactive was performed. It was found that while docking can be used successfully to address the first two questions, it seems to be more difficult to perform the classification. Different scoring functions were included, and the well-characterized water molecule in the active site was included in various ways. Results are compared to experimental data and earlier classification data using machine learning methods. The possibilities and limitations of using structure-based drug design tools for cytochrome P450 1A2 come to light and are discussed.

Optimization of replica exchange molecular dynamics by fast mimicking

We present an approach to mimic replica exchange molecular dynamics simulations (REMD) on a microsecond time scale within a few minutes rather than the years, which would be required for real REMD. The speed of mimicked REMD makes it a useful tool for "testing" the efficiency of different settings for REMD and then to select those settings, that give the highest efficiency. We present an optimization approach with the example of Hamiltonian REMD using soft-core interactions on two model systems, GTP and 8-Br-GTP. The optimization process using REMD mimicking is very fast. Optimization of Hamiltonian-REMD settings of GTP in explicit water took us less than one week. In our study we focus not only on finding the optimal distances between neighboring replicas, but also on finding the proper placement of the highest level of softness. In addition we suggest different REMD simulation settings at this softness level. We allow several replicas to be simulated at the same Hamiltonian simultaneously and reduce the frequency of switching attempts between them. This approach allows for more efficient conversions from one stable conformation to the other.

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.

Characterization of the interaction of hypericin with protein kinase C in U-87 MG human glioma cells

A fluorescence imaging technique was used to monitor intracellular localization of protein kinase C (PKC) in U-87 MG human glioma cells in the presence of hypericin (Hyp) and phorbol 12-myristate-13-acetate (PMA). It is shown that PKC localization, which reflects its activity, is influenced by Hyp and this influence is different from that observed for PMA which acts as PKC activator. Fluorescence binding experiments were used to determine the binding constants of Hyp to several isoforms of PKC. The obtained values of K(d)s (approximately 100 nM) suggest that Hyp binds with high affinity to PKC. Finally, molecular modeling was used to compare structural models of the interaction of C1B domain of PKC (PKC isoforms alpha, delta, gamma) with Hyp and our previously published model of the (C1B domain PKCgamma)/PMA complex. The influence of Hyp on PKC translocation in U-87 MG cells in comparison with PMA, colocalization fluorescence pattern of Hyp and PKC, the higher binding affinity of Hyp to PKC in comparison with known binding constants of phorbol esters, as well as the binding mode of Hyp and PMA to the C1B domain of PKC suggested by molecular modeling, support the idea that Hyp and PMA might competitively bind to the regulatory domain of PKC.

Molecular Interaction Model for the C1B Domain of Protein Kinase C-γ in the Complex with Its Activator Phorbol-12-myristate-13-acetate in Water Solution and Lipid Bilayer

Detailed molecular models of the free C1B domain of protein kinase C-gamma (PKC-gamma) and the C1B domain with its activator phorbol-12-myristate-13-acetate (PMA) in water solution and in the presence of dipalmitoylphosphatidylcholine (DPPC) bilayer are presented. Molecular dynamics of the free C1B domain reveals hydrogen bonds, which are critical for the forming of the diacylglycerols/phorbol esters binding site, and indicates the important role of Gln27 for the geometry of this site. According to the model, PMA interacts with the C1B domain by hydrophobic interactions with Pro11 and Tyr22 and by three persistent hydrogen bonds between the C3 carbonyl group of PMA and Gly23 and between the C20 hydroxyl group of PMA and the Leu21 and Thr12 residues of the C1B domain. The C9 hydroxyl group of PMA does not interact with the C1B domain, but it is involved in the interaction with the DPPC bilayer. Two preferential orientations of the C1B-PMA complex toward the DPPC bilayer resulted from our molecular modeling study.

Toc12, a novel subunit of the intermembrane space preprotein translocon of chloroplasts

Translocation of proteins across membranes is essential for the biogenesis of each cell and is achieved by proteinaceous complexes. We analyzed the translocation complex of the intermembrane space from chloroplasts and identified a 12-kDa protein associated with the Toc machinery. Toc12 is an outer envelope protein exposing a soluble domain into the intermembrane space. Toc12 contains a J-domain and stimulates the ATPase activity of DnaK. The conformational stability and the ability to stimulate Hsp70 are dependent on a disulfide bridge within the loop region of the J-domain, suggesting a redox-regulated activation of the chaperone. Toc12 is associated with Toc64 and Tic22. Its J-domain recruits the Hsp70 of outer envelope membrane to the intermembrane space translocon and facilitates its interaction to the preprotein.

Influence of structure of human, rat, and bovine serum albumins on binding properties of photoactive drug hypericin

Fluorescence binding measurements and molecular modeling were employed to study the interaction of hypericin (Hyp) with human (HSA), rat (RSA), and bovine (BSA) serum albumins. Fluorescence emission data show the solubility of Hyp increasing in the order BSA, HSA, and RSA. Molecular modeling was used to construct the detailed structural models of the complexes and to explain the differences in the binding properties of Hyp. It was shown that the structures of Hyp/HSA and Hyp/RSA complexes are more similar and in some aspects different from those found for the Hyp/BSA complex. The role of the amino acid sequence in the IIA subdomains of HSA, RSA, and BSA is discussed to explain the observed differences.

Interaction of hypericin with serum albumins: surface-enhanced Raman spectroscopy, resonance Raman spectroscopy and molecular modeling study

Surface-enhanced Raman spectroscopy, resonance Raman spectroscopy and molecular modeling were employed to study the interaction of hypericin (Hyp) with human (HSA), rat (RSA) and bovine (BSA) serum albumins. The identification of the binding site of Hyp in serum albumins as well as the structural model for Hyp/HSA complex are presented. The interactions mainly reflect: (1) a change of the strength of H bonding at the N1-H site of Trp; (2) a change of the Trp side-chain conformation; (3) a change of the hydrophobicity of the Trp environment; and (4) a formation of an H-bond between the carbonyl group of Hyp and a proton donor in HSA and RSA which leads to a protonated-like carbonyl in Hyp. Our results indicate that Hyp is rigidly bound in IIA subdomain of HSA close to Trp214 (distance 5.12 A between the centers of masses). In the model presented the carbonyl group of Hyp is hydrogen bonded to Asn458. Two other candidates for hydrogen bonds have been identified between the bay-region hydroxyl group of Hyp and the carbonyl group of the Trp214 peptidic link and between the peri-region hydroxyl group of Hyp and the Asn458 carbonyl group. It is shown that the structures of the Hyp/HSA and Hyp/RSA complexes are similar to, and in some aspects different from, those found for the Hyp/BSA complex. The role of aminoacid sequence in the IIA subdomains of HSA, RSA and BSA is discussed to explain the observed differences.

Interaction of hypericin with serum albumins: surface-enhanced Raman spectroscopy, resonance Raman spectroscopy and molecular modeling study

Surface-enhanced Raman spectroscopy, resonance Raman spectroscopy and molecular modeling were employed to study the interaction of hypericin (Hyp) with human (HSA), rat (RSA) and bovine (BSA) serum albumins. The identification of the binding site of Hyp in serum albumins as well as the structural model for Hyp/HSA complex are presented. The interactions mainly reflect: (1) a change of the strength of H bonding at the N1-H site of Trp; (2) a change of the Trp side-chain conformation; (3) a change of the hydrophobicity of the Trp environment; and (4) a formation of an H-bond between the carbonyl group of Hyp and a proton donor in HSA and RSA which leads to a protonated-like carbonyl in Hyp. Our results indicate that Hyp is rigidly bound in IIA subdomain of HSA close to Trp214 (distance 5.12 A between the centers of masses). In the model presented the carbonyl group of Hyp is hydrogen bonded to Asn458. Two other candidates for hydrogen bonds have been identified between the bay-region hydroxyl group of Hyp and the carbonyl group of the Trp214 peptidic link and between the peri-region hydroxyl group of Hyp and the Asn458 carbonyl group. It is shown that the structures of the Hyp/HSA and Hyp/RSA complexes are similar to, and in some aspects different from, those found for the Hyp/BSA complex. The role of aminoacid sequence in the IIA subdomains of HSA, RSA and BSA is discussed to explain the observed differences.