• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Statistical potentials for 3D structure evaluation:From proteins to RNAs?

    2021-03-11 08:34:44YaLanTan譚雅嵐ChenJieFeng封晨潔XunxunWang王勛勛WenbingZhang張文炳andZhiJieTan譚志杰
    Chinese Physics B 2021年2期

    Ya-Lan Tan(譚雅嵐), Chen-Jie Feng(封晨潔), Xunxun Wang(王勛勛),Wenbing Zhang(張文炳), and Zhi-Jie Tan(譚志杰)

    Department of Physics and Key Laboratory of Artificial Micro&Nano-structures of Education,School of Physics and Technology,Wuhan University,Wuhan 430072,China

    Keywords: statistical potential,3-dimensional structure evaluation,RNA,protein

    1. Introduction

    It is widely known that RNAs can pass the genetic information from DNAs to control the synthesis of proteins and can also serve as the storage of genetic information.[1]Beyond the sequence-level functions, a variety of non-coding RNAs have been discovered in recent three decades, which have extensively biological functions such as catalysis and genetic regulation.[2,3]The functions of noncoding RNAs are generally coupled to their structures or proper structure change. For example,ribozymes can only perform the catalysis of splicing or cleavage when they fold into their native 3-dimensional(3D)structures,[4]and riboswitches can only regulate the genetic expression at translation or transcription level through proper change of their 3D structures upon metabolite binding.[5]Therefore, understanding RNA structures, especially RNA 3D structures, would be critical to understanding their biological functions. Experimental methods have been widely employed to derive RNA 3D structures, such as x-ray crystallography, NMR spectroscopy, and cryo-electron.[6,7]However,until now,RNA 3D structures deposited in the protein data bank(PDB)database[8]are still limited,partially because the experimental methods are generally time-consuming and of high cost. In parallel, various computational models have been developed to predict RNA 3D structures in silico recently,[9–12]including knowledge-based and physics-based models.[13–36]Generally,a predictive model would predict an ensemble of structure candidates, and thus, a reliable scoring function becomes highly required to evaluate the quality of the candidate structures. Furthermore,a reliable scoring function can also be employed to guide RNA structure prediction and to optimize predicted structures.[37,38]

    For proteins,statistical potentials,also called knowledgebased potentials or mean force potentials, have been proven to be an effective and efficient scoring function for structure evaluation in protein and protein-complex 3D structure predictions, along with the protein folding problem from three decades ago.[39–48]In the recent decade,several statistical potentials have been developed for RNA structure evaluation or prediction,[49–56]most of which were built in similar ways to those for proteins due to the prior progress made for proteins and some similar polymer features between proteins and RNAs such as backbone and side chain. Nevertheless, there is inherent difference between proteins and RNAs, such as residue types and driving forces for structure folding. Therefore,to build high-performance statistical potentials for RNAs may still require the involvement of the structure characteristics of RNAs,in addition to utilizing similar modeling experience for proteins.

    In this review, we will firstly introduce the derivation of general formulas for statistical potentials based on fundamental principles. Afterwards, we will give a brief overview on typical traditional statistical potentials for protein 3D structure prediction and evaluation due to the important reference role for RNAs. After that, we will introduce the recent advances in the statistical potentials for RNA 3D structure evaluation.Finally,we will emphasize a perspective on the further development of new statistical potentials with higher performance for RNA 3D structure evaluation.

    2. Fundamental principles for statistical potentials

    There are two fundamental principles for deriving and utilizing a statistical potential for proteins or RNAs: (i) the probability of a conformation of a biomacromolecule in a conformation ensemble obeys the Boltzmann’s law which connects the effective potential energy of the system to its probability;[39,57,58](ii) the native structure of a biomacromolecule corresponds to the conformation with the lowest free energy,which was first demonstrated by Anfinsen in 1973[59]and connects the free energy minimum of a biomacromolecule to its native structure.[58]

    In addition, a native folded structure is actually driven and stabilized by various complex intramolecular and intermolecular interactions, while there would not exist a perfect potential energy function to exactly describe all these interactions for a biomacromolecule system.[58]Practically, a biomacromolecule system can be considered as an equilibrium ensemble of different atom(atom block)pair types with different types of characterizing geometrical parameters such as inter-atom contact, inter-atom distance, inter-atom angle,inter-block orientation,and so on. If a geometrical parameter is chosen,the probability of describing different values of the geometrical parameter can be utilized to build a statistical potential as a function of the geometrical parameter for different atom(or atom block)pair types.

    2.1. Derivation of statistical potential

    In principle,any kind of geometrical parameters,such as distances or angles between atoms which can be utilized to distinguish a native conformation from decoy ones, can be adopted to derive a statistical potential.[48]According to the inverse of Boltzmann’s law, the potential energy for a particular atom pair type in the native state ensemble of biomacromolecules can be expressed by[40]

    where kBand T are the Boltzmann constant and the temperature in Kelvin, respectively. Eobs(s) and Pobs(s) are the energy and observed probability of a particular atom pair type for a geometrical parameter s in the native structure ensemble.Zobsis the summation of Boltzmann factor over geometrical parameter s in the native state ensemble

    However,Eobs(s)includes not only the interactions we focus on but also other interactions from surrounding medium.[60]Therefore, to obtain an effective statistical potential ?E(s), a reference state needs to be involved. Similarly, according to Eq. (1), the inverse of Boltzmann’s law for a reference state can be written as[40]

    where Eref(s) and Pref(s) are the potential energy and the probability of the geometrical parameter s in a reference state,and Zrefis the summation of Boltzmann factor over geometrical parameter s in the reference state ensemble

    Therefore, a net statistical potential is obtained as a function of geometrical parameter s[40]

    and from Eqs. (1) and (3), ?E(s) can be computed from the probability Pobs(s)and Pref(s)by

    In addition, at a fixed temperature, Zobs/Zrefis a constant and independent of the geometrical parameter s for a given amino/nucleic acid sequence,[40]and consequently?kBTln[Zobs/Zref] can be treated as zero point of the potential. Therefore, a general expression of a statistical potential can be obtained as[40]

    As shown in Eq.(7),the geometrical parameters s,which may involve two or more atoms (or atom group) and the reference states,are crucial for building statistical potentials,and the core difference between various statistical potentials is attributed to their choice.

    2.2. Utilization of statistical potential

    After deriving a statistical potential from Eqs. (1)–(7),the total energy ?E(S,C)for a conformation C of a given sequence S relatively to the reference state can be given by[40]

    where the summation is over all atom(atom block)pairs with the additive assumption for statistical potentials.[40]

    3. Statistical potentials for protein 3D structures

    Since the pioneering work of Tanaka and Scheraga in 1976,[39]various strategies have been proposed to develop effective statistical potentials for proteins at atom or residue level.[40,42,43,61–64]These strategies can be roughly classified into the following categories: (i)developing various reference states or circumvent reference states;[43,65–70](ii) classifying short-, medium-, and long-range interactions;[39,57,72–74](iii)considering different geometrical parameters such as multibody[63,75–79]and orientation-dependent interactions.[80–86]The geometrical parameters and major features of existing approaches employed in developing statistical potentials for proteins have been summarized in Table 1.

    Table 1. The mainly existing approaches in developing statistical potentials for proteins.

    3.1. Existing widely used reference states

    The existing reference states can be roughly classified into two types: based on experimental structures in PDB database[43,65,68]and based on statistical physics models.[66,67,69]The former type includes averaging reference state,quasi-chemical approximation reference state,and atomshuffled reference state, and the latter type includes finiteideal-gas reference state, spherical-non-interacting reference state,and random-walk-chain reference state.

    3.1.1. On the basis of experimental structure database

    where Nobs(r)is the number of atom pairs within distance interval of [r,r+dr] from experimental structures in database,and Nobsis the summation of Nobs(r) over all distance intervals.

    3.3. Classification of short-, medium-, and long-range interactions

    Tanaka and Scheraga have classified the inter-residue interactions into short-, medium-, and long-ranged ones for proteins,[39]and the criteria for interaction ranges of proteins are based on the separations between residues along protein sequences. Afterwards,the preference of amino acids to form short-, medium-, and long-range contacts has been analyzed for different types of proteins,[72–74,89]and several potentials have also been proposed based on distinguishing interaction ranges.[57,72,91–94]The classification of interaction ranges may be very necessary since such treatment can involve more accurate information extracted from native structure database in different interaction ranges.

    3.3.1. Backbone torsion-based potentials

    To measure the influence of amino acids on the backbone conformation of neighboring residues along protein sequences, backbone torsion-dependent potentials were proposed by Rooman et al.[91]and Kocher et al.[92]A succession of dihedral angles φ, ψ, and ω of each residue represent the backbone of proteins, and these dihedral angles can be clustered into 7 domains.[91]Based on Eqs. (7) and (8), the total torsional energy for a protein conformation can be expressed as the summation of residue-to-torsion potentials along the protein conformation[92]

    3.4. Extraction of multi-body interactions

    Beyond widely used pairwise two-body potentials,multi-body potentials, especially for three- and fourbody potentials, have been proposed for capturing higherorder interactions.[63,75–79,95–97]In addition, several geometrical techniques, such as Voronoi diagram, Delaunay triangulation,[98]and alpha shape[76,96]were employed to define the neighboring bodies and to identify the multi-body interactions.

    3.4.1. Three-body potentials

    Three-body contacts have been introduced in several studies where physical models were developed with explicitly incorporating a three-body interaction.[95,96,99–101]Based on Delaunay triangulation, Li and Liang used the alpha shape to define three-body interactions for protein structures.[96]Here,we briefly introduce this approach as an example of this class of statistical potential. In the model of Li and Liang, volume overlap was used to identify three-body contacts, and such volume overlap occurs if three atoms from nonbonded residues share a Voronoi vertex.[96]Furthermore,Li and Liang introduced a nonadditive coefficient by measuring the deviation of three-body interactions from three independent pairwise interactions.[96]The nonadditive coefficient v of threebody interactions was defined as[96]

    where P(i,j,k) is the three-body propensity for residues of types i, j, and k in native structure database, and P′(i,j,k)is the propensity if the three-body interaction is the simple consequence of three independent pairwise contact interactions. The expressions of P(i,j,k) and P′(i,j,k) are written as follows:[96]

    Here, P(i,j) is the pairwise contact propensity for residue types i and j, and P(i,j) equals to the odds ratio of the observed probability q(i,j) for two-body contact of residues i and j to the excepting probability p(i,j)for two-body contact of residues i and j. The observed probability q(i,j,k)and expecting p(i,j,k)are for three-body contact,respectively. The details of calculating p(i,j), q(i,j), q(i,j,k), and q(i,j,k)can be found in Ref.[96].Finally,the potential energy ?E(M)for a molecule M modeled as a set of fused hard spheres can be expressed as[96]

    where σij∈κ represents pairwise alpha contacts,and σijk∈κ represents three-body alpha contacts.

    3.4.2. Four-body potentials

    To capture more complicated interactions, several fourbody potentials have been developed,[63,75,77–79,95,97]and Delaunay tessellation algorithms were widely used for finding the nearest neighbors in native structures. Krishnamoorthy and Tropsha have verified that their four-body potentials derived with Delaunay tessellation performed better than pairwise two-body potentials in distinguishing correct sequences/structures and generating Z-scores.[77]Later,considering more detailed interactions between backbones and side chains and including some of the sequential information of proteins,Feng et al. have proposed a new four-body statistical potential based on Delaunay tessellation.[78]Instead of using Delaunay tessellation, Gniewek et al. have used a simple geometric construction to develop a four-body contact potential,and they found that the performance of their optimized coarsegrained contact potentials was comparable to the performance of DFIRE for all proteins larger than 80 amino acids.[79]

    As shown above, due to the huge computational complexity and limited structure data, most existing multi-body potentials are based on local residue contact rather than distance and with the growth of protein structures deposited in database, multi-body distance-dependent statistical potentials can be developed for capturing higher order interactions in protein structures.

    3.5. Involvement of orientation-dependent interaction

    Atoms in proteins generally carry charges and thereby there can exist orientation-dependent correlations between polar atom blocks due to dipole–dipole-like interaction and associated polarization effect.[85]Consequently, orientationdependent correlation between atom blocks has been considered as a complementary contribution to above described statistical potentials. Recently, Wang and Huang have shown that the involvement of orientation contribution can visibly improve the performance of the statistical potential.[86]In geometry, the relative position and orientation between two 3D objects need to be determined by six degrees of freedom, including three translations and three rotations, or one distance and five angular variables.[83]Hence,the definition and calculation approaches for the orientation between two blocks are more diverse than those for the distance of two blocks, and a more adequate native structure database is needed for developing orientation-dependent statistical potentials.

    3.5.1. Definition of relative orientation

    Ideally,for proteins,an amino acids residue can be represented by several blocks,and each block is centered by each of heavy atoms. As shown in Fig.1(a),the relative orientation of a block pair is described by five angles: θi,?i,θj,?j,and ωij.The first four angles are the polar coordinates of the rijvector in the local 3D reference frame of each amino acid and ωijdescribes the relative rotation of Vzvectors of two amino acids along the rijaxis. The local 3D reference frame is defined as

    where r12=r(i1)?r(i)and r13=r(i2)?r(i)are the relative vectors from atom i to atoms i1and i2,respectively;Vzis the normal vector to the plane within which Vxand Vzlie.

    Fig.1. (a)The complete definition of relative position and orientation between two blocks centered by atoms i and j, respectively, in which one scalar distance rij and five scalar angles(θi,Ψi,θj,Ψj,ωij)parameters are used. (b) The definition of relative position and orientation between two blocks i and j for OPUS-PSP,in which two direction vectors ri j and rji and one scalar inter-rotation angle ωij parameters are used. (c) The definition of relative position and orientation between two polar atoms i and j for dDFIRE, in which one scalar distance rij and three scalar angles(θi,θj,θi j)parameters are used. (d)The definition of relative position and orientation between two atoms i and j for ORDER_AVE,in which one scalar distance rij and three scalar angles(θi,θi,ωij)parameters are used.

    Afterwards,by taking the relative position and orientation as geometrical parameters,equation(7)can be rewritten as

    3.5.2. Recent developed orientation-dependent potentials

    A perfect definition of relative orientation for proteins is so detailed that building an orientation-dependent potential would experience the challenge of sparse experimental data and long computational quantity. The existing strategies used to overcome the problem mainly include: (i) reducing representations of proteins; (ii) reducing variables of defining the orientation;and(iii)reducing dimensions of the joint probability distribution. In this subsection, several recent orientationdependent statistical potentials are briefly introduced as follow.

    OPUS-PSPOPUS-PSP potential,proposed by Lu et al.,defines orientations for blocks of atoms bonded rigidly within the same residue.[80]In OPUS-PSP potential, 20 amino acid residues are decomposed into 19 rigid-body block types, in which all heavy atoms are assumed to be in the same plane,and each residue contains several blocks. As shown in Fig.1(b),the relative orientation of blocks i and j respectively with the origins oiand ojis determined by two relative direction vectors rijand rji, and an inter-rotation angle ωij. For the differences among these 19 rigid-body blocks, the definition of local 3D reference frame for each block type and interrotation angle is specified,which can be found in Lu et al.[80]Finally, the orientation-dependent energy function can be expressed by[80]

    In dDFIRE,the finite-ideal-gas reference state was used as the reference state.

    ORDER_AVEORDERAVE is a four-dimensional joint probability distribution with five atoms per residue.[102]As shown in Fig.1(d), for atom pair i and j, two heavy atoms ir and jr that directly connect to i and j, were selected, and the orientations of atoms i and j can be described by four variables, rij, α, β, and γ. Here, rijis the distance between atoms of types i and j; α and β are the scalar angles of the atom groups ir–i–j and i–j–jr,respectively;and γ is the dihedral angle of atom groups ir–i–j–jr. Finally, the orientationdependent energy ?Eijfor the atom pair i and j can be expressed as[102]

    KORPKORP, a side-chain independent potential, was defined by a 6D joint probability that only depends on the relative orientation and position of three backbone atoms.[83]As shown in Fig.1(a),six parameters including a distance rijand five angles (θi,?i,θj,?j,ωij) are used. However, it is quite challenging to obtain a detailed 6D joint probability based on the relatively small number of experimental structures. Thus,a constant z was added to prevent infinite values for very small probabilities and to improve the numerical stability for lowcount statistics,and the orientation-dependent potential can be given by[83]

    The above described involvements of orientationdependent potentials can provide a very valuable contribution to a statistical potential,[80,83,85,102]while due to sparse experimental data and long computational cost, the quality of an orientation-dependent potential may strongly rely on the definition of relative orientation for atom groups in proteins.

    4. Statistical potentials for RNA 3D structures

    In the recent decade,as the number of experimental RNA 3D structures deposited in PDB database[8]increases,several statistical potentials have been derived based on known structures for RNA structure prediction and assessment,[37,49–56]partially with the aid of the experience for proteins. However,the existing statistical potentials for RNAs are significantly less than those for proteins and still do not reach a satisfactory performance. In the following, we will give an overview on the recently developed statistical potentials for RNA 3D structure evaluation and will emphasize the possible development on statistical potentials for RNAs in the near future.

    4.1. Existing statistical potentials for RNAs

    To our knowledge, the existing statistical potentials for RNA 3D structure evaluation[50–56]were developed(partially)with the aid of the experience or those in protein structure prediction and assessment. The existing statistical potentials for RNAs have been organized in Table 2,and will be introduced explicitly in the following.

    4.1.1. Distance- and torsion angle-dependent statistical potentials

    RASPBased on the averaging reference state, Capriotti et al. have built RASP at both coarse-grained and all-atom(23 clustered atom types) levels from a non-redundant training set of 85 RNA structures,and distance between atom pairs was considered as the geometrical parameter.[51]Capriotti et al. have shown that RASP has a better performance than NAST[37]which is a nucleotide-level coarse-grained statistical potential composed of bond,angle,dihedral,and non-bond terms.

    KB potentialUsing the quasi-chemical approximation reference state, Bernauer et al. have derived distancedependent and fully differentiable statistical potentials at both coarse-grained and all-atom (85 atom types) levels from a training set of 77 native RNAs.[50]Owing to the implementation of Dirichlet process mixture models, the statistical potentials are fully differentiable,which makes them applicable for molecular dynamics simulations.

    3dRNAscoreDifferent from RASP and KB potentials, 3dRNAscore is composed of the distance- and torsion angle-dependent potentials based on the averaging reference state.[52]In 3dRNAscore, 85 atom types and 7 torsion angle types were defined, and a set of 317 RNA native structures was used for training the potential.[52]In addition, to appropriately consider the contributions of the two kinds of energy terms, a weight factor w has been optimized by four typical RNA decoys.[52]The test of 3dRNAscore has shown that, the dihedral-dependent potential alone performs visibly worse than the distance-dependent potential averagely,while the combined potential with dihedral-dependent and distance-dependent contributions generally performs better than the distance-dependent potential.[52]Thereby, the dihedral-dependent potential has been verified useful to improve the accuracy of a distance-dependent potential.

    According to the results reported by Wang et al.,[52]overall,the performance of RASP is worse than that of KB potential and 3dRNAscore, and 3dRNAscore has the best performance among these three statistical potentials.[52]The mainly reasons might be that the KB potential of all-atom and 3dRNAscore are based on 85 atom types,while only 23 clustered atom types are involved in the all-atom version of RASP. In addition, the best performance of 3dRNAscore might be attributed to the explicit emphasis on the local torsional structure feature of RNA backbones which can be important for RNA 3D structures.

    4.1.2. Comparisons of different reference states for RNA 3D structure evaluation

    For RNAs, only the averaging and quasi-chemical approximation reference states were used to build statistical potentials until our very recent study.[55]To understand the existing widely used six reference states for RNAs, we performed a comprehensive examination on these six reference states,including the averaging,quasi-chemical approximation,atom-shuffled, finite-ideal-gas, spherical-non-interacting, and random-walk-chain reference states; see also Subsection 3.1.In the work,[55]using the same training set and parameters,e.g., bin width and distance cutoff, we constructed six allatom(85 atom types)distance-dependent statistical potentials corresponding to the six reference states,and examined these potentials against three RNA test sets;see Table 2. Our extensive examinations showed that, overall, for identifying native structures and ranking decoy structures,the finite-ideal-gas[66]and random-walk-chain[69]reference states were slightly better than others, while for identifying near-native structures,there was slight difference among the six reference states.Moreover, compared among RASP, KB potential, and 3dRNAscore,3dRNAscore has the similar performance to the top ones of these six statistical potentials, while RASP and KB potential have visibly lower performance than the others.[55]

    Table 2. The existing statistical potentials for RNA structure evaluation.

    4.1.3. Four-body contact potential

    Masso has developed a four-body contact potential of RAMP, which is the first multibody statistical potential for RNAs.[56]In RAMP, atomic four-body nearest-neighbors were generated by the Delaunary tessellation[98]for each structure in a training set of 85 native RNAs, and each RNA was represented by four atom types of C,N,O,and P.[56]Thus,35 distinct quadruplet types can be produced by 4-letter atomic alphabet. With atomic quadruplet contacts as geometric parameter and the inversed Boltzmann’s law,equation(7)can be rewritten as[56]

    4.1.4. Utilization of 3D convolutional neural network

    Completely diverged from the above-described traditional statistical potentials, Li et al. have employed a 3D convolutional neural network to develop two 3DCNNbased scoring models, named RNA3DCNN_MD and RNA3DCNN_MDMC, for assessing near-native RNA decoys and RNA decoys with large root-mean-square-deviation(RMSD) fluctuation, respectively.[53]Two training sets were generated by molecular dynamics (MD) simulations and Monte Carlo (MC) structure prediction for each of 414 RNAs. RNA3DCNN_MD was trained by the frist set, while RNA3DCNN_MDMC was trained by the training set combined the two training sets. The test showed that, for decoys with RMSD less than 1.0 ?A,the performance of RNA3DCNN was similar to or worse than that of traditional statistical potentials, while for decoys generated by RNA structure predictive models, RNA3DCNN performed obviously superior to other existing statistical potentials in identifying native structures.[53]

    4.2. Limitation of current statistical potentials and perspective for RNAs

    According to Eq. (1), a set of pairwise potentials can be extracted from the distance distribution of atom pairs by the inverse of Boltzmann’s law. However, an individual pairwise potential involves not only the interaction between two atoms but also the influence of the surrounding medium on the interaction.[60]Therefore, reference states and other approaches have been developed aiming to eliminate redundant information.Based on the experiences in protein structure prediction and assessment from which the existing RNA statistical potentials benefit,we will make a brief discussion for developing new statistical potentials for RNA structure evaluation.

    4.2.1. On reference states

    The existing six reference states can be classified into two types: based on experimental RNA 3D structures[43,65,68]and based on physical modeling.[66,67,69]The physical-modelbased reference states are built by various physical models,and they can better describe the conformation space of RNA decoys than the 3D-structure-based reference states.However,the physical-model-based reference states are still deviate significantly from the ideal one due to the ignorance of detailed structure features,and our recent study shows that the performance of statistical potentials based on physical-model-based reference states is still not at a high level for realistic decoys from RNA structure prediction models.[55]For example, the ideal-gas reference state completely neglects the connectivity of RNA chains,and a random-walk reference state ignores the intrinsic rigidity of RNA chains. Moreover,both of the reference states neglect side chains and atom types. Therefore, it is still necessary to propose more realistic reference states for RNAs or circumvent the problem of reference states.

    One possible way to simulate more realistic RNA chains is to involve the chain rigidity of RNAs, e.g., an RNA chain is modeled as a semiflexible polymer rather than a completely random polymer.[103,104]Another possible way is to develop a combined reference state composed of the physical-modelbased and structure-based reference states which may complement each other. A third possible way is to develop an ensemble-based reference state of ideal-gas-like states of different dimension parameters or of random-walk-chain states of difference Kuhn lengths, since backbone atoms and side chain ones may behave very differently. Due to the difficulty in modeling an ideal reference state, an alternative way is to circumvent the problem. Huang and Zou have proposed an iterative method to develop scoring functions for protein–ligand and protein–protein docking.[61,62]This method provides a possible way to circumvent reference states for modeling statistical potentials for RNAs, although the performance may rely strongly on the high-quality decoy training set which should be generated by various structure prediction models.

    4.2.2. On short-,medium-,and long-range interactions

    Tanaka and Scheraga[39]have classified the inter-residue interactions into short-, medium-, and long-ranged ones for proteins, and afterwards, the preference of amino acids to form medium-and long-range contacts between has been analyzed for proteins.[72–74,89]Until now,most statistical potentials for RNAs do not distinguish the contributions of short-,medium-, and long-ranged interactions. However, these different residue-separation-ranged interactions should play different roles in RNA folding process or in stabilizing folded RNA structures. For example,the short-,medium-,and longranged potentials may correspond to the interactions within loops, those stabilizing secondary segments, and those stabilizing tertiary contacts for RNAs. Although the local interactions were split from non-local interactions in RASP by Capriotti et al.,[51]the local potential was built at each nucleotide separation which is too subtle to avoid sparse data problem and might not be able to represent specific topological features. Therefore, the classification of short-, medium-, and long-ranged interactions is in accordance with the folding and structure features of RNAs and may be valuable for developing statistical potentials of higher performance for RNA structure evaluation.

    4.2.3. On geometrical parameters

    Compared with pairwise two-body potential, multi-body potential could capture more complex 3D interactions.[78,96]Up to now,RAMP is the only RNA four-body statistical potential,but its overall performance is not better than that of RASPALL,[56]which may be attributed to that RAMP only involves four atom types to represent RNA structures and RAMP is a (spatially) local contact potential ignoring spatially nonlocal interactions. Therefore, for building a multi-body potential, a proper geometrical parameter may be very important which should not only describe the key geometrical relation between multi-body atoms, but also be conveniently used to build a multi-body potential based on limited native structure database.Another key problem in modeling multi-body potential may be the coarse-graining of RNA atoms which should not only retain the major structure feature of RNAs, but also effectively decrease the complexity of multi-body potentials and the demand for structure data in training set.

    On the other hand, most existing RNA statistical potentials are distance-dependent, expect for 3dRNAscore which is composed of the distance- and the local torsional angledependent terms.[52]The better performance of 3dRNAscore than prior traditional statistical potentials is mainly attributed to the explicit emphasis on the local torsional structure feature of RNA backbone. Hence, other geometrical parameters can be involved into distance-dependent statistical potentials to improve the performance, such as orientation, which has been highlighted important for proteins;[86]see Subsection 3.5. Another geometrical parameter can be the nonlocal torsional angle beyond one nucleotide, which may also improve the description for RNA backbone and be helpful for a distance-dependent potential. Owing to current limited experimental RNAs, it is very important to define the orientation variables,the importance of which has been discussed above.Based on current limited experimental RNAs,it is very important to define the orientation variables and deal with the dependence among the variables. Compared to 20 amino acids with 167 atoms for proteins,RNAs has 4 nucleotides with 85 atom types, which may make it easier to build high-performance orientation-dependent potentials for RNAs.

    4.2.4. On unique feature of RNAs

    From the above, the existing statistical potentials for RNAs were developed mainly according to those for proteins.There are some important similar structure features between RNAs and proteins. For examples, first, RNAs and proteins are both polymers with backbones and side chains. Second,there are distinctively ordered secondary segments in RNA and protein structures, such as helices and loops for RNAs,and α-helices,β-sheets,and loops for proteins;third,the tertiary structures of RNA and proteins are both very compact and approximately sphere-like. Such similarity between RNA and proteins structures has enabled the great benefit of the existing statistical potentials for RNAs from those for proteins.

    However,there are also several distinctive differences between RNAs and proteins. The major forces for stabilizing structures are very different for RNAs and proteins. RNAs are generally stabilized by strong secondary interactions of base-pairing/stacking and moderate tertiary interactions,while proteins are generally stabilized by hydrophilic and hydrophobic interactions. Consequently, the secondary and tertiary interactions for RNAs may be roughly classified while those for proteins are mutually coupled. The separate involvement of secondary and tertiary contributions in a statistical potential may improve its performance. Furthermore, basestacking can promote the helical orientation of RNAs and consequently RNA structures are apparently less spherical than protein structures,[105,106]which may suggest that angular or orientation-dependent contribution is more important for a statistical potential for RNAs.

    Very importantly, RNAs are strongly charged polyanion,[106–109]which is distinctively different from proteins generally with certain positively and negatively charged domains. Generally, for RNAs, metal ions are essentially required to neutralize the negative charges on RNA backbones,and consequently would favor the folding of compact native structures. Furthermore, cations with different valences and concentrations could have very different influences on the deformation of folded RNA structures.[108,109,111–113]For example, due to polyanionic nature, RNA structures can become less compact at lower salt and an RNA can have visibly more compact structure at multivalent salt than at monovalent salt.[32,33,35]Therefore,the involvement of ion-electrostatic interaction in a statistical potential may be helpful for improving its performance for RNA 3D structure evaluation, especially for the NMR structures.[55]Finally, machine and deep learning methods which could dig critical and uneasy to be detected RNA feature information are also able to be applied to build statistical potentials.[53]

    5. Conclusions

    Since various computational RNA 3D structure prediction models have appeared in recent years, scoring functions are essentially needed for evaluating the predicted structure candidates. Because statistical potentials for protein 3D structure prediction and evaluation have been developed from a long time ago,there are lots of experience which can be very helpful for RNAs. In this review, we firstly introduced the fundamental principles of deriving the formulas of statistical potentials. Afterwards, we gave a brief view on the development of statistical potentials for proteins, and we emphasize the choice of geometrical parameters and reference states which are crucial for building different statistical potentials.Furthermore, we introduced the recent advances in statistical potentials for RNAs. Finally, we make a detailed discussion on possible further improvements of statistical potentials for RNAs,which may be helpful for developing high-performance statistical potentials for RNA 3D structure evaluation.

    Acknowledgements

    We are grateful to Professors Shi-Jie Chen(University of Missouri), Xiaoqin Zou(University of Missouri), Jian Zhang(Nanjing University), Yi Xiao (Huazhong University of Science and Technology),and Sheng-You Huang(Huazhong University of Science and Technology)for valuable discussions.

    免费看不卡的av| 午夜精品国产一区二区电影| av网站在线播放免费| 亚洲第一区二区三区不卡| 国产爽快片一区二区三区| 一区二区三区激情视频| 国产免费视频播放在线视频| 亚洲国产欧美在线一区| av在线观看视频网站免费| av网站免费在线观看视频| 99国产综合亚洲精品| 国产毛片在线视频| 欧美av亚洲av综合av国产av | 丰满少妇做爰视频| 人人妻人人澡人人看| 香蕉国产在线看| 久久99精品国语久久久| 欧美精品人与动牲交sv欧美| 午夜福利在线观看免费完整高清在| 久久国产精品大桥未久av| 久久精品国产亚洲av涩爱| 亚洲男人天堂网一区| 精品人妻在线不人妻| 一区二区日韩欧美中文字幕| 亚洲国产看品久久| 欧美成人午夜免费资源| 国产精品国产三级国产专区5o| 中文乱码字字幕精品一区二区三区| 最近手机中文字幕大全| 国产免费一区二区三区四区乱码| 精品少妇内射三级| 国产精品女同一区二区软件| 男女边吃奶边做爰视频| 一本—道久久a久久精品蜜桃钙片| 如日韩欧美国产精品一区二区三区| 色网站视频免费| 久久鲁丝午夜福利片| 欧美激情 高清一区二区三区| 波多野结衣av一区二区av| av在线app专区| 18禁观看日本| 亚洲图色成人| 久久人人97超碰香蕉20202| 在线 av 中文字幕| 国产黄色视频一区二区在线观看| 久久99热这里只频精品6学生| 亚洲欧洲日产国产| 大码成人一级视频| 免费播放大片免费观看视频在线观看| 婷婷色综合大香蕉| 久久免费观看电影| 只有这里有精品99| 亚洲av欧美aⅴ国产| 午夜福利在线免费观看网站| 久久久久久久久久人人人人人人| 国产精品一区二区在线观看99| 成人免费观看视频高清| 久久精品人人爽人人爽视色| 中文字幕人妻丝袜一区二区 | 国产淫语在线视频| 免费少妇av软件| 赤兔流量卡办理| 亚洲国产精品999| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 亚洲成色77777| 国产精品.久久久| 日韩,欧美,国产一区二区三区| 不卡视频在线观看欧美| 国产欧美亚洲国产| 久久久久久久久久人人人人人人| 亚洲人成77777在线视频| 看免费成人av毛片| 亚洲经典国产精华液单| 熟女av电影| 如何舔出高潮| 少妇人妻精品综合一区二区| 亚洲精品国产色婷婷电影| 亚洲情色 制服丝袜| 久久97久久精品| 大片免费播放器 马上看| 丝袜美足系列| 国产片特级美女逼逼视频| 亚洲国产精品成人久久小说| 下体分泌物呈黄色| 丰满饥渴人妻一区二区三| 人妻系列 视频| 18禁动态无遮挡网站| 久久久久久久国产电影| 免费黄网站久久成人精品| 国产亚洲精品第一综合不卡| 亚洲欧美色中文字幕在线| 亚洲久久久国产精品| 波野结衣二区三区在线| 日产精品乱码卡一卡2卡三| 亚洲欧美成人综合另类久久久| 久久影院123| 精品少妇内射三级| 夜夜骑夜夜射夜夜干| 亚洲国产色片| 国产黄色视频一区二区在线观看| 国产乱来视频区| 麻豆乱淫一区二区| 女人被躁到高潮嗷嗷叫费观| 老熟女久久久| 精品国产一区二区三区四区第35| 18禁国产床啪视频网站| 免费黄网站久久成人精品| 亚洲国产成人一精品久久久| 涩涩av久久男人的天堂| 午夜免费男女啪啪视频观看| 91精品三级在线观看| 免费高清在线观看日韩| 天天躁夜夜躁狠狠久久av| 超色免费av| 十分钟在线观看高清视频www| 欧美国产精品一级二级三级| 国产免费视频播放在线视频| 一区二区三区乱码不卡18| 亚洲欧美中文字幕日韩二区| 久久精品久久久久久久性| 伊人久久国产一区二区| 少妇的逼水好多| 免费在线观看完整版高清| 性高湖久久久久久久久免费观看| 国产有黄有色有爽视频| 男女边摸边吃奶| 国产精品无大码| 国产97色在线日韩免费| 成人漫画全彩无遮挡| 久久97久久精品| 色婷婷久久久亚洲欧美| 激情视频va一区二区三区| 久久久久久人妻| 亚洲三级黄色毛片| 久久av网站| 成人免费观看视频高清| 电影成人av| 亚洲伊人久久精品综合| xxx大片免费视频| 久久久精品94久久精品| 天堂8中文在线网| 日韩免费高清中文字幕av| 97精品久久久久久久久久精品| 最近的中文字幕免费完整| 一区二区三区激情视频| 男人舔女人的私密视频| 18+在线观看网站| 美女脱内裤让男人舔精品视频| 免费高清在线观看日韩| 满18在线观看网站| 日韩中字成人| 熟女少妇亚洲综合色aaa.| 最近中文字幕高清免费大全6| 最近手机中文字幕大全| 水蜜桃什么品种好| 五月天丁香电影| 国产精品久久久久久精品古装| 日本91视频免费播放| 少妇的逼水好多| 99精国产麻豆久久婷婷| 视频在线观看一区二区三区| 欧美日韩成人在线一区二区| 18禁裸乳无遮挡动漫免费视频| 搡老乐熟女国产| 99久久精品国产国产毛片| 欧美日韩亚洲高清精品| 美女午夜性视频免费| 久久精品国产亚洲av涩爱| 免费女性裸体啪啪无遮挡网站| 国产黄频视频在线观看| 中国三级夫妇交换| 搡老乐熟女国产| 午夜福利影视在线免费观看| 另类精品久久| 精品视频人人做人人爽| 久久久久国产网址| 高清黄色对白视频在线免费看| 国产精品 国内视频| 欧美精品一区二区大全| 视频在线观看一区二区三区| 国产成人精品久久二区二区91 | 国产av一区二区精品久久| 亚洲久久久国产精品| 在线天堂最新版资源| 青草久久国产| 中文精品一卡2卡3卡4更新| av国产精品久久久久影院| 999精品在线视频| 国产精品久久久av美女十八| 亚洲精品美女久久av网站| 久久人妻熟女aⅴ| 日本wwww免费看| 久久热在线av| 一本大道久久a久久精品| 纯流量卡能插随身wifi吗| 午夜福利一区二区在线看| 精品国产露脸久久av麻豆| 亚洲一区二区三区欧美精品| 十分钟在线观看高清视频www| 91精品国产国语对白视频| 九色亚洲精品在线播放| www.自偷自拍.com| 国产成人一区二区在线| 久久99一区二区三区| 日日撸夜夜添| 久久综合国产亚洲精品| 欧美日韩精品网址| 久久久久久人妻| 亚洲第一青青草原| av网站免费在线观看视频| a级毛片在线看网站| 久久鲁丝午夜福利片| 不卡av一区二区三区| 国产精品国产av在线观看| 国产熟女午夜一区二区三区| 性色av一级| 亚洲av福利一区| 久久97久久精品| 天堂俺去俺来也www色官网| 国产av一区二区精品久久| 九草在线视频观看| 久久99精品国语久久久| 一本—道久久a久久精品蜜桃钙片| 一区二区三区四区激情视频| 亚洲一区中文字幕在线| 丝瓜视频免费看黄片| 人妻一区二区av| 99国产综合亚洲精品| 99香蕉大伊视频| 久久久精品94久久精品| 99久国产av精品国产电影| 免费黄网站久久成人精品| 欧美老熟妇乱子伦牲交| 天天操日日干夜夜撸| 桃花免费在线播放| 看免费av毛片| 成人二区视频| 亚洲第一青青草原| 日韩人妻精品一区2区三区| 爱豆传媒免费全集在线观看| 毛片一级片免费看久久久久| 亚洲精品成人av观看孕妇| 男女午夜视频在线观看| 日韩,欧美,国产一区二区三区| 视频在线观看一区二区三区| 大码成人一级视频| 男男h啪啪无遮挡| 亚洲精品在线美女| 久久av网站| 这个男人来自地球电影免费观看 | 日日撸夜夜添| 午夜91福利影院| 国产亚洲av片在线观看秒播厂| 日本免费在线观看一区| 国产成人免费观看mmmm| 2021少妇久久久久久久久久久| 久久久久久久大尺度免费视频| 国产精品不卡视频一区二区| 亚洲av日韩在线播放| 视频区图区小说| 免费观看无遮挡的男女| 蜜桃国产av成人99| 日韩精品免费视频一区二区三区| 最近中文字幕2019免费版| 黄片无遮挡物在线观看| 国产高清国产精品国产三级| 一级a爱视频在线免费观看| 精品视频人人做人人爽| 国产成人精品久久久久久| 亚洲av男天堂| 色播在线永久视频| 亚洲精品在线美女| 久久毛片免费看一区二区三区| 久久免费观看电影| 日韩三级伦理在线观看| 久久精品国产综合久久久| 婷婷成人精品国产| 亚洲成av片中文字幕在线观看 | 中国国产av一级| 亚洲欧美一区二区三区国产| 欧美97在线视频| 精品国产乱码久久久久久男人| 性少妇av在线| 王馨瑶露胸无遮挡在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产欧美亚洲国产| 观看av在线不卡| 国产精品女同一区二区软件| 国产精品国产三级国产专区5o| 1024香蕉在线观看| 2022亚洲国产成人精品| 国产片内射在线| 最近2019中文字幕mv第一页| 亚洲第一av免费看| 你懂的网址亚洲精品在线观看| 桃花免费在线播放| av国产久精品久网站免费入址| av线在线观看网站| 久久久久久久久久久免费av| 精品国产一区二区久久| 9191精品国产免费久久| 天天躁夜夜躁狠狠躁躁| 中文字幕人妻丝袜一区二区 | 欧美日韩成人在线一区二区| 欧美精品av麻豆av| 侵犯人妻中文字幕一二三四区| 久久婷婷青草| 国产免费福利视频在线观看| 色婷婷av一区二区三区视频| 国产免费一区二区三区四区乱码| 丁香六月天网| 亚洲一码二码三码区别大吗| xxx大片免费视频| 久久久久网色| h视频一区二区三区| 成年美女黄网站色视频大全免费| 熟女少妇亚洲综合色aaa.| 国产乱来视频区| 精品人妻熟女毛片av久久网站| 亚洲美女黄色视频免费看| 在线观看国产h片| 麻豆精品久久久久久蜜桃| 热99久久久久精品小说推荐| 亚洲一区二区三区欧美精品| 久久久久人妻精品一区果冻| 美女视频免费永久观看网站| 亚洲成色77777| 可以免费在线观看a视频的电影网站 | 美女中出高潮动态图| 欧美最新免费一区二区三区| 亚洲美女黄色视频免费看| 国产 精品1| 水蜜桃什么品种好| 久久亚洲国产成人精品v| 大香蕉久久网| 99久久中文字幕三级久久日本| 中文乱码字字幕精品一区二区三区| 久久av网站| 亚洲欧洲国产日韩| 卡戴珊不雅视频在线播放| 男女午夜视频在线观看| 边亲边吃奶的免费视频| 美女xxoo啪啪120秒动态图| 丝袜在线中文字幕| 老司机影院毛片| 777米奇影视久久| 午夜激情av网站| h视频一区二区三区| 欧美日韩成人在线一区二区| 午夜福利视频精品| 国产一区有黄有色的免费视频| 久久精品国产亚洲av高清一级| 亚洲精品国产一区二区精华液| 精品国产超薄肉色丝袜足j| 大香蕉久久网| 又粗又硬又长又爽又黄的视频| 麻豆乱淫一区二区| 亚洲av电影在线进入| xxx大片免费视频| 热re99久久精品国产66热6| 看免费成人av毛片| 99热网站在线观看| 国产成人aa在线观看| 亚洲精品美女久久久久99蜜臀 | 久久午夜福利片| 亚洲三级黄色毛片| 大香蕉久久成人网| 午夜福利影视在线免费观看| 一区福利在线观看| 免费观看a级毛片全部| 久久精品国产综合久久久| 青春草视频在线免费观看| 香蕉国产在线看| 久久午夜福利片| 成年女人在线观看亚洲视频| 美女主播在线视频| 日本av手机在线免费观看| videos熟女内射| 91午夜精品亚洲一区二区三区| 亚洲av电影在线观看一区二区三区| 欧美成人午夜免费资源| 亚洲av日韩在线播放| 观看美女的网站| 日韩视频在线欧美| 街头女战士在线观看网站| 国产精品久久久久久av不卡| 91成人精品电影| 日韩一本色道免费dvd| 日本爱情动作片www.在线观看| 亚洲三级黄色毛片| 哪个播放器可以免费观看大片| 国语对白做爰xxxⅹ性视频网站| 精品人妻在线不人妻| 日本av免费视频播放| 欧美人与性动交α欧美软件| 考比视频在线观看| 日韩制服骚丝袜av| 2021少妇久久久久久久久久久| 久久av网站| 国产精品人妻久久久影院| 性色av一级| 国产av码专区亚洲av| 国产高清国产精品国产三级| 亚洲图色成人| 1024视频免费在线观看| 欧美日韩国产mv在线观看视频| 精品人妻在线不人妻| 日韩制服骚丝袜av| 久久人人爽人人片av| 伦理电影大哥的女人| 免费av中文字幕在线| 欧美日韩一级在线毛片| 久久久久久久精品精品| 日本欧美视频一区| 看非洲黑人一级黄片| 亚洲一区中文字幕在线| 免费大片黄手机在线观看| 午夜激情久久久久久久| 激情视频va一区二区三区| 中文字幕人妻丝袜制服| 亚洲国产精品一区三区| 乱人伦中国视频| 亚洲精品自拍成人| av卡一久久| 国产精品国产av在线观看| 七月丁香在线播放| 国产 精品1| 欧美日韩综合久久久久久| 激情五月婷婷亚洲| 久久久a久久爽久久v久久| 香蕉国产在线看| 国产精品熟女久久久久浪| 国产一区二区三区综合在线观看| 精品午夜福利在线看| 熟妇人妻不卡中文字幕| 午夜日本视频在线| 嫩草影院入口| 久久久欧美国产精品| 国产97色在线日韩免费| 韩国av在线不卡| 欧美xxⅹ黑人| 人人澡人人妻人| 欧美av亚洲av综合av国产av | 黄片无遮挡物在线观看| 久久国内精品自在自线图片| 99re6热这里在线精品视频| 男女边吃奶边做爰视频| 精品亚洲成国产av| 久久精品久久久久久噜噜老黄| 在线观看美女被高潮喷水网站| 少妇猛男粗大的猛烈进出视频| 精品人妻在线不人妻| 久久久久久伊人网av| 一区二区三区激情视频| 国产成人精品久久二区二区91 | 欧美精品一区二区免费开放| 不卡av一区二区三区| 美女高潮到喷水免费观看| 久久ye,这里只有精品| 久久久a久久爽久久v久久| 日韩伦理黄色片| 日韩 亚洲 欧美在线| 人人妻人人澡人人爽人人夜夜| 天堂8中文在线网| 日韩一区二区视频免费看| 一本色道久久久久久精品综合| 晚上一个人看的免费电影| 丝袜人妻中文字幕| 在线 av 中文字幕| 亚洲成色77777| 天堂8中文在线网| 日韩电影二区| 成人国产av品久久久| 这个男人来自地球电影免费观看 | 春色校园在线视频观看| 不卡视频在线观看欧美| 多毛熟女@视频| 久久精品国产a三级三级三级| 日本色播在线视频| 97人妻天天添夜夜摸| 叶爱在线成人免费视频播放| 校园人妻丝袜中文字幕| 久热这里只有精品99| 日韩中文字幕视频在线看片| 999久久久国产精品视频| 国产一区二区三区av在线| 在线免费观看不下载黄p国产| 狠狠婷婷综合久久久久久88av| 亚洲色图综合在线观看| 午夜福利,免费看| www.精华液| 高清在线视频一区二区三区| 国产极品粉嫩免费观看在线| 婷婷成人精品国产| 国产极品天堂在线| 亚洲欧美中文字幕日韩二区| 久久人妻熟女aⅴ| 精品卡一卡二卡四卡免费| 精品人妻熟女毛片av久久网站| 亚洲色图 男人天堂 中文字幕| 我要看黄色一级片免费的| 波野结衣二区三区在线| 一级片'在线观看视频| 女人被躁到高潮嗷嗷叫费观| 久久精品亚洲av国产电影网| 亚洲精品第二区| 中文字幕另类日韩欧美亚洲嫩草| 9191精品国产免费久久| 女性生殖器流出的白浆| 久久青草综合色| 欧美人与善性xxx| 制服人妻中文乱码| 久久久欧美国产精品| 国产精品免费大片| 日韩视频在线欧美| 亚洲精品乱久久久久久| 久久精品国产a三级三级三级| 免费播放大片免费观看视频在线观看| 亚洲成人一二三区av| 国产一区二区三区av在线| 99热全是精品| 最新的欧美精品一区二区| 国产日韩欧美视频二区| 亚洲精品日韩在线中文字幕| 丝袜脚勾引网站| 国产成人精品无人区| 亚洲av免费高清在线观看| 国产成人精品在线电影| av在线播放精品| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 久久精品熟女亚洲av麻豆精品| 一本色道久久久久久精品综合| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 日本91视频免费播放| 十分钟在线观看高清视频www| 老熟女久久久| av片东京热男人的天堂| 乱人伦中国视频| 久久久久久伊人网av| 国产在线一区二区三区精| 亚洲国产av新网站| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 日本91视频免费播放| 成人国产av品久久久| 精品一区二区三区四区五区乱码 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久国内精品自在自线图片| 日本免费在线观看一区| 99久国产av精品国产电影| 制服丝袜香蕉在线| 精品久久久精品久久久| 在线天堂最新版资源| 2021少妇久久久久久久久久久| 日本-黄色视频高清免费观看| 国产在视频线精品| 国产av码专区亚洲av| 亚洲欧美清纯卡通| 狠狠婷婷综合久久久久久88av| 中文字幕精品免费在线观看视频| 国产有黄有色有爽视频| 一级毛片我不卡| 日本-黄色视频高清免费观看| av不卡在线播放| 亚洲精品aⅴ在线观看| 亚洲欧美中文字幕日韩二区| 国产精品 国内视频| 国产 精品1| 亚洲人成电影观看| 黄色视频在线播放观看不卡| 亚洲av欧美aⅴ国产| 免费黄色在线免费观看| 黄片无遮挡物在线观看| 欧美日韩亚洲高清精品| 久久久久精品久久久久真实原创| 高清在线视频一区二区三区| 爱豆传媒免费全集在线观看| 精品视频人人做人人爽| 男女高潮啪啪啪动态图| 亚洲第一青青草原| 啦啦啦中文免费视频观看日本| 在线天堂中文资源库| 精品少妇内射三级| 国产男女超爽视频在线观看| 午夜91福利影院| 满18在线观看网站| 水蜜桃什么品种好| 国产精品国产三级专区第一集| 丰满饥渴人妻一区二区三| 26uuu在线亚洲综合色| 欧美成人午夜免费资源| 亚洲一级一片aⅴ在线观看| 欧美97在线视频| 国产高清国产精品国产三级| 欧美+日韩+精品| 卡戴珊不雅视频在线播放| 丝袜美腿诱惑在线| 曰老女人黄片| 只有这里有精品99| 久久久久久久国产电影| 国产成人免费观看mmmm| 国产精品香港三级国产av潘金莲 | 五月开心婷婷网| 婷婷色综合www| 亚洲av欧美aⅴ国产| 亚洲av综合色区一区| 免费播放大片免费观看视频在线观看| 女人高潮潮喷娇喘18禁视频| 国产男人的电影天堂91| 午夜免费男女啪啪视频观看|