SOD1 exhibits allosteric frustration to facilitate metal binding affinity

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.

Structure of a superoxide dismutase from a tardigrade: Ramazzottius varieornatus strain YOKOZUNA-1

Superoxide dismutase (SOD) is an essential and ubiquitous antioxidant protein that is widely present in biological systems. The anhydrobiotic tardigrades are some of the toughest micro-animals. They have an expanded set of genes for antioxidant proteins such as SODs. These proteins are thought to play an essential role in oxidative stress resistance in critical situations such as desiccation, although their functions at the molecular level have yet to be explored. Here, crystal structures of a copper/zinc-containing SOD (RvSOD15) from an anhydrobiotic tardigrade, Ramazzottius varieornatus strain YOKOZUNA-1, are reported. In RvSOD15, one of the histidine ligands of the catalytic copper center is replaced by a valine (Val87). The crystal structures of the wild type and the V87H mutant show that even though a histidine is placed at position 87, a nearby flexible loop can destabilize the coordination of His87 to the Cu atom. Model structures of other RvSODs were investigated and it was found that some of them are also unusual SODs, with features such as deletion of the electrostatic loop or β3 sheet and unusual metal-binding residues. These studies show that RvSOD15 and some other RvSODs may have evolved to lose the SOD function, suggesting that gene duplications of antioxidant proteins do not solely explain the high stress tolerance of anhydrobiotic tardigrades.

Myricetin: A Potent Anti-Amyloidogenic Polyphenol against Superoxide Dismutase 1 Aggregation

Amyotrophic lateral sclerosis (ALS) is believed to be caused by the aggregation of misfolded or mutated superoxide dismutase 1 (SOD1). As there is currently no treatment, research into aggregation inhibitors continues. Based on docking, molecular dynamics (MD) simulations, and experimental observations, we propose that myricetin, a plant flavonoid, can act as a potent anti-amyloidogenic polyphenol against SOD1 aggregation. Our MD simulation results showed that myricetin stabilizes the protein interface, destabilizes the preformed fibril, and decreases the rate of fibril elongation. Myricetin inhibits the aggregation of SOD1 in a dose-dependent manner as shown by the ThT aggregation kinetics curves. Our transmission electron microscopy, dynamic light scattering, and circular dichroism experiments indicate that fewer shorter fibrils have formed. Fluorescence spectroscopy results predict the involvement of a static quenching mechanism characterized by a strong binding between protein and myricetin. Importantly, size exclusion chromatography revealed the potential of myricetin for fibril destabilization and depolymerization. These experimental observations complement the MD results. Thus, myricetin is a potent SOD1 aggregation inhibitor that can reduce the fibril load. Using the structure of myricetin as a reference, it is possible to design more effective therapeutic inhibitors against ALS that prevent the disease and reverse its effects.

Governing dynamics and preferential binding of the AXH domain influence the aggregation pathway of Ataxin-1

The present state of understanding the mechanism of Spinocerebellar Ataxia-1, a fatal neurodegenerative disease linked to the protein Ataxin-1 (ATXN1), is baffled by a set of self-contradictory, and hence, inconclusive observations. This fallacy poses a bottleneck to the effective designing of curable drugs as the field is currently missing the specific druggable site. To understand the fundamentals of pathogenesis, we tried to decipher the intricacies of the extremely complicated landscape by targeting the relevant species that supposedly dictate the structure-function paradigm. The atomic-level description and characterization of the dynamism of the systems reveal the existence of structural polymorphism in all the leading stakeholders of the overall system. The very existence of conformational heterogeneity in every species creates numerous possible combinations of favorable interactions because of the variability in segmental cross-talks and hence claims its role in the choice of routes between functional activity and dysfunctional disease-causing aggregation. Despite this emergent configurational diversity, there is a common mode of operative intermolecular forces that dictates the extent of stability of all the multimeric complexes due to the localized population of a specific type of residue. The present research proposes a dynamic switch mechanism between aggregability and functional activity, based on the logical interpretation of the estimated variables, which is practically dictated by the effective concentration of the interacting species involved in the cell.

Exploring the Folding Mechanism of Dimeric Superoxide Dismutase

The Cu/Zn Human Superoxide Dismutase (SOD1) is a dimeric metalloenzyme whose genetic mutations are directly related to amyotrophic lateral sclerosis (ALS), so understanding its folding mechanism is of fundamental importance. Currently, the SOD1 dimer formation is studied via molecular dynamics simulations using a simplified structure-based model and an all-atom model. Results from the simplified model reveal a mechanism dependent on distances between monomers, which are limited by constraints to mimic concentration dependence. The stability of intermediates (during the int state) is significantly affected by this distance, as well as by the presence of two folded monomers prior to dimer formation. The kinetics of interface formation are also highly dependent on the separation distance. The folding temperature of the dimer is about 4.2% higher than that of the monomer, a value not too different from experimental data. All-atom simulations on the apo dimer give binding free energy between monomers similar to experimental values. An intermediate state is evident for the apo form at a separation distance between monomers slightly larger than the native distance which has little formed interface between monomers. We have shown that this intermediate is stabilized by non-native intra- and intercontacts.

Amyotrophic lateral sclerosis disease-related mutations disrupt the dimerization of superoxide dismutase 1 - A comparative molecular dynamics simulation study

