Molecular Mechanism of Thioflavin-T Binding to the Surface of β-Rich Peptide Self-Assemblies

A number of small organic molecules have been developed that bind to amyloid fibrils, a subset of which also inhibit fibrillization. Among these, the benzothiol dye Thioflavin-T (ThT) has been used for decades in the diagnosis of protein-misfolding diseases and in kinetic studies of self-assembly (fibrillization). Despite its importance, efforts to characterize the ThT-binding mechanism at the atomic level have been hampered by the inherent insolubility and heterogeneity of peptide self-assemblies. To overcome these challenges, we have developed a minimalist approach to designing a ThT-binding site in a "peptide self-assembly mimic” (PSAM) scaffold. PSAMs are engineered water-soluble proteins that mimic a segment of β-rich peptide self-assembly, and they are amenable to standard biophysical techniques and systematic mutagenesis. The PSAM β-sheet contains rows of repetitive amino acid patterns running perpendicular to the strands (cross-strand ladders) that represent a ubiquitous structural feature of fibril-like surfaces. We successfully designed a ThT-binding site that recapitulates the hallmarks of ThT-fibril interactions by constructing a cross-strand ladder consisting of contiguous tyrosines. The X-ray crystal structures suggest that ThT interacts with the β-sheet by docking onto surfaces formed by a single tyrosine ladder, rather than in the space between adjacent ladders. Systematic mutagenesis further demonstrated that tyrosine surfaces across four or more β-strands formed the minimal binding site for ThT. Our work thus provides structural insights into how this widely used dye recognizes a prominent subset of peptide self-assemblies, and proposes a strategy to elucidate the mechanisms of fibril-ligand interactions.     

High-resolution structure of a self-assembly-competent form of a hydrophobic peptide captured in a soluble beta-sheet scaffold

β-Rich self-assembly is a major structural class of polypeptides, but still little is known about its atomic structures and biophysical properties. Major impediments for structural and biophysical studies of peptide self-assemblies include their insolubility and heterogeneous composition. We have developed a model system, termed peptide self-assembly mimic (PSAM), based on the single-layer β-sheet of Borrelia outer surface protein A. PSAM allows for the capture of a defined number of self-assembly-like peptide repeats within a water-soluble protein, making structural and energetic studies possible. In this work, we extend our PSAM approach to a highly hydrophobic peptide sequence. We show that a penta-Ile peptide (Ile₅), which is insoluble and forms β-rich self-assemblies in aqueous solution, can be captured within the PSAM scaffold in a form capable of self-assembly. The 1.1-Å crystal structure revealed that the Ile₅ stretch forms a highly regular β-strand within this flat β-sheet. Self-assembly models built with multiple copies of the crystal structure of the Ile₅ peptide segment showed no steric conflict, indicating that this conformation represents an assembly-competent form. The PSAM retained high conformational stability, suggesting that the flat β-strand of the Ile₅ stretch primed for self-assembly is a low-energy conformation of the Ile₅ stretch and rationalizing its high propensity for self-assembly. The ability of the PSAM to “solubilize” an otherwise insoluble peptide stretch suggests the potential of the PSAM approach to the characterization of self-assembling peptides.     

Binding Modes of Thioflavin-T to the Single-Layer β-Sheet of the Peptide Self-Assembly Mimics

Although the amyloid dye thioflavin-T (ThT) is among the most widely used tools in the study of amyloid fibrils, the mechanism by which ThT binds to fibrils and other β-rich peptide self-assemblies remains elusive. The development of the water-soluble peptide self-assembly mimic (PSAM) system has provided a set of ideal model proteins for experimentally exploring the properties and minimal dye-binding requirements of amyloid fibrils. PSAMs consist of a single-layer β-sheet (SLB) capped by two globular domains, which capture the flat, extended β-sheet features common among fibril-like surfaces. Recently, a PSAM that binds to ThT with amyloid-like affinity (low micromolar K_d) has been designed, and its crystal structure in the absence of bound ThT was determined. This PSAM thus provides a unique opportunity to examine the interactions of ThT with a β-rich structure. Here, we present molecular dynamics simulations of the binding of ThT to this PSAM β-sheet. We show that the primary binding site for ThT is along a shallow groove formed by adjacent Tyr and Leu residues on the β-sheet surface. These simulations provide an atomic-scale rationale for this PSAM's experimentally determined dye-binding properties. Together, our results suggest that an aromatic-hydrophobic groove spanning across four consecutive β-strands represents a minimal ThT binding site on amyloid fibrils. Grooves formed by aromatic-hydrophobic residues on amyloid fibril surfaces may therefore offer a generic mode of recognition for amyloid dyes.     

Detection of initiation sites in protein folding of the four helix bundle ACBP by chemical shift analysis

