JBC Invitrogen Ultrasensitive Cytokine Assays

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Originally published In Press as doi:10.1074/jbc.M102444200 on May 21, 2001

J. Biol. Chem., Vol. 276, Issue 30, 28380-28387, July 27, 2001
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
276/30/28380    most recent
M102444200v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wei, R. R.
Right arrow Articles by Richardson, J. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wei, R. R.
Right arrow Articles by Richardson, J. P.

Identification of an RNA-binding Site in the ATP Binding Domain of Escherichia coli Rho by H2O2/Fe-EDTA Cleavage Protection Studies*

Ronnie R. WeiDagger and John P. Richardson§

From the Departments of Dagger  Biology and § Chemistry Indiana University, Bloomington, Indiana 47405

Received for publication, March 19, 2001, and in revised form, May 17, 2001

    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Transcription factor Rho is a ring-shaped, homohexameric protein that causes transcript termination through actions on nascent RNAs that are coupled to ATP hydrolysis. The Rho polypeptide has a distinct RNA binding domain of known structure as well as an ATP binding domain for which a structure has been proposed based on homology modeling. Treatment of Rho with H2O2 in the presence of Fe-EDTA caused single-cut cleavage at a number of points that coincide with solvent-exposed loops in both the known and predicted structures, thereby providing support for the validity of the tertiary and quaternary structural models of Rho. The binding of ATP caused one distinct change in the cleavage pattern, a strong protection at a cleavage point in the P-loop of the ATP binding domain. Binding of RNA and single-stranded DNA (poly(dC)) caused strong protection at several accessible parts of the oligosaccharide/oligonucleotide binding (OB) fold in the RNA binding domain. RNA molecules but not DNA molecules also caused a strong, ATP-dependent protection at a cleavage site in the predicted Q-loop of the ATP binding domain. These results suggest that Rho has two distinct binding sites for RNA. Besides the site composed of multiples of the RNA binding domain, to which single-stranded DNA as well as RNA can bind, it has a separate, RNA-specific site on the Q-loop in the ATP binding domain. In the proposed quaternary structure of Rho, the Q-loops from the six subunits form the upper entrance to the hole in the ring-shaped hexamer through which the nascent transcript is translocated by actions coupled to ATP hydrolyses.

    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Transcript termination factor Rho is essential for orderly gene expression in many bacteria (1, 2). In Escherichia coli, Rho functions as a homohexamer of a 419-amino acid polypeptide (3-6). Both sequence and functional analyses have shown that Rho polypeptide has two distinct functional domains; one domain is the N-terminal RNA binding domain of residues 1-130 (RNA-BD),1 and the other is C-terminal adenosine triphosphate binding domain of residues 131-419 (ATP-BD). Sequence comparison of Rho with other proteins indicates that it has ribonucleoprotein-like motifs in the RNA-BD and several motifs that are characteristic of ATPases in the ATP-BD (reviewed in Ref. 7). High resolution structures of the RNA-BD by itself (8, 9) and in a complex with oligo(C9) (10) reveal that it contains a classic OB-fold motif that can form a stable complex with RNA (11, 12). A model for the tertiary structure of the Rho ATP-BD has also been proposed based on the close sequence similarity of that part of Rho with the corresponding part of the alpha  and beta  subunits of the F1-ATPase (13).

A model for the quaternary structure of hexameric Rho has also been proposed (7, 13, 14) based on its sequence similarity and morphological resemblance to the F1-ATPase (5, 6, 15). More recently, a low resolution, three-dimensional model of Rho has been constructed from the analysis of electron microscopic images of Rho stained with uranyl acetate (16). This model places the six RNA-BDs at one end of a ring-shaped structure with a C6 symmetry. Although the precise orientation of the RNA binding clefts of the individual domain is uncertain, they are known to come together in the hexamer to form a single continuous cleft that is large enough to protect 60 nucleotides of poly(C) from cleavage by RNase A (4, 17).

In its function as a termination factor, Rho binds initially to an attachment site, called rut (for Rho utilization), on the nascent RNA (1). This interaction presumably involves a direct interaction between the rut sequence on RNA and the RNA-binding site that extends across several RNA-BDs in the hexamer. Termination then occurs as a result of subsequent interactions of Rho with the RNA. These interactions are coupled to ATP hydrolysis and permit Rho to translocate toward the 3' end of the transcript. Several mechanistic models have been proposed to account for the ATP-driven interaction of Rho with the RNA. In two of the models (18, 19), the interactions between Rho and RNA involve only contacts of the RNA with RNA-BDs in various subunits. The third model (7) involves contacts in the ATP-BD as well.

It has been proposed that Rho has two types of RNA-binding sites (20); one type is a primary site that is able to bind single-stranded polynucleotides (RNA or DNA), and the other type is a secondary site that is RNA-specific and that is coupled to interactions with ATP, the hydrolysis substrate. This proposal was based on enzymatic studies in which various polynucleotides were used to activate ATP hydrolysis (20). However, no direct demonstration has been made yet to identify a distinct RNA-specific site on Rho. The binding of poly(dC) by hexameric Rho is competitive with the binding of poly(C) (17), and oligo(dC) binds to the isolated RBA-BD almost as well as does oligo(C) (21). Thus the proposed primary site is very likely composed of the RNA-BD from multiple subunits. The putative secondary site, however, could also be located in RNA-BDs and coincide with conformational changes in some of the subunits, or alternatively, it could be located elsewhere within the Rho hexamer. The results in this paper address this question.

