Testing the feasibility of targeting a conserved region on the S2 domain of the SARS-CoV-2 spike protein.

The efficacy of vaccines against the SARS-CoV-2 virus significantly declines with the emergence of mutant strains, prompting investigation into the feasibility of targeting highly conserved but often cryptic regions on the S2 domain of spike protein. Using tools from molecular dynamics, we find that exposure of a conserved S2 epitope located in the central helices below the receptor binding domains would require large-scale motion beyond receptor binding domain up-down motion, but, along the reaction coordinates we explored, it is unlikely to be exposed by such large-scale dynamic fluctuations of the S1 domain without any external facilitating factors, despite some previous computational evidence suggesting transient exposure of this region. Furthermore, glycans, particularly those on N165 and N234, hinder S2-exposing opening dynamics, and thus stabilize spike in addition to immunologically shielding the protein surface. Although the S2 epitope region examined here is central to large-scale conformational changes during viral entry, free energy landscape analysis obtained using the path coordinate formalism reveals no inherent "loaded spring" effect, suggesting that a vaccine immunogen would tend to present the epitope in a prefusion-like conformation and may be effective in neutralization. These findings contribute to a deeper understanding of the dynamic origins of the function of the spike protein, as well as further characterizing the feasibility of the S2 epitope as a therapeutic target.

Cyclization Scaffolding for Improved Vaccine Immunogen Stability: Application to Tau Protein in Alzheimer's Disease.

Effective scaffolding of immunogens is crucial for generating conformationally selective antibodies through active immunization, particularly in the treatment of protein misfolding diseases such as Alzheimer's and Parkinson's disease. Previous computational work has revealed that a disorder-prone region of the tau protein, when in a stacked form, is predicted to structurally resemble a small, soluble protofibril, having conformational properties similar to those of experimental in vitro tau oligomers. Such an oligomeric structural mimic has the potential to serve as a vaccine immunogen design for Alzheimer's disease. In this study, we developed a cyclization scaffolding method in Rosetta, in which multiple cyclic peptides are stacked into a protofibril. Cyclization results in significant stabilization of protofibril-like structures by constraining the conformational space. Applying this method to the disorder-prone region of the tau fibril, we evaluated the metastability of the cyclized tau immunogen using molecular dynamics simulations, and we identified sequences of two cyclic constructs having high metastability in the protofibril. We then assessed their thermodynamic stability by computing the free energy required to separate a distal chain from the rest of the stacked structure. Our computational results, based on molecular dynamics simulations and free energy calculations, demonstrate that two cyclized constructs, cyclo-(VKSEKLDFKDRVQSKIFyN) and cyclo-(VKSEKLDFKDRVQSKIYvG) (lowercase letters indicate d-form amino acids), possess significantly increased thermodynamic stability in the protofibril over an uncyclized linear construct VKSEKLDFKDRVQSKI. The cyclization scaffolding approach proposed here holds promise as a means to effectively design immunogens for protein misfolding diseases, particularly those involving liposome-conjugated peptide constructs.

Targeting RACK1 to alleviate TDP-43 and FUS proteinopathy-mediated suppression of protein translation and neurodegeneration.

