![]()
|
|
||||||||
J. Biol. Chem., Vol. 278, Issue 52, 53007-53016, December 26, 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




||
From the
Department of Biochemistry, University of Alberta, Edmonton, Alberta T6G 2H7, Canada, the
Division of Preventative Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02215, and ¶Variagenics/Nuvelo, Inc., Cambridge, Massachusetts 02139
Received for publication, September 12, 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The molecular details of how BRCA1 mutations contribute to the pathogenesis of cancer remain largely unknown. The functional significance of the BRCT region is highlighted by the high degree of sequence conservation within the BRCT regions of among mammalian, Xenopus, and avian BRCA1 homologues (8-10). Several lines of evidence reveal the BRCT is required for tumor suppressor function. A nonsense mutation, which removes 11 C-terminal residues of the second, BRCT (Tyr1853
stop), is associated with early-onset breast cancer (11). Two cancer-linked BRCT missense mutations (12) that destabilize the BRCT fold (13-15), A1708E and M1775R, ablate the double-strand break repair and transcription function of BRCA1 (16) and inhibit BRCT interactions with histone deacetylases (17), BACH1 (18), and the transcriptional co-repressor CtIP (19, 20). Furthermore, mice with homozygous targeted mutations removing the C-terminal half of BRCA1 are viable but develop tumors, suggesting the missing BRCT and/or other domains are expendable for survival, but not for tumor suppression (21).
Although all frameshift or nonsense mutations recorded in the Breast cancer Information Core (BIC) resulting in BRCA1 protein truncation are viewed as functionally deleterious (6, 7), the physiological significance of the majority of missense variants has not been determined due to the absence of a distinctive functional assay for BRCA1. More than 70 missense substitutions have been recorded that alter the primary sequence of the tandem BRCT repeats, but pedigree analysis clarifying the disease linkage of these alleles is available for only eight of these variants (6, 7, 12, 23-27). Many of these amino acid substitutions may be linked with disease but remain as unclassified in the BIC, because the presence of the allele has not been tested in the general population, or the segregation of the allele with disease within a family is unclear (6, 7).
The recent determination of the x-ray crystal structures of the rat and human BRCA1 BRCT repeat domains were important first steps toward understanding tumorigenic BRCT mutations and provide a novel platform for the interpretation of the effects of these alterations in the absence of clinical data (13, 15, 28). In the present study we directly evaluate the consequences missense mutation on the structure of the human BRCA1 BRCT repeats. Using a proteolysis-based assay to probe the BRCT for non-native conformations, we show that the majority of the tested missense and truncations alter the folding state of the BRCT. Cross-validated computational analyses using the BRCT structure and the sequences of proteins homologous to the human BRCT from other organisms further suggest that many of the unclassified BRCT missense mutations are likely to be disease-predisposing and perturb BRCA1 structure/function through protein-destabilizing effects.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
All BRCT single amino acid substitutions were introduced into the BRCT fragment 1646-1859. For missense mutations A1708E, M1775R, and W1837R, mutated BRCA1 coding sequences were used as template for PCR amplification with the FT7 and CD6 primers. All other missense substitutions were engineered using PCR splicing methods (30). Primary PCR mutagenesis reactions used oligonucleotide FT7 with the appropriate reverse (R) mutagenesis oligonucleotide (see below) and CD6 with the appropriate foward (F) mutagenesis oligonucleotide. PCR products from the primary reactions were gel-purified from 1.5% agarose gels using a QIAEX2 kit (Qiagen) and mixed together with oligonucleotides FT7 and CD6 in the secondary PCR splicing reactions to generate mutated PCR products that were subsequently digested with EcoRI/BamHI and ligated to pLM1. The mutagenesis oligonucleotides used were: D1692Y, F-5'-gtt atg aaa aca tat gct gag ttt gtg-3', R-5'-cac aaa ctc agc ata tgt ttt cat aac-3'; F1695L, F-5'-aca gat gct gag ctt gtg tgt gaa cgg-3', R-5'-ccg ttc aca cac aag ctc agc atc tgt-3'; V1696L, F-5'-gat gct gag ttt ttg tgt gaa cgg aca-3', R-5'-tgt ccg ttc aca caa aaa ctc agc atc-3'; C1697R, F-5'-gct gag ttt gtg cgt gaa cgg aca ctg-3', R-5'-cag tgt ccg ttc acg cac aaa ctc agc-3'; R1699W, F-5'-ttt gtg tgt gaa tgg aca ctg aaa tat-3', R-5'-ata ttt cag tgt cca ttc aca cac aaa-3'; R1699Q, F-5'-ttt gtg tgt gaa cag aca ctg aaa tat-3', R-5'-ata ttt cag tgt ctg ttc aca cac aaa-3'; S1715R, F-5'-aaa tgg gta gtt aga tat ttc tgg gtg-3', R-5'-cac cca gaa ata tct aac tac cca ttt-3'; W1718C, F-5'-gtt agc tat ttc tgt gtg acc cag tct-3', R-5'-aga ctg ggt cac aca gaa ata gct aac-3'; T1720A, F-5'-tat ttc tgg gtg gcc cag tct att aaa-3', R-5'-ttt aat aga ctg ggc cac cca gaa ata-3'; G1738E, F-5'-ttt gaa gtc aga gaa gat gtg gtc aat g-3', R-5'-cat tga cca cat ctt ctc tga ctt caa a-3'; G1738R F-5'-ttt gaa gtc aga aga gat gtg gtc aat g-3', R-5'-cat tga cca cat ctc ttc tga ctt caa a-3'; P1749R, F-5'-aac cac caa ggt cgt aag cga gca aga g-3', R-5'-ctc ttg ctc gct tac gac ctt ggt ggt t-3'; R1751Q, F-5'-caa ggt cca aag caa gca aga gaa tcc-3', R-5'-gga ttc tct tgc ttg ctt tgg acc ttg-3'; A1752P, F-5'-ggt cca aag cga cca aga gaa tcc cag-3', R-5'-ctg gga ttc tct tgg tcg ctt tgg acc-3'; I1766S, F-5'-agg ggg cta gaa agc tgt tgc tat ggg-3', R-5'-ccc ata gca aca gct ttc tag ccc cct-3'; M1783T, F-5'-caa ctg gaa tgg acc gta cag ctg tgt g-3', R-5'-cac aca gct gta cgg tcc att cca ggt t; G1788V, F-5'-gta cag ctg tgt gtt gct tct gtg gtg-3', R-5'-cac cac aga agc aac aca cag ctg tac-3'; V1804D, F-5'-ctt ggc aca ggt gac cac cca att gtg-3', R-5'-cac aat tgg gtg gtc acc tgt gcc aag-3'; V1809F, F-5'-cac cca att gtg ttt gtg cag cca gat-3', R-5'-atc tgg ctg cac aaa cac aat tgg gtg-3'; W1837G, F-5'-gtg acc cga gag ggg gtg ttg gac agt g-3', R-5'-cac tgt cca aca ccc cct ctc ggg tca c-3'. All vectors were sequenced to confirm the success of the mutagenesis reactions.
Proteolysis Assays
0.2-0.5 µg of pLM1 plasmid encoding the BRCT variants were used directly as template for protein synthesis reactions with the TNT-Quick in vitro transcription/translation system (Promega). Immediately prior to proteolytic digestion, proteins were translated and labeled with [35S]methionine at 30 °C for 2 h. The reticulocyte lysates were then centrifuged for 2 min at 10,000 x g to remove insoluble material, and 3 µl of the lysate supernatants containing the labeled translation products were added to 12 µl of digestion buffer (150 mM NaCl, 50 mM potassium phosphate, pH 7.5) containing increasing concentrations of trypsin (Sigma) or 1-chloro-3-tosylamido-7-amino-2-heptanone or N
-p-tosyl-L-lysine chloromethyl ketone-treated chymotrypsin (Sigma). After digestion at 20 °C for 12 min, the reactions were stopped with phenylmethylsulfonyl fluoride. Digestion products were electrophoresed on 15% SDS-PAGE gels and visualized with a phosphorimaging plate and a Molecular Dynamics Typhoon scanner. A local average background correction was used during quantification of the reaction products with ImageQuaNT (Amersham Biosciences).
Molecular Graphics
Structural diagrams were created with Bobscript (31, 32) and rendered using Povray (www.povray.org).
Computational Analysis of Risks Associated with Missense Mutations
Method 1: Structure and Sequence-based AnalysisA probability of an effect on function for the missense mutations in BRCT was determined exactly as described using both feature set A and feature set B in Ref. 33. Briefly, the crystal structure (PDB ID: 1JNX
[PDB]
), the multiple sequence alignment for proteins homologous to human BRCT (see Fig. 3), and the chemical nature of the amino acid substitution are used to compute the values of features that are useful for predicting the effects of amino acid substitutions on protein function. For example, the quantitative estimate of solvent accessibility for a residue in a structure or its normalized phylogenetic entropy from a multiple sequence alignment are both features that can be viewed as having a quantitative relationship to the probability of an effect on function for the introduction of a mutant amino acid (33). A probability of an effect on function for a test mutation is estimated by conditional probability as the fraction of training mutations derived from exhaustive mutagenesis of the Lac repressor (34) and T4 lysozyme (35) with an effect on function from among those with feature values are similar to the feature values of the test mutation.
|
| RESULTS |
|---|
|
|
|---|
|
|
-sheet of the C-terminal BRCT in the structures of the human and rat BRCA1-BRCT repeats (Fig. 2D) (13, 28). Hence, the transcriptional activation defects observed for BRCT deletion mutants likely result from destabilization of the protein.
Missense Substitutions Destabilize the BRCT
Similar to the truncation mutants, two cancer predisposing missense mutations, A1708E and M1775R, are destabilizing and exhibit altered BRCT protease susceptibility (13-15). To gain insights into the effects of other patient-derived mutations recorded in the BIC, we generated 23 additional missense variants and tested these proteins for proteolytic sensitivity (see "Experimental Procedures," Fig. 4). 20/25 of the missense mutations tested showed varying degrees of enhanced sensitivity to tryptic digestion at 20 °C (Fig. 4A). Five of six of the mutations that substitute an arginine into the protein (C1697R, S1715R, G1738R, P1749R, and W1837R) also show increased sensitivity to chymotryptic cleavage at 20 °C (Fig. 4B) suggesting that destabilizing effects, rather than the introduction of a new trypsin cleavage site, are responsible for the protease sensitivity. Mutant M1775R is also clearly destabilizing (14) and shows sensitivity to chymotrypsin at elevated temperatures (15).
|
|
|
3-
2 connecting loop (D1692Y and V1696L) confer moderate protease susceptibility to the domain (Fig. 5A). As highlighted by Joo et al. (28), this loop forms a extended
-hairpin structure, participates in the formation of one of the two conserved BRCT surfaces targeted by missense mutations, and may be the primary site of interaction with the BRCA1-associated helicase BACH1. Hence, the BACH1 binding defect reported for mutation F1695L likely results from disruption of a contact site, whereas reduced BACH1 binding for the V1696L mutation may be due to a destabilizing effect, disruption of the contact site, or both. BRCT Interface MutationsThe BRCA1 C-terminal domain consists of two BRCT repeats that pack together head-to-tail via a conserved triple-helical interface, and several of the key residues mediating these BRCT-BRCT contacts are targeted by mutation (Fig. 5B and Table I) (13). Four of five of the tested BRCT interface mutants are destabilizing. Three destabilizing mutants, A1708E, M1775R, and M1783T, likely disrupt the hydrophobic packing between the repeats. The crystal structure of the M1775R variant revealed that mutation-induced structural rearrangements, including flipping of the mutated arginine out of the hydrophobic core of the protein, contribute to fold destabilization (15). The intermediate protease sensitivity of mutation M1783T likely results from combined deleterious effects of protein core cavitation and the burial of a polar hydroxyl group at the BRCT interface. Residue Arg1699 normally participates in a salt bridge between the BRCT repeats. The loss of salt-bridging interactions and steric strain associated with accommodating the tryptophan may contribute to conformational instability of the R1699W mutant. The intermediate stability of R1699W explains the temperature-sensitive transcription phenotype ascribed to this mutation (27). Conversely, R1699Q has little to no effect on BRCT structure and appears to have little effect on transcription activation (25).
BRCA1 Fold MutationsWe have designated a third class of mutations as BRCA1 fold mutants. These substitutions include residues that participate in folding of the BRCT linker region, and residues that do not fall at BRCT fold special positions, but are buried and conserved among BRCA1 homologues (Fig. 5C). All of the tested BRCA1 fold class mutants tested alter the folding of the domain. The majority of these mutants (C1697R, S1715R, G1738R, G1738E, and P1749R) introduce charged residues into the protein core and are highly destabilizing. The A1752P mutant likely disrupts the linker helix and is highly destabilizing. The position of Arg1751 in the crystal structure of human BRCA1 is unclear, but the equivalent residue in the rat structure indicates this residue is involved in salt bridging interactions and the packing of BRCT linker helix (28). This arginine is conserved among all known BRCA1 homologues, and the R1751Q mutation may disrupt similar electrostatic stabilization in the human protein.
BRCT Fold MutationsA conserved hydrophobic clustering signature for the BRCT fold superfamily of proteins was originally identified using sequence-based methods (39, 40). Residues at these positions dictate the fold of an individual BRCT and participate in formation of the BRCT hydrophobic core or are found in turns. As shown in Fig. 5D, 7 of 8 of the mutations tested (W1718C, I1766S, G1788V, V1809F, W1837G, W1837R, and Y1853C) at BRCT fold positions are intermediately or highly destabilizing. The highly destabilizing mutations W1718C, W1837G, and W1837R mutate the invariant BRCT fold tryptophan of both the N- and C-terminal BRCT domains. This critical core residue appears intolerant to both cavitating (W1718C and W1837G) or charge substitution mutations (W1837R) and mediates van der Walls contacts from helix
3 to other secondary structure elements of the BRCT fold, including
-sheet, helix
1, and the 310 helix. Mutation G1788V disrupts the conserved tight turn between
1' and
2' of the C-terminal BRCT. Figs. 2D and 5D highlight the role of Y1853 in positioning the C-terminal BRCT 310 helix that packs against the
-sheet, and substitution of this residue with a cysteine is highly destabilizing. Two of the three tested BRCT fold class mutants that target residues that contribute to intra
-sheet packing, I1766S and V1809F, are destabilizing. The third
-sheet mutant, M1652I, has been classified as a benign polymorphism (41, 42) and does not increase the protease sensitivity of the domain.
Mutations That Destabilize the BRCT Predispose Carriers to Disease
Pedigree analysis clarifying the disease predisposition of BRCA1 alleles is currently available for 8 of the 79 reported BRCA1-BRCT single amino acid substitution variants (Table I). Seven of the mutations, D1692Y, C1697R, R1699W, A1708E, S1715R, P1749R, and M1775R, are destabilizing (Figs. 4 and 5) and are linked to cancer (Table I). In contrast, the frequently recorded BIC polymorphism M1652I exhibits no structural defect, indicating the protease stability assay can successfully discriminate benign mutations from disease-causing variants. The cancer-associated truncation mutants are also protease sensitive (Fig. 2).
Taken together, these results indicate that protease-based detection of altered BRCT stability provides a novel and powerful predictive tool that can be used to assess disease linkage of BRCT mutations in instances where pedigree data is not available. Thus, we suggest that the 20 destabilizing missense mutants and truncations greater than eight amino acids are cancer-predisposing.
Predicting the Structural Consequences of Mutation on the BRCT
The recent development of computational methods that incorporate detailed structural and sequence information to predict the effects of single amino acid substitutions on protein structure/function provides us with alternative tools to study the BRCT mutations (33).2 We have applied two independent methods to predict the potential effects of the 25 missense mutations studied here (Table I) and all known BRCT single amino acid substitutions recorded in the breast cancer information core (www.nhgri.nih.gov/Intramural_research/Lab_transfer/Bic/) (Supplementary Table I).
In the first computational method, a set of quantitative and qualitative features are defined for each amino acid substitution in the BRCT based on its structural disposition, an assessment of the sequence conservation mapped onto the structure, and the potential consequences of introducing the alternative amino acid into the structural environment of the mutated residue. An overall probability of an effect on function is then calculated by comparing the values of the features to the values of features of a large number of mutations with known effects on function (in this case from the Lac repressor and T4 lysozyme (33)). Predictions using two sets of features (A and B in Ref. 33) yielded similar results (see Supplementary Table I). About 51% (35 out of 69) of the mutations are predicted to effect BRCT structure or function with a probability greater than 0.5. Comparison to the experimental analysis of the mutations by proteolysis indicated disagreement for six of 22 of the mutations (0.27), consistent with previous estimates of the cross-validated error described previously (Table I) (33).
The second computational methodology estimates potential effects of mutations on function by extracting as much information as possible from the observed amino acid substitutions in the comparison of the human BRCT sequence to its homologues from other organisms. The method2 relies extensively on Bayesian prior information representing empirically observed amino acid exchangeability in a large number of sequence families (BLOCKS data base (36) and Blocks9 priors (43)) to infer whether an alternative amino acid introduced by mutation is consistent with the biological and chemical character of amino acids found at the corresponding residues of BRCT homologues. This method predicts 14 of 25 of the mutations tested by proteolysis are incompatible with normal BRCT function, representing agreement with the experimental data for 19 of them (0.76). The first method combining structure and sequence is consistent with predictions from method 2 for 18/24 mutations (0.75), and both methods agree with the proteolysis data for 14/24 (0.58) of the mutations.
| DISCUSSION |
|---|
|
|
|---|
Comparison of the transcription and protease-based assay data reveal a striking correlation between destabilizing phenotypes and transcriptional defects (Table I). That is, less stable BRCT variants, including C1697R, R1699W, A1708E, S1715R, G1738E, and M1775R, as well as the truncation mutants, disrupt transactivation function, whereas mutations with no effect on structure (M1652I and R1699Q) are fully active in these assays. However, it has yet to be determined whether BRCT protein misfolding causes BRCA1 tumor suppressor inactivation via BRCA1 transcription function, DNA repair function, or both.
The Protease-based Assay for Ranking BRCT-destabilizing EffectsProteolytic degradation proceeds via an unfolded state for small globular proteins (51, 52) indicating that a correlation between proteolytic resistance and the thermodynamic stability of a protein may exist. This principle forms the basis for phage-based proteolytic selection methods where the evolution of proteins with increased thermodynamic stability closely follows the selection of polypeptides with enhanced resistance to degradation by increasing concentration of protease (53, 54). Thus, the application of a protease-based assay to assess the structural consequences of missense mutations on the BRCT provides a quick, effective, complimentary method to categorize and rank the extent of destabilization of the mutant BRCT proteins.
A recent biophysical assessment of the effects of 8 missense substitutions and the truncation Y1853ter on the thermodynamic stability of the BRCT revealed that four of these missense mutations and the truncation were highly destabilizing and could not be produced as soluble protein in Escherichia coli (14). All four of these missense mutants, A1708E, G1738E, G1788V, and W1837R and the truncation show extreme sensitivity to tryptic digestion (Table I and Figs. 4 and 5). The remaining four (M1775R, M1783T, V1808A, and V1833M) can be produced recombinantly, but destabilize the protein by 3.5-5.5 kcal/mol. Two of these, M1775R and M1783T, show an intermediate sensitivity to proteolysis. Six of the seven other BRCT mutant proteins (V1696L, R1699W, R1751Q, M1783T, V1809F, and Y1853C) with intermediate protease sensitivity are soluble in E. coli when expressed at 20 °C (Table I). Altogether, these data indicate a three-tiered hierarchy of destabilizing effects inferred from the proteolytic data is consistent with results obtained from solubility analysis and direct thermodynamic measurements of BRCT protein stability. Highly destabilizing mutations show sensitivity to low levels of trypsin and tend to be degraded or insoluble when expressed in E. coli. Intermediate thermodynamically destabilizing mutations are sensitive to moderate levels of protease and can be produced in soluble form in E. coli.
The remaining set of mutations may not affect the folding detected by the proteolysis assay and yet still affect the functional properties of human BRCT. The computational methods we have explored represent a first attempt to identify alternative correlates within this class of disease predisposing substitutions. The computational methods were largely consistent with the proteolysis data, whether or not there was an effect on protein stability. The purely sequence-based computational methodology was more consistent with the experimental evidence than the structure- and sequence-based approach. Whether the discrepancy can be interpreted or not remains to be seen through further studies. Methods for predicting the biological consequences of amino acid substitutions is an area of active research, especially because the genome initiatives are discovering too large a number of amino acid altering genetic variants with potential effects on biological function for experimental analysis (for example, see Refs. 55 and 56).
Detection of BRCT MutationsThe observations that the BRCA1 BRCT domains form a proteolytically resistant domain and that cancer-predisposing BRCT variants (but not benign single amino acid substitutions) have compromised stability indicate that a protease-based screen for mutant BRCT conformations could be incorporated into routine BRCA1 screening protocols. The protease-based assay we have used tests the stability of BRCT domains expressed in commercially available reticulocyte lysates. These lysates are employed in protein truncation test genetic screens that have been used for detecting BRCA1 mutations (5, 57, 58). In this method, patient-isolated RNA is reverse-transcribed to generate a cDNA that is then amplified using oligonucleotides that target a coding region of the protein that is to be tested for protein truncating effects. The cDNA message is then transcribed/translated and radiolabeled within the lysate, and the presence of truncated protein product is visualized by SDS-PAGE. Although effective at detecting frameshifts, nonsense and deletion mutations that lead to truncation of the expressed protein message, this technique is incapable of identifying missense substitutions.
An adaptation of the protein truncation test, where a protease digestion step is added could be appropriate for the detection of the large majority of cancer-associated BRCT mutations (Fig. 6A). Here, oligonucleotides would be specifically designed to amplify BRCT coding sequence (amino acids 1646-1863) from patient samples, and the translation step would be followed by a trypsinolysis series. This test would have the distinct advantage of sensing the protein destabilizing effects of both missense and truncation mutations. Conservative estimates indicate it could detect as much as 80% of the cancer-associated mutations that fall within the BRCT coding region. Alternatively, for cases where a BRCT missense mutation has already been detected by sequencing, the mutant BRCT coding sequence could be produced by PCR (Fig. 6B). Direct transcription/translation from the PCR product, followed by protein digestion, would provide a quick, relatively inexpensive test for mutant BRCT conformations.
|
| CONCLUSIONS |
|---|
|
|
|---|
-sandwich and zinc binding mutations that unfold the core DNA binding domain of the p53 tumor suppressor (reviewed in Ref. 22). Conversely, the p53 core domain mutation data base is largely populated by mutations that have little effect on stability, but directly target residues involved in sequence-specific DNA binding. The identification of analogous cancer-associated mutations that are not destabilizing but disrupt specific BRCT protein-protein or protein-DNA binding would provide strong support for the role of these interactions in mediating BRCA1 tumor suppressor function. To this end, two patient-derived mutations (F1695L and V1696L) on the surface of the BRCT that affect BACH1 binding have recently been identified (28). Further biochemical and structural characterization of these interactions will be necessary to confirm the nature of these defects. We have established a set of probability and protease-based criteria on which we can define the structural effects of mutation on the BRCT at the protein level. The early identification of carriers of potentially deleterious BRCA1 alleles is an essential component of breast and ovarian cancer screening programs that facilitates detection, surveillance, and prevention of tumor growth. Further development of complimentary methods that test the destabilizing and biological repercussions of missense variants will provide clinicians and researchers with important tools to unravel BRCA1 function and misfunction.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Table S1. ![]()
|| To whom correspondence should be addressed: Dept. of Biochemistry, University of Alberta, 437 Medical Sciences Bldg., Edmonton, Alberta T6G 2H7, Canada. Tel.: 780-492-2136; Fax: 780-492-0886; E-mail: mark.glover{at}ualberta.ca.
1 The abbreviations used are: BRCT, BRCA1 C-terminal domain; BIC, Breast Cancer Information Core Database. ![]()
2 A. Y. Lau and D. I. Chasman, submitted for publication. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Kumar, W. S. Joo, G. Meinke, S. Moine, E. N. Naumova, and P. A. Bullock Evidence for a Structural Relationship between BRCT Domains and the Helicase Domains of the Replication Initiators Encoded by the Polyomaviridae and Papillomaviridae Families of DNA Tumor Viruses J. Virol., September 1, 2008; 82(17): 8849 - 8862. [Abstract] [Full Text] [PDF] |
||||