基因組規(guī)模代謝網絡模型自動化修正

吳曉紅， 薛 衛(wèi)， 張 梁*， 石貴陽

（1.糧食發(fā)酵工藝與技術國家工程實驗室，江南大學，江蘇 無錫 214122；2.南京農業(yè)大學 信息科學技術學院 江蘇 南京 210095）

基因組規(guī)模代謝網絡模型自動化修正

吳曉紅1， 薛 衛(wèi)2， 張 梁*1， 石貴陽1

（1.糧食發(fā)酵工藝與技術國家工程實驗室，江南大學，江蘇 無錫 214122；2.南京農業(yè)大學 信息科學技術學院 江蘇 南京 210095）

基于KEGG在線數(shù)據(jù)庫以及6個蛋白質區(qū)間預測數(shù)據(jù)庫，對基因組規(guī)模代謝網絡模型進行了自動化修正。作者提出了蛋白質區(qū)間預測結果的權重打分機制，同時利用圖像處理算法確定可信度高的特異性反應。上述修正的研究均在Spathaspora passalidarum NRRL Y-27907基因組規(guī)模代謝網絡精煉過程中得到運用實施，對于提高模型構建效率意義重大。

基因組規(guī)模；代謝網絡；斷點補齊；圖像處理；區(qū)間預測

隨著基因組高通量測序數(shù)據(jù)的涌現(xiàn)以及大量的生物學數(shù)據(jù)的產生，代謝網絡模型構建成為研究生物信息學的熱點之一。代謝網絡構建是一個耗時費力的過程，因此許多自動化構建的工具隨之應運而生。通常這些自動化工具側重關注代謝網絡粗模型的構建如 metaSHARK[1]和 AUTOGRAPH[2]，其次關注代謝網絡模型的模擬過程，如CellNetAnalyzer[3]、OptFlux[4]和 COBRA Toolbox[5]，只有少量的自動化工具是針對代謝網絡模型的精煉過程。目前能夠提供代謝網絡模型自動化精煉過程的工具有Model SEED、Pathway Tools、RAVEN 和 SuBliMinaL。

代謝網絡的模型構建包括粗模型的構建、模型的精煉、數(shù)學模型的轉換、模型的預測驗證四個過程。一個高質量的代謝網絡模型，應達到模型模擬結果和生物實際生長表型一致，否則要不斷的重復精煉修正過程，直到模擬與表型一致。模型的精煉修正無疑是代謝網絡模型構建過程中最耗時耗力的過程，現(xiàn)有模型精煉工具并不能真正實現(xiàn)真菌代謝網絡模型精煉過程的自動化。模型的精煉過程必須包括漏洞代謝的填補、反應區(qū)間定位等。Model SEED[6]和Pathway Tools[7]只能提供原核生物的代謝網絡模型的精煉自動化過程，不能提供反應區(qū)間的定位。RAVEN[8]和SuBliMinaL[9]是基于Wolf PSORT蛋白質區(qū)間預測數(shù)據(jù)庫實現(xiàn)自動化定位區(qū)間的程序。但是Wolf PSORT[10]只是基于氨基酸組成特征的在線預測數(shù)據(jù)庫。研究表明，基于氨基酸組成、二肽和物理化學三種綜合特征的蛋白質區(qū)間定位預測結果更為準確[11]。

利用作者所在實驗室自動化構建全基因組代謝網絡模型的程序，自動構建了Spathasporapassalidarum NRRL Y-27907全基因組規(guī)模代謝的粗模型。以S.passalidarum NRRL Y-27907的基因組規(guī)模代謝網絡模型的精煉過程為例，以簡單、面向對象的Java語言為基礎，對精煉過程中人工冗雜的斷點補齊的方法進行了研究，提出了一種基于KEGG[12]在線數(shù)據(jù)庫自動化填補漏洞反應的方法，并利用權重打分機制分析，6個真菌蛋白質定位數(shù)據(jù)庫預測S.passalidarum NRRL Y-27907的結果，在保證模型中反應的物種特異性的同時，實現(xiàn)了真菌代謝網絡模型精煉的自動化。自動化修正的流程見圖1。圖中進程g、進程n、進程o為一個小的流程循環(huán)。進程g中判斷反應包含斷點，則進入進程h，查找該反應在注釋圖譜中對應的坐標，并在進程i中讀取此坐標，在進程j中判斷此坐標是否為特異性反應，如果是，則在進程p中記錄該反應。如果不是，則在進程l中判斷此坐標是否為最后一個坐標，如果是最后一個坐標，則進入進程n，即進入進程g、進程n、進程o該流程循環(huán)。如果不是最后一個坐標，則進入進程m，讀取下一個坐標，判斷此坐標是否為特異性反應，重復此循環(huán)直至將所有的特異性反應都被找出，進入進程q，進行模型修正。在進程r中判斷模型中是否已經包含此反應，若已經包含，則回到進程n，即進入進程g、進程n、進程o該流程循環(huán)，檢查下一條反應。若不包含此反應，則進入進程s，將此反應加入到模型中。

圖1 自動補齊斷點流程Fig.1 Process of the auto-refinement of gap

1 自動填補網絡漏洞

采用柴文平[13]等人的方法構建了S.passalidarum NRRL Y-27907代謝網絡粗模型。構建的代謝網絡粗模型需要進一步精細化與修正，最終完成一個高質量的基因組規(guī)模代謝網絡模型。

1.1 代謝網絡漏洞查找

模型導入到裝有COBRA工具包和GLPK線性規(guī)劃器的Matlab中，將模型轉化為計算機可讀的格式 （SBML）才能進行代謝網絡漏洞查找。通過xls2model程序將模型Excel表讀取為計量學S矩陣。S矩陣（828×984）表示該模型由828個代謝物和984個反應組成。同時通過GapFind程序完成代謝漏洞的查找，其中上游漏洞代謝物有為44個，下游漏洞代謝物有128個。

1.2 基于KEGG網絡爬蟲反應

KEGG是代謝網絡構建常用數(shù)據(jù)庫，含有多個在線子數(shù)據(jù)庫，其中REACTION數(shù)據(jù)庫包含迄今為止發(fā)現(xiàn)的所有生化反應。各個子數(shù)據(jù)庫的網頁數(shù)據(jù)格式比較統(tǒng)一明確，方便人們進行遠程服務器訪問。但是，KEGG數(shù)據(jù)庫更新頻繁，各個子數(shù)據(jù)庫不能夠免費下載，需要付費使用。而在基因組代謝網絡斷點補齊過程中，因為數(shù)據(jù)信息量浩大，頻繁訪問遠程服務器比較耗時耗力。因此，實現(xiàn)一種批量在線獲取并存取數(shù)據(jù)的方法意義重大。

