Part:BBa_K316012

TEV protease S219P autocatalysis resistant variant

TEV protease S219P autocatalysis resistant variant. This part had been reversed for the 3' strand in order to reduce any read-through that may be caused by upstream elements.

Introduction :

This is the nuclear inclusion protease, endogenous to Tobacco Etch Virus and is used in the late lifecycle to cleave polyprotein precursors. The recognition sequence is ENLYFQG/S 1 between QG or QSDue to it’s stringent sequence specificity, TEV is commonly used to cleave genetically engineered proteins.

Uses:

TEV proteinase is used to cleave fusion proteins. It is useful due to its high degree of specificity1 and potential to be used in vivo or in vitro applications.

Auto-inactivation:

Wild type TEV protease also cleaves itself at Met 218 and Ser 2192. This leads to auto-inactivation of the TEV protease and progressive loss of activity of the protein. The rate of inactivation is proportional to the concentration of protease. More stable Mutants have been produced by single amino acid substitutions S219V (AGC(serine) to GTG(valine) and S219P (AGC(serine) to CCG(proline)3.

Table I.

Kinetic parameters for wild-type and mutant TEV proteases with the peptide substrate TENLYFQSGTRR-NH2. From original paper by Kapust et.al. 20013

Enzyme	Km (mM)	kcat (s-1)	kcat /Km (mM-1 s-1)
Wild type	0.061 ± 0.010	0.16 ± 0.01	2.62 ± 0.46
S219V	0.041± 0.010	0.19 ± 0.01	4.63 ± 1.16
S219P	0.066 ± 0.008	0.09 ±0.01	1.36 ± 0.22

S219V* - retains same activity as wild type

S219P* - virtually imperivious to autocatalysis

Sequence and Features

Assembly Compatibility:

• 10
  COMPATIBLE WITH RFC[10]
• 12
  COMPATIBLE WITH RFC[12]
• 21
  COMPATIBLE WITH RFC[21]
• 23
  COMPATIBLE WITH RFC[23]
• 25
  INCOMPATIBLE WITH RFC[25]

  Illegal AgeI site found at 32
  
• 1000
  COMPATIBLE WITH RFC[1000]

References

<biblio>

1. 1 pmid=8179197
2. 2 pmid=7793070
3. 3 pmid=2047602

</biblio>

For more information about our project please visit our wiki.

Contribution

• Group: TecCEM, iGEM 2024
• Authors: Giovana Andrea Osorio León, Ana Laura Torres Huerta, Aurora Antonio Pérez, Lorena Gallegos Solís
• Summary: The sequences of TEV proteases S219P and 1LVM were compared to determine if TEV S219P exhibits greater efficiency and interaction with a substrate. Sequence alignment was performed using ESPript3, and both enzymes were modeled with AlphaFold for structural analysis. Additionally, heatmaps of protein stability were generated using Protein-Sol and molecular docking was carried out with HADDOCK 2.4 and BIOVIA Discovery Studio Visualizer to compare their interactions.

Documentation

The Tobacco etch virus (TEV) protease is a 27 kDa catalytic domain of the nuclear inclusion polyprotein a (NIa) in TEV. It specifically recognizes the amino acid sequence ENLYFQG/S and cleaves between the Q and G/S residues. TEV belongs to the Potyviridae family of viruses, which includes other positive-strand RNA viruses. The TEV genome is translated into a single large polyprotein that is subsequently processed by virus-encoded proteins with proteolytic activity. Despite its substrate specificity, the use of TEV protease is limited due to its self-inactivation through autocleavage and its low solubility during purification, caused by its high hydrophobicity .

Studies have shown that TEV protease undergoes self-cleavage near its C terminus after residue 218, at a site that does not follow the canonical sequence, resulting in a truncated protein with reduced activity . According to UniProt, its active sites are in 214, 223, 256, 649, 722, 2083, 2118, and 2188 .

In this study, the sequences of mutant TEV protease S219P and wild type 1LVM were compared to determine if TEV S219P exhibits greater efficiency and interaction with its recognition site or synthetic substrate (ENLYFQG). Structure-based sequence alignment was performed using ESPript3, and both enzymes were modeled with AlphaFold for structural analysis. Additionally, heatmaps of protein stability were generated using Protein-Sol and molecular docking was carried out with HADDOCK 2.4 and BIOVIA Discovery Studio Visualizer to compare their interactions. Finally, to corroborate the information obtained from modeling and docking we made our cloning design for TEV S219P in the plasmid pET-24b(+) using Benchling.

Initially, the sequences of TEV protease S219P and TEV 1LVM were analyzed using Clustal Omega to generate a protein alignment in ALN format, which was subsequently processed with ESPript3. ESPript3 provides a detailed representation of sequence similarities and secondary structure information for the aligned sequences. As shown in Figure 1, both sequences exhibit a high degree of similarity. However, the TEV S219P sequence consists of 237 amino acids, while the TEV 1LVM sequence contains 229 amino acids, indicating that TEV S219P is longer.

In panel A, differences are observed between amino acids 1-8 of TEV 1LVM and TEV S219P, attributed to the presence of a His-tag, as well as at position 10. Additionally, amino acids 229-237 in TEV S219P represent seven extra residues not present in TEV 1LVM. The differing amino acids between the two proteases are highlighted in yellow. The catalytic sites for both proteases are located at amino acids 214 and 223, as indicated by green triangles in panel A of Figure 1. The peptidase C4 domain, marked in dark blue, spans residues 1-218, corresponding to the coding sequence shared by both proteins.

In panel B, the domains of chain A of TEV 1LVM, selected for comparison with TEV S219P, are illustrated. The domains of 1lvmA01 are shown in aqua, while those of 1lvmA02 are highlighted in pink. Table 1 in Figure 1 provides a comprehensive description of these domains.

Figure 1. Proteín alignment of TEV 1LVM with TEV S219P and characteristics of TEV 1LVM domains in chain A. The structure of TEV S219P obtained from the iGEM registry (BBa_K316012) does not exhibit the S219P mutation as stated. Therefore, the subsequent analyses were conducted using the specified mutation. Image A retrieved from ESPript. [6]

To evaluate the impact of the differences in the sequence of the mutant with respect to the wild-type protein, modeling and simulation of the interaction with the target peptide were performed. Afterwards, the program ChimeraX was used to obtain the 3D model of both TEV proteases using the AlphaFold tool for structure prediction, as shown in Figure 2. The image labeled with letter A is the model obtained for the TEV protease S219P, and the image labeled with letter B is the model for chain A of TEV 1LVM. Both proteases were overlapped to compare them and observe similarities, as shown in image C. The overlap is highlighted in two colors: light blue for the TEV protease S219P and purple for TEV 1LVM, showing areas where the structures align or differ. It can be seen that most of the chains of both proteins align, as well as the main beta sheets and alpha helices, indicating that the S219P mutation does not significantly affect the native structure of the enzyme. However, there are regions at the ends where conformational differences between TEV S219P and TEV 1LVM are observed. These differences may influence the function or stability of the protein.

Figure 2. 3D models of TEV proteases obtained from AlphaFold. A) 3D model of TEV S219P. B) 3D model of TEV 1LVM. C) Overlap of both enzymes.

