Molecular evolution of dihydrouridine synthases
© Kasprzak et al.; licensee BioMed Central Ltd. 2012
Received: 17 June 2011
Accepted: 24 May 2012
Published: 28 June 2012
Dihydrouridine (D) is a modified base found in conserved positions in the D-loop of tRNA in Bacteria, Eukaryota, and some Archaea. Despite the abundant occurrence of D, little is known about its biochemical roles in mediating tRNA function. It is assumed that D may destabilize the structure of tRNA and thus enhance its conformational flexibility. D is generated post-transcriptionally by the reduction of the 5,6-double bond of a uridine residue in RNA transcripts. The reaction is carried out by dihydrouridine synthases (DUS). DUS constitute a conserved family of enzymes encoded by the orthologous gene family COG0042. In protein sequence databases, members of COG0042 are typically annotated as “predicted TIM-barrel enzymes, possibly dehydrogenases, nifR3 family”.
To elucidate sequence-structure-function relationships in the DUS family, a comprehensive bioinformatic analysis was carried out. We performed extensive database searches to identify all members of the currently known DUS family, followed by clustering analysis to subdivide it into subfamilies of closely related sequences. We analyzed phylogenetic distributions of all members of the DUS family and inferred the evolutionary tree, which suggested a scenario for the evolutionary origin of dihydrouridine-forming enzymes. For a human representative of the DUS family, the hDus2 protein suggested as a potential drug target in cancer, we generated a homology model. While this article was under review, a crystal structure of a DUS representative has been published, giving us an opportunity to validate the model.
We compared sequences and phylogenetic distributions of all members of the DUS family and inferred the phylogenetic tree, which provides a framework to study the functional differences among these proteins and suggests a scenario for the evolutionary origin of dihydrouridine formation. Our evolutionary and structural classification of the DUS family provides a background to study functional differences among these proteins that will guide experimental analyses.
KeywordsDihydrouridine synthases Protein structure prediction Fold recognition Remote homology RNA modification Molecular evolution Enzymes acting on RNA
D is introduced in the RNA by dihydrouridine synthases (DUS). They constitute a conserved family of enzymes encoded by the orthologous gene family COG0042. In protein sequence databases, members of COG0042 are typically annotated as “predicted TIM-barrel enzymes, possibly dehydrogenases, nifR3 family”. They are found in all free-living organisms whose genomes have been sequenced. For instance, in E. coli three DUS have been identified (DusA, DusB, DusC), while Saccharomyces cerevisiae has four members of this family (Dus1p, Dus2p, Dus3p, Dus4p)  (see the most up to date list of enzymes in the MODOMICS database ). All yeast DUS are substrate specific and can modify only one or two positions in a tRNA molecule. It was shown  that Dus1p and Dus4p are able to modify two different positions in tRNA, U16/U17 and U20a/U20b, respectively. The other two yeast DUS, Dus2p and Dus3p, catalyze modification of unique positions in tRNA: 20 and 47, respectively. The targets of Dus1p, Dus2p, and Dus4p are within the D-loop of tRNA, while Dus3p acts on residues within a variable loop [6, 9]. Among the three DUS identified in E. coli, it has been reported that DusA (formerly YjbN) acts on position 21 of tRNA . It has been noticed that while the Mycoplasma mycoides genome encodes only a single DUS, its tRNAs contain D at many different positions .
Dihydropirimidine dehydrogenase (DHPDH)  and dihydroorotate dehydrogenase (DHODH)  are enzymes with activities similar to DUS. DHPDH catalyzes NADPH-dependent reduction of uracil to 5,6-dihydrouracil, while DHODH catalyzes oxidation of dihydroorotate to orotate in the pirimidine biosynthesis pathway. Both enzymes are homologous to DUS, their crystallographic structures have been determined and the mechanism of action is known. The availability of these structures prompted us to use comparative modeling to elucidate the structure of the DUS active site, its interactions with the substrate, and the mechanism of action.
Studies on DUS enzymes may have an applied edge, as it has been found that a human member of the DUS family, hDus2, has been implicated in cancer. In particular, the Nakamura group discovered that silencing the HDUS2 gene decreases the abundance of D residues in tRNA molecules, reduces the effectiveness of translation, and in consequence blocks the growth of cancer cells . This finding suggests that hDus2 may be a target for anticancer therapy. In particular, selective inhibitors able to block the DUS activity of hDus2 by interactions with its active site or tRNA binding site may be interesting leads for new anti-cancer drugs.
At the time of the original submission of this article (21st June 2011) the structure of DUS enzymes was unknown. To elucidate sequence-structure-function relationships in the DUS family we carried out a comprehensive bioinformatic analysis, including sequence searches, clustering and phylogenetic analyses, and protein structure prediction. We have also built a homology model of hDus2 and predicted enzyme-substrate interactions. Based on these analyses, we identified a putative active site and substrate-binding residues and proposed the mechanisms of DUS activity and its potential inhibition. While the manuscript was under review, a crystal structure of another representative of the DUS family has been determined in complex with tRNA , giving us an opportunity to validate theoretical predictions.
Results and discussion
Sequence database searches and retrieval of members of the DUS family
To identify a complete set of DUS sequences, we used full-length sequences of experimentally characterized DUS from E. coli and S. cerevisiae to carry out exhaustive PSI-BLAST  searches of the (nr) database (until reaching convergence) and retrieved all sequences reported with e-value better than 1e-25. We removed identical proteins retrieved in different searches and, as a result, we obtained a set of about 11000 sequences.
Subdivision of DUS sequences into closely related groups
Clustering of all full-length DUS sequences was performed based on their pair-wise BLAST similarity scores, using CLANS . We had experimentally found that for this particular dataset the P-value threshold of 1e-7 produced qualitatively best results (more stringent values caused disconnection of the most diverged families, while more permissive values caused over-compaction of the whole dataset into a single cluster with only a few outliers). This clustering revealed that from all retrieved sequences only members of COG0042 and KOG2333, KOG2334 and KOG2335 grouped together with the genuine DUS proteins. All other gathered sequences showed only slight similarity to DUS and while they should be considered homologous, they are likely to exhibit different enzymatic activities, possibly based on a generally similar mechanism. Importantly, the clustering confirmed sequence similarity of KOG1799 (dihydropirimidine dehydrogenase) and KOG1436 (dihydroorotate dehydrogenase) to DUS (data not shown). This analysis has also confirmed that the TM0096 protein from Thermatoga maritima, annotated as “a putative flavin oxidoreductase”, whose crystal structure has been solved in complex with a flavin mononucleotide by the Structural GenomiX consortium (Protein Data Bank ID - 1vhn) , is indeed a member of the DUS family.
Based on the results of preliminary clustering, we extracted members of individual subfamilies, calculated family-specific multiple sequence alignments using ClustalX  and adjusted them manually (as described in Methods) to remove truncated sequences and redundant, nearly identical versions of the same protein, and to improve the placement of insertions and deletions. We have submitted these alignments, as well as individual sequences of representatives from E. coli and S. cerevisiae to the GeneSilico metaserver , to predict a number of structural and functional features, such as: secondary structure, intrinsic disorder, RNA-binding residues, and to identify matches to related proteins with experimentally determined structures, in particular to TM0096.
A multiple sequence alignment of representative DUS sequences
The multiple sequence alignment also reveals residues that are specific to individual subfamilies, and hence may contribute to differences in substrate specificity between subfamilies. For example the eukaryotic Dus1 and Dus4 subfamilies possess common conserved residues (T/S)PM (see positions 50–52 in Additional file 1), which correspond to a very similar motif (T/S)EM present in all bacterial DUS. This region is not conserved in archaeal Dus, Dus2 and Dus3 subfamilies. DusA members possess a subfamily-specific conserved insertion GLSPKENREIPPLD (positions 292 to 336 in Additional file 1). Sequences of Dus3 members are distinguished by a highly conserved motif LTTVGNPPFRR (positions 22–32 in Additional file 1). Further, a common feature of Dus4, DusA, DusC and archaeal Dus is the lack of a conserved GA motif (positions 37–38 in Additional file 1). There are also several regions in the alignment corresponding to fragments with high sequence variability that are enriched in predicted intrinsic disorder and in predicted RNA-binding residues (e.g. positions 55–115, 260–270 in Additional file 1). While the presence or absence or length of these regions may be characteristic for individual subfamilies, their sequences are not highly conserved, which suggests that they may be responsible for sequence-nonspecific interactions with the RNA (e.g. via the backbone). Such regions may for instance influence the subfamily-specific recognition of entire regions in the RNA substrates and restrict the action of enzymes to particular loops in the RNA. Such regions are however typical for eukaryotic sequences, and are relatively scarce in bacterial subfamilies. Nonetheless, we predict that it is not the variability of residues in conserved motifs, but rather the length and sequence composition of loops immediately following the conserved motifs (and shaping the surface around the active site) that may contribute most to the substrate specificity of different DUS members. This prediction remains to be tested experimentally.
Domain architecture of DUS members
DusB and DusC consist only of the catalytic DUS domain. However, many members of the DUS family possess additional domains or extensions fused to N- or C-termini of the catalytic domain. Dus1p, Dus2p, Dus4p and DusA have additional disordered regions in at least one of the termini. Dus3p contains two long disordered fragments and additionally one zinc finger domain in the N-terminus. Human Dus2 (hDus2) contains a dsRBD (double stranded RNA-binding motif domain) in the C-terminus. The dsRBD domain is also present in Dus2 from other animals, including mammals (e.g. cow, mouse, rat), amphibia (e.g. western clawed frog), flatworms (e.g. blood fluke), nematodes (e.g. filarial nematode worm), and insects (e.g. african malaria mosquito, fruit fly). On the other hand it is absent in Dus2 from fungi and plants. It will be interesting to determine whether the presence of the dsRBD domain may influences the range of substrates modified by enzymes in the Dus2 family (e.g. in the human vs the yeast enzyme).
Eukaryotic proteins often differ from their bacterial counterparts by possessing additional regions that exhibit intrinsic structural disorder, and according to our predictions, DUS proteins are no exception to this rule. Among DUS from E. coli, only DusA contains a short N-terminal extension that is predicted to be disordered. Dus1p from S. cerevisiae has two regions predicted to be disordered, spanning residues 1–15 and 324–423. Dus2p possesses a C-terminal disordered region spanning residues 330–382. Dus3p has very long disordered fragments in both termini, spanning residues 1–280 and 601–668. Both termini of Dus4p are predicted to be disordered, spanning residues 1–36 and 300–367. Disordered regions are commonly involved in interactions with other molecules, and in DUS they may play a role in tRNA binding and substrate specificity. Computational prediction of RNA binding sites for all DUS representatives from S. cerevisiae and E. coli (See Methods for details) indicate that disordered regions in Dus1p (residues 356–423), Dus2p (345–382), Dus3p (1–30, 65–90, 145–165, and 655–668), Dus4p (1–36, 335–367), and DusA (1–28) are rich in RNA binding residues.
Phylogenetic distribution of members of the DUS family
Phylogenetic tree of DUS family
Based on the sequence alignment of the structurally conserved (and therefore reliably alignable) region of DUS sequences (Additional file 1), we attempted to infer the phylogenetic tree of the family. Traditional methods for phylogenetic reconstruction based on sequence data, including neighbor-joining (NJ), minimum evolution (ME), and maximum parsimony (MP) failed to produce a confident tree with well-resolved branches (data not shown). This failure was most likely caused by high sequence divergence and uneven rates of evolution between and within different DUS subfamilies. On the other hand, the number of sequences in the alignment is relatively large for the computationally demanding maximum likelihood (ML) method. Therefore, we decided to use the Bayesian approach, which combines relative reliability of ML with fast scanning of the parameter landscape by the Markov chain Monte Carlo (MCMC) approach implemented in MrBayes3 program . The family trees were calculated for an alignment comprising only 107 representatives selected from gathered COG and KOG sequences. Only well aligned regions were used for the calculations.
A speculative scenario of the evolutionary history of the DUS family
The order of branching of eukaryotic DUS subfamilies suggests that Dus3 may be the ancestral eukaryotic enzyme, from which others have been derived by gene duplication, starting with Dus2, and then Dus1 and Dus4. The absence of Dus4 in Viridiplantae can be explained by gene loss. Dus2p and Dus3p are able to modify only a single specific position in tRNA molecules (20 and 47 respectively), in contrast to Dus4p and Dus1p that can act on two neighboring positions in the same tRNA region. This suggests that the specificity towards one target residue may be evolutionally older than the ability to modify several residues in one region, at least in yeast. Thus, the evolution of these enzymes has followed a trend in which new members of the family showed relaxed specificity compared to their ancestors.
The phylogeny of bacterial DUS families is more complex. DusB appears in nearly all Bacteria, which suggest that it might have been a bacterial DUS ancestor. The ancestor of the DusC subfamily probably appeared as a result of DusB duplication shortly after the divergence of main Proteobacteria groups, in a branch leading to β and γ Proteobacteria. The DusA subfamily contains proteins highly similar to each other, but relatively distinct from other subfamilies. We suggest that DusA appeared by duplication of DusB in an ancestor of Proteobacteria. A single representative of DusA in Cyanobacteria was probably transferred there horizontally.
An interesting feature in the DusB subfamily is the appearance of additional paralogous copies of DUS (Figure 4) in Bacillus subtilis, Clostridium acetobutylicum and Nostoc sp. Functions of these duplicated genes are unknown. It is possible that they acquired specificities for new positions or that they specialized in modifying the same position in a subset of substrates. Unfortunately, little is known about the position of D residues in tRNAs from these organisms. Definitely, the determination of tRNA sequences from additional organisms would dramatically help in the evolutionary studies of RNA modification enzymes, in particular with respect to their sequence specificities.
The DusA subfamily has members not only in Bacteria but also in plants (Arabidopsis thaliana; two proteins: At3g63510 and At5g47970) and fungi (Encephalitozoon cuniculi). These genes are most likely products of horizontal gene transfers from endocellular bacterial symbionts: from the alpha-proteobacterial ancestor of a mitochondrion to an ancestor of all contemporary Eukaryota, and from the cyanobacterial ancestor of a chloroplast to the ancestor of green plants. The absence of DusA-like members in many eukaryotic genomes can be simply explained by the gene loss following the initial transfer.
For At3g63510 sequence localization predictions according to the methods used (Sherlock , WolfPSORT , Plant-Ploc , Protein-Powler , TargetP  and Cello ) indicated chloroplast localization, which confirmed our hypothesis that this protein has a cyanobacterial origin. Predictions for At5g47970 sequence, however, failed to identify one preferred localization, and predictions for the DusA-like protein from Encephalitozoon cuniculi indicated cytoplasmic localization, hence the site of action of these proteins remains to be characterized experimentally.
Structural model for the hDus2 protein
In the absence of experimentally determined structure for functionally characterized DUS members (as of 21st June 2011), we constructed a comparative model of a human DUS2 enzyme, to provide a structural platform for the investigation of sequence-structure-function relationships in this family. However, while this article was under review, a crystal structure of dihydrouridine synthase from Thermus thermophilus (TthDus) and its complex with tRNA have been determined (PDB code 3b0v) . Hence, we present only a short description of structural model of human Dus in the Results section, and we compare our structural and functional predictions with the independently obtained experimental data in the Discussion section.
We mapped information about conserved residues from the multiple sequence alignment of hDus2 homologs containing the DSRM domain onto the structure of hDus2 model (Figure 7C) and predicted that the substrate-binding site of hDus2 is most likely formed by residues that are highly conserved among all DUS family members, which correspond to N113, C116, K155, R157, H183, and R185 in hDus2. In the dsRBD the conserved residues included Y388, T390, D396, R397, F399, and S401, and the fragment with residues 419–429. Such conservation suggests that the dsRBD may be functionally important, and the most obvious hypothesis is that its function is related to RNA binding. At present it is not clear how the dsRBD does participate in tRNA binding and recognition by Dus2 members, in particular compared to those members of the Dus2 family that lack the dsRBD.
Prediction of ligand binding
As mentioned in the introduction, due to the sequence and structure similarity between DHODs, DHPDHs, and TM0096, it has been suggested that these enzymes may use a similar mechanism of catalysis . The DHPDH mechanism of action has been studied in detail . Pyrimidine binding triggers a conformational change of a flexible active-site loop (residues 666–683, contains conserved residues Leu-Asn-Leu-Ser-Cys-Pro) in the barrel domain, resulting in the placement of a catalytically indispensable cysteine residue (C671) close to the bound substrate. The flexible loop closure is an absolute prerequisite for the catalytic activity, as it not only excludes the surrounding solvent from the active site, but most importantly, also places C671 at the location required to enable the proton transfer to the pirimidine C5 atom. From FMN N5, a hydride is transferred to the uracil C6, and the C5 atom receives a proton from the C671 SH group.
The availability of the TthDus-tRNA complex allowed us to test also the results of the ligand binding prediction. Indeed, as we predicted the position, orientation, and conformation of ligands in DHODH and DHPDH are similar to those of the target nucleotide U20 in the crystal structure of the TthDus-tRNA complex. Thus, our prediction of the protein-ligand complex structure and the mechanism of action of DUS enzymes can be regarded as highly accurate.
Structural comparison of DUS representatives reveals differences that may be responsible for different substrate specificities
An important question in the comparative analysis of enzymes is the identification of the molecular basis of their different substrate specificities. For DUS, such studies are hampered by the paucity of information about substrate specificities; in particular it is not known whether substrate specificity is conserved within subfamilies. Nonetheless, the comparison of predicted DUS structures, in particular those for members from yeast, together with the analysis of subfamily-specific conserved residues and prediction of RNA-binding residues has revealed potential specificity-determining elements.
In order to identify structural differences between DUS enzymes that may be responsible for their different substrate specificities, we have built homology models for all seven representatives of DUS from yeast (Dus1p, Dus2p, Dus3p, and Dus4p), and E. coli (DusA, DusB, and DusC), using the same methodology as for modeling of hDus2. Superposition of these models with annotated conserved and predicted RNA-binding residues (Additional file 3) revealed that while the catalytic cores of these enzymes are very similar, they exhibit variations in internal loops surrounding the active site, and they have very different N and C termini. The terminal extensions are particularly different between eukaryotic enzymes.
According to predictions RNA binding sites might be located in nine regions: I32-R45 (for the purpose of the article called binding region 1, BR1), H62-R73 (BR2), R87-D97 (BR3), D116-L135 (BR4), N147-R162 (BR5), T186-W204 (BR6), P215-Y228 (BR7), D242-P254 (BR8), and A291-F300 (BR9) (numbering of residues as in Dus1p). Dus1p and Dus4p and all bacterial DUS possess predicted RNA-binding residues in all these regions (BR1-BR9) and they cover large part of the protein surfaces. In both Dus1p and Dus4p, predicted RNA-binding surfaces are of a triangular shape with the top vertex formed by a loop (62–73 in Dus1p). In Dus2p and Dus3p, regions corresponding to BR2 and BR3 are not predicted as involved in RNA binding. In Dus3p, BR1 and BR7 are also devoid of predicted RNA-binding residues. In addition, the yeast DUS exhibit different terminal extensions predicted both to be largely disordered (at least in the absence of the RNA substrate) and to contain RNA-binding residues. We predict that it is not the variability of residues in conserved motifs, but rather the length and sequence composition of loops immediately following the conserved motifs (and shaping the surface around the active site) that may contribute most to the substrate specificity of different DUS members. In addition, the highly divergent termini are likely to be involved in sequence-nonspecific interactions with the RNA substrates in such a way that they facilitate the binding of tRNA in some orientations (e.g. by electrostatic interactions) and/or prevent others (e.g. by steric exclusion). These predictions remain to be tested experimentally.
Evolution of DUS
We have carried out extensive bioinformatics analyses aimed at comprehensive classification of the DUS family and making functional predictions for the previously uncharacterized proteins. Our results revealed all members of the DUS family, which can be subdivided into 8 subfamilies (DusA, DusB, DusC, Dus1, Dus2, Dus3, Dus4 and archaeal Dus). The results of DUS family clustering combined with the phylogenetic analysis suggest that members of KOG2335 encode proteins from three different subfamilies: DusA, Dus1 and Dus4 (3, 7, and 5 members, respectively). Similarly, COG0042 consists of members from four subfamilies: DusA, archaeal Dus, DusC, and DusB (18, 6, 9, and 41 members, respectively). Furthermore, members of KOG2333 that have been annotated as “uncharacterized conserved proteins of unknown function” belong to the Dus3 subfamily, whereas members of KOG2334 correspond to the Dus2 family. Additionally, based on clustering, a new subfamily of archaeal Dus have been classified which should be from now on included in DUS classification and analysis. The results of our analysis suggest that the groups defined earlier in the COG/KOG database may require revision. We have also analyzed the phylogenetic distributions of DUS family members and inferred their evolutionary tree. The most important conclusion of this part of our work is that LUCA possessed one ancestral DUS enzyme, which underwent independent duplications in Eukaryota and Bacteria, but not in Archaea.
The refined grouping of DUS enzymes into orthologous and paralogous branches provides a framework to study the functional differences among these proteins, in particular their different substrate specificities. Our results will enable easier classification of new DUS members identified in the future.
Validation of the predictions in the light of the independently determined crystal structure
Proposed mechanism of DUS inhibition
The human hDus2 enzyme has been suggested as a possible target for cancer therapy. The knowledge of DUS-ligand interactions and of the catalytic mechanism may help in the development of inhibitors of hDus2. Our homology model of hDus2 enzyme may aid in the efforts towards development of inhibitors for this enzyme. It must be emphasized that DUS-specific inhibitors should be inactive against DHPDH and other related enzymes with a similar structure and mechanism of action in order to reduce their toxicity against human cells. It was shown that 5-iodouracil (5 IU) is an inhibitor of DHPDH activity . NADPH-dependent reduction of 5 IU leads to the formation of 5-iodo-5,6-dihydrouracil, which is released from the active site, but then re-binds, and causes covalent modification of C671, which ultimately blocks the enzyme activity . Derivatives of 5 IU, in particular variants with 5′ and/or 3′ extensions that would prevent binding to DHPDH and other related human enzymes and would e.g. mimic phosphate groups and increase the affinity towards DUS, might be thus potential candidates for inhibitors against hDus2.
Sequence database searches, multiple alignments and comparison of sequence-structure profile HMMs
A set of known members of the DUS family (DusA, NCBI Gene Identification [GI] number: 85676801, DusB GI: 85676051, DusC GI: 85675254 from Escherichia coli and Dus1p GI: 6323560, Dus2p GI: 6324342, hDus2 GI: 8923374, Dus3p GI: 85666121, Dus4p GI: 6323437 from Saccharomyces cerevisiae) were used as queries in PSI-BLAST  searches of the non-redundant (nr) version of the NCBI sequence database (as of 2011) with the expectation (e) value threshold for the retrieval of related sequences set to 10-25.
To visualize pairwise similarities between and within protein of DUS homologs we used CLANS (CLuster ANalysis of Sequences), a Java program that applies version of the Fruchterman-Reingold graph layout algorithm . CLANS uses the P-values of high-scoring segment pairs (HSPs) obtained from an N x N BLAST search, to compute attractive and repulsive forces between each sequence pair in a user-defined dataset. 3D or 2D representation is achieved by having points corresponding to sequences randomly positioned in space. The points are then moved within this environment according to the force vectors resulting from all pairwise interactions, the forces are recalculated after the move, and the process is repeated to convergence.
All sequences classified as members of the DUS family were aligned using ClustalX . Incomplete sequences were discarded (if the deletion spanned > 30% of the alignment) or repaired using amino acid sequences predicted from the available DNA sequences of the corresponding genes. Manual adjustments were introduced into the alignment to preserve the continuity of secondary structure elements, either observed in crystal structures or predicted computationally (see Protein Structure Prediction).
The phylogenetic tree of the DUS family was inferred for 98 representative members of COG0042 and KOGs 2333, 2334 and 2335 using the alignment of a complete catalytic domain. A Minimum Evolution (ME) analysis carried out with MEGA 4  (with pairwise gaps deletion and either Dayhoff or JTT matrices) were not sufficient to infer a tree with pre-defined subfamilies grouped into monophyletic branches, thus the Bayesian analysis was attempted. Phylogenetic trees were calculated, using the JTT model of substitutions and pairwise deletions, with the initial tree calculated by the NJ method and the Closest Neighbor Search option set to level = 2. The stability of individual nodes was calculated using the bootstrap test (1000 replicates) and additionally confirmed by the interior branch test – in all trees, for all branches with bootstrap support >50%, the ITB support was equal or higher.
MrBayes MPI version 3.1.2  was used to carry out a Bayesian analysis of sequences data. All analyses were performed using the sequence alignment for family representatives, with a gamma distribution of substitution rates, using the default approximation of four rate classes for each. Preliminary runs with MrBayes using a mixture of model priors demonstrated conclusively that substitution rates from the WAG/BLOSUM matrices  provided best fit to the sequence data. Therefore, the WAG/BLOSUM model was used to provide substitution priors for the sequence partition of the data. A Metropolis-coupled Markov-chain Monte-Carlo analysis was performed with 2 million generations, two runs and eight chains (four per run). The Markov chain was sampled every 100 generations. Convergence of runs was confirmed by average standard split deviation factor that falls under the recommended value of 0.01. The final tree was obtained after removing the first 25% of samples.
Protein structure prediction
Protein sequences were submitted to the GeneSilico metaserver, which is a gateway to a large number of third-party methods that facilitates comparison and interpretation of predictions made by different algorithms . In particular, the metaserver was used for secondary structure prediction and for protein fold-recognition (i.e. alignment of protein sequence to proteins with experimentally determined structures that can be used as templates for modeling). Fold-recognition alignments reported by primary methods were compared, evaluated, and ranked by the PCONS method . PCONS score >1 in general indicates estimation that the protein fold has been correctly guessed by FR methods. However, lower scores do not necessarily exclude correct predictions, in particular for folds with strongly diverged members. In such cases, a good estimator of prediction quality is the number of occurrences of a given fold at the top positions of the PCONS ranking. PCONS has traditionally performed very well in all editions of CASP and other benchmarks .
The fold-recognition alignments of hDus2 and the top-scoring templates (PDB codes 1vhn and 1whn for N-terminal and C-terminal hDus2 domains, respectively) were used as a starting point for modeling of hDus2 tertiary structure comprising cycles of model building by Modeller , evaluation by MetaMQAP , realignment in poorly scored regions as long as manual alignment changes does not improve model quality. Uncertain regions (residues 1–11, 62–78, 253–266) were modeled de novo using ROSETTA  in the context of ‘frozen’ remainder of the hDus2 model. Briefly, fragment selection based on profile-profile and secondary structure comparison with the ROSETTA database was performed and 3 and 9 amino acids fragments lists were generated for the remodeled regions. The fragment assembly was performed with default options and medium level of side chains rotamers optimization. The set of 5000 preliminary models (decoys) of was clustered and representatives of 3 largest clusters were selected as the final structures for evaluation.
Protein model evaluation
For evaluation of models we used two Model Quality Assessment Programs (MQAPs): MetaMQAP  and PROQ . It must be emphasized that MQAP scores only predict the deviation of a model from the real structure (the real deviation can be calculated only by comparison to the real structure, which of course is not available). Thus, the scores reported in this work that indicate e.g. ‘very good models’, must be interpreted as estimations or predictions that our models are ‘very good’, and not as ultimate validation of the model quality. However, it should be mentioned that both PROQ and MetaMQAP performed very well in CASP and in independent benchmarks and can be regarded as robust predictors.
Prediction of intrinsically disordered residues
Predictions of intrinsically disordered residues were made using MetaDisorder (http://iimcb.genesilico.pl/metadisorder/; , a meta-method which combines the predictions of the following primary methods: DisEMBL , DISPROT(VSL2) , GlobPlot , IUPred , PDISORDER (SoftBerry, http://linux1.softberry.com/berry.phtml), POODLE-S , POODLE-L , PrDOS , Prosat (http://glaros.dtc.umn.edu/gkhome/node/456) and RONN .
Prediction of RNA-binding residues
Prediction of RNA-binding residues for protein sequences was made using a specialized meta-predictor  based on three sequence-based primary predictors that ranked highest in our tests (PiRaNhA , PPRInt , and BindN + ).
Authors would like to thank Kaja Milanowska, Dominik Kasprzak, Katarzyna H. Kamińska, Iga Korneta and Grzegorz Papaj for critical reading of the manuscript and Katarzyna H. Kamińska and Urszula Sobczyńska, Irina Tuszyńska and Tomasz Puton for their help during the project. This analysis was funded by the Foundation for Polish Science (grant TEAM/2009-4/2 to J.M.B. and grant POMOST_C/58 to J.M.K.). JMB has been also supported by the 7th Framework Programme of the European Commission (EC FP7, grant HEALTHPROT, contract number 229676). A.C. and J.M.K. have been supported by the Polish Ministry of Science and Higher Education (grants 0067/P01/2010/70 and 0083/IP1/2011/71 to AC and N N301 123138 to JMK). JMK has been a scholarship-holder of Adam Mickiewicz University Foundation for 2011.
- Sprinzl M, Horn C, Brown M, Ioudovitch A, Steinberg S: Compilation of tRNA sequences and sequences of tRNA genes. Nucleic Acids Res 1998, 26(1):148–153. 10.1093/nar/26.1.148PubMed CentralView ArticlePubMedGoogle Scholar
- Kowalak JA, Bruenger E, McCloskey JA: Posttranscriptional modification of the central loop of domain V in Escherichia coli 23 S ribosomal RNA. J Biol Chem 1995, 270(30):17758–17764. 10.1074/jbc.270.30.17758View ArticlePubMedGoogle Scholar
- Dalluge JJ, Hashizume T, Sopchik AE, McCloskey JA, Davis DR: Conformational flexibility in RNA: the role of dihydrouridine. Nucleic Acids Res 1996, 24(6):1073–1079. 10.1093/nar/24.6.1073PubMed CentralView ArticlePubMedGoogle Scholar
- Dalluge JJ, Hamamoto T, Horikoshi K, Morita RY, Stetter KO, McCloskey JA: Posttranscriptional modification of tRNA in psychrophilic bacteria. J Bacteriol 1997, 179(6):1918–1923.PubMed CentralPubMedGoogle Scholar
- Kato T, Daigo Y, Hayama S, Ishikawa N, Yamabuki T, Ito T, Miyamoto M, Kondo S, Nakamura Y: A novel human tRNA-dihydrouridine synthase involved in pulmonary carcinogenesis. Cancer Res 2005, 65(13):5638–5646. 10.1158/0008-5472.CAN-05-0600View ArticlePubMedGoogle Scholar
- Bishop AC, Xu J, Johnson RC, Schimmel P, de Crecy-Lagard V: Identification of the tRNA-dihydrouridine synthase family. J Biol Chem 2002, 277(28):25090–25095. 10.1074/jbc.M203208200View ArticlePubMedGoogle Scholar
- Czerwoniec A, Dunin-Horkawicz S, Purta E, Kaminska KH, Kasprzak JM, Bujnicki JM, Grosjean H, Rother K: MODOMICS: a database of RNA modification pathways. 2008 update. Nucleic Acids Res 2009, 37(Database issue):D118-D121.PubMed CentralView ArticlePubMedGoogle Scholar
- Xing F, Hiley SL, Hughes TR, Phizicky EM: The specificities of four yeast dihydrouridine synthases for cytoplasmic tRNAs. J Biol Chem 2004, 279(17):17850–17860. 10.1074/jbc.M401221200View ArticlePubMedGoogle Scholar
- Xing F, Martzen MR, Phizicky EM: A conserved family of Saccharomyces cerevisiae synthases effects dihydrouridine modification of tRNA. RNA 2002, 8(3):370–381. 10.1017/S1355838202029825PubMed CentralView ArticlePubMedGoogle Scholar
- Savage DF, de Crecy-Lagard V, Bishop AC: Molecular determinants of dihydrouridine synthase activity. FEBS Lett 2006, 580(22):5198–5202. 10.1016/j.febslet.2006.08.062View ArticlePubMedGoogle Scholar
- Dobritzsch D, Ricagno S, Schneider G, Schnackerz KD, Lindqvist Y: Crystal structure of the productive ternary complex of dihydropyrimidine dehydrogenase with NADPH and 5-iodouracil. Implications for mechanism of inhibition and electron transfer. J Biol Chem 2002, 277(15):13155–13166. 10.1074/jbc.M111877200View ArticlePubMedGoogle Scholar
- Rowland P, Bjornberg O, Nielsen FS, Jensen KF, Larsen S: The crystal structure of Lactococcus lactis dihydroorotate dehydrogenase A complexed with the enzyme reaction product throws light on its enzymatic function. Protein Sci 1998, 7(6):1269–1279. 10.1002/pro.5560070601PubMed CentralView ArticlePubMedGoogle Scholar
- Yu F, Tanaka Y, Yamashita K, Suzuki T, Nakamura A, Hirano N, Suzuki T, Yao M, Tanaka I: Molecular basis of dihydrouridine formation on tRNA. Proc Natl Acad Sci U S A 2011, 108(49):19593–19598. 10.1073/pnas.1112352108PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Frickey T, Lupas A: CLANS: a Java application for visualizing protein families based on pairwise similarity. Bioinformatics 2004, 20(18):3702–3704. 10.1093/bioinformatics/bth444View ArticlePubMedGoogle Scholar
- Park F, Gajiwala K, Noland B, Wu L, He D, Molinari J, Loomis K, Pagarigan B, Kearins P, Christopher J, et al.: The 1.59 A resolution crystal structure of TM0096, a flavin mononucleotide binding protein from Thermotoga maritima. Proteins 2004, 55(3):772–774. 10.1002/prot.20086View ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. NucleicAcids Res 1997, 25(24):4876–4882. 10.1093/nar/25.24.4876View ArticleGoogle Scholar
- Kurowski MA, Bujnicki JM: GeneSilico protein structure prediction meta-server. Nucleic Acids Res 2003, 31(13):3305–3307. 10.1093/nar/gkg557PubMed CentralView ArticlePubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, et al.: The COG database: an updated version includes eukaryotes. BMC Bioinforma 2003, 4(1):41. 10.1186/1471-2105-4-41View ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19(12):1572–1574. 10.1093/bioinformatics/btg180View ArticlePubMedGoogle Scholar
- Shatkay H, Hoglund A, Brady S, Blum T, Donnes P, Kohlbacher O: SherLoc: high-accuracy prediction of protein subcellular localization by integrating text and protein sequence data. Bioinformatics 2007, 23(11):1410–1417. 10.1093/bioinformatics/btm115View ArticlePubMedGoogle Scholar
- Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K: WoLF PSORT: protein localization predictor. Nucleic Acids Res 2007, 35(Web Server issue):W585-W587.PubMed CentralView ArticlePubMedGoogle Scholar
- Chou KC, Shen HB: Euk-mPLoc: a fusion classifier for large-scale eukaryotic protein subcellular location prediction by incorporating multiple sites. J Proteome Res 2007, 6(5):1728–1734. 10.1021/pr060635iView ArticlePubMedGoogle Scholar
- Hawkins J, Davis L, Boden M: Predicting nuclear localization. J Proteome Res 2007, 6(4):1402–1409. 10.1021/pr060564nView ArticlePubMedGoogle Scholar
- Emanuelsson O, Brunak S, von Heijne G, Nielsen H: Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc 2007, 2(4):953–971. 10.1038/nprot.2007.131View ArticlePubMedGoogle Scholar
- Yu CS, Lin CJ, Hwang JK: Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci 2004, 13(5):1402–1406. 10.1110/ps.03479604PubMed CentralView ArticlePubMedGoogle Scholar
- Glaser F, Pupko T, Paz I, Bell RE, Bechor-Shental D, Martz E, Ben-Tal N: ConSurf: identification of functional regions in proteins by surface-mapping of phylogenetic information. Bioinformatics 2003, 19(1):163–164. 10.1093/bioinformatics/19.1.163View ArticlePubMedGoogle Scholar
- Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA: Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci U S A 2001, 98(18):10037–10041. 10.1073/pnas.181342398PubMed CentralView ArticlePubMedGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. DeLano Scientific, Palo Alto, CA, USA; 2002.Google Scholar
- Rider LW, Ottosen MB, Gattis SG, Palfey BA: Mechanism of dihydrouridine synthase 2 from yeast and the importance of modifications for efficient tRNA reduction. J Biol Chem 2009, 284(16):10324–10333. 10.1074/jbc.M806137200PubMed CentralView ArticlePubMedGoogle Scholar
- Porter DJ, Chestnut WG, Taylor LC, Merrill BM, Spector T: Inactivation of dihydropyrimidine dehydrogenase by 5-iodouracil. J Biol Chem 1991, 266(30):19988–19994.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol 2007, 24(8):1596–1599. 10.1093/molbev/msm092View ArticlePubMedGoogle Scholar
- Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol 2001, 18(5):691–699. 10.1093/oxfordjournals.molbev.a003851View ArticlePubMedGoogle Scholar
- Lundstrom J, Rychlewski L, Bujnicki J, Elofsson A: Pcons: a neural-network-based consensus predictor that improves fold recognition. Protein Sci 2001, 10(11):2354–2362.PubMed CentralView ArticlePubMedGoogle Scholar
- Wallner B, Elofsson A: Prediction of global and local model quality in CASP7 using Pcons and ProQ. Proteins 2007, 69(Suppl 8):184–193.View ArticlePubMedGoogle Scholar
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol 1993, 234(3):779–815. 10.1006/jmbi.1993.1626View ArticlePubMedGoogle Scholar
- Pawlowski M, Gajda MJ, Matlak R, Bujnicki JM: MetaMQAP: a meta-server for the quality assessment of protein models. BMC Bioinforma 2008, 9(1):403. 10.1186/1471-2105-9-403View ArticleGoogle Scholar
- Misura KM, Chivian D, Rohl CA, Kim DE, Baker D: Physically realistic homology models built with ROSETTA can be more accurate than their templates. Proc Natl Acad Sci U S A 2006, 103(14):5361–5366. 10.1073/pnas.0509355103PubMed CentralView ArticlePubMedGoogle Scholar
- Wallner B, Fang H, Elofsson A: Automatic consensus-based fold recognition using Pcons, ProQ, and Pmodeller. Proteins 2003, 53(Suppl 6):534–541.View ArticlePubMedGoogle Scholar
- Kozlowski LP, Bujnicki JM: MetaDisorder: a meta-server for the prediction of intrinsic disorder in proteins. BMC Bioinformatics 2012, 13(13):111–1459.PubMed CentralView ArticlePubMedGoogle Scholar
- Linding R, Jensen LJ, Diella F, Bork P, Gibson TJ, Russell RB: Protein disorder prediction: implications for structural proteomics. Structure 2003, 11(11):1453–1459. 10.1016/j.str.2003.10.002View ArticlePubMedGoogle Scholar
- Dunker AK, Brown CJ, Lawson JD, Iakoucheva LM, Obradovic Z: Intrinsic disorder and protein function. Biochemistry 2002, 41(21):6573–6582. 10.1021/bi012159+View ArticlePubMedGoogle Scholar
- Obradovic Z, Peng K, Vucetic S, Radivojac P, Dunker AK: Exploiting heterogeneous sequence properties improves prediction of protein disorder. Proteins 2005, 61(Suppl 7):176–182.View ArticlePubMedGoogle Scholar
- Dosztanyi Z, Csizmok V, Tompa P, Simon I: IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics 2005, 21(16):3433–3434. 10.1093/bioinformatics/bti541View ArticlePubMedGoogle Scholar
- Shimizu K, Hirose S, Noguchi T: POODLE-S: web application for predicting protein disorder by using physicochemical features and reduced amino acid set of a position-specific scoring matrix. Bioinformatics 2007, 23(17):2337–2338. 10.1093/bioinformatics/btm330View ArticlePubMedGoogle Scholar
- Hirose S, Shimizu K, Kanai S, Kuroda Y, Noguchi T: POODLE-L: a two-level SVM prediction system for reliably predicting long disordered regions. Bioinformatics 2007, 23(16):2046–2053. 10.1093/bioinformatics/btm302View ArticlePubMedGoogle Scholar
- Ishida T, Kinoshita K: PrDOS: prediction of disordered protein regions from amino acid sequence. Nucleic Acids Res 2007, 35(Web Server issue):W460-W464.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang ZR, Thomson R, McNeil P, Esnouf RM: RONN: the bio-basis function neural network technique applied to the detection of natively disordered regions in proteins. Bioinformatics 2005, 21(16):3369–3376. 10.1093/bioinformatics/bti534View ArticlePubMedGoogle Scholar
- Puton T, Kozlowski L, Tuszynska I, Rother K, Bujnicki JM: Computational methods for prediction of protein-RNA interactions. J Struct Biol 2011, ():-. 10.1016/j.jsb.2011.10.001PubMedGoogle Scholar
- Murakami Y, Spriggs RV, Nakamura H, Jones S, et al.: PiRaNhA: a server for the computational prediction of RNA-binding residues in protein sequences. Nucleic Acids Res 2010, 38(Web Server issue):W412-W416.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar M: Gromiha MM. Prediction of RNA binding sites in a protein using SVM and PSSM profile. Proteins, Raghava GP; 2007.Google Scholar
- Wang L, Huang C, Yang MQ, Yang JY: BindN + for accurate prediction of DNA and RNA-binding residues from protein sequence features. BMC Syst Biol 2010, 4(Suppl 1):S3. 10.1186/1752-0509-4-S1-S3PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.