CN103235902A - Prediction method for ribose nucleic acid (RNA) structure comprising false knots - Google Patents

Prediction method for ribose nucleic acid (RNA) structure comprising false knots Download PDF

Info

Publication number
CN103235902A
CN103235902A CN2013101363864A CN201310136386A CN103235902A CN 103235902 A CN103235902 A CN 103235902A CN 2013101363864 A CN2013101363864 A CN 2013101363864A CN 201310136386 A CN201310136386 A CN 201310136386A CN 103235902 A CN103235902 A CN 103235902A
Authority
CN
China
Prior art keywords
rna
structural unit
rna structure
energy
false knot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013101363864A
Other languages
Chinese (zh)
Other versions
CN103235902B (en
Inventor
刘振栋
张鹏
崔巍
张志军
李跃军
柳楠
徐功文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Jianzhu University
Original Assignee
Shandong Jianzhu University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong Jianzhu University filed Critical Shandong Jianzhu University
Priority to CN201310136386.4A priority Critical patent/CN103235902B/en
Publication of CN103235902A publication Critical patent/CN103235902A/en
Application granted granted Critical
Publication of CN103235902B publication Critical patent/CN103235902B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

The invention discloses a prediction method for a ribose nucleic acid (RNA) structure comprising false knots. The method comprises the steps of determining all structural units in a to-be-predicted RNA sequence including the false knots, placing all known existing structural units in a structural unit pool S0={s1, s2, s3,...sn}, enabling n to be a structural unit sum, and enabling sn to represent an nth structural unit; determining U={U1, U2,..., Ur, ..., UR} by iteration, enabling Ur to represent an RNA structure obtained by iteration in an rth time with small RNA structure energy, and enabling R to be total iteration times; respectively determining similar values of every element in Ur and an actual RNA structure according to free energy of every element in Ur and the sum of occurrence frequency of every element in all RNA structures; and enabling an element with a high similar value in U to be predicted as the RNA structure of the to-be-predicted RNA sequence. According to the prediction method, prediction time of the RNA structure and space complexity can be reduced, and prediction sensitivity and specificity are improved.

Description

The RNA structure prediction method that comprises false knot
Technical field
The invention belongs to biological information engineering field, relate to a kind of forecast method of carrying out to RNA (ribonucleic acid) (hereinafter, abbreviating RNA as), relate in particular to the RNA structure prediction method that comprises false knot.
Background technology
False knot (pseudoknot also claims pseudo-knot) is the nucleic acid tertiary structure that comprises at least two loop-stem structures, and wherein, half of one of two stem rings is inserted between the two halves of another stem ring.Nineteen eighty-two has been found false knot first in Brassica 2 et 4 (turnip yellow mosaic virus).False knot is folded into the 3 D stereo conformation of knot shape, but is not real topology knot.
In fact, in advance the general considerations of the minimum free energy RNA structure of measuring tape false knot has been proved to be np complete problem.Yet many important biological methods depend on the prediction to the RNA structure of band false knot.For example, telomerase RNA component (Telomerase RNA component is with reference to figure 1) comprises its active vital false knot.Many viruses are used the false knot structure to form similar tRNA motif (tRNA-like motif) and are penetrated into host cell.Having widely, the RNA molecule of tertiary structure often has a large amount of false knots.
Yet, because the context sensitivity (context-sensitivity) of false knot structure or the characteristic of " overlapping " are difficult to that it is carried out biological computation and detect.The base pairing of false knot is well not nested, in other words, and the base-pair appearance that in sequence, overlaps each other.This makes the standard method (standard method of dynamic programming) of existing dynamic programming be difficult to the false knot sequence in the default RNA (ribonucleic acid) (hereinafter, abbreviating RNA as).The newer stochastic context method (method of stochastic context-free grammars) that has nothing to do has also run into same problem.In secondary structure prediction methods such as more popular Mfold and Pfold, the false knot structure that exists in the RNA sequence that even can not go to predict.
Therefore, how as far as possible raising is a science difficult problem that needs to be resolved hurrily to susceptibility and the specificity of the prediction of the RNA structural series of band false knot.
Summary of the invention
One of technical matters to be solved by this invention is to provide a kind of RNA structure to comprising false knot that can reduce time and space complexity to carry out forecast method.
In order to solve the problems of the technologies described above, the invention provides a kind of RNA structure prediction method that comprises false knot.This method comprises: step S10, determine all structural units in the RNA sequence to be predicted, and comprise false knot, the structural unit of all known existence is put into structural unit pond S 0={ s 1, s 2, s 3... s n, n is the structural unit sum, s nRepresent n structural unit; Step S20 based on all structural units in the RNA sequence to be predicted, determines U={ U by iteration 1, U 2..., U r..., U R, U rRepresent the less RNA structure of RNA structural energy that the r time iteration obtains, R is total iterations; Step S30 is according to U rIn each element free energy and in the RNA structure summation of the structural unit frequency of occurrences, determine U respectively rIn the similar value of each element and actual RNA structure; Step S40, the element that similar value among the U is high is predicted as the RNA structure of this RNA sequence to be predicted.
Wherein, described based on all structural units in the RNA sequence to be predicted, determine U={ U by iteration 1, U 2..., U r..., U RStep, further comprise: substep S201, to structural unit U rCompose null value, U rRepresent the folding RNA structure that obtains the r time; K is initialized as 1 with iterations; Substep S202 carries out K time by substep S203 to substep S204 and folds processing; Substep S203 judges the structural unit s among the S one by one iWhether can make current RNA structure U rFree energy reduce and s iWith U rIn structural unit not overlapping, if be judged as and be, K=K+1 then is with s iIncorporate U into r, return substep S202 then folding next time, otherwise, if be judged as not, then continue to judge the s among the S I+1Whether can make current RNA structure U rFree energy reduce and s iWith U rIn structural unit not overlapping, up to all elements among the S is judged; Substep S204 if all elements among the S is judged, has not existed among the expression S and can reduce RNA structure U rThe structural unit of free energy, then this iteration finishes.
Further, described substep S203 utilizes following expression formula to calculate current RNA structure U rFree energy:
E=E pseudo+E nest
E Nest=E Multiloop+ E Stem+ E Bulge+ E Interloop+ E Hairpin+ E A, wherein,
E represents the free energy of RNA structure; E NestThe energy of nested structure in the expression RNA structure, E MultiloopBe the energy of multi-branched ring in the RNA structure, E StemBe the energy in stem district in the RNA structure, E BulgeBe the energy of RNA structure protrusions, E InterloopBe the energy of interior ring in the RNA structure, E HairpinBe the energy of hair fastener ring in the RNA structure, E AThe compensation that is each structural unit connects parameter; E PseudoThe energy of false knot in the expression RNA structure.
Further, calculate RNA structure U by following expression rThe energy of middle false knot:
E pseudo=A 1S b+A 2P b+A 3N b+A 4+A 5
A 1Be the energy value that produces a false knot structure neutron false knot, S bBe the number of false knot structure neutron false knot, P bBe the pairing base-pair number on the false knot inner boundary, N bBe the inner unpaired base number of false knot, A 2, A 3Be respectively P bAnd N bEnergy value, can be measured A by experiment 4Be the coaxial stacking energy weights in the false knot, A 5Be the connection parameter between the RNA structural unit, connect parameter and represent the modified value that is connected between RNA false knot and the RNA structural unit.
Preferably, among the step S30, calculate the similar value of each element and actual RNA structure among the U according to following formula:
Function(U r)=E(U r)×a/n+P(U r)×b+c,r∈[1,k]
Wherein,
Function (U r) expression U rSimilar value with actual RNA structure;
E (U r) expression structural unit U rFree energy;
A, b, c, k are constants;
N represents RNA base sequence length;
P (U r) expression U rIn the summation of structural unit frequency of occurrences in the RNA structure.
Preferably, a ∈ [5,15], b ∈ [1,8], c ∈ [1,10], k ∈ [90,110].Optimally, a=10, b=3, c=5, k=100.
Compared with prior art, one or more embodiment of the present invention can have following advantage:
In the present invention, proposed to comprise the Forecasting Methodology of the RNA structure of false knot, this method is by utilizing hereinafter with the approximate data that describes in detail, can improve the RNA structure prediction accuracy that comprises false knot preferably, improve specificity and the susceptibility of prediction, prediction RNA structure is conformed to or close with practical structures.
In addition, this invention can not only be predicted the RNA secondary structure, even also has specificity and susceptibility preferably aspect the prediction RNA tertiary structure.
Other advantages of the present invention, target, to set forth in the following description to a certain extent with feature, and to a certain extent, based on being apparent to those skilled in the art to investigating hereinafter, perhaps can obtain instruction from the practice of the present invention.Target of the present invention and other advantages can be passed through following instructions, claims, and the specifically noted structure realizes and obtains in the accompanying drawing.
Description of drawings
Accompanying drawing is used to provide further understanding of the present invention, and constitutes the part of instructions, is used for jointly explaining the present invention with embodiments of the invention, is not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is telomerase RNA component structure synoptic diagram;
Fig. 2 A illustrates the process flow diagram of structure that method is according to an embodiment of the invention determined the RNA of band false knot;
Fig. 2 B illustrates the process flow diagram of the detailed step of step S20 among Fig. 2 A;
Fig. 3 is time complexity and the space complexity comparison diagram that the inventive method and Pknots technology are preset these several RNA structures of mottle virus, mosaic virus, ingspot virus and sapiens mRNA;
Fig. 4 be the inventive method and ILM technology, MWM technology, Pknots technology respectively to multiple RNA sequence (HIV-1-RT, DR7741, TYMV, BMV1, TMV-3 '-up, TMV-3 '-down, HDV, Anti-HDV) when predicting, the prediction effect comparison diagram aspect susceptibility, specificity two.
Embodiment
Describe embodiments of the present invention in detail below with reference to drawings and Examples, how the application technology means solve technical matters to the present invention whereby, and the implementation procedure of reaching technique effect can fully understand and implements according to this.Need to prove that only otherwise constitute conflict, each embodiment among the present invention and each feature among each embodiment can mutually combine, formed technical scheme is all within protection scope of the present invention.
In addition, can in the computer system such as one group of computer executable instructions, carry out in the step shown in the process flow diagram of accompanying drawing, and, though there is shown logical order in flow process, but in some cases, can carry out step shown or that describe with the order that is different from herein.
Illustrate that at first some are about the concept in false knot and stem district etc.False knot (pseudoknot) is as by intersecting the structure that base-pair constitutes, and is tertiary structure unit the most widely among the RNA, is complicated but stable RNA structure, and the plane false knot is false knot subclass the most widely.By contrast, the stem district is made of adjacent continuous base-pair.
Present embodiment utilizes following formula to determine the energy of a false knot:
E pseudo=A 1S b+A 2P b+A 3N b+A 4+A 5
A wherein 1Be the energy value that produces a false knot structure neutron false knot, can be measured by experiment.S bIt is the number of false knot structure neutron false knot; P bIt is the pairing base-pair number on the false knot inner boundary; N bIt is the inner unpaired base number of false knot; A 2, A 3Be respectively P bAnd N bEnergy value, can be measured by experiment; A 4Be coaxial stacking energy (coaxial stacking) weights in the false knot, A 5Be the connection parameter that (comprises for example hair fastener ring, multi-branched ring, interior ring, projection, stem district etc.) between the RNA structural unit, measured by experiment.Connect parameter and represent the modified value that is connected between RNA false knot and the RNA structural unit, different RNA false knot structures may connect also difference of parameter.
Present embodiment can utilize following expression to determine the energy of nested structure:
E nest=E multiloop+E stem+E bulge+E interloop+E hairpin+E A
E wherein MultiloopBe the energy of multi-branched ring in the RNA structure, E StemBe the energy in stem district in the RNA structure, E BulgeBe the energy of RNA structure protrusions, E InterloopBe the energy of interior ring in the RNA structure, E HairpinBe the energy of hair fastener ring in the RNA structure, E AThe compensation that is each structural unit connects parameter, can be determined by experiment.
Like this, present embodiment can utilize following expression to determine the gross energy (free energy) of a RNA structure:
E=E pseudo+E nest
Use which kind of expression formula to determine the false knot energy, be directly connected to susceptibility and the specificity of Forecasting Methodology.Above-mentioned expression formula in conjunction with following each step that will illustrate, can be improved susceptibility and the specificity of band false knot RNA prediction preferably.
Describe each step of structure of determining the RNA of band false knot according to present embodiment in detail below in conjunction with Fig. 2 A.
Step S10 determines all structural units in the RNA sequence to be predicted, comprises false knot, and the structural unit of all known existence is put into structural unit pond S 0={ s 1, s 2, s 3... s n, n is the structural unit sum, s nRepresent n structural unit.
Step S20 based on all structural units in the RNA sequence to be predicted, determines U={ U by iteration 1, U 2..., U r..., U R, U rRepresent the less RNA structure of RNA structural energy that the r time iteration obtains, R is total iterations, U 0=Φ;
Step S30, according to the free energy of each element among the U and in all RNA structures the summation of the structural unit frequency of occurrences, determine the similar value of each element and actual RNA structure among the U respectively;
Step S40, the element that similar value among the U is high is predicted as the RNA structure of this RNA sequence to be predicted.
Describe each step of the r time iteration among the step S20 in detail below with reference to Fig. 2 B.
Step S201, initialization structure (tax empty set) U r← Φ, U rRepresent the folding RNA structure that obtains the r time; Iterations K is initialized as 1, K ← 1.
Step S202 carries out folding the processing K time by step S203 to S204, and the value of folding number of times K is different because of the different RNA base sequence, can arrange flexibly according to requirements such as susceptibility and time complexities.Carrying out K folding meaning is: reach steady state (SS) by folding RNA base sequence.Judge that whether folding times is above K time.
Step S203 judges the structural unit s among the S one by one iWhether can make current RNA structure U rFree energy reduce and s iWith U rIn structural unit not overlapping, namely whether make E (U r∪ { s i)<E (U r) and If be judged as and be, K=K+1 then is with s iIncorporate U into r, i.e. U r← U r∪ { s i, return step S202 then folding next time, otherwise, if be judged as not, then continue to judge the s among the S I+1Whether can make current RNA structure U rFree energy reduce and s I+1With U rIn structural unit not overlapping up to all elements among the S is judged.
Substep S204 if all elements among the S is judged, has not existed among the expression S and can reduce RNA structure U rThe structural unit of free energy, then this iteration finishes.
Describe step S30 below in detail, according to the free energy of each element among the U and in all RNA structures the summation of the frequency of occurrences, determine the similar value of each element and actual RNA structure among the U respectively.
For example, can calculate the similar value of each element and actual RNA structure among the U according to following formula:
Function (U r)=E (U r) * a/n+P (U r) * b+c, r ∈ [1, k], Function (U r) expression U rWith the similar value of actual RNA structure, above-mentioned expression formula can make the RNA structure U that the more high expression of functional value is pre- rMore near practical structures.Wherein:
U rRepresent the folding RNA structure that obtains the r time
E (U r) expression structural unit U rFree energy
A, b, c, k are constants, can be preferably a ∈ [5,15] among the present invention by measuring, b ∈ [1,8], c ∈ [1,10], k ∈ [90,110], more preferably.In Fig. 4 experiment, a=10; B=3; C=5; K=100.
N represents RNA base sequence length;
P (U r) expression U rIn the summation of structural unit frequency of occurrences in the RNA structure.
Above-mentioned expression formula formula can make when selecting final RNA structure, both considered the factor of free energy minimum, also consider simultaneously the frequency of occurrences of structural unit. in other words, if a structural unit occurs, can think that then this structural unit also probably also can appear in the final RNA structure in a plurality of Local Search results.
The method that we propose is studied various RNA structures as a kind of structural unit unification, according to the contribution of the probability that adds free energy that these structural units reduce and appearance in the RNA folding process as the RNA structure, predict the RNA structure.Through RNA structures such as for example mottle virus, the mosaic virus that comprise false knot, ringspot virus, sapiens mRNA, HIV-1-RT, DR7741 are predicted, higher forecasting efficiency and accuracy (with reference to figure 3 and Fig. 4) have been obtained.
Fig. 3 is time complexity and the space complexity comparison diagram that the inventive method and Pknots technology are preset these several RNA structures of mottle virus, mosaic virus, ingspot virus and sapiens mRNA, and n represents that n is the number of RNA sequence base among Fig. 3.As can be seen from Figure 3, the present invention has clear superiority at time complexity.Fig. 4 be the inventive method and ILM technology, MWM technology, Pknots technology respectively to multiple RNA sequence (HIV-1-RT, DR7741, TYMV, BMV1, TMV-3 '-up, TMV-3 '-down, HDV, Anti-HDV) when predicting, the prediction effect comparison diagram aspect susceptibility, specificity two.As can be seen from Figure 4, to compare ILM technology, MWM technology and Pknots technology aspect susceptibility and specificity all relative more excellent for method of the present invention.
Though the disclosed embodiment of the present invention as above, the embodiment that described content just adopts for the ease of understanding the present invention is not in order to limit the present invention.Technician in any the technical field of the invention; under the prerequisite that does not break away from the disclosed spirit and scope of the present invention; can do any modification and variation in the details of implementing that reaches in form; but scope of patent protection of the present invention still must be as the criterion with the scope that appending claims was defined.

Claims (7)

1. a RNA structure prediction method that comprises false knot is characterized in that, comprising:
Step S10 determines all structural units in the RNA sequence to be predicted, comprises false knot, and the structural unit of all known existence is put into structural unit pond S 0={ s 1, s 2, s 3... s n, n is the structural unit sum, s nRepresent n structural unit;
Step S20 based on all structural units in the RNA sequence to be predicted, determines U={ U by iteration 1, U 2..., U r..., U R, U rRepresent the less RNA structure of RNA structural energy that the r time iteration obtains, R is total iterations;
Step S30 is according to U rIn each element free energy and in the RNA structure summation of the structural unit frequency of occurrences, determine U respectively rIn the similar value of each element and actual RNA structure;
Step S40, the element that similar value among the U is high is predicted as the RNA structure of this RNA sequence to be predicted.
2. method according to claim 1 is characterized in that, and is described based on all structural units in the RNA sequence to be predicted, determines U={ U by iteration 1, U 2..., U r..., U RStep, further comprise:
Substep S201 is to structural unit U rCompose null value, U rRepresent the folding RNA structure that obtains the r time; K is initialized as 1 with iterations;
Substep S202 carries out K time by substep S203 to substep S204 and folds processing;
Substep S203 judges the structural unit s among the S one by one iWhether can make current RNA structure U rFree energy reduce and s iWith U rIn structural unit not overlapping, if be judged as and be, K=K+1 then is with s iIncorporate U into r, return substep S202 then folding next time, otherwise, if be judged as not, then continue to judge the s among the S I+1Whether can make current RNA structure U rFree energy reduce and s iWith U rIn structural unit not overlapping, up to all elements among the S is judged;
Substep S204 if all elements among the S is judged, has not existed among the expression S and can reduce RNA structure U rThe structural unit of free energy, then this iteration finishes.
3. method according to claim 2 is characterized in that, described substep S203 utilizes following expression formula to calculate current RNA structure U rFree energy:
E=E pseudo+E nest
E Nest=E Multiloop+ E Stem+ E Bulge+ E Interloop+ E Hairpin+ E A, wherein,
E represents the free energy of RNA structure;
E NestThe energy of nested structure in the expression RNA structure, E MultiloopBe the energy of multi-branched ring in the RNA structure, E StemBe the energy in stem district in the RNA structure, E BulgeBe the energy of RNA structure protrusions, E InterloopBe the energy of interior ring in the RNA structure, E HairpinBe the energy of hair fastener ring in the RNA structure, E AThe compensation that is each structural unit connects parameter;
E PseudoThe energy of false knot in the expression RNA structure.
4. method according to claim 3 is characterized in that, calculates RNA structure U by following expression rThe energy of middle false knot:
E pseudo=A 1S b+A 2P b+A 3N b+A 4+A 5
A 1Be the energy value that produces a false knot structure neutron false knot, S bBe the number of false knot structure neutron false knot, P bBe the pairing base-pair number on the false knot inner boundary, N bBe the inner unpaired base number of false knot, A 2, A 3Be respectively P bAnd N bEnergy value, can be measured A by experiment 4Be the coaxial stacking energy weights in the false knot, A 5Be the connection parameter between the RNA structural unit, connect parameter and represent the modified value that is connected between RNA false knot and the RNA structural unit.
5. according to each described method in the claim 1 to 4, it is characterized in that, among the step S30, calculate the similar value of each element and actual RNA structure among the U according to following formula:
Function(U r)=E(U r)×a/n+P(U r)×b+c,r∈[1,k]
Wherein,
Function (U r) expression U rSimilar value with actual RNA structure;
E (U r) expression structural unit U rFree energy;
A, b, c, k are constants;
N represents RNA base sequence length;
P (U r) expression U rIn the summation of structural unit frequency of occurrences in the RNA structure.
6. method according to claim 5 is characterized in that, a ∈ [5,15], b ∈ [1,8], c ∈ [1,10], k ∈ [90,110].
7. method according to claim 6 is characterized in that, a=10, b=3, c=5, k=100.
CN201310136386.4A 2013-04-18 2013-04-18 Comprise the RNA Structure Prediction Methods of false knot Expired - Fee Related CN103235902B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310136386.4A CN103235902B (en) 2013-04-18 2013-04-18 Comprise the RNA Structure Prediction Methods of false knot

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310136386.4A CN103235902B (en) 2013-04-18 2013-04-18 Comprise the RNA Structure Prediction Methods of false knot

Publications (2)

Publication Number Publication Date
CN103235902A true CN103235902A (en) 2013-08-07
CN103235902B CN103235902B (en) 2016-03-09

Family

ID=48883942

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310136386.4A Expired - Fee Related CN103235902B (en) 2013-04-18 2013-04-18 Comprise the RNA Structure Prediction Methods of false knot

Country Status (1)

Country Link
CN (1) CN103235902B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104298894A (en) * 2014-09-17 2015-01-21 山东建筑大学 Prediction method and device for RNA false knot structure based on k stems
CN104765983A (en) * 2015-04-23 2015-07-08 山东建筑大学 Predicting method and device of ribonucleic pseudoknot structure based on half-extension structure
CN116092575A (en) * 2023-02-03 2023-05-09 中国科学院地理科学与资源研究所 G DNA structure discrimination method based on GMNS rule

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040030111A1 (en) * 2002-06-19 2004-02-12 Turner Douglas H. Oligonucleotide directed misfolding of RNA
US20060269911A1 (en) * 2004-09-16 2006-11-30 Avi Biopharma, Inc. Antisense antiviral compound and method for treating ssRNA viral infection
CN101717822A (en) * 2009-12-16 2010-06-02 吉林大学 RNA sequence secondary structure prediction method based on base fragment coding and ant colony optimization
CN101908102A (en) * 2010-08-13 2010-12-08 山东建筑大学 Ribosomal stalk based predicting method and device of RNA (Ribonucleic Acid) secondary structure
CN102880811A (en) * 2012-10-24 2013-01-16 吉林大学 Method for predicting secondary structure of ribonucleic acid (RNA) sequence based on complex programmable logic device (CPLD) base fragment encoding and ant colony algorithm

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040030111A1 (en) * 2002-06-19 2004-02-12 Turner Douglas H. Oligonucleotide directed misfolding of RNA
US20060269911A1 (en) * 2004-09-16 2006-11-30 Avi Biopharma, Inc. Antisense antiviral compound and method for treating ssRNA viral infection
CN101717822A (en) * 2009-12-16 2010-06-02 吉林大学 RNA sequence secondary structure prediction method based on base fragment coding and ant colony optimization
CN101908102A (en) * 2010-08-13 2010-12-08 山东建筑大学 Ribosomal stalk based predicting method and device of RNA (Ribonucleic Acid) secondary structure
CN102880811A (en) * 2012-10-24 2013-01-16 吉林大学 Method for predicting secondary structure of ribonucleic acid (RNA) sequence based on complex programmable logic device (CPLD) base fragment encoding and ant colony algorithm

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ZHENDONG LIU: "《New Heuristic Algorithm of RNA Structure Prediction including Pseudoknots》", 《JOURNAL OF COMPUTERS》 *
陈翔,卜东波, 张法, 高文: "基于局部茎搜索的RNA二级结构预测算法", 《生物化学与生物物理进展》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104298894A (en) * 2014-09-17 2015-01-21 山东建筑大学 Prediction method and device for RNA false knot structure based on k stems
CN104298894B (en) * 2014-09-17 2017-10-20 山东建筑大学 The Forecasting Methodology and device of ribonucleic acid pseudoknot structure based on k stems
CN104765983A (en) * 2015-04-23 2015-07-08 山东建筑大学 Predicting method and device of ribonucleic pseudoknot structure based on half-extension structure
CN116092575A (en) * 2023-02-03 2023-05-09 中国科学院地理科学与资源研究所 G DNA structure discrimination method based on GMNS rule

Also Published As

Publication number Publication date
CN103235902B (en) 2016-03-09

Similar Documents

Publication Publication Date Title
EP3061022B1 (en) Systems and methods for using paired-end data in directed acyclic structure
Nazareth Conjugate gradient method
Guo et al. A parallel attractor finding algorithm based on Boolean satisfiability for genetic regulatory networks
Zheng et al. GenRev: exploring functional relevance of genes in molecular networks
Shih et al. A single source k-shortest paths algorithm to infer regulatory pathways in a gene network
CN103500250B (en) Automatic layout method for digraph created GMF graph editor
CN104199852A (en) Label propagation community structure mining method based on node membership degree
Li et al. The normalized Laplacian spectrum of quadrilateral graphs and its applications
Kováč Complexity of the path avoiding forbidden pairs problem revisited
CN103235902A (en) Prediction method for ribose nucleic acid (RNA) structure comprising false knots
Willson Regular networks can be uniquely constructed from their trees
Tache General sum-connectivity index with α≥ 1 for bicyclic graphs
Brankovic et al. Linear-time superbubble identification algorithm for genome assembly
Zmazek et al. The obnoxious center problem on weighted cactus graphs
Yerneni et al. IAS: Interaction specific GO term associations for predicting Protein-Protein Interaction Networks
Ali et al. A note on chemical trees with minimum Wiener polarity index
CN108052743B (en) Method and system for determining step approach centrality
US9619362B2 (en) Event sequence construction of event-driven software by combinational computations
Pham et al. Graph-based information diffusion method for prioritizing functionally related genes in protein-protein interaction networks
Liu et al. Predicting scheme of RNA folding structure including pseudoknots
Liu et al. An ant colony optimization based algorithm for identifying gene regulatory elements
Yang et al. Finding K shortest looping paths with waiting time in a time–window network
Wiedenhoeft et al. The plexus model for the inference of ancestral multidomain proteins
Backofen et al. Fast detection of common sequence structure patterns in RNAs
Junior et al. FB-APSP: A new efficient algorithm for computing all-pairs shortest-paths

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20130807

Assignee: Nuoxin medical technology (Shandong) Co.,Ltd.

Assignor: SHANDONG JIANZHU University

Contract record no.: X2020370000011

Denomination of invention: RNA structure prediction method with false knots

Granted publication date: 20160309

License type: Common License

Record date: 20200821

EE01 Entry into force of recordation of patent licensing contract
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160309

CF01 Termination of patent right due to non-payment of annual fee