Converging experimental and computational views of the knotting mechanism of a small knotted protein

Discrepant Effects of Hydrated and Neat Reline on the Conformational Stability of a Knotted Protein

Although knotted proteins are rare in number, their peculiar topology has long intrigued the scientific community. In this study, we have explored the conformational stability of a trefoil-knotted protein, YbeA, in reline (choline chloride:urea in a 1:2 ratio), a well-characterized deep eutectic solvent, using classical molecular dynamics simulation. Deep eutectic solvents (DESs) are explored as a reliable alternative to conventional solvents, effectively altering a protein's structural stability and activity, either stabilizing its native state or disordering its conformation depending on the relevant interactions involved. Here, using pure and hydrated concentrations of reline, we observe the conflicting effect of the DES on the knotted protein's stability. Our studies at room temperature and elevated temperatures show that in pure reline, the protein is conformationally stable and rigid. In contrast, the protein tends to lose its structural integrity in hydrated reline. The stable knotted topology also gets untied as the protein, solvated in hydrated reline, is exposed to an elevated temperature. Using Minimum Distance Distribution Functions and Kirkwood-Buff Integrals, we analyzed the solvation pattern of the DES constituents around the protein. We expect that this study will lead to more effective strategies in developing tailored solvent systems for comprehending the conformational behavior of knotted proteins.

Effect of a Knot on the Thermal Stability of Protein MJ0366: Insights into Molecular Dynamics and Monte Carlo Simulations

Protein MJ0366 is a hypothetical protein from Methanocaldococcus jannaschii that has a rare and complex knot in its structure. The knot is a right-handed trefoil knot that involves about half of the protein's residues. In this work, we investigate the thermal stability of protein MJ0366 using numerical simulations based on molecular dynamics and Monte Carlo methods. We compare the results with those of a similar unknotted protein and analyze the effects of the knot on the folding and unfolding processes. We show that the knot in protein MJ0366 increases its thermal stability by creating a topological barrier that prevents the protein from unfolding at high temperatures.

The Critical Role of the Variable Domain in Driving Proteotoxicity and Aggregation in Full-length Light Chains

Light chain (AL) amyloidosis is the most common systemic amyloid disease characterized by abnormal accumulation of amyloid fibrils derived from immunoglobulin light chains (LCs). Both full-length (FL) LCs and their isolated variable (VL) and constant (CL) domains contribute to amyloid deposits in multiple organs, with the VL domain predominantly forming the fibril core. However, the role and interplay of these domains in amyloid aggregation and toxicity are poorly understood. Characterizing the amyloidogenic λ6-LC AL55, this study explores the properties of both FL and isolated domains and compares them with the available patient-derived data. FL AL55 biophysical features result from the interplay between its VL and CL domains where the limited VL-CL interface might play a major role. Slow refolding kinetic of FL confirms the unfolded VL domain as a kinetic trap possibly shifting the process towards misfolding. The X-ray structure of FL AL55 shows that VL domains may detach from the native dimeric assembly and establish non-native interdimeric interfaces. Additionally, isolated VL domains display significantly lower soluble toxicity compared to FL and do not form fibrils similar to those found ex vivo. Thus the data obtained in this work allowed us to draw a molecular sketch of the aggregation pathway for amyloidogenic LCs.

The mechanism of amyloid fibril growth from Φ-value analysis

Amyloid fibrils are highly stable misfolded protein assemblies that play an important role in several neurodegenerative and systemic diseases. Although structural information of the amyloid state is now abundant, mechanistic details about the misfolding process remain elusive. Inspired by the Φ-value analysis of protein folding, we combined experiments and molecular simulations to resolve amino-acid contacts and determine the structure of the transition-state ensemble-the rate-limiting step-for fibril elongation of PI3K-SH3 amyloid fibrils. The ensemble was validated experimentally by Tanford β analysis and computationally by free energy calculations. Although protein folding proceeds on funnel-shaped landscapes, here we find that the energy landscape for the misfolding reaction consists of a large 'golf course' region, defined by a single energy barrier and transition state, accessing a sharply funnelled region. Thus, misfolding occurs by rare, successful monomer-fibril end collisions interspersed by numerous unsuccessful binding attempts. Taken together, these insights provide a quantitative and highly resolved description of a protein misfolding reaction.

Salt-bridge mediated conformational dynamics in the figure-of-eight knotted ketol acid reductoisomerase (KARI)