More than 150 genes are involved in amyotrophic lateral sclerosis (ALS), with superoxide dismutase 1 (SOD1) being one of the most studied. Mutations in SOD1 gene, which encodes the enzyme SOD1 is the second most prevalent and studied cause of familial ALS. SOD1 is a ubiquitous, homodimeric metalloenzyme that forms a critical component of the cellular defense against reactive oxygen species. Several mutations in the SOD1 enzyme cause misfolding, dimerization instability, and increased aggregate formation in ALS. However, there is a lack of information on the dimerization of SOD1 monomers and the mechanistic underpinnings on how the pathogenic mutations disrupt the dimerization mechanism. Here, we presented microsecond-scale molecular dynamics (MD) simulations to unravel how interface-based mutations compromise SOD1 dimerization and provide mechanistic understanding into the corresponding process using WT and three interface-based mutant systems (A4V, T54R, and I113T). Structural stability analysis showed that the mutant systems displayed disparate variations in the catalytic sites which may directly alter the stability and activity of the SOD1 enzyme. Based on the dynamic network analysis and principal component analysis, it has been identified that the mutations weakened the correlated motions along the dimer interface and altered the protein conformational behavior, thus weakening the stability of dimer formation. Moreover, the simulation results identified crucial residues such as G51, D52, G114, I151, and Q153 in establishing the dimerization interaction network, which were weakened or absent in the presence of interfacial mutants. Surface potential analysis on mutant systems also displayed changes in the dimerization potential, thus showing the unfavorable dimer formation. Furthermore, network analysis identified the hotspot residues necessary for SOD1 signal transduction which were surprisingly found in the catalytic sites rather than the anticipated dimerization interface.

Effects of Protein Crowders and Charge on the Folding of Superoxide Dismutase 1 Variants: A Computational Study

The neurodegenerative disease amyotrophic lateral sclerosis (ALS) is associated with the misfolding and aggregation of the metalloenzyme protein superoxide dismutase 1 (SOD1) via mutations that destabilize the monomer-dimer interface. In a cellular environment, crowding and electrostatic screening play essential roles in the folding and aggregation of the SOD1 monomers. Despite numerous studies on the effects of mutations on SOD1 folding, a clear understanding of the interplay between crowding, folding, and aggregation in vivo remains lacking. Using a structure-based minimal model for molecular dynamics simulations, we investigate the role of self-crowding and charge on the folding stability of SOD1 and the G41D mutant where experimentalists were intrigued by an alteration of the folding mechanism by a single point mutation from glycine to charged aspartic acid. We show that unfolded SOD1 configurations are significantly affected by charge and crowding, a finding that would be extremely costly to achieve with all-atom simulations, while the native state is not significantly altered. The mutation at residue 41 alters the interactions between proteins in the unfolded states instead of those within a protein. This paper suggests electrostatics may play an important role in the folding pathway of SOD1 and modifying the charge via mutation and ion concentration may change the dominant interactions between proteins, with potential impacts for aggregation of the mutants. This work provides a plausible reason for the alteration of the unfolded states to address why the mutant G41D causes the changes to the folding mechanism of SOD1 that have intrigued experimentalists.

Misfolding-Associated Exposure of Natively Buried Residues in Mutant SOD1 Facilitates Binding to TRAF6

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease primarily impacting motor neurons. Mutations in superoxide dismutase 1 (SOD1) are the second most common cause of familial ALS. Several of these mutations lead to misfolding or toxic gain of function in the SOD1 protein. Recently, we reported that misfolded SOD1 interacts with TNF receptor-associated factor 6 (TRAF6) in the SOD1^G93A rat model of ALS. Further, we showed in cultured cells that several mutant SOD1 proteins, but not wildtype SOD1 protein, interact with TRAF6 via the MATH domain. Here, we sought to uncover the structural details of this interaction through molecular dynamics (MD) simulations of a dimeric model system, coarse grained using the AWSEM force field. We used direct MD simulations to identify buried residues, and predict binding poses by clustering frames from the trajectories. Metadynamics simulations were also used to deduce preferred binding regions on the protein surfaces from the potential of the mean force in orientation space. Well-folded SOD1 was found to bind TRAF6 via co-option of its native homodimer interface. However, if loops IV and VII of SOD1 were disordered, as typically occurs in the absence of stabilizing Zn²⁺ ion binding, these disordered loops now participated in novel interactions with TRAF6. On TRAF6, multiple interaction hot-spots were distributed around the equatorial region of the MATH domain beta barrel. Expression of TRAF6 variants with mutations in this region in cultured cells demonstrated that TRAF6^T475 facilitates interaction with different SOD1 mutants. These findings contribute to our understanding of the disease mechanism and uncover potential targets for the development of therapeutics.

Probing the polyphenolic flavonoid, morin as a highly efficacious inhibitor against amyloid(A4V) mutant SOD1 in fatal amyotrophic lateral sclerosis

Deposition of misfolded protein aggregates in key areas of human brain is the quintessential trait of various pertinent neurodegenerative disorders including amyotrophic lateral sclerosis (ALS). Genetic point mutations in Cu/Zn superoxide dismutase (SOD1) are found to be the most important contributing factor behind familial ALS. Especially, single nucleotide polymorphism (SNP) A4V is the most nocuous since it substantially decreases life expectancy of patients. Besides, the use of naturally occurring polyphenolic flavonoids is profoundly being advocated for palliating amyloidogenic behavior of proteopathic proteins. In the present analysis, through proficient computational tools, we have attempted to ascertain a pharmacodynamically promising flavonoid compound that effectively curbs the pathogenic behavior of A4V SOD1 mutant. Initial screening of flavonoids that exhibit potency against amyloids identified morin, myricetin and epigallocatechin gallate as promising leads. Further, with the help of feasible and yet adept protein-ligand interaction studies and stalwart molecular simulation analyses, we were able to observe that aforementioned flavonoids were able to considerably divert mutant A4V SOD1 from its distinct pathogenic behavior. Among which, morin showed the most curative potential against A4V SOD1. Therefore, morin holds a great therapeutic potential in contriving highly efficacious inhibitors in mitigating fatal and insuperable ALS.

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.

