Difference between revisions of "Team:NCTU Formosa/Model"

Line 607: Line 607:
 
         <h2 class="content-1" id="titleI" style="color:#FFFF82">II. The Life History of Moth</h2>
 
         <h2 class="content-1" id="titleI" style="color:#FFFF82">II. The Life History of Moth</h2>
 
         <div class="modelingPartContent" id="partI">
 
         <div class="modelingPartContent" id="partI">
−
           <p class="content">After understanding the life history of the moths from Dr. Huang in Taiwan Agriculture Research Institute(see more), we can know that the time of pupa becoming moth ranges from six days to 14 days. So the moths we caught in 20 days must be the larvae in our farm right now. Therefore, we utilized the data of the accumulated number of moth in the next 20 days as the target of our prediction system, and we can use the output to know the time that requires spraying Pantide for prevention.</p>
+
           <p class="content">After understanding the life history of the moths from Dr. Huang in Taiwan Agriculture Research Institute, we can know that the time of pupa becoming moth ranges from six days to 14 days. So the moths we caught in 20 days must be the larvae in our farm right now. Therefore, we utilized the data of the accumulated number of moth in the next 20 days as the target of our prediction system, and we can use the output to know the time that requires spraying Pantide for prevention.</p>
 
         </div>
 
         </div>
 
     </div>
 
     </div>
Line 639: Line 639:
 
         <h2 class="content-1" id="titleL" style="color:#FFFF82">I. Summary</h2>
 
         <h2 class="content-1" id="titleL" style="color:#FFFF82">I. Summary</h2>
 
         <div class="modelingPartContent" id="partL">
 
         <div class="modelingPartContent" id="partL">
−
           <p class="content"></p>
+
           <p class="content">The aim of this modeling was to predict and simulate the degradation rate of Pantide. Practically, the results would be integrated into our device to promote automatic control system. Once Pantide degraded below the effective level, it will spray the solution to replenish.</p>
 +
          <p class="content">After discussion and experimental verification, we show that UV radiolytic oxidation contributes to the degradation of Pantide. We tried to build a model for simulating the relationship between applied UV light and degradation rate and then combined with UV intensity sensor to calculate the period of spraying Pantide.</p>
 
         </div>
 
         </div>
 
     </div>
 
     </div>
Line 648: Line 649:
 
         <h2 class="content-1" id="titleM" style="color:#FFFF82">II. Pantide degradation process</h2>
 
         <h2 class="content-1" id="titleM" style="color:#FFFF82">II. Pantide degradation process</h2>
 
         <div class="modelingPartContent" id="partM">
 
         <div class="modelingPartContent" id="partM">
−
           <p class="content"></p>
+
           <p class="content">For the proteins used as Pantide, the inhibitor cystine knot (ICK) is significant to their function. Pantide with several disulfide bonds is often more stable in solution. If “native protein” is denatured, it loses disulfide bonds and becomes a less stable form, normally called “linear protein,” which is easily degradable.</p>
 +
          <p class="content">Recent research shows that the spider toxin proteins containing ICK structure, for example, ω-hexatoxin-Hv1a (Hv1a), have high stability against temperature, pH, solvents and protease. In contrast, when Hv1a is denatured to linear form, it loses its stability and then degrades rapidly. [1]</p>
 +
          <p class="content">There are many possible processes of Pantide degradation we discussed below. (Figure 1) Pantide may have a chance to be reduced to a linear form by reductants or reductases. For both native form and linear form proteins, it may suffer hydrolysis and proteolysis, resulting in denaturing or amino acid cleavage. Also, UV light of sun also leads to Pantide degradation. Though the energy of UV light may not be not enough to break the covalent bonds efficiently, proteins still could undergo radiolytic oxidation.</p>
 +
    <div>
 +
<img src=”” class=”picture”>
 +
<p class=”Figure 1. Pantide degradation process”></p>
 +