A simple alternative method for obtaining "random coil" chemical shifts by intrinsic referencing using the protein's own peptide sequence is presented. These intrinsic random coil backbone shifts were then used to calculate secondary chemical shifts, that provide important information on the residual secondary structure elements in the acid-denatured state of an acyl-coenzyme A binding protein. This method reveals a clear correlation between the carbon secondary chemical shifts and the amide secondary chemical shifts 3-5 residues away in the primary sequence. These findings strongly suggest transient formation of short helix-like segments, and identify unique sequence segments important for protein folding.     

Alanine Scan of Core Positions in Ubiquitin Reveals Links between Dynamics,Stability, and Function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.     

Global changes in local protein dynamics reduce the entropic cost of carbohydrate binding in the arabinose-binding protein

Protein dynamics make important but poorly understood contributions to molecular recognition phenomena. To address this, we measure changes in fast protein dynamics that accompany the interaction of the arabinose-binding protein (ABP) with its ligand, d-galactose, using NMR relaxation and molecular dynamics simulation. These two approaches present an entirely consistent view of the dynamic changes that occur in the protein backbone upon ligand binding. Increases in the amplitude of motions are observed throughout the protein, with the exception of a few residues in the binding site, which show restriction of dynamics. These counter-intuitive results imply that a localised binding event causes a global increase in the extent of protein dynamics on the pico- to nanosecond timescale. This global dynamic change constitutes a substantial favourable entropic contribution to the free energy of ligand binding. These results suggest that the structure and dynamics of ABP may be adapted to exploit dynamic changes to reduce the entropic costs of binding.     

Residual structure, backbone dynamics, and interactions within the synuclein family

The human synuclein protein family includes alpha-synuclein, which has been linked to both familial and sporadic Parkinson's disease, and the highly homologous beta and gamma-synuclein. Mutations in alpha-synuclein cause autosomal dominant early onset Parkinson's, and the protein is found deposited in a fibrillar form in hereditary and idiopathic forms of the disease. No genetic link between beta and gamma-synuclein, and any neurodegenerative disease has been established, and it is generally considered that these proteins are not highly pathogenic. In addition, beta and gamma-synuclein are reported to aggregate less readily than alpha-synuclein in vitro. Indeed, beta-synuclein has been reported to protect against alpha-synuclein aggregation in vitro, as well as alpha-synuclein-mediated toxicity in vivo. Earlier, we compared the structural properties of the highly helical states adopted by all three synucleins in association with detergent micelles in an attempt to delineate the basis for functional differences between the three proteins. Here, we report a comparison of the structural and dynamic properties of the free states of all three proteins in order to shed light on differences that may help to explain their different propensities to aggregate, which in turn may underlie their differing contributions to the etiology of Parkinson's disease. We find that gamma-synuclein closely resembles alpha-synuclein in its free-state residual secondary structure, consistent with the more similar propensities of the two proteins to aggregate in vitro. beta-Synuclein, however, differs significantly from alpha-synuclein, exhibiting a lower predisposition towards helical structure in the second half of its lipid-binding domain, and a higher preference for extended structures in its C-terminal tail. Both beta and gamma-synuclein show less extensive transient long-range structure than that observed in alpha-synuclein. These results raise questions regarding the role of secondary structure propensities and transient long-range contacts in directing synuclein aggregation reactions.     

An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding–unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.     

The HIV fusion peptide adopts intermolecular parallel beta-sheet structure in membranes when stabilized by the adjacent N-terminal heptad repeat: a 13C FTIR study

The HIV gp41 protein mediates fusion with target host cells. The region primarily involved in directing fusion, the fusion peptide (FP), is poorly understood at the level of structure and function due to its toxic effect in expression systems. To overcome this, we used a synthetic approach to generate the N70 construct, whereby the FP is stabilized in context of the adjacent auto oligomerization domain. The amide I profile of unlabeled N70 in membranes reveals prominent alpha-helical contribution, along with significant beta-structure. By truncating the N terminus (FP region) of N70, beta-structure is eliminated, suggesting that the FP adopts a beta-structure in membranes. To assess this directly, (13)C Fourier-transformed infra-red analysis was carried out to map secondary structure of the 16 N-terminal hydrophobic residues of the fusion peptide (FP16). The (13)C isotope shifted absorbance of the FP was filtered from the global secondary structure of the 70 residue construct (N70). On the basis of the peak shift induced by the (13)C-labeled residues of FP16, we directly assign beta-sheet structure in ordered membranes. A differential labeling scheme in FP16 allows us to distinguish the type of beta-sheet structure as parallel. Dilution of each FP16-labeled N70 peptide, by mixing with unlabeled N70, shows directly that the FP16 beta-strand region self-assembles. We discuss our structural findings in the context of the prevailing gp41 fusion paradigm. Specifically, we address the role of the FP region in organizing supramolecular gp41 assembly, and we also discuss the mechanism by which exogenous, free FP constructs inhibit gp41-induced fusion.     

