Folding Pathways of a Knotted Protein with a Realistic Atomistic Force Field

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native -sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting.     

Knotted fusion proteins reveal unexpected possibilities in protein folding

Proteins that contain a distinct knot in their native structure are impressive examples of biological self-organization. Although this topological complexity does not appear to cause a folding problem, the mechanisms by which such knotted proteins form are unknown. We found that the fusion of an additional protein domain to either the amino terminus, the carboxy terminus, or to both termini of two small knotted proteins did not affect their ability to knot. The multidomain constructs remained able to fold to structures previously thought unfeasible, some representing the deepest protein knots known. By examining the folding kinetics of these fusion proteins, we found evidence to suggest that knotting is not rate limiting during folding, but instead occurs in a denatured-like state. These studies offer experimental insights into when knot formation occurs in natural proteins and demonstrate that early folding events can lead to diverse and sometimes unexpected protein topologies.     

Aromatic residues engineered into the beta-turn nucleation site of ubiquitin lead to a complex folding landscape, non-native side-chain interactions, and kinetic traps

The fast folding of small proteins is likely to be the product of evolutionary pressures that balance the search for native-like contacts in the transition state with the minimum number of stable non-native interactions that could lead to partially folded states prone to aggregation and amyloid formation. We have investigated the effects of non-native interactions on the folding landscape of yeast ubiquitin by introducing aromatic substitutions into the beta-turn region of the N-terminal beta-hairpin, using both the native G-bulged type I turn sequence (TXTGK) as well as an engineered 2:2 XNGK type I' turn sequence. The N-terminal beta-hairpin is a recognized folding nucleation site in ubiquitin. The folding kinetics for wt-Ub (TLTGK) and the type I' turn mutant (TNGK) reveal only a weakly populated intermediate, however, substitution with X = Phe or Trp in either context results in a high propensity to form a stable compact intermediate where the initial U-->I collapse is visible as a distinct kinetic phase. The introduction of Trp into either of the two host turn sequences results in either complex multiphase kinetics with the possibility of parallel folding pathways, or formation of a highly compact I-state stabilized by non-native interactions that must unfold before refolding. Sequence substitutions with aromatic residues within a localized beta-turn capable of forming non-native hydrophobic contacts in both the native state and partially folded states has the undesirable consequence that folding is frustrated by the formation of stable compact intermediates that evolutionary pressures at the sequence level may have largely eliminated.     

A circumventing role for the non-native intermediate in the folding of β-lactoglobulin

Folding experiments have suggested that some proteins have kinetic intermediates with a non-native structure. A simple G ?o model does not explain such non-native intermediates. Therefore, the folding energy landscape of proteins with non-native intermediates should have characteristic properties. To identify such properties, we investigated the folding of bovine β-lactoglobulin (βLG). This protein has an intermediate with a non-native α-helical structure, although its native form is predominantly composed of β-structure. In this study, we prepared mutants whose α-helical and β-sheet propensities are modified and observed their folding using a stopped-flow circular dichroism apparatus. One interesting finding was that E44L, whose β-sheet propensity was increased, showed a folding intermediate with an amount of β-structure similar to that of the wild type, though its folding took longer. Thus, the intermediate seems to be a trapped intermediate. The high α-helical propensity of the wild-type sequence likely causes the folding pathway to circumvent such time-consuming intermediates. We propose that the role of the non-native intermediate is to control the pathway at the beginning of the folding reaction.     

Yeast epiarginase regulation,an enzyme-enzyme activity control: identification of residues of ornithine carbamoyltransferase and arginase responsible for enzyme catalytic and regulatory activities

In the presence of ornithine and arginine, ornithine carbamoyltransferase (OTCase) and arginase form a one-to-one enzyme complex in which the activity of OTCase is inhibited whereas arginase remains catalytically active. The mechanism by which these nonallosteric enzymes form a stable complex triggered by the binding of their respective substrates raises the question of how such a cooperative association is induced. Analyses of mutations in both enzymes identify residues that are required for their association, some of them being important for catalysis. In arginase, two cysteines at the C terminus of the protein are crucial for its epiarginase function but not for its catalytic activity and trimeric structure. In OTCase, mutations of putative ornithine binding residues, Asp-182, Asn-184, Asn-185, Cys-289, and Glu-256 greatly reduced the affinity for ornithine and impaired the interaction with arginase. The four lysine residues located in the SMG loop, Lys-260, Lys-263, Lys-265, and Lys-268, also play an important role in mediating the sensitivity of OTCase to ornithine and to arginase and appear to be involved in transducing and enhancing the signal given by ornithine for the closure of the catalytic domain.     