In this study, we used a protein footprinting method similar to the one developed by Heyduk and Heyduk (22) to determine whether the binding of various ligands alters the portions of the Rho polypeptide that are exposed to solvent. In this method, mixtures of H2O2 with Fe-EDTA were used to cause partial chemical cleavage of the protein backbone. Because these reagents cleave proteins preferentially at solvent-accessible sites, this reaction was first used to test for predictions of the structural models. This method was then used to identify possible points of contact with various ligands. Both RNA and DNA caused protection at exposed parts of the OB fold in the RNA-BD, whereas only RNA caused the exposed Q-loop in the ATP-BD to be protected. These observations in conjunction with measurements of ATP hydrolysis suggest that the secondary binding site is at or near the Q-loop.

    EXPERIMENTAL PROCEDURES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Materials-- All nucleotides were purchased from Roche Molecular Biochemicals. Poly(C), poly(U), and poly(A) were purchased from Miles Laboratories, Inc. Poly(dC) was purchased from Amersham Pharmacia Biotech. Eight-residue-long oligocytidine ((Cp)7C) was purchased from Dharmacon Research, Inc. Oligonucleotides (U) and (A) ((Np)7N) were prepared previously (20), and their lengths were confirmed on 20% acrylamide-denaturing gel. Endoproteinase Lys-C and 3,3-diaminobenzidine were products of Sigma. E. coli Rho protein was overexpressed from plasmid pCB111 and purified according to the previously published procedure (23). The peptide corresponding to Rho N-terminal sequence, MNLTELKNTPVSELIT, was synthesized in the laboratory of Dr. Roger Roeske, Department of Biochemistry, Indiana University, Indianapolis, IN, as previously described (24). The peptide corresponding to the C-terminal sequence, KTNDDFFEMMKRS, was synthesized by Research Genetics, Inc. Each peptide had an additional cysteine at the C-terminal end. They were reacted to Imject® maleimide-activated keyhole limpet hemocyanin (Pierce) following the manufacture's instructions. Rabbit antisera were raised against the two conjugated peptides by Cocalico Biologicals, Inc. The collected sera were lyophilized for storage and redesolved in 50% glycerol before use. The secondary antibody, the horseradish peroxidase-conjugated goat anti-rabbit IgG, was purchased from Jackson ImmunoResearch Laboratories, Inc.

The different lambda  cro RNA derivatives were synthesized in vitro by transcription of plasmid pIF2 using T7 RNA polymerase (25). Three DNA templates were prepared by linearizing plasmid pIF2 with either TaqI, BstXI, or BglII restriction enzymes and purified as described in Burgess and Richardson (26). RNA transcripts from these DNA templates were then synthesized and purified as described by Gan and Richardson (27).

H2O2/Fe-EDTA Protein Footprinting-- Each 20-µl reaction mixture contained 0.58 µM Rho hexamer (3.2 µg), 40 mM Hepes-NaOH (pH 7.9), 150 mM potassium acetate. The concentrations of magnesium acetate, adenine nucleotides, and polynucleotides were specified in each experiment. The cleavage was initiated by adding freshly prepared 10× (NH4)2Fe(II)(SO4)2, EDTA (pH 8.0), and H2O2 to reach a final concentration of 1, 2, and 15 mM, respectively. Reactions were carried out at 37 °C for 5 min and terminated by adding 6 µl of 4× SDS loading buffer (200 mM Tris-Cl (pH 8.0), 8% sodium dodecyl sulfate, 0.4% bromphenol blue, 40% glycerol, 0.4 M 2-mercaptoethanol, and 0.2 M thiolurea). The products were then separated by electrophoresis on a polyacrylamide gel composed of a 5% stacking gel, a 10% spacer gel, and a 16.5% separation gel in a Tricine/SDS system (28).

Internal Molecular Weight Marker-- Rho peptide fragments cleaved at Met residues were prepared by digestion of 2 µg of Rho in 20 µl of 10 mg/ml CNBr in 70% formic acid for 40 min at 22 °C. The reaction was terminated by adding 1 ml of deionized distilled water. The products were lyophilized and redesolved in 1× SDS loading buffer before use. Fragments cleaved at Lys residues were prepared by incubation of 2 µg of Rho with 25 ng of Lys-C in 20 µl of 50 mM Tris-Cl (pH 8.0) and M urea at 37 °C for 15 min.