Molecular interactions between C9ORF72 and SMCR8: A local energetic frustration perspective

The hexanucleotide repeat expansion in C9orf72 represents a major cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). C9orf72, together with SMCR8 and WDR41, can form a stable complex that regulates autophagy and membrane trafficking. Very recently, the cryo-EM structure of C9orf72-SMCR8-WDR41 helps in understanding the structure-function relationship of C9orf72. This protein complex is indispensable to several cellular processes and is strongly linked to familial ALS and FTD. Understanding the molecular basis of the C9orf72-SMCR8 protein-protein interaction is thus important to comprehend their function. To establish a basis for understanding the relationships between sequence, structure, and function of the C9orf72, this study reports a local frustration analysis on the C9orf72-SMCR8 complex structure. An analysis of local frustration profiles indicated that (1) the structural domains in C9orf72 are minimally-frustrated and relatively conserved, (2) high frustration patches on the protein-protein interface (3) increased frustration in the C-terminal helices involved in the dimerization of C9orf72 structures.

Free energy calculations of ALS-causing SOD1 mutants reveal common perturbations to stability and dynamics along the maturation pathway

With over 150 heritable mutations identified as disease-causative, superoxide dismutase 1 (SOD1) has been a main target of amyotrophic lateral sclerosis (ALS) research and therapeutic efforts. However, recent evidence has suggested that neither loss of function nor protein aggregation is responsible for promoting neurotoxicity. Furthermore, there is no clear pattern to the nature or the location of these mutations that could suggest a molecular mechanism behind SOD1-linked ALS. Here, we utilize reliable and accurate computational techniques to predict the perturbations of 10 such mutations to the free energy changes of SOD1 as it matures from apo monomer to metallated dimer. We find that the free energy perturbations caused by these mutations strongly depend on maturational progress, indicating the need for state-specific therapeutic targeting. We also find that many mutations exhibit similar patterns of perturbation to native and non-native maturation, indicating strong thermodynamic coupling between the dynamics at various sites of maturation within SOD1. These results suggest the presence of an allosteric network in SOD1 which is vulnerable to disruption by these mutations. Analysis of these perturbations may contribute to uncovering a unifying molecular mechanism which explains SOD1-linked ALS and help to guide future therapeutic efforts.

A pH Switch Controls Zinc Binding in Tomato Copper-Zinc Superoxide Dismutase

Copper-zinc superoxide dismutase (SOD1) is a major antioxidant metalloenzyme that protects cells from oxidative damage by superoxide anions (O2-). Structural, biophysical, and other characteristics have in the past been compiled for mammalian SOD1s and for the highly homologous fungal and bovine SOD1s. Here, we characterize the biophysical properties of a plant SOD1 from tomato chloroplasts and present several of its crystal structures. The most unusual of these structures is a structure at low pH in which tSOD1 harbors zinc in the copper-binding site but contains no metal in the zinc-binding site. The side chain of D83, normally a zinc ligand, adopts an alternate rotameric conformation to form an unusual bidentate hydrogen bond with the side chain of D124, precluding metal binding in the zinc-binding site. This alternate conformation of D83 appears to be responsible for the previously observed pH-dependent loss of zinc from the zinc-binding site of SOD1. Titrations of cobalt into apo tSOD1 at a similar pH support the lack of an intact zinc-binding site. Further characterization of tSOD1 reveals that it is a weaker dimer relative to human SOD1 and that it can be activated in vivo through a copper chaperone for the SOD1-independent mechanism.

Structural fluctuations and mechanical stabilities of the metamorphic protein RfaH

RfaH is a compact two-domain bacterial transcription factor that functions both as a regulator of transcription and an enhancer of translation. Underpinning the dual functional roles of RfaH is a partial but dramatic fold switch, which completely transforms the ~50-amino acid C-terminal domain (CTD) from an all-α state to an all-β state. The fold switch of the CTD occurs when RfaH binds to RNA polymerase (RNAP), however, the details of how this structural transformation is triggered is not well understood. Here we use all-atom Monte Carlo simulations to characterize structural fluctuations and mechanical stability properties of the full-length RfaH and the CTD as an isolated fragment. In agreement with experiments, we find that interdomain contacts are crucial for maintaining a stable, all-α CTD in free RfaH. To probe mechanical properties, we use pulling simulations to measure the work required to inflict local deformations at different positions along the chain. The resulting mechanical stability profile reveals that free RfaH can be divided into a "rigid" part and a "soft" part, with a boundary that nearly coincides with the boundary between the two domains. We discuss the potential role of this feature for how fold switching may be triggered by interaction with RNAP.