TAR DNA-binding protein 43 (TDP-43) and Fused in Sarcoma/Translocated in Sarcoma (FUS) are ribonucleoproteins associated with pathogenesis of amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). Under physiological conditions, TDP-43 and FUS are predominantly localized in the nucleus, where they participate in transcriptional regulation, RNA splicing and metabolism. In disease, however, they are typically mislocalized to the cytoplasm where they form aggregated inclusions. A number of shared cellular pathways have been identified that contribute to TDP-43 and FUS toxicity in neurodegeneration. In the present study, we report a novel pathogenic mechanism shared by these two proteins. We found that pathological FUS co-aggregates with a ribosomal protein, the Receptor for Activated C-Kinase 1 (RACK1), in the cytoplasm of spinal cord motor neurons of ALS, as previously reported for pathological TDP-43. In HEK293T cells transiently transfected with TDP-43 or FUS mutant lacking a functional nuclear localization signal (NLS; TDP-43ΔNLS and FUSΔNLS), cytoplasmic TDP-43 and FUS induced co-aggregation with endogenous RACK1. These co-aggregates sequestered the translational machinery through interaction with the polyribosome, accompanied by a significant reduction of global protein translation. RACK1 knockdown decreased cytoplasmic aggregation of TDP-43ΔNLS or FUSΔNLS and alleviated associated global translational suppression. Surprisingly, RACK1 knockdown also led to partial nuclear localization of TDP-43ΔNLS and FUSΔNLS in some transfected cells, despite the absence of NLS. In vivo, RACK1 knockdown alleviated retinal neuronal degeneration in transgenic Drosophila melanogaster expressing hTDP-43WT or hTDP-43Q331K and improved motor function of hTDP-43WT flies, with no observed adverse effects on neuronal health in control knockdown flies. In conclusion, our results revealed a novel shared mechanism of pathogenesis for misfolded aggregates of TDP-43 and FUS mediated by interference with protein translation in a RACK1-dependent manner. We provide proof-of-concept evidence for targeting RACK1 as a potential therapeutic approach for TDP-43 or FUS proteinopathy associated with ALS and FTLD.

De Novo Design of a ß-Helix Tau Protein Scaffold: An Oligomer-Selective Vaccine Immunogen Candidate for Alzheimer's Disease.

Tau pathology is associated with many neurodegenerative disorders, including Alzheimer's disease (AD), where the spatio-temporal pattern of tau neurofibrillary tangles strongly correlates with disease progression, which motivates therapeutics selective for misfolded tau. Here, we introduce a new avidity-enhanced, multi-epitope approach for protein-misfolding immunogen design, which is predicted to mimic the conformational state of an exposed epitope in toxic tau oligomers. A predicted oligomer-selective tau epitope 343KLDFK347 was scaffolded by designing a ß-helix structure that incorporated multiple instances of the 16-residue tau fragment 339VKSEKLDFKDRVQSKI354. Large-scale conformational ensemble analyses involving Jensen-Shannon Divergence and the embedding depth D showed that the multi-epitope scaffolding approach, employed in designing the ß-helix scaffold, was predicted to better discriminate toxic tau oligomers than other "monovalent" strategies utilizing a single instance of an epitope for vaccine immunogen design. Using Rosetta, 10,000 sequences were designed and screened for the linker portions of the ß-helix scaffold, along with a C-terminal stabilizing α-helix that interacts with the linkers, to optimize the folded structure and stability of the scaffold. Structures were ranked by energy, and the lowest 1% (82 unique sequences) were verified using AlphaFold. Several selection criteria involving AlphaFold are implemented to obtain a lead-designed sequence. The structure was further predicted to have free energetic stability by using Hamiltonian replica exchange molecular dynamics (MD) simulations. The synthesized ß-helix scaffold showed direct binding in surface plasmon resonance (SPR) experiments to several antibodies that were raised to the structured epitope using a designed cyclic peptide. Moreover, the strength of binding of these antibodies to in vitro tau oligomers correlated with the strength of binding to the ß-helix construct, suggesting that the construct presents an oligomer-like conformation and may thus constitute an effective oligomer-selective immunogen.

PROTHON: A Local Order Parameter-Based Method for Efficient Comparison of Protein Ensembles.

