BACKGROUND
protein-rna interactions play fundamental roles in many biological processes, such as regulation of gene expression, rna splicing, protein synthesis, replication of viral rnas, and virus assembly .
defects of protein-rna interactions are responsible for many human diseases ranging from neurological disorders to cancer  <cit> .
the understanding of these processes improves as new structures of protein-rna complexes are solved and the molecular details of interactions analyzed.
unfortunately, the experimental determination of protein-rna complexes is a slow and difficult process  <cit> .
the ability to predict structures of protein-rna complexes computationally would greatly help us study protein-rna interactions.
however, while there is a growing number of methods for modeling protein and rna structures , the number of methods for modeling protein-rna complexes remains scarce.
docking methods are widely used to predict structures of macromolecular complexes, starting from structures of the individual components  <cit> .
all docking methods face two main challenges: to search the space of possible orientations and conformations  of the components and to identify near-native structures among the alternative complex models  generated.
an ideal docking method should be able to reliably reconstruct a native complex from its 'bound' components, and score it significantly better than any non-native decoys.
in real life, the structure of the complex is unknown, and while the structures of binding partners are solved in isolation, the task of the 'unbound' docking experiment is not only to assemble them into a complex, but also to take into account possible conformational changes upon binding.
conformational changes are either modeled explicitly at the atomic level , or a certain level of 'fuzziness' is introduced e.g.
by allowing for some extent of steric overlaps between atoms or by 'coarse-graining' of the representation i.e.
by neglecting some atoms or grouping them into 'united atoms' to be considered jointly .
one interesting and frequently neglected aspect of rna structure and interactions is the presence of posttranscriptional modifications, which increase the basic set of four nucleotides  to more than  <dig> variants with altered base and/or ribose moieties  <cit> .
modified residues in rna are involved in many processes, including rna folding and rna-rna interactions , but also specific rna-protein recognition and binding  <cit> .
to our knowledge, among freely available macromolecular docking methods only three are suitable for handling post-transcriptional nucleotide modifications in rna without the need of 'demodification'.
haddock accepts rna as a part of the complex to be predicted, but the user is required to provide force field parameters for all modified nucleotides  <cit> .
gramm  <cit>  and hex  <cit>  can also perform protein-rna docking for rna structures with modified nucleotides, but the scoring functions of these programs have been developed to evaluate protein-protein complexes, and while they can generate the poses for protein-rna complexes, they are unable to identify near-native variants from a set of decoys.
a useful extension of the latter kind of methods would be the development of a potential for scoring protein-rna complexes.
recently, new statistical potentials for scoring protein-rna complexes have been proposed: a distance-dependent all-atom potential developed by the varani group  <cit> , and a residue-level potential developed by the fernandez group  <cit> .
the varani potential performs well in discriminating models of protein-rna complexes that are very close to the native structure, i.e.
with the root mean square deviation  <  <dig> Å  <cit> .
however, during a real  docking experiment it is difficult to obtain many decoys with rmsd <  <dig> Å. the fernandez potential was designed to improve the discriminative power of the ftdock potential and is not available as a standalone program.
the ftdock program  <cit>  was developed for protein-protein and protein-dna docking, but it accepts only conventional rna molecules .
in this article we introduce two new, medium-resolution, knowledge-based potentials for scoring protein-rna models: the quasi-chemical potential  and the decoys as the reference state  <cit>  potential .
these potentials are based on a reduced representation of protein and rna, use the same mathematical base but differ in their reference state.
we compare the discriminative power of our new statistical potentials to the varani and fernandez potentials.
for the bound docking test set our potentials discriminated native-like  structures of protein-rna complexes, the potential developed by the varani group recognized structures very close to the native , whereas the fernandez potential recognized near native structures only for some protein-rna complexes.
for the unbound docking test set, our potentials have the highest discriminative ability of alternative models.
our new knowledge-based potentials are a useful tool for scoring protein-rna complexes generated by macromolecular docking methods, such as gramm or hex.
implementation
docking
to search the space of possible orientations and conformations of the components, we employed the gramm method for medium- to low-resolution docking  <cit> .
as opposed to high-resolution methods that typically operate in the continuous space, gramm discretizes the system  by projecting the macromolecular structures on a grid and allows for imprecise fit by 'softening' the van der waals interactions and permitting some degree of steric conflicts.
one of the components of a binary complex, referred to as the "receptor", is fixed, while the other component, referred to as the "ligand", is rotated and translated around the receptor to obtain geometrically compatible poses.
the van der waals radius was used as a projection of an atom.
the value of the grid step was set to  <dig>  Å for complexes from the bound docking set and to the minimal value allowed by the program for complexes from the unbound docking set, a repulsion parameter to  <dig>  and attraction double range to  <dig>  ligand structures were rotated by 10° angle intervals.
if significant steric clashes were observed in a large fraction of models obtained in preliminary gramm runs, as assessed by visual inspection, the repulsion parameter was increased to decrease the volume of steric conflict, until docking decoys reached physically reasonable geometry .
we defined "native-like" poses as those with the ligand rmsd <  <dig> Å from the native complex structure.
according to our experience, this value is appropriate for consideration of medium- and low-resolution docking experiments.
the distance between the experimentally determined and decoy complex structure was calculated as the rmsd of all heavy atoms of ligands, following optimal superposition of the receptor structures.
the resulting sets of decoys are supposed to approximate a broad distribution of structures that exhibit relatively good spatial complementarity between the receptor and the ligand, regardless of other interactions that contribute to the observed strength of binding.
such sets may be enriched into native-like decoys comparably to the completely random models , but "spatially complementary" models are much better approximation of the real-life sampling.
the false models were not selected by the gramm scores, as we found that the gramm score alone is a poor predictor of the complex native-likeness .
we treated all  <dig> decoys per complex as equal, both for training of our potential, and for evaluation with all four potentials analyzed in this work.
statistical potentials used
the distance and orientation-dependent knowledge-based potentials described in this article  were generated using reverse boltzmann statistics, where the interaction energy ε between the united atom type from the protein i and the united atom type from the rna j is calculated according to formula
ϵi,j,d=-rtlnnobsi,j,dnexpi,j,d
where r is the gas constant, t is the temperature, nobs  is the number of contacts between atom types i and j observed in the training set in a distance/angular bin d, while nexp  is the expected number of contacts at the same distance/angle in the reference state.
there are two types of bins, connected with different terms of the potential: a distance bin  used in the distance-dependent term, and an angular bin  used in the angular-dependent term.
the energy is calculated for each pair of protein-rna united atoms that are within the distance of  <dig> Å from each other.
this threshold parameter was chosen based on the analysis carried out by the shakhnovich group  <cit> , who identified  <dig> Å as the optimal threshold for evaluating protein-nucleic acid interactions with all-atom multiple bin potentials.
for our reduced representation potentials we have tested five distance thresholds  by calculating correlation coefficients and a position of the native structure in the decoys ranking.
the best results was for  <dig> Å for both dars and quasi potentials , and it was not statistically distinguishable from  <dig> Å and  <dig> Å .
our two statistical potentials calculate the reference state in different ways.
the reference state for the quasi-rnp potential is calculated using mole fractions of residues:
nexpi,j,d=xi*xj*nobsd,
where xi and xj are the mole fractions of atom type i and j respectively, while nobs  is the total number of contacts in bin d. on the other hand, in the dars approach  <cit> nexp  is a normalized number of contacts between atom types i and j in a set of decoy structures that are considered as random models.
<dig> decoys were generated for each native structure of protein-rna complex from the training set by using the gramm docking program  <cit> , with the following parameters: values for repulsion and attraction were  <dig> and  <dig>  respectively, while a grid step used the minimal value allowed by the program, depending on the size of the components to be docked.
in a few cases where default values led to generation of decoys with significant steric clashes, the repulsion value was increased stepwise, until physically realistic decoys were obtained.
it must be emphasized that these decoys maximize the geometric fit between the protein and rna molecule, but do not take any interaction energies into account.
both statistical potentials comprise a distance-dependent energy term , an angular-dependent energy term , a site-dependent energy term , and a penalty for steric clashes  :
e=ed+ea+es+ep
the site-dependent term assesses the probability of interaction of amino acid residues with edges of nucleotide residues: watson-crick, sugar and hoogsteen, as defined by leontis at.al  <cit>  .
all four terms of the energy function exhibit comparable values .
hence, we assigned equal weights to all four terms and have not optimized them further in any way.
additional file  <dig>  figures s <dig>  s <dig> and s <dig> show selected graphs of nobs common for dars-rnp and quasi-rnp, nexp for dars-rnp, and nexp for quasi-rnp, respectively, as a function of distance, angle and nucleotide edges, while additional files  <dig>   <dig> and  <dig> include the list of values for all cases.
representation of molecules
to calculate the parameters of quasi-rnp and dars-rnp potentials, the all-atom representation of all macromolecular structures was transformed into a reduced representation.
the atoms of each amino acid or nucleotide residue were replaced by a number of united atoms depending on the residue type.
for amino acid residues, we used the representation used in the refiner program  <cit> , which involves one to three united atoms per residue, depending on the size of the residue.
for nucleotide residues we employed the representation used in the redrna method, currently under development in our laboratory .
the rna backbone was represented by two united atoms, one for the phosphate group  and one for the ribose , while the pyrimidine and purine rings were represented by one and two atoms, respectively .
all united atoms and alpha carbons from the reduced representation in every residue were considered for the potentials as separate atom types, e.g.
the alpha carbon of alanine and alpha carbon of lysine had different types, because they represented different type of residues.
training set for developing the statistical potentials
in order to compare our quasi-rnp and dars-rnp potentials to the previously published varani potential we used the same training set as the varani group  <cit> .
there are  <dig> protein-rna complexes in the training set taken from crystal structures of protein-rna complexes downloaded from the protein data bank   <cit> .
one of these is the 50s ribosomal subunit, which comprises of  <dig> individual peptide chains in complex with rna.
due to the limited number of protein-rna complexes in this set, we performed a "leave one out" cross-validation, in which the potential was recalculated for testing of each structure, based on a training set with the tested structure excluded.
our software uses an in-house modified variant of pdbparser from biopython  to parse pdb files and ignore information about e.g.
atom occupancy.
testing statistical potentials for protein-rna docking
to test the discriminatory power of the quasi-rnp and dars-rnp knowledge-based potentials we compared them with two existing statistical potentials developed by the varani  <cit>  and fernandez  <cit>  groups.
software to calculate the varani potential has been kindly provided by the author .
since the fernandez potential is not available as a standalone software, the authors have kindly provided raw statistical data for each amino acid-nucleotide pair calculated based on their training set , which we used to calculate our local implementation of their scoring function, independent of the ftdock program.
two types of test sets were used, based on bound structures , or on unbound structures .
bound docking test sets
we used one set of molecules with unmodified  conformations of entire rna molecules and protein backbones, but with optimized protein side chains.
this set was developed by the varani group to perform decoy-discrimination tests of their all-atom potentials  <cit> .
their decoys were obtained by modifying five native protein-rna complexes  using the docking module of rosetta, through the use of the protein side chain repacking algorithm  <cit> .
for each protein complex they have generated  <dig> structures with the rmsd from the native complex structure ranging from  <dig>  Å to about  <dig> Å. the rmsd of two complexes was calculated based on ligand heavy atoms, after superposition of receptors.
the second set contains all molecules in unmodified  conformations, and was generated by ourselves using the high-resolution mode of the gramm docking program.
for each of the five protein-rna complexes from the varani test set  <dig>  alternative docking decoys were generated according to the procedure described above in the docking section.
unbound docking test set
the unbound docking test set was based on twelve native protein-rna complexes , previously used by the fernandez group to test their potential  <cit> .
for each component of these complexes, at least one independently solved 3d structure per complex is available.
the gramm program was used to generate  <dig>  docking decoys for each complex.
we used the same parameters as in the bound docking procedure.
with these settings, the gramm program generated at least five native-like structures  for eight out of twelve protein-rna complexes.
only these eight decoy sets with native-like structures were considered.
clustering the best scored decoys
critical assessment of protein structure predictions  has demonstrated that the scoring functions alone may not be the best discriminators of native-like structures, and better results may be achieved by clustering well-scored suboptimal structures  <cit> .
we have applied this approach, in particular using the clustering algorithm proposed by baker and coworkers  <cit> , which has worked very well in protein structure prediction.
first, an rmsd is calculated for all pairs of structures and stored in a distance matrix.
then, the row of the distance matrix with the largest number of rmsd values smaller than a cutoff  is selected.
structures in that row with rmsd below the cutoff value are assigned to one cluster and excluded from the matrix.
the process is iterated until all structures with rmsd smaller than the cutoff are assigned to clusters.
three biggest clusters are then considered as candidates for groups of potentially native-like structures, from which then the lowest-energy decoys are identified.
as a rule of thumb, the larger the cluster and the better score of its decoys - the higher chance that it contains a native-like structure.
RESULTS
in order to compare the ability to identify near-native protein-rna complex structures by our two potentials , the potential developed by varani et al.
<cit> , and the potential developed by fernandez et al.
<cit> , we used test sets, obtained by bound and unbound docking.
to give all methods equal chances in finding structures close to the native one, we included the clustering procedure implemented in our potentials as an additional step in the scoring performed by both the varani and fernandez potentials.
we have calculated correlation coefficients for decoys with rmsd in the range of 0- <dig> Å, 0- <dig> Å and 0- <dig> Å. the rmsd range <  <dig> Å was chosen because the varani potential has a recognition area of near-native structures around  <dig> Å. the rmsd threshold  <dig> Å is a "golden standard" of near to native decoys definition used in the capri  experiment.
the generous  <dig> Å threshold was used because we observed that for large structures there exist decoys with a biologically reasonable match of the binding sites and a number of native-like protein-rna contacts, despite the global rmsd of the ligand in the range of 10- <dig> Å. such deviations are typically due to a rotation of the ligand that retains contacts in the binding site, but moves around atoms that are far away from the binding site, as shown e.g.
in additional file  <dig>  figures s <dig> and s <dig>
decoy discrimination for the bound docking test set
we tested the dars-rnp, quasi-rnp, varani, and fernandez potentials for scoring of rna-protein complexes on two bound test sets, and examined which potential gives the best correlation coefficient between scores and rmsd of the corresponding decoys.
for the varani test set, the dars-rnp and quasi-rnp potentials recognized most structures with ligand rmsd <  <dig> Å. we found strong correlation between dars-rnp scores and rmsd, as well between quasi-rnp scores and rmsd for models up to  <dig> Å from the native structure .
the varani potential discriminates models with ligand rmsd <  <dig> Å, but generally fails to distinguish structures with 5- <dig> Å rmsd from the native and random structures with larger rmsd values .
the fernandez potential fails to recognize individual near-native decoys from the varani test set; and its scores exhibit positive correlation coefficients with rmsds only for decoys from three complexes in five .
the clustering procedure improves the results obtained by the fernendez potential, as the biggest clusters contain structures with rmsd <  <dig> Å for all complexes considered.
nonetheless, the native structures and structures with the smallest rmsd are scored too poorly by the fernandez potential to be included in the top clusters for two of five complexes - 1urn and 1ec <dig> .
correlation coefficients were calculated for scatter plots from figure  <dig>
in the bound docking test set generated by the gramm program, there is a smaller number of near-native structures than in the varani set, and some of them exhibit steric clashes.
for this test set, our dars-rnp and quasi-rnp potentials exhibit lower values of the correlation coefficient than for the varani test set .
the varani and fernanez potentials discriminate gramm decoys better than varani decoys for three complexes, and worse for two complexes   .
correlation coefficients were calculated for scatter plots from figure  <dig>
summarizing, the dars-rnp and quasi-rnp scores exhibit the highest correlation coefficients for all cutoffs in the varani test set, and for  <dig> Å and  <dig> Å thresholds in the gramm test set .
thus, dars-rnp and quasi-rnp potentials can be declared as "winners" of the bound docking test, except for the structures that are very close to the native complexes, where they are outperformed by the varani potential.
decoy discrimination for the unbound docking set
in an analogous way, we examined the discriminatory power of the dars-rnp, quasi-rnp, varani, and fernandez potentials for decoys of the unbound docking test set.
the assessment of unbound docking results reveals, expectedly, that all potentials exhibit worse results than for the bound docking set .
here, the best results are obtained by our dars-rnp potential, followed closely by the quasi-rnp potential.
these potentials recognized native-like structures for four out of eight complexes from the unbound test set , while the varani and fernandez potentials recognized native-like structures only for one complex in this set with the default options of clustering, and two and three complexes respectively after increasing both number of clustering structures and rmsd threshold to  <dig> and  <dig> Å .
the correlation coefficients between the model score and rmsd are relatively low for the varani potential  and for the fernandez potential , while the dars-rnp/quasi-rnp potentials exhibit correlation coefficients of  <dig> / <dig> ,  <dig> / <dig> , and  <dig> / <dig>  for decoys with rmsd from the native structure lesser than  <dig>   <dig>  and  <dig> Å, respectively  and .
rna
correlation coefficients were calculated for scatter plots from figure  <dig>
clustering of the best-scored decoys
the application of clustering to identify groups of similar structures among the top-scored decoys improves the predictive power of all statistical potentials analyzed in this work.
it helps fernandez potential to recognize near-native structures in the varani set for bound docking with optimization of side chains , our dars-rnp and quasi-rnp potentials in the gramm bound docking test set , and all potentials in the unbound docking test set .
figure  <dig> shows examples of four complexes where the native-structure was found owing to the clustering of well-scored models identified by the dars-rnp potential.
discussion
the quasi-rnp and dars-rnp potentials described in this work exhibit the highest discriminatory power for the bound varani set, where there are many near-native structures without steric clashes.
likewise, our methods performed well for another set of decoys generated for the same rna-protein complex structures with the gramm method.
our potentials failed to recognize native-like structures generated by gramm only for the 1jid structure .
both dars-rnp and quasi-rnp favored a structure that is very different from the native complex, even though they were able to recognize native-like structures for the same complex generated by varani.
the visual examination of models for the 1jid complex that were best-scored by our potentials revealed structures, in which the rna backbone has entered a very deep and narrow groove on the protein structure far from the true rna-binding site, leading to decoys with very big area of protein-rna interactions, and hence with significantly more contacts than in the native structure .
one way to avoid such situations is to identify  the rna-binding site on the protein surface and filter the initial decoys  to remove those with rna away from the binding site.
it is worth to mention that the five structures in the bound docking test sets were excluded from the training set only for the quasi-rnp and dars-rnp potentials.
the training set of the fernandez potential contained three out of five assessed complexes, and the varani training set contained all five complexes.
therefore, the ability of the varani and fernandez potentials to discriminate native-like models for these structures may be overestimated.
in particular, the varani potential has the best results for those decoys from the gramm-generated unbound set that are very close to native structure .
there, the varani potential easily recognizes the structures close to those in its training set.
interestingly, for decoys from the same set, with rmsd up to  <dig> or  <dig> Å, our potentials still exhibit better results than the varani potential, suggesting that they do have a power to discriminate between these medium-quality decoys and those that are totally native-unlike.
the results of tests for the unbound docking set are more objective due to a complete separation of training and testing data, and because they simulated the predictive power of the potentials in a real docking experiment, where the bound conformations of components are unknown.
among the four methods tested, quasi-rnp and dars-rnp potentials have the biggest average correlation coefficients between the scores and the rmsd of the model from the native structure.
however, it must be emphasized that even these "winners" of our benchmark were able to identify native-like structures only for four out of eight tested cases .
as expected, our potentials recognized near native structures only for these complexes, whose components exhibit relatively small structural changes  during complex formation .
still, in our hands, the fernandez and varani potentials recognized near-native structures for one complex only with the default options of clustering , and three and two complexes, respectively, after increasing the number of clustered best-scored decoys to  <dig> and the rmsd threshold to  <dig> Å .
such "relaxed" clusters are of course more heterogeneous.
for two complexes  all four potentials considered in this work have recognized native-like structures.
for both of these complexes only one component of the complex was solved in isolation from the other  and that 'unbound' component underwent only a very minor conformational change with respect to the bound form  .
hence, these two targets must be considered as very easy.
among the four knowledge-based potentials tested here for their ability to identify native-like protein-rna docking models, the high-resolution varani potential exhibited the best ability to recognize models closest to the native structure .
this potential appears as the method of choice for high-accuracy docking methods that are able to generate structures very close to the native ones.
it must be emphasized, however, that there are no computational tools, with which to reliably predict conformational changes upon binding.
therefore high-resolution docking is only applicable in situations, where the receptor and ligand structures exhibit very little conformational change between the bound and unbound forms.
unfortunately, whether the conformational change occurs or not cannot be reliably predicted.
in most cases of protein-rna binding, moderate conformational changes of protein and/or rna molecules occur upon complex formation.
there, low resolution methods that apply a coarse-grained energy model to a coarse-grained representation  have a chance to be practically useful.
among the four potentials tested in this work, the fernandez potential has the weakest discriminatory power for identification of near-native structures, but we believe this potential may perform much better when combined with the ftdock scoring function, as originally intended by the fernandez group.
however, the combination of fernandez and ftdock potentials is possible only for ftdock decoys, because the ftdock potential is not available as a standalone program, hence we could not apply it to the decoys generated in our study.
ftdock is also unable to deal with modified residues in the rna, which precludes it from applicability to many rna-protein complexes, where rna modifications play a critical role .
the main difference between all potentials considered in this article is the set of atoms taken into consideration.
the varani potential considers interactions between all atoms, with chemically similar atoms  treated in the same way.
it contains only a distance-dependent multiple bin term.
the fernandez potential calculates interactions between entire residues represented as single interaction centers, using only one bin .
our potentials use multiple bins for distance as well as orientation, hence they take into account more information about the possible arrangements of amino acid and nucleotide residues, even though they use less atoms than the varani method.
in our study we have used both bound and unbound conformations for docking .
the shapes of score-rmsd dependence plots show differences associated not only with the methods, but also with the type of docking.
as expected, all potentials exhibit best results for bound docking experiment, where there are many near-native structures without steric clashes.
this observation underlines the influence of decoy generation method on the ability to successfully identify native-like decoys in the generated dataset.
our study allows for direct comparison of the decoy-based and quasi-chemical approaches for calculating statistical potentials.
the small, but significant average advantage of the dars-rnp potential over the quasi-rnp potential can be explained by the more realistic treatment of "random" protein-rna interactions.
in the dars-based approach, the statistics of amino acid-nucleotide contacts are inferred from geometrically plausible, but biologically irrelevant decoys , while the quasi-chemical approach predicts the occurrence of such contacts based on the frequency of individual residues.
we expect this advantage of the dars approach to hold also for other types of docking.
the calculation of a dars-based potential requires, however, the calculation of a large number of decoys for each complex in the training set, hence it requires considerably bigger computational effort, which may be prohibitive in case of large training sets.
by definition, none of the rigid-body docking methods analyzed here is capable of predicting the structures of complexes that involve large conformational changes.
we also found that the presence of extensive steric clashes in decoys deteriorates the discriminatory power of all potentials tested in our benchmark.
thus, we propose that the next step in the development of methods for modeling of protein-rna complexes should be taken towards algorithms that enable simultaneous docking and folding of protein and rna components.
recently, a number of methods for modeling of rna 3d structures have been reported that utilize very similar methodology to that used for protein modeling  <cit> .
this suggests that the combination of comparative modeling and "de novo" folding should be possible not only for proteins and rna separately, but also as components of the same molecular system.
CONCLUSIONS
among the four potentials tested in this work the quasi-rnp and dars-rnp potentials exhibit the highest discriminatory power for both bound and unbound docking experiments.
the small average advantage of the dars-rnp potential over the quasi-rnp potential can be caused by the more realistic treatment of "random" protein-rna interactions.
none of the rigid-body docking methods analyzed here is capable of predicting the structures of complexes that involve large conformational changes.
our potentials recognized near native structures only for these complexes, whose components exhibit relatively small structural changes  during formation of the protein-rna complex.
availability and requirements
project name: statistical potentials for protein-rna docking; project home page: http://iimcb.genesilico.pl/rnp/; operating system: linux; programming language: python; other requirements: biopython and numpy libraries; licence: gnu gpl; no restrictions to use by non-academics.
list of abbreviations
rmsd: root mean square deviation; rnp: rna-protein complexes; quasi-rnp: quasi-chemical potential for rna-protein complexes; dars-rnp: decoys as reference state potential for rna-protein complexes.
authors' contributions
it developed statistical potentials, performed testing procedures, drafted the manuscript and prepared the figures.
jmb conceived of the project and edited the manuscript.
both authors analyzed and interpreted the data.
all authors read and approved the final manuscript.
supplementary material
additional file 1
additional_figures_and_table.pdf.
files and tables with additional data illustrating the details of implementation, ordered according to their appearance in the text.
click here for file
additional file 2
energy_dars.pdf.
energy for each distance, angle, and site bin, for each pair wise interaction in the dars-rnp potential.
click here for file
additional file 3
energy_quasi.pdf.
energy for each distance, angle, and site bin, for each pair wise interaction in the quasi-rnp potential.
click here for file
additional file 4
observed.pdf.
observed number of contacts in each distance, angle, and site bin, for each pair wise interaction.
click here for file
additional file 5
expected_dars.pdf.
expected number of contacts in each distance, angle, and site bin, for each pair wise interaction in the dars-rnp potential.
click here for file
additional file 6
expected_quasi.pdf.
expected number of contacts in each distance, angle, and site bin, for each pair wise interaction in the quasi-rnp potential.
click here for file
acknowledgements and funding
we thank michal boniecki and marcin pawlowski for inspiration and discussions and michal pietal for sharing brutepdbparser.
we also wish to thank kristian rother, grzegorz chojnowski, anna philips, and lukasz kozlowski for critical reading of the manuscript.
we thank gabriele varani and juan fernandez-recio for providing us with their software and data for comparison with our methods.
this project was funded by the foundation for polish science .
i.t.
was supported by the 6th framework programme of the european commission .
j.m.b.
was supported by the european research council  and by the "ideas for poland" fellowship from the foundation of polish science .

 protein-rna interactions play fundamental roles in many biological processes, such as regulation of gene expression, rna splicing, protein synthesis, replication of viral rnas, and virus assembly. defects of protein-rna interactions are responsible for many human diseases ranging from neurological disorders to cancer. the understanding of these processes improves as new structures of protein-rna complexes are solved and the molecular details of interactions analyzed. unfortunately, the experimental determination of protein-rna complexes is a slow and difficult process. the ability to predict structures of protein-rna complexes computationally would greatly help us study protein-rna interactions. however, while there is a growing number of methods for modeling protein and rna structures, the number of methods for modeling protein-rna complexes remains scarce. docking methods are widely used to predict structures of macromolecular complexes, starting from structures of the individual components. all docking methods face two main challenges: to search the space of possible orientations and conformations of the components and to identify near-native structures among the alternative complex models generated. an ideal docking method should be able to reliably reconstruct a native complex from its'bound' components, and score it significantly better than any non-native decoys. in real life, the structure of the complex is unknown, and while the structures of binding partners are solved in isolation, the task of the'unbound' docking experiment is not only to assemble them into a complex, but also to take into account possible conformational changes upon binding. conformational changes are either modeled explicitly at the atomic level, or a certain level of'fuzziness' is introduced e.g. by allowing for some extent of steric overlaps between atoms or by'coarse-graining' of the representation i.e. by neglecting some atoms or grouping them into'united atoms' to be considered jointly. one interesting and frequently neglected aspect of rna structure and interactions is the presence of posttranscriptional modifications, which increase the basic set of four nucleotides to more than variants with altered base and/or ribose moieties. modified residues in rna are involved in many processes, including rna folding and rna-rna interactions, but also specific rna-protein recognition and binding. to our knowledge, among freely available macromolecular docking methods only three are suitable for handling post-transcriptional nucleotide modifications in rna without the need of'demodification'. haddock accepts rna as a part of the complex to be predicted, but the user is required to provide force field parameters for all modified nucleotides. gramm and hex can also perform protein-rna docking for rna structures with modified nucleotides, but the scoring functions of these programs have been developed to evaluate protein-protein complexes, and while they can generate the poses for protein-rna complexes, they are unable to identify near-native variants from a set of decoys. a useful extension of the latter kind of methods would be the development of a potential for scoring protein-rna complexes. recently, new statistical potentials for scoring protein-rna complexes have been proposed: a distance-dependent all-atom potential developed by the varani group, and a residue-level potential developed by the fernandez group. the varani potential performs well in discriminating models of protein-rna complexes that are very close to the native structure, i.e. with the root mean square deviation < Å. however, during a real docking experiment it is difficult to obtain many decoys with rmsd < Å. the fernandez potential was designed to improve the discriminative power of the ftdock potential and is not available as a standalone program. the ftdock program was developed for protein-protein and protein-dna docking, but it accepts only conventional rna molecules. in this article we introduce two new, medium-resolution, knowledge-based potentials for scoring protein-rna models: the quasi-chemical potential and the decoys as the reference state potential. these potentials are based on a reduced representation of protein and rna, use the same mathematical base but differ in their reference state. we compare the discriminative power of our new statistical potentials to the varani and fernandez potentials. for the bound docking test set our potentials discriminated native-like structures of protein-rna complexes, the potential developed by the varani group recognized structures very close to the native, whereas the fernandez potential recognized near native structures only for some protein-rna complexes. for the unbound docking test set, our potentials have the highest discriminative ability of alternative models. our new knowledge-based potentials are a useful tool for scoring protein-rna complexes generated by macromolecular docking methods, such as gramm or hex. implementation docking to search the space of possible orientations and conformations of the components, we employed the gramm method for medium- to low-resolution docking. as opposed to high-resolution methods that typically operate in the continuous space, gramm discretizes the system by projecting the macromolecular structures on a grid and allows for imprecise fit by'softening' the van der waals interactions and permitting some degree of steric conflicts. one of the components of a binary complex, referred to as the "receptor", is fixed, while the other component, referred to as the "ligand", is rotated and translated around the receptor to obtain geometrically compatible poses. the van der waals radius was used as a projection of an atom. the value of the grid step was set to Å for complexes from the bound docking set and to the minimal value allowed by the program for complexes from the unbound docking set, a repulsion parameter to and attraction double range to ligand structures were rotated by 10° angle intervals. if significant steric clashes were observed in a large fraction of models obtained in preliminary gramm runs, as assessed by visual inspection, the repulsion parameter was increased to decrease the volume of steric conflict, until docking decoys reached physically reasonable geometry. we defined "native-like" poses as those with the ligand rmsd < Å from the native complex structure. according to our experience, this value is appropriate for consideration of medium- and low-resolution docking experiments. the distance between the experimentally determined and decoy complex structure was calculated as the rmsd of all heavy atoms of ligands, following optimal superposition of the receptor structures. the resulting sets of decoys are supposed to approximate a broad distribution of structures that exhibit relatively good spatial complementarity between the receptor and the ligand, regardless of other interactions that contribute to the observed strength of binding. such sets may be enriched into native-like decoys comparably to the completely random models, but "spatially complementary" models are much better approximation of the real-life sampling. the false models were not selected by the gramm scores, as we found that the gramm score alone is a poor predictor of the complex native-likeness. we treated all decoys per complex as equal, both for training of our potential, and for evaluation with all four potentials analyzed in this work. statistical potentials used the distance and orientation-dependent knowledge-based potentials described in this article were generated using reverse boltzmann statistics, where the interaction energy ε between the united atom type from the protein i and the united atom type from the rna j is calculated according to formula ϵi,j,d=-rtlnnobsi,j,dnexpi,j,d where r is the gas constant, t is the temperature, nobs is the number of contacts between atom types i and j observed in the training set in a distance/angular bin d, while nexp is the expected number of contacts at the same distance/angle in the reference state. there are two types of bins, connected with different terms of the potential: a distance bin used in the distance-dependent term, and an angular bin used in the angular-dependent term. the energy is calculated for each pair of protein-rna united atoms that are within the distance of Å from each other. this threshold parameter was chosen based on the analysis carried out by the shakhnovich group, who identified Å as the optimal threshold for evaluating protein-nucleic acid interactions with all-atom multiple bin potentials. for our reduced representation potentials we have tested five distance thresholds by calculating correlation coefficients and a position of the native structure in the decoys ranking. the best results was for Å for both dars and quasi potentials, and it was not statistically distinguishable from Å and Å. our two statistical potentials calculate the reference state in different ways. the reference state for the quasi-rnp potential is calculated using mole fractions of residues: nexpi,j,d=xi*xj*nobsd, where xi and xj are the mole fractions of atom type i and j respectively, while nobs is the total number of contacts in bin d. on the other hand, in the dars approach nexp is a normalized number of contacts between atom types i and j in a set of decoy structures that are considered as random models. decoys were generated for each native structure of protein-rna complex from the training set by using the gramm docking program, with the following parameters: values for repulsion and attraction were and respectively, while a grid step used the minimal value allowed by the program, depending on the size of the components to be docked. in a few cases where default values led to generation of decoys with significant steric clashes, the repulsion value was increased stepwise, until physically realistic decoys were obtained. it must be emphasized that these decoys maximize the geometric fit between the protein and rna molecule, but do not take any interaction energies into account. both statistical potentials comprise a distance-dependent energy term, an angular-dependent energy term, a site-dependent energy term, and a penalty for steric clashes: e=ed+ea+es+ep the site-dependent term assesses the probability of interaction of amino acid residues with edges of nucleotide residues: watson-crick, sugar and hoogsteen, as defined by leontis at.al. all four terms of the energy function exhibit comparable values. hence, we assigned equal weights to all four terms and have not optimized them further in any way. additional file figures s s and s show selected graphs of nobs common for dars-rnp and quasi-rnp, nexp for dars-rnp, and nexp for quasi-rnp, respectively, as a function of distance, angle and nucleotide edges, while additional files and include the list of values for all cases. representation of molecules to calculate the parameters of quasi-rnp and dars-rnp potentials, the all-atom representation of all macromolecular structures was transformed into a reduced representation. the atoms of each amino acid or nucleotide residue were replaced by a number of united atoms depending on the residue type. for amino acid residues, we used the representation used in the refiner program, which involves one to three united atoms per residue, depending on the size of the residue. for nucleotide residues we employed the representation used in the redrna method, currently under development in our laboratory. the rna backbone was represented by two united atoms, one for the phosphate group and one for the ribose, while the pyrimidine and purine rings were represented by one and two atoms, respectively. all united atoms and alpha carbons from the reduced representation in every residue were considered for the potentials as separate atom types, e.g. the alpha carbon of alanine and alpha carbon of lysine had different types, because they represented different type of residues. training set for developing the statistical potentials in order to compare our quasi-rnp and dars-rnp potentials to the previously published varani potential we used the same training set as the varani group. there are protein-rna complexes in the training set taken from crystal structures of protein-rna complexes downloaded from the protein data bank. one of these is the 50s ribosomal subunit, which comprises of individual peptide chains in complex with rna. due to the limited number of protein-rna complexes in this set, we performed a "leave one out" cross-validation, in which the potential was recalculated for testing of each structure, based on a training set with the tested structure excluded. our software uses an in-house modified variant of pdbparser from biopython to parse pdb files and ignore information about e.g. atom occupancy. testing statistical potentials for protein-rna docking to test the discriminatory power of the quasi-rnp and dars-rnp knowledge-based potentials we compared them with two existing statistical potentials developed by the varani and fernandez groups. software to calculate the varani potential has been kindly provided by the author. since the fernandez potential is not available as a standalone software, the authors have kindly provided raw statistical data for each amino acid-nucleotide pair calculated based on their training set, which we used to calculate our local implementation of their scoring function, independent of the ftdock program. two types of test sets were used, based on bound structures, or on unbound structures. bound docking test sets we used one set of molecules with unmodified conformations of entire rna molecules and protein backbones, but with optimized protein side chains. this set was developed by the varani group to perform decoy-discrimination tests of their all-atom potentials. their decoys were obtained by modifying five native protein-rna complexes using the docking module of rosetta, through the use of the protein side chain repacking algorithm. for each protein complex they have generated structures with the rmsd from the native complex structure ranging from Å to about Å. the rmsd of two complexes was calculated based on ligand heavy atoms, after superposition of receptors. the second set contains all molecules in unmodified conformations, and was generated by ourselves using the high-resolution mode of the gramm docking program. for each of the five protein-rna complexes from the varani test set alternative docking decoys were generated according to the procedure described above in the docking section. unbound docking test set the unbound docking test set was based on twelve native protein-rna complexes, previously used by the fernandez group to test their potential. for each component of these complexes, at least one independently solved 3d structure per complex is available. the gramm program was used to generate docking decoys for each complex. we used the same parameters as in the bound docking procedure. with these settings, the gramm program generated at least five native-like structures for eight out of twelve protein-rna complexes. only these eight decoy sets with native-like structures were considered. clustering the best scored decoys critical assessment of protein structure predictions has demonstrated that the scoring functions alone may not be the best discriminators of native-like structures, and better results may be achieved by clustering well-scored suboptimal structures. we have applied this approach, in particular using the clustering algorithm proposed by baker and coworkers, which has worked very well in protein structure prediction. first, an rmsd is calculated for all pairs of structures and stored in a distance matrix. then, the row of the distance matrix with the largest number of rmsd values smaller than a cutoff is selected. structures in that row with rmsd below the cutoff value are assigned to one cluster and excluded from the matrix. the process is iterated until all structures with rmsd smaller than the cutoff are assigned to clusters. three biggest clusters are then considered as candidates for groups of potentially native-like structures, from which then the lowest-energy decoys are identified. as a rule of thumb, the larger the cluster and the better score of its decoys - the higher chance that it contains a native-like structure.  in order to compare the ability to identify near-native protein-rna complex structures by our two potentials, the potential developed by varani et al., and the potential developed by fernandez et al., we used test sets, obtained by bound and unbound docking. to give all methods equal chances in finding structures close to the native one, we included the clustering procedure implemented in our potentials as an additional step in the scoring performed by both the varani and fernandez potentials. we have calculated correlation coefficients for decoys with rmsd in the range of 0- Å, 0- Å and 0- Å. the rmsd range < Å was chosen because the varani potential has a recognition area of near-native structures around Å. the rmsd threshold Å is a "golden standard" of near to native decoys definition used in the capri experiment. the generous Å threshold was used because we observed that for large structures there exist decoys with a biologically reasonable match of the binding sites and a number of native-like protein-rna contacts, despite the global rmsd of the ligand in the range of 10- Å. such deviations are typically due to a rotation of the ligand that retains contacts in the binding site, but moves around atoms that are far away from the binding site, as shown e.g. in additional file figures s and s decoy discrimination for the bound docking test set we tested the dars-rnp, quasi-rnp, varani, and fernandez potentials for scoring of rna-protein complexes on two bound test sets, and examined which potential gives the best correlation coefficient between scores and rmsd of the corresponding decoys. for the varani test set, the dars-rnp and quasi-rnp potentials recognized most structures with ligand rmsd < Å. we found strong correlation between dars-rnp scores and rmsd, as well between quasi-rnp scores and rmsd for models up to Å from the native structure. the varani potential discriminates models with ligand rmsd < Å, but generally fails to distinguish structures with 5- Å rmsd from the native and random structures with larger rmsd values. the fernandez potential fails to recognize individual near-native decoys from the varani test set; and its scores exhibit positive correlation coefficients with rmsds only for decoys from three complexes in five. the clustering procedure improves the results obtained by the fernendez potential, as the biggest clusters contain structures with rmsd < Å for all complexes considered. nonetheless, the native structures and structures with the smallest rmsd are scored too poorly by the fernandez potential to be included in the top clusters for two of five complexes - 1urn and 1ec. correlation coefficients were calculated for scatter plots from figure in the bound docking test set generated by the gramm program, there is a smaller number of near-native structures than in the varani set, and some of them exhibit steric clashes. for this test set, our dars-rnp and quasi-rnp potentials exhibit lower values of the correlation coefficient than for the varani test set. the varani and fernanez potentials discriminate gramm decoys better than varani decoys for three complexes, and worse for two complexes. correlation coefficients were calculated for scatter plots from figure summarizing, the dars-rnp and quasi-rnp scores exhibit the highest correlation coefficients for all cutoffs in the varani test set, and for Å and Å thresholds in the gramm test set. thus, dars-rnp and quasi-rnp potentials can be declared as "winners" of the bound docking test, except for the structures that are very close to the native complexes, where they are outperformed by the varani potential. decoy discrimination for the unbound docking set in an analogous way, we examined the discriminatory power of the dars-rnp, quasi-rnp, varani, and fernandez potentials for decoys of the unbound docking test set. the assessment of unbound docking results reveals, expectedly, that all potentials exhibit worse results than for the bound docking set. here, the best results are obtained by our dars-rnp potential, followed closely by the quasi-rnp potential. these potentials recognized native-like structures for four out of eight complexes from the unbound test set, while the varani and fernandez potentials recognized native-like structures only for one complex in this set with the default options of clustering, and two and three complexes respectively after increasing both number of clustering structures and rmsd threshold to and Å. the correlation coefficients between the model score and rmsd are relatively low for the varani potential and for the fernandez potential, while the dars-rnp/quasi-rnp potentials exhibit correlation coefficients of /, /, and / for decoys with rmsd from the native structure lesser than and Å, respectively and. rna correlation coefficients were calculated for scatter plots from figure clustering of the best-scored decoys the application of clustering to identify groups of similar structures among the top-scored decoys improves the predictive power of all statistical potentials analyzed in this work. it helps fernandez potential to recognize near-native structures in the varani set for bound docking with optimization of side chains, our dars-rnp and quasi-rnp potentials in the gramm bound docking test set, and all potentials in the unbound docking test set. figure shows examples of four complexes where the native-structure was found owing to the clustering of well-scored models identified by the dars-rnp potential. discussion the quasi-rnp and dars-rnp potentials described in this work exhibit the highest discriminatory power for the bound varani set, where there are many near-native structures without steric clashes. likewise, our methods performed well for another set of decoys generated for the same rna-protein complex structures with the gramm method. our potentials failed to recognize native-like structures generated by gramm only for the 1jid structure. both dars-rnp and quasi-rnp favored a structure that is very different from the native complex, even though they were able to recognize native-like structures for the same complex generated by varani. the visual examination of models for the 1jid complex that were best-scored by our potentials revealed structures, in which the rna backbone has entered a very deep and narrow groove on the protein structure far from the true rna-binding site, leading to decoys with very big area of protein-rna interactions, and hence with significantly more contacts than in the native structure. one way to avoid such situations is to identify the rna-binding site on the protein surface and filter the initial decoys to remove those with rna away from the binding site. it is worth to mention that the five structures in the bound docking test sets were excluded from the training set only for the quasi-rnp and dars-rnp potentials. the training set of the fernandez potential contained three out of five assessed complexes, and the varani training set contained all five complexes. therefore, the ability of the varani and fernandez potentials to discriminate native-like models for these structures may be overestimated. in particular, the varani potential has the best results for those decoys from the gramm-generated unbound set that are very close to native structure. there, the varani potential easily recognizes the structures close to those in its training set. interestingly, for decoys from the same set, with rmsd up to or Å, our potentials still exhibit better results than the varani potential, suggesting that they do have a power to discriminate between these medium-quality decoys and those that are totally native-unlike. the results of tests for the unbound docking set are more objective due to a complete separation of training and testing data, and because they simulated the predictive power of the potentials in a real docking experiment, where the bound conformations of components are unknown. among the four methods tested, quasi-rnp and dars-rnp potentials have the biggest average correlation coefficients between the scores and the rmsd of the model from the native structure. however, it must be emphasized that even these "winners" of our benchmark were able to identify native-like structures only for four out of eight tested cases. as expected, our potentials recognized near native structures only for these complexes, whose components exhibit relatively small structural changes during complex formation. still, in our hands, the fernandez and varani potentials recognized near-native structures for one complex only with the default options of clustering, and three and two complexes, respectively, after increasing the number of clustered best-scored decoys to and the rmsd threshold to Å. such "relaxed" clusters are of course more heterogeneous. for two complexes all four potentials considered in this work have recognized native-like structures. for both of these complexes only one component of the complex was solved in isolation from the other and that'unbound' component underwent only a very minor conformational change with respect to the bound form. hence, these two targets must be considered as very easy. among the four knowledge-based potentials tested here for their ability to identify native-like protein-rna docking models, the high-resolution varani potential exhibited the best ability to recognize models closest to the native structure. this potential appears as the method of choice for high-accuracy docking methods that are able to generate structures very close to the native ones. it must be emphasized, however, that there are no computational tools, with which to reliably predict conformational changes upon binding. therefore high-resolution docking is only applicable in situations, where the receptor and ligand structures exhibit very little conformational change between the bound and unbound forms. unfortunately, whether the conformational change occurs or not cannot be reliably predicted. in most cases of protein-rna binding, moderate conformational changes of protein and/or rna molecules occur upon complex formation. there, low resolution methods that apply a coarse-grained energy model to a coarse-grained representation have a chance to be practically useful. among the four potentials tested in this work, the fernandez potential has the weakest discriminatory power for identification of near-native structures, but we believe this potential may perform much better when combined with the ftdock scoring function, as originally intended by the fernandez group. however, the combination of fernandez and ftdock potentials is possible only for ftdock decoys, because the ftdock potential is not available as a standalone program, hence we could not apply it to the decoys generated in our study. ftdock is also unable to deal with modified residues in the rna, which precludes it from applicability to many rna-protein complexes, where rna modifications play a critical role. the main difference between all potentials considered in this article is the set of atoms taken into consideration. the varani potential considers interactions between all atoms, with chemically similar atoms treated in the same way. it contains only a distance-dependent multiple bin term. the fernandez potential calculates interactions between entire residues represented as single interaction centers, using only one bin. our potentials use multiple bins for distance as well as orientation, hence they take into account more information about the possible arrangements of amino acid and nucleotide residues, even though they use less atoms than the varani method. in our study we have used both bound and unbound conformations for docking. the shapes of score-rmsd dependence plots show differences associated not only with the methods, but also with the type of docking. as expected, all potentials exhibit best results for bound docking experiment, where there are many near-native structures without steric clashes. this observation underlines the influence of decoy generation method on the ability to successfully identify native-like decoys in the generated dataset. our study allows for direct comparison of the decoy-based and quasi-chemical approaches for calculating statistical potentials. the small, but significant average advantage of the dars-rnp potential over the quasi-rnp potential can be explained by the more realistic treatment of "random" protein-rna interactions. in the dars-based approach, the statistics of amino acid-nucleotide contacts are inferred from geometrically plausible, but biologically irrelevant decoys, while the quasi-chemical approach predicts the occurrence of such contacts based on the frequency of individual residues. we expect this advantage of the dars approach to hold also for other types of docking. the calculation of a dars-based potential requires, however, the calculation of a large number of decoys for each complex in the training set, hence it requires considerably bigger computational effort, which may be prohibitive in case of large training sets. by definition, none of the rigid-body docking methods analyzed here is capable of predicting the structures of complexes that involve large conformational changes. we also found that the presence of extensive steric clashes in decoys deteriorates the discriminatory power of all potentials tested in our benchmark. thus, we propose that the next step in the development of methods for modeling of protein-rna complexes should be taken towards algorithms that enable simultaneous docking and folding of protein and rna components. recently, a number of methods for modeling of rna 3d structures have been reported that utilize very similar methodology to that used for protein modeling. this suggests that the combination of comparative modeling and "de novo" folding should be possible not only for proteins and rna separately, but also as components of the same molecular system.  among the four potentials tested in this work the quasi-rnp and dars-rnp potentials exhibit the highest discriminatory power for both bound and unbound docking experiments. the small average advantage of the dars-rnp potential over the quasi-rnp potential can be caused by the more realistic treatment of "random" protein-rna interactions. none of the rigid-body docking methods analyzed here is capable of predicting the structures of complexes that involve large conformational changes. our potentials recognized near native structures only for these complexes, whose components exhibit relatively small structural changes during formation of the protein-rna complex. availability and requirements project name: statistical potentials for protein-rna docking; project home page: http://iimcb.genesilico.pl/rnp/; operating system: linux; programming language: python; other requirements: biopython and numpy libraries; licence: gnu gpl; no restrictions to use by non-academics. list of abbreviations rmsd: root mean square deviation; rnp: rna-protein complexes; quasi-rnp: quasi-chemical potential for rna-protein complexes; dars-rnp: decoys as reference state potential for rna-protein complexes. authors' contributions it developed statistical potentials, performed testing procedures, drafted the manuscript and prepared the figures. jmb conceived of the project and edited the manuscript. both authors analyzed and interpreted the data. all authors read and approved the final manuscript. supplementary material additional file 1 additional_figures_and_table.pdf. files and tables with additional data illustrating the details of implementation, ordered according to their appearance in the text. click here for file additional file 2 energy_dars.pdf. energy for each distance, angle, and site bin, for each pair wise interaction in the dars-rnp potential. click here for file additional file 3 energy_quasi.pdf. energy for each distance, angle, and site bin, for each pair wise interaction in the quasi-rnp potential. click here for file additional file 4 observed.pdf. observed number of contacts in each distance, angle, and site bin, for each pair wise interaction. click here for file additional file 5 expected_dars.pdf. expected number of contacts in each distance, angle, and site bin, for each pair wise interaction in the dars-rnp potential. click here for file additional file 6 expected_quasi.pdf. expected number of contacts in each distance, angle, and site bin, for each pair wise interaction in the quasi-rnp potential. click here for file acknowledgements and funding we thank michal boniecki and marcin pawlowski for inspiration and discussions and michal pietal for sharing brutepdbparser. we also wish to thank kristian rother, grzegorz chojnowski, anna philips, and lukasz kozlowski for critical reading of the manuscript. we thank gabriele varani and juan fernandez-recio for providing us with their software and data for comparison with our methods. this project was funded by the foundation for polish science. i.t. was supported by the 6th framework programme of the european commission. j.m.b. was supported by the european research council and by the "ideas for poland" fellowship from the foundation of polish science.