Conformational dynamics of superoxide dismutase (SOD1) in osmolytes: a molecular dynamics simulation study

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease caused by the misfolding of Cu, Zn superoxide dismutase (SOD1). Several earlier studies have shown that monomeric apo SOD1 undergoes significant local unfolding dynamics and is the predecessor for aggregation. Here, we have employed atomistic molecular dynamics (MD) simulations to study the structure and dynamics of monomeric apo and holo SOD1 in water, aqueous urea and aqueous urea-TMAO (trimethylamine oxide) solutions. Loop IV (zinc-binding loop) and loop VII (electrostatic loop) of holo SOD1 are considered as functionally important loops as they are responsible for the structural stability of holo SOD1. We found larger local unfolding of loop IV and VII of apo SOD1 as compared to holo SOD1 in water. Urea induced more unfolding in holo SOD1 than apo SOD1, whereas the stabilization of both the form of SOD1 was observed in ternary solution (i.e. water/urea/TMAO solution) but the extent of stabilization was higher in holo SOD1 than apo SOD1. The partially unfolded structures of apo SOD1 in water, urea and holo SOD1 in urea were identified by the exposure of the hydrophobic cores, which are highly dynamic and these may be the initial events of aggregation in SOD1. Our simulation studies support the formation of aggregates by means of the local unfolding of monomeric apo SOD1 as compared to holo SOD1 in water.

Local frustration determines loop opening during the catalytic cycle of an oxidoreductase

Local structural frustration, the existence of mutually exclusive competing interactions, may explain why some proteins are dynamic while others are rigid. Frustration is thought to underpin biomolecular recognition and the flexibility of protein-binding sites. Here, we show how a small chemical modification, the oxidation of two cysteine thiols to a disulfide bond, during the catalytic cycle of the N-terminal domain of the key bacterial oxidoreductase DsbD (nDsbD), introduces frustration ultimately influencing protein function. In oxidized nDsbD, local frustration disrupts the packing of the protective cap-loop region against the active site allowing loop opening. By contrast, in reduced nDsbD the cap loop is rigid, always protecting the active-site thiols from the oxidizing environment of the periplasm. Our results point toward an intricate coupling between the dynamics of the active-site cysteines and of the cap loop which modulates the association reactions of nDsbD with its partners resulting in optimized protein function.

Role of substrate-product frustration on enzyme functional dynamics

Natural enzymes often have enormous catalytic power developed by evolution. Revealing the underlying physical strategy used by enzymes to achieve high catalysis efficiency is one of the central focuses in the field of biological physics. Our recent work demonstrated that multisubstrate enzymes can utilize steric frustration encountered in the substrate-product cobound complex to overcome the bottleneck of the enzymatic cycle [W. Li et al., Phys. Rev. Lett. 122, 238102 (2019)10.1103/PhysRevLett.122.238102]. However, the key atomic-level interactions by which the steric frustration contributes to the enzymatic cycle remain elusive. In this work we study the microscopic mechanism for the role of the substrate-product frustration on the key physical steps in the enzymatic cycle of adenylate kinase (AdK), a multisubstrate enzyme catalyzing the reversible phosphoryl transfer reaction ATP+AMP⇋ADP+ADP. By using atomistic molecular dynamics simulations with enhanced sampling, we showed that the competitive interactions from the phosphate groups of the substrate ATP and product ADP in the ATP-ADP cobound complex of the AdK lead to local frustration in the binding pockets. Such local frustration disrupts the hydrogen bond network around the binding pockets, which causes lowered barrier height for the opening of the enzyme conformations and expedited release of the bottleneck product ADP. Our results directly demonstrated from the atomistic level that the local frustration in the active sites of the enzyme can be utilized to facilitate the key physical steps of the enzymatic cycle, providing numerical evidence to the predictions of the previous theoretical work.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

Stability and Local Unfolding of SOD1 in the Presence of Protein Crowders

Using NMR and Monte Carlo (MC) methods, we investigate the stability and dynamics of superoxide dismutase 1 (SOD1) in homogeneous crowding environments, where either bovine pancreatic trypsin inhibitor (BPTI) or the B1 domain of streptococcal protein G (PGB1) serves as a crowding agent. By NMR, we show that both crowders, and especially BPTI, cause a drastic loss in the overall stability of SOD1 in its apo monomeric form. Additionally, we determine chemical shift perturbations indicating that SOD1 interacts with the crowder proteins in a residue-specific manner that further depends on the identity of the crowding protein. Furthermore, the specificity of SOD1-crowder interactions is reciprocal: chemical shift perturbations on BPTI and PGB1 identify regions that interact preferentially with SOD1. By MC simulations, we investigate the local unfolding of SOD1 in the absence and presence of the crowders. We find that the crowders primarily interact with the long flexible loops of the folded SOD1 monomer. The basic mechanisms by which the SOD1 β-barrel core unfolds remain unchanged when adding the crowders. In particular, both with and without the crowders, the second β-sheet of the barrel is more dynamic and unfolding-prone than the first. Notably, the MC simulations (exploring the early stages of SOD1 unfolding) and the NMR experiments (under equilibrium conditions) identify largely the same set of PGB1 and BPTI residues as prone to form SOD1 contacts. Thus, contacts stabilizing the unfolded state of SOD1 in many cases appear to form early in the unfolding reaction.

Local frustration around enzyme active sites

Conflicting biological goals often meet in the specification of protein sequences for structure and function. Overall, strong energetic conflicts are minimized in folded native states according to the principle of minimal frustration, so that a sequence can spontaneously fold, but local violations of this principle open up the possibility to encode the complex energy landscapes that are required for active biological functions. We survey the local energetic frustration patterns of all protein enzymes with known structures and experimentally annotated catalytic residues. In agreement with previous hypotheses, the catalytic sites themselves are often highly frustrated regardless of the protein oligomeric state, overall topology, and enzymatic class. At the same time a secondary shell of more weakly frustrated interactions surrounds the catalytic site itself. We evaluate the conservation of these energetic signatures in various family members of major enzyme classes, showing that local frustration is evolutionarily more conserved than the primary structure itself.