The comparison of protein conformational ensembles is of central importance in structural biology. However, there are few computational methods for ensemble comparison, and those that are readily available, such as ENCORE, utilize methods that are sufficiently computationally expensive to be prohibitive for large ensembles. Here, a new method is presented for efficient representation and comparison of protein conformational ensembles. The method is based on the representation of a protein ensemble as a vector of probability distribution functions (pdfs), with each pdf representing the distribution of a local structural property such as the number of contacts between Cß atoms. Dissimilarity between two conformational ensembles is quantified by the Jensen-Shannon distance between the corresponding set of probability distribution functions. The method is validated for conformational ensembles generated by molecular dynamics simulations of ubiquitin, as well as experimentally derived conformational ensembles of a 130 amino acid truncated form of human tau protein. In the ubiquitin ensemble data set, the method was up to 88 times faster than the existing ENCORE software, while simultaneously utilizing 48 times fewer computing cores. We make the method available as a Python package, called PROTHON, and provide a GitHub page with the Python source code at https://github.com/PlotkinLab/Prothon.

Ensemble Generation for Linear and Cyclic Peptides Using a Reservoir Replica Exchange Molecular Dynamics Implementation in GROMACS.

The profile of shapes presented by a cyclic peptide modulates its therapeutic efficacy and is represented by the ensemble of its sampled conformations. Although some algorithms excel at creating a diverse ensemble of cyclic peptide conformations, they seldom address the entropic contribution of flexible conformations and often have significant practical difficulty producing an ensemble with converged and reliable thermodynamic properties. In this study, an accelerated molecular dynamics (MD) method, namely, reservoir replica exchange MD (R-REMD or Res-REMD), was implemented in GROMACS ver. 4.6.7 and benchmarked on two small cyclic peptide model systems: a cyclized furin cleavage site of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (cyclo-(CGPRRARSG)) and oxytocin (disulfide-bonded CYIQNCPLG). Additionally, we also benchmarked Res-REMD on alanine dipeptide and Trpzip2 to demonstrate its validity and efficiency over REMD. For Trpzip2, Res-REMD coupled with an umbrella-sampling-derived reservoir generated similar folded fractions as regular REMD but on a much faster time scale. For cyclic peptides, Res-REMD appeared to be marginally faster than REMD in ensemble generation. Finally, Res-REMD was more effective in sampling rare events such as trans to cis peptide bond isomerization. We provide a GitHub page with the modified GROMACS source code for running Res-REMD at https://github.com/PlotkinLab/Reservoir-REMD.

Rational Generation of Monoclonal Antibodies Selective for Pathogenic Forms of Alpha-Synuclein.

Misfolded toxic forms of alpha-synuclein (α-Syn) have been implicated in the pathogenesis of synucleinopathies, including Parkinson's disease (PD), dementia with Lewy bodies (DLB), and multiple system atrophy (MSA). The α-Syn oligomers and soluble fibrils have been shown to mediate neurotoxicity and cell-to-cell propagation of pathology. To generate antibodies capable of selectively targeting pathogenic forms of α-Syn, computational modeling was used to predict conformational epitopes likely to become exposed on oligomers and small soluble fibrils, but not on monomers or fully formed insoluble fibrils. Cyclic peptide scaffolds reproducing these conformational epitopes exhibited neurotoxicity and seeding activity, indicating their biological relevance. Immunization with the conformational epitopes gave rise to monoclonal antibodies (mAbs) with the desired binding profile showing selectivity for toxic α-Syn oligomers and soluble fibrils, with little or no reactivity with monomers, physiologic tetramers, or Lewy bodies. Recognition of naturally occurring soluble α-Syn aggregates in brain extracts from DLB and MSA patients was confirmed by surface plasmon resonance (SPR). In addition, the mAbs inhibited the seeding activity of sonicated pre-formed fibrils (PFFs) in a thioflavin-T fluorescence-based aggregation assay. In neuronal cultures, the mAbs protected primary rat neurons from toxic α-Syn oligomers, reduced the uptake of PFFs, and inhibited the induction of pathogenic phosphorylated aggregates of endogenous α-Syn. Protective antibodies selective for pathogenic species of α-Syn, as opposed to pan α-Syn reactivity, are expected to provide enhanced safety and therapeutic potency by preserving normal α-Syn function and minimizing the diversion of active antibody from the target by the more abundant non-toxic forms of α-Syn in the circulation and central nervous system.