Subsequently, the Protein-Sol program was used to generate energy and charge heatmaps, providing insights into the stability of the proteases (TEV 1LVM and TEV S219P) under varying pH and ionic strength (salt concentration) conditions. The Debye-Hückel (DH) method was employed to model interactions between ionizable groups and to calculate pKa values. Each heatmap consists of 91 combinations of pH and ionic strength. For qualitative comparison with the experiment, the CH2 and CH3 domains of the IgG1 structure (PDB 1HZH) were used, as pH and ionic strength variations have been reported to affect the stability of these domains in IgG.

TEV 1LVM

The heatmap results for TEV 1LVM are shown in Figure 3. The color scale ranges from red (positive values) to green (negative values). Positive values (ranging from red to orange) suggest conditions where the protease has higher energy, potentially indicating instability or less favorable folding. Conversely, negative values (green) indicate more stable conditions with lower energy.

At a pH between 2.0 and 4.0 and high ionic strengths (0.1–0.3 M), the TEV 1LVM protease displays higher energy values, which may indicate enzyme instability. As the pH increases beyond 5.0 and the ionic strength decreases, the energy values drop, suggesting that these conditions favor a more stable folded state of the protease. However, the optimal conditions for this enzyme are observed at a pH of 6.5 to 8.0 and low ionic strengths (0.005–0.1 M), which promote the stability of the TEV 1LVM protease.