Tryptophan 32 mediates SOD1 toxicity in a in vivo motor neuron model of ALS and is a promising target for small molecule therapeutics

SOD1 misfolding, toxic gain of function, and spread are proposed as a pathological basis of amyotrophic lateral sclerosis (ALS), but the nature of SOD1 toxicity has been difficult to elucidate. Uniquely in SOD1 proteins from humans and other primates, and rarely in other species, a tryptophan residue at position 32 (W32) is predicted to be solvent exposed and to participate in SOD1 misfolding. We hypothesized that W32 is influential in SOD1 acquiring toxicity, as it is known to be important in template-directed misfolding. We tested if W32 contributes to SOD1 cytotoxicity and if it is an appropriate drug target to ameliorate ALS-like neuromuscular deficits in a zebrafish model of motor neuron axon morphology and function (swimming). Embryos injected with human SOD1 variant with W32 substituted for a serine (SOD1^W32S) had reduced motor neuron axonopathy and motor deficits compared to those injected with wildtype or disease-associated SOD1. A library of FDA-approved small molecules was ranked with virtual screening based on predicted binding to W32, and subsequently filtered for analogues using a pharmacophore model based on molecular features of the uracil moiety of a small molecule previously predicted to interact with W32 (5'-fluorouridine or 5'-FUrd). Along with testing 5'-FUrd and uridine, a lead candidate from this list was selected based on its lower toxicity and improved blood brain barrier penetrance; telbivudine significantly rescued SOD1 toxicity in a dose-dependent manner. The mechanisms whereby the small molecules ameliorated motor neuron phenotypes were specifically mediated through human SOD1 and its residue W32, because these therapeutics had no measurable impact on the effects of UBQLN4^D90A, EtOH, or tryptophan-deficient human SOD1^W32S. By substituting W32 for a more evolutionarily conserved residue (serine), we confirmed the significant influence of W32 on human SOD1 toxicity to motor neuron morphology and function; further, we performed pharmaceutical targeting of the W32 residue for rescuing SOD1 toxicity. This unique residue offers future novel insights into SOD1 stability and toxic gain of function, and therefore poses an potential target for drug therapy.

TFE-induced local unfolding and fibrillation of SOD1: bridging the experiment and simulation studies

Misfolding and aggregation of Cu, Zn Superoxide dismutase (SOD1) is involved in the neurodegenerative disease, amyotrophic lateral sclerosis. Many studies have shown that metal-depleted, monomeric form of SOD1 displays substantial local unfolding dynamics and is the precursor for aggregation. Here, we have studied the structure and dynamics of different apo monomeric SOD1 variants associated with unfolding and aggregation in aqueous trifluoroethanol (TFE) through experiments and simulation. TFE induces partially unfolded β-sheet-rich extended conformations in these SOD1 variants, which subsequently develops aggregates with fibril-like characteristics. Fibrillation was achieved more easily in disulfide-reduced monomeric SOD1 when compared with wild-type and mutant monomeric SOD1. At higher concentrations of TFE, a native-like structure with the increase in α-helical content was observed. The molecular dynamics simulation results illustrate distinct structural dynamics for different regions of SOD1 variants and show uniform local unfolding of β-strands. The strands protected by the zinc-binding and electrostatic loops were found to unfold first in 20% (v/v) TFE, leading to a partial unfolding of β-strands 4, 5, and 6 which are prone to aggregation. Our results thus shed light on the role of local unfolding and conformational dynamics in SOD1 misfolding and aggregation.

Phylogenetic spread of sequence data affects fitness of SOD1 consensus enzymes: Insights from sequence statistics and structural analyses

Non-natural protein sequences with native-like structures and functions can be constructed successfully using consensus design. This design strategy is relatively well understood in repeat proteins with simple binding function, however detailed studies are lacking in globular enzymes. The SOD1 family is a good model for such studies due to the availability of large amount of sequence and structure data motivated by involvement of human SOD1 in the fatal motor neuron disease amyotrophic lateral sclerosis (ALS). We constructed two consensus SOD1 enzymes from multiple sequence alignments from all organisms and eukaryotic organisms. A significant difference in their catalytic activities shows that the phylogenetic spread of the sequences used affects the fitness of the construct obtained. A mutation in an electrostatic loop and overall design incompatibilities between bacterial and eukaryotic sequences were implicated in this disparity. Based on this analysis, a bioinformatics approach was used to classify mutations thought to cause familial ALS providing a unique high level view of the physical basis of disease-causing aggregation of human SOD1.

The Role of Metal Binding in the Amyotrophic Lateral Sclerosis-Related Aggregation of Copper-Zinc Superoxide Dismutase

Protein misfolding and conformational changes are common hallmarks in many neurodegenerative diseases involving formation and deposition of toxic protein aggregates. Although many players are involved in the in vivo protein aggregation, physiological factors such as labile metal ions within the cellular environment are likely to play a key role. In this review, we elucidate the role of metal binding in the aggregation process of copper-zinc superoxide dismutase (SOD1) associated to amyotrophic lateral sclerosis (ALS). SOD1 is an extremely stable Cu-Zn metalloprotein in which metal binding is crucial for folding, enzymatic activity and maintenance of the native conformation. Indeed, demetalation in SOD1 is known to induce misfolding and aggregation in physiological conditions in vitro suggesting that metal binding could play a key role in the pathological aggregation of SOD1. In addition, this study includes recent advances on the role of aberrant metal coordination in promoting SOD1 aggregation, highlighting the influence of metal ion homeostasis in pathologic aggregation processes.