Metabolic channeling of carbamoyl phosphate, a thermolabile intermediate: evidence for physical interaction between carbamate kinase-like carbamoyl-phosphate synthetase and ornithine carbamoyltransferase from the hyperthermophile Pyrococcus furiosus

Two different approaches provided evidence for a physical interaction between the carbamate kinase-like carbamoyl-phosphate synthetase (CKase) and ornithine carbamoyltransferase (OTCase) from the hyperthermophilic archaeon Pyrococcus furiosus. Affinity electrophoresis indicated that CKase and OTCase associate into a multienzyme cluster. Further evidence for a biologically significant interaction between CKase and OTCase was obtained by co-immunoprecipitation combined with formaldehyde cross-linking experiments. These experiments support the hypothesis that CKase and OTCase form an efficient channeling cluster for carbamoyl phosphate, an extremely thermolabile and potentially toxic metabolic intermediate. Therefore, by physically interacting with each other, CKase and OTCase prevent the thermodenaturation of carbamoyl phosphate in the aqueous cytoplasmic environment.     

Duplex DNA knots produced by Escherichia coli topoisomerase I. Structure and requirements for formation

We investigated systematically the knotting of nicked circular duplex DNA by Escherichia coli topoisomerase I. Agarose gel electrophoresis of knots forms a ladder of DNA bands. Each rung is made up of a variety of knots with the same number of nodes, or segment crossings; knots in adjacent rungs differ by one node. We extended the technique of electron microscopy of recA protein-coated DNA to the visualization of the complex knots tied by topoisomerase I. The striking result is that the enzyme produces every knot theoretically possible. The requirement for excess enzyme to form complex knots suggests a role for topoisomerase I in contorting the DNA in addition to promoting strand passage. We conclude that nodes formed are equally likely to be positive or negative and that topoisomerase I can pass DNA strands through a transient enzyme-generated break without regard to orientation of the passing strand. The results are interpreted in terms of a formulation for the topological requirements for knotting.     

Evolutionary divergence of genes for ornithine and aspartate carbamoyl-transferases--complete sequence and mode of regulation of the Escherichia coli argF gene; comparison of argF with argI and pyrB.

The complete nucleotide sequence of argF is presented, together with that of an operator-constitutive mutant. ArgF is compared with the other gene coding for ornithine carbamoyltransferase (OTCase) in E. coli K-12, argI, and with pyrB, encoding the catalytic monomer of aspartate carbamoyltransferase (ATCase). ArgF and argI appear very closely related having emerged from a relatively recent ancestor gene. The relationship between OTCase and ATCase appears more distant. Nevertheless, the homology observed between the two proteins (mainly in the polar domain) suggests a common origin.     

Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state

Comparatively little is known about the role of non-native interactions in protein folding and their role in both folding and stability is controversial. We demonstrate that non-native electrostatic interactions involving specific residues in the denatured state can have a significant effect upon protein stability and can persist in the transition state for folding. Mutation of a single surface exposed residue, Lys12 to Met, in the N-terminal domain of the ribosomal protein L9 (NTL9), significantly increased the stability of the protein and led to faster folding. Structural and energetic studies of the wild-type and K12M mutant show that the 1.9 kcal mol(-1) increase in stability is not due to native state effects, but rather is caused by modulation of specific non-native electrostatic interactions in the denatured state. pH dependent stability measurements confirm that the increased stability of the K12M is due to the elimination of favorable non-native interactions in the denatured state. Kinetic studies show that the non-native electrostatic interactions involving K12 persist in the transition state. The analysis demonstrates that canonical Phi-values can arise from the disruption of non-native interactions as well as from the development of native interactions.     

Statistics of knots, geometry of conformations, and evolution of proteins

