Computational modeling suggests dimerization of equine infectious anemia virus Rev is required for RNA binding
Retrovirology volume 11, Article number: 115 (2014)
The lentiviral Rev protein mediates nuclear export of intron-containing viral RNAs that encode structural proteins or serve as the viral genome. Following translation, HIV-1 Rev localizes to the nucleus and binds its cognate sequence, termed the Rev-responsive element (RRE), in incompletely spliced viral RNA. Rev subsequently multimerizes along the viral RNA and associates with the cellular Crm1 export machinery to translocate the RNA-protein complex to the cytoplasm. Equine infectious anemia virus (EIAV) Rev is functionally homologous to HIV-1 Rev, but shares very little sequence similarity and differs in domain organization. EIAV Rev also contains a bipartite RNA binding domain comprising two short arginine-rich motifs (designated ARM-1 and ARM-2) spaced 79 residues apart in the amino acid sequence. To gain insight into the topology of the bipartite RNA binding domain, a computational approach was used to model the tertiary structure of EIAV Rev.
The tertiary structure of EIAV Rev was modeled using several protein structure prediction and model quality assessment servers. Two types of structures were predicted: an elongated structure with an extended central alpha helix, and a globular structure with a central bundle of helices. Assessment of models on the basis of biophysical properties indicated they were of average quality. In almost all models, ARM-1 and ARM-2 were spatially separated by >15 Å, suggesting that they do not form a single RNA binding interface on the monomer. A highly conserved canonical coiled-coil motif was identified in the central region of EIAV Rev, suggesting that an RNA binding interface could be formed through dimerization of Rev and juxtaposition of ARM-1 and ARM-2. In support of this, purified Rev protein migrated as a dimer in Blue native gels, and mutation of a residue predicted to form a key coiled-coil contact disrupted dimerization and abrogated RNA binding. In contrast, mutation of residues outside the predicted coiled-coil interface had no effect on dimerization or RNA binding.
Our results suggest that EIAV Rev binding to the RRE requires dimerization via a coiled-coil motif to juxtapose two RNA binding motifs, ARM-1 and ARM-2.
The Rev protein of lentiviruses mediates nuclear export of singly spliced and unspliced viral RNA transcripts. HIV-1 Rev binds its target RNA, the Rev responsive element (RRE), as a monomer, then multimerizes along the RRE RNA before being shuttled to the cytoplasm through association with the Crm1 host export factor. Several well-characterized motifs mediate known functions of HIV-1 Rev: a nuclear localization signal (NLS), which overlaps an RNA-binding arginine-rich motif (ARM); a nuclear export signal; and a pair of oligomerization domains that flank the ARM (reviewed in ).
RNA recognition by HIV-1 Rev is mediated by a 17-residue long ARM that adopts an alpha-helical conformation and initially docks into the major groove of a highly structured region in the HIV-1 RRE, termed stem loop IIB (SLIIB) -. Biochemical and biophysical studies have revealed that HIV-1 Rev oligomerizes , and that monomeric, dimeric, and higher-order oligomeric forms associate with the RRE -. These and subsequent studies  have shown that HIV-1 Rev binding to the RRE is a stepwise process: initial binding of Rev to SLIIB acts as a nucleation event that drives further oligomerization of additional copies of Rev along the RRE (reviewed in ). Although, monomeric HIV-1 Rev has been shown to bind the RRE in gel shift, filter binding, and single molecule fluorescence spectroscopy assays ,,,, studies of the tertiary structure of HIV-1 Rev and the RRE suggest that the “fundamental building block” for RRE binding is a Rev dimer ,. Both dimeric (head-to-head) contacts and higher-order oligomeric (tail-to-tail) intermolecular contacts are critical for Rev-mediated RNA export ,-.
NMR studies of HIV-1 Rev revealed that the C-terminal half of the protein is intrinsically disordered ; however, crystal structures of the N-terminal half of HIV-1 Rev, including the ARM and oligomerization motifs, have provided valuable insights into the structural basis of RNA binding and multimerization ,,. In the crystal structures, the N-terminal half of the Rev monomer adopts a helix-loop-helix structure with hydrophobic patches on opposite surfaces. Hydrophobic patches on one surface contain residues that drive dimerization, whereas hydrophobic patches on the opposite surface contain residues that mediate oligomerization (reviewed in ,). Dimerization of HIV-1 Rev orients monomers in a ‘V’ shape with an angle of 120-140° and a distance of ~55 Å between the distal ends ,. Recent SAXS analysis  indicates that the HIV-1 RRE adopts an unusual topology resembling the letter ‘A’, with the ‘legs’ forming Rev binding tracks. The ‘legs’ are spaced ~55 Å apart and appear to match the distance between the ARMs in an HIV-1 Rev dimer ,,. Although, Rev monomers can bind the HIV-1 RRE, it is believed that the specific structural arrangement of Rev dimers combined with the complementary topology of the RRE dictates Rev-RRE binding specificity and aids recognition of cognate RNA substrate from among an abundant pool of host RNAs ,,,.
Equine infectious anemia virus (EIAV) Rev is functionally homologous to HIV-1 Rev, but shares very little sequence similarity and differs in domain organization (reviewed in ). We previously identified a bipartite RNA binding domain in EIAV Rev that contains two short arginine-rich motifs (designated ARM-1 and ARM-2) spaced 79 amino acids apart in the primary sequence . ARM-1 is located in the central region of the protein while ARM-2 resides at the C-terminus and also functions as an NLS. It is possible that ARM-1 and ARM-2 are in close proximity in the Rev monomer, forming a single RNA-binding interface. Alternately, ARM-1 and ARM-2 could each bind different sites on the RRE RNA. RNA footprinting and chemical modification experiments have shown that the RRE target of EIAV Rev contains two Rev binding regions (designated RBR-1 and RBR-2) that undergo conformational changes in the presence of EIAV Rev . RBR-1 encompasses the minimal RRE sequence, which overlaps a characterized exonic splicing enhancer -, while RBR-2 is necessary for high-affinity Rev binding in vitro .
Insight as to how the bipartite RNA binding domain interacts with the RRE target requires knowledge of the tertiary structure of EIAV Rev and relative positioning of ARM-1 and ARM-2 in the folded protein. Obtaining high-resolution structures of Rev proteins has proven very challenging due to the tendency of Rev to spontaneously aggregate into insoluble filaments in solution -. Computational modeling of the EIAV Rev structure has been challenging, also, because the amino acid sequence similarity between HIV-1 Rev and EIAV Rev is almost undetectable . Thus, it is not possible to simply use the available experimental structures of HIV-1 Rev as templates for homology modeling of non-primate Rev proteins. Recent progress in ab initio and threading methods for structure prediction, however, has provided a viable platform for modeling structures of proteins that have proven difficult to characterize experimentally -.
The first proposed structural model for EIAV Rev suggested that ARM-1 and ARM-2 are juxtaposed to form a single RNA binding interface on the monomer structure . In the present study, newer and more accurate structural modeling approaches were employed to predict the topology and relative orientation of ARM-1 and ARM-2 within the overall structure of EIAV Rev. Our results suggest that ARM-1 and ARM-2 do not form a single RNA binding interface within a single Rev monomer. Instead, our computational analyses, supported by experimental data, suggest that dimerization of Rev is a prerequisite for RNA binding. Thus, dimerization of EIAV Rev may be required to juxtapose ARMs from two Rev monomers so that they form a single functional RNA binding domain that recognizes the EIAV RRE.
Generation of Rev structural models
A computational approach was employed to model the tertiary structure of EIAV Rev and obtain insight into the topology of the bipartite RNA binding domain. Because EIAV Rev is highly variable in sequence and contains non-essential regions (reviewed in ), our analyses included deletion mutants as well as a divergent Rev variant (Figure 1). Rev165 is the full-length sequence of the R1 variant derived from the EIAVWyo2078 field isolate ; Rev135 contains an N-terminal deletion of R1 encompassing all of exon 1, while Rev∆HVR contains a 13-amino-acid deletion in the hypervariable region, located in the C-terminal half of the protein ,. Rev135 and Rev∆HVR are functionally equivalent to Rev165 in in vitro assays of nuclear export activity . RevFDD is from the fetal donkey dermal cell-adapted Chinese isolate EIAVFDD-10 . FDD Rev and Rev165 differ at 54 positions across the length of the protein (Figure 1B), demonstrating marked variation in primary amino acid sequence.
A total of 235 computational models were generated using several state-of-the-art protein structure prediction servers that implement different algorithms (Table 1). These include the QUARK server, which implements an ab initio algorithm ,; the ITASSER and LOMETS servers, which implement threading-based algorithms ,; and the PROTINFO server, which implements an homology modeling algorithm . Examination of the predicted models revealed significant variation in their overall shapes (Table 1): 52% of the Rev165 models had an elongated topology, 34% had a globular topology, and 14% were either unfolded and/or truncated. Models generated for Rev135, Rev∆HVR, and RevFDD also showed elongated and globular topologies, in proportions comparable to those for Rev165. The topology of models obtained depended, in part, on the method of structure prediction: the QUARK ab initio server predominately yielded elongated topologies, whereas homology servers generated exclusively globular topologies. A mixture of elongated and globular topologies was generated by the ITASSER and LOMETS threading servers.
Assessment and structural features of Rev models
The quality of generated models was assessed using the QMEAN and ProQ2 model quality assessment programs - (see Additional file 1). These programs evaluate the physicochemical and structural features of a given model by comparison with those of known experimental structures. Both programs generate a score in the range 0–1, with 1 signifying the highest possible quality score. The calculated QMEAN score for each model was plotted against its corresponding ProQ2 score and the distribution of model quality scores for each Rev sequence analyzed is shown in Figure 2. The majority of models have QMEAN and ProQ2 scores ~0.5, indicating most are of average quality. Although the elongated models generally scored higher than globular models, the overlap in quality scores precluded selection of a single preferred tertiary topology.
Key structural features in EIAV Rev elongated and globular models were identified by visual inspection using PyMol software . Figure 3A shows representative examples of the top ranking elongated and globular models for Rev165, Rev135, Rev∆HVR, and RevFDD. A distinguishing structural feature common to virtually all elongated models was an extended alpha helix in the central region of the protein. In globular models, this central helix was disrupted by a kink (Figure 3A, black arrows), resulting in a compact bundle of helices. In all models, ARM-1, ARM-2, and exon1 adopted alpha-helical conformations (Figure 3A). The NES formed a short alpha helix flanked on both sides by flexible loops, which usually formed a helix-turn-helix motif with the adjacent helix. ARM-1 was always positioned at the N-terminus of the central region. The other functional motifs were separated by flexible regions, and the positioning of these motifs relative to the central region was the major difference among the various models in both elongated and globular structures. The greatest variability was observed in the position of exon1 (Figure 3A), which is rich in solvent-exposed, hydrophilic residues.
Relative positioning of the bipartite RNA binding domain in Rev models
PyMol software  was used to inspect the relative positioning of ARM-1 and ARM-2 on the surface of the predicted structures. The distance separating the two closest atoms in ARM-1 and ARM-2 was calculated for each model (Additional file 1). In 204 of 235 generated models, ARM-1 and ARM-2 were separated by ≥15 Å on the monomer surface. In addition, ARM-1 and ARM-2 were positioned on opposite faces of the monomer in many of the top structures (Figure 3B). Although, electrostatic views show that the two ARMs could be bridged by a continuous stretch of positive charge in some models (Figure 3C), the bridging region consisted of positively charged residues from exon1, which can be deleted with no loss of Rev activity in vitro . These results strongly suggest that ARM-1 and ARM-2 do not form a single RNA binding interface in the Rev monomer.
A coiled-coil motif in EIAV Rev may promote dimerization
Given that ARM-1 and ARM-2 are not predicted to form a single RNA binding interface on the Rev monomer, two scenarios for RNA binding are possible: i) ARM-1 and ARM-2 form two distinct RNA binding interfaces, or ii) dimerization of Rev juxtaposes ARM-1 from one monomer with ARM-2 from a second monomer to form a single RNA binding interface. The latter scenario predicts that EIAV Rev dimerizes, and that dimerization is essential for RNA binding. Therefore, the primary sequence of EIAV Rev was computationally analyzed for oligomerization motifs -. Results identified a canonical coiled-coil motif, spanning residues 82–109, within an extended alpha helix predicted in the central region of Rev (Figure 4A). The predicted coiled-coil motif displayed characteristics typical for an oligomerization domain ,,,, with hydrophobic residues predominantly occupying ‘a’ and ‘d’ registers of the coiled-coil and charged residues preferentially occurring in ‘e’ and ‘g’ registers.
A helical wheel projection of the predicted coiled-coil motif (Figure 4B) shows that the ‘a’ and ‘d’ residues (Leu, Ile, Val, Ala) constitute the hydrophobic face of an amphipathic helix and are well positioned to mediate dimerization. Additionally, a bulky Trp residue is predicted to reside on the opposite side of the interhelical interface. Docking of predicted coiled-coil structures using the ClusPro server - resulted in formation of a head-to-tail dimer, with residues in the ‘a’ and ‘d’ registers forming an interhelical interface and the bulky Trp residue segregating to the opposite face, in a position where it could mediate further oligomerization (Figure 4C). Although both head-to-head and head-to-tail orientations were obtained by docking, the head-to-tail orientation resulted in a larger number of contacts between hydrophobic ‘a’ and ‘d’ residues and a more energetically favorable dimer structure. Fewer interactions between ‘a’ and ‘d’ residues of the coiled-coil were observed when docking full-length elongated structures, in either the head-to-head or head-to-tail orientation (not shown).
The coiled-coil motif is highly conserved among EIAV Rev variants
There is high degree of genetic variation in EIAV Rev sequences (reviewed in ), and it was of interest to examine conservation of residues in the predicted coiled-coil motif. Accordingly, 200 distinct Rev amino acid sequences encompassing phylogenetically diverse isolates were retrieved from GenBank (Additional file 2), aligned, and analyzed using the WebLogo server . These analyses revealed that a large number of residues in the predicted coiled-coil region are, in fact, invariant (Figure 4D). More importantly, residues in the ‘a’ and ‘d’ positions are either completely conserved, or were substituted only with similarly hydrophobic residues. The high degree of conservation suggests that the predicted coiled-coil motif contributes an essential function in Rev activity. In support of this, mutation of hydrophobic residues located in the predicted interhelical interface (L95D, L109D) abrogated Rev activity, whereas mutation of hydrophobic residues that lie outside the predicted interface (e.g., V112D) retained wild-type Rev activity ,.
Dimerization is required for RNA binding in EIAV Rev
Coiled-coil motifs generally mediate intermolecular interactions and, less frequently, intramolecular interactions -. In the context of our predicted EIAV Rev structures, intermolecular coiled-coil interactions would be favored by elongated structures whereas intramolecular coiled-coil interactions would occur in globular structures. In the elongated structures, interactions between coiled-coil motifs in two different Rev monomers could juxtapose ARM-1 and ARM-2 to form a single RNA binding interface in a dimeric structure. In the globular models, intramolecular interactions between the two smaller helices in the center of the monomer could be important for structural stability . To test the hypothesis that the Rev coiled-coil motif mediates dimerization, purified MPB-Rev fusion proteins containing mutations in the predicted coiled-coil interface were analyzed by Blue native PAGE (Figure 5). Because Rev aggregates readily in solution, samples were resuspended in Blue native PAGE loading buffer supplemented with 0.2% SDS. MBP-Rev165 and MBP-Rev135 samples migrated as monomeric, dimeric and higher oligomeric forms (Figure 5B). In contrast, only the monomeric form was present in samples of MBP-Rev145-165, which contains only the C-terminal 21 amino acids of Rev. Aspartic acid substitution of Leu 95, which is predicted to be critical for mediating intermolecular coiled-coil interactions, resulted in loss of dimerization, whereas alanine substitutions of non-interfacial residues in the coiled-coil motif (ERLE to AALA) had no effect on dimerization (Figure 5B). These data support the hypothesis that the EIAV Rev coiled-coil motif mediates intermolecular interactions between Rev monomers, resulting in formation of dimers.
To explore the importance of dimerization for RNA binding, we re-examined previous studies that mapped determinants of EIAV Rev required for RNA binding . In UV-crosslinking experiments, the L95D mutation, which abolishes dimerization, resulted in markedly reduced RNA binding activity, whereas the AALA mutant in which dimerization is not affected retained wild-type binding activity (compare Figures 5B and C). Thus, loss of Rev dimerization is correlated with loss of RNA binding activity. Furthermore, mutations within ARM-1 and ARM-2 that disrupted RNA binding did not affect dimerization (Figure 5B,C), indicating that dimerization and RNA binding are distinct and separable functions of Rev. Taken together, these results indicate that a coiled-coil motif mediates dimerization of EIAV Rev, and that dimerization is a prerequisite for Rev binding to the RRE.
The Rev protein of EIAV contains a bipartite RNA binding domain, containing two arginine-rich motifs, designated ARM-1 and ARM-2, which are separated by 79 residues in the amino acid sequence. In this study, computational models were generated and evaluated in an effort to determine the relative positioning of ARM-1 and ARM-2 on the tertiary structure of EIAV Rev. Two overall topologies for the Rev monomer were predicted: an elongated structure with an extended central alpha helix, and a globular structure with a kink in the central helix, resulting in a bundle of helices. In 204 of 235 generated models, ARM-1 and ARM-2 were well separated on the tertiary structure, strongly suggesting that a single RNA binding interface is not formed on the Rev monomer. A highly conserved coiled-coil motif was identified in the central region of EIAV Rev and was found to mediate dimerization of Rev monomers in vitro. Mutation of residues predicted to form key intermolecular coiled-coil contacts abolished dimerization and also disrupted RNA binding. In contrast, mutation of residues predicted to lie outside the coiled-coil interface had no effect on dimerization or RNA binding activity. Taken together, our results suggest that the EIAV Rev monomer adopts an elongated structure that dimerizes through intermolecular interactions mediated by a highly conserved coiled-coil motif in the central region of the protein. Dimerization is predicted to juxtapose ARM-1 from one monomer with ARM-2 from a second monomer to form a single RNA binding interface.
The central region of Rev is known to be sensitive to mutation , but a specific role for this region in the Rev nuclear export pathway has not been identified. The presence of a highly conserved coiled-coil motif in the central region suggests that it is required for intermolecular and/or intramolecular interactions essential for Rev activity. In elongated structural models, the coiled-coil is positioned to meditate intermolecular interactions required for dimerization and RNA binding; in the globular models, the coiled-coil motif would mediate intramolecular interactions that contribute to protein stability. Our data are most consistent with an elongated topology wherein the coiled-coil motif mediates formation of a Rev dimer. In this scenario, coiled-coil intermolecular interactions that stabilize the EIAV Rev dimer are maximized in an antiparallel orientation, suggesting that EIAV Rev binds RNA as a head-to-tail dimer. In support of this model, a series of trans-complementation experiments reported by Harris et al.,  showed that co-transfection of ARM-1 and ARM-2 mutants, each deficient for RNA export, restored Rev activity. In contrast, trans-complementation was abolished by mutation of residues that correspond to key contacts in the coiled-coil interface. In total, the computational models and experiments reported here, combined with previous experimental results, indicate that a coiled-coil motif in the central region of EIAV Rev mediates dimerization of Rev, which in turn, plays an essential role in RNA binding and Rev activity.
The predicted overall fold of EIAV Rev reported here shows both similarities and differences compared with the crystal structure of the HIV-1 Rev monomer ,. In both Rev proteins, the ARM motifs adopt an alpha-helical conformation. Oligomerization domains are found in both proteins and play an essential role in Rev function. The oligomerization domain of EIAV Rev contains the strong signature of a coiled-coil motif, which is required for dimerization and binding to the RRE. A corresponding canonical coiled-coil motif is not found in HIV-1 Rev. Instead, hydrophobic residues flanking the ARM mediate oligomerization of Rev on the RRE ,,. One difference between the two lentiviral Rev proteins is that dimerization is required for RNA binding of EIAV, but not HIV-1, Rev in vitro. In both cases, however, dimerization may be the biologically relevant configuration that determines RNA-binding specificity and formation of a functional nuclear export complex in vivo.
Our study highlights the value of employing computational methods to gain insight into structure-function relationships of Rev proteins, which have proven extremely difficult to characterize experimentally. In particular, recent advances in ab initio and threading based modeling has resulted in increased power and accuracy in predicting protein structure. Ab initio methods have the advantage of not requiring a structural template that shares sequence homology to that of the protein of interest; current ab initio methods can reliably predict tertiary structures of proteins ≤ 200 amino acids in length -,. Model quality assessment has also improved significantly in recent years and provides a quantitative measure of confidence in the quality of predicted protein structures ,. Due to the low level of sequence identity between EIAV and HIV-1 Rev and the lack of other homologous templates, the “average” scores of our predicted models were not unexpected. The quality score of a given model depends in part, on whether the overall fold of the model is consistent with predictions of secondary structure generated by independent methods -; therefore, models of average quality can yield useful information on general topology and spatial features of a protein. The elongated topology is most consistent with secondary structure predictions, in which the central region of EIAV Rev adopts an extended alpha helical conformation (Figure 4A). This explains, in part, why the elongated models generally scored higher than globular models, especially those generated by ab initio servers. Although we were unable to select a single topology based on computational predictions alone, both the globular and elongated models indicate that ARM-1 and ARM-2 do not form a single RNA binding interface, a finding that motivated the search for an oligomerization motif in EIAV Rev. It will be of interest to determine whether coiled-coil motifs are found in other retroviral Rev or Rev-like proteins where they may contribute to oligomerization and nuclear export activity.
This study provides computational and experimental data indicating that dimerization of EIAV Rev is required for RNA binding. Our results suggest dimerization is mediated by a coiled-coil motif in the central region of Rev. This work illustrates that computational modeling, combined with a molecular genetics approach, can be a valuable tool for interrogating the tertiary structure of Rev proteins and generating testable hypotheses regarding the mechanisms by which lentiviral Rev proteins recognize and bind their cognate RNA targets.
Generation of EIAV Rev structural models
EIAV Rev R1 [GenBank:AAG53100] was used as the reference amino acid sequence for generating full-length EIAV Rev165 structural models. R1 was isolated from a pony experimentally infected with EIAVWyo2078, a highly virulent strain of EIAV . Additional Rev sequence variants included: R1 Rev135, which lacks the first 30 amino acids encoded by exon1; R1 Rev∆HVR, in which the hypervariable region (residues 131–143) is deleted ,; and RevFDD , the full-length Rev sequence from the Chinese isolate EIAVFDD-10 [GenBank:ADK35837].
The QUARK, ITASSER, LOMETS, and PROTINFO protein structure prediction servers were used for automated modeling of Rev and are described in ,,,,,. Default settings were used for the QUARK, ITASSER, and LOMETS servers. The “generate comparative models” option was used for PROTINFO. For the ITASSER server, in addition to default settings, Rev was modeled using an HIV Rev crystal structure (PDB:3lph)  as the specified template with two different parameterized settings: i) the “specify template without an alignment” mode; and ii) the “specify template with alignment” mode. Pairwise alignment of R1 and HIV-1 Rev 3lph:A amino acid sequences was generated with the T-Coffee webserver . All models were manually inspected and models with an unfolded topology or those missing C-terminal residues encompassing ARM-2 were excluded from further analysis.
Quality assessment of EIAV Rev structural models
The QMEAN , and ProQ2  servers were used to evaluate models for consistency with known protein structural features. These are among the top performing model quality assessment servers, routinely outperforming other assessment programs in recent CASP competitions ,,,. To discriminate between high and low quality models, QMEAN uses a composite scoring function based on four geometrical features: i) local geometry, ii) long-range interactions, iii) all-atom potential, and iv) solvation energy of residues ,. The output score ranges from 0 to 1, where 1 is the highest score. The mean scores of high, medium and low quality models are 0.68, 0.58, and 0.40, respectively ,. ProQ2 predicts both local and global “correctness” of models using a support vector machine algorithm that considers the following features of a given model: i) atom-atom and residue-residue contacts, ii) solvent accessibility, iii) predicted secondary structure, iv) predicted surface area exposure, and v) evolutionary information . The quality of a model predicted by ProQ2 is consistent with predictions by QMEAN . All models generated in this study were evaluated with both servers, using default parameters.
Alignment of EIAV Rev protein sequences
Pairwise protein alignments were performed with the T-Coffee webserver , using the default settings of the T-Coffee mode. Multiple sequence alignments were performed with MacVector software, using the Gonnet substitution matrix with default settings ,.
Prediction of coiled-coil motifs
Coiled-coil motif prediction was performed using the COILS ,, PAIRCOIL , and CCHMMPROF  servers. For the COILS server, the following parameters were used: a 28-residue window width, the MTDIK matrix, and the 2.5 fold weighting of positions ‘a’ and ‘d’. For the PARCOIL server, a 28-residue window width and a p-score cut-off of 0.05 were used. Default settings were used for CCHMMPROF. The DrawCoil 1.0 server  was used to generate helical wheel representations of predicted coiled-coils.
Analysis of sequence conservation
Two hundred distinct EIAV Rev amino acid sequences from the US, Ireland, and China, were retrieved from the NCBI GenBank protein database (see Additional file 2). A multiple sequence alignment of the central region of Rev was generated and a sequence logo corresponding to the coiled-coil motif (a.a. 82–109) was derived using the WebLogo server . Sequence logos generated by WebLogo summarize the overall conservation of residues at each column position in a sequence alignment by depicting stacks of residues at each position: the height of each residue indicates its relative frequency. Relative frequencies are expressed in terms of information content, or bits, on the y-axis.
Prediction of protein secondary structures
Secondary structure predictions for Rev proteins were obtained from the PSIPRED , ITASSER , and QUARK , webservers and manually aligned to generate a consensus secondary structure.
The central region of Rev165 encompassing the predicted coiled-coiled motif (amino acids 82–109) was modeled with ITASSER. The ClusPro 2.0 docking server was used to generate dimeric structures, using default parameters -.
Expression and purification of EIAV Rev
MBP-Rev fusion proteins were cloned and expressed in E. coli strain Rosetta Gami in NZY media as described previously . Following expression, cells were pelleted and resuspended in lysis buffer containing 25 mM HEPES pH 7.5, 200 mM NaCl, 2 mM beta-mercaptoethanol (BME), supplemented with 2 mM phenylmethylsulfonyl fluoride (PMSF) and Roche cOmplete® protease inhibitor cocktail tablet, according to the manufacturer’s protocol. The suspension was incubated with 1 mg/ml lysozyme on ice for 20 min and subjected to 10 cycles of freeze-thaw and 20 cycles of sonication. The suspension was clarified by centrifugation and mixed by rocking with Ni-NTA beads equilibrated in 50 mM Tris pH 8.0, 2 M NaCl, 2 mM BME, 0.1% Tween-20, 10 mM imidazole. After overnight incubation at 4°C, resin was rinsed with 5 sample volumes of equilibration buffer, washed with 5 sample volumes of wash buffer (50 mM Tris pH 8.0, 250 mM NaCl, 2 mM BME, 10 mM imidazole), and MBP-Rev fusion proteins were eluted in 50 mM Tris pH 8.0, 250 mM NaCl, 2 mM BME, 250 mM imidazole. Eluted protein samples were dialyzed against 50 mM Tris pH 8.0, 200 mM NaCl, 2 mM BME, 10% glycerol. The purity of all proteins preparations was confirmed by SDS-PAGE analysis.
Blue native PAGE assay
Purified MBP-Rev protein samples were added to 6X Blue native sample loading buffer (12 mM EDTA, 120 mM NaCl, 120 mM Bis-Tris pH 7.0, 60% glycerol, and 0.5% Coomassie brilliant blue G-250 manufactured by Thermo Scientific, Waltham, MA) supplemented with 0.2% SDS. Samples were analyzed by electrophoresis in 8% Blue native polyacrylamide gels with 50 mM Bis-Tris pH 7.0 anode buffer and 50 mM Tricine, 15 mM Bis-Tris pH 7.0, 0.002% Coomassie brilliant blue G-250 cathode buffer.
RNA binding assays
UV-crosslinking RNA binding assays were described previously . Briefly, 2–4 μg purified MBP-Rev was incubated with 104 cpm of 32P-labeled EIAV RRE RNA in binding buffer (10 mM HEPES-KOH, pH 7.5, 100 mM KCL, 1 mM MgCl2, 0.5 mM EDTA, 1 mM dithiothreitol, 50 μg/ml E. coli tRNA and 10% glycerol) for 20 min at room temperature. Following incubation, samples were UV-irradiated with 3×105 μJ at 254 nm for 7 min, followed by treatment with 0.1 mg/ml RNase A at 37°C for 2 min. Samples were boiled in an equal volume of SDS for 5 min and separated in 12% SDS-PAGE in Tris-glycine buffer. Gels were fixed in 50% methanol-10% acetic acid, dried, and exposed to phosphorimager screens overnight. UV cross-linked complexes were detected using a PersonalFX scanner and Quality One software (Bio-Rad, Hercules, CA).
Pollard VW, Malim MH: The HIV-1 Rev protein. Annu Rev Microbiol. 1998, 52: 491-532. 10.1146/annurev.micro.52.1.491.
Tan R, Chen L, Buettner JA, Hudson D, Frankel AD: RNA recognition by an isolated alpha helix. Cell. 1993, 73: 1031-1040. 10.1016/0092-8674(93)90280-4.
Tan R, Frankel AD: Costabilization of peptide and RNA structure in an HIV Rev peptide-RRE complex. Biochemistry. 1994, 33: 14579-14585. 10.1021/bi00252a025.
Battiste JL, Mao H, Rao NS, Tan R, Muhandiram DR, Kay LE, Frankel AD, Williamson JR: Alpha helix-RNA major groove recognition in an HIV-1 rev peptide-RRE RNA complex. Science. 1996, 273: 1547-1551. 10.1126/science.273.5281.1547.
Jain C, Belasco JG: A structural model for the HIV-1 Rev-RRE complex deduced from altered-specificity rev variants isolated by a rapid genetic strategy. Cell. 1996, 87: 115-125. 10.1016/S0092-8674(00)81328-8.
Grate D, Wilson C: Role REVersal: understanding how RRE RNA binds its peptide ligand. Structure. 1997, 5: 7-11. 10.1016/S0969-2126(97)00161-5.
Wingfield PT, Stahl SJ, Payton MA, Venkatesan S, Misra M, Steven AC: HIV-1 Rev expressed in recombinant Escherichia coli: purification, polymerization, and conformational properties. Biochemistry. 1991, 30: 7527-7534. 10.1021/bi00244a023.
Cole JL, Gehman JD, Shafer JA, Kuo LC: Solution oligomerization of the rev protein of HIV-1: implications for function. Biochemistry. 1993, 32: 11769-11775. 10.1021/bi00095a004.
Daly TJ, Doten RC, Rennert P, Auer M, Jaksche H, Donner A, Fisk G, Rusche JR: Biochemical characterization of binding of multiple HIV-1 Rev monomeric proteins to the Rev responsive element. Biochemistry. 1993, 32: 10497-10505. 10.1021/bi00090a028.
Mann DA, Mikaelian I, Zemmel RW, Green SM, Lowe AD, Kimura T, Singh M, Butler PJ, Gait MJ, Karn J: A molecular rheostat. Co-operative rev binding to stem I of the rev-response element modulates human immunodeficiency virus type-1 late gene expression. J Mol Biol. 1994, 241: 193-207. 10.1006/jmbi.1994.1488.
Thomas SL, Oft M, Jaksche H, Casari G, Heger P, Dobrovnik M, Bevec D, Hauber J: Functional analysis of the human immunodeficiency virus type 1 Rev protein oligomerization interface. J Virol. 1998, 72: 2935-2944.
Jain C, Belasco JG: Structural model for the cooperative assembly of HIV-1 Rev multimers on the RRE as deduced from analysis of assembly-defective mutants. Mol Cell. 2001, 7: 603-614. 10.1016/S1097-2765(01)00207-6.
Daugherty MD, Booth DS, Jayaraman B, Cheng Y, Frankel AD: HIV Rev response element (RRE) directs assembly of the Rev homooligomer into discrete asymmetric complexes. Proc Natl Acad Sci U S A. 2010, 107: 12481-12486. 10.1073/pnas.1007022107.
Fernandes J, Jayaraman B, Frankel A: The HIV-1 Rev response element: an RNA scaffold that directs the cooperative assembly of a homo-oligomeric ribonucleoprotein complex. RNA Biol. 2012, 9: 6-11. 10.4161/rna.9.1.18178.
Daugherty MD, Liu B, Frankel AD: Structural basis for cooperative RNA binding and export complex assembly by HIV Rev. Nat Struct Mol Biol. 2010, 17: 1337-1342. 10.1038/nsmb.1902.
Pond SJ, Ridgeway WK, Robertson R, Wang J, Millar DP: HIV-1 Rev protein assembles on viral RNA one molecule at a time. Proc Natl Acad Sci U S A. 2009, 106: 1404-1408. 10.1073/pnas.0807388106.
Fang X, Wang J, O'Carroll IP, Mitchell M, Zuo X, Wang Y, Yu P, Liu Y, Rausch JW, Dyba MA, Kjems J, Schwieters CD, Seifert S, Winans RE, Watts NR, Stahl SJ, Wingfield PT, Byrd RA, Le Grice SFJ, Rein A, Wang Y: An unusual topological structure of the HIV-1 Rev response element. Cell. 2013, 155: 594-605. 10.1016/j.cell.2013.10.008.
Malim MH, Cullen BR: HIV-1 structural gene expression requires the binding of multiple Rev monomers to the viral RRE: implications for HIV-1 latency. Cell. 1991, 65: 241-248. 10.1016/0092-8674(91)90158-U.
Daugherty MD, D'Orso I, Frankel AD: A solution to limited genomic capacity: using adaptable binding surfaces to assemble the functional HIV Rev oligomer on RNA. Mol Cell. 2008, 31: 824-834. 10.1016/j.molcel.2008.07.016.
Edgcomb SP, Aschrafi A, Kompfner E, Williamson JR, Gerace L, Hennig M: Protein structure and oligomerization are important for the formation of export-competent HIV-1 Rev-RRE complexes. Protein Sci. 2008, 17: 420-430. 10.1110/ps.073246608.
Hoffmann D, Schwarck D, Banning C, Brenner M, Mariyanna L, Krepstakies M, Schindler M, Millar DP, Hauber J: Formation of trans-activation competent HIV-1 Rev:RRE complexes requires the recruitment of multiple protein activation domains. PLoS One. 2012, 7: e38305-10.1371/journal.pone.0038305.
Vercruysse T, Daelemans D: HIV-1 Rev multimerization: mechanism and insights. Curr HIV Res. 2013, 11: 623-634. 10.2174/1570162X12666140307094603.
DiMattia MA, Watts NR, Stahl SJ, Rader C, Wingfield PT, Stuart DI, Steven AC, Grimes JM: Implications of the HIV-1 Rev dimer structure at 3.2 A resolution for multimeric binding to the Rev response element. Proc Natl Acad Sci U S A. 2010, 107: 5810-5814. 10.1073/pnas.0914946107.
Hammarskjold MH, Rekosh D: A long-awaited structure is rev-ealed. Viruses. 2011, 3: 484-492. 10.3390/v3050484.
Carpenter S, Dobbs D: Molecular and biological characterization of equine infectious anemia virus Rev. Curr HIV Res. 2010, 8: 87-93. 10.2174/157016210790416424.
Lee JH, Murphy SC, Belshan M, Sparks WO, Wannemuehler Y, Liu S, Hope TJ, Dobbs D, Carpenter S: Characterization of functional domains of equine infectious anemia virus Rev suggests a bipartite RNA-binding domain. J Virol. 2006, 80: 3844-3852. 10.1128/JVI.80.8.3844-3852.2006.
Lee JH, Culver G, Carpenter S, Dobbs D: Analysis of the EIAV Rev-responsive element (RRE) reveals a conserved RNA motif required for high affinity Rev binding in both HIV-1 and EIAV. PLoS One. 2008, 3: e2272-10.1371/journal.pone.0002272.
Gontarek RR, Derse D: Interactions among SR proteins, an exonic splicing enhancer, and a lentivirus Rev protein regulate alternative splicing. Mol Cell Biol. 1996, 16: 2325-2331.
Belshan M, Park GS, Bilodeau P, Stoltzfus CM, Carpenter S: Binding of equine infectious anemia virus rev to an exon splicing enhancer mediates alternative splicing and nuclear export of viral mRNAs. Mol Cell Biol. 2000, 20: 3550-3557. 10.1128/MCB.20.10.3550-3557.2000.
Chung H, Derse D: Binding sites for Rev and ASF/SF2 map to a 55-nucleotide purine-rich exonic element in equine infectious anemia virus RNA. J Biol Chem. 2001, 276: 18960-18967. 10.1074/jbc.M008996200.
Ihm Y, Sparks WO, Lee JH, Cao H, Carpenter S, Wang CZ, Ho KM, Dobbs D: Structural model of the Rev regulatory protein from equine infectious anemia virus. PLoS One. 2009, 4: e4178-10.1371/journal.pone.0004178.
Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008, 9: 40-10.1186/1471-2105-9-40.
Zhang Y: I-TASSER: fully automated protein structure prediction in CASP8. Proteins. 2009, 77 (Suppl 9): 100-113. 10.1002/prot.22588.
Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5: 725-738. 10.1038/nprot.2010.5.
Xu D, Zhang J, Roy A, Zhang Y: Automated protein structure modeling in CASP9 by I-TASSER pipeline combined with QUARK-based ab initio folding and FG-MD-based structure refinement. Proteins. 2011, 79 (Suppl 10): 147-160. 10.1002/prot.23111.
Xu D, Zhang Y: Ab initio protein structure assembly using continuous structure fragments and optimized knowledge-based force field. Proteins. 2012, 80: 1715-1735. 10.1002/prot.24105.
Xu D, Zhang Y: Toward optimal fragment generations for ab initio protein structure assembly. Proteins. 2013, 81: 229-239. 10.1002/prot.24179.
Xu D, Zhang Y: Ab Initio structure prediction for Escherichia coli: towards genome-wide protein structure modeling and fold assignment. Sci Rep. 2013, 3: 1895-
Zhang Y: Interplay of I-TASSER and QUARK for template-based and ab initio protein structure prediction in CASP10. Proteins. 2014, 82 (Suppl 2): 175-187. 10.1002/prot.24341.
Tai CH, Bai H, Taylor TJ, Lee B: Assessment of template-free modeling in CASP10 and ROLL. Proteins. 2014, 82 (Suppl 2): 57-83. 10.1002/prot.24470.
Carpenter S, Chen WC, Dorman KS: Rev variation during persistent lentivirus infection. Viruses. 2011, 3: 1-11. 10.3390/v3010001.
Belshan M, Baccam P, Oaks JL, Sponseller BA, Murphy SC, Cornette J, Carpenter S: Genetic and biological variation in equine infectious anemia virus Rev correlates with variable stages of clinical disease in an experimentally infected pony. Virology. 2001, 279: 185-200. 10.1006/viro.2000.0696.
Sparks WO, Dorman KS, Liu S, Carpenter S: Naturally arising point mutations in non-essential domains of equine infectious anemia virus Rev alter Rev-dependent nuclear-export activity. J Gen Virol. 2008, 89: 1043-1048. 10.1099/vir.0.83195-0.
Wang X, Wang S, Lin Y, Jiang C, Ma J, Zhao L, Lv X, Wang F, Shen R, Kong X, Zhou J: Genomic comparison between attenuated Chinese equine infectious anemia virus vaccine strains and their parental virulent strains. Arch Virol. 2011, 156: 353-357. 10.1007/s00705-010-0877-8.
Wu S, Zhang Y: LOMETS: a local meta-threading-server for protein structure prediction. Nucleic Acids Res. 2007, 35: 3375-3382. 10.1093/nar/gkm251.
Hung LH, Ngan SC, Liu T, Samudrala R: PROTINFO: new algorithms for enhanced protein structure predictions. Nucleic Acids Res. 2005, 33: W77-W80. 10.1093/nar/gki403.
Benkert P, Tosatto SC, Schomburg D: QMEAN: a comprehensive scoring function for model quality assessment. Proteins. 2008, 71: 261-277. 10.1002/prot.21715.
Benkert P, Biasini M, Schwede T: Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011, 27: 343-350. 10.1093/bioinformatics/btq662.
Ray A, Lindahl E, Wallner B: Improved model quality assessment using ProQ2. BMC Bioinformatics. 2012, 13: 224-10.1186/1471-2105-13-224.
Delano WL: The PyMOL Molecular Graphics System. 2002, DeLano Scientific, San Carlos, CA, USA
Lupas A, Van Dyke M, Stock J: Predicting coiled coils from protein sequences. Science. 1991, 252: 1162-1164. 10.1126/science.252.5009.1162.
Lupas A: Prediction and analysis of coiled-coil structures. Methods Enzymol. 1996, 266: 513-525. 10.1016/S0076-6879(96)66032-7.
McDonnell AV, Jiang T, Keating AE, Berger B: Paircoil2: improved prediction of coiled coils from sequence. Bioinformatics. 2006, 22: 356-358. 10.1093/bioinformatics/bti797.
Bartoli L, Fariselli P, Krogh A, Casadio R: CCHMM_PROF: a HMM-based coiled-coil predictor with evolutionary information. Bioinformatics. 2009, 25: 2757-2763. 10.1093/bioinformatics/btp539.
Lupas A: Predicting coiled-coil regions in proteins. Curr Opin Struct Biol. 1997, 7: 388-393. 10.1016/S0959-440X(97)80056-5.
Lupas AN, Gruber M: The structure of alpha-helical coiled coils. Adv Protein Chem. 2005, 70: 37-78. 10.1016/S0065-3233(05)70003-6.
Beck K, Brodsky B: Supercoiled protein motifs: the collagen triple-helix and the alpha-helical coiled coil. J Struct Biol. 1998, 122: 17-29. 10.1006/jsbi.1998.3965.
Grigoryan G, Keating AE: Structural specificity in coiled-coil interactions. Curr Opin Struct Biol. 2008, 18: 477-483. 10.1016/j.sbi.2008.04.008.
Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14: 1188-1190. 10.1101/gr.849004.
Kozakov D, Brenke R, Comeau SR, Vajda S: PIPER: an FFT-based protein docking program with pairwise potentials. Proteins. 2006, 65: 392-406. 10.1002/prot.21117.
Kozakov D, Beglov D, Bohnuud T, Mottarella SE, Xia B, Hall DR, Vajda S: How good is automated protein docking?. Proteins. 2013, 81: 2159-2166. 10.1002/prot.24403.
Comeau SR, Gatchell DW, Vajda S, Camacho CJ: ClusPro: an automated docking and discrimination method for the prediction of protein complexes. Bioinformatics. 2004, 20: 45-50. 10.1093/bioinformatics/btg371.
Comeau SR, Gatchell DW, Vajda S, Camacho CJ: ClusPro: a fully automated algorithm for protein-protein docking. Nucleic Acids Res. 2004, 32: W96-W99. 10.1093/nar/gkh354.
Harris ME, Gontarek RR, Derse D, Hope TJ: Differential requirements for alternative splicing and nuclear export functions of equine infectious anemia virus Rev protein. Mol Cell Biol. 1998, 18: 3889-3899.
Kryshtafovych A, Fidelis K, Tramontano A: Evaluation of model quality predictions in CASP9. Proteins. 2011, 79 (Suppl 10): 91-106. 10.1002/prot.23180.
Kryshtafovych A, Barbato A, Fidelis K, Monastyrskyy B, Schwede T, Tramontano A: Assessment of the assessment: evaluation of the model quality estimates in CASP10. Proteins. 2014, 82 (Suppl 2): 112-126. 10.1002/prot.24347.
Notredame C, Higgins DG, Heringa J: T-Coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-217. 10.1006/jmbi.2000.4042.
Olson SA: MacVector: an integrated sequence analysis program for the Macintosh. Methods Mol Biol. 1994, 25: 195-201.
Rastogi PA: MacVector: integrated sequence analysis for the Macintosh. Methods Mol Biol. 2000, 132: 47-69.
McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics. 2000, 16: 404-405. 10.1093/bioinformatics/16.4.404.
We wish to thank Marit-Nilsen Hamilton for helpful input during the course of the study. This work was supported in part by grants from National Institutes of Health CA128568 (SC), National Science Foundation DBI 0923827 (DD); Iowa State University's Center for Integrated Animal Genomics (DD) and ISU Presidential Initiative for Innovative Research (DD and SC).
The authors declare that they have no competing interests.
SC, CNU, and DD conceived, designed, and implemented the study. HL and KC assisted in protein purification and characterization. JRC contributed towards conception and design of the study. SC supervised the study and provided advice throughout. CNU, SC, and DD contributed equally to writing the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Quality Assessment of Rev Models. Complete list of elongated and globular models for all four sequences modeled (Revs 165, 135, ∆HVR, and FDD). The server from which each model was generated is listed. Each model’s quality scores, calculated using QMEAN and ProQ2, and the average of the two quality scores are also listed. The last column lists the calculated distance between ARM-1 and ARM-2. aE: elongated; G: globular. b: Servers used for protein prediction included QUARK ,; ITASSER ,; LOMETS ; and PROTINFO . cThe QMEAN model quality assessment server is described in ,. dThe ProQ2 model quality assessment server is described in . eDistance (in Ångstrom) was calculated for the pair of closest atoms between ARM-1 and ARM-2. (XLSX 56 KB)
Additional file 2: GenBank Accession Codes and Sequences of the EIAV Rev central region. Complete list of sequences and Genbank accession numbers for the EIAV Rev central region (residues 76–120, based on EIAV R1) used for generating WebLogo of the coiled-coil motif. (PDF 87 KB)
About this article
Cite this article
Umunnakwe, C.N., Loyd, H., Cornick, K. et al. Computational modeling suggests dimerization of equine infectious anemia virus Rev is required for RNA binding. Retrovirology 11, 115 (2014). https://doi.org/10.1186/s12977-014-0115-7
- Bipartite RNA binding domain
- Coiled-coil motif
- Arginine-rich motif