The unfolding mechanism of monomeric mutant SOD1 by simulated force spectroscopy

Mechanical unfolding of mutated apo, disulfide-reduced, monomeric superoxide dismutase 1 protein (SOD1) has been simulated via force spectroscopy techniques, using both an all-atom (AA), explicit solvent model and a coarse-grained heavy-atom Gō (HA-Gō) model. The HA-Gō model was implemented at two different pulling speeds for comparison. The most-common sequence of unfolding in the AA model agrees well with the most-common unfolding sequence of the HA-Gō model, when the same normalized pulling rate was used. Clustering of partially-native structures as the protein unfolds shows that the AA and HA-Gō models both exhibit a dominant pathway for early unfolding, which eventually bifurcates repeatedly to multiple branches after the protein is about half-unfolded. The force-extension curve exhibits multiple force drops, which are concomitant with jumps in the local interaction potential energy between specific β-strands in the protein. These sudden jumps in the potential energy coincide with the dissociation of specific pairs of β-strands, and thus intermediate unfolding events. The most common sequence of β-strand dissociation in the unfolding pathway of the AA model is β-strands 5, 4, 8, 7, 1, 2, then finally β-strands 3 and 6. The observation that β-strand 5 is among the first to unfold here, but the last to unfold in simulations of loop-truncated SOD1, could imply the existence of an evolutionary compensation mechanism, which would stabilize β-strands flanking long loops against their entropic penalty by strengthening intramolecular interactions. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Incorporation of a Reporter Peptide in FPOP Compensates for Adventitious Scavengers and Permits Time-Dependent Measurements

Incorporation of a reporter peptide in solutions submitted to fast photochemical oxidation of proteins (FPOP) allows for the correction of adventitious scavengers and enables the normalization and comparison of time-dependent results. Reporters will also be useful in differential experiments to control for the inclusion of a radical-reactive species. This incorporation provides a simple and quick check of radical dosage and allows comparison of FPOP results from day-to-day and lab-to-lab. Use of a reporter peptide in the FPOP workflow requires no additional measurements or spectrometers while building a more quantitative FPOP platform. It requires only measurement of the extent of reporter-peptide modification in a LC/MS/MS run, which is performed by using either data-dependent scanning or an inclusion list. Graphical Abstract ᅟ.

As Simple As Possible, but Not Simpler: Exploring the Fidelity of Coarse-Grained Protein Models for Simulated Force Spectroscopy

Mechanical unfolding of a single domain of loop-truncated superoxide dismutase protein has been simulated via force spectroscopy techniques with both all-atom (AA) models and several coarse-grained models having different levels of resolution: A Gō model containing all heavy atoms in the protein (HA-Gō), the associative memory, water mediated, structure and energy model (AWSEM) which has 3 interaction sites per amino acid, and a Gō model containing only one interaction site per amino acid at the Cα position (Cα-Gō). To systematically compare results across models, the scales of time, energy, and force had to be suitably renormalized in each model. Surprisingly, the HA-Gō model gives the softest protein, exhibiting much smaller force peaks than all other models after the above renormalization. Clustering to render a structural taxonomy as the protein unfolds showed that the AA, HA-Gō, and Cα-Gō models exhibit a single pathway for early unfolding, which eventually bifurcates repeatedly to multiple branches only after the protein is about half-unfolded. The AWSEM model shows a single dominant unfolding pathway over the whole range of unfolding, in contrast to all other models. TM alignment, clustering analysis, and native contact maps show that the AWSEM pathway has however the most structural similarity to the AA model at high nativeness, but the least structural similarity to the AA model at low nativeness. In comparison to the AA model, the sequence of native contact breakage is best predicted by the HA-Gō model. All models consistently predict a similar unfolding mechanism for early force-induced unfolding events, but diverge in their predictions for late stage unfolding events when the protein is more significantly disordered.

Superoxide Dismutase 1 (SOD1)-Derived Peptide Inhibits Amyloid Aggregation of Familial Amyotrophic Lateral Sclerosis SOD1 Mutants

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder that leads to the death of the upper and lower motor neurons. Superoxide dismutase 1 (SOD1) is an ALS pathogenic protein, whose misfolding results in the formation of amyloid aggregates. The mechanism underlying SOD1 pathogenesis in ALS remains obscure, but one possible mechanism involves gain-of-interaction, in which the misfolded soluble SOD1 forms abnormal protein-protein interactions (PPIs) with various cellular proteins, including with other SOD1 molecules, thereby interfering with their function. The structural basis of this gain-of-interaction mechanism is unknown. Here, we characterized the backbone dynamics landscape of misfolded SOD1 to pinpoint surface areas predisposed to aberrant PPIs. This analysis enabled us to formulate a working hypothesis for the mechanism of the gain-of-function of misfolded SOD1, according to which an abnormal PPI potential results from the increased mobility of the SOD1 surface backbone. Guided by the backbone dynamics landscape, we have identified a SOD1-derived peptide that can bind SOD1 proteins and divert the typical amyloid aggregation of ALS-related SOD1 mutants toward a potentially less toxic amorphous aggregation pathway.

A mechanism for propagated SOD1 misfolding from frustration analysis of a G85R mutant protein assembly

Application of landscape theory and the dehydron hypothesis to a crystal structure of a G85R mutant superoxide dismutase (SOD1) tetrameric complex allows for the description of a prion-like hypothesis that serves to explain propagated SOD1 misfolding. We have developed two conformational-change scenarios, one local to the ESL at the complex interface, and a second displacement at the ESL of the otherdimeric subunit. When taken together these provide for a prion-like mechanism that can serve to explain the observed conversion of wtSOD1 to a misfolded form by the G85R mutant.