Figure 3. Energy heatmap for TEV 1LVM. Image obtained from Protein-Sol.[4]

Figure 4 displays the net charge per amino acid (e/aa) under various pH and ionic strength conditions. The color scale ranges from blue (positive charges) to red (negative charges). At low pH values (between 2.0 and 4.0), the protease exhibits a positive charge per amino acid. As the pH increases (between 5.0 and 8.0), the charges approach zero, particularly in the pH range of 6.0 to 8.0. This suggests that within these pH ranges, the protease is near its isoelectric point, where the net charge is essentially neutral. Figure 4 indicates that the TEV protease (1LVM) is more positively charged under acidic conditions (low pH) and approaches a neutral charge under more neutral to basic conditions (pH 6.0–8.0) .

Figure 4. Charge heatmap for TEV 1LVM. Image obtained from Protein-Sol [4].

TEV Protease S219P

Figure 5 shows that TEV S219P exhibits increased stability (negative energy) under pH conditions between 6.5 and 8.0, particularly at low ionic strengths (0.005–0.1 M). In contrast, at low pH (2.0–4.0), the enzyme displays very high energy values (deep red color) across the entire range of ionic strengths, indicating significant instability or a propensity for denaturation. This instability is attributed to the S219P mutation, where a serine (AGC) was replaced with a proline (CCG). Proline often has negative effects on protein structure because it is the only amino acid whose side chain is covalently bonded to the alpha amine, forming a rigid cyclic structure. This structure introduces rigid turns in the peptide chain, stabilizing the edges of beta-sheets and alpha-helices, thereby facilitating protein folding. Additionally, peptides with consecutive proline residues can fold into a characteristic polyproline helix (PII helix), a common motif in protein-protein interactions . However, this effect was not observed in the 3D model generated by AlphaFold, as both sequences are identical except for a single amino acid change and the addition of a histidine tag. Furthermore, homology modeling might not capture such stability changes. Several studies on peptide bond formation with proline suggest that proline can hinder the rate of protein synthesis by increasing the entropy of peptide bond formation. In general, the substitution of any amino acid can destabilize the native conformation of the protein , leading to its instability.

Figure 5. Energy heatmap for TEV S219P . Image obtained from Protein-Sol [4].

Figure 6 displays the net charge map for TEV S219P. Instability at low pH is shown, with high energies in this range. However, the magnitude of the energy values ​​in this TEV is lower compared to TEV 1LVM, which could indicate that TEV S219 is slightly more tolerant to acidic conditions. The charge of the protein decreases both with increasing pH and ionic strength, which is a typical behavior for many proteins. The smaller change in charge at higher pH suggests that residues affecting the net charge (probably basic groups) reach a point where their protonation is minimal .

Figure 6. Charge heatmap for TEV S219P . Image obtained from Protein-Sol [4].

Molecular docking to evaluate binding to the TEV recognition site

Molecular docking is a computational technique used to predict the interaction between two molecules, typically a small molecule (ligand) and a protein (receptor). The primary objective of molecular docking is to simulate the binding of a ligand to a specific active site of a protein and to predict the binding affinity and position. In this study, we evaluated the binding affinity of TEV protease S219P and TEV 1LVM with their synthetic substrate or recognition site (ENLYFQG). The docking simulations were performed using HADDOCK 2.4, resulting in the identification of the top 10 clusters, to observe the interaction probability in different random sites of the receptor. Subsequently, the PDB files of the best 5 clusters of each enzyme were analyzed using BIOVIA Discovery Studio Visualizer, yielding a total of 60 images from the 10 clusters of both TEV proteases. However, only the images corresponding to the first clusters were discussed, as outlined below .

