Outline
We designed a system combined antimicrobial peptides and microneedles, hoping to cure skin diseases at an early stage. After cycles of brainstorming and modifying, we figure out the best way to reach our expectations. In fact, the engineering cycle helps us a lot to put our ideas into practice and the cycle can be found in every part of our project.
Fig.1 shows that our experiment went through five stages and the Model has four stages in the project. In the end, we combined them to make our project much better.
Fig.1 From the left, experiment , model engineering cycle
Experiment
Our project has two parts: AMPC and MN. This part mainly describes the evolution of AMPC, and Fig.2 is the evolution diagram of the microneedle. For the detailed evolution process, please visit Hardware.
Fig.2 The evolution of microneedles
Design
Aiming at specific detection and sterilization of S. aureus, we designed an antimicrobial peptide complex using the V8 protease cleavage site as a linker between our antimicrobial peptides. V8 protease is an extracellular protease of S. aureus that cleaves peptide bonds exclusively on the aspartate (D) and glutamic acid (E) amino residues. In our antimicrobial peptide complex, we combine one kind of antimicrobial peptide. The main antimicrobial peptide we planned to use the human AMPs, such as hBD-3, HNP-1, or LL-37.
Redesign
However, IDT company can’t produce repetitive sequences. Interestingly, we found that two different AMPs have the synergistic effect for better sterilization in the literature, so we redesign that with two different AMPs. Besides, we selected an AMP called lysostaphin to be the main AMP component, which can hydrolyze the cell wall of S. aureus. We also combined it with other AMP, at this time, not only hBD-3, HNP-1, and LL-37, but Histatin 5 and Ranalexin became the AMP component. Different from the first design, we tried to use different AMP and nucleic acid to increase the randomness, making our sequences work.
Fig.3 The evolution of AMPC design (a) first designed AMPC (b) redesigned AMPC
Build & Test
We put the bacterium into 37°C incubator. However, we found that the protein expression was not obvious, so the sample was changed to be incubated at 16°C to facilitate stable expression of the protein. According to the SDS-PAGE analysis results, we divided protein expression into three categories, soluble , insoluble and not induced. Therefore, we try to purify the protein by the method of solubilization of inclusion body to break this barrier, so that it could reach the bacteria and function.
Learn
At the purification part, we tried hard to get the proteins and do the disc diffusion assays. We have tried eight proteins and found that protein No.3 and No.7 being the most effective. Next, we focus on protein No.3 and No.7, cultivating and purifying them repeatedly. After that, we adjusted their concentration and redesigned the experiment in order to enhance the antimicrobial effect.
See Result - Microneedle page and Result - AMPC page for detail information.
Model
Design
Our model purpose was to predict the effective time of our system. That is, we’re trying to see how long the AMP will be effective after the AMPC release from the microneedle and be cleaved by V8 protease. We’ve tried to get the reaction rate constant of the reaction of that AMPC to become single AMPs, the minimum inhibitory concentration of different single AMP, and assume a value to the degradation rate of a single AMP according to experiments. Last, we transfer the chemical reaction to a differentiated equation and use python code to draw the drug concentration-time curve.
We designed the formula and assumed two different reaction rate constants for the different site. we assumed k_1 is the first V8 protease cleavage site, and k_2 is the second one.
Preference sites are represented EL, ED, EF, and other sites are EG, EA, and EK.
(’’EL’’ means E and L are P1 and P1’ sites respectively.)
Next, we would need two more variables for our model. One is the MIC values of each AMP, the other one is the degradation rate of a single AMP.
Bulid
MIC Values
We found the minimum inhibitory concentration of each AMP, we check whether the MIC is reached to judge whether our drug is efficient.
AMP | Molecular Weight (Da) | Minimum Inhibitory Concentraion (MIC) |
---|---|---|
IcCCL28-25 | 8500 | 6.25-50 μg/mL(0.735-5.88μM) |
HNP1 | 3442 | 2.2 - 7.9 μg/mL(0.639-2.29μM) |
Ranalexin | 2618 | 4 - 16 mg/L(1.9-7.6μM) |
LL37 | 4493 | 2.9 - 3.6 μg/mL(0.645-0.801μM) |
hBD3 | 5157 | 100 μg/mL(19.39μM) |
Lysostaphin | 27000 | 0.003 - 2 μg/mL(0.0001-0.074μM) |
Histatin5 | 3037 | >400 μg/mL(>131.7μM) |
Table.1 AMPs’ molecular weight and minimum inhibitory concentration
Degradation Rate
Because we couldn’t find more related research about degradation rate, we refer to the research paper that LL37’s half-life time is 1 hour in cell as the AMPs half-life time. Then, set the parameter in the differential equation to draw the concentration-time curve.
Python Code
In fact, we’ve completed the code and tried to set the default value randomly based on the 2021 CCU_Taiwan team’s result and some papers before we get the real data of each reaction rate from the experiment to see if the result is reasonable or not. However, we think there’s another way to simulate realistic circumstances. We then talked with seniors and redesigned our model. Full source code can be found in our GitLab and more design details can be seen on the Model page.
Learn
Che-Kang Chang who is NTHU 2016 member advised us to apply the enzyme kinetics in our model, making our model more persuasive. We looked into some tutorials and examples to start redesigning our work.
Redesign
We modified our model design and applied the Enzyme Kinetics formula.
Enzyme Kinetics has four hypotheses:
- [ES] is constant. So k_1[E][S] = k_2[ES] + k_3[ES] always holds.
- All enzyme concentration [E_t] = exist alone [E_f] + [ES]
- The reaction rate is decided by k_3[ES].
- The maximum reaction rate V_max is the situation in which all the enzymes are turned into [ES]. V_max = k_3[E_t]
Keeping the MIC values and the assumption of degradation rate in the original version, we modified more in the reaction rate constant part among the three main variables in our model.
Bulid
We referred to the paper about V8 protease to set the k_3. Because [S] is much larger than [E], we consider that k_3 is much larger than k_2 and assume the k_2 value is 0. Therefore, we can get the k_1 value by definition K_max = (k2+k3)/k1.
Reaction rate const | Value(hr^-1) |
---|---|
k_1 | 0.037 |
k_2 | 0 |
k_3 | 1044 |
Table.2 Reaction rate constant we assumed.
In the new version of python code, the skeleton was basically the same but was changed the equations and had three reaction rate constants this time. Please check out the full source code in our GitLab.
Conclusion
According to the above, we knew the reaction rate is controlled by k_3[ES], and we got the k_3 value from the research paper. In order to reach our purpose, we hoped our system could release the drug slowly and stably. Therefore, we expect our system can reach the MIC of AMPs keeping for a few hours. By assuming the enzyme concentration, we estimate how much initial concentration of the drug we need to achieve the minimum inhibitory concentration, and estimate how much concentration of the drug we need.
See Model page for detail information.
Reference
- Domhan C, Uhl P, Kleist C, Zimmermann S, Umstätter F, Leotta K, Mier W, Wink M. Replacement of l-Amino Acids by d-Amino Acids in the Antimicrobial Peptide Ranalexin and Its Consequences for Antimicrobial Activity and Biodistribution. Molecules. 2019 Aug 17;24(16):2987. doi: 10.3390/molecules24162987. PMID: 31426494; PMCID: PMC6720431.
- Yang XY, Li CR, Lou RH, Wang YM, Zhang WX, Chen HZ, Huang QS, Han YX, Jiang JD, You XF. In vitro activity of recombinant lysostaphin against Staphylococcus aureus isolates from hospitals in Beijing, China. J Med Microbiol. 2007 Jan;56(Pt 1):71-76. doi: 10.1099/jmm.0.46788-0. PMID: 17172520.
- Singh D, Vaughan R, Kao CC. LL-37 peptide enhancement of signal transduction by Toll-like receptor 3 is regulated by pH: identification of a peptide antagonist of LL-37. J Biol Chem. 2014 Oct 3;289(40):27614-24. doi: 10.1074/jbc.M114.582973. Epub 2014 Aug 4. PMID: 25092290; PMCID: PMC4183800.
- Su J, Li H, Hu J, Wang D, Zhang F, Fu Z, Han F. LcCCL28-25, Derived from Piscine Chemokine, Exhibits Antimicrobial Activity against Gram-Negative and Gram-Positive Bacteria In Vitro and In Vivo. Microbiol Spectr. 2022 Jun 29;10(3):e0251521. doi: 10.1128/spectrum.02515-21. Epub 2022 May 26. PMID: 35616397; PMCID: PMC9241943.
- Wanmakok M, Orrapin S, Intorasoot A, Intorasoot S. Expression in Escherichia coli of novel recombinant hybrid antimicrobial peptide AL32-P113 with enhanced antimicrobial activity in vitro. Gene. 2018 Sep 10;671:1-9. doi: 10.1016/j.gene.2018.05.106. Epub 2018 May 30. PMID: 29859288.
- Komatsu T, Salih E, Helmerhorst EJ, Offner GD, Oppenheim FG. Influence of histatin 5 on Candida albicans mitochondrial protein expression assessed by quantitative mass spectrometry. J Proteome Res. 2011 Feb 4;10(2):646-55. doi: 10.1021/pr100861k. Epub 2010 Dec 6. PMID: 21080726; PMCID: PMC3033980.
- Midorikawa K, Ouhara K, Komatsuzawa H, Kawai T, Yamada S, Fujiwara T, Yamazaki K, Sayama K, Taubman MA, Kurihara H, Hashimoto K, Sugai M. Staphylococcus aureus susceptibility to innate antimicrobial peptides, beta-defensins and CAP18, expressed by human keratinocytes. Infect Immun. 2003 Jul;71(7):3730-9. doi: 10.1128/IAI.71.7.3730-3739.2003. PMID: 12819054; PMCID: PMC162002.