Like shoelaces, the backbones of proteins may get entangled and form knots. However, only a few knots in native proteins have been identified so far. To more quantitatively assess the rarity of knots in proteins, we make an explicit comparison between the knotting probabilities in native proteins and in random compact loops. We identify knots in proteins statistically, applying the mathematics of knot invariants to the loops obtained by complementing the protein backbone with an ensemble of random closures, and assigning a certain knot type to a given protein if and only if this knot dominates the closure statistics (which tells us that the knot is determined by the protein and not by a particular method of closure). We also examine the local fractal or geometrical properties of proteins via computational measurements of the end-to-end distance and the degree of interpenetration of its subchains. Although we did identify some rather complex knots, we show that native conformations of proteins have statistically fewer knots than random compact loops, and that the local geometrical properties, such as the crumpled character of the conformations at a certain range of scales, are consistent with the rarity of knots. From these, we may conclude that the known “protein universe” (set of native conformations) avoids knots. However, the precise reason for this is unknown—for instance, if knots were removed by evolution due to their unfavorable effect on protein folding or function or due to some other unidentified property of protein evolution.     

A single mutation in the active site swaps the substrate specificity of N-acetyl-L-ornithine transcarbamylase and N-succinyl-L-ornithine transcarbamylase

Transcarbamylases catalyze the transfer of the carbamyl group from carbamyl phosphate (CP) to an amino group of a second substrate such as aspartate, ornithine, or putrescine. Previously, structural determination of a transcarbamylase from Xanthomonas campestris led to the discovery of a novel N-acetylornithine transcarbamylase (AOTCase) that catalyzes the carbamylation of N-acetylornithine. Recently, a novel N-succinylornithine transcarbamylase (SOTCase) from Bacteroides fragilis was identified. Structural comparisons of AOTCase from X. campestris and SOTCase from B. fragilis revealed that residue Glu92 (X. campestris numbering) plays a critical role in distinguishing AOTCase from SOTCase. Enzymatic assays of E92P, E92S, E92V, and E92A mutants of AOTCase demonstrate that each of these mutations converts the AOTCase to an SOTCase. Similarly, the P90E mutation in B. fragilis SOTCase (equivalent to E92 in X. campestris AOTCase) converts the SOTCase to AOTCase. Hence, a single amino acid substitution is sufficient to swap the substrate specificities of AOTCase and SOTCase. X-ray crystal structures of these mutants in complexes with CP and N-acetyl-L-norvaline (an analog of N-acetyl-L-ornithine) or N-succinyl-L-norvaline (an analog of N-succinyl-L-ornithine) substantiate this conversion. In addition to Glu92 (X. campestris numbering), other residues such as Asn185 and Lys30 in AOTCase, which are involved in binding substrates through bridging water molecules, help to define the substrate specificity of AOTCase. These results provide the correct annotation (AOTCase or SOTCase) for a set of the transcarbamylase-like proteins that have been erroneously annotated as ornithine transcarbamylase (OTCase, EC 2.1.3.3).     

Atomic-level structure characterization of an ultrafast folding mini-protein denatured state

Atomic-level analyses of non-native protein ensembles constitute an important aspect of protein folding studies to reach a more complete understanding of how proteins attain their native form exhibiting biological activity. Previously, formation of hydrophobic clusters in the 6 M urea-denatured state of an ultrafast folding mini-protein known as TC5b from both photo-CIDNP NOE transfer studies and FCS measurements was observed. Here, we elucidate the structural properties of this mini-protein denatured in 6 M urea performing (15)N NMR relaxation studies together with a thorough NOE analysis. Even though our results demonstrate that no elements of secondary structure persist in the denatured state, the heterogeneous distribution of R(2) rate constants together with observing pronounced heteronuclear NOEs along the peptide backbone reveals specific regions of urea-denatured TC5b exhibiting a high degree of structural rigidity more frequently observed for native proteins. The data are complemented with studies on two TC5b point mutants to verify the importance of hydrophobic interactions for fast folding. Our results corroborate earlier findings of a hydrophobic cluster present in urea-denatured TC5b comprising both native and non-native contacts underscoring their importance for ultra rapid folding. The data assist in finding ways of interpreting the effects of pre-existing native and/or non-native interactions on the ultrafast folding of proteins; a fact, which might have to be considered when defining the starting conditions for molecular dynamics simulation studies of protein folding.     