The utility of knotted proteins in biological activities has been ambiguous since their discovery. From their evolutionary significance to their functionality in stabilizing the native protein structure, a unilateral conclusion hasn't been achieved yet. While most studies have been performed to understand the stabilizing effect of the knotted fold on the protein chain, more ideas are yet to emerge regarding the interactions in stabilizing the knot. Using classical molecular dynamics (MD) simulations, we have explored the dynamics of the figure-of-eight knotted domain present in ketol acid reductoisomerase (KARI). Our main focus was on the presence of a salt bridge network evident within the knotted region and its role in shaping the conformational dynamics of the knotted chain. Through the potential of mean forces (PMFs) calculation, we have also marked the specific salt bridges that are pivotal in stabilizing the knotted structure. The correlated motions have been further monitored with the help of principal component analysis (PCA) and dynamic cross-correlation maps (DCCM). Furthermore, mutation of the specific salt bridges led to a change in their conformational stability, vindicating their importance.

Is There a Functional Role for the Knotted Topology in Protein UCH-L1?

Knotted proteins are present in nature, but there is still an open issue regarding the existence of a universal role for these remarkable structures. To address this question, we used classical molecular dynamics (MD) simulations combined with in vitro experiments to investigate the role of the Gordian knot in the catalytic activity of UCH-L1. To create an unknotted form of UCH-L1, we modified its amino acid sequence by truncating several residues from its N-terminus. Remarkably, we find that deleting the first two N-terminal residues leads to a partial loss of enzyme activity with conservation of secondary structural content and knotted topological state. This happens because the integrity of the N-terminus is critical to ensure the correct alignment of the catalytic triad. However, the removal of five residues from the N-terminus, which significantly disrupts the native structure and the topological state, leads to a complete loss of enzymatic activity. Overall, our findings indicate that UCH-L1's catalytic activity depends critically on the integrity of the N-terminus and the secondary structure content, with the latter being strongly coupled with the knotted topological state.

Exploring the electronic structure of knotted proteins: the case of two ornithine transcarbamylase family

CONTEXT

Geometrical knots are rare structural arrangements in proteins in which the polypeptide chain ties itself into a knot, which is very intriguing due to the uncertainty of their impact on the protein properties. Presently, classical molecular dynamics is the most employed technique in the few studies found on this topic, so any information on how the presence of knots affects the reactivity and electronic properties of proteins is even scarcer. Using the electronic structure methods and quantum chemical descriptors analysis, we found that the same amino-acid residues in the knot core have statistically larger values for the unknotted protein, for both hard-hard and soft-soft interaction descriptors. In addition, we present a computationally feasible protocol, where we show it is possible to separate the contribution of the geometrical knot to the reactivity and other electronic structure properties.

METHODS

In order to investigate these systems, we used PRIMoRDiA, a new software developed by our research group, to explore the electronic structure of biological macromolecules. We evaluated several local quantum chemical descriptors to unveil relevant patterns potentially originating from the presence of the geometrical knot in two proteins, belonging to the ornithine transcarbamylase family. We compared several sampled structures from these two enzymes that are highly similar in both tertiary structure and function, but one of them has a knot whereas the other does not. The sampling was carried out through molecular dynamics simulations using ff14SB force field along 50 ns, and the semiempirical convergence was performed with PM7 Hamiltonian.

Tying a true topological protein knot by cyclization

Knotted proteins are fascinating to biophysicists because of their robust ability to fold into intricately defined three-dimensional structures with complex and topologically knotted arrangements. Exploring the biophysical properties of the knotted proteins is of significant interest, as they could offer enhanced chemical, thermal, and mechanostabilities. A true mathematical knot requires a closed path; in contrast, knotted protein structures have open N- and C-termini. To address the question of how a truly knotted protein differs from the naturally occurring counterpart, we enzymatically cyclized a 31 knotted YibK protein from Haemophilus influenza (HiYibK) to investigate the impact of path closure on its structure-function relationship and folding stability. Through the use of a multitude of structural and biophysical tools, including X-ray crystallography, NMR spectroscopy, small angle X-ray scattering, differential scanning calorimetry, and isothermal calorimetry, we showed that the path closure minimally perturbs the native structure and ligand binding of HiYibK. Nevertheless, the cyclization did alter the folding stability and mechanism according to chemical and thermal unfolding analysis. These molecular insights contribute to our fundamental understanding of protein folding and knotting that could have implications in the protein design with higher stabilities.

