Team:XMU Software/Project/protein

<html>
<head>
	<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
	<script src="http://ajax.googleapis.com/ajax/libs/jquery/1.8.0/jquery.min.js"></script>

	<script type="text/javascript" async="" src="https://2013.igem.org/wiki/index.php?title=Team:XMU_Software/Project/js&action=raw&amp;ctype=text/javascript"></script>
	<link rel="stylesheet" type="text/css" href="https://2013.igem.org/wiki/index.php?title=Team:XMU_Software/css&action=raw&amp;ctype=text/css"/>
	<link rel="stylesheet" type="text/css" href="https://2013.igem.org/wiki/index.php?title=Team:XMU_Software/Project/css&action=raw&amp;ctype=text/css"/>
</head>
<body>
	 <div id="xmu_body">	
		<div id="xmu_software">
			<div id="xmu_software_right_matrix"></div>
			<!--project-->
			<div id="project"><a name="_project" id="_project"></a>
				<div class="navigation">
					<div>
						<ul id="jsddm">
							<li><a href="https://2013.igem.org/Team:XMU_Software">Home</a></li>
							<li><a href="https://2013.igem.org/Team:XMU_Software/Team">Team</a></li>
							<li><a href="https://2013.igem.org/Team:XMU_Software/Project/promoter">Project</a></li>
							<li><a href="https://2013.igem.org/Team:XMU_Software/Requirements">Requirements</a></li>
							<li><a href="https://2013.igem.org/Team:XMU_Software/Outreach">Outreach</a></li>
							<li><a href="https://2013.igem.org/Team:XMU_Software/Notebook">Notebook</a></li>			
						</ul>	
					</div>
					<div class="home_s3l1"></div>
				</div>
				<div id="project_s1">
					<canvas id="myCanvas" width="1347px" height="395px"></canvas>
					<div class="project_title">PROJECT</div>
					<div class="project_introduce">Our project includes 2 independent software tools-the brick worker and E' NOTE. The former is a software suit for the evaluation and optimization of biobricks, i.e., promoter, RBS, protein coding sequences and terminator. E' NOTE is a web application serving as an assistant for experiments. Its useful functions such as experiments recording and experimental template customization make experimental process easier and more enjoyable.</div>
				</div>
				<div class="navi_text">						
					<a id="project_navi0" class="navi0" href="https://2013.igem.org/Team:XMU_Software/Project/promoter"  onmouseover="navi0();" onmouseout="navi01();"   style="color:white">Promoter-decoder</a>
					<a id="project_navi1" class="navi1" href="https://2013.igem.org/Team:XMU_Software/Project/rbs"       onmouseover="navi1();" onmouseout="navi11();"   style="color:white">RBS-decoder</a>
					<a id="project_navi2" class="navi2" href="https://2013.igem.org/Team:XMU_Software/Project/terminator" onmouseover="navi2();" onmouseout="navi21();"   style="color:white">Terminator</a>
					<a id="project_navi3" class="navi3" href="https://2013.igem.org/Team:XMU_Software/Project/protein"    onmouseover="navi3();" onmouseout="navi31();"   style="color:white">SynoProteiner</a>
					<a id="project_navi4" class="navi4" href="#https://2013.igem.org/Team:XMU_Software/Project/enote"     onmouseover="navi4();" onmouseout="navi41();"   style="color:white">E' NOTE</a>
					<a id="project_navi5" class="navi5" href="#https://2013.igem.org/Team:XMU_Software/Project/achieve"   onmouseover="navi5();" onmouseout="navi51();"   style="color:white">Achievement</a>
					<a id="project_navi6" class="navi6" href="javascript:void(0)"                                 onmouseover="navi6();" onmouseout="navi61();"   style="color:white">Back To Top</a>
				</div>
					
				<div id="protein"><a name="_protein" id="_protein"></a>
						<div id="protein_m1">
							<img src="https://static.igem.org/mediawiki/2013/5/50/%E6%A0%87%E9%A2%98%E6%A0%8F.png" width="1348px" class="sep_photo">
							<div class="sep_title">SynoProteiner</div>
						</div>
						<div id="protein_m2">
							</br></br>
							<!--Abstract-->
							<h1>Abstract</h1>
							<p>Our team mainly focuses on programming the software by genetic algorithm, evaluating both optimization of single codon and codon pair and hence determining the fittest optimized sequences for expression in heterologous host cell.</p>
							<p>In addition, Synoproteiner protects enzyme cutting site. Users can choose corresponding cutting sites, not substituted during the optimization, so that the power of chosen enzyme cutting site won’t fail.</p>
							<p>Apart from the optimization, we have two additional functions. One is the statistics analysis, which provides the numbers and the proportion of the codon in the original and optimized sequences, making the optimization easier to understand. The other is the prediction of the protein folding rate. The purpose of the prediction is to seek the law of the folding rate in general, computing a relatively accurate folding rate value of the optimized sequences for the users.</p>
							</br>
							
							<!--Background-->
							<h1>Background</h1>
							<h4>Synonymous codons and the efficiency</h4>
							<p>Except methionine and tryptophan, all amino acids can be encoded by two to six synonymous codons, resulting from the degeneracy of the genetic code.<sup>1</sup>  However, unequal utilization of the synonymous condons leads to the phenomenon of codon usage bias, which is mainly due to natural selection, mutation and genetic drift.<sup>2</sup> According to related studies, codon usage bias has certain connection with gene expression level.<sup>3</sup> The larger the value of codon usage bias is, the higher gene expression will be. So the problem, how to substitute the synonymous codons aimed at raising the efficiency of gene expression and thus increasing the production of recombination protein in heterologous host cell, is expected to be addressed.</p>
							<h4>Protein folding rate</h4>
							<p>Protein is an important class of biological macromolecules. It is the main bearer of life activities and occupies a special position in vivo. Each protein has its own unique amino acid composition and sequences. Only when the amino acid chain is folded into the correct three-dimensional structure, will the protein have normal biological functions. Misfolded ones will not only lose its biological function but also even cause diseases such as mad cow disease, Alzheimer's syndrome, etc. The protein folding problem, an important biological question that the central dogma of molecular biology has not solved yet, has been listed as an important topic in twenty-first century. The folding mechanism of the protein is a challenging task, one of which is to determine factor influencing the folding rate. Although the answer can be found in a variety of biological experiments, such as various spectroscopy, mass spectrometry and nuclear magnetic resonance, these methods are time-consuming and costly. With the development of physics, mathematics, especially the progress of computer technology, how to apply a fast and accurate calculation method to predict protein folding rate attracts more and more attention.<sup>4</sup></p>
							</br>
							
							<!--Introduction-->
							<h1>Introduction</h1>
							<h4>Balance with single codon and codon pair</h4>
							<p>Individual codon usage optimization has been attached importance to, taking Codon optimizer,<sup>5</sup> Gene Designer,<sup>6</sup> OPTIMIZER<sup>7</sup> for example. Subsequently, people found the effect of gene expression optimization cannot be perfect just by single codon optimization. Codon pair, namely the pair of k-th and (k+1)-th codons from the 5’ to 3’ end, is another crucial factor. Due to potential tRNA-tRNA steric interaction within the ribosomes,<sup>8</sup> the usage of rare condon pairs, which correlate with translation elongation, decrease protein translation rates.<sup>9</sup> Optimization of individual codon has an influence on the corresponding codon pair resulting in maybe-not-the-best codon pair optimization. In the same way, optimizing codon pair merely contributes to maybe-not-the-best single codon optimization. Therefore, it is a challenging way for us to apply a method considering and weighing the effects of single codon and codon pair optimization and thus make the whole best.</p>
							<p>Our team focuses on evaluating both optimization of single codon and codon pair and thus selecting the best sequences for expression in heterologous host cell.</p>
							<h4>Host Cell</h4>
							<p>Considering <i>E. coli</i> and <i>S. cerevisiae</i> are the ideal hosts for recombinant proteinexpression, and Gram-positive bacterium <i>L. lactis</i> and methylotrophic yeast <i>P. pastoris</i> are also promising candidates for expressing recombinant proteins,<sup>10</sup> we attached importance to selecting these four kinds of bacterium as host cell to optimize the sequences.</p>
							<h4>Enzyme cutting site protection</h4>
							<p>To avoid deconstruction of enzyme cutting site and therefore incapability during substitution of the base, Synoproteiner succeeds in exchanging <a href="https://2013.igem.org/Team:XMU_Software/Project_protein/cutting_site" target="_blank">14 kinds of commonly used cutting sites </a>such as EcoRI, XbaI and so on.</p>
							


							<p>Users can multi-select the enzyme cutting site that they want to stay unaltered and hence get an optimized sequence without changing selected region.</p>
							<h4>Method of prediction</h4>
							<p>In recent years, many researchers have made great efforts to explore the determinants of the folding rate, and various forecasting methods have been proposed. The existed prediction methods can be roughly divided into three categories.<sup>11-12</sup> The first one is based on the tertiary structure.<sup>13-19</sup> However, it takes lots of molecular experiments, expensive and in long period, to acquire the information of the tertiary structure, which fails to meet the demand of rapid prediction. The second category is based on the secondary structure.<sup>20-24</sup> This kind of method requires information of the secondary structure, similarly obtained by molecular experiments, or from the primary sequences prediction, but it will be limited by accuracy of the secondary structure prediction method. The last one is based on the primary structure,<sup>25-34</sup> which predicts the folding rate from amino acid sequences without most structure information.<sup>4</sup> And our prediction of the protein folding rate focuses on the last method.</p>
							</br>
							
							<!--Algorithm-->
							<h1>Algorithm</h1>
							<h4>Part I—Optimization</h4>
							<h3>Basic Table</h3>
							<p>Based on the synonymous codons table, we calculate, we calculate function of single codon, function of codon pair  and the function of multi-objective codon optimization. The method aims at make the optimization of whole best by calculating the relative effect of single codon and codon pair.</p>
							
							
							
							<h4>Fitness calculation<sup>36</sup></h4>
							<h3>Fitness function：</h3>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/0/08/Xmusoftware_Protein1.png"></div>
							<p>In the function,</p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/b/b7/Xmusoftware_Protein2.png"></div>
							<p><i>cpi</i> is a value larger than zero, ranging from 10<sup>-4</sup> to 0.5,<i>fit<sub>cp</sub></i> (<i>g</i>) is the fitness function of the codon pair，<i>fit<sub>sc</sub></i> (<i>g</i>) is the fitness function of the single codon，<i>w</i> ( (<i>c</i> (<i>k</i>),<i>c</i> (<i>k+1</i>)) is the weight of codon pairs in sequences <i>g</i>,|<i>g</i>| is the length of encoding sequences, <i>c</i> (<i>k</i>) is k-th codon in the sequences,<img src="https://static.igem.org/mediawiki/2013/6/6c/Xmusoftware_Protein3.png"> is the target ratio of k-th codon, <img src="https://static.igem.org/mediawiki/2013/0/0a/Xmusoftware_Protein4.png"> is the actual ratio of k-th codon in the sequences，the best value of <i>cpi</i> is 0.2 in the software.</p>
							<p>In the function, the target ratio of k-th codon can be approximated by the equation below:</p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/6/66/Xmusoftware_Protein5.png"></div>
							<p>In the function, weight can be calculated by the equation below:</p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/e/e5/Xmusoftware_Protein6.png"></div>
							<p><img src="https://static.igem.org/mediawiki/2013/2/26/Xmusoftware_Protein8.png">stands for the ratio of single codon c<sub>k</sub> in the complete genome'<img src="https://static.igem.org/mediawiki/2013/f/f6/Xmusoftware_Protein9.png">is the number of pair ( <i>c<sub>i</sub>,c<sub>j</sub></i> )  in high-expression genes, and high-expression genes are genes whose copy numbers of mRNA can be detected at least 20 per cell. </p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/1/1f/Xmusoftware_Protein10.png"></div>
							<p>syn (c<sub>k</sub>) stands for the synonymous codon set related to <i>c<sub>k</sub></i>,<img src="https://static.igem.org/mediawiki/2013/a/ae/Xmusoftware_Protein11.png">equals to the number of amino acid encoded by <i>c<sub>i</sub></i> in the whole protein set. </p>
							
							
							<h3>Calculation of Fitness</h3>
							<p>NSGA-II algorithm applied<sup>35</sup>:</p>
							<p>1.	Randomly initialize a population of coding sequences for target protein.</p>
							<p>2.	Evaluate <i>fit<sub>sc</sub></i> (<i>g</i>) and <i>fit<sub>sp</sub></i> (<i>g</i>) fitness of each sequences in the population.</p>
							<p>3.	Group the sequences into nondominated sets and rank the sets.</p>
							<p>4.	Check termination criterion.</p>
							<p>5.	If termination criterion is not satisfied, select the “fittest” sequences (top 50% of the population) as the parents for creation of offsprings via recombination and mutation.</p>
							<p>6.	Combine the parents and offsprings to form a new population.</p>
							<p>7.	Repeat steps 2 to 5 until termination criterion is satisfied.</p>
							<p>The identification and ranking of nondominated sets in step 3 is performed via pair-wise comparison of the sequences' <i>fit<sub>sc</sub></i> (<i>g</i>) and <i>fit<sub>sp</sub></i> (<i>g</i>) fitness. For a given pair of sequences with fitness values expressed as (<i>fit<sup>1</sup><sub>cp</sub></i> (<i>g</i>),<i>fit<sup>1</sup><sub>sc</sub></i> (<i>g</i>)) and (<i>fit<sup>2</sup><sub>cp</sub></i> (<i>g</i>),<i>fit<sup>2</sup><sub>sc</sub></i> (<i>g</i>)), the domination status can be evaluated as follows:</p>
							<p>•	If (<i>fit<sup>1</sup><sub>sc</sub></i> (<i>g</i>)&gt;<i>fit<sup>2</sup><sub>sc</sub></i> (<i>g</i>))and (<i>fit<sup>1</sup><sub>cp</sub></i> (<i>g</i>)&gt;=<i>fit<sup>2</sup><sub>cp</sub></i> (<i>g</i>)), sequences 1 dominates sequences 2.</p>
							<p>•	If (<i>fit<sup>1</sup><sub>sc</sub></i> (<i>g</i>)&gt;=<i>fit<sup>2</sup><sub>sc</sub></i> (<i>g</i>)) and (<i>fit<sup>1</sup><sub>cp</sub></i> (<i>g</i>)&gt;<i>fit<sup>2</sup><sub>cp</sub></i> (<i>g</i>)), sequences 1 dominates sequences 2.</p>
							<p>•	If  (<i>fit<sup>1</sup><sub>sc</sub></i> (<i>g</i>)&lt;<i>fit<sup>2</sup><sub>sc</sub></i> (<i>g</i>))and (<i>fit<sup>1</sup><sub>cp</sub></i> (<i>g</i>)&lt;=<i>fit<sup>2</sup><sub>cp</sub></i> (<i>g</i>)), sequences 2 dominates sequences 1.</p>
							<p>•	If  (<i>fit<sup>1</sup><sub>sc</sub></i> (<i>g</i>)&lt;=<i>fit<sup>2</sup><sub>sc</sub></i> (<i>g</i>))and (<i>fit<sup>1</sup><sub>cp</sub></i> (<i>g</i>)&lt;<i>fit<sup>2</sup><sub>cp</sub></i> (<i>g</i>)), sequences 2 dominates sequences 1.</p>
							<p>The process is showed in the figure below:</p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/7/7b/Xmusoftware_Protein30.jpg"></div>
							<div class="body_photo">Figure 1</div>
							<h3>Result:</h3>
							<p>By this method, there are enough experimental data to prove the sequences optimized works. Xylose isomerase in <i>Bacillus stearothermophilus</i>, Xylose isomerase in <i>Streptomyces olivochromogenes</i> and L-arabinose isomerase in <i>Thermoanaerobacter mathranii</i>  all, the optimized ones, were highly expressed  in <i>Bacillus subtilis</i>. In addition, the activity of the optimized Aspergillusniger fungal amylase was enhanced to 400% compared with the original sequences in <i>A. niger</i>.<sup>36</sup></p>
							<p><a href="https://2013.igem.org/Team:XMU_Software/Project_protein/Example_of_optimization" target="_blank">Link to the page of Example of optimization</a></p>
							<h4>Part II—Prediction of protein folding rate</h4>
							<p>In order to illustrate protein folding rate quantitatively, we determine the folding rate of 60 kinds of proteins as an experimental data set from literature and database<sup>37</sup>,  and information of the sequences comes from PBD and NCBI.</p>
							<table border="1" cellspacing="0" class="protein_table2">
								<tr>
									<th>&nbsp;&nbsp;&nbsp;protein&nbsp;&nbsp;&nbsp;</th>
									<th>Logarithm of the folding rate Ln(k<sub>f</sub>)</th>
									<th>&nbsp;&nbsp;&nbsp;protein&nbsp;&nbsp;&nbsp;</th>
									<th>Logarithm of the folding rate  Ln(k<sub>f</sub>)</th>
									<th>&nbsp;&nbsp;&nbsp;protein&nbsp;&nbsp;&nbsp;</th>
									<th>Logarithm of the folding rate  Ln(k<sub>f</sub>)</th>
								</tr>
								<tr>
									<td>2PDD</td>
									<td>9.8</td>
									<td>1FKB</td>
									<td>1.5</td>
									<td>1RA9</td>
									<td>-2.5</td>
								</tr>
								<tr>
									<td>2ABD</td><td>6.6</td><td>2CI2</td><td>3.9</td><td>1QOP</td><td>-6.9</td>
								</tr>
								<tr>
									<td>256B</td><td>12.2</td><td>1URN</td><td>5.8</td><td>1PHP</td><td>2.3</td>
								</tr>
								<tr>
									<td>1IMQ</td><td>7.3</td><td>1APS</td><td>-1.5</td><td>1PHP</td><td>-3.5</td>
								</tr>
								<tr>
									<td>1LMB</td><td>8.5</td><td>1RIS</td><td>5.9</td><td>1BNI</td><td>2.6</td>
								</tr>
								<tr>
									<td>1WIT</td><td>0.4</td><td>1POH</td><td>2.7</td><td>2LZM</td><td>4.1</td>
								</tr>
								<tr>
									<td>1TEN</td><td>1.1</td><td>1DIV</td><td>6.1</td><td>1UBQ</td><td>5.9</td>
								</tr>
								<tr>
									<td>1SHG</td><td>1.4</td><td>2VIK</td><td>6.8</td><td>1SCE</td><td>4.2</td>
								</tr>
								<tr>
									<td>1SRL</td><td>4</td><td>1A6N</td><td>1.1</td><td>1YCC</td><td>9.62</td>
								</tr>
								<tr>
									<td>1PNJ</td><td>-1.1</td><td>1CEI</td><td>5.8</td><td>1VII</td><td>11.52</td>
								</tr>
								<tr>
									<td>1SHF</td><td>4.5</td><td>2CRO</td><td>3.7</td><td>1NYF</td><td>4.54</td>
								</tr>
								<tr>
									<td>1PSF</td><td>3.2</td><td>2A5E</td><td>3.5</td><td>2AIT</td><td>4.2</td>
								</tr>
								<tr>
									<td>1CSP</td><td>7</td><td>1IFC</td><td>3.4</td><td>1PIN</td><td>9.44</td>
								</tr>
								<tr>
									<td>1C9O</td><td>7.2</td><td>1EAL</td><td>1.3</td><td>1C8C</td><td>6.91</td>
								</tr>
								<tr>
									<td>1G6P</td><td>6.3</td><td>1OPA</td><td>1.4</td><td>1BRS</td><td>3.4</td>
								</tr>
								<tr>
									<td>1MJC</td><td>5.3</td><td>1CBI</td><td>-3.2</td><td>1UBQ</td><td>5.9</td>
								</tr>
								<tr>
									<td>1LOP</td><td>6.6</td><td>1QOP</td><td>-2.5</td><td>3CHY</td><td>1</td>
								</tr>
								<tr>
									<td>1C8C</td><td>7</td><td>1BRS</td><td>3.4</td><td>1BIN</td><td>2.6</td>
								</tr>
								<tr>
									<td>1HZ6</td><td>4.1</td><td>3CHY</td><td>1</td><td>1SCE</td><td>4.2</td>
								</tr>
								<tr>
									<td>1PGB</td><td>6</td><td>2RN2</td><td>0.1</td><td>1GXT</td><td>4.38</td>
								</tr>
							</table>
							<p>In order that the characteristic factors of the folding rate can be extracted from protein sequences, we introduced the Chou's pseudo amino acid composition concept.<sup>38</sup> According to the pseudo amino acid composition principle, the position information of protein sequences can be, to some extent, reflected by a group of serial correlation factors <i>θ<sub>1</sub>，θ<sub>2</sub> ，θ<sub>3</sub>……，θ<sub>n</sub></i> ,which is defined as follows:</p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/f/fb/Xmusoftware_Protein31.png"></div>
							<p>in the function, <i>θ<sub>1</sub></i> is called the first-tier correlation factor that reflects the sequences order correlation between all the most contiguous residues along a protein chain (Fig. 2a), <i>θ<sub>2</sub></i> the second-tier correlation factor that reflects the sequences order correlation between all the second most contiguous residues (Fig.2b), <i>θ<sub>3</sub></i> the third-tier correlation factor that reflects the sequences order correlation between all the 3rd most contiguous residues (Fig.2c), and so forth.<sup>38</sup></p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/f/f8/Xmusoftware_Protein32.jpg"></div>
							<div class="body_photo">Figure 2</div>
							<p>the correlation function is given by<sup>4</sup>:</p>
							<div class="body_photo">Θ(R<sub>i</sub>,R<sub>j</sub>)=|<i>H</i>(<i>R<sub>j</sub></i>)-<i>H</i>(<i>R<sub>i</sub></i>)|</div>
							<p>where <i>H1(R<sub>i</sub>)), H2(R<sub>i</sub>), and M(R<sub>i</sub>)</i> are, respectively, the hydrophobicity value. Studies have shown that λ=10 will be the best predictor.<sup>39</sup> But there will be a large amount of calculation considering all possible situations—the 30 factors. We should select factors that can obtain the best prediction accuracy in least calculation. For that reason, we drew lessons from the literature<sup>4</sup> by using the method of Monte Carlo simulation and then 14 optimal characteristic factor were obtained. Other studies have indicated that the logarithm of the sequences length has a good correlation with folding rate, so Ln (L) will be the fifteenth factors. We apply SPSS software to calculate the coefficient of 15 factor by multivariate linear regression, and this will be the forecast formula of the rate of protein folding. We compared the experimental data and the predicted data and the results are as follows:</p>
							<div class="body_photo"><img src="https://static.igem.org/mediawiki/2013/2/25/Xmu-software_synofig3.jpg"></div>
							<div class="body_photo">Figure 3 The comparision of the experimental and the predicted results</div>
							<p>Through the test, our software succeeded in showing a relatively accurate folding rate value.</p>
							</br>
							
							<!--Future work-->
							<h1>Future work</h1>
							<p>First of all, we will modify our software by advancing the program and the framework to improve its ability of concurrent computation and shorten the computing time.</p>
							<p>Secondly, to accelerate the calculation, we may simplify the function of calculation by neglecting some term in our equations. However, considering the time spent on running program was extremely little, we will pay more attention on how to modify the equations for increasing the accuracy which maybe dramatically progress optimization result.</p>
							<p>Thirdly, enriching the database is other way to improve our software. According to time-space tradeoff law, we could pre-process a bunch of sequences in common use to optimized one and save the result into our database. By assessing our data, investigators could select the optimized sequences for their synthesis. Then, users are required to feedback their result. When it collects enough information, our app will learn users’ bias therefore modify our optimizing function by some methods, like NSGA-II algorithm.</p>
							<p>The specific points are listed as following:</p>
							<p>1.	Shortening the computing time of the software.</p>
							<p>2.	Expanding the range of the host cells.</p>
							<p>3.	Improving bacterium's resistance to toxic molecule.</p>
							<p>4.	Advancing existed paths of synthetic biology by the method.</p>
							<p>5.	Designing new paths of synthetic biology by the method.</p>
							<p>6.	Increasing the output of recombinant protein.</p>
							<p>7.	Predicting the expression of heterologous gene in a new host cell.</p>
							<p>8.	Considering more factors such as spiral structure in folding which influence the folding rate and thereby obtaining more accurate prediction rate.</p>
							<p>9.	Providing a set of software tools for protein folding, especially in molecular dynamics simulation of protein folding.</p>
							</br>
							
							<!--References-->
							<h1>References</h1>
							<div class="body_references">[1] Grantham, R.; Gautier, C.; Gouy, M.; Mercier, R.; Pave, A., Codon catalog usage and the genome hypothesis. <i>Nucleic acids research  </i> <b>1980</b>, <i>8</i> (1),  197-197.</div>
							<div class="body_references">[2] Hershberg, R.; Petrov, D. A., Selection on codon bias. <i> Annual review of genetics </i> <b>2008</b>, <i> 42,  </i>287-299.</div>
							<div class="body_references">[3] Gouy, M.; Gautier, C., Codon usage in bacteria: correlation with gene expressivity. <i> Nucleic acids research </i> <b>1982</b>,  <i>10</i> (22), 7055-7074.</div>
							<div class="body_references">[4] 郭建秀，饶妮妮, 刘广雄, 李杰, &amp; 王云鹤. 从氨基酸序列预测蛋白质折叠速率. 生物化学与生物物理进展 <i>Progress in Biochemistry and Biophysics</i> <b>2010</b>, <i>37</i>(12): 1331~1338</div>
							<div class="body_references">[5] Fuglsang, A., Codon optimizer: a freeware tool for codon optimization.  <i>Protein expression and purification  </i><b>2003</b>,  <i>31</i> (2),  247-249.</div>
							<div class="body_references">[6] Villalobos, A.; Ness, J. E.; Gustafsson, C.; Minshull, J.; Govindarajan, S., Gene Designer: a synthetic biology tool for constructing artificial DNA segments. <i> Bmc Bioinformatics </i> <b>2006</b>,  <i>7</i> (1),  285.</div>
							<div class="body_references">[7]	Puigbò, P.; Guzmán, E.; Romeu, A.; Garcia-Vallvé, S., OPTIMIZER: a web server for optimizing the codon usage of DNA sequences.  <i>Nucleic acids research </i> <b>2007</b>, <i> 35</i> (suppl 2),  W126-W131.</div>
							<div class="body_references">[8] Smith, D.; Yarus, M., tRNA-tRNA interactions within cellular ribosomes.  <i>Proceedings of the National Academy of Sciences  </i><b>1989</b>,  <i>86 </i>  (12),4397-4401.</div>
							<div class="body_references">[9] Coleman, J. R.; Papamichail, D.; Skiena, S.; Futcher, B.; Wimmer, E.; Mueller, S., Virus attenuation by genome-scale changes in codon pair bias.  <i>Science </i> <b>2008</b>, <i> 320</i> (5884),  1784-1787.</div>
							<div class="body_references">[10] (a) Wildt, S.; Gerngross, T. U., The humanization of N-glycosylation pathways in yeast. Nature Reviews Microbiology <b>2005</b>, <i>3</i> (2), 119-128; (b) Morello, E.; Bermudez-Humaran, L.; Llull, D.; Sole, V.; Miraglio, N.; Langella, P.; Poquet, I., Lactococcus lactis, an efficient cell factory for recombinant protein production and secretion.  <i>Journal of molecular microbiology and biotechnology  </i><b>2007</b>,  <i>14</i> (1-3),  48-58.</div>
							<div class="body_references">[11] 郭建秀, 马彬广, 张红雨. 蛋白质折叠速率预测研究进展. 生物物理学报, <b>2006</b>, <i>22</i>(2):89-95 Guo J X, Ma B G, Zhang H Y. <i>Acta Biophys Sin</i>, <b>2006</b>, <i>22</i> (2):89-95.</div>
							<div class="body_references">[12] Gromiha M M, Selvaraj S. Bioinformatics approaches for understanding and predicting protein folding rates. <i>Current Bioinformatics</i>, <b>2008</b>, <i>3</i>(1): 1-9</div>
							<div class="body_references">[13] Plaxco K W, Simons K T, Baker D. Contact order, transition state placement and the refolding rates of single domain proteins. <i>J MolBiol</i>, <b>1998</b>, <i>277</i>(4): 985-994</div>
							<div class="body_references">[14] Gromiha M M, Selvaraj S. Comparison between long-range interactions and contact order in determining the folding rate of two-state proteins: application of long-range order to folding rate prediction. <i>J Mol Biol</i>, <b>2001</b>, <i>310</i>(1): 27-32</div>
							<div class="body_references">[15] Zhou H, Zhou Y. Folding rate prediction using total contact distance. <i>Biophys J</i>, <b>2002</b>, <i>82</i>(1): 458-463</div>				 
							<div class="body_references"> [16] Nölting B, Schälike W, Hampel P, <i>et al</i>. Structural determinants of the rate of protein folding.<i> J Theor Biol</i>, <b>2003</b>, <i>223</i>(3): 299-307</div>
							<div class="body_references"> [17] Weikl T R, Dill K A. Folding kinetics of two-state proteins: Effect of circularization, permutation, and crosslinks. <i>J Mol Biol</i>, <b>2003</b>,<i>332</i>(4): 953-963</div>
							<div class="body_references"> [18] Ivankov D N, Garbuzynskiy S O, Alm E, et al. Contact order revisited: influence of protein size on the folding rate. <i>Protein Sci</i>,<b>2003</b>, <i>12</i>(9): 2057-2062</div>
							<div class="body_references"> [19] Mirny L, Shakhnovich E. Protein folding theory: from lattice to all-atom models. <i>Annu Rev Biophys Biomol Struct</i>, <b>2001</b>, <i>30</i> (1):361-396</div>
							<div class="body_references"> [20] Gong H, Isom D G, Srinivasan R, et al. Local secondary structure content predicts folding rates for simple two-state proteins. <i>J MolBiol</i>, <b>2003</b>, <i>327</i>(5): 1149-1154</div>
							<div class="body_references"> [21] Ivankov D N, Finkelstein A V. Prediction of protein folding rates from the amino acid sequences-predicted secondary structure. <i>Proc Nat Acad Sci USA</i>, <b>2004</b>, <i>101</i>(24): 8942-8944</div>
							<div class="body_references"> [22] Fleming P J, Gong H P, Rose G D. Secondary structure determines protein topology. <i>Protein Sci</i>, <b>2006</b>, <i>15</i>(8): 1829-1834</div>
							<div class="body_references"> [23] Huang J T, Cheng J P, Chen H. Secondary structure length as a determinant of folding rate of proteins with two- and three-state kinetics. <i>Proteins</i>, <b>2007</b>, <i>67</i>(1): 12-17</div>
							<div class="body_references"> [24] Prabhu N P, Bhuyan A K. Prediction of folding rates of small proteins: empirical relations based on length, secondary structure content, residue type, and stability. <i>Biochemistry</i>,<b> 2006</b>, <i>45</i> (11):3805-3812</div>
							<div class="body_references"> [25] Shao H, Peng Y, Zeng Z H. A simple parameter relating sequenceswith folding rates of small helical proteins. <i>Protein Pept Lett</i>, <b>2003</b>,<i>10</i>(3): 277-280</div>
							<div class="body_references"> [26] Galzitskaya O V, Garbuzynskiy S O, Ivankov D N, et al. Chainlength is the main determinant of the folding rate for proteins withthree-state folding kinetics. <i>Proteins</i>, <b>2003</b>, <i>51</i>(2): 162-166</div>
							<div class="body_references"> [27] Huang J T, Jing T. Amino acid sequences predicts folding rate for middle-size two-state proteins. <i>Proteins</i>, <b>2006</b>, <i>63</i>(3): 551-554</div>
							<div class="body_references"> [28] Gromiha M M. A statistical model for predicting protein folding rates from amino acid sequences with structural class information.<i>J Chem Inf Model</i>, <b>2005</b>, <i>45</i>(2): 494-501</div>
							<div class="body_references"> [29] Ma B G, Guo J X, Zhang H Y. Direct correlation between proteins'folding rates and their amino acid compositions: an ab initio foldingrate prediction. <i>Proteins</i>, <b>2006</b>, <i>65</i>(2): 362-372</div>
							<div class="body_references"> [30] Gromiha M M, Thangakani A M, Selvaraj S. FOLD-RATE:prediction of protein folding rates from amino acid sequences.<i>Nucleic Acids Res</i>, <b>2006</b>, <i>34</i>(suppl_2): 70-74.</div>
							<div class="body_references"> [31] OuYang Z, Liang J. Predicting protein folding rates from geometric contact and amino acid sequences. <i>Protein Sci</i>, <b>2008</b>, <i>17</i>(7): 1256-1263</div>
							<div class="body_references"> [32] Huang L T, Gromiha M M. Analysis and prediction of proteinfolding rates using quadratic responde surface models. <i>J ComputChem</i>, <b>2008</b>, <i>29</i>(10): 1675-1683</div>
							<div class="body_references"> [33] Shen H B, Song J N, Chou K C. Prediction of protein folding ratesfrom primary sequences by fusing multiple sequential features.<i>J Biomedical Science and Engineering</i>, <b>2009</b>,<i> 2</i>(3): 136-143</div>
							<div class="body_references"> [34] Jiang Y, Iglinski P, Kurgan L. Prediction of protein folding ratesfrom primary sequences using hybrid sequences representation.<i>J Comput Chem</i>, <b>2009</b>, <i>30</i>(5): 772-783</div>
							<div class="body_references">[35] Chung, B.; Lee, D.-Y., Computational codon optimization of synthetic gene for protein expression.<i> BMC systems biology </i> <b>2012</b>,  <i>6</i> (1),  134.</div>
							<div class="body_references">[36] 	帝斯曼知识产权资产管理有限. 公司实现改进的多肽表达的方法: 中国, 200780024670.5[P]. 2009-07-22</div>
							<div class="body_references"> [37]	Gromiha, M. M.; Thangakani, A. M.; Selvaraj, S., FOLD-RATE: prediction of protein folding rates from amino acid sequences.<i> Nucleic acids research</i> <b>2006</b>, <i>34 </i>(suppl 2), W70-W74. </div>
							<div class="body_references"> [38] Chou, K. C., Prediction of protein cellular attributes using pseudo‐amino acid composition. <i>Proteins: Structure, Function, and Bioinformatics</i> <b>2001</b>, <i>43</i> (3), 246-255. </div>
							<div class="body_references"> [39]	Galzitskaya, O. V.; Garbuzynskiy, S. O.; Ivankov, D. N.; Finkelstein, A. V., Chain length is the main determinant of the folding rate for proteins with three‐state folding kinetics. <i>Proteins: Structure, Function, and Bioinformatics</i> <b>2003</b>, <i>51 </i>(2), 162-166. </div>
							</br></br></br></br></br></br></br></br>	
						</div>
					</div>
				
			</div>		
		</div>
	</div>
