Loading

L o a d i n g ,   p l e a s e   w a i t . . .

Engineering

 Overview

In our project, we conducted two engineering efforts: we investigated the optimal reaction conditions and physiological properties of PET05 and 40 new enzymes, and completed the modification of PET05 for enhancing its enzymatic activity and thermostability.

 Engineering 1 ——
Optimization of Reaction Conditions for PET05 and 40 New Enzymes

1.1 Optimization of Reaction Conditions for PET05

Cycle 1:

Design:

First, we designed the following orthogonal experiment:

pH NaCl(M) Enzyme Loading (nM) Temperature(℃)
8 9 10 0 1 2 3 4 5 25 50 100 30 50
- + - + - - - - - - + - - +
- + - - + - - - - - + - - +
- + - - - + - - - - + - - +
- + - - - - + - - - + - - +
- + - - - - - + - - + - - +
- + - - - - - - + - + - - +
+ - - - + - - - - - + - - +
- - + - + - - - - - + - - +
- + - - + - - - - + - - - +
- + - - + - - - - - - + - +
- + - - + - - - - - + - + -

Build:

The plasmid vector we used was pET-32a(+), and proteins expressed from this vector carry a TrxA tag at the C-terminus. Based on preliminary experiments, we found that PET05 with the TrxA tag had no enzymatic activity. Additionally, it is known from literature that PET05 is highly soluble; therefore, we performed tag removal on PET05. We successfully purified PET05 without the TrxA tag and used NanoDrop to determine the concentration of the purified protein solution. Using the Expasy - ProtParam website, we calculated that the relative molecular mass of PET05 is 27724.34 Da, the extinction coefficient is 1.431, and the isoelectric point is 5.52.
According to the formula:

Mass (ng) = Enzyme Loading (nM) × Reaction Volume (μL) × Relative Molecular Mass of the Enzyme (g/mol) × 1000

we can calculate the mass of PET05 to be added in a 600 μL reaction system.According to the formula:

Volume of Enzyme Solution to be Added (μL) = Mass (ng) × Extinction Coefficient / Concentration of the Enzyme Solution (ng/μL)

we can calculate the volume of the PET05 enzyme solution to be added.

Test:

We completed 11 groups of orthogonal experiments with 3 parallel replicates in each group, and determined the release amounts of the products TPA and MHET after 18 hours of reaction via ultra performance liquid chromatography (UPLC). The detailed experimental data results can be found in here.

Learn:

We found that PET05 only exhibits enzymatic activity in buffer solutions with a NaCl concentration of 4-5 M.

Cycle 2:

Design:

Next, we increased the overall NaCl concentration of the reaction system and designed the following second round of orthogonal experiment:

pH NaCl(M) Enzyme Loading (nM) Temperature
8 9 10 4 4.5 5 45 50 55 60 50 500 1000
+ - - - - + - + - - - + -
- + - - - + - + - - - + -
- - + - - + - + - - - + -
- + - + - - - + - - - + -
- + - - + - - + - - - + -
- + - - - + + - - - - + -
- + - - - + - - + - - + -
- + - - - + - - - + - + -
- + - - - + - + - - + - -
- + - - - + - + - - - - +

Build:

Same as the first round.

Test:

We completed 10 groups of orthogonal experiments with 3 parallel replicates in each group. The treatment method is the same as the first round. The detailed experimental data results can be found in here.

Learn:

After two rounds of orthogonal experiments, we concluded that the optimal reaction conditions for PET05 are as follows:the optimal pH is 9, the optimal NaCl concentration is 5 M, the optimal reaction temperature is 55 °C, and saturation of enzyme amount is achieved when the Enzyme Loading is 50 nM.

1.2 Testing and Optimization of Reaction Conditions for 40 New Enzymes

Design:

To ensure that a comparison of properties can be made between our 40 new enzymes and PET05, we first designed an orthogonal experiment identical to the first round for PET05:

pH NaCl(M) Enzyme Loading (nM) Temperature
8 9 10 0 1 2 3 4 5 25 50 100 30 50
- + - + - - - - - - + - - +
- + - - + - - - - - + - - +
- + - - - + - - - - + - - +
- + - - - - + - - - + - - +
- + - - - - - + - - + - - +
- + - - - - - - + - + - - +
+ - - - + - - - - - + - - +
- - + - + - - - - - + - - +
- + - - + - - - - + - - - +
- + - - + - - - - - - + - +
- + - - + - - - - - + - + -

