Thermotoga Maritima Research Articles

Biosynthesis of Lacto-N-biose I from starch and N-acetylglucosamine via an in vitro synthetic enzymatic biosystem

Human milk oligosaccharides (HMOs) are very distinctive components in human milk and are beneficial for infant health. Lacto-N-biose I (LNB) is the core structural unit of HMOs, which could be used for the synthesis of other HMOs. In this study, an ATP-free in vitro synthetic enzymatic biosystem contained four thermostable enzymes (alpha-glucan phosphorylase from Thermococcus kodakarensis, UDP-glucose-hexose-1-phosphate uridylyltransferase from Thermotoga maritima, UDP-glucose 4-epimerase from T. maritima, lacto-N-biose phosphorylase from Clostridium thermobutyricum) were constructed for the biosynthesis of LNB from starch and N-acetylglucosamine (GlcNAc). Under the optimal conditions, 0.75 g/L and 2.23 g/L LNB was produced from 1.1 g/L and 4.4 g/L GlcNAc and excess starch, with the molar yield of 39% and 29% based on the GlcNAc concentration, respectively, confirming the feasibility of this in vitro synthetic enzymatic biosystem for LNB synthesis and shedding light on the biosynthesis of other HMOs using LNB as the core structural unit from low cost polysaccharides.

The Role of a Loop in the Non-catalytic Domain B on the Hydrolysis/Transglycosylation Specificity of the 4-α-Glucanotransferase from Thermotoga maritima.

The mechanism by which glycoside hydrolases control the reaction specificity through hydrolysis or transglycosylation is a key element embedded in their chemical structures. The determinants of reaction specificity seem to be complex. We looked for structural differences in domain B between the 4-α-glucanotransferase from Thermotoga maritima (TmGTase) and the α-amylase from Thermotoga petrophila (TpAmylase) and found a longer loop in the former that extends towards the active site carrying a W residue at its tip. Based on these differences we constructed the variants W131G and the partial deletion of the loop at residues 120-124/128-131, which showed a 11.6 and 11.4-fold increased hydrolysis/transglycosylation (H/T) ratio relative to WT protein, respectively. These variants had a reduction in the maximum velocity of the transglycosylation reaction, while their affinity for maltose as the acceptor was not substantially affected. Molecular dynamics simulations allow us to rationalize the increase in H/T ratio in terms of the flexibility near the active site and the conformations of the catalytic acid residues and their associated pKas.

Alpha-1,4-transglycosylation Activity of GH57 Glycogen Branching Enzymes Is Higher in the Absence of a Flexible Loop with a Conserved Tyrosine Residue.

Starch-like polymers can be created through the use of enzymatic modification with glycogen branching enzymes (GBEs). GBEs are categorized in the glycoside hydrolase (GH) family 13 and 57. Both GH13 and GH57 GBEs exhibit branching and hydrolytic activity. While GH13 GBEs are also capable of α-1,4-transglycosylation, it is yet unknown whether GH57 share this capability. Among the four crystal structures of GH57 GBEs that have been solved, a flexible loop with a conserved tyrosine was identified to play a role in the branching activity. However, it remains unclear whether this flexible loop is also involved in α-1,4-transglycosylation activity. We hypothesize that GH57 GBEs with the flexible loop and tyrosine are also capable of α-1,4-transglycosylation, similar to GH13 GBEs. The aim of the present study was to characterize the activity of GH57 GBEs to investigate a possible α-1,4-transglycosylation activity. Three GH57 GBEs were selected, one from Thermococcus kodakarensis with the flexible loop and two beta-strands; one from Thermotoga maritima, missing the flexible loop and beta-strands; and one from Meiothermus sp., missing the flexible loop but with the two beta-strands. The analysis of chain length distribution over time of modified maltooctadecaose, revealed, for the first time, that all three GH57 GBEs can generate chains longer than the substrate itself, showing that α-1,4-transglycosylation activity is generally present in GH57 GBEs.

Structure of reverse gyrase with a minimal latch that supports ATP-dependent positive supercoiling without specific interactions with the topoisomerase domain.