Interaction of the 106-126 prion peptide with lipid membranes and potential implication for neurotoxicity

Prion diseases are fatal neurodegenerative disorders characterized by the accumulation in the brain of an abnormally misfolded, protease-resistant, and beta-sheet rich pathogenic isoform (PrP(SC)) of the cellular prion protein (PrP(C)). In the present work, we were interested to study the mode of prion protein interaction with the membrane using the 106-126 peptide and small unilamellar lipid vesicles as model. As previously demonstrated, we showed by MTS assay that PrP 106-126 induces alterations in the human neuroblastoma SH-SY5Y cell line. We demonstrated for the first time by lipid-mixing assay and by the liposome vesicle leakage test that PrP 106-126, a non-tilted peptide, induces liposome fusion thus a potential cell membrane destabilization, as supported by membrane integrity assay (LDH). By circular dichroism (CD) analysis we showed that the fusogenic property of PrP 106-126 in the presence of liposome is associated with a predominantly beta-sheet structure. These data suggest that the fusogenic property associated with a predominant beta-sheet structure exhibited by the prion peptides contributes to the neurotoxicity of these peptides by destabilizing cellular membranes. The latter might be attached at the membrane surface in a parallel orientation as shown by molecular modeling.     

The cytoplasmic domain of the chloride channel ClC-0: structural and dynamic characterization of flexible regions

Eukaryotic members of the ClC family of chloride channels and transporters are composed of a transmembrane ion transport domain followed by a cytoplasmic domain, which is believed to be involved in the modulation of ClC function. In some family members this putative regulatory domain contains next to a well-folded structured part, long sequence stretches with low sequence complexity. These regions, a 96 residue long linker connecting two structured sub-domains, and 35 residues on the C teminus of the domain were found disordered in a recent crystal structure of this domain in ClC-0. Both regions have a large influence on the modulation of channel function in closely related family members. Here we describe a NMR study to characterize the structural and dynamic properties of these putatively unstructured stretches. Our study reveals that the two regions indeed show large conformational flexibility with dynamics on the nanosecond timescale. However, small islands of secondary structure are found interdispersed between the unfolded regions. This study characterizes for the first time the biophysical properties of these protein segments, which may become important for the understanding of novel regulatory mechanisms within the ClC family.     

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.     

Solution structure determination of the two DNA-binding domains in the Schizosaccharomyces pombe Abp1 protein by a combination of dipolar coupling and diffusion anisotropy restraints

We have solved the solution structure of the N-terminal region of the fission yeast centromere protein, Abp1, bound to a 21-base pair DNA fragment bearing its recognition site (Mw = 30 kDa). Although the two DNA-binding domains in the Abp1 protein were defined well by a conventional NOE-based NMR methodology, the overall structure of the Abp1 protein was poorly defined, due to the lack of interdomain distance restraints. Therefore, we additionally used residual dipolar couplings measured in a weakly aligned state, and rotational diffusion anisotropies. Neither the NH residual dipolar couplings nor the backbone ¹⁵N T ₁/T ₂ data were sufficient to determine the overall structure of the Abp1 protein, due to spectral overlap. We used a combination of these two orientational restraints (residual dipolar coupling and rotational diffusion anisotropy), which significantly improved the convergence of the overall structures. The range of the observed T ₁/T ₂ ratios was wider (20–50 for the secondary structure regions of Abp1) than the previously reported data for several globular proteins, indicating that the overall shape of the Abp1DNA complex is ellipsoid. This extended form would facilitate the recognition of the two separate sites in the relatively long DNA sequence by the DNA-binding domains of Apb1.     

Breaking symmetry in the structure determination of (large) symmetric protein dimers

We demonstrate a novel methodology to disrupt the symmetry in the NMR spectra of homodimers. A paramagnetic probe is introduced sub-stoichiometrically to create an asymmetric system with the paramagnetic probe residing on only one monomer within the dimer. This creates sufficient magnetic anisotropy for resolution of symmetry-related overlapped resonances and, consequently, detection of pseudocontact shifts and residual dipolar couplings specific to each monomeric component. These pseudocontact shifts can be readily incorporated into existing structure refinement calculations and enable determination of monomer orientation within the dimeric protein. This methodology can be widely used for solution structure determination of symmetric dimers.     

Three-dimensional structure of the weakly associated protein homodimer SeR13 using RDCs and paramagnetic surface mapping