In general terms, for both TEV proteases, the images labeled with letter A shows the interaction scheme between the ligand (polypeptide substrate, shown in aqua blue) and the receptor (TEV 1LVM in multiple colors according to secondary structure: cyan for beta-sheets, red for alpha-helices, and white for random coil). The amino acids involved in the interaction with the ligand are specifically colored (yellow-green/olive green). In images labeled with letter B, the interaction site was zoomed in. The substrate is shown in dark aqua blue, while the amino acids are depicted in olive green/yellowish green. The interactions are represented by dashed (or dotted) lines and may have different color codes, depending on the type of interaction established

On the other hand, in images labeled with letter C, the surface analysis of the receptor based on aromatic residues is shown. A color scale is used in the surface scheme, referencing 'nuclear' interactions referred to as 'Edge' or 'Face.' This does not pertain to a 'force' or 'magnitude,' but rather to whether the interaction occurs on the 'Face' of the aromatic ring plane, or on the 'Edge,' where one of the ring's carbons establishes the interaction. In images labeled with the letter D, the hydrogen bond interaction scheme on the receptor's surface is shown. The color scale indicates the detection and which amino acid of the receptor participates as a donor or acceptor. In images labeled with letter E, the surface charge interaction scheme is shown. In this scheme, the scale corresponds to both charge and magnitude. Finally, in images labeled with letter F, the surface interaction scheme of the receptor through hydrophobic interactions is shown. In this case, the scale is related to magnitude. The more negative (red/orange), the stronger the interaction.

TEV Protease 1LVM

In Figure 7, the results obtained in cluster 1 for TEV protease 1LVM for different types of schemes (A, B, C, D, E and F) can be seen. This cluster was chosen because it was the one that obtained the most interactions between ligand and receptor, giving a total of 36 interactions, as mentioned in Table 2. The more interactions there are between the ligand and the receptor, the greater the possibility of reaction. In panels A and B, these interactions can be seen, the green dotted line shows the conventional hydrogen bonds, while the orange dotted line shows an electrostatic bond.

On the other hand, panel C shows the surface analysis of the receptor by means of aromatic residues. The receptor is shown in different colors, but the most predominant are blue (Edge) and white. The "Edge" areas (blue) indicate that there are residues that interact more peripherally with the ligand, which could influence the orientation and conformation that the ligand adopts when binding to the receptor. Panel D shows the hydrogen bond interaction scheme on the surface of the receptor. The combination of colors (pink and green) shows that there are several areas of both hydrogen donation and acceptance, suggesting a network of hydrogen bond interactions, as shown in Table 2, since this type of bond is the most predominant in the interaction. This is essential to maintain the precise orientation of the ligand within the active site of the TEV protease. These interactions contribute to the affinity and specificity of the ligand binding to the protease.

The distribution of charges on the receptor surface and their complementarity with the ligand charges play a crucial role in molecular binding and recognition, as shown in panel E. The presence of oppositely charged regions (blue and red) can create an electrostatic environment that guides the ligand to its correct position within the binding site. The blue areas (positive charge) could be aligned to interact with negatively charged regions of the ligand, while the red areas (negative charge) would interact with positively charged regions of the ligand. It can be seen that neutral charges (0) and charges up to 0.100 predominate.

In panel F, the orange areas of the receptor indicate a high hydrophobicity, which tends to interact favorably with the hydrophobic regions of the ligand, contributing to the stability of the complex. The scale also allows to identify hydrophilic areas (blue), which could interact with polar or charged groups of the ligand, complementing the hydrophobic interactions.

Figure 7.Interactions between ligand and receptor of Cluster 1 TEV 1LVM, Image created in BIOVIA Discovery Studio Visualizer. A) Schematic of the interaction between ligand (polypeptide substrate) and receptor (TEV 1LVM). B) Zoom of the interaction site between ligand and receptor. C) Schematic of the surface analysis of the receptor by means of aromatic residues. D) Schematic of the interaction by hydrogen bonding on the surface of the receptor. E) Schematic of the interaction by charges on the surface of the receptor. F) Schematic of the interaction of the surface of the receptor by hydrophobic interactions.

TEV Protease S219P