Western Blot-- After the electrophoresis, the separated polypeptides were transferred to nitrocellulose membranes (Schleicher & Schuell, BA 83, 0.2 µm). The electrophoretic transfer was carried out using a Genie electrophoretic blotter (Idea Scientific) at 12 V for 2 h in the electroblotting buffer at 22 °C (25 mM Tris, 192 mM glycine, and 20% methanol, pre-chilled to 4 °C). All the subsequent procedures were done at room temperature. The membranes were air-dried and then treated for 30 min with blocking buffer (5% nonfat dry milk in phosphate-buffered saline (137 mM NaCl, 2.7 mM KCl, 4.3 mM Na2HPO4, and 1.4 mM KH2PO4). The membranes were incubated for 1.5 h with agitation in 30 ml of a solution containing antiserum raised against either the Rho N-terminal or C-terminal peptide conjugate diluted 1:250 in blocking buffer. The blots were then washed three times for 10 min each with phosphate-buffered saline and then incubated for another 1.5 h in a 1:500 dilution of the horseradish peroxidase-conjugated secondary antibody in the blocking buffer. After washing three times as before, the blots were developed in 100 ml of phosphate-buffered saline containing 50 mg 3,3-diaminobenzidine, 0.025% CoCl2, 0.02% NiCl2, and 0.03% H2O2.

Data Analysis-- The blots were scanned and analyzed with IMAGEQUANT (Molecular Dynamics) to obtain the integrated intensity of each cleavage result. These data were then transferred to Microsoft EXCEL and quantified further. The differences in sample loading and cleavage efficiency were corrected based on methods described previously (29, 30). The percentage differences in cleavage intensity (protection) were calculated using the equation (Icomplex - IRho)/Irho × 100, where Icomplex is the corrected intensity for Rho in complex with other ligands, and IRho is the corrected intensity for Rho alone or as indicated in each experiment. The percentage protection values presented were generated by averaging the results from at least three independent experiments.

ATPase Assay-- The ATPase activities of Rho induced by the binding of various nucleic acid ligands were determined as described by Nowatzke and Richardson (31). 80 ng of Rho was incubated with 1 mM ATP at 37 °C for 5 min in 100 µl of the footprinting buffer (40 mM Hepes-NaOH (pH 7.9), 150 mM potassium acetate, 2 mM magnesium acetate) and RNA as indicated for each reaction.

    RESULTS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Cleavage of Rho by Hydrogen Peroxide and Fe-EDTA-- Treatment of proteins with a solution containing H2O2 and Fe(II)-EDTA causes partial cleavage of the polypeptide chain at solvent-exposed segments (22, 29, 30, 32, 33). The procedure is thus useful for probing the structure of proteins and for locating exposed regions that become blocked upon binding of ligands. To apply this approach to the structure and function of Rho factor, we first tested several conditions to find a reaction procedure that allowed partial cleavage of Rho during a relatively short incubation. Under our standard conditions, ~25% of the polypeptides were cleaved after 5 min (Fig. 1). This extent was low enough to avoid more than one cut per molecule (22, 34) but high enough to generate an array of distinct products. We found in the preliminary experiments that ascorbic acid, which was used in previously published procedures (22, 29, 30), was not needed and that reactions depended strongly on the H2O2 concentration. Also, although Fe(II) was used in all the experiments presented here, we obtained identical results with Fe(III). Thus we conclude that the cleavages are caused by actions of the iron-activated H2O2 on the protein and probably not from hydroxyl radicals.


View larger version (9K):
[in this window]
[in a new window]
 
Fig. 1.   Time-dependent cleavage of Rho with H2O2/Fe-EDTA. The cleavage reaction was performed as described under "Experimental Procedures." Samples of the reaction were quenched at the indicated times. The cleavage products were resolved by electrophoresis on a polyacrylamide gel. After staining with 0.1% Coomassie Brilliant Blue, the gel was scanned, and the relative intensity of the stained intact protein was determined using IMAGEQUANT.

To locate the positions of the cleavage sites, we used an indirect immunological detection method (Western blot) to visualize the fragments ending at the N terminus or those ending at the C terminus, depending on the antiserum used (Fig. 2, A and B). One antiserum was raised against a synthetic conjugate containing the N-terminal 16 residues of Rho attached to a carrier protein (keyhole limpet hemocyanin). The other was raised against a similar conjugate with the C-terminal 13 residues of Rho. One drawback of this Western blot approach is that peptides smaller than 150 residues are not retained on the nitrocellulose membranes and are thus not detected. However, we could detect fragments produced by cleavages in the N-terminal two-thirds of Rho by using the C-terminal antiserum (Fig. 2A) and fragments produced by cleavages in the C-terminal two-thirds of Rho by using the N-terminal antiserum (Fig. 2B). Thus, by running the fragments on two gels simultaneously and using both antisera, we were able to locate cleavage sites all across the Rho polypeptide except for the 30 residues closest to the two termini. Cleaved peptides from these two regions could not be resolved from the major uncut protein.


View larger version (56K):
[in this window]
[in a new window]
 
Fig. 2.   Positions of the preferred cleavage points in native Rho protein. Rho was treated with H2O2 and Fe-EDTA, and the cleavage products were separated by gel electrophoresis as described under "Experimental Procedures. A, a blot of separated peptides stained using the anti-C-terminal serum. B, blot stained using the anti-N-terminal serum. Molecular weight markers were generated by residue-specific partial cleavage of Rho by CNBr, which targets methionine residues, or by Lys-C, which cleaves on the carboxyl group side of lysine residues. These residues are shown in the lanes labeled CNBr and LysC, with their approximate positions listed on the left of panels A and B. C, combined Fe-EDTA cleavage profile of Rho polypeptide. Arbitrary intensity was plotted against the residue number. The profile on the left of the dashed line was from panel A and that on the right was from panel B. D, a diagram showing the secondary structural motifs of Rho. The motifs for the RNA-BD were from Ref. 9. The motifs of the model of ATP-BD (13) were determined by Secondary  Render using Protein Data Bank classification in Insight II, BIOSYM/Molecular Simulations. Dotted, filled, and open boxes represent alpha  helices, beta  strands, and loops, respectively. The arrow indicates the trypsin cleavage point that defines the border between the RNA-BD and ATP-BD. Solid bars on top of the diagram mark the positions of the peaks in C.

Molecular weight markers were generated by residue-specific partial digestion (lanes CNBr and LysC in Fig. 2, A and B) and were used to create standard curves that could relate the electrophoretic mobility to the position of cleavage within a few residues by interpolation. A cleavage profile (Fig. 2C) for most of Rho polypeptide was generated by combining the two separate profiles of the Rho-alone lanes in Fig. 2, A and B. This profile revealed that Rho was cleaved with H2O2/Fe-EDTA complex at several distinct points throughout the polypeptide. However, there were also large stretches of residues that were quite resistant to cleavage. Because the H2O2/Fe-EDTA complex cleaves preferentially at solvent-exposed regions of a protein, this procedure allowed us to identify parts of Rho that are solvent-exposed. Fig. 2D is a diagram showing the positions of the known alpha  helices and beta  strands of the RNA-BD as well as the positions of the predicted alpha  helices and beta  strands in the model of the ATP-BD. The positions of cleavage are placed above the diagram as black bars. This alignment shows that nearly all of the major cleavage points are in regions connecting known or predicted secondary structural elements.

In the Rho RNA-BD, the segments of Rho most susceptible to the cleavage were near 61 on loop L12, which connects beta  strands 1 and 2, near residue 70 on loop L23, which connects beta  strands 2 and 3, near residue 85 on the alpha  helix alpha 34 (part of loop L34) that connects beta  strands 3 and 4, and near residue 105 on loop L45, which connects beta  strands 4 and 5 (8, 9). In the ATP-BD, the parts of Rho that were accessible to the chemical cleavage include the P-loop near residue 179, the Q-loop near residue 285, and the R-loop near residue 324.

Fig. 2D also shows that not all known or predicted loops connecting secondary structure elements were cleaved to the same extent. The fact that there was very little cleavage in the middle part of Rho suggests that this region is structurally very compact. Also, many of its loops may be involved in subunit-subunit interaction. With the protein concentration and buffer conditions used in the cleavage reactions, Rho subunits associate predominantly into a hexameric form. This was confirmed by measuring the sedimentation coefficient of Rho at the same concentration and ionic strength as was used in the footprinting assays (data not shown).

ATP Protects the P-loop of Rho-- Rho binds three ATP molecules per hexamer with high affinity (Kd approx  0.2 µM) and perhaps three more ATPs with a lower affinity (35, 36). Stitt and Xu show that the three tightly bound ATPs on each Rho hexamer are rapidly hydrolyzed when RNA is added (37).

To determine whether the interaction of ATP and other adenosine nucleotides has any effect on the cleavage pattern of Rho by the H2O2/Fe-EDTA complex, we conducted separate footprinting reactions in the presence of 1 mM ATP, ATPgamma S, ADP, or AMPPNP. Fig. 3 shows the cleavage products detected with the N-terminal-specific antiserum. The only significant change in cleavage intensity observed was near residue 179, which is on the highly conserved P-loop, a part of the Rho polypeptide that is expected to have a close contact with the beta  and gamma  phosphate of ATP (15, 38, 39). The decrease in the cleavage intensity at this site upon ATP binding was also detected using the C-terminal-specific antiserum (data not shown). No significant change in cleavage intensity in other parts of Rho was evident with either antiserum. ATP and ATPgamma S gave much stronger protection than did ADP and AMPPNP. The differences in the cleavage intensity that resulted from ligand binding compared with that of Rho alone were calculated. The averaged value of the percentage protection ((Icomplex - IRho)/Irho × 100) at residue 179 by each ligand is summarized in Fig. 3. These results provide direct evidence that ATP and ADP interacts with the P-loop of Rho without RNA present and that ATP protects the P-loop from the cleavage by the H2O2/Fe-EDTA complex better than does ADP.


View larger version (72K):
[in this window]
[in a new window]
 
Fig. 3.   Adenine nucleotides protect the cleavage site in the P-loop. A blot of separated cleaved proteins was stained using the anti-N-terminal serum. Each 20-µl reaction mixture contained 3.2 µg of Rho (0.58 µM hexamer), 2 mM magnesium acetate, and 1 mM adenine nucleotide (lane 1, Rho alone; lane 2, ATP; lane 3, ATPgamma S; lane 4, ADP; lane 5, AMPPNP; lane 0, Rho alone without cleavage reagent added). The arrow points to the position of the fragment from cleavage at residue 179 (in the P-loop). The averaged percentage of protection at the P-loop by each ligand (lanes 2-5) was calculated using equation (Icomplex - IRho)/Irho × 100 (Icomplex is the corrected intensity for Rho in complex with ligands, and IRho is the corrected intensity for Rho alone). The values are listed under the image of the blot.

Protection of Rho RNA-BD by Various Polynucleotide Ligands-- We next used partial cleavages with H2O2/Fe-EDTA to probe for changes in Rho upon the addition of various synthetic RNA and DNA molecules. We first used the C-terminal-specific antiserum to detect the fragments arising from cleavage in the RNA-BD. Rho is known to bind to single-stranded C-rich RNA and DNA molecules (17, 21). A major part of this interaction is believed to occur in an extensive cleft formed by an arrangement of the RNA-BDs in hexameric Rho. When Rho and ATPgamma S were mixed with saturating amounts of poly(C) or poly(dC) and the mixture was treated with H2O2 and Fe-EDTA, there was a significant decrease in the cleavage of several of the readily accessible parts of the RNA-BD (Fig. 4A; compare lanes 2 and 3 with lane 1). These cleavage points that became protected were at or near residues 62, 85, and 105 on the L12, alpha 34 (part of L34), and L45, respectively. The percentage of protection at each site is summarized in Fig. 4A. Poly(dC) was just as effective as poly(C) in the extent of protection at the three sites, suggesting that the two polymers bind similarly in this domain. The extent of protection at each of the three positions by poly(C) and poly(dC) was about the same, suggesting a nearly uniform interaction at each of these different points. Poly(U), which binds to Rho much less tightly than Poly(C) (17), gave partial protection at the three points, whereas Poly(A), which has a very weak interaction with Rho (17), did not give significant protection (Fig. 4A, lanes 5 and 6).


View larger version (32K):
[in this window]
[in a new window]
 
Fig. 4.   Both RNA and DNA protect multiple cleavage sites in the RNA-BD. Blots of separated cleaved proteins were stained using the anti-C-terminal serum. A, protection with synthetic homopolynucleotides as ligands. Each reaction contained 3.2 µg of Rho, 3.2 µg of homopolymer (as indicated), 11 mM magnesium acetate, 10 mM ATPgamma S for samples in lanes 1-6, and where indicated, with C8, 40 µM oligo(C8). The arrows point to the positions of the fragments from cleavage at the indicated residue numbers. The percentage of protection at each of the three sites by each ligand is listed under the blot. The protections in lanes 2-6 were calculated with respect to the cleavage in lane 1. The protections in lanes 7 and 8 were calculated with respect to the cleavage of Rho alone without ATPgamma S (not shown). B, diagram showing the segments of the cro RNA derivatives (named for the restriction enzymes used to prepare the DNA templates used for their synthesis) used in C and in Table I. The shadowed boxes indicate the approximate rut site region. C, protections with lambda  cro derivatives. Each reaction contained 3.2 µg of Rho, 1.2 µM indicated RNA, 11 mM magnesium acetate, and 10 mM ATP (lanes 2-4) or no ATP (lanes 5 and 6). The arrows point to the positions of the fragments from cleavage at the indicated residue numbers. The percentage protection at each position is listed under the blot. The protections in lanes 2-4 were calculated with respect to the cleavage in lane 1. The protections in lanes 7 and 8 were calculated with respect to the cleavage of Rho alone without ATP (not shown). nt, nucleotides.

Strong protection at these same elements in the RNA-BD was also observed when Rho was mixed with a full-length lambda  cro gene transcript, a natural RNA that is terminated by the action of Rho (Fig. 4C). Identical protection was achieved with a shorter lambda  cro transcript that extended only just through most of the cro gene rut site (Fig. 4, B and C), but no protection was found with a partial transcript lacking the rut sequence (Fig. 4, B and C). Thus, these results with natural transcripts as well as those with synthetic polymers showed that extensive protection of loop regions in the RNA-BD correlates with high affinity binding of these polymers to Rho (17, 27).

Although poly(dC) binds as tightly as poly(C) to Rho, it does not activate ATPase hydrolysis. However, Rho will hydrolyze ATP when it is mixed with poly(dC) and oligo(C8) (20), (Table I). This result has been interpreted as an indication that in addition to the primary binding site, which is not specific for RNA, there is a secondary binding site that is specific for RNA (20). To determine whether this possible RNA-specific site is within the RNA-BD, the extent of protection was examined for a mixture of poly(dC) with oligo(C8). The result showed that the presence of oligo(C8) did not alter the pattern of protection in the RNA-RBD when compared with the pattern with poly(dC) alone (Fig. 4A, compare lane 4 with lane 3). Thus, this experiment failed to reveal the presence of a possible RNA-specific component in the RNA-BD of Rho. This result essentially confirms the lack of differences in the protection by poly(C) and by poly(dC) in this region of Rho.

                              
View this table:
[in this window]
[in a new window]
 
Table I
ATP hydrolysis of Rho activated by various RNA and DNA polymers
Each reaction mixture contained 80 ng of Rho, 80 ng of RNA or DNA homopolyers, 40 µM oligonucleotide and 0.6 pmol of cro derivatives where indicated.

Another characteristic of the protection seen here in the RNA-BD was the absence of dependence on ATP or ATPgamma S. No significant difference in the protection patterns afforded by all polynucleotides tested in the RNA-BD was observed when ATPgamma S or ATP was omitted from the reaction mixture (compare lanes 7 and 8 with lanes 2 and 3 in Fig. 4A and lanes 5 and 6 with lanes 2 and 3 in Fig. 4C). The only difference seen was the protection at P-loop region afforded by ATP or ATPgamma S as described earlier in Fig. 3.

RNA-specific Protection of Rho at the Q-loop Region in the ATP-BD-- To determine whether RNA-specific protection could be detected in other regions of Rho, the gel-separated, partially cleaved complexes were visualized with the N-terminal-specific antiserum. When Rho was cleaved in the presence of various synthetic polymers, either poly(C) alone or a mixture of poly(dC) with oligo(C8) gave extensive protection at a prominent cleavage point near residue 285 (Fig. 5A, lanes 2 and 6). However, cleavage at that point was not protected by poly(dC) alone (Fig. 5A, lane 5). Thus, unlike the protections seen in the RNA-BD, the protection at residue 285 is specific for RNA. In the tertiary and the quaternary structural models of Rho, residue 285 is predicted to be on a solvent-exposed loop called the Q-loop (14).


View larger version (54K):
[in this window]
[in a new window]
 
Fig. 5.   RNA-specific protection of Rho at the Q-loop region in the ATP-BD. Blots of separated cleaved proteins were stained using the anti-N-terminal serum. A, protection with synthetic homopolymers. Each reaction contained 3.2 µg of Rho, 3.2 µg of homopolymer (as indicated), 11 mM magnesium acetate, 10 mM ATPgamma S for samples in lanes 1-8, and where indicated, with (N)8, 40 µM oligo(N8). The arrow points to the position of the fragment from cleavage at residue 285 (in the Q-loop). The percentage protection by each ligand or combination of ligands is listed under the blot. The protections in lanes 2-8 were calculated with respect to the cleavage in lane 1. The protections in lanes 9 and 10 were calculated with respect to the cleavage of Rho alone without ATPgamma S (as in Fig. 2B). C, protection with cro RNA derivatives. Each reaction contained 3.2 µg of Rho, 1.2 µM RNA, 11 mM magnesium acetate, 10 mM ATP for lane 1-4, and 10 mM ATPgamma S for lane 5-8. The percentage protection at Q-loop by each cro derivative is listed under the blot. The protections in lanes 2-4 were calculated with respect to the cleavage in lane 1. The protections in lanes 6-8 were calculated with respect to the cleavage in lane 5.

Some of the other synthetic polymers also gave significant cleavage protection at residue 285, albeit less extensively than with poly(C) or the poly(dC)-oligo(C8) mixture. Such protection was seen with poly(U) or with a mixture of poly(dC) and oligo(U8) (Fig. 5A, lanes 3 and 7). When the effectiveness of these various polymers in activating ATP hydrolysis was measured (Table I), only those RNAs or combinations of DNA with oligonucleotides that activated ATP hydrolysis gave protection.

When a series of cro RNA derivatives was used as ligands in the footprinting assay under similar conditions, we could not observe significant protection at residue 285 when using ATPgamma S as the cofactor (Fig. 5B, lane 6-8) or when no adenine nucleotide was present (data not shown). However, we were able to detect a significant protection by cro-Taq RNA, the full-length cro transcript, when ATP was present (Fig. 5B, lane 2). Cro-BstX RNA, the cro derivative that ends near the 3' end of the rut site, gave only borderline protection at residue 285, whereas cro-Bgl RNA gave no protection at this region under these conditions (Fig. 5B, lane 3 and 4). Again, the extent of protection with these cro RNA derivatives correlated with their abilities to activate ATP hydrolysis (Table I). Therefore, the segment of transcript that is 3' to the rut site is needed to protect the Q-loop under these conditions, and the presence of ATP is critical to this observed protection pattern.

Although the protection of residue 285 by cro RNA is dependent upon the presence of ATP, this protection can be achieved by poly(C) alone (i.e. with no ATP or ATPgamma S; lane 9, Fig. 5A). The extent of protection with poly(C) was about the same with and without ATPgamma S. However, the protection of residue 285 by the mixture of oligo(C8) with poly(dC) was decreased from 37 to 11% in the absence of ATPgamma S (compare lane 6 to lane 11, Fig. 5A). Thus the protection in this instance correlates with the presence of RNA but not with ongoing ATP hydrolysis. The positions where adenine nucleotides and RNA ligands gave protection are summarized in the proposed Rho quaternary structure model (26) (Fig. 6).


View larger version (49K):
[in this window]
[in a new window]
 
Fig. 6.   Summary of the footprinting results on the model of the quaternary structure of Rho. The two-subunit representation of the proposed quaternary structure is the same as the one shown in Burgess and Richardson (26). The model shows two subunits rotated 180° with respect to each other. RNA-BD is green. ATP-BD is cyan. The P-loop region (residue 175-183) that is protected by adenine nucleotides is yellow. Segments of Rho RNA-BD are protected by RNA and are shown in pink (residue 59-65 on L12), blue (residue 80-90 on alpha 34), and reddish orange (residue 100-110 on L45). The Q-loop region (residue 278-288) in the ATP-BD that is protected by RNA is red.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

In this study we used a protein footprinting technique to locate solvent-exposed regions on termination factor Rho and to determine which of them are the potential binding sites for various ligands. The cleavage profile showed that the positions at which the Rho polypeptide was most susceptible to the H2O2/Fe-EDTA digestion were primarily on parts that are predicted to be solvent-exposed loops in the proposed structural model (Fig. 2D). This evidence has thus provided further support of the validity of that model.

Residues on the P-loop of Rho have been implicated in the binding of ATP through cross-linking studies using ATP analogs (40, 41) and by the functional properties of mutants with changes in conserved P-loop residues (40). In the structure of the complex of ATP with T7 gp4, a DNA helicase that is a structural homolog of Rho, Ser-314, which is its P-loop residue that is equivalent to Lys-181 in Rho, makes a hydrogen bond contact with an oxygen on the gamma -phosphate of ATP (42). The formation of a similar hydrogen bond between Lys-181 in Rho and ATP or ATPgamma S could account for the change in conformation of the part of the P-loop that protects residues near Lys-181 from the cleavage with H2O2/Fe-EDTA. The absence of the gamma -phosphate in ADP is likely to be the reason it gave much less protection than did ATP. Since the concentration of ADP used (1 mM) was well above the Kd for its complex with Rho (35), the lower protection with it cannot be a consequence of incomplete site occupancy. AMPPNP does have a gamma -phosphate group, but it only gave partial protection. Again, because the concentration used was well above the Kd value for the Rho-AMPPNP complex (35), its binding site must have been fully occupied. Thus, the lack of protection must be a result of steric differences, possibly because of small differences in the bond lengths and geometries associated with the phosphoramidate link (43).

We observed that three segments of the Rho RNA-BD were protected by homopolynucleotides and by cro RNAs containing the rut site (Fig. 4). All three segments are located in the OB-fold region. The first protected region is on the loop that connects beta  strands 1 and 2 (residues 58-62). This loop contains several residues that are important for the strength and specificity of RNA binding. Martinez et al. (44) finds that F62S Rho was very defective in binding to RNA, whereas D60G Rho bound RNA more avidly than did wild type Rho but also with less discrimination. In the crystal structure of Rho RNA-BD with oligo(C9), residue Leu-58 and Phe-62 create a hydrophobic pocket for the ribose of the oligonucleotide (10), and their interactions are consistent with our chemical cleavage protection results.

The second protected part is on the loop that connects beta  strands 3 and 4. In some other OB-fold proteins, residues on this loop are also involved in polynucleotide binding (9, 45). However, this region was not found to be involved in the interaction with the RNA ligand in the structure of the crystallized complex of RNA-BD with oligo(C9) (10). This inconsistency may be a consequence of the special characteristics of the complex in the crystal. In that complex, two RNA-BDs share a single oligo(C9) and, thus, has an average site occupancy of 4.5 C residues per RNA-BD. In hexameric Rho, the site occupancy is 10-13 C residues per RNA-BD (17, 46). Thus, the oligo(C9)-RNA-BD complex may have formed without full binding site coverage in each of the subunits. On the other hand, the cleavage protection results are consistent with recent evidence that Arg-88 and Phe-89 display chemical shifts in the NMR spectrum of the isolated RNA-BD upon binding of C-rich oligonucleotides.2

The third protected region is on the loop connecting beta  strands 4 and 5 (residues 103-110). It has been suggested (10) that the side chains of Glu-108, Arg-109, and Tyr-110 on this loop are engaged in Van der Waals interactions with the base of a cytidine, whereas the main chain of this loop forms hydrogen bonds with the bases. The size of the cavity enclosed by L45 and beta  strand 3 is relatively restricted, which can easily shelter a pyrimidine but would be too small for a purine (9, 10). This could explain the preference Rho shows for binding pyrimidine nucleotides and was confirmed by the lack of protection by poly(A) in the footprinting results.

The three sites protected by RNA in the RNA-BD are at disparate locations (Fig. 6), indicating that the binding site extends across all subunits. This result complements the finding that Rho protects a continuous 60-nucleotide segment of poly(C) from digestion with RNase A (4, 17). This RNaseA protection could come about from the RNA being bound in an extensive cleft made from the RNA-BDs of the six subunits arranged in a ring structure. Because the protection of cleavage by H2O2/Fe-EDTA in Rho RNA-BD was essentially the same with poly(dC) as with poly(C) and was the same in the presence of ATP or ATPgamma S as in their absence, this cleft must be the primary polynucleotide binding site.

Contrary to the relatively well defined primary RNA binding site, the location of the proposed secondary binding site has remained elusive. However, the protection by RNA at residue 285 on the Q-loop has several of the characteristics that are expected for the secondary binding site (20). First, this protection was RNA-specific. Both poly(C) and the combination of poly(dC) plus oligo(C8) gave strong protection at the Q-loop, whereas poly(dC) alone did not. Second, this protection showed the following base preference: C U > A. This preference correlates with the ability of oligonucleotides with these residues to activate ATP hydrolysis in conjunction with poly(dC) (Table I). Third, the extent of protection at the Q-loop by some RNAs was affected by the presence of ATP. The observation that the presence of ATPgamma S increases the protection by oligo(C8) in combination with poly(dC) is consistent with the earlier finding that the Km for oligo(C8) in activating ATP hydrolysis decreased when the concentration of ATP was increased (20). The protection of the Q-loop with cro RNAs was even more strictly dependent on ATP (Fig. 5B). On the other hand, poly(C) protected the Q-loop in the absence of ATP (Fig. 5A). This could be a consequence of its extremely high affinity for both the primary and secondary sites of Rho. Possibly, the binding of ATP may lower the stringency for an inherent cytidine requirement at the secondary site. This would be important if the secondary site is the one responsible for translocation, in which case it should be compatible with most RNA sequences.

Protection of the cleavage sites in the Q-loop could be caused by a direct physical contact with RNA or by a conformational change in the loop as a result of the interaction of RNA at another position. We favor the first interpretation for the following reason. In the Rho structural model, the Q-loops of the six subunits form the narrowest part of the central channel in the ring shaped hexamer (Fig. 6), and each Q-loop is right above another loop, called the R-loop, which can be cross-linked to RNA by a 20-Å photo-activable cross-linker (26). This cross-linking result indicates that RNA can pass through the channel. Hence, RNA is likely to make close contact with Q-loop residues because they are in the narrowest part of the channel.

The sequence alignment of Rho polypeptides encoded by the rho genes from several bacteria shows that the Q-loop region has several highly conserved residues (47). None of the known rho mutants changes residues in the Q-loop. We predict that mutational changes of some of these residues should affect the interactions of Rho with RNA that are coupled to ATP hydrolysis.

There have been two hypotheses on the nature of the two types of RNA binding sites on Rho. The first proposes that the two sites are different conformational states assumed by the RNA-BDs on different subunits (18, 19). The other proposes that the primary binding site is composed of multiple RNA-BDs and that there is an additional interaction site for RNA outside of the traditionally defined RNA-BD (7). Our findings in this study support the latter hypothesis and suggest that the secondary binding occurs in the ATP-BD and is closely coupled to ATP hydrolysis.

    ACKNOWLEDGEMENTS

We thank Dr. T. Heyduk and Dr. N. Loizos for gracious help in the footprinting data analysis. We are also grateful for the insight from Dr. B. Stitt concerning the kinetics of ATP hydrolysis.

    FOOTNOTES

* This research was supported by National Institutes of Health Grant GM 56095.The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

To whom correspondence should be addressed: Dept. of Chemistry, Indiana University, 800 E. Kirkwood Ave, Bloomington, IN 47405. Tel.: 812-855-1520; Fax: 812-855-8300; E-mail: richardj@indiana.edu.

Published, JBC Papers in Press, May 21, 2001, DOI 10.1074/jbc.M102444200

2 T. K. Hitchins and G. S. Rule, personal communication.

    ABBREVIATIONS

The abbreviations used are: BD, binding domain; OB, oligosaccharide/oligonucleotide binding; Tricine, N-[2-hydroxy-1,1-bis(hydroxymethyl)ethyl]glycine; ATPgamma S, adenosine 5'-O-(thiotriphosphate); AMPPNP, adenylyl imidodiphosphate.

    REFERENCES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

1. Richardson, J. P., and Greenblatt, J. L. (1996) in Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology (Neidhardt, F. C. , Curtiss, R., III , Ingraham, J. L. , Lin, E. C. C. , Low, K. B. , Magasanik, B. , Reznidoff, W. S. , Riley, M. , Schaechter, M. , and Umbarger, H. E., eds) , pp. 822-848, American Society for Microbiology, Washington, D. C.
2. Richardson, J. P. (1990) Biochim. Biophys. Acta 1048, 127-138
3. Pinkham, J. L., and Platt, T. (1983) Nucleic Acids Res. 11, 3531-3545
4. Bear, D. G., Hicks, P. S., Escudero, K. W., Andrews, C. L., McSwiggen, J. A., and von Hippel, P. H. (1988) J. Mol. Biol. 199, 623-635
5. Oda, T., and Takanami, M. (1972) J. Mol. Biol. 71, 799-802
6. Gogol, E. P., Seifried, S. E., and von Hippel, P. H. (1991) J. Mol. Biol. 221, 1127-1138
7. Richardson, J. P. (1996) J. Biol. Chem. 271, 1251-1254
8. Allison, T. J., Wood, T. C., Briercheck, D. M., Rastinejad, F., Richardson, J. P., and Rule, G. S. (1998) Nat. Struct. Biol. 5, 352-356
9. Briercheck, D. M., Wood, T. C., Allison, T. J., Richardson, J. P., and Rule, G. S. (1998) Nat. Struct. Biol. 5, 393-399
10. Bogden, C. E., Fass, D., Gergman, N., Nichols, M. D., and Berger, J. M. (1999) Mol. Cell 3, 487-493
11. Murzin, A. G. (1993) EMBO J. 12, 861-867
12. Murzin, A. G., Brenner, S. E., Hubbard, T., and Chothia, C. (1995) J. Mol. Biol. 247, 536-540
13. Wood, T. C. (1999) Theory and Application of Protein Homology.Doctoral dissertation , University of Virginia
14. Miwa, Y., Horiguchi, T., and Shigesada, K. (1995) J. Mol. Biol. 254, 815-837
15. Abrahams, J. P., Leslie, A. G. W., Lutter, R., and Walker, J. E. (1994) Nature 370, 621-628
16. Yu, X., Horiguchi, T., Shigesada, K., and Egelman, E. H. (2000) J. Mol. Biol. 299, 1279-1287
17. Galluppi, G. R., and Richardson, J. P. (1980) J. Mol. Biol. 138, 513-539
18. Platt, T. (1994) Mol. Microbiol. 11, 983-990
19. Geiselmann, J., Wang, Y., Seifried, S. E., and von Hippel, P. H. (1993) Proc. Natl. Acad. Sci. U. S. A. 90, 7754-7758
20. Richardson, J. P. (1982) J. Biol. Chem. 257, 5760-5766
21. Briercheck, D. M., Allison, T. J., Richardson, J. P., Ellena, J. F., Wood, T. C., and Rule, G. S. (1996) J. Biomol. NMR 8, 429-444
22. Heyduk, E., and Heyduk, T. (1994) Biochemistry 33, 9643-9650
23. Nowatzke, W. L., Richardson, L. V., and Richardson, J. P. (1996) Methods Enzymol. 274, 353-363
24. Dolan, J. W., Marshall, N. F., and Richardson, J. P. (1990) J. Biol. Chem. 265, 5747-5754
25. Faus, I., and Richardson, J. P. (1989) Biochemistry 28, 3510-3517
26. Burgess, B. R., and Richardson, J. P. (2001) J. Biol. Chem. 276, 4182-4189
27. Gan, E., and Richardson, J. P. (1999) Biochemistry 38, 16882-16888
28. Schägger, H., and von Jagow, G. (1987) Anal. Biochem. 166, 368-379
29. Heyduk, T., Heduk, E., Severinov, K., Tang, H., and Ebright, R. H. (1996) Proc. Natl. Acad. Sci. U. S. A. 93, 10162-10166
30. Loizos, N., and Darst, S. A. (1999) J. Biol. Chem. 274, 23378-23386
31. Nowatzke, W. L., and Richardson, J. P. (1996) J. Biol. Chem. 271, 742-747
32. Rana, T. M., and Meares, C. F. (1991) Proc. Natl. Acad. Sci. U. S. A. 88, 10578-10582
33. Platis, I. E., Ermacora, M. R., and Fox, R. O. (1993) Biochemistry 32, 12761-12767
34. Bernowitz, M., Senear, D., Shea, M. A., and Ackers, G. K. (1986) Methods Enzymol. 30, 133-181
35. Stitt, B. L. (1988) J. Biol. Chem. 263, 11130-11137
36. Geiselmann, J., and von Hippel, P. H. (1992) Protein Sci. 1, 850-860
37. Stitt, B. L., and Xu, Y. (1998) J. Biol. Chem. 273, 26477-26486
38. Story, R. M., and Steitz, T. A. (1992) Nature 355, 374-376
39. Gulick, A. M., Bauer, C. B, Thoden, J. B., and Rayment, I. (1997) Biochemistry 36, 11619-11628
40. Dombroski, A. J., LaDine, J. R., Cross, R. L., and Platt, T. (1988) J. Biol. Chem. 263, 18810-18815
41. O, I., and Stitt, B. L. (1994) J. Biol. Chem. 269, 5009-5015
42. Sawaya, M. R., Guo, S., Tabor, S., Richardson, C. C., and Ellenberger, T. (1999) Cell 99, 167-177
43. Yount, R. G., Babcock, D., Ballantyne, W., and Ojala, D. (1971) Biochemistry 10, 2484-2480
44. Martinez, A. A., Burns, C. M., and Richardson, J. P. (1996) J. Mol. Biol. 257, 909-918
45. Raghunathan, S., Kozlov, A. G., Lohman, T. M., and Waksman, G. (2000) Nat. Struct. Biol. 7, 648-652
46. McSwiggen, J. A., Bear, D. G., and von Hippel, P. H. (1988) J. Mol. Biol. 199, 609-622
47. Opperman, T., and Richardson, J. P. (1994) J. Bacteriol. 176, 5033-5043


Copyright © 2001 by The American Society for Biochemistry and Molecular Biology, Inc.


This article has been cited by other articles:


Home page
J. Biol. Chem.Home page
A. Schwartz, E. Margeat, A. R. Rahmouni, and M. Boudvillain
Transcription Termination Factor Rho Can Displace Streptavidin from Biotinylated RNA
J. Biol. Chem., October 26, 2007; 282(43): 31469 - 31476.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
T. K. Hitchens, Y. Zhan, L. V. Richardson, J. P. Richardson, and G. S. Rule
Sequence-specific Interactions in the RNA-binding Domain of Escherichia coli Transcription Termination Factor Rho
J. Biol. Chem., November 3, 2006; 281(44): 33697 - 33703.
[Abstract] [Full Text] [PDF]


Home page
MicrobiologyHome page
M. S. Ciampi
Rho-dependent terminators and transcription termination.
Microbiology, September 1, 2006; 152(Pt 9): 2515 - 2528.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
Y.-J. Jeong, D.-E. Kim, and S. S. Patel
Nucleotide Binding Induces Conformational Changes in Escherichia coli Transcription Termination Factor Rho
J. Biol. Chem., April 30, 2004; 279(18): 18370 - 18376.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
X. Chen and B. L. Stitt
The Binding of C10 Oligomers to Escherichia coli Transcription Termination Factor Rho
J. Biol. Chem., April 16, 2004; 279(16): 16301 - 16310.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
A. Brevet, J. Chen, S. Commans, C. Lazennec, S. Blanquet, and P. Plateau
Anticodon Recognition in Evolution: SWITCHING tRNA SPECIFICITY OF AN AMINOACYL-tRNA SYNTHETASE BY SITE-DIRECTED PEPTIDE TRANSPLANTATION
J. Biol. Chem., August 15, 2003; 278(33): 30927 - 30935.
[Abstract] [Full Text] [PDF]


Home page
J. Cell Sci.Home page
X. Q. Wang, A. Y. Xiao, C. Sheline, K. Hyrc, A. Yang, M. P. Goldberg, D. W. Choi, and S. Ping Yu
Apoptotic insults impair Na+, K+-ATPase activity as a mechanism of neuronal death mediated by concurrent ATP deficiency and oxidant stress
J. Cell Sci., May 15, 2003; 116(10): 2099 - 2110.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
Y. Xu, H. Kohn, and W. R. Widger
Mutations in the Rho Transcription Termination Factor That Affect RNA Tracking
J. Biol. Chem., August 9, 2002; 277(33): 30023 - 30030.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
276/30/28380    most recent
M102444200v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wei, R. R.
Right arrow Articles by Richardson, J. P.
Right arrow Search for Related Content