Build:

We successfully purified the 40 new enzymes and constructed the reaction system following the same process as that used for PET05.

Test:

We completed 11 groups of orthogonal experiments with 3 parallel replicates in each group. The treatment method is the same as that for PET05's rounds. The detailed experimental data results can be found in here.

Learn:

Among all the enzymes, we identified four of them that exhibit PET-degrading enzyme activity:

GOMC ID Numbers and Optimal Conditions

Part Name Name/Short Description Part Type ID Number in GOMC Description
Optimum pH Optimal NaCl Concentration/M Optimum Temperature/℃
BBa_25ERRSY9 B6 Coding >SRR8062318_GL0674422 8 4 50
BBa_25NLL6X5 B7 Coding >SRR7479623_GL0026367 8 4 50
BBa_25A3MGA6 B8 Coding >SRR100298_GL0002528 3
BBa_25P0HHP3 B14 Coding >SRR6207735_GL6219842 8 3 50

* Due to its relatively low melting temperature (Tm), the enzymatic activity of B8 decreases sharply with temperature changes. Consequently, its optimal NaCl concentration and optimum temperature were not determined in the experiment.

In Engineering 1, we found that among all the marine-derived PET-degrading enzymes we studied, PET05 exhibits the best properties. Both its enzymatic activity and thermostability are higher than those of the newly identified enzymes. Thus, we decided to select PET05 as the target for engineering modification.

 Engineering 2 ——
Engineering Modification of PET05

2.1 Mutation of the Enzyme Active Site Pocket Based on Structural Homology and Molecular Docking Results

Cycle 1:Single-point Mutation

Design:

To directly enhance the catalytic activity of 05PET enzyme, we successively utilized structural homology and molecular docking to locate its pocket, selected hot spot amino acids and conducted mutation and wet experiments for verification.

We first used AutoDock Vina to perform molecular docking of 05PET enzyme with trimeric PET long-chain molecules. By comparing with the structure of LCC, we excluded some non-active pocket areas. Based on this, we delineated a more reliable docking region and re-performed molecular docking. After identifying the sequence sites corresponding to the active pocket space, we took 05PET enzyme as the basis and conducted mutations based on the sequences of other PET enzymes, determining the following single-point mutation schemes:

G55A T56A T56Y T56F T57A T58G T62S T62A Y63I Y63L Y63A H123A H123W S124A Q125M Q125A Y148A Y148W I150A I150T I150L I150S G151A G152A D174A I176A I186V L199A L199I S200N S200A S200D S200R G201A G201N A202G S203A S203T H204A F205A F205L F205S E206A E206C E206V V208N V208A G209A G209I G209T G209S S200A S200N S200G S200P A211N A211D G212A G212T G212K G212S D213A D213I D213N D213T D234A D242A D246A D248A D248E

Build:

We designed forward and reverse primers for each mutant, and performed PCR using 32a-PET05 (i.e., the wild-type PET05 gene constructed in the pET-32a(+) plasmid vector) as the PCR template. The PCR products were subjected to agarose gel electrophoresis and gel extraction. The gel-extracted products were then subjected to homologous recombination. The homologous recombination products were used for slow transformation of E. coli DH5α. After transformation, the Ampicillin-resistant (Amp-resistant) plate medium was incubated at 37°C for 12 hours, and then sent to a sequencing company for sequencing. The sequencing and assembly results were checked. The detailed construction procedure can be found in notebook here.

Plasmids with correct sequencing results were returned to the laboratory by the company. We transformed the plasmids into E. coli BL21(DE3), followed by bacterial culture and protein purification. We determined the concentration (mg/mL) of the purified protein solution and the A260/280 ratio.

Test:

We constructed the same reaction system for both PET05 mutants and PET05 as follows:
Reaction buffer (600 μL): 50 mM Glycine, 5 M NaCl, pH = 9
Reaction time: 18 hours
Reaction temperature: 55 ℃
Rotation speed: 300 rpm
Reactant: flaky lcPET (with a uniform diameter)