In Figure 8, the results obtained in cluster 2 for TEV protease S219P for different types of schemes (A, B, C, D, E and F) can be seen. This cluster was chosen because it was the one that obtained the most interactions between ligand and receptor, giving a total of 18 interactions, as mentioned in Table 3, as well as a better binding strength (Gibbs energy). In panels A and B, the interactions of the ligand and the receptor are shown, compared to TEV 1LVM, it can be observed that there are mostly conventional hydrogen bonds, as well as an orange dotted line representing an electrostatic bond and the pink dotted line representing a hydrophobic bond (Pi-Pi in the form of a T). In panel C, it can be observed that the blue color predominates on the protease, which means that it has more "Edge" interactions with the aromatic rings of the substrate used. In panel D, it can be seen that the protease has more donor activity than acceptor activity, since it has more pink and white areas in its structure. Panel E shows the interpolar charges that exist within the protease; neutral charges (0) and up to 0.100 predominate in its structure, since white and blue regions are shown. Finally, in panel F, it can be seen that TEV S219P has as many hydrophilic (blue) or hydrophobic (orange) regions.

Figure 8. Interactions between ligand and receptor of Cluster 2 TEV protease S219P. Image created in BIOVIA Discovery Studio Visualizer.A) Schematic of the interaction between ligand (polypeptide substrate) and receptor (TEV S219P). B) Zoom of the interaction site between ligand and receptor. C) Schematic of the surface analysis of the receptor by means of aromatic residues. D) Schematic of the interaction by hydrogen bonding on the surface of the receptor. E) Schematic of the interaction by charges on the surface of the receptor. F) Schematic of the interaction of the surface of the receptor by hydrophobic interactions.

The interactions for the 5 clusters of each TEV protease are summarized in Table 2 and Table 3. The data indicate that TEV 1LVM exhibited the highest number of interactions, with a total of 105, in contrast to TEV S219P, which displayed a maximum of 64 interactions. This observation suggests that TEV 1LVM has a higher likelihood of achieving greater stability, as these interactions are crucial for maintaining the ligand in the correct orientation within the binding site. Notably, these interactions were concentrated in the catalytic sites of both enzymes and extended to adjacent residues. Similarly, TEV 1LVM, having more interactions, is more likely to allow the ligand to bind with higher affinity to the protein. Higher affinity indicates that the ligand can remain bound to the protein for a longer period, which can be essential for its functional efficacy. In the case of TEV S219P, a lower number of interactions results in weaker binding, potentially leading to a less efficient formation of the enzyme-substrate complex. For enzymes, strong and specific binding to their substrates is crucial for catalysis. The reduction in the number of interactions in the S219P variant suggests a lower binding affinity for its substrates, which could result in decreased enzymatic efficiency. The mutation in TEV S219P indicates that it may lead to a less stable binding and, consequently, a reduction in enzymatic activity compared to TEV 1LVM. In both Table 2 and Table 3, the most frequent type of interactions were conventional hydrogen bonds. In TEV 1LVM, a greater interaction was observed at active site 214 (HIS), as well as at nearby residues such as 213 (GLY) and 215 (LYS). Likewise, these residues within the receptor mostly behaved as H-Donor. Likewise, TEV S219P showed a low interaction at its active site 214 (GLY), but a greater affinity for residues such as 215 (HIS) and 216 (LYS), behaving as H-Donor.

Proposal for in silico Cloning Design in Expression Vector

To corroborate the information obtained from modeling and docking, experiments are necessary. We planned the cloning and expression of this mutant protease, however due to lack of time this was not possible. We show and make available our cloning design in case someone wants to clone this biobrick.

Benchling was used for in silico cloning of the mutant TEV protease previously mentioned, to insert its sequence into the pET-24b(+) vector (resistant to kanamycin). The steps followed for cloning are mentioned below.

First, the file of the plasmid pET-24b(+) was downloaded from Snapgene and then inserted into Benchling for its further modification. Also, the biobrick sequence was downloaded from the iGEM registry and copied into a new file in Benchling. Then, in the same biobrick file, the restriction sites for the enzymes NdeI and XhoI were added, as well as the extra base pairs that each enzyme required. For NdeI, in the 5’ extreme, we added 3 bases and also for XhoI, in the 3’ extreme, like it seems in Figure 9 and Figure 10. On the other hand, an additional 21 extra base pairs were added on each side of the sequence (5' and 3') to facilitate in silico cloning into the pET-24b(+) vector, as well as to increase the sequence length and qualify for the IDT sponsorship. Afterward, the modified TEV protease sequence was forward translated.