Topological Reaction Coordinate Captures the Folding Transition State Ensemble in a Pierced Lasso Protein

Proteins with a pierced lasso topology (PLT) have a covalent loop created by a disulfide bond, and the backbone circles back to thread the loop. This threaded topology has unique features compared to knotted topologies; notably, the topology is controlled by the chemical environment and the covalent loop remains intact even when denatured. In this work, we use the hormone leptin as our model PLT system and study its folding using molecular dynamics simulations that employ a structure-based (Go̅-like) model. We find that the reduced protein has a two-state folding mechanism with a transition state ensemble (TSE) that can be characterized by the reaction coordinate Q, the fraction of native contacts formed. In contrast, the oxidized protein, which must thread part of the polypeptide chain through a covalent loop, has a folding process that is poorly characterized by Q. Instead, we find that a topological coordinate that monitors the residue crossing the loop can identify the TSE of oxidized leptin. By precisely identifying the predicted TSE, one may now reliably calculate theoretical phi-values for the PLT protein, thereby enabling a comparison with experimental measurements. We find the loop-threading constraint leads to noncanonical phi-values that are uniformly small because this PLT protein has a flat energy landscape through the TSE.

Structure, dynamics, and stability of the smallest and most complex 71 protein knot

Proteins can spontaneously tie a variety of intricate topological knots through twisting and threading of the polypeptide chains. Recently developed artificial intelligence algorithms have predicted several new classes of topological knotted proteins, but the predictions remain to be authenticated experimentally. Here, we showed by X-ray crystallography and solution-state NMR spectroscopy that Q9PR55, an 89-residue protein from Ureaplasma urealyticum, possesses a novel 71 knotted topology that is accurately predicted by AlphaFold 2, except for the flexible N terminus. Q9PR55 is monomeric in solution, making it the smallest and most complex knotted protein known to date. In addition to its exceptional chemical stability against urea-induced unfolding, Q9PR55 is remarkably robust to resist the mechanical unfolding-coupled proteolysis by a bacterial proteasome, ClpXP. Our results suggest that the mechanical resistance against pulling-induced unfolding is determined by the complexity of the knotted topology rather than the size of the molecule.

Folding and functions of knotted proteins

Topologically knotted proteins have entangled structural elements within their native structures that cannot be disentangled simply by pulling from the N- and C-termini. Systematic surveys have identified different types of knotted protein structures, constituting as much as 1% of the total entries within the Protein Data Bank. Many knotted proteins rely on their knotted structural elements to carry out evolutionarily conserved biological functions. Being knotted may also provide mechanical stability to withstand unfolding-coupled proteolysis. Reconfiguring a knotted protein topology by circular permutation or cyclization provides insights into the importance of being knotted in the context of folding and functions. With the explosion of predicted protein structures by artificial intelligence, we are now entering a new era of exploring the entangled protein universe.

De novo design of knotted tandem repeat proteins

De novo protein design methods can create proteins with folds not yet seen in nature. These methods largely focus on optimizing the compatibility between the designed sequence and the intended conformation, without explicit consideration of protein folding pathways. Deeply knotted proteins, whose topologies may introduce substantial barriers to folding, thus represent an interesting test case for protein design. Here we report our attempts to design proteins with trefoil (31) and pentafoil (51) knotted topologies. We extended previously described algorithms for tandem repeat protein design in order to construct deeply knotted backbones and matching designed repeat sequences (N = 3 repeats for the trefoil and N = 5 for the pentafoil). We confirmed the intended conformation for the trefoil design by X ray crystallography, and we report here on this protein's structure, stability, and folding behaviour. The pentafoil design misfolded into an asymmetric structure (despite a 5-fold symmetric sequence); two of the four repeat-repeat units matched the designed backbone while the other two diverged to form local contacts, leading to a trefoil rather than pentafoil knotted topology. Our results also provide insights into the folding of knotted proteins.

Effects of sequence-dependent non-native interactions in equilibrium and kinetic folding properties of knotted proteins