According to the formula:
Mass (ng) = Enzyme Loading (nM) × Reaction Volume (μL) × Relative Molecular Mass of the Enzyme (g/mol) × 1000
we can calculate the mass of the enzyme solution to be added in a 600 μL reaction system.
According to the formula:
Volume of the Enzyme Solution to be Added (μL) = Mass (ng) × Extinction Coefficient / Concentration of the Enzyme Solution (ng/μL)
we can calculate the volume of the enzyme solution to be added.

We completed experiments for several mutants and PET05, with 3 parallel replicates in each group. The release amounts of the products TPA and MHET after 18 hours of reaction were determined via ultra performance liquid chromatography (UPLC). The detailed experimental data results can be found in here.

Learn:

After testing, we identified 10 mutants with higher enzymatic activity than the wild-type:

I150A I150T I150L I150S I176A L199A S200A S200N S203T F205A

We will combine these mutations for the next round of engineering cycle.

Cycle 2:Double-Point Combinatorial Mutation

Design:

Based on the 70 sequences in section 2.3, we selected 10 mutants with enzymatic activity higher than that of the wild-type PET05. These 10 mutants were combined in pairs to form 38 types of double-point combinatorial mutants. We then designed the PCR templates and primers required for constructing these mutants:

Mutants, Templates, and Corresponding Primers

Name Template Mutant primer (F) Mutant primer (R) Mutant
D1 32a-PET05-I150A I176A-F I176A-R I150A-I176A
D2 32a-PET05-I150A L199A-F L199A-R I150A-L199A
D3 32a-PET05-I150A S200A-F S200A-R I150A-S200A
D4 32a-PET05-I150A S200N-F S200N-R I150A-S200N
D5 32a-PET05-I150A S203T-F S203T-R I150A-S203T
D6 32a-PET05-I150A F205A-F F205A-R I150A-F205A
D7 32a-PET05-I150T I176A-F I176A-R I150T-I176A
D8 32a-PET05-I150T L199A-F L199A-R I150T-L199A
D9 32a-PET05-I150T S200A-F S200A-R I150T-S200A
D10 32a-PET05-I150T S200N-F S200N-R I150T-S200N
D11 32a-PET05-I150T S203T-F S203T-R I150T-S203T
D12 32a-PET05-I150T F205A-F F205A-R I150T-F205A
D13 32a-PET05-I150L I176A-F I176A-R I150L-I176A
D14 32a-PET05-I150L L199A-F L199A-R I150L-L199A
D15 32a-PET05-I150L S200A-F S200A-R I150L-S200A
D16 32a-PET05-I150L S200N-F S200N-R I150L-S200N
D17 32a-PET05-I150L S203T-F S203T-R I150L-S203T
D18 32a-PET05-I150L F205A-F F205A-R I150L-F205A
D19 32a-PET05-I150S I176A-F I176A-R I150S-I176A
D20 32a-PET05-I150S L199A-F L199A-R I150S-L199A
D21 32a-PET05-I150S S200A-F S200A-R I150S-S200A
D22 32a-PET05-I150S S200N-F S200N-R I150S-S200N
D23 32a-PET05-I150S S203T-F S203T-R I150S-S203T
D24 32a-PET05-I150S F205A-F F205A-R I150S-F205A
D25 32a-PET05-I176A L199A-F L199A-R I176A-L199A
D26 32a-PET05 L199A-S200A-F L199A-S200A-R L199A-S200A
D27 32a-PET05 L199A-S200N-F L199A-S200N-R L199A-S200N
D28 32a-PET05 L199A-S203T-F L199A-S203T-R L199A-S203T
D29 32a-PET05 L199A-F205A-F L199A-F205A-R L199A-F205A
D30 32a-PET05-I176A S200A-F S200A-R I176A-S200A
D31 32a-PET05 S200A-S203T-F S200A-S203T-R S200A-S203T
D32 32a-PET05 S200A-F205A-F S200A-F205A-R S200A-F205A
D33 32a-PET05-I176A S200N-F S200N-R I176A-S200N
D34 32a-PET05 S200N-S203T-F S200N-S203T-R S200N-S203T
D35 32a-PET05 S200N-F205A-F S200N-F205A-R S200N-F205A
D36 32a-PET05-I176A S203T-F S203T-R I176A-S203T
D37 32a-PET05 S203T-F205A-F S203T-F205A-R S203T-F205A
D38 32a-PET05-I176A F205A-F F205A-R I176A-F205A

