| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Journal of Bacteriology, June 2007, p. 4038-4045, Vol. 189, No. 11
0021-9193/07/$08.00+0 doi:10.1128/JB.01950-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Cheryl J. Whiting,2 and
Lindsay D. Eltis1,2*
Department of Biochemistry, Université Laval, Québec City, Québec, Canada,1 Department of Microbiology and Immunology, Life Sciences Institute, University of British Columbia, Vancouver, BC, Canada,2 Department of Molecular and Cellular Biology, University of Guelph, Guelph, ON N1G 2W1, Canada,3 Department of Microbiology, Immunology and Molecular Genetics, University of Kentucky, Lexington, Kentucky 405364
Received 22 December 2006/ Accepted 6 March 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
|
MCP hydrolases are members of the
/ß-hydrolase superfamily (28). The hydrolyzed CC bond lies between a dienoate (C-1 to C-5) (Fig. 1) and a carbonyl (C-6). Hydrolysis is mediated by conserved Ser-His-Asp residues reminiscent of the catalytic triad of proteases. The proposed catalytic mechanism involves His-mediated ketonization of the dienoate (18, 23), followed by Ser-mediated hydrolysis at C-6. It is unclear whether the catalytic Ser attacks C-6 directly to form an acyl-enzyme intermediate or whether it activates a water molecule which then attacks C-6 (9, 18, 23, 24). Crystal structures have been reported for five MCP hydrolases: BphDLB400 (18), BphDRHA1 from Rhodococcus sp. strain RHA1 (26), CumDIP01 from Pseudomonas fluorescens IP01, involved in cumene degradation (13), CarCJ3 from Janthinobacterium sp. strain J3, involved in carbazole degradation (14), and MhpC (7). All comprise a "core domain," consisting of an eight-stranded ß-sheet flanked by
-helices, and a "lid domain," which occurs as an insertion into the core domain. The active site, located between the core and lid domains, includes the conserved Ser, His, and Asp residues and is divided into two subsites: P (polar) and NP (nonpolar), located on each side of the catalytic serine. The P subsite binds the dienoate moiety, and the NP subsite accommodates the C-6 substituent, a phenyl ring in the case of HOPDA.
Sphingomonas wittichii RW1, one of the best-characterized dibenzofuran-degrading bacteria (37), is able to grow on either dibenzofuran or dibenzo-p-dioxin as a sole organic substrate (36). In addition, the strain partially transforms 2-Cl, 3-Cl, 4-Cl, and 2,3-diCl dibenzofurans to the corresponding chlorinated hydroxysalicylates (36). It is not entirely clear which enzymes are responsible for these transformations, as the bacterium possesses multiple homologues of key catabolic enzymes. For example, Armengaud and coworkers isolated four genes that encode
subunits of ring-hydroxylating dioxygenases and two genes that encode ß subunits (2). The bacterium also contains at least two extradiol dioxygenases and three MCP hydrolases (2, 4). Two of the MCP hydrolases, DxnB and DxnB2 (previously H1 and H2), were purified and shown to hydrolyze HOPDA and 8-hydroxy HOPDA (4). In addition, the gene encoding a third potential MCP hydrolase was cloned (2). DxnB2 has been reported to be more related to MCP hydrolases involved in the catabolism of monocyclic aromatic compounds than to enzymes such as BphD, which are involved in the degradation of bicyclic aromatics (15). This result, derived from amino acid sequence comparisons, is somewhat surprising as the substrates of DxnB and BphD differ by a single hydroxyl group on the phenyl ring. Better characterization of the DxnB enzymes may clarify the phylogenetic relationships of MCP hydrolases and provide insight into why S. wittichii RW1 contains multiple isozymes.
Herein, we report the cloning, sequencing, and heterologous expression of dxnB2. Heterologously produced DxnB2 was purified, and steady-state kinetic studies revealed that the enzyme possesses useful activities towards chlorinated HOPDAs. We generated a structure-based sequence alignment of MCP hydrolases and used this to perform phylogenetic analyses. The results identified a previously unrecognized class of MCP hydrolases and revealed possible structural determinants of substrate specificity in the various classes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Bacterial strains and plasmids. S. wittichii RW1 was purchased from the German Collection of Microorganisms and Cell Cultures (DSMZ, Braunschweig, Germany). This bacterium was grown at 30°C in mineral medium (DSMZ medium 457) with crystalline dibenzofuran or 5 mM sodium salicylate as a carbon source. The vectors and bacterial strains used to test the expression of dxnB2 were the same as those described for bphDP6 and bphDLB400 (30, 32).
Recombinant DNA techniques. DNA was manipulated, and genomic DNA of S. wittichii RW1 was isolated using established protocols (29). DNA amplification reaction mixtures typically contained 20 to 100 ng of DNA template, 0.3 µM of each primer, 0.25 mM deoxynucleoside triphosphates, 1.5 mM MgCl2, and 1.75 U Pfx high-fidelity DNA polymerase (Boehringer Mannheim) in a total volume of 25 µl. PCR was performed using 10 to 30 of the following cycles: denaturation at 94°C for 1 min, annealing at 48°C to 56°C (depending on primers) for 1 min, and extension at 72°C for 2 min. Amplicons were purified using a QIAquick PCR purification kit (QIAGEN).
The dxnB2 gene was amplified from S. wittichii RW1 genomic DNA using a three-step PCR procedure. First, a fragment of dxnB2 was amplified using degenerate primers. The sense primer was GCGGAATTCATATGTTYGARCARTTYGAR, in which the underlined sequence encodes the first six amino acids, MFEQFE, derived from the N-terminal sequence of purified H2 (4). The antisense primer was CCGCCCGGGATCCTGIRCCCARTGICCRCA, in which the underlined sequence corresponds to a highly conserved amino acid sequence, CGHW, found near the C terminus of MCP hydrolases and the bold sequence corresponds to partially conserved codons encoding alanine and glutamine. Restriction sites (EcoRI and NdeI on the sense primer and XmaI and BamHI on the antisense primer) were introduced to facilitate insertion of the resultant amplicon into pT7-7 (33). Amplification yielded a single fragment of the expected 0.8-kbp size. The inserts of two independently obtained clones were sequenced. In the second step, the flanking dxnB2 DNA was amplified using inverse PCR. In this reaction, the template consisted of circularized fragments of genomic DNA that had been generated through PstI digestion and religation using T4 DNA ligase. The primers for inverse PCR were based on the terminal sequences of the 0.8-kbp amplicon: GGGGGAATTCAATGTAATGGGTACG and GCTAGGATCCCGCAGGCGATC. The resultant 1.6-kbp amplicon was inserted into pT7-7 by using EcoRI and BamHI. In the third and final step, dxnB2 was amplified in its entirety together with the upstream Shine-Dalgarno sequence from genomic DNA, using primers based on the nucleotide sequences of two clones containing the amplicon obtained by inverse PCR: GGGGTCTAGAATTGCCGGTCGGTGAATATTTG and GGGGCTGCAGCGATCAATCCAGCTTTCGCGAAG. This third PCR yielded a 0.84-kbp amplicon that was cloned into pT7-7 by using XbaI and PstI, yielding pT7DXN1. The nucleotide sequences of cloned inserts were determined, and the insert was subcloned into pEMBL18 (6) and pVLT31 to yield pEMDXN1 and pVLTDXN1, respectively. To prevent the production of DxnB2 as a fusion protein with the LacZ
-peptide by using pEMDXN1, the pEMBL18 construct was digested with BamHI (upstream of the XbaI insertion site of dxnB2), the 5' overhangs were blunted with Klenow DNA polymerase, and the construct was recircularized by ligation. The production of DxnB2 was tested in three systems: Escherichia coli BL21(DE3) containing pT7DXN1, E. coli DH5
containing pEMDXN1, and Pseudomonas putida KT2442 containing pVLTDXN1.
Protein purification.
DxnB2 was purified from E. coli DH5
containing pEMDXN1. Freshly transformed cells were grown in LB broth at 37°C. When the optical density at 600 nm of the culture reached 0.6, IPTG (isopropyl-ß-D-thiogalactopyranoside) was added to a concentration of 0.5 mM. Incubation was continued overnight, after which the cells were harvested by centrifugation. The cell pellet obtained from 4 liters of culture (16 g [wet weight]) was resuspended in 15 ml of buffer A (20 mM sodium HEPES, pH 7.5) and disrupted by double passage through a French press operated at 20,000 lb/in2. The cell debris was removed by ultracentrifugation (150,000 x g, 60 min), and the supernatant was filtered through a 0.45-µm filter. The filtered supernatant (
20 ml, 0.085 U/mg, 570 mg) was divided into two portions, each of which was loaded onto a Source 15Q (GE Healthcare) anion-exchange column (2 by 9 cm) equilibrated with buffer A. Proteins were eluted using a gradient of 70 to 220 mM NaCl over 10 column volumes. Activity-containing fractions eluting at
140 mM NaCl were pooled, concentrated to about 10 ml by ultrafiltration (Amicon YM 10 membrane), and added to an equal volume of buffer A containing 1 M ammonium sulfate. This mixture (0.29 U/mg, 160 mg) was divided into two portions, each of which was loaded onto a 1- by 9-cm column of phenyl-Sepharose (GE Healthcare) that had been equilibrated with buffer A containing 20% saturation of ammonium sulfate. Proteins were eluted with a gradient of 20% to 0% saturation of ammonium sulfate over 4 column volumes. Activity-containing fractions which eluted at
10% saturation of ammonium sulfate were pooled, concentrated, and exchanged into buffer A by ultrafiltration. The protein preparation (
11 ml, 0.36 U/mg, 61 mg), was then loaded onto a MonoQ HR10/10 (GE Healthcare) column equilibrated with buffer B (20 mM Tris-sulfate, pH 7.5). The proteins were eluted using a gradient of buffer B containing 30 to 150 mM sodium sulfate over 15 column volumes. Activity-containing fractions eluting at
60 mM sodium sulfate were concentrated to about 1 ml (0.44 U/mg, 39 mg), flash frozen as beads in liquid nitrogen, and stored at 80°C.
Recombinant BphDLB400, BphDP6, and 2,3-dihydroxybiphenyl dioxygenase (DHBD) of B. xenovorans LB400 were purified as described previously (30, 32).
Determination of protein concentration, purity, and molecular mass. Protein concentrations were determined using the Bradford method (3a), using bovine serum albumin as the standard. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was performed using a Bio-Rad Mini-Protean II apparatus and stained with Coomassie blue according to established procedures (29). Molecular weight markers of the range 14,000 to 200,000 were used (Bio-Rad).
Steady-state kinetics. Enzymatic activity was measured by following the consumption of the yellow substrate on a Varian Cary 3 spectrophotometer equipped with a thermojacketed cuvette holder. Standard activity assays were performed at 25°C using 0.1 M phosphate, pH 7.5, containing 10 µM HOPDA as previously described (32). One unit of enzymatic activity is defined as the quantity of enzyme required to consume 1 µmol of HOPDA or substituted HOPDA per minute. Extinction coefficients for HOPDA and monochlorinated and dichlorinated HOPDAs were reported previously (30). Extinction coefficients for trichlorinated HOPDAs were determined by measuring the absorbance of solutions containing weighed amounts of trichlorinated DHB (0.1 M phosphate, pH 7.5, 25°C) to which an excess amount of DHBD had been added. Steady-state kinetic parameters were evaluated using 0.1 M phosphate, pH 7.5, at 25°C essentially as described previously (32). Initial velocities were evaluated using substrate concentrations that typically ranged from 0.2 to 10 times the Km value for that substrate. Appropriate equations were fitted to the initial velocities determined at different substrate concentrations by using the least-squares and dynamic weighting options of LEONORA (5). Best-fit parameters were calculated using a minimum of 10 independent data points. The validity of the kinetic models was evaluated using residual plots.
The enzyme-catalyzed hydrolysis of 8-OH HOPDA was monitored at 438 nm using 1 mM potassium phosphate, pH 7.5, at 8°C and an extinction coefficient of 16.2 mM1 cm1, determined as described above. The extinction coefficient for HODPA under these conditions was unchanged from the previously reported value (30). Due to its low stability, 8-OH HOPDA was generated in situ by adding an excess of DHBD to assay mixtures containing 5 or 30 µM of THB. The rates of enzyme-catalyzed reactions were calculated from the first 30 s of the progress curves and were corrected for nonenzymatic transformation of 8-OH HODPA.
Sequence analyses.
All sequences were obtained from the Entrez server at the National Center for Biotechnology Information. Crystal structures were obtained from the Protein Data Bank (http://www.pdb.org/) (accession numbers in parentheses) for MhpCW3110 (1U2E), BphD1RHA1 (1C4X), CumDIP01 (1IUP), CarCJ3 (1J1I), and BphDLB400 (2OG1). A structure-guided alignment of eight MCP hydrolases was obtained in a two-step process. First, the sequences of eight MCP hydrolases of interest were aligned in an automated fashion using the program ClustalX version 1.83 (34) and its default parameters. Second, this alignment was modified using information obtained by superposing the C
backbones of the structures of MhpCW3110, BphD1RHA1, CumDIP01, and CarCJ3 on that of BphDLB400, using the least-squares superposition commands of the program O with default parameters (19). The initial transformation relating each of the four pairs of structures was obtained by superposing the two structures, using the Lsq_exp command and explicitly defined pairs of C
atoms (i.e., those of the conserved catalytic residues). The Lsq_imp command was then used to improve the superposition by locating the longest matching fragments between the two molecules, as determined by the superposition of C
atoms (20). The ClustalX-generated alignment was visually checked against the well-superimposed segments in the structures to eliminate any program-generated errors, such as misaligned residues or segments. Well-superimposed segments included all
-helices and ß-sheets as well as some of the loops. An alignment containing more sequences was obtained by using ClustalX and defining the structure-guided eight-sequence alignment as a profile. The final alignment was inspected to ensure that it was consistent with the structure-guided alignment. Phylogenetic analyses were performed using the maximum-likelihood program of PHYLIP (8). Bootstrap analyses were performed using 100 data sets. Trees were calculated using default parameters except for the jumbling parameter, J, which was set at 10.
| RESULTS |
|---|
|
|
|---|
subunit of salicylate hydroxylase from Pseudomonas aeruginosa JB2 [16]). This result is consistent with DxnB2 being involved in dibenzofuran catabolism, as salicylate would be one of the products of the hydrolase-catalyzed reaction in such a pathway. The nucleotide sequence of the 2.2-kbp fragment was submitted to NCBI (accession no. DQ975235). Sequence and phylogenetic analyses. The amino acid sequences of MCP hydrolases were aligned with the help of information derived from the crystallographic structures of BphDLB400 (18), MhpCW3110 (7), CarCJ3 (14), CumDIP01 (13), and BphDRHA1 (26). The alignment of eight representative sequences, including DxnB2, is presented in Fig. 2. Conserved residues in these enzymes include the catalytic triad (Ser105, His257, and Asp227 of DxnB2) (Fig. 2) and the residues that comprise the P subsite, which binds the dienoate moiety shared by all substrates. These residues include those that interact with the C-1 carboxylate and C-2 hydroxyl of the dienoate moiety: Gly35, Asn43, Asn104, Arg180, His257, and Trp258. The residues comprising the NP subsite, which binds the C-6 substituent of the substrate, occur at positions 106, 131, 132, 135, 146, 149, 150, 196, 200, 203, 206, 229, and 230 in DxnB2 (Fig. 2) and are less well conserved.
|
|
Production and purification of DxnB2.
Among three tested expression systems, DxnB2 was produced to the highest level using E. coli DH5
containing pEMDXN1 based on SDS-PAGE analyses. The heterologously produced DxnB2 was purified using anion-exchange and hydrophobic interaction chromatography. Thirty-nine milligrams of purified protein with a specific activity of 0.44 U/mg was obtained from 4 liters of cell culture. The protein was of >99% purity as judged from SDS-PAGE. The molecular weight of the purified protein (
31,000) corresponds to that of DxnB2 partially purified from its native host (4) as well as to the molecular mass deduced from the nucleotide sequence of the gene (30.4 kDa).
Steady-state kinetic analysis.
The steady-state kinetic parameters of DxnB2 for 8-OH HOPDA, its likely natural substrate, could not be reliably determined under the standard assay conditions due to the compound's rapid nonenzymatic transformation (half-life [t1/2],
18 s), consistent with previous reports (21). Reducing the concentration of phosphate to 1 mM increased the t1/2 of 8-OH HOPDA an order of magnitude (t1/2,
3 min), and lowering the temperature to 8°C resulted in a further fivefold increase (t1/2,
15 min). Under these conditions, DxnB2 catalyzed the hydrolysis of 5 µM 8-OH HOPDA at a rate of 0.30 ± 0.09 U/mg (mean ± standard error). The rate at 30 µM 8-OH HOPDA was similar, indicating that the enzyme was saturated, and the kcat value was 0.15 ± 0.04 s1 (1 mM potassium phosphate, pH 7.5, 8°C). The kcat value for HOPDA under these conditions was approximately sixfold lower (0.026 ± 0.001 s1). By contrast, BphDLB400 catalyzed the hydrolysis of HOPDA (kcat,
1.4 ± 0.1 s1) an order of magnitude more rapidly than for 8-OH HOPDA (kcat,
0.11 ± 0.01 s1) under these conditions.
Under the standard assay conditions, DxnB2 catalyzed the efficient hydrolysis of HOPDA: the specificity constant (kcat/Km) was within a factor of two of those previously determined for BphDLB400 and BphDP6, whose physiological substrate is HOPDA (Table 1). Interestingly, the specificity of DxnB2 for 8-Cl HOPDA, a steric analog of the enzyme's physiological substrate 8-OH HOPDA, was hydrolyzed with a specificity constant that was over an order of magnitude lower than for HOPDA. As observed for the BphDs, other HOPDAs that were chlorinated on the phenyl ring tended to be better substrates of DxnB2 than HOPDA. However, the specificity of DxnB2 appears to be distinct in that it preferentially hydrolyzed 9-Cl HOPDA over 10-Cl HOPDA, in contrast to the BphDs in which this specificity is reversed.
|
Another distinguishing feature of the DxnB2 steady-state hydrolysis of HOPDAs was the occurrence of substrate inhibition. This was observed for all substrates except for 8-Cl HOPDA and 9,11-Cl HOPDA. By contrast, substrate inhibition was not observed for either BphD. In DxnB2, the Ksi values were in the millimolar range while the Km values were in the micromolar range.
| DISCUSSION |
|---|
|
|
|---|
In a previous phylogenetic analysis, DxnB-type MCP hydrolases were classified together with MCP hydrolases involved in the degradation of alkylbenzenes (15). The current analysis is based on a structure-guided sequence alignment. Moreover, the previous classification scheme included
/ß hydrolases that are not MCP hydrolases, such as dehalogenases and epoxide hydrolases (15), which were omitted from our analyses. Nevertheless, it is interesting that the DxnB-type hydrolases appear to be more distinct from the BphD-type hydrolases than are the HsaD-type enzymes, considering that the substrates of the first two enzymes are structurally much more similar. The current analysis reveals several instances of the convergent evolution of substrate specificity. For example, BphD1RHA1 from the PCB-degrading Rhodococcus sp. strain RHA1 is evolutionarily divergent from the class I BphDs, such as BphDLB400, clustering with ThnDTFA from the tetralin catabolic pathway of Sphingomonas macrogoltabida TFA. However, they have functionally converged to transform MCPs of bicyclic compounds. Sequence alignment suggests that the key determinants for recognizing the C-6 phenyl ring are conserved in these two enzymes. Similarly, the three MCP hydrolases that act on 6-propionate HODA, MhpCW3110, HppCPWD1 from Rhodococcus globerulus PWD1, and OhpCRHA1 from Rhodococcus sp. strain RHA1, appear to be evolutionarily divergent despite possessing similar substrate specificities. We expect our understanding of the phylogeny of MCP hydrolases to improve as additional structural and functional data for these enzymes become available, particularly with respect to the enzymes that do not belong to one of the three identified classes. Functional characterization of OrfG7 would be especially interesting in light of the residues that constitute this enzyme's P subsite. In particular, the replacement of Arg180 by Ile in OrfG7 suggests that this enzyme may not be an MCP hydrolase. Arg180 is absolutely conserved among MCP hydrolases, neutralizing the negatively charged C-1 carboxylate of the substrate (7, 18).
The sequence alignment provides useful insight into the structural determinant of substrate specificity in MCP hydrolases. More particularly, the residues comprising the NP subsite to which the C-6 substituent of the substrate binds (positions 106, 131, 132, 135, 146, 149, 150, 196, 200, 203, 206, 229, and 230) appear to correlate with the specificities of the class I, II, and III MCP hydrolases. In the class II enzymes, which transform substrates with smaller alkyl C-6 substituents, the NP subsite residues tend to be bulky hydrophobic residues, such are Phe106, Phe135, Trp150, and Phe197. In contrast, in the class I and class III enzymes, which transform HOPDA, these tend to be smaller aliphatic residues, presumably creating a larger substrate binding pocket to accommodate the C-6 phenyl substituent. Of particular interest are the residues at positions 196 and 203. Residue 196 is conserved as Thr for class III enzymes, except BphD_JF8, and may interact with the ortho-hydroxyl or amine group of the substrate's phenyl ring. In class II enzymes, this residue is Met196 and is part of a conserved motif: MFPXP. Within this motif, Phe197 interacts with two of the key determinants of class II enzymes noted above: Trp150 and Phe106. Residue 203 is notable in that it is lysine or arginine in the divergent hydrolases that utilize substrates with a 6-propionate HODA (MhpC, OhpC, and HppC), suggesting that the positive charge interacts with the negatively charged carboxylate of the substrate. In class III enzymes, residue 203 is asparagine or serine and may form a hydrogen bond with the 8-hydroxyl group of the substrate produced in dibenzofuran degradation. In the other MCP hydrolases, this residue is hydrophobic.
Considering that 8-OH HOPDA is produced during the catabolism of dibenzofuran (Fig. 1) and is the likely physiological substrate of DxnB2, it is somewhat surprising that DxnB2 has a sevenfold-lower specificity for 8-Cl HOPDA than for HOPDA. Structure-based sequence alignment suggests that this may be due to Thr196 and/or Asn203 of DxnB2, which appears to be a significant determinant of substrate specificity in class III MCP hydrolases as discussed above. Thus, the hydrogen bonding that we predict to occur between these residues and the 8-OH group of 8-OH HOPDA could not occur with the chloro substituent of 8-Cl HOPDA. Alternatively, it is possible that the physiological substrate of DxnB2 is substituted 8-OH HODPA. For example, it has been reported that S. wittichii RW1 metabolizes 4-chlorodibenzofuran and 2,7-dichlorodibenzo-p-dioxin to form 9-Cl-, 10-Cl-, or 11-Cl-substituted 8-OH HOPDA (1, 17). Consistent with this possibility, DxnB2 possessed higher specificity for 9-Cl and 10-Cl HOPDAs than for unchlorinated HOPDA (Table 1).
Although DxnB2 transforms HOPDAs substituted on the dienoate moiety better than either BphDLB400 or BphDP6, the effects of such substituents are similar in the three enzymes. For example, the kcat/Km for DxnB2 for 3,9,11-triCl HOPDA was
300-fold lower than for 9,11-diCl HOPDA. This is consistent with the proposed catalytic mechanism of MCP hydrolases, which involves ketonization of the dienoate moiety prior to the hydrolytic step (22). Substituents on the dienoate moiety could interfere with ketonization due to steric and/or electronic effects. Nevertheless, despite the relatively low specificities of the three enzymes for HOPDAs substituted on the dienoate moiety, DxnB2 has a higher specificity for such compounds than either BphDLB400 or BphDP6. For example, DxnB2 has at least a 10-fold-higher specificity towards 3-Cl and 4-OH HOPDAs than either BphDLB400 or BphDP6. Moreover, the Km of DxnB2 for 3-Cl HOPDA is 30-fold higher than for HOPDA, implying that this enzyme is not as strongly competitively inhibited by this compound as BphDLB400. This complementarity in substrate specificity between DxnB2 and BphDLB400 could be exploited to improve the biodegradation of mixtures of PCBs. In light of the current data, it will be interesting to test whether appropriate derivatives of 4,9-DSHA inhibit HsaDH37Rv.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published ahead of print on 6 April 2007. ![]()
Present address: Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA 02115. ![]()
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Appl. Environ. Microbiol. | Infect. Immun. | Eukaryot. Cell |
|---|---|---|
| Mol. Cell. Biol. | J. Virol. | Microbiol. Mol. Biol. Rev. |
| ALL ASM JOURNALS |