A residue in helical conformation in the native state adopts a β-strand conformation in the folding transition state despite its high and canonical Φ-value

Delineating structures of the transition states in protein folding reactions has provided great insight into the mechanisms by which proteins fold. The most common method for obtaining this information is Φ-value analysis, which is carried out by measuring the changes in the folding and unfolding rates caused by single amino acid substitutions at various positions within a given protein. Canonical Φ-values range between 0 and 1, and residues displaying high values within this range are interpreted to be important in stabilizing the transition state structure, and to elicit this stabilization through native-like interactions. Although very successful in defining the general features of transition state structures, Φ-value analysis can be confounded when non-native interactions stabilize this state. In addition, direct information on backbone conformation within the transition state is not provided. In the work described here, we have investigated structure formation at a conserved β-bulge (with helical conformation) in the Fyn SH3 domain by characterizing the effects of substituting all natural amino acids at one position within this structural motif. By comparing the effects on folding rates of these substitutions with database-derived local structure propensity values, we have determined that this position adopts a non-native backbone conformation in the folding transition state. This result is surprising because this position displays a high and canonical Φ-value of 0.7. This work emphasizes the potential role of non-native conformations in folding pathways and demonstrates that even positions displaying high and canonical Φ-values may, nevertheless, adopt a non-native conformation in the transition state.     

Effective stiffening of DNA due to nematic ordering causes DNA molecules packed in phage capsids to preferentially form torus knots

Observation that DNA molecules in bacteriophage capsids preferentially form torus type of knots provided a sensitive gauge to evaluate various models of DNA arrangement in phage heads. Only models resulting in a preponderance of torus knots could be considered as close to reality. Recent studies revealed that experimentally observed enrichment of torus knots can be qualitatively reproduced in numerical simulations that include a potential inducing nematic arrangement of tightly packed DNA molecules within phage capsids. Here, we investigate what aspects of the nematic arrangement are crucial for inducing formation of torus knots. Our results indicate that the effective stiffening of DNA by the nematic arrangement not only promotes knotting in general but is also the decisive factor in promoting formation of DNA torus knots in phage capsids.     

Structure of a loose dimer: an intermediate in nitric oxide synthase assembly

Cooperativity among ligand binding, subunit association, and protein folding has implications for enzyme regulation as well as protein aggregation events associated with disease. The binding of substrate l-arginine or cofactor tetrahydrobiopterin converts nitric oxide synthases (NOSs) from a "loose dimer", with an exposed active center and higher sensitivity to proteolysis, to a "tight dimer" competent for catalysis. The crystallographic structure of the Bacillus subtilis NOS loose dimer shows an altered association state with severely destabilized subdomains. Ligand binding or heme reduction converts loose dimers to tight dimers in solution and crystals. Mutations at key positions in the dimer interface that distinguish prokaryotic from eukaryotic NOSs affect the propensity to form loose dimers. The loose dimer structure indicates that non-native interactions can mediate subunit association in NOS.     

The effect of increasing the stability of non-native interactions on the folding landscape of the bacterial immunity protein Im9

How stabilising non-native interactions influence protein folding energy landscapes is currently not well understood: such interactions could speed folding by reducing the conformational search to the native state, or could slow folding by increasing ruggedness. Here, we examine the influence of non-native interactions in the folding process of the bacterial immunity protein Im9, by exploiting our ability to manipulate the stability of the intermediate and rate-limiting transition state (TS) in the folding of this protein by minor alteration of its sequence or changes in solvent conditions. By analysing the properties of these species using Phi-value analysis, and exploration of the structural properties of the TS ensemble using molecular dynamics simulations, we demonstrate the importance of non-native interactions in immunity protein folding and demonstrate that the rate-limiting step involves partial reorganisation of these interactions as the TS ensemble is traversed. Moreover, we show that increasing the contribution to stability made by non-native interactions results in an increase in Phi-values of the TS ensemble without altering its structural properties or solvent-accessible surface area. The data suggest that the immunity proteins fold on multiple, but closely related, micropathways, resulting in a heterogeneous TS ensemble that responds subtly to mutation or changes in the solvent conditions. Thus, altering the relative strength of native and non-native interactions influences the search to the native state by restricting the pathways through the folding energy landscape.     

The mechanism of beta-hairpin formation