</div>
 +
          <p class="content">According to this degradation process, we attempted to build up our model. The degradation rate of Pantide is contributed by hydrolysis, proteolysis, UV radiolytic oxidation and the reduction to linear form. We considered that the rate Pantide transformed to linear form is a part of degradation rate because linear form protein loses its function and it is also likely to be degraded.</p>
 +
          <p class="content">However, since the whole process is too complex to verify, we divided the experiment into three parts, hydrolysis test, proteolysis test and UV radiolytic oxidation test, and further summarized the results to conclude. The purpose is to find out the degradation rate of Pantide and to verify the less stable linear form protein has.</p>
 +
          <p class="content">On the other hand, because ICK structure domain mainly contributes the stability of these proteins; we could assume that the degradation processes of three target proteins are roughly the same, for the reason that, we chose Hv1a and Hv1a with lectin (Hv1a-lectin) to demonstrate.</p>
 
         </div>
 
         </div>
 
     </div>
 
     </div>
 +
 +
 
<!--3-->
 
<!--3-->
 
     <div class="modelingPart">
 
     <div class="modelingPart">
 
         <h2 class="content-1" id="titleN" style="color:#FFFF82">III. Hydrolysis test</h2>
 
         <h2 class="content-1" id="titleN" style="color:#FFFF82">III. Hydrolysis test</h2>
 
         <div class="modelingPartContent" id="partN">
 
         <div class="modelingPartContent" id="partN">
 +
          <p class="content-1" style="color:#00E600"> i. Theory</p>
 +
<p class="content">Like the reverse direction of polypeptide formation, proteins can be hydrolyzed into their constituent amino acids. The mechanism undergoes an E2 elimination since a nucleotide attach to sp2 hybridized acyl carbon (amide). [2] Therefore, the reaction rate depends on the concentration of both nucleotides and the protein: </p>
 +
    <div>
 +
<img src=”” class=”picture”>
 +
<p class=”(1)”></p>
 +
</div>
 +
          <p class="content">Where, </p>
 +
          <p class="content">[<i>P</i>] is the concentation of Pantide, <i>K<sub>h</sub></i> is the reaction constant of hydrolysis, and <i>k<sub>A</sub></i>, <i>k<sub>N</sub></i>, <i>k<sub>B</sub></i> are the three contributing components of overall hydorlysis which are dependent to condition.</p>
 +
          <p class="content"> For Pantide, we would apply the protein solution in the farm at a constant pH. Thus, the rate of hydrolysis is proportional to the concentration of Pantide by a constant <i>K<sub>h</sub></i>, and the concentration of Pantide undergoes an exponential decay as time goes on.</p>
 +
  <p class="content-1" style="color:#00E600"> ii. Experimental proof</p>
 +
          <p class="content">We tested the chemical stability of both native and linear types of Hv1a and Hv1a-lectin in neutral PBS solvent (phosphate buffered saline, pH=7.4) in 4 &#8451; for one day and seven days. SDS-PAGE showed the remained protein concentration and calculated by software ImageJ. (Figure 2, Figure 3)</p>
 +
          <p class="content"></p>
 
           <p class="content"></p>
 
           <p class="content"></p>
 
         </div>
 
         </div>

Revision as of 23:29, 19 October 2016

Software—Toxin selection

I. Purpose

To prove the concept of Pantide, we wanted to select three existing distinct spider toxin peptides with probable oral toxicity against the testee－Spodoptera litura(Tobacco cutworms). For the actual application of Pantide, we needed some more knowledge base of peptides which have different molecular targets to promote Pantide applying to other orders of insects, and a different toxic mechanism to regularly alternate so as to avoid drug resistance.

To date, about 1500 toxin peptides from 97 spider species have been studied, though the number of spider toxin peptides is conservatively estimated up to 10 million. [1] So, our purpose is to establish a database collecting the information of those peptides, such as molecular target, taxon, toxicity, sequence. According to the database, if we first choose a target insect, then we can easily find out groups of suitable peptides used as Pantide. Therefore, we also need to create a method to select peptides from the database.

II. Method

The method of toxin selection can be separated into three part: crawler, filter, and selection.

  • Toxin Collection—we planned to collect information of toxin peptides to establish our own database for Pantide from protein databases and some research results like taxon and toxicity from published papers.
  • Toxin Filtering—based on background knowledge of toxin peptides, we set up some conditions to filter out those unsuitable to use as Pantide.
  • Toxin Processing—we used online protein analytic tools to classify the remained peptides into groups by their similarity. Finally, we select out three distinct peptides from different groups to proof concept of Pantide.