1.2.1 方法概述 利用超文本轉移協(xié)議和Java控件HttpClient相結合，實現(xiàn)對網頁中特定信息的抓取KEGG提供物種特異性基因組信息以及所有反應式信息查詢網頁，通過一定的URL（Uniform Resource Locator，統(tǒng)一資源定位符）格式地址發(fā)送HTTP請求并獲取網頁中的基因信息。在漏洞填補的過程中需要訪問大量不同的網絡資源，獲取相關的基因信息，由于數(shù)據(jù)量較大且人工操作比較繁瑣，這里利用Java控件HttpClient實現(xiàn)爬蟲技術，抓去符合特定條件的網絡資源。HttpClient是Apache Jakarta Common下的子項目，可以用來提供高效的、最新的、功能豐富的支持HTTP協(xié)議的客戶端編程工具包，并且它支持HTTP協(xié)議最新的版本和建議。利用HttpClient訪問具體的URL地址，獲取服務器端返回的獲取html內容，html內容由標題、js代碼、正文、相關鏈接、聲明等區(qū)域組成，而有用信息只出現(xiàn)在正文中的各種html標簽標記內，分析html標簽并獲取特定的網頁信息。

1.2.2 漏洞填補算法實現(xiàn)

1）獲取注釋圖譜：提交物種基因組蛋白質序列至KAAS自動注釋服務器，獲取注釋信息，下載html和text格式。

2）查找包含斷點的注釋圖譜：根據(jù)Matlab軟件中GapFind程序返回的漏洞代謝物列表，在代謝網絡模型Excel格式中確定代謝物的反應途徑，依據(jù)KASS注釋返回的途徑圖譜找到包含漏洞代謝物的所有反應。

注釋返回的KEGG代謝途徑為包含糖代謝等在內的110個途徑。查找包含斷點的代謝圖譜的流程見圖2。具體思路和偽代碼步驟如下：

A：獲取斷點化合物所對應的Subsystem信息，記為sub。

B：向注釋查詢網頁URL地址發(fā)送HTTP請求。

C：如果服務器端響應代碼為HTTPStatus.SC_OK則正常響應，否則繼續(xù)請求，獲取html正文內容。

D：分析html內容，設i為行號，由第一行開始遍歷標簽對中的每一行，

For i from 1 to n

｛

if（該行中第二個標簽中的內容與sub相等）

｛

提取對應的第一個標簽中的內容，記為KO；

｝

else

忽略該行，遍歷下一行；

｝

E：根據(jù)D中的KO號得到滿足條件圖譜的URL地址，向URL地址發(fā)送HTTP請求得到服務器端響應的網頁圖片記為T1，T1即為整個網絡結構圖，其中綠色酶號表示包含斷點的特異性反應。

F：點擊T1左上角途徑方框，進去包含所有反應頁面page1，網頁中每一個EC號對應圖譜中的一個具體反應，它的URL地址指向具體的反應方程式。

G：獲取page1中所有EC號對應的反應，設ec_num為每一個EC號，從第一個開始

EC_K_Break.txt保存包含斷點化合物的EC，K號的信息。3）查找EC_K_Break.txt中每個K對應的坐標根據(jù)K號獲取其在T1中對應的坐標，判斷特異性反應。

圖2 斷點代謝途徑定位Fig.2 Orientation of gap metabolic pathway

1.3 判斷特異性反應

KEGG所有的反應都包含在通路數(shù)據(jù)庫（PATHWAY database）中，PATHWAY 圖譜上有顏色標記的酶號是指這個物種特定的基因或酶，只有有顏色標記的酶號表示的反應才是具有該物種特異性的反應，也才能添加到代謝網絡模型中。在代謝網絡模型中添加非特異性的反應會改變整個代謝途徑和代謝物流量，進而使模型模擬的結果偏離實驗數(shù)據(jù)，影響模型的準確性和可信度。

構建代謝網絡模型需要提取代謝途徑中的特異性反應，圖中特異性反應對應的酶號所在的方形框有顏色標記。因此通過網絡爬蟲技術獲得方形框的位置列表，定位到某酶號所在的方形框后需要選取框內的像素點，讀取其顏色值，如果顏色分量RGB均為0或255，則沒有顏色標記，反之則有。代謝網絡特異性反應獲取流程見圖3。

基本思路為：

根據(jù)得到的position坐標讀取T1對應點的RGB色彩值。

Picture（Key：酶號；Value：代謝網絡圖中所有方形框的坐標向量集｛V1，V2，……，Vn｝）

For i from 1 to n