Optimizing Epitope Conformational Ensembles Using α-Synuclein Cyclic Peptide "Glycindel" Scaffolds: A Customized Immunogen Method for Generating Oligomer-Selective Antibodies for Parkinson's Disease.

Effectively presenting epitopes on immunogens, in order to raise conformationally selective antibodies through active immunization, is a central problem in treating protein misfolding diseases, particularly neurodegenerative diseases such as Alzheimer's disease or Parkinson's disease. We seek to selectively target conformations enriched in toxic, oligomeric propagating species while sparing the healthy forms of the protein that are often more abundant. To this end, we computationally modeled scaffolded epitopes in cyclic peptides by inserting/deleting a variable number of flanking glycines ("glycindels") to best mimic a misfolding-specific conformation of an epitope of α-synuclein enriched in the oligomer ensemble, as characterized by a region most readily disordered and solvent-exposed in a stressed, partially denatured protofibril. We screen and rank the cyclic peptide scaffolds of α-synuclein in silico based on their ensemble overlap properties with the fibril, oligomer-model and isolated monomer ensembles. We present experimental data of seeded aggregation that support nucleation rates consistent with computationally predicted cyclic peptide conformational similarity. We also introduce a method for screening against structured off-pathway targets in the human proteome by selecting scaffolds with minimal conformational similarity between their epitope and the same solvent-exposed primary sequence in structured human proteins. Different cyclic peptide scaffolds with variable numbers of glycines are predicted computationally to have markedly different conformational ensembles. Ensemble comparison and overlap were quantified by the Jensen-Shannon divergence and a new measure introduced here, the embedding depth, which determines the extent to which a given ensemble is subsumed by another ensemble and which may be a more useful measure in developing immunogens that confer conformational selectivity to an antibody.

First Principles Calculation of Protein-Protein Dimer Affinities of ALS-Associated SOD1 Mutants.

Cu,Zn superoxide dismutase (SOD1) is a 32 kDa homodimer that converts toxic oxygen radicals in neurons to less harmful species. The dimerization of SOD1 is essential to the stability of the protein. Monomerization increases the likelihood of SOD1 misfolding into conformations associated with aggregation, cellular toxicity, and neuronal death in familial amyotrophic lateral sclerosis (fALS). The ubiquity of disease-associated mutations throughout the primary sequence of SOD1 suggests an important role of physicochemical processes, including monomerization of SOD1, in the pathology of the disease. Herein, we use a first-principles statistical mechanics method to systematically calculate the free energy of dimer binding for SOD1 using molecular dynamics, which involves sequentially computing conformational, orientational, and separation distance contributions to the binding free energy. We consider the effects of two ALS-associated mutations in SOD1 protein on dimer stability, A4V and D101N, as well as the role of metal binding and disulfide bond formation. We find that the penalty for dimer formation arising from the conformational entropy of disordered loops in SOD1 is significantly larger than that for other protein-protein interactions previously considered. In the case of the disulfide-reduced protein, this leads to a bound complex whose formation is energetically disfavored. Somewhat surprisingly, the loop free energy penalty upon dimerization is still significant for the holoprotein, despite the increased structural order induced by the bound metal cations. This resulted in a surprisingly modest increase in dimer binding free energy of only about 1.5 kcal/mol upon metalation of the protein, suggesting that the most significant stabilizing effects of metalation are on folding stability rather than dimer binding stability. The mutant A4V has an unstable dimer due to weakened monomer-monomer interactions, which are manifested in the calculation by a separation free energy surface with a lower barrier. The mutant D101N has a stable dimer partially due to an unusually rigid ß-barrel in the free monomer. D101N also exhibits anticooperativity in loop folding upon dimerization. These computational calculations are, to our knowledge, the most quantitatively accurate calculations of dimer binding stability in SOD1 to date.