Reverse gyrase is the only topoisomerase that introduces positive supercoils into DNA in an ATP-dependent reaction. Positive DNA supercoiling becomes possible through the functional cooperation of the N-terminal helicase domain of reverse gyrase with its C-terminal type IA topoisomerase domain. This cooperation is mediated by a reverse-gyrase-specific insertion into the helicase domain termed the `latch'. The latch consists of a globular domain inserted at the top of a β-bulge loop that connects this globular part to the helicase domain. While the globular domain shows little conservation in sequence and length andis dispensable for DNA supercoiling, the β-bulge loop is required for supercoiling activity. It has previously been shown that the β-bulge loop constitutes a minimal latch that couples ATP-dependent processes in the helicase domain to DNA processing by the topoisomerase domain. Here, the crystal structure of Thermotoga maritima reverse gyrase with such a β-bulge loop as a minimal latch is reported. It is shown that the β-bulge loop supports ATP-dependent DNA supercoiling of reverse gyrase without engaging in specific interactions with the topoisomerase domain. When only a small latch or no latch is present, a helix in the nearby helicase domain of T. maritima reverse gyrase partially unfolds. Comparison of the sequences and predicted structures of latch regions in other reverse gyrases shows that neither sequence nor structure are decisive factors for latch functionality; instead, the decisive factors are likely to be electrostatics and plain steric bulk.

Single-stranded DNA recruitment mechanism in replication origin unwinding by DnaA initiator protein and HU, an evolutionary ubiquitous nucleoid protein.

The Escherichia coli replication origin oriC contains the initiator ATP-DnaA-Oligomerization Region (DOR) and its flanking duplex unwinding element (DUE). In the Left-DOR subregion, ATP-DnaA forms a pentamer by binding to R1, R5Mand three other DnaA boxes. The DNA-bending protein IHF binds sequence-specifically to the interspace between R1 and R5M boxes, promoting DUE unwinding, which is sustained predominantly by binding of R1/R5M-bound DnaAs to the single-stranded DUE (ssDUE). The present study describes DUE unwinding mechanisms promoted by DnaA and IHF-structural homolog HU, a ubiquitous protein in eubacterial species that binds DNA sequence-non-specifically, preferring bent DNA. Similar to IHF, HU promoted DUE unwinding dependent on ssDUE binding of R1/R5M-bound DnaAs. Unlike IHF, HU strictly required R1/R5M-bound DnaAs and interactions between the two DnaAs. Notably, HU site-specifically bound the R1-R5M interspace in a manner stimulated by ATP-DnaA and ssDUE. These findings suggest a model that interactions between the two DnaAs trigger DNA bending within the R1/R5M-interspace and initial DUE unwinding, which promotes site-specific HU binding that stabilizes the overall complex and DUE unwinding. Moreover, HU site-specifically bound the replication origin of the ancestral bacterium Thermotoga maritima depending on the cognate ATP-DnaA. The ssDUE recruitment mechanism could be evolutionarily conserved in eubacteria.

The cell envelope of Thermotogae suggests a mechanism for outer membrane biogenesis

The presence of a cell membrane is one of the major structural components defining life. Recent phylogenomic analyses have supported the hypothesis that the last universal common ancestor (LUCA) was likely a diderm. Yet, the mechanisms that guided outer membrane (OM) biogenesis remain unknown. Thermotogae is an early-branching phylum with a unique OM, the toga. Here, we use cryo-electron tomography to characterize the insitu cell envelope architecture of Thermotoga maritima and show that the toga is made of extended sheaths of β-barrel trimers supporting small (~200nm) membrane patches. Lipidomic analyses identified the same major lipid species in the inner membrane (IM) and toga, including the rare to bacteria membrane-spanning ether-bound diabolic acids (DAs). Proteomic analyses revealed that the toga was composed of multiple SLH-domain containing Ompα and novel β-barrel proteins, and homology searches detected variable conservations of these proteins across the phylum. These results highlight that, in contrast to the SlpA/OmpM superfamily of proteins, Thermotoga possess a highly diverse bipartite OM-tethering system. We discuss the implications of our findings with respect to other early-branching phyla and propose that a toga-like intermediate may have facilitated monoderm-to-diderm cell envelope transitions.

Screening the Thermotoga maritima genome for new wide-spectrum nucleoside and nucleotide kinases

Enzymes from thermophilic organisms are interesting biocatalysts for a wide variety of applications in organic synthesis, biotechnology, and molecular biology. Next to an increased stability at elevated temperatures, they were described to show a wider substrate spectrum than their mesophilic counterparts. To identify thermostable biocatalysts for the synthesis of nucleotide analogs, we performed a database search on the carbohydrate and nucleotide metabolism of Thermotoga maritima. After expression and purification of 13 enzyme candidates involved in nucleotide synthesis, these enzymes were screened for their substrate scope. We found that the synthesis of 2'-deoxynucleoside 5'-monophosphates (dNMPs) and uridine 5'-monophosphate from nucleosides was catalyzed by the already known wide-spectrum thymidine kinase and the ribokinase. In contrast, no NMP-forming activity was detected for adenosine-specific kinase, uridine kinase, or nucleotidase. The NMP kinases (NMPKs) and the pyruvate-phosphate-dikinase of T.maritima exhibited a rather specific substrate spectrum for the phosphorylation of NMPs, while pyruvate kinase, acetate kinase, and three of the NMPKs showed a broad substrate scope with (2'-deoxy)nucleoside 5'-diphosphates as substrates. Based on these promising results, TmNMPKs were applied in enzymatic cascade reactions for nucleoside 5'-triphosphate synthesis using four modified pyrimidine nucleosides and four purine NMPs as substrates, and we determined that base- and sugar-modified substrates were accepted. In summary, besides the already reported TmTK, NMPKs of T.maritima were identified to be interesting enzyme candidates for the enzymatic production of modified nucleotides.

Orientation of Cel5A and Xyn10B in a fusion construct is important in facilitating synergistic degradation of plant biomass polysaccharides

One approach to achieve efficient and economical saccharification of plant biomass would be using thermostable and multifunctional enzymes from hyperthermophiles such as Thermotoga maritima. Thus, the bifunctional constructs, Cel5A-Xyn10B and Xyn10B-Cel5A, were produced by fusing cellulase Cel5A at the N- and C-terminals of xylanase Xyn10B, respectively. The Cel5A-Xyn10B fusion construct showed cellulase activity of 1483 Uμmol-1 against carboxymethyl cellulose, which was nearly the same as that of Cel5A in the free form. However, xylanase activity of this construct increased by 2-fold against beechwood xylan as compared to that of Xyn10B in free form. The synergistic effect between Cel5A and Xyn10B in the form of Cel5A-Xyn10B fusion resulted an overall increase in the release of reducing sugars. However, Xyn10B-Cel5A showed about 60% decrease in activities of both the component enzymes as compared to their activities in the free form. Both the fusion constructs were active in a wide range of pH from 4.0 to 9.0 and temperatures from 50 to 90°C. Nearly 80% of cellulase and xylanase activities were retained in Cel5A-Xyn10B fusion after incubation at 60°C for 1h. Secondary structures of the component enzymes were retained in the Cel5A-Xyn10B fusion as observed by circular dichroism spectroscopy. Docking and simulation studies suggested that the enhanced xylanase activity in Cel5A-Xyn10B was due to the high binding energy, favorable orientation of the active sites, as well as relative positioning of the active site residues of Cel5A and Xyn10B in closer vicinity, which facilitated the substrate channeling.

Engineering TM1459 for Stabilisation against Inactivation by Amino Acid Oxidation

AbstractOxidative alkene cleavage is a highly interesting reaction to obtain aldehydes and ketones. The Mn‐dependent protein TM1459 from Thermotoga maritima can catalyse alkene cleavage of styrene derivatives in the presence of tert‐butyl hydroperoxide. Despite the high thermal stability of the enzyme, it gets inactivated during the reaction. The data reported here indicate that auto‐oxidation is responsible for the low stability of TM1459 in the oxidative environment required for the alkene cleavage reaction. By targeting the exchange of residues prone to oxidation, this phenomenon was successfully prevented. Importantly, the stability to oxidation conveyed by the amino acid exchanges led to increased enzyme activity. However, the exchanges resulted in slightly modified positions of two of the four metal‐binding amino acids, thereby strongly impacting metal binding.

Co-expression of thermophilic pectinases in a single host for cost-effective pectin bioconversion into D-galacturonic acid

Co-expression of enzymes allow to produce multiple enzymes in a single host, representing a cost-effective alternative in biocatalytic processes which can be used for pectin bioconversion. Pectin-rich biomass is an abundant by-product from the fruit and sugar industries that is usually disposed in landfill or sold as a low value feedstock. The aim of this work was to co-express a thermophilic pectin methyl esterase (PME) and exo-polygalacturonases (exo-PGs) in a single host for pectin bioconversion into D-galacturonic acid (GalA) using different pectic substrates such as apple, citrus and sugar beet pectin. To achieve this, a PME from Bacillus licheniformis (BLI09) with either an exo-PG from Thermotoga maritima (TMA01) or from Bacillus licheniformis (BLI04) were cloned in pETDuet-1 and co-expressed in E. coli BL21 (DE3). Four co-expression plasmids containing both pectinases were constructed and factors such as the effect of the genes’ cloning order and their expression were evaluated. Co-expression constructs 3 and 4 (pETDuet-TMA01-BLI09 and pETDuet-BLI04-BLI09, respectively) showed better expression of both pectinases compared to co-expression constructs 1 and 2 (pETDuet-BLI09-TMA01 and pETDuet-BLI09-BLI04, respectively). Co-expression constructs 3 and 4 were the most efficient for pectin bioconversion into GalA reaching 3 and 2.5 mM GalA, respectively from apple and citrus pectin after 4 h reaction. In conclusion, this work demonstrates that the co-expression of pectinases can potentially contribute to reduce the cost associated to their production and purification as well as to increase their applicability for exploiting pectin-rich biomass to obtain bio-based chemicals.

Regulatory mechanisms of c-di-AMP synthase from Mycobacterium smegmatis revealed by a structure: Function analysis.

Cyclic-di-nucleotide-based secondary messengers regulate various physiological functions, including stress responses in bacteria. Cyclic diadenosine monophosphate (c-di-AMP) has recently emerged as a crucial second messenger with implications in processes including osmoregulation, antibiotic resistance, biofilm formation, virulence, DNA repair, ion homeostasis, and sporulation, and has potential therapeutic applications. The contrasting activities of the enzymes diadenylate cyclase (DAC) and phosphodiesterase (PDE) determine the equilibrium levels of c-di-AMP. Although c-di-AMP is suspected of playing an essential role in the pathophysiology of bacterial infections and in regulating host-pathogen interactions, the mechanisms of its regulation remain relatively unexplored in mycobacteria. In this report, we biochemically and structurally characterize the c-di-AMP synthase (MsDisA) from Mycobacterium smegmatis. The enzyme activity is regulated by pH and substrate concentration; conditions of significance in the homoeostasis of c-di-AMP levels. Substrate binding stimulates conformational changes in the protein, and pApA and ppApA are synthetic intermediates detectable when enzyme efficiency is low. Unlike the orthologous Bacillus subtilis enzyme, MsDisA does not bind to, and its activity is not influenced in the presence of DNA. Furthermore, we have determined the cryo-EM structure of MsDisA, revealing asymmetry in its structure in contrast to the symmetric crystal structure of Thermotoga maritima DisA. We also demonstrate that the N-terminal minimal region alone is sufficient and essential for oligomerization and catalytic activity. Our data shed light on the regulation of mycobacterial DisA and possible future directions to pursue.

Design of an artificial phage-display library based on a new scaffold improved for average stability of the randomized proteins

Scaffold-based protein libraries are designed to be both diverse and rich in functional/folded proteins. However, introducing an extended diversity while preserving stability of the initial scaffold remains a challenge. Here we developed an original approach to select the ensemble of folded proteins from an initial library. The thermostable CheY protein from Thermotoga maritima was chosen as scaffold. Four loops of CheY were diversified to create a new binding surface. The subset of the library giving rise to folded proteins was first selected using a natural protein partner of the template scaffold. Then, a gene shuffling approach based on a single restriction enzyme was used to recombine DNA sequences encoding these filtrated variants. Taken together, the filtration strategy and the shuffling of the filtrated sequences were shown to enrich the library in folded and stable sequences while maintaining a large diversity in the final library (Lib-Cheytins 2.1). Binders of the Oplophorus luciferase Kaz domain were then selected by phage display from the final library, showing affinities in the μM range. One of the best variants induced a loss of 92% of luminescent activity, suggesting that this Cheytin preferentially binds to the Kaz active site.

Reduction of Nitric Oxide to Nitrous Oxide in Flavodiiron Proteins: Catalytic Mechanism and Plausible Intermediates

The flavin dependent nonheme diiron proteins comprise a family of enzymes, which can act as scavengers for both molecular oxygen and nitric oxide. The reduction of nitric oxide to nitrous oxide and water in flavodiiron proteins (FDPs) has been studied both experimentally and computationally, but the reaction mechanism is far from well understood. From experiments, it is known that two NO molecules can bind to the reduced active site, forming an observable diferrous dinitrosyl complex. A main question has been whether nitrous oxide can be formed directly from the diferrous dinitrosyl complex or if further reduction and/or protonation is needed to make this step feasible. Experiments have shown that nitrous oxide can be formed in a deflavinated form of the enzyme, indicating that further reduction is not needed. In the present study, hybrid density functional theory calculations are performed on a cluster model of the Thermotoga maritima FDP active site. We show that nitric oxide can be reduced to nitrous oxide and water using a direct coupling mechanism, i.e., without further additions to the reduced active site. The diferrous dinitrosyl complex can form an unstable N–N bridging hyponitrite intermediate, which can rotate into an N–O bond bridging hyponitrite with a low barrier. From this intermediate, the N–O bond cleavage leading to release of nitrous oxide is energetically feasible. An energy profile for the entire catalytic cycle of such a direct coupling mechanism is presented, and it is shown that the suggested mechanism agrees with data on FDP variants. Finally, an energy profile for the entire process starting with the fully reduced enzyme turning over four NO equivalents is constructed. This energy profile suggests explanations to experimentally observed states, such as the dihydroxyl form of the fully oxidized diferric state, and the difference with respect to returning to the original oxidized state after NO reduction between the flavinated and the deflavinated form of the enzyme.

Redox tuning of the H-cluster by second coordination sphere amino acids in the sensory [FeFe] hydrogenase from Thermotoga maritima.

[FeFe] hydrogenases are exceptionally active catalysts for the interconversion of molecular hydrogen with protons and electrons. Their active site, the H-cluster, is composed of a [4Fe-4S] cluster covalently linked to a unique [2Fe] subcluster. These enzymes have been extensively studied to understand how the protein environment tunes the properties of the Fe ions for efficient catalysis. The sensory [FeFe] hydrogenase (HydS) from Thermotoga maritima has low activity and displays a very positive redox potential for the [2Fe] subcluster compared to that of the highly active prototypical enzymes. Using site directed mutagenesis, we investigate how second coordination sphere interactions of the protein environment with the H-cluster in HydS influence the catalytic, spectroscopic and redox properties of the H-cluster. In particular, mutation of the non-conserved serine 267, situated between the [4Fe-4S] and [2Fe] subclusters, to methionine (conserved in prototypical catalytic enzymes) gave a dramatic decrease in activity. Infra-red (IR) spectroelectrochemistry revealed a 50 mV lower redox potential for the [4Fe-4S] subcluster in the S267M variant. We speculate that this serine forms a hydrogen bond to the [4Fe-4S] subcluster, increasing its redox potential. These results demonstrate the importance of the secondary coordination sphere in tuning the catalytic properties of the H-cluster in [FeFe] hydrogenases and reveal a particularly important role for amino acids interacting with the [4Fe-4S] subcluster.

Identification of Matriglycan by Dual Exoglycosidase Digestion of α-Dystroglycan.

Matriglycan is a linear polysaccharide of alternating xylose and glucuronic acid units [-Xyl-α1,3-GlcA-β1,3]n that is uniquely synthesized on α-dystroglycan (α-DG) and is essential for neuromuscular function and brain development. It binds several extracellular matrix proteins that contain laminin-globular domains and is a receptor for Old World arenaviruses such as Lassa Fever virus. Monoclonal antibodies such as IIH6 are commonly used to detect matriglycan on α-DG. However, endogenous expression levels are not sufficient to detect and analyze matriglycan by mass spectrometry approaches. Thus, there is a growing need to independently confirm the presence of matriglycan on α-DG and possibly other proteins. We used an enzymatic approach to detect matriglycan, which involved digesting it with two thermophilic exoglycosidases: β-Glucuronidase from Thermotoga maritima and α-xylosidase from Sulfolobus solfataricus. This allowed us to identify and categorize matriglycan on α-DG by studying post-digestion changes in the molecular weight of α-DG using SDS-PAGE followed by western blotting with anti-matriglycan antibodies, anti-core α-DG antibodies, and/or laminin binding assay. In some tissues, matriglycan is capped by a sulfate group, which renders it resistant to digestion by these dual exoglycosidases. Thus, this method can be used to determine the capping status of matriglycan. To date, matriglycan has only been identified on vertebrate α-DG. We anticipate that this method will facilitate the discovery of matriglycan on α-DG in other species and possibly on other proteins. Key features • Analysis of endogenous matriglycan on dystroglycan from any animal tissue. • Matriglycan is digested using thermophilic enzymes, which require optimum thermophilic conditions. • Western blotting is used to assay the success and extent of digestion. • Freshly purified enzymes work best to digest matriglycan.

Modification and Production of Encapsulin.

Encapsulins are a class of protein nanocages that are found in bacteria, which are easy to produce and engineer in E. coli expression systems. The encapsulin from Thermotoga maritima (Tm) is well studied, its structure is available, and without modification it is barely taken up by cells, making it promising candidates for targeted drug delivery. In recent years, encapsulins are engineered and studied for potential use as drug delivery carriers, imaging agents, and as nanoreactors. Consequently, it is important to be able to modify the surface of these encapsulins, for example, by inserting a peptide sequence for targeting or other functions. Ideally, this is combined with high production yields and straightforward purification methods. In this chapter, we describe a method to genetically modify the surface of Tm and Brevibacterium linens (Bl) encapsulins, as model systems, to purify them and characterize the obtain nanocages.

Abstract 2570: Structures and kinetics of Thermotoga maritima MetY reveal insights into the predominant route for methionine biosynthesis in bacteria

Abstract 2570: Structures and kinetics of Thermotoga maritima MetY reveal insights into the predominant route for methionine biosynthesis in bacteria

Topochemical Design of Cellulose-Based Carriers for Immobilization of Endoxylanase

Xylooligosaccharides (XOSs) gained much attention for their use in food and animal feed, attributed to their prebiotic function. These short-chained carbohydrates can be enzymatically produced from xylan, one of the most prevalent forms of hemicellulose. In this work, endo-1,4-β-xylanase from Thermotoga maritima was immobilized on cellulose-based beads with the goal of producing xylooligosaccharides with degrees of polymerization (DPs) in the range of 4-6 monomeric units. More specifically, the impact of different spacer arms, tethers connecting the enzyme with the particle, on the expressed enzymatic activity and oligosaccharide yield was investigated. After surface functionalization of the cellulose beads, the presence of amines was confirmed with time of flight secondary ion mass spectrometry (TOF-SIMS), and the influence of different spacer arms on xylanase activity was established. Furthermore, XOSs (DPs 2-6) with up to 58.27 mg/g xylan were obtained, which were greatly enriched in longer oligosaccharides. Approximately 80% of these XOSs displayed DPs between 4 and 6. These findings highlight the importance of topochemical engineering of carriers to influence enzyme activity, and the work puts forward an enzymatic system focusing on the production of longer xylooligosaccharides.

Solubility and Thermal Stability of Thermotoga maritima MreB.

The basis of MreB research is the study of the MreB protein from the Thermotoga maritima species, since it was the first one whose crystal structure was described. Since MreB proteins from different bacterial species show different polymerisation properties in terms of nucleotide and salt dependence, we conducted our research in this direction. For this, we performed measurements based on tryptophan emission, which were supplemented with temperature-dependent and chemical denaturation experiments. The role of nucleotide binding was studied through the fluorescent analogue TNP-ATP. These experiments show that Thermotoga maritima MreB is stabilised in the presence of low salt buffer and ATP. In the course of our work, we developed a new expression and purification procedure that allows us to obtain a large amount of pure, functional protein.

Knotting Optimization and Folding Pathways of a Go-Model with a Deep Knot.

Formation of protein knots is an intriguing offshoot of the protein folding problem. Since experimental resolution on knot formation is limited, theoretical methods currently provide the most detailed insights into the knotting process. While suitable for shallow knots, molecular dynamics simulations have faced challenges capturing the formation of deep knots in proteins such as the minimally tied trefoil α/β methyltransferase from Thermotoga maritima (MTTTM). To improve the efficiency of MTTTM knotting in Cα Go-model simulations, mutant variants of the MTTTM Go-model were investigated. Through a structure-based analysis of knotted and unknotted states, four residues (K71, R72, E75, V76) were identified to increase the knotting efficiency from 2% to 83% when their contact energies were doubled and dihedral strength around the knot loop increased. The key features of this model are (i) a C-terminal slipknot intermediate that threads the knot in a highly unstructured intermediate, (ii) the inability to knot in native-like intermediate states, and (iii) a minor population in a long-lived trap that cannot knot. Examination of residue 71-76 contacts provides a small set of potential mutants that can directly test the model's validity. In addition, the knotting optimization process developed here has broad applicability in generating knotting-efficient models of other knotted proteins.