{

If（某酶號所在的方形框）

｛

沿方形框的長邊內側逐一選取像素點，讀取其顏色值；

If顏色分量RGB均為0或255 then沒有顏色標記

else有顏色標記；

If有顏色標記then該酶號對應的是特異性反應

do將反應加入菌的代謝網絡模型中；

else舍棄該酶號對應的反應。

｝｝

反應式漏洞填補

遍歷new_rec.TXT中每一個反應，查看模型中是否存在，存在則不處理，否則添加。

A：讀取new_rec.TXT中每行反應記為new_rec，i為行號

For i from 1 to n

｛

if（模型中不包含 new_rec）

｛

將new_rec添加到模型中；

｝

else

忽略該反應，查找下一條反應；｝

圖3 特異性反應獲取流程Fig.3 Process of getting the pecificreaction

2 獲取反應區(qū)間定位

細胞是生命活動的基本單位，它由執(zhí)行不同機體功能的稱為亞細胞的各部分組成，如細胞膜、細胞核、線粒體、高爾基體、內質網等。亞細胞功能是由位于其中的蛋白質執(zhí)行的，蛋白質所在的亞細胞稱為蛋白質的亞細胞位置[14]。蛋白質必須轉運到其應在的亞細胞位置上才能正確行使其功能，否則就會出現(xiàn)機體功能紊亂，正確合理的蛋白區(qū)間定位是高質量模型構建的基礎，見表1。

表1 真菌蛋白質亞細胞預測數(shù)據(jù)庫Table 1 Database for subcellular localization of fungal proteins

確定一條蛋白質的亞細胞位置稱為蛋白質亞細胞定位[15]。蛋白質亞細胞定位的傳統(tǒng)方法是通過生物化學實驗，如射線晶體衍射電子顯微鏡核磁共振等方法進行測定[16]。實驗方法精確度高，但費時耗力代價昂貴，而且對難于結晶的蛋白質來說，實驗方法不再有效。借助于先進高效的計算機自動化數(shù)據(jù)處理技術，出現(xiàn)了一些蛋白質定位預測網站。結合Spathasporapassalidarum NRRL Y-27907的生理生化性質和蛋白質特征提取方法、算法和準確性等，選取了6個真菌生物蛋白質區(qū)間預測網站，自動化提取分析網站的預測結果，在權重打分機制的基礎上得到最佳的蛋白質定位區(qū)間。這6個網站是基于蛋白質的氨基酸組成、偽氨基酸組成、二肽、生物化學特征或是四種特征的綜合。

2.1 區(qū)間定位算法實現(xiàn)

A：對每條反應獲取對應的KO號。

B：將A中的KO號在KASS注釋結果中查找基因號，并在本地下載Spathasporapassalidarum NRRL Y-27907蛋白質序列庫提取其對應的蛋白質序列。

C：將蛋白質序列提交到對應網站的表單中，獲取返回的定位信息。

D：獲取定位區(qū)間的信息并填入反應式中。

在獲取具體反應的區(qū)間信息過程中，需要將反應所對應的蛋白質序列提交到網頁的表單中，提交后返回具體的區(qū)間定位信息，此處會遇到兩個問題：1）表單提交過程中不支持大量蛋白質序列自動提交。由于模型中蛋白質序列數(shù)量較大，在有的網站中獲取定位信息時不支持大量序列的一次性提交而只能分別提交單個序列獲得定位信息，在提交過程中任務量大且人工耗費時間長。2）大量蛋白質序列提交耗費時間長，在網站中提交多個序列后等待服務器端反饋的定位信息耗費時間太長，甚至會發(fā)生無響應等問題，見圖4。

圖4 反應亞細胞定位流程Fig.4 Process of subcellular localization

HttpClient支持訪問特定的URL地址，獲取服務器端返回的html信息，并且能夠分析html中form表單中的信息，實現(xiàn)內容的自動提交。由于涉及到的定位頁面所有的表單提交方式都是POST提交，利用HttpClient中的PostMethod方法實現(xiàn)post提交。表單中的元素賦值過程：獲取表單中需要賦值的元素標簽，以蛋白質序列元素賦值標簽為例，標簽為，標簽名稱為 input，利用PostMethod中的addParameter方法為具體標簽元素賦值，形式如下PostMethod.addParameter（“input”，protein）；其中 protein 為蛋白質序列。 提交過程：調用Httpclient的executeMethod方法將賦值后的表單提交到服務器并獲取服務器返回的html信息，其中具體的定位信息包含在html文本中，分析提取即可。</p><h3>2.2 權重打分機制分析蛋白質區(qū)間定位</h3><p>6個蛋白質區(qū)間預測網站的訓練數(shù)據(jù)集不同，不能單純依據(jù)文獻里公布的預測結果準確率來預測Spathasporapassalidarum NRRL Y-27907各個蛋白質區(qū)間的的準確率。同時這6個網站各個區(qū)間預測的準確性也有區(qū)別。因此我們整合了RH2427[24]和PK7579[25]數(shù)據(jù)集，組成一個包含12個蛋白質區(qū)間，每個區(qū)間包含100條蛋白質序列的真菌蛋白質數(shù)據(jù)集，計算出6個網站基于新的真菌數(shù)據(jù)集上的權重。</p><p>權重是一個相對的概念，針對某一指標而言。某一指標的權重是指該指標在整體評價中的相對重要程度。研究6個蛋白質預測網站指標體系權重計算，反映各個蛋白質預測網站在預測結果中的重要性程度的數(shù)量。具體計算步驟如下：</p><p>第一步，統(tǒng)計每個預測網站各個區(qū)間預測正確的蛋白質序列個數(shù)。</p><p>第二步，計算出每個預測網站的平均識別正確數(shù)量，設X{X1，X2…X12}為每個預測網站12個區(qū)間的正確預測區(qū)間個數(shù)，則每個網站平均識別正確數(shù)量為：D=（X1+X2+…X12）/12， 計算結果見表 2 最后第二列。</p><p>第三步，計算6個預測網站的權重，計算結果見表 2最后一列。cello權重=89.3/（89.3+62.4+72.7+62.0+56.4+85.1） =0.208；PsortⅡ 權 重 =62.4/427.9=0.146；Epiloc 權重=72.7/427.9=0.170；Bacello權重=62.0/427.9=0.145；SLPFA權重=56.4/427.9=0.132；Euloc權重=85.1/427.9=0.199。</p><p>根據(jù)各個蛋白質區(qū)間預測網站的權重，采用加權投票的方式計算每個蛋白質序列最佳的蛋白質區(qū)間定位。本研究中采用的權值是每一個蛋白質區(qū)間預測網站的準確率，而不是對預測結果中區(qū)間票數(shù)的簡單加和，這樣就充分考慮到每一個蛋白質區(qū)間預測網站的準確率，有區(qū)別的對待了每一個預測網站，更符合實際預測結果。加權投票方式可以采用以下公式[26]表示：</p><p><img src="https://cimg.fx361.com/images/2023/0407/4370326f2cd1b5eeb8de4e2a5b9d20aa9677776e.webp"/></p><p><img src="https://cimg.fx361.com/images/2023/0407/03d6c8bb40dfea18576b40dbcfce2794cad0f15f.webp"/></p><p>表2 數(shù)據(jù)庫的權重Table 2 Weight of database</p><p>其中，Vi表示第i條蛋白序列的判決區(qū)間結果；wn為第n個蛋白區(qū)間預測網站的權重，其中保持表示第i條蛋白質序列在第n個區(qū)間預測網站上的預測結果，其為m，c，n，e等12個區(qū)間中的一個值（實際預測結果中待預測蛋白質可能只對應這12個區(qū)間值中的某幾個），N表示選取的真菌蛋白質區(qū)間預測網站個數(shù)，c表示所要預測的蛋白質序列數(shù)。當對輸入的待測蛋白質序列做判決時，把預測蛋白質區(qū)間在每一類區(qū)間的得票量排序，把待測蛋白質序列劃分到得票量最大的區(qū)間所在的類。當加權計算后，預測結果中對應多個得票量最大的區(qū)間時，則認定該蛋白質對應多個蛋白質預測區(qū)間。</p><p>下面通過具體的實例說明一下采用加權方式計算每個蛋白質序列最佳蛋白質區(qū)間定位的步驟。如SPAPADRAFT_69954編碼的蛋白質序列上傳到6個預測網站，返回的預測結果是cello、Epiloc、SLPFA預測的區(qū)間為cytos，PsortⅡ預測的區(qū)間為er，Bacello 預 測 區(qū) 間 為 mito，Euloc 預 測 結 果 為golgi。</p><p>第一步那么根據(jù)各個網站的權重和公式計算得</p><p>y=0.208cytos+0.146er+0.170cytos+0.145mito+0.132cytos+0.199golgi=0.51cytos+0.146er+0.145mito+0.199golgi</p><p>第二步，預測結果在每一類區(qū)間得票量排序，0.51>0.199>0.146>0.145，蛋白質序列在各個區(qū)間的概率為cytos>golgi>er>mito，則最佳蛋白質預測區(qū)間為cytos。</p><h2>3 討 論</h2><p>根據(jù)TCDB和TransportDB數(shù)據(jù)庫返回的轉運反應信息，添加不同細胞器間的代謝物轉運反應，使整個代謝物網絡模型連接起來。將精煉后的模型轉化為計算機可讀的格式（SBML）進行模擬分析。通過xls2model程序將模型Excel表讀取為計量學S矩陣。模型由包含828個代謝物和984個反應的粗模型，擴充到包含873個代謝物和1 243個反應的精細模型。</p><p>在補充斷點的過程中，對于模型新添加的代謝物信息，自動化從網站上抓取，節(jié)省了大量的時間和勞力。代謝物信息則包括代謝物縮寫（mets）、代謝物全稱（metNames）、代謝物分子式（metFormulas）、帶 電 荷 數(shù) （metCharge）、 代 謝 物 分 區(qū)（metCompartment）、 在 KEGG數(shù)據(jù)庫中的編號（metKEGGID）、在 PubChem數(shù)據(jù)庫中的編號（metPubChemID）、 在 EBI數(shù) 據(jù) 庫 中 的 編 號（metChEBIID）等。</p><p>基于Matlab平臺的COBRA工具，打開GLPK線 性 規(guī) 劃 器 ， 輸 入 [allGaps，rootGaps，downstreamGaps]=GapFind （model，false，false） 命令，補齊后的模型，不存在任何斷點，說明此普適性自動化程序的可行性。</p><h2>4 結語</h2><p>國內外雖然對代謝網絡自動化修正也做了多方面的研究，但仍然不能夠實現(xiàn)完全自動化修正真菌代謝網絡，代謝網絡模型構建過程中仍然需要手工添加模型特異性反應和反應區(qū)間。作者提出基于KEGG數(shù)據(jù)庫的網絡爬蟲和自動化提取蛋白質區(qū)間預測結果，并加權計算分析結果，能夠在修正一個相對精煉模型的過程中實現(xiàn)計算機技術與代謝網絡模型修正的最大化結合，減少了修正過程中大量的勞力與時間，提高了代謝網絡構建的效率與精確性。</p><p>[1]PINNEY J W，SHIRLEY M W，MCCONKEY G A，et al.MetaSHARK：software for automated metabolic network prediction from DNA sequence and itsapplication to the genomesof Plasmodium falciparum and Eimeriatenella[J].Nucleic Acids Research，2005，33（4）：1399-1409.</p><p>[2]NOTEBAART R A，van Enckevort F H，F(xiàn)RANCKE C，et al.Accelerating the reconstruction of genome-scale metabolic networks[J].BMC Bioinformatics，2006，13（7）：296.</p><p>[3]STEFFEN K，JULIO S R，ERNST D G.Structural and functional analysis of cellular networks with CellNetAnalyzer[J].BMC Systems Biology，2007，1（1）：1-13.</p><p>[4]ROCHA I，MAIA P，EVANGELISTA P，et al.OptFlux：an open-source software platform for in silico metabolic engineering[J].BMC Systems Biology，2010，l 4（1）：45-57.</p><p>[5]SCHELLENBERGER J，QUE R，F(xiàn)LEMING R M T，et al.Quantitative prediction of cellular metabolism with constraint-based models：the COBRA Toolbox v2.0[J].Nature Protocols，2011，6（9）：1290-1307.</p><p>[6]HENRY C S，DEJONGH M，BEST A A，et al.High-throughput generation，optimization and analysis of genome-scale metabolic models[J].Nature Biotechnology，2010，28（9）：977-982.</p><p>[7]KARP P D，PALEY S M，ROMERO P.The pathway tools software[J].BMC Bioinformatics，2002，18：S225-S232.</p><p>[8]AGREN R，LIU L M，SHOAIE S，et al.The RAVEN Toolbox and its use for generating a genome-scale metabolic model for Penicillium chrysogenum[J].PLoS Computational Biology，2013，9（3）：1-16.</p><p>[9]SWAINSTON N，SMALLBONE K，MENDES P，et al.The suBliMinaL toolbox：automating steps in the reconstruction of metabolic networks[J].J IntegrBioinform，2011，8（2）：186.</p><p>[10]PAUL H，KEUN-JOON P，TAKESHI O，et al.WoLF PSORT：protein localization predictor[J].Nucleic Acids Research，2007，35：W585-W587.</p><p>[11]LU Z，SZAFRON D，GREINER R，et al.Predicting subcellular localization of proteins using machine-learned classifiers[J].BMC Bioinformatics，2004，20（4）：547-556.</p><p>[12]KANEHISA M，GOTO S.KEGG：Kyoto encyclopedia of genes and genomes[J].Nucleic Acids Research，2002，28（1）：27-30.</p><p>[13]CHAI Wenping，XUE Wei，ZHANG Liang，et al.Research on the auto-reconstruction of genome-scale metabolic network model[J].Journal of Food Science and Biotechnology，2014，33（9）：957-967.（in Chinese）</p><p>[14]韓榕.細胞生物學[M].北京：科學出版社，2011：55-106.</p><p>[15]葉子弘.生物信息學[M].杭州：浙江大學出版社，2011：179-223.</p><p>[16]LIU Liyuan，CHEN Yuehui，MA Bingxian，et al.Prediction of protein subnuclear location using evolutionary fuzzy K-nearest neighbors and its ensemble[J].Journal of University of JINAN，2010，24（4）：376-379.（in Chinese）</p><p>[17]NAKAI K，HORTON P.PSORT：a program for detecting sorting signals in proteins and predicting their subcellular localization[J].Trends in Biochemical Sciences，1999，24（1）：34-36.</p><p>[18]ANDREA P，PIER L M，F(xiàn)ARISELLI P，et al.BaCelLo：a balanced subcellular localization predictor[J].BMC Bioinformatics，2006，22（14）：408-416.</p><p>[19]YU C S，LIN C J，HWANG J K.Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions[J].Protein Science，2004，13（5）：1402-1406.</p><p>[20]YU C S，CHEN Y C，LU C H.Prediction of protein subcellular localization[J].Proteins：Structure，F(xiàn)unction and Genetics，2006，64（3）：643-651.</p><p>[21]BRADY S，SHATKAY H.EpiLoc：a （working） text-based system forpredictingprotein subcelluarlocation[J].Pacific Symposium on Biocomputing，2008（13）：604-615.</p><p>[22]TAMURA T，AKUTSU T.Subcellular location prediction of proteins using support vector machines with alignment of block sequences utilizing amino acid composition[J].BMC Bioinformatics，2007，8（1）：466-478.</p><p>[23]CHANG TH1，WU L C，LEE T Y.EuLoc：a web-server for accurately predict protein subcellular localization in eukaryotes by incorporating various features of sequence segments into the general form of Chou's PseAAC[J].Journal of Computer-Aided Molecular Design，2013，27（1）：91-103.</p><p>[24]GARG A，RAGHAVA P S.ESLpred2：Improved method for predicting subcellular localization of eukaryotic proteins[J].BMC Bioinformatics，2008，9（1）：1-10.</p><p>[25]CHOU K C，SHEN H B.Euk-mPLoc：A fusion classifier for large-scale eukaryotic protein subcellular location prediction by incorporating multiple sites[J].Journal of Proteome Research，2007，6（5）：1728-1734.</p><p>Auto-Refinement of Genome-Scale Metabolic Network Model</p><p>WU Xiaohong1， XUE Wei2， ZHANG Liang*1， SHI Guiyang1<br/>(1.National Engineering Laboratory for Cereal Fermention Technology，Jiangnan University，Wuxi 214122，China；2.College of Information Science and Technology，Nanjing Agricultural University，Nanjing210095，China）</p><p>KEGGonline database and six subcellular prediction databaseshave been studied for the process of auto-refinement.The weighted scoring mechanism was proposed to analyze the results of subcellular prediction databases，using image processing algorithm to determine high credibility specific reaction.As an illustration example，all of the automatic methods were implemented in the process of genome-scale metabolic network refinement of Spathasporapassalidarum NRRLY-27907，whichconfirmed that these methods can improve the efficiency of model reconstruction.</p><p>genome-scale，metabolicnetworks，gapsupplement，imageprocessing，subcellularprediction</p><p>TP 391；Q 939</p><p>A</p><p>1673—1689（2017）09—0982—08</p><p>2015-03-02</p><p>江蘇省自然科學基金項目（BK2012363，BK2011153）。</p><p>*通信作者：張 梁（1978—），男，江蘇無錫人，工學博士，教授，博士研究生導師，主要從事代謝工程方面的研究。</p><p>E-mail：zhangl@jiangnan.edu.cn</p><p>吳曉紅，薛衛(wèi)，張梁，等.基因組規(guī)模代謝網絡模型自動化修正[J].食品與生物技術學報，2017，36（09）：982-989.</p></p> <!-- <div id="g0gggggg" class="article_pdf"><a >查看pdf文檔請下載app</a></div>--><div id="g0gggggg" class="article_love"> <div id="g0gggggg" class="title">猜你喜歡</div> <div id="g0gggggg" class="article_love_keyword"><span><a href="/tags/8/0/a5d4f41f098c85c8/1.html" target="_blank">模型</a></span></div> <div id="g0gggggg" class="article_love_news"><dd><a class="txt_title" href="/page/2024/0619/24295180.shtml" target="_blank" title="一半模型">一半模型</a><div id="g0gggggg" class="rsorc"><a href="/bk/thwgqmljtl/20245.html" class="ly" title="童話王國·奇妙邏輯推理(2024年5期)">童話王國·奇妙邏輯推理(2024年5期)</a><span id="g0gggggg" class="txt">2024-06-19 16:03:38</span></div></dd><dd><a class="txt_title" href="/page/2022/0829/12456558.shtml" target="_blank" title="一種去中心化的域名服務本地化模型">一種去中心化的域名服務本地化模型</a><div id="g0gggggg" class="rsorc"><a href="/bk/wlaqysjgl/20221.html" class="ly" title="網絡安全與數(shù)據(jù)管理(2022年1期)">網絡安全與數(shù)據(jù)管理(2022年1期)</a><span id="g0gggggg" class="txt">2022-08-29 03:15:20</span></div></dd><dd><a class="txt_title" href="/page/2022/0815/10688729.shtml" target="_blank" title="適用于BDS-3 PPP的隨機模型">適用于BDS-3 PPP的隨機模型</a><div id="g0gggggg" class="rsorc"><a href="/bk/dhdwxb/20224.html" class="ly" title="導航定位學報(2022年4期)">導航定位學報(2022年4期)</a><span id="g0gggggg" class="txt">2022-08-15 08:27:00</span></div></dd><dd><a class="txt_title" href="/page/2022/0614/13840119.shtml" target="_blank" title="提煉模型 突破難點">提煉模型 突破難點</a><div id="g0gggggg" class="rsorc"><a href="/bk/zxsslhzkb/20228.html" class="ly" title="中學生數(shù)理化·中考版(2022年8期)">中學生數(shù)理化·中考版(2022年8期)</a><span id="g0gggggg" class="txt">2022-06-14 06:55:24</span></div></dd><dd><a class="txt_title" href="/page/2021/1124/13632982.shtml" target="_blank" title="函數(shù)模型及應用">函數(shù)模型及應用</a><div id="g0gggggg" class="rsorc"><a href="/bk/xsjznsxbk/20219.html" class="ly" title="新世紀智能(數(shù)學備考)(2021年9期)">新世紀智能(數(shù)學備考)(2021年9期)</a><span id="g0gggggg" class="txt">2021-11-24 01:14:36</span></div></dd><dd><a class="txt_title" href="/page/2021/0719/12335968.shtml" target="_blank" title="p150Glued在帕金森病模型中的表達及分布">p150Glued在帕金森病模型中的表達及分布</a><div id="g0gggggg" class="rsorc"><a href="/bk/cdyxyxb/20212.html" class="ly" title="成都醫(yī)學院學報(2021年2期)">成都醫(yī)學院學報(2021年2期)</a><span id="g0gggggg" class="txt">2021-07-19 08:35:14</span></div></dd><dd><a class="txt_title" href="/page/2021/0104/13703595.shtml" target="_blank" title="函數(shù)模型及應用">函數(shù)模型及應用</a><div id="g0gggggg" class="rsorc"><a href="/bk/xsjznsxbk/20209.html" class="ly" title="新世紀智能(數(shù)學備考)(2020年9期)">新世紀智能(數(shù)學備考)(2020年9期)</a><span id="g0gggggg" class="txt">2021-01-04 00:25:14</span></div></dd><dd><a class="txt_title" href="/page/2020/1126/13925649.shtml" target="_blank" title="重要模型『一線三等角』">重要模型『一線三等角』</a><div id="g0gggggg" class="rsorc"><a href="/bk/zxsslhqnjsxrjb/202010.html" class="ly" title="中學生數(shù)理化·七年級數(shù)學人教版(2020年10期)">中學生數(shù)理化·七年級數(shù)學人教版(2020年10期)</a><span id="g0gggggg" class="txt">2020-11-26 08:24:50</span></div></dd><dd><a class="txt_title" href="/page/2020/0602/17732655.shtml" target="_blank" title="重尾非線性自回歸模型自加權M-估計的漸近分布">重尾非線性自回歸模型自加權M-估計的漸近分布</a><div id="g0gggggg" class="rsorc"><a href="/bk/sxwlxb/20202.html" class="ly" title="數(shù)學物理學報(2020年2期)">數(shù)學物理學報(2020年2期)</a><span id="g0gggggg" class="txt">2020-06-02 11:29:24</span></div></dd><dd><a class="txt_title" href="/page/2016/1107/18790514.shtml" target="_blank" title="3D打印中的模型分割與打包">3D打印中的模型分割與打包</a><div id="g0gggggg" class="rsorc"><a href="/bk/gxjmgc/20166.html" class="ly" title="光學精密工程(2016年6期)">光學精密工程(2016年6期)</a><span id="g0gggggg" class="txt">2016-11-07 09:07:19</span></div></dd></div> </div><div id="g0gggggg" class="other_pel mt80"> <p class="fl"><a href="/bk/spyswjsb/20179.html" target="_blank"><img src="https://cimg.fx361.com/images/2023/0407/26b03e6fba8e6ec9e4922e17b2ea7b03e206cec3.webp" alt=""></a><span id="g0gggggg" class="p1"><a href="/bk/spyswjsb/" target="_blank">食品與生物技術學報</a></span><span id="g0gggggg" class="p2"><a href="/bk/spyswjsb/20179.html" target="_blank">2017年9期</a></span></p> <dl class="fl"><dt>食品與生物技術學報的其它文章</dt><dd><a href="/page/2017/1101/20647410.shtml" title="牛羊乳蛋白質粒度對其熱處理沉淀率的影響">牛羊乳蛋白質粒度對其熱處理沉淀率的影響</a></dd><dd><a href="/page/2017/1101/20647389.shtml" title="復合酶解法優(yōu)化黃精多糖提取工藝">復合酶解法優(yōu)化黃精多糖提取工藝</a></dd><dd><a href="/page/2017/1101/20647356.shtml" title="艾蒿多糖抑菌活性及穩(wěn)定性">艾蒿多糖抑菌活性及穩(wěn)定性</a></dd><dd><a href="/page/2017/1101/20647306.shtml" title="黑曲霉葡萄糖氧化酶基因在畢赤酵母SMD1168中的表達">黑曲霉葡萄糖氧化酶基因在畢赤酵母SMD1168中的表達</a></dd><dd><a href="/page/2017/1101/20647261.shtml" title="油炸兔肉酶法預處理及改良劑配方優(yōu)化">油炸兔肉酶法預處理及改良劑配方優(yōu)化</a></dd><dd><a href="/page/2017/1101/20647149.shtml" title="鐵皮石斛莖部和葉部多糖的性質和活性">鐵皮石斛莖部和葉部多糖的性質和活性</a></dd></dl> </div></div> </div> </div> <div id="g0gggggg" class="sidebarR"> <!-- tab選項卡 --> <div id="g0gggggg" class="tab01 mb20"><div id="g0gggggg" class="tabArrow"></div><div id="g0gggggg" class="tabItem"><div id="g0gggggg" class="tabTit"><a href="#">雜志排行</a></div> <div id="g0gggggg" class="tabCont"><ol><li><p class="row01"><span id="g0gggggg" class="topNum">1</span><a href="/bk/sdjy/202410.html" class="row01a">《師道·教研》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/sdjy/202410.html">2024年10期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">2</span><a href="/bk/swyzhsby/202411.html" class="row01a">《思維與智慧·上半月》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/swyzhsby/202411.html">2024年11期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">3</span><a href="/bk/xdgyjjhxxh/20242.html" class="row01a">《現(xiàn)代工業(yè)經濟和信息化》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/xdgyjjhxxh/20242.html">2024年2期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">4</span><a href="/bk/wxxsyb/202410.html" class="row01a">《微型小說月報》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/wxxsyb/202410.html">2024年10期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">5</span><a href="/bk/gywsw/20241.html" class="row01a">《工業(yè)微生物》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/gywsw/20241.html">2024年1期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">6</span><a href="/bk/xl/20249.html" class="row01a">《雪蓮》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/xl/20249.html">2024年9期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">7</span><a href="/bk/sjbl/202421.html" class="row01a">《世界博覽》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/sjbl/202421.html">2024年21期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">8</span><a href="/bk/zxqyglykj/20246.html" class="row01a">《中小企業(yè)管理與科技》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/zxqyglykj/20246.html">2024年6期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">9</span><a href="/bk/xdsp/20244.html" class="row01a">《現(xiàn)代食品》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/xdsp/20244.html">2024年4期</a></span></p></li><li><p class="row01"><span id="g0gggggg" class="topNum">10</span><a href="/bk/wszyjy/202410.html" class="row01a">《衛(wèi)生職業(yè)教育》</a><span id="g0gggggg" class="row01_fr"><a href="/bk/wszyjy/202410.html">2024年10期</a></span></p></li></ol> </div></div> </div> </div> <div id="g0gggggg" class="clr"></div> </div> </div> <!--div class="advertisement"> </div--> <div id="g0gggggg" class="footer"> <p><a href="/aboutus/index.html">關于參考網</a></p> </div> <!-- <script> if ('serviceWorker' in navigator) { window.onload = function () { navigator.serviceWorker.register('/sw.js'); }; } </script> --> <script type="text/javascript" src="https://s1.pstatp.com/cdn/expire-1-M/jquery/3.4.0/jquery.min.js"></script> <script type="text/javascript" src="https://s1.pstatp.com/cdn/expire-1-M/sticky-kit/1.1.3/sticky-kit.min.js"></script> <script type="text/javascript"> document.write('<script src="https://cimg.fx361.com/cdn/w/index.js"><\/script>'); </script> <footer> <div class="friendship-link"> <p>感谢您访问我们的网站，您可能还对以下资源感兴趣：自贡牌麓投资有限公司</p> <a href="http://m.handmsg.com"title=999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看">999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看</a> <a href="/sitemap.xml">网站地图</a> <div style="position:fixed;left:-9000px;top:-9000px;"></div> <div class="friend-links"> </div> </div> </footer> 主站蜘蛛池模板： <a href="http://www.123shula.com" target="_blank">91免费观看视频</a>| <a href="http://www.cglgo.com" target="_blank">亚洲国产精品成人久久综合影院</a>| <a href="http://www.foleyflyer.com" target="_blank">亚洲一区二区日韩欧美gif</a>| <a href="http://www.hzwinjet.com" target="_blank">美女国产在线</a>| <a href="http://www.ymmob.com" target="_blank">a色毛片免费视频</a>| <a href="http://www.checkcashingloans.net " target="_blank">亚洲成人动漫在线观看</a>| <a href="http://www.shouxindc.com" target="_blank">囯产av无码片毛片一级</a>| <a href="http://www.jiejiong.com" target="_blank">国产成人欧美</a>| <a href="http://www.gangjinhui.com" target="_blank">亚洲精品免费网站</a>| <a href="http://www.51lbxx.net" target="_blank">久久a级片</a>| <a href="http://www.xiangshanspa.com" target="_blank">亚洲最新在线</a>| <a href="http://www.collect200.net" target="_blank">精品国产美女福到在线不卡f</a>| <a href="http://www.qdlikai.com" target="_blank">国产欧美在线观看一区</a>| <a href="http://www.quxiaoqi.com" target="_blank">国产91熟女高潮一区二区</a>| <a href="http://www.pzh4x.com" target="_blank">亚洲AV电影不卡在线观看</a>| <a href="http://www.jlhangkong.com" target="_blank">国产亚洲现在一区二区中文</a>| <a href="http://www.wasds.cc" target="_blank">97精品国产高清久久久久蜜芽</a>| <a href="http://www.lanh1.com" target="_blank">99成人在线观看</a>| <a href="http://www.9chi.net" target="_blank">四虎AV麻豆</a>| <a href="http://www.wlyxa.com" target="_blank">国内精自视频品线一二区</a>| <a href="http://www.jxdnzs.com" target="_blank">精品福利视频网</a>| <a href="http://www.9e555.com" target="_blank">亚洲黄色成人</a>| <a href="http://www.jdhyqd.com" target="_blank">又黄又湿又爽的视频</a>| <a href="http://www.inhyderabad.net" target="_blank">国产精品无码影视久久久久久久</a>| <a href="http://www.tthcp05.com" target="_blank">国产无码在线调教</a>| <a href="http://www.jiuzhaigo.com" target="_blank">久久99热这里只有精品免费看</a>| <a href="http://www.jieyali.com" target="_blank">伊人久综合</a>| <a href="http://www.yakamz.com" target="_blank">国产精品福利在线观看无码卡</a>| <a href="http://www.chuanjiuwang.com" target="_blank">中文字幕永久视频</a>| <a href="http://www.sxszgsyp.com" target="_blank">伊人AV天堂</a>| <a href="http://www.hgxsxx.com" target="_blank">午夜福利在线观看成人</a>| <a href="http://www.hsswzny.com" target="_blank">美女国产在线</a>| <a href="http://www.hgjinzun.com" target="_blank">亚洲 欧美 中文 AⅤ在线视频</a>| <a href="http://www.lysshop.cn" target="_blank">青草视频久久</a>| <a href="http://www.bnegpx.cn" target="_blank">欧美色图久久</a>| <a href="http://www.cyhzf.com" target="_blank">在线观看av永久</a>| <a href="http://www.al-pacabags.net" target="_blank">99在线观看精品视频</a>| <a href="http://www.halojk.com" target="_blank">三级国产在线观看</a>| <a href="http://www.yuyouke.com" target="_blank">精品国产成人国产在线</a>| <a href="http://www.gxsoben.com" target="_blank">国产精品久久久久无码网站</a>| <a href="http://www.act4health.net" target="_blank">日韩美女福利视频</a>| <a href="http://www.wanjuche.net" target="_blank">www欧美在线观看</a>| <a href="http://www.hzhugo.com" target="_blank">欧类av怡春院</a>| <a href="http://www.allerbind.net" target="_blank">2020久久国产综合精品swag</a>| <a href="http://www.aspnetxx.com" target="_blank">区国产精品搜索视频</a>| <a href="http://www.tswin.net" target="_blank">亚洲an第二区国产精品</a>| <a href="http://www.xmjjw.net" target="_blank">久久久91人妻无码精品蜜桃HD</a>| <a href="http://www.reddragonfly.net" target="_blank">亚洲人成网站在线播放2019</a>| <a href="http://www.linlinjia.com" target="_blank">白浆免费视频国产精品视频</a>| <a href="http://www.cbs258.com" target="_blank">国产精品一老牛影视频</a>| <a href="http://www.sds678.com" target="_blank">国产成人精品免费av</a>| <a href="http://www.yhz987.com" target="_blank">在线日韩日本国产亚洲</a>| <a href="http://www.starrylai.com" target="_blank">她的性爱视频</a>| <a href="http://www.njrycc.com" target="_blank">亚洲视屏在线观看</a>| <a href="http://www.612717.com" target="_blank">中文字幕亚洲电影</a>| <a href="http://www.adtinstall.net" target="_blank">国产成人精品男人的天堂下载 </a>| <a href="http://www.shuyousup.com" target="_blank">欧美激情第一欧美在线</a>| <a href="http://www.bdccz.com" target="_blank">99视频全部免费</a>| <a href="http://www.zfynzs.com" target="_blank">在线观看免费人成视频色快速</a>| <a href="http://www.i138.net" target="_blank">中美日韩在线网免费毛片视频</a>| <a href="http://www.adtinstall.net" target="_blank">国产亚洲精品精品精品</a>| <a href="http://www.nmcxtx.com" target="_blank">美女国内精品自产拍在线播放</a>| <a href="http://www.llelee.com" target="_blank">狠狠色噜噜狠狠狠狠色综合久</a>| <a href="http://www.ae-trade.com" target="_blank">91人妻在线视频</a>| <a href="http://www.guoranxing.com" target="_blank">国产大片黄在线观看</a>| <a href="http://www.mytastegr.com" target="_blank">丰满的少妇人妻无码区</a>| <a href="http://www.yyxhmp.com" target="_blank">亚洲av日韩综合一区尤物</a>| <a href="http://www.495898.com" target="_blank">波多野结衣一区二区三区AV</a>| <a href="http://www.jscyjj.com" target="_blank">a毛片在线播放</a>| <a href="http://www.mapfont.net" target="_blank">欧美日韩动态图</a>| <a href="http://www.968012.com" target="_blank">91啦中文字幕</a>| <a href="http://www.hk1717.com" target="_blank">亚洲高清中文字幕在线看不卡</a>| <a href="http://www.sztyad.com" target="_blank">免费人成网站在线观看欧美</a>| <a href="http://www.meireyu.com" target="_blank">日本人真淫视频一区二区三区</a>| <a href="http://www.lcbj18.com" target="_blank">一区二区欧美日韩高清免费</a>| <a href="http://www.aiyan123.com" target="_blank">国产亚洲一区二区三区在线</a>| <a href="http://www.absolutelyportable.com" target="_blank">久久精品国产电影</a>| <a href="http://www.71stock.com" target="_blank">国产永久免费视频m3u8</a>| <a href="http://www.xawbdwyy.com" target="_blank">色视频国产</a>| <a href="http://www.oveetrade.com" target="_blank">亚洲狼网站狼狼鲁亚洲下载</a>| <a href="http://www.un-fair.net" target="_blank">午夜欧美理论2019理论</a>| <a href="http://www.cb5000.com" target="_blank">色亚洲激情综合精品无码视频</a>| <script> (function(){ var bp = document.createElement('script'); var curProtocol = window.location.protocol.split(':')[0]; if (curProtocol === 'https') { bp.src = 'https://zz.bdstatic.com/linksubmit/push.js'; } else { bp.src = 'http://push.zhanzhang.baidu.com/push.js'; } var s = document.getElementsByTagName("script")[0]; s.parentNode.insertBefore(bp, s); })(); </script> </body><div id="4qqqq" class="pl_css_ganrao" style="display: none;"><dd id="4qqqq"></dd><ul id="4qqqq"></ul><source id="4qqqq"></source><tbody id="4qqqq"><xmp id="4qqqq"><cite id="4qqqq"></cite></xmp></tbody><center id="4qqqq"></center><center id="4qqqq"></center><noframes id="4qqqq"></noframes><nav id="4qqqq"><cite id="4qqqq"><option id="4qqqq"></option></cite></nav><wbr id="4qqqq"><noframes id="4qqqq"><small id="4qqqq"></small></noframes></wbr><tfoot id="4qqqq"></tfoot><pre id="4qqqq"></pre><noscript id="4qqqq"><optgroup id="4qqqq"><strong id="4qqqq"></strong></optgroup></noscript><menu id="4qqqq"><del id="4qqqq"><button id="4qqqq"></button></del></menu><noframes id="4qqqq"><noscript id="4qqqq"><optgroup id="4qqqq"></optgroup></noscript></noframes><cite id="4qqqq"></cite><source id="4qqqq"></source><dl id="4qqqq"></dl><kbd id="4qqqq"></kbd><noscript id="4qqqq"></noscript><pre id="4qqqq"></pre><strong id="4qqqq"></strong><dd id="4qqqq"><dfn id="4qqqq"><samp id="4qqqq"></samp></dfn></dd><em id="4qqqq"></em><s id="4qqqq"></s><center id="4qqqq"></center><input id="4qqqq"></input><optgroup id="4qqqq"></optgroup><input id="4qqqq"></input><tbody id="4qqqq"><nav id="4qqqq"><code id="4qqqq"></code></nav></tbody><cite id="4qqqq"></cite><kbd id="4qqqq"></kbd><center id="4qqqq"></center><source id="4qqqq"></source><del id="4qqqq"></del><samp id="4qqqq"></samp><pre id="4qqqq"><th id="4qqqq"><abbr id="4qqqq"></abbr></th></pre><delect id="4qqqq"></delect><li id="4qqqq"><table id="4qqqq"><small id="4qqqq"></small></table></li><menu id="4qqqq"></menu><rt id="4qqqq"></rt><th id="4qqqq"></th><pre id="4qqqq"></pre><li id="4qqqq"></li><wbr id="4qqqq"><noframes id="4qqqq"><small id="4qqqq"></small></noframes></wbr><code id="4qqqq"></code><acronym id="4qqqq"><dl id="4qqqq"><wbr id="4qqqq"></wbr></dl></acronym><input id="4qqqq"><acronym id="4qqqq"><ul id="4qqqq"></ul></acronym></input><tr id="4qqqq"><blockquote id="4qqqq"><center id="4qqqq"></center></blockquote></tr><dd id="4qqqq"><dfn id="4qqqq"><samp id="4qqqq"></samp></dfn></dd><del id="4qqqq"></del><small id="4qqqq"><menu id="4qqqq"><strong id="4qqqq"></strong></menu></small><wbr id="4qqqq"></wbr><pre id="4qqqq"><abbr id="4qqqq"><abbr id="4qqqq"></abbr></abbr></pre><s id="4qqqq"></s><kbd id="4qqqq"></kbd><em id="4qqqq"><tbody id="4qqqq"><nav id="4qqqq"></nav></tbody></em><tr id="4qqqq"></tr><dfn id="4qqqq"></dfn><optgroup id="4qqqq"><strong id="4qqqq"><s id="4qqqq"></s></strong></optgroup><table id="4qqqq"><code id="4qqqq"><object id="4qqqq"></object></code></table><dl id="4qqqq"></dl><wbr id="4qqqq"></wbr><sup id="4qqqq"><ul id="4qqqq"><option id="4qqqq"></option></ul></sup><noscript id="4qqqq"></noscript><xmp id="4qqqq"></xmp><strong id="4qqqq"></strong><s id="4qqqq"></s><optgroup id="4qqqq"></optgroup><tbody id="4qqqq"></tbody><dd id="4qqqq"></dd><pre id="4qqqq"><abbr id="4qqqq"><abbr id="4qqqq"></abbr></abbr></pre><optgroup id="4qqqq"><pre id="4qqqq"><del id="4qqqq"></del></pre></optgroup><table id="4qqqq"></table><tbody id="4qqqq"><nav id="4qqqq"><cite id="4qqqq"></cite></nav></tbody><kbd id="4qqqq"></kbd><blockquote id="4qqqq"></blockquote><option id="4qqqq"><tfoot id="4qqqq"><dd id="4qqqq"></dd></tfoot></option><center id="4qqqq"></center><tr id="4qqqq"></tr><pre id="4qqqq"><del id="4qqqq"><dl id="4qqqq"></dl></del></pre><strike id="4qqqq"></strike><option id="4qqqq"></option><xmp id="4qqqq"><nav id="4qqqq"><sup id="4qqqq"></sup></nav></xmp><option id="4qqqq"><noscript id="4qqqq"><optgroup id="4qqqq"></optgroup></noscript></option><samp id="4qqqq"></samp><sup id="4qqqq"><delect id="4qqqq"><pre id="4qqqq"></pre></delect></sup><th id="4qqqq"></th><th id="4qqqq"></th><wbr id="4qqqq"></wbr><small id="4qqqq"><menu id="4qqqq"><fieldset id="4qqqq"></fieldset></menu></small><ul id="4qqqq"></ul><ul id="4qqqq"><option id="4qqqq"><abbr id="4qqqq"></abbr></option></ul><tr id="4qqqq"><blockquote id="4qqqq"><center id="4qqqq"></center></blockquote></tr><li id="4qqqq"></li><pre id="4qqqq"><source id="4qqqq"><tbody id="4qqqq"></tbody></source></pre><ul id="4qqqq"></ul><small id="4qqqq"></small><wbr id="4qqqq"><noframes id="4qqqq"><small id="4qqqq"></small></noframes></wbr><delect id="4qqqq"><pre id="4qqqq"><kbd id="4qqqq"></kbd></pre></delect><rt id="4qqqq"></rt><th id="4qqqq"><abbr id="4qqqq"><td id="4qqqq"></td></abbr></th><fieldset id="4qqqq"></fieldset><center id="4qqqq"></center><fieldset id="4qqqq"><source id="4qqqq"><em id="4qqqq"></em></source></fieldset><abbr id="4qqqq"></abbr><delect id="4qqqq"></delect><tr id="4qqqq"></tr><ul id="4qqqq"><rt id="4qqqq"><th id="4qqqq"></th></rt></ul><em id="4qqqq"></em><li id="4qqqq"></li><pre id="4qqqq"></pre><dl id="4qqqq"><wbr id="4qqqq"><noframes id="4qqqq"></noframes></wbr></dl><strike id="4qqqq"></strike><sup id="4qqqq"></sup><small id="4qqqq"></small><th id="4qqqq"><abbr id="4qqqq"><td id="4qqqq"></td></abbr></th><samp id="4qqqq"><acronym id="4qqqq"><dl id="4qqqq"></dl></acronym></samp><blockquote id="4qqqq"></blockquote><tr id="4qqqq"></tr><noframes id="4qqqq"></noframes><dfn id="4qqqq"><bdo id="4qqqq"><input id="4qqqq"></input></bdo></dfn><strong id="4qqqq"></strong><center id="4qqqq"><del id="4qqqq"><s id="4qqqq"></s></del></center><dl id="4qqqq"></dl><optgroup id="4qqqq"></optgroup><acronym id="4qqqq"></acronym><ul id="4qqqq"></ul><object id="4qqqq"></object><center id="4qqqq"><fieldset id="4qqqq"><source id="4qqqq"></source></fieldset></center><cite id="4qqqq"><rt id="4qqqq"><option id="4qqqq"></option></rt></cite><dd id="4qqqq"><dfn id="4qqqq"><bdo id="4qqqq"></bdo></dfn></dd><source id="4qqqq"></source><small id="4qqqq"></small><menu id="4qqqq"></menu><noscript id="4qqqq"></noscript><bdo id="4qqqq"></bdo><acronym id="4qqqq"></acronym><abbr id="4qqqq"></abbr><xmp id="4qqqq"><sup id="4qqqq"><ul id="4qqqq"></ul></sup></xmp><input id="4qqqq"></input><input id="4qqqq"><acronym id="4qqqq"><ul id="4qqqq"></ul></acronym></input><del id="4qqqq"></del><center id="4qqqq"></center><del id="4qqqq"></del><dl id="4qqqq"></dl><center id="4qqqq"></center><wbr id="4qqqq"></wbr><optgroup id="4qqqq"></optgroup><tr id="4qqqq"></tr><pre id="4qqqq"></pre></div> </html>