Protein Frustratometer 2: a tool to localize energetic frustration in protein molecules, now with electrostatics

The protein frustratometer is an energy landscape theory-inspired algorithm that aims at localizing and quantifying the energetic frustration present in protein molecules. Frustration is a useful concept for analyzing proteins’ biological behavior. It compares the energy distributions of the native state with respect to structural decoys. The network of minimally frustrated interactions encompasses the folding core of the molecule. Sites of high local frustration often correlate with functional regions such as binding sites and regions involved in allosteric transitions. We present here an upgraded version of a webserver that measures local frustration. The new implementation that allows the inclusion of electrostatic energy terms, important to the interactions with nucleic acids, is significantly faster than the previous version enabling the analysis of large macromolecular complexes within a user-friendly interface. The webserver is freely available at URL: http://frustratometer.qb.fcen.uba.ar.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.

Omega-3 fatty acids differentially modulate enzymatic anti-oxidant systems in skeletal muscle cells

During physical activity, increased reactive oxygen species production occurs, which can lead to cell damage and in a decline of individual's performance and health. The use of omega-3 polyunsaturated fatty acids as a supplement to protect the immune system has been increasing; however, their possible benefit to the anti-oxidant system is not well described. Thus, the aim of this study was to evaluate whether the omega-3 fatty acids (docosahexaenoic acid and eicosapentaenoic acid) can be beneficial to the anti-oxidant system in cultured skeletal muscle cells. C2C12 myocytes were differentiated and treated with either eicosapentaenoic acid or docosahexaenoic acid for 24 h. Superoxide content was quantified using the dihydroethidine oxidation method and superoxide dismutase, catalase, and glutathione peroxidase activity, and expression was quantified. We observed that the docosahexaenoic fatty acids caused an increase in superoxide production. Eicosapentaenoic acid induced catalase activity, while docosahexaenoic acid suppressed superoxide dismutase activity. In addition, we found an increased protein expression of the total manganese superoxide dismutase and catalase enzymes when cells were treated with eicosapentaenoic acid. Taken together, these data indicate that the use of eicosapentaenoic acid may present both acute and chronic benefits; however, the treatment with DHA may not be beneficial to muscle cells.

Molecular interactions between photosystem I and ferredoxin: an integrated energy frustration and experimental model

The stromal domain (PsaC, PsaD, and PsaE) of photosystem I (PSI) reduces transiently bound ferredoxin (Fd) or flavodoxin. Experimental structures exist for all of these protein partners individually, but no experimental structure of the PSI/Fd or PSI/flavodoxin complexes is presently available. Molecular models of Fd docked onto the stromal domain of the cyanobacterial PSI site are constructed here utilizing X-ray and NMR structures of PSI and Fd, respectively. Predictions of potential protein-protein interaction regions are based on experimental site-directed mutagenesis and cross-linking studies to guide rigid body docking calculations of Fd into PSI, complemented by energy landscape theory to bring together regions of high energetic frustration on each of the interacting proteins. The results identify two regions of high localized frustration on the surface of Fd that contain negatively charged Asp and Glu residues. This study predicts that these regions interact predominantly with regions of high localized frustration on the PsaC, PsaD, and PsaE chains of PSI, which include several residues predicted by previous experimental studies.

Insights into the role of the unusual disulfide bond in copper-zinc superoxide dismutase

The functional and structural significance of the intrasubunit disulfide bond in copper-zinc superoxide dismutase (SOD1) was studied by characterizing mutant forms of human SOD1 (hSOD) and yeast SOD1 lacking the disulfide bond. We determined x-ray crystal structures of metal-bound and metal-deficient hC57S SOD1. C57S hSOD1 isolated from yeast contained four zinc ions per protein dimer and was structurally very similar to wild type. The addition of copper to this four-zinc protein gave properly reconstituted 2Cu,2Zn C57S hSOD, and its spectroscopic properties indicated that the coordination geometry of the copper was remarkably similar to that of holo wild type hSOD1. In contrast, the addition of copper and zinc ions to apo C57S human SOD1 failed to give proper reconstitution. Using pulse radiolysis, we determined SOD activities of yeast and human SOD1s lacking disulfide bonds and found that they were enzymatically active at ∼10% of the wild type rate. These results are contrary to earlier reports that the intrasubunit disulfide bonds in SOD1 are essential for SOD activity. Kinetic studies revealed further that the yeast mutant SOD1 had less ionic attraction for superoxide, possibly explaining the lower rates. Saccharomyces cerevisiae cells lacking the sod1 gene do not grow aerobically in the absence of lysine, but expression of C57S SOD1 increased growth to 30-50% of the growth of cells expressing wild type SOD1, supporting that C57S SOD1 retained a significant amount of activity.

Free-radical first responders: the characterization of CuZnSOD and MnSOD regulation during freezing of the freeze-tolerant North American wood frog, Rana sylvatica

BACKGROUND

The North American wood frog, Rana sylvatica, is able to overcome subzero conditions through overwintering in a frozen state. Freezing imposes ischemic and oxidative stress on cells as a result of cessation of blood flow. Superoxide dismutases (SODs) catalyze the redox reaction involving the dismutation of superoxide (O(2)(-)) to molecular oxygen and hydrogen peroxide.

METHODS

The present study investigated the regulation of CuZnSOD and MnSOD kinetics as well as the transcript, protein and phosphorylation levels of purified enzyme from the muscle of control and frozen R. sylvatica.

RESULTS