Figure 9. Addition of the restriction site for NdeI and the 24 extra pairs of bases for the in silico cloning.

Figure 10. Addition of the restriction site for Xho and the 24 extra pairs of bases for the in silico cloning.

Later, a forward and reverse primer were designed, with the following characteristics mentioned in Table 4.

Table 4. Characteristics of the forward and reverse primers.

	Sequence	GC (%)	Length (bp)	Min ΔG Homodimer (kcal)	Min ΔG Monomer (kcal)	Tm (°C)
Primer Forward	agccatatgaatgggaag	44.44	18	-6.4	-2.4	54.3
Primer Reverse	cacagttaatgaacdtcg	44.44	18	-3.9	-0.3	54.4

After we designed both primers, they were analyzed in the OligoAnalyzer tool for searching hetero-dimers, like it seems in Figure 11.

Figure 11. Hetero-dimers analyzer for both primers.

Next, an in silico digestion was performed with these enzymes in pET-24b(+) and in the file that contained the TEV protease S219P sequence. For both pET-24b(+) plasmid and construct, the buffer 2.1, 3.1, and 4/CS had an efficiency of 100%. After we performed both digests, with the assembly wizard tool, we made a ligation of pET-24b(+) with the TEV protease S219P sequence. Finally we confirmed that our cloning was successful by making a reverse translation on the open reading frame after the lac operator, where the amino acids of the modified TEV protease sequence were correctly translated and the histidine tag that will help us to purify this protease was conserved like the original plasmid before the cloning, like it seems in Figure 12.

Figure 12. Successful cloning file of TEV S219P in pET-24b(+).

The Benchling files and other details of the cloning proposal is explained in Table 5.

Table 5. Cloning proposal of TEV S219P in pET-24b(+).

Proposal	Description and link	Sequence	Length (bp)
TEV protease sequence ready to clone. TEV-Protease-BB	The cutting sites for the restriction enzymes XhoI and NdeI were added to the BioBrick sequence (BBa_K316012), as well as the extra base pairs that each enzyme occupies to cut correctly. In total, three extra bases were added to both the 5' and 3' extremes. Additionally, another 21 extra pairs of bases were added to facilitate in silico cloning into the pET-24b(+) vector, as well as to increase the sequence length and qualify for the IDT sponsorship. The insert is compatible with pET-24 b(+) and pET-29 a(+). Also, a forward and reverse primer were created for amplifying the TEV sequence by PCR. - Link to PDF of Original BioBrick in Benchling - Link to PDF of Modified BioBrick ready to clone in Benchling - Link to TEV protease cloned in pET24b(+) in Benchling	acgtggggcaaatttccctat agccatatga atgggaagatcgttgtttaaaggaccacgtgac tataatccgattagctcgactatt tgccatctga ccaatgagagtgatggtcataccactagcttgtat ggcattggctttgggccattcatcatcacgaacaa acacctctttaggcgcaataatggtacact gttgg tacaatcccttcatggagtctttaaggtcaaaaac acaacgacgcttcagcaacatct gatagatggaag agacatgattatcattcgaatgccgaaagactttc caccgtttcctcagaaactcaagtttcgcgaacct cagcgtgaagaacggatctgct tagtcacaacaaa ctttcagaccaaatctatgtcctcaatggtatcag acactagctgtacattccctag ctctgatggcatc ttttggaagcattggattcagacaaaagatgggca atgtggctctcctcttgtgtca acacgggatgggt ttattgtgggcatacactctgcgtcaaacttcacc aatacgaataattactttacgagtgttcccaagaa cttcatggagttactgacgaac caagaagctcaac aatgggtttcaggctggagactgaatgcagattcc gttctttggggaggtcacaaag tgttcatggataa accggaagaaccgtttcaaccggttaaagaggcca cacagttaatgaacCTCGAGagc ATCGCTAGGCAC ACGTCGCCG	772

The previous steps on how to make an in silico cloning of the TEV protease S219P into the plasmid or destiny vector pET-24b(+) are summarized in Figure 13.

Figure 13. Diagram of the steps followed to make an in silico cloning of the TEV protease S219P into the plasmid or destiny vector pET-24b(+).