![]()
|
|
||||||||
J. Biol. Chem., Vol. 278, Issue 47, 46556-46564, November 21, 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



From the Department of Biochemistry, Kansas State University, Manhattan, Kansas 66506
Received for publication, September 2, 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
-1,3-glucan, initiating activation of a series of proteinases. Phenoloxidase (PO),1 the final enzyme in this cascade, exists in hemolymph as an inactive zymogen, prophenoloxidase (pro-PO), which is activated by a specific proteinase. Active PO catalyzes the hydroxylation of monophenols to o-diphenols and oxidation of o-diphenols to quinones that can polymerize to form melanin. Cytotoxic molecules produced in this process, including quinones and reactive oxygen intermediates, may kill the invading microorganisms that are trapped by melanin (4).
Our understanding of the pro-PO activation cascade is still incomplete, and the number of proteolytic steps in the pathway is unknown. However, serine proteinases that directly activate pro-PO, known as pro-PO activating enzyme (PPAE) or pro-PO-activating proteinase (PAP), have been isolated and cloned from several arthropod species, including the silkworm, Bombyx mori (5, 6), the tobacco hornworm, Manduca sexta (79), a beetle, Holotrichia diomphalia (10, 11), and a crayfish, Pacifastacus leniusculus (12, 13). These enzymes are composed of a carboxyl-terminal serine proteinase domain and at least one clip domain at their amino terminus. A specific proteolytic cleavage between the domains activates the proteinase zymogens, but the two resulting chains remain covalently linked by a disulfide bond (14). The activated proteinase cleaves pro-PO at a conserved Arg-Phe bond at approximately residue 50 and leads to its activation (613).
Three PAPs from M. sexta have been identified and cloned as cDNAs. PAP-1 contains one clip domain (7), whereas PAP-2 and PAP-3 each have two clip domains (8, 9). PAP-1 and PAP-3 are expressed as zymogens at a low constitutive level in naive larvae. Bacterial challenge induces synthesis of all three of the pro-PAPs, whose concentrations increase dramatically in the hemolymph (79). PAP-2 and PAP-3 also accumulate in the hemolymph of prepupae, and PAP-1 is present in the cuticle at the same stage.
The PAPs are synthesized as inactive zymogens and require activation by as yet unidentified proteinases upstream in the pro-PO cascade pathway. Once the PAPs are active, it is likely that they are regulated by specific inhibitors, because the quinones and reactive oxygen species generated by uncontrolled PO activation would be harmful to the host insect (1, 4). Potential candidates for regulators of PAPs are members of the serpin superfamily, which include proteins known to regulate serine proteinase cascades involved in mammalian inflammatory responses, coagulation, complement activation, and fibrinolysis (15, 16). Serpins, present in hemolymph of insects, including M. sexta (17), are proteins of
45 kDa with a solvent-exposed reactive site loop (RSL) near their carboxyl terminus. In an inhibition reaction, hydrolysis of the P1P1' peptide bond is accompanied by formation of a stable acyl intermediate, with the P1 residue of the serpin linked to the catalytic Ser residue of the proteinase. The cleaved serpin undergoes a large conformational rearrangement, in which the RSL inserts into a
-sheet in the serpin (16).
In M. sexta, two serpin genes have been identified (18, 19). The serpin-1 gene contains 10 exons, with 12 alternative forms of exon 9. Alternative mRNA splicing generates 12 serpin variants identical in their amino-terminal 336 residues, encoded by the first eight exons, and a variable carboxyl terminus encoded by one of the exon 9 forms (20). Because exon 9 encodes the reactive center region, serpin-1 variants with different reactive site sequences have different inhibitory selectivity (21). The serpin-1 gene is constitutively expressed in fat body at a high level and, to lesser degree, in the granular hemocytes (22). The serpin-1 proteins are secreted into the plasma. M. sexta serpin-2 is expressed in cytoplasm of granular hemocytes in response to bacterial infection, but its function is not yet understood (19). There is also evidence for the existence of additional serpins in M. sexta, which are possibly encoded by different genes and may regulate the cuticular pro-PO activation (23).
Here we report the molecular cloning and analysis of an immune-responsive serpin (serpin-3) from M. sexta. Serpin-3, whose expression is stimulated by microbial infection, was shown to block the pro-PO activation cascade and to inhibit M. sexta prophenoloxidase-activating proteinases.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Collection of Hemolymph, Hemocytes, and Fat Body from M. sexta LarvaeFor determining the induction pattern of serpin-3 protein in plasma, day 2 fifth instar larvae were injected with 10 µl of filter-sterilized saline (0.85% NaCl), 10 µg/µl Micrococcus luteus (Sigma), 2 x 106 cells/µl Escherichia coli, or 5 x 105 cells/µl Saccharomyces cerevisiae. Hemolymph was collected into microcentrifuge tubes at 0, 1, 2, 6, 12, 24, 32, and 48 h after injection by cutting the dorsal horn. A few crystals of 1-phenyl-2-thiourea were added to each tube to prevent hemolymph melanization. Hemocytes were removed by centrifugation at 12,000 x g for 10 min at 4 °C. For total RNA preparation, hemocytes and fat body were collected from larvae at 24 h after injection as described previously (25).
Molecular Cloning and Sequence AnalysisA serpin-3 cDNA fragment, recovered from clone M5 (GenBankTM accession number BI262710
[GenBank]
) (26), was labeled with [
-32P]dCTP and used as a probe to screen a bacteria-induced M. sexta fat body cDNA library in Uni-Zap XR (Stratagene) (27). The positive plaques were purified and subcloned by in vivo excision. The nucleotide sequence of the longest cDNA insert was determined from both strands using a BigDye Terminator Cycle Sequencing Ready Reaction Kit (PE Applied Biosystems). Sequences were edited with IBI PUSTELL programs. The amino acid sequence deduced from the cDNA was used to search the non-redundant protein data base with the BLAST program (28). The most similar sequences were retrieved and aligned using the PILEUP program from the GCG Sequence Analysis software package (29). The signal peptide was predicted using the program of Nielsen et al. (30).
Detection of an Intron in the serpin-3 GeneM. sexta genomic DNA (0.2 µg) was amplified by PCR in 30 cycles of 95 °C for 30 s, 53 °C for 30 s, and 72 °C for 2 min. The forward primer (5'-TCA AAG GTG TTC CAG AGA AC-3') and reverse primer (5'-ATT GTT TTG CGT GTG GCT TTC-3') were located in the serpin-3 cDNA corresponding to 3' end of exon 8 and the 5' end of exon 10 in the serpin-1 gene, respectively. The PCR product was cloned and sequenced to reveal the sequence of the corresponding region of the serpin-3 gene.
Expression of Recombinant Serpin-3 in E. coliThe serpin-3 cDNA encoding the mature protein from residue 1 to 435 was amplified by PCR using a forward primer (5'-TCT CCC ATG GAC GAC GTC GAC-3') that contains an NcoI site and a reverse primer (5'-TGC GCT GCA GCT CGT TCT AAA T-3') that contains a PstI site. The PCR product digested with NcoI and PstI was purified by agarose gel electrophoresis and cloned into the same sites in H6pQE60 (31). Serpin-3 was expressed in a 1-liter culture of E. coli strain XL1blue carrying the recombinant plasmid. Serpin-3 was expressed in an insoluble form and therefore was purified under denaturing conditions by nickel-affinity chromatography (32). An aliquot of the purified recombinant serpin-3 (1.5 mg) was further resolved by preparative SDS-PAGE and used to produce a polyclonal rabbit antiserum.
Production of Soluble Serpin-3 in a Baculovirus Expression System The BamHI-XhoI cDNA fragment, containing the complete open reading frame for serpin-3, was inserted into pFastBac1 (Invitrogen). A recombinant baculovirus was generated from the resulting plasmid and amplified according to the manufacturer's instructions. For producing soluble serpin-3, Spodoptera frugiperda Sf9 cells (2x 106 cells/ml) in 250 ml of Sf-900 II serum-free medium were infected with the virus at a multiplicity of infection of 4 and incubated at 27 °C with shaking at 140 rpm for 72 h. The conditioned medium was clarified by centrifugation at 500 x g for 10 min, and the supernatant was placed in dialysis tubing and concentrated to 28 ml by covering the tubing with polyethylene glycol 8000 and then dialyzed twice against 1.0 liter of buffer A (20 mM Tris-HCl, 10 mM NaCl, pH 7.5) overnight at 4 °C. The protein solution was applied to a 12-ml Q-Sepharose column and washed with 60 ml of buffer A. The bound proteins were eluted with linear gradients of 10220 mM NaCl (150 ml) and then 220500 mM NaCl (50 ml) in buffer A. Fractions (2 ml each) were analyzed by 10% SDS-PAGE (33) and immunoblotting using serpin-3 antiserum as the first antibody. Concentration of the purified protein was determined by SDS-PAGE analysis along with a series of diluted ovalbumin standards. Intensity of the protein bands was digitized using Kodak Digital Science one-dimensional gel analysis software. The amino-terminal sequence of serpin-3 was determined by automated Edman degradation as described before (7).
Northern Blot Analysis and RT-PCRTotal RNA samples were isolated from hemocytes or fat body using a Rapid Total RNA Isolation kit (Eppendorf-5 Prime, Inc). Northern blot analysis was performed using 20 µg of RNA for each sample (8). The expression profile of serpin-3 was examined by RT-PCR using SuperScript One-step RT-PCR kit and Platinum Taq DNA polymerase (Invitrogen). At 0.5, 1, 2, 4, 6, 8, 12, and 24 h after injection with 100 µg of M. luteus, fat body from four larvae in each group was dissected and pooled for total RNA isolation. One µg of each RNA sample was treated with 2 units of DNase I to remove trace amounts of genomic DNA. Reverse transcription was performed at 45 °C for 30 min, followed by PCR amplification with 20 cycles of 94 °C for 30 s, 50 °C for 30 s, and 72 °C for 90 s using gene-specific primers for serpin-3 (5'-GAT TCC TCG CGA TTC GAT GC-3' and 5'-CAT TTA CGT CAT TAA GTT TCA TG-3'). Primers for ribosomal protein S3 (5'-CGC GAG TTG ACT TCG GT-3' and 5'-GCC GTT CTT GCC CTG TT-3') were used in the RT-PCR as an internal control. PCR products were analyzed by agarose gel electrophoresis with detection by ethidium bromide staining.
Proteinase Inhibition AssaysM. sexta pro-PAP-1 and pro-PAP-3 produced as recombinant proteins by baculovirus expression (34) contained spontaneously activated PAP-1 and PAP-3, which were separated from the zymogen forms by anion exchange chromatography (data not shown). Concentrations of the active PAPs were determined as described above for serpin-3. The purified PAP-1 or PAP-3, each at 50 nM, was preincubated with serpin-3 at different molar ratios (0.4:1 to 4:1) in 10 µl of buffer B (0.1 M Tris-HCl, pH 8.0, 0.03% bovine serum albumin) for 10 min at room temperature. Residual amidase activity was assayed by adding 190 µl of 50 µM N-acetyl-Ile-Glu-Ala-Arg-p-nitroanilide (IEAR-pNA) in buffer C (50 mM Tris-HCl, pH 7.8, 50 mM NaCl) (8).
Determination of the Association Rate Constants (ka)The ka value for PAP-1 and serpin-3 was determined under second order conditions according to Beatty et al. (35). Briefly, recombinant PAP-1 and serpin-3, each at 16 nM, were incubated in 18 µl of buffer B at room temperature. At each designated time (t) (1590 s), 182 µlof50 µM IEARpNa in buffer C was added to determine the residual amidase activity (v). Half-life of the inhibition reaction (t0.5) was determined by plotting 1/v against t, and the ka was calculated according to ka = 1/t0.5x [E0] (35). The ka value for PAP-3 and serpin-3 was also determined under second order conditions with 8.8 nM PAP-3 and 38.6 nM serpin-3. The ka was calculated using the standard Equation 1 (36),
![]() |
Regulation of Pro-PO Activation in Hemolymph by Serpin-3A preparation of hemolymph proteins that contains the components of the pro-PO activation cascade (pro-PO activation fraction) was produced as described previously (21). Ten µg of M. luteus in 1-µl sterile saline was added to 100 µl of pro-PO activating fraction (diluted 1:5 in water). At different incubation times (060 min), 10-µl aliquots of the reaction mixture were added to 700 µl of phenoloxidase substrate solution (2 mM 3-hydroxytyramine in 50 mM sodium phosphate, pH 6.8) to measure PO activity (37). To test the inhibition of pro-PO activation, serpin-3 at final concentrations of 061 µg/ml was added to diluted pro-PO activation fraction (10 µl), along with 1 µg of M. luteus. PO activity was determined after incubation for 10 min at room temperature.
For detecting the formation of serpin-PAP complex, purified serpin-3 was reacted with the purified PAP-1 or PAP-3 in buffer B at a molar ratio of 1:1 or 1:2 (serpin/proteinase). In control samples, diisopropyl fluorophosphate (DFP) at a final concentration of 5 mM was mixed with PAPs prior to the addition of serpin-3. After incubation at room temperature for 10 min, the reaction mixtures were treated with SDS sample buffer at 95 °C for 5 min and resolved by 10% SDS-PAGE. Proteins transferred onto a nitrocellulose membrane were subjected to immunoblot analysis using 1:3000 diluted antisera against PAP-1 (7), PAP-3 (9), or serpin-3 as the first antibody.
Detection of Serpin-PAP Complexes in Hemolymph by Affinity ChromatographyA serpin-3 immunoaffinity column was prepared according to Harlow and Lane (38). Briefly, 2.4 ml of protein A-Sepharose beads (Sigma) were incubated with 4.8 ml of serpin-3 antiserum in 19.2 ml of phosphate-buffered saline (0.8% NaCl, 0.02% KCl, 0.144% NaH2PO4, 0.024% KH2PO4, pH 7.4) for 1 h at room temperature to allow binding of the antibodies to protein A. Following a washing step, covalent coupling of antibodies to protein A was carried out by adding 20 mM dimethylpimelimidate in 0.2 M sodium borate, pH 9.0. After incubation at room temperature for 2 h, the reaction was stopped by incubating the beads with 0.2 M ethanolamine at pH 8.0.
Hemolymph from day 2 fifth instar larvae was collected and centrifuged to remove hemocytes. The plasma was divided into 2 aliquots (2 ml each). DFP was immediately added to 1 aliquot at a final concentration of 5 mM. M. luteus (100 µg) was added to the other aliquot, which was incubated for 5 min before addition of DFP. The two samples were then separately incubated with 0.6 ml of affinity beads in 3.4 ml of phosphate-buffered saline overnight at 4 °C. The resin was packed into two columns and washed with 20 ml of 10 mM sodium phosphate, pH 8.0. The bound proteins were eluted by sequential applications of 0.3 ml of 100 mM glycine-HCl, pH 2.5, with collection of the eluents into 60 µl, 1 M sodium phosphate, pH 8.0. The fractions were subjected to immunoblot analysis as described above.
| RESULTS |
|---|
|
|
|---|
|
Sequence ComparisonsSerpin-3 was most similar in amino acid sequence to a bacteria-induced serpin from another lepidopteran insect, Hyphantria cunea (39), with 60.3% identity. The second most similar sequence was Drosophila serpin 27A (40, 41), which shares 39.6% identity with serpin-3. Serpin-3 is only 29.8% identical to M. sexta serpin-1. A multiple sequence alignment of these insect serpins indicates that M. sexta serpin-3 and H. cunea serpin contain an amino-terminal extension of 4550 amino acid residues that is absent from M. sexta serpin-1 and most other members of the serpin superfamily (Fig. 2). Drosophila serpin 27A has a shorter amino-terminal extension of
30 residues. The putative reactive center loop of serpin-3 is from Gly374 to Tyr397, including Asn-Lys-Phe-Gly at the predicted P2 to P2' positions. The same sequence is present in the H. cunea serpin and Drosophila serpin 27A. With a Lys residue as its putative P1 residue, serpin-3 probably inhibits serine proteinases with trypsin-like specificity.
|
150 mM NaCl at pH 7.5. The purified serpin-3 from the baculovirus expression system (Fig. 3B) was used for characterization and functional analysis in later experiments. It has an apparent molecular mass of 56.8 kDa as estimated by SDS-PAGE, which is larger than the value calculated from the deduced amino acid sequence (48.9 kDa). This difference may be a result of glycosylation. The first five amino acid residues of the purified serpin-3 were determined to be Asp-Asp-Val-Asp-Pro, confirming the prediction that the signal peptide ends at residue 20 (Fig. 1).
|
|
|
|
|
|
90-kDa complex under reducing conditions, which was absent from the controls lanes with only PAP-1 or serpin-3. This complex was recognized by both serpin-3 antibody and PAP-1 antibody (Fig. 8, A and B), indicating that it was composed of these two proteins. The combined apparent molecular mass of serpin-3 (residues 1388) (51.5 kDa) and the PAP-1 catalytic domain (34 kDa) is close to the size of the complex (
90 kDa). When the molar ratio of PAP-1 to serpin-3 was increased from 1:1 to 2:1, more of the serpin-3 was converted to the complex form (Fig. 8B, lanes 2 and 3). A band corresponding to a lower molecular weight protein (51.5 kDa) was also detected by serpin antibody. Because this band was not present in the control of serpin-3 alone, it probably represents the cleaved serpin without the 5.3-kDa carboxyl terminus (residues 389435). This is consistent with the observation that more PAP-1 led to more of the cleaved form of serpin-3. The complex formation was blocked when an irreversible serine proteinase inhibitor, DFP, was preincubated with PAP-1 before addition of serpin-3 to the reaction mixture (Fig. 8B, lane 4). This result indicates that active PAP-1 is required for forming the complex. PAP-1 is composed of an amino-terminal clip domain and a carboxyl-terminal catalytic domain, which are linked by an interchain disulfide bond (7). Under reducing conditions, the clip domain is separated from the catalytic domain. The intensity of the clip domain band recognized by PAP-1 antibodies was similar in each sample (Fig. 9A), whereas the catalytic domain band became much weaker after reaction with serpin-3, further indicating that linkage between serpin-3 and the active site of the catalytic domain resulted in complex formation. Similar results were obtained for the inhibition of PAP-3 by serpin-3. A complex at
90 kDa was formed between serpin-3 (51.5 kDa) and the PAP-3 catalytic domain (39 kDa) (Fig. 8, C and D).
|
|
90-kDa band, which is present at a much lower intensity in the control, untreated sample (Fig. 9A). Similarly, PAP-3 antibodies detected a band at about the same size in M. luteus-treated hemolymph, which was hardly detected in the control (Fig. 9B). This result demonstrated that in plasma, the activated PAP-1 and PAP-3 resulting from the bacteria-triggered activation of the pro-PO cascade can be inhibited by serpin-3, with consequent formation of a stable serpin-3-PAP complex. The complex in naive hemolymph may be because of a trace amount of enzyme that was spontaneously activated during hemolymph collection. The same samples were also examined using serpin-3 antibodies, which revealed two bands at the size range of serpin-3-proteinase complexes (Fig. 9C). We have not yet determined from this experiment which band represents which proteinase complex. | DISCUSSION |
|---|
|
|
|---|
It is not surprising that PAP-1 and PAP-3 share a regulator because the enzymes are similar and they cleave the same substrates. They have an overall amino acid sequence identity of 42%, and many key residues are conserved. PAP-2 and PAP-3 cleave the two subunits of pro-PO at Arg51 within the sequence Asn-Arg-Phe-Gly (7, 8, 37, 48). Serpin-3 contains Asn-Lys-Phe-Gly at its predicted P2-P2' position, which is strikingly similar to the pro-PO cleavage site. This sequence similarity may explain the structural basis for the highly specific and efficient binding between the PAPs and serpin-3.
The cleavage site of pro-POs from different insect species and a crayfish are highly conserved with an occasional substitution of (P2) Asn to Ser or Thr and (P2') Gly to Ser (39). Serpins with a P2-P2' sequence similar to the pro-PO activation site, perhaps homologs of serpin-3, may be present in other insect species to regulate pro-PO activation. Indeed, an immune-responsive H. cunea serpin, which is 60% identical to serpin-3, was shown to inhibit pro-PO activation (39). Its target protease was not identified. Drosophila serpin 27A, which is 40% identical to serpin-3 and which has the same P2-P2' residues, is also immune-inducible and inhibits activation of pro-PO (40, 41). M. sexta serpin-3, H. cunea serpin, and Drosophila serpin 27A all have an extended amino-terminal sequence that is not present in most other serpins. Spn27A mutants have higher PO activity than wild type flies; they develop melanized tissues, and they react to injury with uncontrolled melanization (40, 41). The natural target of Spn27A is unknown, but the serpin inhibits a prophenoloxidase-activating enzyme (PPAE) purified from the beetle H. diomphalia, suggesting that Spn27A may inhibit a Drosophila PPAE (40). M. sexta serpin-1J is another serpin that blocks pro-PO activation (21) and forms complexes with PAP-1, -2, and -3 (9).3 The reactive site of serpin-1J is DRCC, which is less similar to the pro-PO cleavage sequence than the serpin-3 P2-P2' sequence. Once kinetic data for serpin-1J inhibition of the PAPs are available, we can compare the inhibitory efficiency and specificity of serpin-1J and serpin-3.
Serpin-3 is an immune-responsive protein. It is constitutively expressed but rapidly up-regulated at the mRNA level after bacterial injection. PAP-1 and PAP-3 are also constitutively expressed at a low level and are immune-inducible (9, 32). The similar induction patterns among serpin-3 and the PAPs are consistent with regulation of the PAPs by serpin-3. In healthy larvae, serpin-3 may function as a scavenger to inactivate any spontaneously activated PAPs to avoid unnecessary melanization and, thus, protect the host from the cytotoxic molecules produced during this process. The small amount of complex formed between serpin-3 and PAP-1 and PAP-3 in naive larvae and the fast kinetics of the PAP-serpin-3 reactions support this hypothesis. In infected larvae, inhibition of the PAPs by serpin-3 may be necessary to prevent systemic melanization. We expect that serpin-3 acts to limit the duration of the melanization reaction and restrict its location to the infection site.
| FOOTNOTES |
|---|
* This work was supported by National Institutes of Health Grants GM41247 (to M. R. K.) and GM58634 (to H. J.). This is contribution 04-022-J from the Kansas Agricultural Experiment Station. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
Present address: Dept. of Entomology and Plant Pathology, Oklahoma State University, Stillwater, OK 74078. ![]()
To whom correspondence should be addressed: Dept. of Biochemistry, Kansas State University, Manhattan, KS 66506. Tel.: 785-532-6964; Fax: 785-532-7278; E-mail: kanost{at}ksu.edu.
1 The abbreviations used are: PO, phenoloxidase; pro-PO, prophenoloxidase; PAP, prophenoloxidase-activating proteinase; PPAE, prophenoloxidase-activating enzyme; DFP, diisopropyl fluorophosphate; IEARpNa, acetyl-Ile-Glu-Ala-Arg-p-nitroanilide; RT, reverse transcriptase; RSL, reactive site loop. ![]()
2 H. Jiang and Y. Wang, unpublished results. ![]()
3 Y. Wang, H. Jiang, and M. R. Kanost, unpublished results. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Z. Lu, M. H. Beck, Y. Wang, H. Jiang, and M. R. Strand The Viral Protein Egf1.0 Is a Dual Activity Inhibitor of Prophenoloxidase-activating Proteinases 1 and 3 from Manduca sexta J. Biol. Chem., August 1, 2008; 283(31): 21325 - 21333. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Gorman, Y. Wang, H. Jiang, and M. R. Kanost Manduca sexta Hemolymph Proteinase 21 Activates Prophenoloxidase-activating Proteinase 3 in an Insect Innate Immune Response Proteinase Cascade J. Biol. Chem., April 20, 2007; 282(16): 11742 - 11749. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Michel, C. Suwanchaichinda, I. Morlais, L. Lambrechts, A. Cohuet, P. H. Awono-Ambene, F. Simard, D. Fontenille, M. R. Kanost, and F. C. Kafatos Increased melanizing activity in Anopheles gambiae does not affect development of Plasmodium falciparum PNAS, November 7, 2006; 103(45): 16858 - 16863. [Abstract] [Full Text] [PDF] |
||||
![]() |
|