</body>
	<script>
		var c=document.getElementById("myCanvas");
		var cxt=c.getContext("2d");
		var grd=cxt.createLinearGradient(20,10,1348,395);
		grd.addColorStop(0,"#becbad");
		grd.addColorStop(0.5,"#e4d7aa");
		grd.addColorStop(1,"#becbad");
		cxt.fillStyle=grd;
		cxt.fillRect(0,0,1348,395);
	</script>
	<script>
		//定位1366*768
		
					var x = window.innerWidth-1366;
					if(x>0){
						x = window.innerWidth-1348;
						x = x / 2;
						document.getElementById("xmu_software").style.left= x+"px";
						$("#search-controls").get(0).style.right= x+"px";
						$(".right-menu").get(0).style.right= x+"px";
						$(".left-menu").get(0).style.left= x+"px";
						
						
					}
					else{
						x = window.innerWidth-1348;
						document.getElementById("xmu_software").style.left="0px";
						$("#search-controls").get(0).style.right=x+"px";
						$(".right-menu").get(0).style.right=x+"px";
						$(".left-menu").get(0).style.left="0px";
						
						
					}

					$("#xmu_software_right_matrix").mouseover(function(){$("#search-controls").show();$(".right-menu").show();$("#mw-searchButton").show();$("#searchGoButton").show();});
					 $("#search-controls").mouseover(function(){$("#search-controls").show();$(".right-menu").show();$("#mw-searchButton").show();$("#searchGoButton").show();});
					$(".right-menu").mouseover(function(){$("#search-controls").show();$(".right-menu").show();$("#mw-searchButton").show();$("#searchGoButton").show();});


					 $("#xmu_software_right_matrix").mouseout(function(){$("#search-controls").hide();$(".right-menu").hide();$("#mw-searchButton").hide();$("#searchGoButton").hide();});
				
					if($(".left-menu").find("li").size()==7)
					{
						$(".left-menu").find("li").eq(6).html("");  
						$("#menubar").get(0).style.width="330px";
					}
					else
					{
						$(".left-menu").find("li").eq(4).html("");  
						$("#menubar").get(0).style.width="280px";	
						$("#pt-login").get(0).style.right="50px";
					}
					
					
					$("#xmu_body").scroll(function(){ movetotop();$("#search-controls").hide();$(".right-menu").hide();$("#mw-searchButton").hide();$("#searchGoButton").hide();});
					
					color1=$("#notebook_s3p2").css("background-color");
					color2=$("#notebook_s3p3").css("background-color");

	</script>		

</html>