The traditional NMR‐based method for determining oligomeric protein structure usually involves distinguishing and assigning intra‐ and intersubunit NOEs. This task becomes challenging when determining symmetric homo‐dimer structures because NOE cross‐peaks from a given pair of protons occur at the same position whether intra‐ or intersubunit in origin. While there are isotope‐filtering strategies for distinguishing intra from intermolecular NOE interactions in these cases, they are laborious and often prove ineffectual in cases of weak dimers, where observation of intermolecular NOEs is rare. Here, we present an efficient procedure for weak dimer structure determination based on residual dipolar couplings (RDCs), chemical shift changes upon dilution, and paramagnetic surface perturbations. This procedure is applied to the Northeast Structural Genomics Consortium protein target, SeR13, a negatively charged Staphylococcus epidermidis dimeric protein (K_d 3.4 ± 1.4 mM) composed of 86 amino acids. A structure determination for the monomeric form using traditional NMR methods is presented, followed by a dimer structure determination using docking under orientation constraints from RDCs data, and scoring under residue pair potentials and shape‐based predictions of RDCs. Validation using paramagnetic surface perturbation and chemical shift perturbation data acquired on sample dilution is also presented. The general utility of the dimer structure determination procedure and the possible relevance of SeR13 dimer formation are discussed.     

Measurements of residual dipolar couplings in peptide inhibitors weakly aligned by transient binding to peptide amyloid fibrils**

In this communication, we suggest that transferred residual dipolar couplings (trRDCs) can be employed to restrain the structure of peptide inhibitors transiently binding to beta-amyloid fibrils. The effect is based on the spontaneous alignment of amyloid fibrils with the fibril axis parallel to the magnetic field. This alignment is transferred to the transiently binding peptide inhibitor and is reflected in the size of the trRDCs. We find that the peptide inhibitor adopts a beta-sheet conformation with the backbone N-H and C-H dipolar vectors aligned preferentially parallel and perpendicular, respectively, to the fibril axis.     

Structural Characterization of Protein-Protein Complexes by Integrating Computational Docking with Small-angle Scattering Data

X-ray crystallography and NMR can provide detailed structural information of protein-protein complexes, but technical problems make their application challenging in the high-throughput regime. Other methods such as small-angle X-ray scattering (SAXS) are more promising for large-scale application, but at the cost of lower resolution, which is a problem that can be solved by complementing SAXS data with theoretical simulations. Here, we propose a novel strategy that combines SAXS data and accurate protein-protein docking simulations. The approach has been benchmarked on a large pool of known structures with synthetic SAXS data, and on three experimental examples. The combined approach (pyDockSAXS) provided a significantly better success rate (43% for the top 10 predictions) than either of the two methods alone. Further analysis of the influence of different docking parameters made it possible to increase the success rates for specific cases, and to define guidelines for improving the data-driven protein-protein docking protocols.     

The solution structure of antigen MPT64 from Mycobacterium tuberculosis defines a new family of beta-grasp proteins

The MPT64 protein and its homologs form a highly conserved family of secreted proteins with unknown function that are found within the pathogenic Mycobacteria genus. The founding member of this family from Mycobacterium tuberculosis (MPT64 or protein Rv1980c) is expressed only when Mycobacteria cells are actively dividing. By virtue of this relatively unique expression profile, Rv1980c is currently under phase III clinical trials to evaluate its potential to replace tuberculin, or purified protein derivative, as the rapid diagnostic of choice for detection of active tuberculosis infection. We describe here the NMR solution structure of Rv1980c. This structure reveals a previously undescribed fold that is based upon a variation of a beta-grasp motif most commonly found in protein-protein interaction domains. Examination of this structure in conjunction with multiple sequence alignments of MPT64 homologs identifies a candidate ligand-binding site, which may help guide future studies of Rv1980c function. The work presented here also suggests structure-based approaches for increasing the antigenic potency of a Rv1980c-based diagnostic.     

Solution structure of (gamma)S-crystallin by molecular fragment replacement NMR

The solution structure of murine gammaS-crystallin (gammaS) has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of dipolar couplings, recorded in different alignment media, and supplemented by a small number of NOE distance restraints. gammaS consists of two topologically similar domains, arranged with an approximate twofold symmetry, and each domain shows close structural homology to closely related (approximately 50% sequence identity) domains found in other members of the gamma-crystallin family. Each domain consists of two four-strand "Greek key" beta-sheets. Although the domains are tightly anchored to one another by the hydrophobic surfaces of the two inner Greek key motifs, the N-arm, the interdomain linker and several turn regions show unexpected flexibility and disorder in solution. This may contribute entropic stabilization to the protein in solution, but may also indicate nucleation sites for unfolding or other structural transitions. The method used for solving the gammaS structure relies on the recently introduced molecular fragment replacement method, which capitalizes on the large database of protein structures previously solved by X-ray crystallography and NMR.