Search

International audience; 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 216PTPP219, 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.

The microtubule associated protein tau forms disease-specific filamentous aggregates in several different neurodegenerative diseases. In order to understand how tau undergoes misfolding into a specific filament type and to be able to control this process for drug development purposes, it is crucial to studyin vitrotau aggregation methods and investigate the structures of the obtained filaments at atomic level. Here, we used the tau fragment dGAE(297-391), which aggregates spontaneously, to seed formation of fulllength tau filaments. The structures of dGAE and full-length tau filaments were investigated by solid-state MAS NMR showing that dGAE allows propagating a chronic traumatic encephalopathy (CTE)-like fold to full-length tau. This work demonstrates thatin vitropreparation of disease-specific types of full-length tau filaments is feasible.

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

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.

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.

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 K

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.

Quantum mechanics (QM) calculations are applied to examine 1H, 13C, 15N, and 31P chemical shifts of two phosphorylation sites in an intrinsically disordered protein region. The QM calculations employ a combination of (1) structural ensembles generated by molecular dynamics, (2) a fragmentation technique based on the adjustable density matrix assembler, and (3) density functional methods. The combined computational approach is used to obtain chemical shifts (i) in the S19 and S40 residues of the nonphosphorylated and (ii) in the pS19 and pS40 residues of the doubly phosphorylated human tyrosine hydroxylase 1 as the system of interest. We study the effects of conformational averaging and explicit solvent sampling as well as the effects of phosphorylation on the computed chemical shifts. Good to great quantitative agreement with the experiment is achieved for all nuclei, provided that the systematic error cancellation is optimized by the choice of a suitable NMR standard. The effect of the standard reference on the computed 15N and 31P chemical shifts is demonstrated by employing three different referencing methods. Error bars associated with the statistical averaging of the computed 31P chemical shifts are larger than the difference between the 31P chemical shift of pS19 and pS40. The sequence trend of 31P shifts therefore could not be reliably reproduced. On the contrary, the calculations correctly predict the change of the 13C chemical shift for CB induced by the phosphorylation of the serine residues. The present work demonstrates that QM calculations coupled with molecular dynamics simulations and fragmentation techniques can be used as an alternative to empirical prediction tools in the structure characterization of intrinsically disordered proteins.

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.