Build:

The method is the same as Cycle 1.

Test:

The method is the same as Cycle 1.

Learn:

We have completed the construction and testing of a total of 16 double mutants:The mutants D2, D3, D4, D5, D6, D9, D11, D14, D16, D17, D21, D25, D30, D33, D36, and D38 exhibit significantly higher enzymatic activity than the wild-type. Among them, the D5 mutant achieves a product release amount of 1.28 μM within 18 hours, representing a 62% increase compared to wild-type PET05.



2.2 Enzyme Thermostability Engineering Based on Deep Learning and Physicochemical Algorithms

Cycle 1 :Single-point mutation

Design:

Overall process: ThermoMPNN->AlphaFold->FoldX

Therefore, we first performed structure prediction on the 05PET enzyme using AlphaFold, and then ran ThermoMPNN with the predicted structure of PET05 as input. After obtaining a 20×257 heatmap and its corresponding CSV file, we screened out mutations with a model-predicted ddG (delta delta Gibbs free energy) greater than 0, and used AlphaFold to conduct structure prediction for these single-point mutations. Subsequently, FoldX, a software capable of rapidly calculating protein stability, was employed to sequentially perform "Repair" operations and predict the folding free energy of both PET05 and all its mutants. Through this process, ten single-point mutations with a ddG of less than -6 kcal/mol were identified.

To prevent the mutants from affecting the catalysis of the enzyme's active pocket, we first determined the active pocket position of LCC (a reference enzyme) and then identified the catalytic pocket of PET05 by leveraging homology analysis. Molecular docking within the catalytic pocket of the 05PET enzyme was carried out using AutoDock Vina, and the single-point mutation sites identified in the previous step were further verified based on the set of amino acid residues involved in the interaction with the substrate. Finally, we screened out 9 single-point mutations:

T142L V208L D213R D234K Q146F D234S D242R D248M D248K

Disulfide bonds represent a potential structural feature capable of enhancing enzyme stability. To further improve the stability of the PET05 enzyme, we employed Disulfide by Design 2.0 to engineer disulfide bonds into PET05. Disulfide by Design 2.0 is an algorithm that utilizes protein structures to design disulfide bonds. We executed the algorithm using the structure of PET05 as input. The program output possible disulfide bond mutations along with associated metrics. By applying the threshold criteria recommended in the literature, we selected two triple-disulfide bond variants and utilized AlphaFold for structural prediction to validate the feasibility of forming three disulfide bonds. The simultaneous presence of all three disulfide bonds was observed in one structural prediction outcome, confirming the viability of the disulfide bond design. The three disulfide bonds are as follows:

A58C-T62C Q159C-R189C A121C-A143C

Build:

The method is the same as Cycle 1 in Section 2.1

Test:

The method is the same as Cycle 1 in Section 2.1

Learn:

After testing, we excluded the Q146F and T142L mutants, whose activity was significantly reduced. We then determined the Tm values of the other mutants to rule out any impact on thermostability.



Cycle 2 :Composite mutant

Design:

After ruling out potential impacts on the active pocket, we combined the single-point mutations, predicted their structures following the same protocol as described in section 2.1, performed repairs, and calculated the folding free energies. Finally, we integrated both the single-point mutations and the three disulfide bonds into the sequences simultaneously, generating multiple sequence variants for subsequent wet laboratory experiments.

Combination Mutants and Their Designations

Name Mutant
C1 A58C_T62C_Q159C_R189C_
C2 q159w_v208l_
C3 v208l_d213r_
C4 d213r_d234k_
C5 d234k_d242r_
C6 d242r_d248k_
C7 q159w_v208l_d213r_d234k_d242r_d248k_A58C_T62C_Q159C_R189C_

Build:

The method is the same as Cycle 1 in Section 2.1

Test:

The method is the same as Cycle 1 in Section 2.1

Learn:

Email copied! Paste into your email client: tjusls_2025china@163.com