The protein folding rate and the geometry and topology of the native state

Folds from fold: Exploring topological isoforms of a single-domain protein

Expanding the protein fold space beyond linear chains is of fundamental significance, yet remains largely unexplored. Herein, we report the creation of seven topological isoforms (i.e., linear, cyclic, knot, lasso, pseudorotaxane, and catenane) from a single protein fold precursor by rewiring the connectivity of secondary structure elements of the SpyTag-SpyCatcher complex and mutating the reactive residue on SpyTag to abolish the isopeptide bonding. These topological isoforms can be directly expressed in cells. Their topologies were confirmed by combined techniques of proteolytic digestion, fluorescence correlation spectroscopy (FCS), size-exclusion chromatography (SEC), and topological transformation. To study the effects of topology on their structures and properties, their biophysical properties were characterized by differential scanning calorimetry (DSC), heteronuclear single quantum coherence nuclear magnetic resonance spectroscopy (HSQC-NMR), and circular dichroism (CD) spectroscopy. Molecular dynamics (MD) simulations were further performed to reveal the atomic details of structural changes upon unfolding. Both experimental and simulation results suggest that they share a similar, well-folded hydrophobic core but exhibit distinct folding/unfolding dynamic behaviors. These results shed light onto the folding landscape of topological isoforms derived from the same protein fold. As a model system, this work improves our understanding of protein structure and dynamics beyond linear chains and suggests that protein folds are highly amenable to topological variation.

Mathematical topology and geometry-based classification of tauopathies

Neurodegenerative diseases, like Alzheimer's, are associated with the presence of neurofibrillary lesions formed by tau protein filaments in the cerebral cortex. While it is known that different morphologies of tau filaments characterize different neurodegenerative diseases, there are few metrics of global and local structure complexity that enable to quantify their structural diversity rigorously. In this manuscript, we employ for the first time mathematical topology and geometry to classify neurodegenerative diseases by using cryo-electron microscopy structures of tau filaments that are available in the Protein Data Bank. By employing mathematical topology metrics (Gauss linking integral, writhe and second Vassiliev measure) we achieve a consistent, but more refined classification of tauopathies, than what was previously observed through visual inspection. Our results reveal a hierarchy of classification from global to local topology and geometry characteristics. In particular, we find that tauopathies can be classified with respect to the handedness of their global conformations and the handedness of the relative orientations of their repeats. Progressive supranuclear palsy is identified as an outlier, with a more complex structure than the rest, reflected by a small, but observable knotoid structure (a diagrammatic structure representing non-trivial topology). This topological characteristic can be attributed to a pattern in the beginning of the R3 repeat that is present in all tauopathies but at different extent. Moreover, by comparing single filament to paired filament structures within tauopathies we find a consistent change in the side-chain orientations with respect to the alpha carbon atoms at the area of interaction.

Determining suitable surfactant concentration ranges to avoid protein unfolding in pharmaceutical formulations using UV analysis

Protein stability is fundamental to maintain pharmaceutical efficacy in the nascent field of biologics. One particular property that is essential for therapeutic effect is retention of the folded 3-dimensional conformation, i.e. once unfolding has occurred the biologic is often rendered inactive. In this work we propose a modified form of a recently published UV spectroscopic method that identifies protein unfolding. In this study we determine concentration limits to avoid protein unfolding of two model surfactants, namely polysorbate 20 and polysorbate 80, by correlating surfactant concentration with percentage 'unfolded' for three model proteins. For each scenario two distinct regions were observed, firstly surfactant concentrations at which no unfolding had occurred, followed by a second region whereby unfolding steadily increased with surfactant concentration. In general for the combinations analysed in this study, this second region began to appear around ten times below the critical micellar concentration of each surfactant, regardless of the protein or polysorbate chosen. It is therefore proposed that this adapted method could be used by researchers in the early stages of formulation development as a convenient and simple screening tool to confirm the 'onset of unfolding' concentration for protein-surfactant formulations, thus helping to optimise surfactant concentration selection in pharmaceutical formulations to maintain the benefits of surfactants yet avoid inadvertent unfolding.

Pathfinder: Protein folding pathway prediction based on conformational sampling

The study of protein folding mechanism is a challenge in molecular biology, which is of great significance for revealing the movement rules of biological macromolecules, understanding the pathogenic mechanism of folding diseases, and designing protein engineering materials. Based on the hypothesis that the conformational sampling trajectory contain the information of folding pathway, we propose a protein folding pathway prediction algorithm named Pathfinder. Firstly, Pathfinder performs large-scale sampling of the conformational space and clusters the decoys obtained in the sampling. The heterogeneous conformations obtained by clustering are named seed states. Then, a resampling algorithm that is not constrained by the local energy basin is designed to obtain the transition probabilities of seed states. Finally, protein folding pathways are inferred from the maximum transition probabilities of seed states. The proposed Pathfinder is tested on our developed test set (34 proteins). For 11 widely studied proteins, we correctly predicted their folding pathways and specifically analyzed 5 of them. For 13 proteins, we predicted their folding pathways to be further verified by biological experiments. For 6 proteins, we analyzed the reasons for the low prediction accuracy. For the other 4 proteins without biological experiment results, potential folding pathways were predicted to provide new insights into protein folding mechanism. The results reveal that structural analogs may have different folding pathways to express different biological functions, homologous proteins may contain common folding pathways, and α-helices may be more prone to early protein folding than β-strands.

The Jones polynomial of collections of open curves in 3-space

Measuring the entanglement complexity of collections of open curves in 3-space has been an intractable, yet pressing mathematical problem, relevant to a plethora of physical systems, such as in polymers and biopolymers. In this manuscript, we give a novel definition of the Jones polynomial that generalizes the classic Jones polynomial to collections of open curves in 3-space. More precisely, first we provide a novel definition of the Jones polynomial of linkoids (open link diagrams) and show that this is a well-defined single variable polynomial that is a topological invariant, which, for link-type linkoids, coincides with that of the corresponding link. Using the framework introduced in (Panagiotou E, Kauffman L. 2020 Proc. R. Soc. A 476 , 20200124. (( doi:10.1098/rspa.2020.0124 )), this enables us to define the Jones polynomial of collections of open and closed curves in 3-space. For collections of open curves in 3-space, the Jones polynomial has real coefficients and it is a continuous function of the curves’ coordinates. As the endpoints of the curves tend to coincide, the Jones polynomial tends to that of the resultant link. We demonstrate with numerical examples that the novel Jones polynomial enables us to characterize the topological/geometrical complexity of collections of open curves in 3-space for the first time.

The Local Topological Free Energy of the SARS-CoV-2 Spike Protein

The novel coronavirus SARS-CoV-2 infects human cells using a mechanism that involves binding and structural rearrangement of its Spike protein. Understanding protein rearrangement and identifying specific amino acids where mutations affect protein rearrangement has attracted much attention for drug development. In this manuscript, we use a mathematical method to characterize the local topology/geometry of the SARS-CoV-2 Spike protein backbone. Our results show that local conformational changes in the FP, HR1, and CH domains are associated with global conformational changes in the RBD domain. The SARS-CoV-2 variants analyzed in this manuscript (alpha, beta, gamma, delta Mink, G614, N501) show differences in the local conformations of the FP, HR1, and CH domains as well. Finally, most mutations of concern are either in or in the vicinity of high local topological free energy conformations, suggesting that high local topological free energy conformations could be targets for mutations with significant impact of protein function. Namely, the residues 484, 570, 614, 796, and 969, which are present in variants of concern and are targeted as important in protein function, are predicted as such from our model.