The Ca2+-dependent clip-and-link activity of large repeat-in-toxin (RTX) proteins is an exceptional posttranslational process in which an internal domain called a self-processing module (SPM) mediates Ca2+-dependent processing of a highly specific aspartate-proline (Asp-Pro) peptide bond and covalent linkage of the released aspartyl to an adjacent lysine residue through an isopeptide bond. Here, we report the solution structures of the Ca2+-loaded SPM (Ca-SPM) defining the mechanism of the autocatalytic cleavage of the Asp414-Pro415 peptide bond of the Neisseria meningitidis FrpC exoprotein. Moreover, deletion of the SPM domain in the ApxIVA protein, the FrpC homolog of Actinobacillus pleuropneumoniae, resulted in attenuation of virulence of the bacterium in a pig infection model, indicating that the Ca2+-dependent clip-and-link activity plays a role in the virulence of Gram-negative pathogens., The posttranslational Ca2+-dependent “clip-and-link” activity of large repeat-in-toxin (RTX) proteins starts by Ca2+-dependent structural rearrangement of a highly conserved self-processing module (SPM). Subsequently, an internal aspartate-proline (Asp-Pro) peptide bond at the N-terminal end of SPM breaks, and the liberated C-terminal aspartyl residue can react with a free ε-amino group of an adjacent lysine residue to form a new isopeptide bond. Here, we report a solution structure of the calcium-loaded SPM (Ca-SPM) derived from the FrpC protein of Neisseria meningitidis. The Ca-SPM structure defines a unique protein architecture and provides structural insight into the autocatalytic cleavage of the Asp-Pro peptide bond through a “twisted-amide” activation. Furthermore, in-frame deletion of the SPM domain from the ApxIVA protein of Actinobacillus pleuropneumoniae attenuated the virulence of this porcine pathogen in a pig respiratory challenge model. We hypothesize that the Ca2+-dependent clip-and-link activity represents an unconventional strategy for Gram-negative pathogens to adhere to the host target cell surface.

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.

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 (CKI1RD) from Arabidopsis thaliana. We observed that the crystal structures of free, Mg2+-bound, and beryllofluoridated CKI1RD (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 CKI1RD variants exhibited different conformational dynamics in solution. NMR studies revealed that Mg2+ 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 CKI1RD was supported by a comparison with another A. thaliana histidine kinase, ETR1. In contrast to the highly dynamic β3–α3 loop of CKI1RD, the corresponding loop of the ETR1 receiver domain (ETR1RD) exhibited little conformational exchange and adopted a different orientation in crystals. Biochemical data indicated that ETR1RD is involved in phosphorylation-independent signaling, implying a direct link between conformational behavior and the ability of eukaryotic receiver domains to participate in MSP.

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.

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.

International audience; 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.

Computational simulations are used to study the structural and dynamics properties of biomoleculer systems at atomistic resolution. Morover, the simuations also allow to access the energetics of studied systems that can be applied in free energy calculations. Free energy determines the thermodynamic stability and solubility of biomoleclacues in given solution, their affinities towards another biomolecules and their populations in available conformational states. Traditional molecular dynamics simulations of biomacromolecules in explicit water solvent technique are currently restricted to the microseconds time scale but this limitation can be overcome by variety of enhanced sampling computational simulations.

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 naive 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.

The FtsZ protein is a self-polymerizing GTPase that plays a central role in bacterial cell division. Several C8-substituted GTP analogs are known to inhibit the polymerization of FtsZ by competing for the same binding site as its endogenous activating ligand GTP. Free energy calculations of the relative binding affinities to FtsZ for a set of five C8-substituted GTP analogs were performed. The calculated values agree well with the available experimental data, and the main contribution to the free energy differences is determined to be the conformational restriction of the ligands. The dihedral angle distributions around the glycosidic bond of these compounds in water are known to vary considerably depending on the physicochemical properties of the substituent at C8. However, within the FtsZ protein, this substitution has a negligible influence on the dihedral angle distributions, which fall within the narrow range of −140° to −90° for all investigated compounds. The corresponding ensemble average of the coupling constants 3 J(C4,H1′) is calculated to be 2.95 ± 0.1 Hz. The contribution of the conformational selection of the GTP analogs upon binding was quantified from the corresponding populations. The obtained restraining free energy values follow the same trend as the relative binding affinities to FtsZ, indicating their dominant contribution.

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.

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 byNMRspectroscopy 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. © 2009 American Chemical Society.

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.

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.

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.

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., Author Summary 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.

Recently, we compared amino acid sequences of the E2 glycoprotein of natural North American eastern equine encephalitis virus (NA-EEEV) isolates and demonstrated that naturally circulating viruses interact with heparan sulfate (HS) and that this interaction contributes to the extreme neurovirulence of EEEV (C. L. Gardner, G. D. Ebel, K. D. Ryman, and W. B. Klimstra, Proc. Natl. Acad. Sci. U. S. A., 108:16026–16031, 2011). In the current study, we have examined the contribution to HS binding of each of three lysine residues in the E2 71-to-77 region that comprise the primary HS binding site of wild-type (WT) NA-EEEV viruses. We also report that the original sequence comparison identified five virus isolates, each with one of three amino acid differences in the E2 71-to-77 region, including mutations in residues critical for HS binding by the WT virus. The natural variant viruses, which possessed either a mutation from lysine to glutamine at E2 71, a mutation from lysine to threonine at E2 71, or a mutation from threonine to lysine at E2 72, exhibited altered interactions with heparan sulfate and cell surfaces and altered virulence in a mouse model of EEEV disease. An electrostatic map of the EEEV E1/E2 heterotrimer based upon the recent Chikungunya virus crystal structure (J. E. Voss, M. C. Vaney, S. Duquerroy, C. Vonrhein, C. Girard-Blanc, E. Crublet, A. Thompson, G. Bricogne, and F. A. Rey, Nature, 468:709–712, 2010) showed the HS binding site to be at the apical surface of E2, with variants affecting the electrochemical nature of the binding site. Together, these results suggest that natural variation in the EEEV HS binding domain may arise during EEEV sylvatic cycles and that this variation may influence receptor interaction and the severity of EEEV disease.

Human APOBEC3A is a single-stranded DNA cytidine deaminase that restricts viral pathogens and endogenous retrotransposons, and has a role in the innate immune response. Furthermore, its potential to act as a genomic DNA mutator has implications for a role in carcinogenesis. A deeper understanding of APOBEC3A's deaminase and nucleic acid-binding properties, which is central to its biological activities, has been limited by the lack of structural information. Here we report the nuclear magnetic resonance solution structure of APOBEC3A and show that the critical interface for interaction with single-stranded DNA substrates includes residues extending beyond the catalytic centre. Importantly, by monitoring deaminase activity in real time, we find that A3A displays similar catalytic activity on APOBEC3A-specific TTCA- or A3G-specific CCCA-containing substrates, involving key determinants immediately 5' of the reactive C. Our results afford novel mechanistic insights into APOBEC3A-mediated deamination and provide the structural basis for further molecular studies.

As the most important phase I drug metabolizing enzymes, the human Cytochromes P450 display an enormous versatility in the molecular structures of possible substrates. Individual isoforms may preferentially bind specific classes of molecules, but also within these classes, some isoforms show remarkable levels of promiscuity. In this work, we try to link this promiscuity to the versatility and malleability of the active site at the hand of examples from our own work. Mainly focusing on the flexibility of protein structures and the presence or absence of water molecules, we establish molecular reasons for observed promiscuity, determine the relevant factors to take into account when predicting binding poses and rationalize the role of individual interactions in the process of ligand binding. A high level of active site flexibility does not only allow for the binding of a large variety of substrates and inhibitors, but also appears to be important to facilitate ligand binding and unbinding.

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. © 2010 American Chemical Society.

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.

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. © 2008 American Chemical Society.

To overcome the problem of insufficient conformational sampling within biomolecular simulations, we have developed a novel Hamiltonian replica exchange molecular dynamics (H-REMD) scheme that uses soft-core interactions between those parts of the system that contribute most to high energy barriers. The advantage of this approach over other H-REMD schemes is the possibility to use a relatively small number of replicas with locally larger differences between the individual Hamiltonians. Because soft-core potentials are almost the same as regular ones at longer distances, most of the interactions between atoms of perturbed parts will only be slightly changed. Rather, the strong repulsion between atoms that are close in space, which in many cases results in high energy barriers, is weakened within higher replicas of our proposed scheme. In addition to the soft-core interactions, we proposed to include multiple replicas using the same Hamiltonian/level of softness. We have tested the new protocol on the GTP and 8-Br-GTP molecules, which are known to have high energy barriers between the anti and syn conformation of the base with respect to the sugar moiety. During two 25 ns MD simulations of both systems the transition from the more stable to the less stable (but still experimentally observed) conformation is not seen at all. Also temperature REMD over 50 replicas for 1 ns did not show any transition at room temperature. On the other hand, more than 20 of such transitions are observed in H-REMD using six replicas (at three different Hamiltonians) during 6.8 ns per replica for GTP and 12 replicas (at six different Hamiltonians) during 8.7 ns per replica for 8-Br-GTP. The large increase in sampling efficiency was obtained from an optimized H-REMD scheme involving soft-core potentials, with multiple simulations using the same level of softness. The optimization of the scheme was performed by fast mimicking [J. Hritz and C. Oostenbrink, J. Chem. Phys. 127, 204104 (2007)]. © 2008 American Institute of Physics.

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. © 2007 American Institute of Physics.

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 −f0.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 (−f0.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. Proteins 2006. © 2006 Wiley-Liss, Inc.

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.

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.

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.

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 31P 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.