III. Step 1: Crawler

In the beginning, we searched on UniProtKB/Swiss-Prot. It is a freely accessible database of protein sequence and functional information that is the manually annotated and reviewed section. (http://www.uniprot.org/) By searching the keyword “insecticidal NOT crystal” we wanted to find all the proteins that have insecticidal activity excluding those crystal proteins of Bacillus thuringiensis, and we got 216 proteins as results.

Using the result, we established our Pantide database by crawling 11 entries of the protein information from UniProt. The entries are as follows.

  • The name of the protein
  • The description of protein function
  • The organisms/source of the protein sequence
  • The length of amino acids
  • The number of disulfides bonds
  • Propeptide & signal peptide—If the proteins have an N-terminal signal peptide and propeptide, a part of protein will be cleaved during maturation or activation.
  • Uniprot entry & Arachnoserver id—the accession number of protein in UniProtKB and ArachnoServer*.

*ArachnoServer is a manually curated database for protein toxins derived from spider venom.(http://www.arachnoserver.org/).

We also crawled other seven entries of protein toxicity recorded by Arachnoserver—molecular target, taxon, ED50, LD50, PD50, qualitative information, protein sequence from Arachnoserver. The term, Molecular target, is the effect site of toxin peptides, such as voltage-gated ion channels, GABA receptors and so on. Taxon, ED50, LD50, PD50, and the qualitative information are the toxicity against taxon that had been tested by experiments. The protein sequence from two databases is entirely the same.

We utilized BeautifulSoup 4.4.0, sqlite3 and gevent modules in Python 3.5 to develop our crawler. Moreover, we have submitted the code to GitHub.
(Link:https://github.com/chengchingwen/iGEM/blob/master/crawler.py)

IV. Step 2: Filter

After crawling the data, we used DB Browser for SQLite software to browse and used SQL to process our Pantide database. We tried to build a filter to find out peptides suitable to use as Pantide.

According to the previous articles, we knew that around 90% of spider venom toxin peptides contain ICK structure which is the most important domain that reacts with the voltage-gated ion channels of insects and some other receptors specifically. [2]

Therefore, to find these spider venom toxin peptides from Pantide database, we could start from searching for ICK structure, whose mass is among 1-10 kDa containing at least three disulfide bonds. [2] So we set a filter with three conditions.

  • The organism we choose must be spiders or tarantulas.
  • The length of the a.a. sequences are between 27 and 271 base pairs (1 kDa of protein has averagely nine amino acids, encoded by 27 base pairs)
  • The number of disulfide bonds is greater or equal to 3. After filtering with the three conditions, 113 peptides remained. Next, we set another filter to find out insecticidal peptides.
  • Molecular target contains “invertebrate,” but we also remain peptides without data.
          The reason why we keep the peptides without data was that they have the probability to be effective. In this stage, we got 63 candidates.

    • For efficacy experiment of Pantide, we choose our testee－Spodoptera litura as target insect. While there are 14 kinds of distinct Taxon in our database, including 4 Lepidoptera genus. Thus, we also set the other filter to find out peptides against Lepidoptera:

  • Taxon contains at least one of Spodoptera litura, Heliothis virescens, Manduca sexta and Spodoptera exigua, but we also remain peptides without data
          On the other hand, because we designed to produce Pantide by E.coli, that is difficult to express proteins containing disulfide bonds. We had chosen E.coli Rosetta-gami strain for enhanced disulfide bond formation, but to express a protein with more than four disulfide bonds is still a heavy load. So we finally filtered out those peptides containing too much disulfide bonds.
  • The number of disulfide bonds is less than or equal to four.
          The result was that we got 46 peptides which have the possibility to use as Pantide in proof concept experiment, and all of them is targeted to insects’ voltage-gated ion channels (excluding NULL).

V. Step 3: Selection

In this step, we tried to find three peptides that have different molecular target or mechanism from filtering result to do the test experiment. The method we used was to classify the remained peptides into groups by their structure similarity.

We used online analytic tools on NCBI to process those peptides.

We started with using Protein BLAST (Basic Local Alignment Search Tool) to search from the whole protein database for the similar query protein sequences related to all the 46 peptides and put those related peptides into groups.

The next was using COBALT (Constraint-based Multiple Alignment Tool) to align the sequence between groups to find out whether or not the two groups have the similar structure while they were not got together on the last step because of side chains and other factors. At last, we separated 46 peptides into four groups, containing 27, 12, 3, 2 peptides, and two alone.

Then we chose the three larger groups and used Conserved Domains Search, and found out that they belonged to the three conserved protein domain family. There are Omega-toxin Superfamily (cl05707), Toxin_28 Superfamily (cl06928) and Toxin_20 Superfamily (cl06915). The strings in brackets are unique ID of superfamilies in the conserved protein domain family database. Finally, we selected the representative peptides from each superfamily and got these three peptides, ω-hexatoxin-Hv1a, μ-segestritoxin-Sf1a and Orally active insecticidal peptide (OAIP).

VI. Future

To promote the applicability of Pantide, we still need to extend our database. The next step is to integrate with other toxin peptide databases, such as scorpions or cone snails, collect more peptides’ information from research results, and even combine with bioinformatics to build a new scoring system, and search for new potential peptides.

Reference

[1] King, G.F.; Gentz, M.C.; Escoubas, P.; Nicholson, G.M. A rational nomenclature for naming peptide toxins from spiders and other venomous animals. Toxicon 2008, 52, 264–276.

[2] Monique J. Windley, Volker Herzig, Sławomir A. Dziemborowicz, Margaret C. Hardy, Glenn F. King and Graham M. Nicholson (2012). Spider-Venom Peptides as Bioinsecticides. Toxins, 4, 191-227.

The Pest Prediction System

I. The Factors of Prediction

We built a prediction model for our device to support our spraying system. The model will predict the number of the pest, so we can know when our farm will be in the pest threat in the future. Thus we can open our spraying system spraying PANTIDE to protect our farmland. To know the number of pests in the future, we used seven weather data in the past 20 days including air pressure, the highest temperature, the lowest temperature, average temperature, humidity, precipitation, wind velocity, and also the accumulated number of bugs from the two periods including twenty to ten days ago and ten to one days ago. Because Pantide only influences on the larvae, we cannot easily predict the number of moths tomorrow or two or three days later. We needed to know the number of larvae in the future.

II. The Life History of Moth

After understanding the life history of the moths from Dr. Huang in Taiwan Agriculture Research Institute, we can know that the time of pupa becoming moth ranges from six days to 14 days. So the moths we caught in 20 days must be the larvae in our farm right now. Therefore, we utilized the data of the accumulated number of moth in the next 20 days as the target of our prediction system, and we can use the output to know the time that requires spraying Pantide for prevention.

III.The Software Design of Prediction Model

The method we used was called neural network, one of the popular machine learning solutions, or another well-known name called deep learning. The model we used was a combination of two kinds of neural network, Recurrent Neural Network (RNN) and Artificial Neural Network (ANN). Becuase we got seven features of 20 days, equals 140 features in total, we cannot simply put it into ANN to train the model. Therefore, we used RNN which was good for this kind of time series data to compress the 140 features into seven compressed data and got the seven feature and the other two feature which we have previously mentioned. Therefore, we only had nine features in total and then put them into the ANN as the input to get the final prediction. Then, it will compute the errors between the actual answer and the output and use gradient descent and backpropagation to modify the weight in each network each neuron. After a bunch of train steps, we got our model with about 80% of accuracy.

IV. The Quality of the Model

The prediction model is also a part of our device system. After using the device to collect data in a farm, we can not only predict the number of moth with the model we already trained but also retrain our model with the data we collect in that farm. So, now we can modify the model according to each farm to get a specific model for it.

The programming language we use is python3.5, and the deep learning framework is the one developed by Google called tendorflow.(See on GitHub)

Degradation Rate

I. Summary

The aim of this modeling was to predict and simulate the degradation rate of Pantide. Practically, the results would be integrated into our device to promote automatic control system. Once Pantide degraded below the effective level, it will spray the solution to replenish.

After discussion and experimental verification, we show that UV radiolytic oxidation contributes to the degradation of Pantide. We tried to build a model for simulating the relationship between applied UV light and degradation rate and then combined with UV intensity sensor to calculate the period of spraying Pantide.

II. Pantide degradation process

For the proteins used as Pantide, the inhibitor cystine knot (ICK) is significant to their function. Pantide with several disulfide bonds is often more stable in solution. If “native protein” is denatured, it loses disulfide bonds and becomes a less stable form, normally called “linear protein,” which is easily degradable.

Recent research shows that the spider toxin proteins containing ICK structure, for example, ω-hexatoxin-Hv1a (Hv1a), have high stability against temperature, pH, solvents and protease. In contrast, when Hv1a is denatured to linear form, it loses its stability and then degrades rapidly. [1]

There are many possible processes of Pantide degradation we discussed below. (Figure 1) Pantide may have a chance to be reduced to a linear form by reductants or reductases. For both native form and linear form proteins, it may suffer hydrolysis and proteolysis, resulting in denaturing or amino acid cleavage. Also, UV light of sun also leads to Pantide degradation. Though the energy of UV light may not be not enough to break the covalent bonds efficiently, proteins still could undergo radiolytic oxidation.

According to this degradation process, we attempted to build up our model. The degradation rate of Pantide is contributed by hydrolysis, proteolysis, UV radiolytic oxidation and the reduction to linear form. We considered that the rate Pantide transformed to linear form is a part of degradation rate because linear form protein loses its function and it is also likely to be degraded.

However, since the whole process is too complex to verify, we divided the experiment into three parts, hydrolysis test, proteolysis test and UV radiolytic oxidation test, and further summarized the results to conclude. The purpose is to find out the degradation rate of Pantide and to verify the less stable linear form protein has.

On the other hand, because ICK structure domain mainly contributes the stability of these proteins; we could assume that the degradation processes of three target proteins are roughly the same, for the reason that, we chose Hv1a and Hv1a with lectin (Hv1a-lectin) to demonstrate.

III. Hydrolysis test

i. Theory

Like the reverse direction of polypeptide formation, proteins can be hydrolyzed into their constituent amino acids. The mechanism undergoes an E2 elimination since a nucleotide attach to sp2 hybridized acyl carbon (amide). [2] Therefore, the reaction rate depends on the concentration of both nucleotides and the protein:

Where,

[P] is the concentation of Pantide, Kh is the reaction constant of hydrolysis, and kA, kN, kB are the three contributing components of overall hydorlysis which are dependent to condition.

For Pantide, we would apply the protein solution in the farm at a constant pH. Thus, the rate of hydrolysis is proportional to the concentration of Pantide by a constant Kh, and the concentration of Pantide undergoes an exponential decay as time goes on.

ii. Experimental proof

We tested the chemical stability of both native and linear types of Hv1a and Hv1a-lectin in neutral PBS solvent (phosphate buffered saline, pH=7.4) in 4 ℃ for one day and seven days. SDS-PAGE showed the remained protein concentration and calculated by software ImageJ. (Figure 2, Figure 3)

IV. Proteolysis test

V. UV radiolytic oxidation test

VI. Summarization

Reference

[1] Volker Herzig and Glenn F. King (2015). The Cystine Knot Is Responsible for the Exceptional Stability of the Insecticidal Spider Toxin ω-Hexatoxin-Hv1a. Toxins, 7, 4366-4380.

[2] Anonymous. (n.d.). HYDROLYSIS. Retrieved October 16, 2016 , from https://zh.scribd.com/document/79207692/Hydrolysis-2006

[3] Hedstrom, L. (2002, May 14). Serine Protease Mechanism and Specificity. Chem. Rev. 2002, 102, 4501-4523.

[4] Guozhong Xu & Mark R. Chance (2005). Radiolytic Modification of Sulfur-Containing Amino Acid Residues in Model Peptides: Fundamental Studies for Protein Footprinting. Anal. Chem, 77, 2437-2449.

[5] Mary E. Dzaugis, Arthur J. Spivack, Steven D'Hondt (2015, April 10). A quantitative model of water radiolysis and chemical production rates near radionuclide-containing solids. Radiation Physics and Chemistry, 115, 127-134.

[6] Bachari, T. S. (n.d.). Theoretical Investigation on The Kinetics of Free Radical Reactions of Styrene Emulsion Polymerization . Retrieved from http://www.iasj.net/iasj?func=fulltext&aId=13996