Hailey M.
Wallace‡
a,
Hyunjun
Yang‡
ab,
Sophia
Tan
a,
Henry S.
Pan
a,
Rose
Yang
a,
Junyi
Xu
a,
Hyunil
Jo
a,
Carlo
Condello
bc,
Nicholas F.
Polizzi
de and
William F.
DeGrado
*ab
aDepartment of Pharmaceutical Chemistry, The Cardiovascular Research Institution, University of California, San Francisco, CA 94158, USA. E-mail: Bill.DeGrado@ucsf.edu; Hailey.Wallace@ucsf.edu; Hyunjun.Yang@ucsf.edu; Sophia.Tan@ucsf.edu; Henry.Pan@ucsf.edu; Rose.Yang@ucsf.edu; Junyixu@berkeley.edu; Hyunil.Jo@ucsf.edu; Carlo.Condello@ucsf.edu; NicholasF_Polizzi@DFCI.harvard.edu
bInstitute for Neurodegenerative Diseases, University of California, San Francisco, CA 94143, USA
cDepartment of Neurology, University of California, San Francisco, CA 94143, USA
dDana Farber Cancer Institute, Harvard Medical School, Boston, MA 02215, USA
eDepartment of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA 02215, USA
First published on 30th March 2024
Insoluble amyloids rich in cross-β fibrils are observed in a number of neurodegenerative diseases. Depending on the clinicopathology, the amyloids can adopt distinct supramolecular assemblies, termed conformational strains. However, rapid methods to study amyloids in a conformationally specific manner are lacking. We introduce a novel computational method for de novo design of peptides that tile the surface of α-synuclein fibrils in a conformationally specific manner. Our method begins by identifying surfaces that are unique to the conformational strain of interest, which becomes a “target backbone” for the design of a peptide binder. Next, we interrogate structures in the PDB with high geometric complementarity to the target. Then, we identify secondary structural motifs that interact with this target backbone in a favorable, highly occurring geometry. This method produces monomeric helical motifs with a favorable geometry for interaction with the strands of the underlying amyloid. Each motif is then symmetrically replicated to form a monolayer that tiles the amyloid surface. Finally, amino acid sequences of the peptide binders are computed to provide a sequence with high geometric and physicochemical complementarity to the target amyloid. This method was applied to a conformational strain of α-synuclein fibrils, resulting in a peptide with high specificity for the target relative to other amyloids formed by α-synuclein, tau, or Aβ40. This designed peptide also markedly slowed the formation of α-synuclein amyloids. Overall, this method offers a new tool for examining conformational strains of amyloid proteins.
With the advancement of solid-state nuclear magnetic resonance (ssNMR) and cryo-electron microscopy (cryo-EM), the structures of many fibrils at the molecular resolution are now available.9–12 Depending on clinicopathology, each amyloidogenic protein can adopt multiple, distinct supramolecular assemblies, termed conformational strains.13–15 This conformational strain polymorphism of each amyloidogenic protein adds complexity to study amyloid-related disorders and, in some cases, multiple conformational strains are associated with a single disease (e.g. two conformers of tau amyloids are observed in Alzheimer's disease patients16). A method that could enable the identification and association of specific conformational strains at low cost (relative to, e.g., cryo-electron microscopy) would provide valuable insight into the nature of amyloidopathies. The development of amyloid conformational-specific binders will advance our fundamental understanding of the molecular recognition of amyloid structure for therapeutic strategies. Such binders will also enable the identification of polymorphs in vitro and in animal models, providing valuable reagents for histological examination of brain tissue. Ultimately, the knowledge gained in this process can contribute to the design of molecular therapeutics for the recognition and/or degradation of fibrils in a cellular context.
Current approaches have used small molecules, peptides, and antibodies as tools to recognize fibrils and disrupt their aggregation kinetics.17–26 However, these approaches exhibit varying degrees of conformational specificity. There are also molecular dyes that recognize differences in conformational strains and shift excitation/emission profiles; however, they often bind to multiple types of fibrils.27–29 Progress has been made in the field of protein design to engineer high-affinity fibril capping peptides that inhibit amyloid growth.23–26 However, these capping sites are limited in a growing fibril when compared to the immense fibril surface, providing only limited surface area for specific molecular recognition. These limitations underscore the feasibility of de novo designed binders targeting the surface of fibrils.
In this work, we have developed a computational method to design peptides that bind amyloids in a conformational strain-specific manner. We focused on in vitro fibrils for this model study because they are available in quantities for bulk studies and do not require biosafety requirements associated with handling human tissues. The designed peptides interact in a structure-specific manner with amyloids, and amyloid growth. Finally, our new computational approach should, in principle, allow for the development of new tools for the structure-based design of proteins that bind to other repeating protein surfaces.
Fig. 1 Cryo-EM structures of WT (green) and E46K (cyan) αSyn fibrils. The most discriminating residues Gln62, Thr64, and Val66 project outward in PDB 6cu7 but inward in PDB 6ufr. E46 and K46 are buried in each structure. |
To construct the target binding surface, we stacked the 5-residue fragment of the distinguishable region along the Z-axis, using the screw symmetry of the fibril cryo-EM structure (Fig. 2). We then extracted three neighboring strands and used the program MASTER33 to search for structures in the Protein Data Bank (PBD) that contained similar backbone coordinates to this motif. Because 3-stacked β-strands are highly prevalent in structures in the PDB, a large number of matches were retrieved (15395 matches at 1.0 Å RMSD). We then extracted the secondary structural motifs that are in the vicinity of the 3-strand motif (i.e., the target binding surface). In this process, interacting β-sheet motifs were filtered out because we expected that any hits conforming to the β-sheet conformation of an amyloid would probably self-assemble into an amyloid, even in the absence of aSyn fibrils. This resulted in a large number of helical poses, which were clustered based on the geometry of the interacting helix, its orientation relative to the target binding surface, and its frequency of occurrence (designability).
Having identified helical poses that could interact favorably with the underlying amyloid surface, we next used the screw symmetry of the fibril to create the overall assembly and design helix–helix interactions. The repeating geometry of the underlying fibrils dictated the stoichiometry and topology of the peptides relative to the target amyloid surface. The spacing of the β-strands is approximately 5 Å, restrained by the formation of inter-strand mainchain hydrogen bonds. On the other hand, the interhelical distance of helices within water-soluble proteins is approximately 8–11 Å, depending on the bulk of the residues lining the interface.34 These values dictated a stoichiometry of one helix for every two β-strands. Application of the fibril symmetry created a helical array with 9.6 Å spacing (twice the 4.8 Å spacing of the β-strands in the fibril). Moreover, the symmetry dictated a parallel orientation between the interacting helices.
Parallel packing of helices with a 10 Å spacing occurs quite frequently in proteins,34 suggesting that this should be a highly designable interface and may facilitate cooperative peptide binding on the amyloid surface. However, other interhelical parameters such as the α-helical phase and the precise interhelical crossing angle should also be favorable in the design.35 Thus, to identify helices suitable for helix–helix stacking along the amyloid surface, we again used MASTER to query two neighboring helices. The number of matches determined the feasibility of each helix–helix geometry, resulting in 75 helical backbones (>200 matches at 1.0 Å RMSD) that were well suited for the joint optimization of helix–strand and helix–helix interfaces in the next step of design.
We then ran a flexible backbone sequence design protocol using symmetry files for symmetric helix–helix packing within Rosetta.36–38 Rosetta symmetry files enable sequence design considering the symmetry of the fibril and the peptides, engendering favorable helix–strand and helix–helix interactions as a joint optimization problem. A total of 75000 designs were generated using flexible backbone design. To rank the designs, we used Rosetta's REF2015 score function. We took the energy score of [six β-strands·one helix] and subtracted the sum of [six β-strands] and [one helix] to estimate the helix–strand interaction. The helix–helix interaction score was similarly computed as the contribution from interactions in a symmetrical complex (Fig. S1†). During the selection process, designs were further filtered to minimize unsatisfied hydrogen bonds, ensuring only the most promising candidates were retained.
Peptide 1 | SWEEIIKALAELMRKVLEALSKL | Peptide 1* | FITC-Aβ-SWEEIIKALAELMRKVLEALSKL |
Peptide 2 | SMLEKVVEAGQKIVAK | Peptide 2* | FITC-Aβ-SMLEKVVEAGQKIVAK |
Peptide 3 | KMWELWARMAKEAFKMAESAYKS | Peptide 3* | FITC-Aβ-KMWELWARMAKEAFKMAESAYKS |
Peptide 4 | KMWELWARMAKEAFKMAESAAKS | Peptide 4* | FITC-Aβ-KMWELWARMAKEAFKMAESAAKS |
Scramble | ASKAKMAAEKRWMESLWEAMKAF | Scramble* | FITC-Aβ-ASKAKMAAEKRWMESLWEAMKAF |
FITC-β-Ala* | FITC-Aβ |
We examined the binding of our peptides to WT αSyn fibrils using confocal microscopy. Peptides 1*–4* were dissolved and diluted in PBS buffer at concentrations ranging from 10 nM to 0.1 nM with a final volume of 50 μL in a 384-well plate (Fig. 3A). To each well, 1 μL aliquot of 300 μM WT αSyn fibril was added and incubated for 30 min. The plate was centrifuged to pellet the fibrils for confocal microscopy. At 10 nM concentration, Peptide 4* displayed strong intensity indicating good binding. Even at concentrations as low as 1 nM, binding was still detectable with a good signal-to-noise ratio. The other three peptides 1*–3*, however, did not stain the fibrils brightly. It is interesting to note that peptide 3* was not bright despite having only a single alanine to tyrosine mutation from peptide 4*. Therefore, we focused additional characterization on peptide 4*.
To confirm the specificity of peptide 4* binding to WT αSyn fibrils, we compared the fluorescence intensity with E46K αSyn (Fig. 3B). Since the targeting region comprised of Gln62, Thr64 and Val66 is buried in the E46K αSyn fibril structure, no bright signal was expected. Indeed, we observed bright signal from the WT αSyn fibrils whereas only weak fluorescence from E46K αSyn which matched the intensity from the scramble* peptide. These data suggested that peptide 4* bound WT αSyn fibrils with good affinity and did not bind E46K αSyn fibrils validating our peptide design method to bind structure-defined amyloid conformer of interest. To rule out any nonspecific binding of the FITC dye, we also synthesized FITC-β-Ala, and no signal was observed with αSyn fibrils. We then examined if binding occurred with other types of amyloids such as Aβ40 and 0N3R tau fibrils (Fig. 3C). Peptide 4* showed dim fluorescence for Aβ40 and a low signal for 0N3R tau fibrils, in comparison to bright WT αSyn. Collectively, these data suggest that peptide 4* recognizes the targeted WT αSyn fibril conformational strain over E46K as well as other amyloid fibrils.
We then investigated whether peptide 4* could distinguish conformational strains of αSyn based solely on the structural differences, where the conformational strains each share identical amino-acid sequences. Three different wildtype αSyn fibril conformers (named WT, ribbon, and fibril) are produced in three different buffer conditions.40 Peptide 4* showed good intensity for the targeted WT αSyn conformer, dim signal for αSyn Ribbon, and low signal for αSyn fibril, respectively (Fig. 3D). These data clearly demonstrate that peptide 4* can effectively distinguish structural differences between conformational strains of the same amyloid sequence, highlighting the power of our peptide-design pipeline.
To gain structural insight into the peptide-amyloid complex, we first turned to protein–protein docking protocols from MOE and Hdock41 to predict the potential binding poses of the peptide bound to the target amyloid structure (Fig. 4A). When three copies of the designed helix were rigidly docked onto six β-strands of the fibril, the two methods were in good agreement with the design. A backbone RMSD difference of only 0.8 Å (relative to the computational design) was observed for peptide 4 when the helices were docked using Hdock and only slightly larger using MOE.
Fig. 4 (A) Docking prediction of the designs using MOE dock and Hdock. The backbone RMSD for design vs. prediction from MOE dock was 2.0 Å and 0.8 Å for Hdock. (B) MD simulation of design and the WT αSyn fibril (PDB 6cu7) with a plot showing the Cα RMSD calculation over 500 ns simulation. |
We next used 500 ns molecular dynamics (MD) simulations to assess the stability of the designed complex (Fig. 4B). Four copies of designed helices in an array were simulated in complex with ten β-strand chains for 500 ns and their helical backbone RMSDs were calculated to monitor stability of the initial configuration on this time scale. Peptides 1 and 2, which did not show appreciable binding in the confocal microscopy study, showed the largest deviations from the initial model. The functional peptide 4 was quite stable by this metric, indicating that MD was able to differentiate this binder from non-functional peptides. However, the Cα RMSD of peptide 3 was similar to peptide 4, indicating that stability on the 500 ns time scale was a necessary but not sufficient metric predictor of binding.
To determine whether the binding interaction was stoichiometric, we determined the binding affinity between peptide 4 and WT αSyn fibril using a ligand depletion sedimentation assay (Fig. 5A). WT αSyn fibrils were titrated into a fixed concentration of peptide 4*, and the amount of peptide remaining in solution was determined after ultracentrifugation. Our results displayed a saturation binding curve, indicating a plateau in binding. The data fits well to a binding isotherm with a KD = 1.3 μM and a stoichiometry of 1 peptide per 2 αSyn monomer, which matched the designed stoichiometry. It is, however, still possible that the peptide initially binds with higher affinity at low peptide/surface ratios and becomes weaker as saturation is approached. In a control experiment with the scramble* peptide, no appreciable binding to WT αSyn fibril was observed.
We next assessed the effect of peptide 4 on the rate of fibrillization of αSyn in the presence of thioflavin T (ThT), a fluorogenic probe for amyloid formation (Fig. 5B). We prepared fibrilization kinetics experiment with WT αSyn monomer in the presence and absence of peptide 4, each in of eight replicates. We noticed an increase in the fluorescence signal for WT αSyn at around 500 min, indicating αSyn fibrillization. This was compared to αSyn fibrilization with peptide 4 where a delayed increase in fluorescence is observed at around 2500 min and a slower increase of fluorescence signal. These data suggest that peptide 4 can interfere with the fibrilization kinetics. We suspect that peptide 4 tiles the solvent exposed fibril surface and prevents the secondary nucleation event.42 It is also possible that peptide 4 competes with ThT for the amyloid binding site, or peptide 4 favored one conformational strain over the other polymorphs. Collectively, peptide 4 slowed the fibrilization kinetics and did not act as an inducer.
More generally, this work presents a comprehensive approach for designing peptide or protein assemblies along any structured, repeating biological surface. In this method, the specificity of interaction is generated through both peptide–strand interactions and peptide–peptide interactions. The interaction of a single peptide is expected to be relatively weak, but by forming highly favorable peptide–peptide repeating interactions, the overall free energy of assembly becomes significantly favorable. Moreover, the specific interaction that each adjacent peptide forms with its neighbors stabilizes the geometry required for peptide-strand interactions. Altogether, a highly specific assembly is created, and might be widely useful as such.
Footnotes |
† Electronic supplementary information (ESI) available: Procedures for the design pipeline; preparation of αSyn fibrils (WT, E46K, rod, and fibril), and Aβ40 and 0N3R tau fibrils; peptide synthesis; ligand depletion sedimentation assay; ThT fluorescence assays; confocal imaging; characterization data (ESI-MS and HPLC traces) for peptide 4, peptide 4*, and scramble*. See DOI: https://doi.org/10.1039/d3sc06245g |
‡ H. M. W. and H. Y. contributed equally. |
This journal is © The Royal Society of Chemistry 2024 |