CuZnSOD from frozen muscle showed a significantly higher V(max) (1.52 fold) in comparison to CuZnSOD from the muscle of control frogs. MnSOD from frozen muscle showed a significantly lower Km for O(2)(-) (0.66 fold) in comparison to CuZnSOD from control frogs. MnSOD from frozen frogs showed higher phosphorylation of serine (2.36 fold) and tyrosine (1.27 fold) residues in comparison to MnSOD from control animals. Susceptibility to digestion via thermolysin after incubation with increasing amount of urea (C(m)) was tested, resulting in no significant changes for CuZnSOD, whereas a significant change in MnSOD stability was observed between control (2.53 M urea) and frozen (2.92 M urea) frogs. Expressions of CuZnSOD and MnSOD were quantified at both mRNA and protein levels in frog muscle, but were not significantly different.

CONCLUSION

The physiological consequence of freeze-induced SOD modification appears to adjust SOD function in freezing frogs.

GENERAL SIGNIFICANCE

Augmented SOD activity may increase the ability of R. sylvatica to overcome oxidative stress associated with ischemia.

Insights into SOD1-linked amyotrophic lateral sclerosis from NMR studies of Ni(2+)- and other metal-ion-substituted wild-type copper-zinc superoxide dismutases

The dimeric Cu-Zn superoxide dismutase (SOD1) is a particularly interesting system for biological inorganic chemical studies because substitutions of the native Cu and/or Zn ions by a nonnative metal ion cause minimal structural changes and result in high enzymatic activity for those derivatives with Cu remaining in the Cu site. The pioneering NMR studies of the magnetically coupled derivative Cu2Co2SOD1 by Ivano Bertini and coworkers are of particular importance in this regard. In addition to Co(2+), Ni(2+) is a versatile metal ion for substitution into SOD1, showing very little disturbance of the structure in Cu2Ni2SOD1 and acting as a very good mimic of the native Cu ion in Ni2Zn2SOD1. The NMR studies presented here were inspired by and are indebted to Ivano Bertini's paramagnetic NMR pursuits of metalloproteins. We report Ni(2+) binding to apo wild-type SOD1 and a time-dependent Ni(2+) migration from the Zn site to the Cu site, and the preparation and characterization of Ni2Ni2SOD1, which shows coordination properties similar to those of Cu2Cu2SOD1, namely, an anion-binding property different from that of the wild type and a possibly broken bridging His. Mutations in the human SOD1 gene can cause familial amyotrophic lateral sclerosis (ALS), and mutant SOD1 proteins with significantly altered metal-binding behaviors are implicated in causing the disease. We conclude by discussing the effects of the ALS mutations on the remarkable stabilities and metal-binding properties of wild-type SOD1 proteins and the implications concerning the causes of SOD1-linked ALS.

Mechanical resistance in unstructured proteins

Single-molecule pulling experiments on unstructured proteins linked to neurodegenerative diseases have measured rupture forces comparable to those for stable folded proteins. To investigate the structural mechanisms of this unexpected force resistance, we perform pulling simulations of the amyloid β-peptide (Aβ) and α-synuclein (αS), starting from simulated conformational ensembles for the free monomers. For both proteins, the simulations yield a set of rupture events that agree well with the experimental data. By analyzing the conformations occurring shortly before rupture in each event, we find that the mechanically resistant structures share a common architecture, with similarities to the folds adopted by Aβ and αS in amyloid fibrils. The disease-linked Arctic mutation of Aβ is found to increase the occurrence of highly force-resistant structures. Our study suggests that the high rupture forces observed in Aβ and αS pulling experiments are caused by structures that might have a key role in amyloid formation.

Local unfolding and aggregation mechanisms of SOD1: a Monte Carlo exploration

Copper, zinc superoxide dismutase 1 (SOD1) is a ubiquitous homodimeric enzyme, whose misfolding and aggregation play a potentially key role in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). SOD1 aggregation is thought to be preceded by dimer dissociation and metal loss, but the mechanisms by which the metal-free monomer aggregates remain incompletely understood. Here we use implicit solvent all-atom Monte Carlo (MC) methods to investigate the local unfolding dynamics of the β-barrel-forming SOD1 monomer. Although event-to-event variations are large, on average, we find clear differences in dynamics among the eight strands forming the β-barrel. Most dynamic is the eighth strand, β8, which is located in the dimer interface of native SOD1. For the four strands in or near the dimer interface (β1, β2, β7, and β8), we perform aggregation simulations to assess the propensity of these chain segments to self-associate. We find that β1 and β2 readily self-associate to form intermolecular parallel β-sheets, whereas β8 shows a very low aggregation propensity.

Hypothesis: The unfolding power of protein dielectricity

A hypothesis is proposed on a potential role of protein dielectricity as an unfolding factor in protein-protein interactions. It is suggested that large protein complexes and aggregation seeds can unfold target proteins by virtue of their effect on the dielectric properties of water at the protein-solvent interface. Here, similar to the effect of membrane surfaces, protein surface can cause decrease in the local dielectric constant of solvent and thereby induce structural changes in a target protein approaching this surface. Some potential implementations of this hypothetical mechanism are also discussed.

Polymer uncrossing and knotting in protein folding, and their role in minimal folding pathways

We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold [Formula: see text], [Formula: see text], [Formula: see text], and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation "alignment". The consensus minimal pathway is constructed and shown schematically for representative cases of an [Formula: see text], [Formula: see text], and knotted protein. An overlap parameter is defined between pathways; we find that [Formula: see text] proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and [Formula: see text] proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.