we developed two medium-resolution, knowledge-based potentials for scoring protein-rna models obtained by docking: the quasi-chemical potential and the decoys as the reference state potential. we compared the discriminative power of dars-rnp and quasi-rnp for selecting rigid-body docking poses with the potentials previously developed by the varani and fernandez groups. protein-rna interactions play fundamental roles in many biological processes. both potentials use a coarse-grained representation for both rna and protein molecules and are capable of dealing with rna structures with posttranscriptionally modified residues. in both bound and unbound docking tests, dars-rnp showed the highest ability to identify native-like structures. understanding the molecular mechanism of protein-rna recognition and formation of protein-rna complexes is a major challenge in structural biology. unfortunately, the experimental determination of protein-rna complexes is tedious and difficult, both by x-ray crystallography and nmr. for many interacting proteins and rnas the individual structures are available, enabling computational prediction of complex structures by computational docking.
 protein-rna interactions play fundamental roles in many biological processes. understanding the molecular mechanism of protein-rna recognition and formation of protein-rna complexes is a major challenge in structural biology. unfortunately, the experimental determination of protein-rna complexes is tedious and difficult, both by x-ray crystallography and nmr. for many interacting proteins and rnas the individual structures are available, enabling computational prediction of complex structures by computational docking. however, methods for protein-rna docking remain scarce, in particular in comparison to the numerous methods for protein-protein docking.  we developed two medium-resolution, knowledge-based potentials for scoring protein-rna models obtained by docking: the quasi-chemical potential and the decoys as the reference state potential. both potentials use a coarse-grained representation for both rna and protein molecules and are capable of dealing with rna structures with posttranscriptionally modified residues. we compared the discriminative power of dars-rnp and quasi-rnp for selecting rigid-body docking poses with the potentials previously developed by the varani and fernandez groups.  in both bound and unbound docking tests, dars-rnp showed the highest ability to identify native-like structures. python implementations of dars-rnp and quasi-rnp are freely available for download at http://iimcb.genesilico.pl/rnp/ rnaproteinrnpmacromolecular dockingcomplex modelingstructural bioinformatics 