Determining the role of non-native interactions in folding dynamics, kinetics, and mechanisms is a classic problem in protein folding. More recently, this question has witnessed a renewed interest in light of the hypothesis that knotted proteins require the assistance of non-native interactions to fold efficiently. Here, we conduct extensive equilibrium and kinetic Monte Carlo simulations of a simple off-lattice C-alpha model to explore the role of non-native interactions in the thermodynamics and kinetics of three proteins embedding a trefoil knot in their native structure. We find that equilibrium knotted conformations are stabilized by non-native interactions that are non-local, and proximal to native ones, thus enhancing them. Additionally, non-native interactions increase the knotting frequency at high temperatures, and in partially folded conformations below the transition temperatures. Although non-native interactions clearly enhance the efficiency of transition from an unfolded conformation to a partially folded knotted one, they are not required to efficiently fold a knotted protein. Indeed, a native-centric interaction potential drives the most efficient folding transition, provided that the simulation temperature is well below the transition temperature of the considered model system.

Long-lived metastable knots in polyampholyte chains

Knots in proteins and DNA are known to have significant effect on their equilibrium and dynamic properties as well as on their function. While knot dynamics and thermodynamics in electrically neutral and uniformly charged polymer chains are relatively well understood, proteins are generally polyampholytes, with varied charge distributions along their backbones. Here we use simulations of knotted polymer chains to show that variation in the charge distribution on a polyampholyte chain with zero net charge leads to significant variation in the resulting knot dynamics, with some charge distributions resulting in long-lived metastable knots that escape the (open-ended) chain on a timescale that is much longer than that for knots in electrically neutral chains. The knot dynamics in such systems can be described, quantitatively, using a simple one-dimensional model where the knot undergoes biased Brownian motion along a "reaction coordinate", equal to the knot size, in the presence of a potential of mean force. In this picture, long-lived knots result from charge sequences that create large electrostatic barriers to knot escape. This model allows us to predict knot lifetimes even when those times are not directly accessible by simulations.

Computational Design of Extraordinarily Stable Peptide Structures through Side-Chain-Locked Knots

Extraordinarily stable protein and peptide structures are critically demanded in many applications. Typical approaches to enhance protein and peptide stability are strengthening certain interactions. Here, we develop a very different approach: stabilizing peptide structures through side-chain-locked knots. More specifically, a peptide core consists of a knot, which is prevented from unknotting and unfolding by large side chains of amino acids at knot boundaries. These side chains impose free energy barriers for unknotting. The free energy barriers are quantified using all-atom and coarse-grained simulations. The barriers become infinitely high for large side chains and tight knot cores, resulting in stable peptide structures, which never unfold unless one chemical bond is broken. The extraordinary stability is essentially kinetic stability. Our new approach lifts the thermodynamic restriction in designing peptide structures, provides extra freedom in selecting sequence and structural motifs that are thermodynamically unstable, and should expand the functionality of peptides. This work also provides a bottom-up understanding of how knotting enhances protein stability.

In silico method for identifying the key residues in a knotted protein: with MJ0366 as an example

A simple in silico method for predicting the key residues for knotting and unknotting a knotted protein is put forward, with the residues ranked by the relevance to knotting and unknotting in the annealing molecular dynamics simulations.

Tying up the Loose Ends: A Mathematically Knotted Protein

Knots have attracted scientists in mathematics, physics, biology, and engineering. Long flexible thin strings easily knot and tangle as experienced in our daily life. Similarly, long polymer chains inevitably tend to get trapped into knots. Little is known about their formation or function in proteins despite >1,000 knotted proteins identified in nature. However, these protein knots are not mathematical knots with their backbone polypeptide chains because of their open termini, and the presence of a “knot” depends on the algorithm used to create path closure. Furthermore, it is generally not possible to control the topology of the unfolded states of proteins, therefore making it challenging to characterize functional and physicochemical properties of knotting in any polymer. Covalently linking the amino and carboxyl termini of the deeply trefoil-knotted YibK from Pseudomonas aeruginosa allowed us to create the truly backbone knotted protein by enzymatic peptide ligation. Moreover, we produced and investigated backbone cyclized YibK without any knotted structure. Thus, we could directly probe the effect of the backbone knot and the decrease in conformational entropy on protein folding. The backbone cyclization did not perturb the native structure and its cofactor binding affinity, but it substantially increased the thermal stability and reduced the aggregation propensity. The enhanced stability of a backbone knotted YibK could be mainly originated from an increased ruggedness of its free energy landscape and the destabilization of the denatured state by backbone cyclization with little contribution from a knot structure. Despite the heterogeneity in the side-chain compositions, the chemically unfolded cyclized YibK exhibited several macroscopic physico-chemical attributes that agree with theoretical predictions derived from polymer physics.