Beta-hairpins constitute an important class of connecting protein secondary structures. Several groups have postulated that such structures form early in the folding process and serve to nucleate the formation of extended beta-sheet structures. Despite the importance of beta-hairpins in protein folding, little is known about the mechanism of formation of these structures. While it is well established that there is a complex interplay between the stability of a beta-hairpin and loop conformational propensity, loop length, and the formation of stabilizing cross-strand interactions (H-bonds and hydrophobic interactions), the influence of these factors on the folding rate is poorly understood. Peptide models provide a simple framework for exploring the molecular details of the formation of beta-hairpin structures. We have explored the fundamental processes of folding in two linear peptides that form beta-hairpin structures, having a stabilizing hydrophobic cluster connected by loops of differing lengths. This approach allows us to evaluate existing models of the mechanism of beta-hairpin formation. We find a substantial acceleration of the folding rate when the connecting loop is made shorter (i.e., the hydrophobic cluster is moved closer to the turn). Analysis of the folding kinetics of these two peptides reveals that this acceleration is a direct consequence of the reduced entropic cost of the smaller loop search.     

AlphaFold predicts the most complex protein knot and composite protein knots

The computer artificial intelligence system AlphaFold has recently predicted previously unknown three‐dimensional structures of thousands of proteins. Focusing on the subset with high‐confidence scores, we algorithmically analyze these predictions for cases where the protein backbone exhibits rare topological complexity, that is, knotting. Amongst others, we discovered a 7₁‐knot, the most topologically complex knot ever found in a protein, as well several six‐crossing composite knots comprised of two methyltransferase or carbonic anhydrase domains, each containing a simple trefoil knot. These deeply embedded composite knots occur evidently by gene duplication and interconnection of knotted dimers. Finally, we report two new five‐crossing knots including the first 5₁‐knot. Our list of analyzed structures forms the basis for future experimental studies to confirm these novel‐knotted topologies and to explore their complex folding mechanisms.     

Purification and characterization of Arabidopsis ornithine transcarbamoylase (OTCase), a member of a distinct and evolutionarily-conserved group of plant OTCases

We affinity-purified an ornithine transcarbamoylase (carbamoyl phosphate:L-ornithine carbamoyltransferase; EC 2.1.3.3) 676-fold to near homogeneity from leaf tissues of Arabidopsis thaliana L. cv. Columbia. The purified OTCase protein exhibited a molecular mass of 37 kDa on SDS-PAGE gels and exhibited a pI = 6.8. A 41-kDa polypeptide was immunoprecipitated from Arabidopsis leaf poly(A)⁺ RNA in vitro translation products by pea OTCase antiserum. This precursor OTCase (pOTCase) is the predicted size (41 170 Da) for a polypeptide encoded by an Arabidopsis OTCase cDNA. Characteristics of N-terminal residues of the deduced amino acid sequence of this pOTCase suggest that it is a chloroplast-targeted protein. The sequences of plant OTCases suggest that they represent a distinct and evolutionarily-conserved group of OTCases. No evidence was found for OTCase isoenzymes in Arabidopsis leaf tissues. The Arabidopsis pOTCase was poorly-expressed in Escherichia coli strain TB-2, an OTCase-deficient mutant, and did not complement the mutant on arginine-minus selection medium.     

Engineering a compact non-native state of intestinal fatty acid-binding protein

The last three C-terminal residues (129-131) of intestinal fatty acid-binding protein (IFABP) participate in four main-chain hydrogen bonds and two electrostatic interactions to sequentially distant backbone and side-chain atoms. To assess if these interactions are involved in the final adjustment of the tertiary structure during folding, we engineered an IFABP variant truncated at residue 128. An additional mutation, Trp-6-->Phe, was introduced to simplify the conformational analysis by optical methods. Although the changes were limited to a small region of the protein surface, they resulted in an IFABP with altered secondary and tertiary structure. Truncated IFABP retains some cooperativity, is monomeric, highly compact, and has the molecular dimensions and shape of the native protein. Our results indicated that residues 129-131 are part of a crucial conformational determinant in which several long-range interactions, essential for the acquisition of the native state, are established. This work suggests that carefully controlled truncation can populate equilibrium non-native states under physiological conditions. These non-native states hold a great promise as experimental models for protein folding.