Search

The NKX3.1 transcription factor is an NK family homeodomain protein and a tumor suppressor gene that is haploinsufficient and down-regulated in the early phases of prostate cancer. Like its cardiac homolog, NKX2.5, NKX3.1 acts synergistically with serum response factor (SRF) to activate expression from the smooth muscle gamma-actin (SMGA) gene promoter. Using NMR spectroscopy, three conserved motifs in a construct containing the N-terminal region and homeodomain of NKX3.1 were observed to interact with the MADS box domain of SRF. These motifs interacted both in the absence of DNA and when both proteins were bound to a SMGA promoter DNA sequence. No significant interaction was seen between the homeodomain and SRF MADS box. One of the SRF-interacting regions was the tinman (TN) or engrailed homology-1 motif (EH-1), residues 29-35 (FLIQDIL), which for other NK proteins is the site of interaction with the repressor protein Groucho. A second hydrophobic interacting region was designated the SRF-interacting (SI) motif and included residues 99-105 (LGSYLLD). A third interacting motif was the acidic region adjacent to the SI motif including residues 88-96 (ETLAETEPE). The acidic domain (AD) motif signals also showed strengthening upon the NKX3.1 homeodomain binding to DNA in the absence of SRF, consistent with the acidic region weakly interacting with the homeodomain in the unbound state. The importance of these linear motifs in the transcriptional interaction of NKX3.1 and SRF was demonstrated by targeted mutagenesis of an NKX3.1 expression vector in a SMGA reporter assay. The results implicate the NKX3.1 N-terminal region in regulation of transcriptional activity of this tumor suppressor.

The cardiac-specific Nkx2.5 homeodomain has been expressed as a 79-residue protein with the oxidizable Cys(56) replaced with Ser. The Nkx2.5 or Nkx2.5(C56S) homeodomain is 73% identical in sequence to and has the same NMR structure as the vnd (ventral nervous system defective)/NK-2 homeodomain of Drosophila when bound to the same specific DNA. The thermal unfolding of Nkx2.5(C56S) at pH 6.0 or 7.4 is a reversible, two-state process with unit cooperativity, as measured by differential scanning calorimetry (DSC) and far-UV circular dichroism. Adding 100 mM NaCl to Nkx2.5(C56S) at pH 7.4 increases T(m) from 44 to 54 +/- 0.2 degrees C and DeltaH from 34 to 45 +/- 2 kcal/mol (giving a DeltaC(p) of approximately 1.2 kcal K(-)(1) mol(-)(1) for homeodomain unfolding). DSC profiles of Nkx2.5 indicate fluctuating nativelike structures at <37 degrees C. Titrations of specific 18 bp DNA with Nkx2.5(C56S) in buffer at pH 7.4 with 100 mM NaCl yield binding constants of 2-6 x 10(8) M(-)(1) from 10 to 37 degrees C and a stoichiometry of 1:1 for homeodomain binding DNA, using isothermal titration calorimetry. The DNA binding reaction of Nkx2.5 is enthalpically controlled, and the temperature dependence of DeltaH gives a DeltaC(p) of -0.18 +/- 0.01 kcal K(-)(1) mol(-)(1). This corresponds to 648 +/- 36 A(2) of buried apolar surface upon Nkx2.5(C56S) binding duplex B-DNA. Thermodynamic parameters differ for Nkx2.5 and vnd/NK-2 homeodomains binding specific DNA. Unbound NK-2 is more flexible than Nkx2.5.

The three-dimensional solution structure obtained by NMR of the A35T mutant vnd/NK-2 homeodomain bound to the vnd/NK-2 consensus 16 bp DNA sequence was determined. This mutation to threonine from alanine in position 35 in helix II of the vnd/NK-2 homeodomain is associated with early embryonic lethality in Drosophila melanogaster. Although the unbound mutant protein is not structured, in the DNA-bound state it adopts the three-helix fold characteristic of all known homeodomains, but with alterations relative to the structure of the wild-type analogue. These structural modifications occur, and are accompanied by a 50-fold reduction in the DNA binding affinity, even though most of the protein-DNA interactions originally seen for the wild-type homeodomain are found likewise in the threonine analogue. Alterations include torsional angle changes in the loop between helix I and helix II, and in the turn between helix II and helix III, as well as in a distortion of the usual antiparallel orientation of helix I with respect to helix II. The alteration of the position of leucine 40 in the A35T mutant is proposed to explain the observed 1.27 ppm upfield shift of the corresponding amide proton resonance relative to the value observed for the wild-type analogue. A detailed comparison of the structures of the mutant A35T and wild-type vnd/NK-2 homeodomains bound to the cognate DNA is presented. The consequences of the structural alteration of the DNA-bound A35T mutant vnd/NK-2 protein may constitute the basis of the observed early embryonic lethality.

The importance in downstream target regulation of tertiary structure and DNA binding specificity of the protein encoded by the vndNK-2 homeobox gene is analyzed. The ectopic expression patterns of WT and four mutant vndNK-2 genes are analyzed together with expression of two downstream target genes, ind and msh, which are down-regulated by vndNK-2. Three mutants are deletions of conserved regions (i.e., tinman motif, acidic motif, and NK-2 box), and the fourth, Y54M vndNK-2, corresponds to a single amino acid residue replacement in the homeodomain. Of the four ectopically expressed mutant genes examined, only the Y54M mutation inactivates the ability of the vndNK-2 homeodomain protein to repress ind and msh. The acidic motif deletion mutant slightly reduced the ability of the protein to repress ind and msh. By contrast, both tinman and NK-2 box deletion mutants behaved as functional vndNK-2 genes in their ability to repress ind and msh. The NMR-determined tertiary structures of the Y54M vndNK-2 homeodomain, both free and bound to DNA, are compared with the WT analog. The only structural difference observed for the mutant homeodomain is in the complex with DNA and involved closer interaction of the methionine-54 with A2, rather than with C3 of the (-) strand of the DNA. This subtle change in the homeodomain-DNA complex resulted in modifications of binding affinities to DNA. These changes resulting from a single amino acid residue replacement constitute the molecular basis for the phenotypic alterations observed on ectopic expression of the Y54M vndNK-2 gene during embryogenesis.

The (15)N relaxation behavior and heteronuclear Overhauser effect data for the wild type and an H52R/T56W double mutant protein that encompasses the vnd/NK-2 homeodomain from Drosophila melanogaster were used to characterize and describe the protein backbone dynamics. This investigation, which includes a description of a model structure for the H52R/T56W double mutant vnd/NK-2 homeodomain, was carried out for the two proteins in both the free and DNA-bound states. The double residue replacement at positions 52 and 56 within the DNA recognition helix of vnd/NK-2 has been shown to lead to a significant secondary structural modification resulting in an increase in the length of the recognition helix for the unbound protein. These structural changes are accompanied by corresponding changes in the T(1) and T(1)(rho) relaxation times as well as in the heteronuclear Overhauser effect (XNOE) values that show that the structural stability of the protein is enhanced by the two residue replacements. The values of the rotational anisotropy, D(parallel)/D(perpendicular), derived from analysis of the (15)N T(1) and T(1)(rho) relaxation values are small (1.189 for the unbound homeodomain and 1.110 for the bound homeodomain; both analyzed as prolate ellipsoids of revolution). A comparison of the T(2) values of the wild type and double mutant homeodomain reveals the presence of a low-frequency exchange contribution for the wild type analogue. These relaxation studies show that the motional behavior of the protein primarily reflects the tertiary structure and stability of the homeodomain backbone as well as the respective changes induced upon site-directed residue replacement or DNA binding.

The conformational stabilities of the vnd (ventral nervous system defective)/NK-2 homeodomain [HD(wt); residues 1-80 that encompass the 60-residue homeodomain] and those harboring mutations in helix III of the DNA recognition site [HD(H52R) and HD(H52R/T56W)] have been investigated by differential scanning calorimetry (DSC) and ellipticity changes at 222 nm. Thermal unfolding reactions at pH 7.4 are reversible and repeatable in the presence of 50-500 mM NaCl with DeltaC(p) = 0.52 +/- 0.04 kcal K(-1) mol(-1). A substantial stabilization of HD(wt) is produced by 50 mM phosphate or by the addition of 100-500 mM NaCl to 50 mM Hepes, pH 7.4, buffer (from T(m) = 35.5 degrees C to T(m) 43-51 degrees C; DeltaH(vH) congruent with 47 +/- 5 kcal mol(-1)). The order of stability is HD(H52R/T56W) > HD(H52R) > HD(wt), irrespective of the anions present. Progress curves for ellipticity changes at 222 nm as a function of increasing temperature are fitted well by a two-state unfolding model, and the cooperativity of secondary structure changes is greater for mutant homeodomains than for HD(wt) and also is increased by adding 100 mM NaCl to Hepes buffer. A 33% quench of the intrinsic tryptophanyl residue fluorescence of HD(wt) by phosphate binding (K(D)' = 2.6 +/- 0.3 mM phosphate) is reversed approximately 60% by DNA binding. Thermodynamic parameters for vnd/NK-2 homeodomain proteins binding sequence-specific 18 bp DNA have been determined by isothermal titration calorimetry (10-30 degrees C). Values of DeltaC(p) are +0.25, -0.17, and -0.10 +/- 0.04 kcal K(-1) mol(-1) for HD(wt), HD(H52R), and HD(H52R/T56W) binding duplex DNA, respectively. Interactions of homeodomains with DNA are enthalpically controlled at 298 K and pH 7.4 with corresponding DeltaH values of -6.6 +/- 0.5, -10.8 +/- 0.1, and -9.0 +/- 0.6 kcal mol(-1) and DeltaG' values of -11.0 +/- 0.1, -11.0 +/- 0.1, and -11.3 +/- 0.3 kcal mol(-1) with a binding stoichiometry of 1.0 +/- 0.1. Thermodynamic parameters for DNA binding are not predicted from homeodomain structural changes that occur upon complexing to DNA and must reflect also solvent and possibly DNA rearrangements.

Purpose: The NKX-3.1gene is an androgen regulated prostate specific homeobox gene that is believed to have a vital role in normal prostate development. In mice the homologue NKx3.1 is exclusively expressed in prostate epithelium. In humans NKX3.1 expression is also restricted to the prostate but to our knowledge the cellular location has not been described. Furthermore, since NKX3.1 maps to chromosomal band 8p21, a region with high loss of heterozygosity in prostate cancer, the gene has been proposed to have tumor suppressor function. In this study we demonstrate that in human prostates NKX3.1 is expressed exclusively in secretory epithelial cells and the level of NKX3.1 expression remains invariant in normal tissue and in tissue showing various grades of prostate cancer. In the 19 cases examined the DNA sequences of the NKX3.1 gene were identical and no mutation was detected., Materials and Methods: Frozen tissue from patients who underwent radical prostatectomy was used for this study. For in situ hybridization experiments a 377 bp fragment corresponding to a portion of the 3' untranslated region of the NKX3.1 gene was amplified by polymerase chain reaction and cloned into the pCRII plasmid vector Invitrogen. Antisense or sense [33P] uridine triphosphate labeled RNA probes were generated with SP6 or T7 RNA polymerase and hybridized to the tissue sections. Slides were exposed to photographic emulsion and visualized on autoradiography. Laser capture microdissection was performed to procure pure populations of malignant epithelium. DNA was isolated by digesting samples in proteinase K buffer. Polymerase chain reaction and direct sequencing was performed using standard protocols., Results: In vitro hybridization showed that NKX3.1 expression was restricted to secretory epithelial cells within benign prostate glands. No expression was detected in stroma or infiltrating lymphocytes. NKX3.1 was expressed in all grades of malignant epithelium in all 25 cases examined. Direct sequencing of the coding region of NKX3.1 revealed the wild-type sequence in all 18 microdissected cancers analyzed., Conclusions: Based on our studies we propose that NKX3.1 gene expression is restricted to benign and malignant secretory epithelium within the prostate but NKX3.1 does not appear to be a classic tumor suppressor gene responsible for prostate cancer initiation. These findings are consistent with the role of NKX3.1 in the development of normal prostate epithelium and maintenance of normal secretory function. Thus, NKX3.1 may represent a useful molecular marker for benign and malignant prostate epithelium.

The mode coupling diffusion theory is applied to the derivation of local dynamics in proteins in solution. The rotational dynamics of the bonds along the protein sequence are calculated and compared to the experimentally measured nmr (15)N spin-lattice relaxation time T(1), at 36.5, 60.8, and 81.1 MHz of the vnd/NK-2 homeodomain from Drosophila melanogaster. The starting point for the calculations is the experimental three-dimensional solution structure of the homeodomain determined by multidimensional nmr spectroscopy. The higher order mode-coupling computations are compared also with the recently published first-order approximation calculations. The more accurate calculations improve substantially the first-order ORZLD calculations and show that the role of the strength of the hydrodynamic interactions becomes crucial to fix the order of magnitude of the rotational dynanics for these very compact molecules characterized by partial screening of the internal atoms to water. However, the relative mobility of the bonds along the sequence and the differential fluctuations depend only weakly on the hydrodynamic strength but strongly on the geometry of the three-dimensional structure and on the statistics incorporated into the theory. Both rigid and fluctuating dynamic models are examined, with fluctuations evaluated using molecular dynamics simulations. The comparison with nmr data shows that mode coupling diffusion accounts for the T(1) relaxation pattern at low frequency where the rotational tumbling dominates. An important contribution of internal motions in the nanosecond time scale is seen at high frequencies and is discussed in terms of diffusive concepts., (Copyright 2000 John Wiley & Sons, Inc.)

A novel implementation of the water flipback technique employing a 45 degrees flip-angle water-selective pulse is presented. The use of this water flipback technique is shown to significantly enhance signal in 3D 15N-edited ROESY in a 20 kDa complex of the vnd/NK-2 homeodomain bound to DNA. The enhancement is seen relative to the same experiment using weak water presaturation during the recovery delay. This enhancement is observed for the signals from both labile and nonlabile protons. ROESY and NOESY pulse sequences with 45 degrees water flipback are presented using both HMQC and HSQC for the 15N dimension. The 45 degrees flipback pulse is followed by a gradient, a water selective 180 degrees pulse, and another gradient to remove quadrature images and crosspeak phase distortion near the water frequency. Radiation damping of the water magnetization during the t1 and t2 evolution periods is suppressed using gradients. Water resonance planes from NOESY-HMQC and NOESY-HSQC spectra show that the HMQC version of the pulse sequences can provide stronger signal for very fast exchanging protons. The HSQC versions of the ROESY and NOESY pulse sequences are designed for the quantitative determination of protein-water crossrelaxation rates, with no water-selective pulses during the mixing time and with phase cycling and other measures for reducing axial artifacts in the water signal.

The three-dimensional solution structure obtained by NMR of the complex formed between the uniformly singly15N and doubly13C/15N-labeled vnd/NK-2 homeodomain and its consensus 16 base-pair DNA binding sequence was determined. This work was carried out using the accepted repertoire of experiments augmented with a novel implementation of the water flipback technique to enhance signals from exchangeable amide protons. The results using this new technique confirm the existence of hydrogen bonding between the invariant Asn51 and the second adenine of the DNA binding sequence, as seen in crystal structures of other homeodomain-DNA complexes, but never before detected by NMR. Hydrogen bonding by Arg5 and Lys3 in the minor groove of the DNA appears to be responsible for two unusually upfield-shifted ribose H1' resonances. The DNA duplex is nearly straight and its structure is primarily that of B -DNA. A detailed comparison is presented for all available homeodomain-DNA structures including the vnd/NK-2 DNA complex, which demonstrates that homology is maintained in the protein structure, whereas for the orientation of the homeodomain relative to DNA, small but significant variations are observed. Interactions are described involving certain residues in specific positions of the homeodomain, namely Leu7, Thr41, and Gln50 of vnd/NK-2, where single amino acid residue mutations lead to dramatic developmental alterations. The availability of our previously determined three- dimensional structure of the vnd/NK-2 homeodomain in the absence of DNA allows us to assess structural changes in the homeodomain induced by DNA binding., (Copyright 1999 Academic Press.)

This work is the first in a series devoted to applying mode coupling diffusion theory to the derivation of local dynamics properties of proteins in solution. The first-order mode-coupling approximation, or optimized Rouse-Zimm local dynamics (ORZLD), is applied here to derive the rotational dynamics of the bonds and compare the calculated with the experimental nmr 15N spin-lattice relaxation time behavior of the vnd/NK-2 homeodomain from Drosophila melanogaster. The starting point for the calculations is the experimental three-dimensional structure of the homeodomain determined by multidimensional nmr spectroscopy. The results of the computations are compared with experimentally measured 15N spin-lattice relaxation times T1, at 34.5 and 60.8 MHz, to check the first-order approximation. To estimate the relative importance of internal and overall rotation, both rigid and fluctuating dynamic models are examined, with fluctuations evaluated using molecular dynamics (MD) simulations. The correlation times for the fundamental bond vector time correlation function and for the second-order bond orientational TCF are obtained as a function of the residue number for vnd/NK-2. The stability of the corresponding local dynamics pattern for the fluctuating structure as a function of the length of the MD trajectory is presented. Diffusive dynamics, which is essentially free of model parameters even at first order in the mode-coupling diffusion approach, confirm that local dynamics of proteins can be described in terms of rotational diffusion of a fluctuating quasi-rigid structure. The comparison with the nmr data shows that the first-order mode coupling diffusion approximation accounts for the correct order of magnitude of the results and of important qualitative aspects of the data sensitive to conformational changes. Indications are obtained from this study to efficiently extend the theory to higher order in the mode-coupling expansion. These results demonstrate the promise of the mode-coupling approach, where the local dynamics of proteins is described in terms of rotational diffusion of a fluctuating quasi-rigid structure, to analyze nmr spin-lattice relaxation behavior.

Two pulse sequences are described that employ a modified water flipback technique to enhance the signal intensity of weak side chain resonances at the protein-DNA interface of the vnd/NK-2 homeodomain/DNA complex in an 15N-edited three-dimensional NOESY-HMQC spectrum. The pulse sequences presented employ water flipback pulses at the beginning of the NOESY mixing time, optimizing the direct NOE transfer of magnetization from the water to the protein by maximizing the z-component of the water magnetization. In one of the pulse sequences, radiation damping during the the indirect 1H and 15N evolution times is suppressed. A modified version of the WATERGATE water suppression technique is employed during the HMQC portion of the experiment. The signal enhancement is demonstrated for the resonances of the side chain amide of Asn51, an invariant homeodomain residue whose contact with the DNA is critical for binding. An ancillary advantage of the experiment is the ability to observe NOE transfer of magnetization from water. The information present in the water resonance plane of the three-dimensional spectrum is illustrated in a comparison with the corresponding HMQC spectrum of the protein/DNA complex., (Copyright 1998 Academic Press.)

The structural and DNA binding behavior is described for an analog of the vnd/NK-2 homeodomain, which contains a single amino acid residue alanine to threonine replacement in position 35 of the homeodomain. Multidimensional nuclear magnetic resonance, circular dichroism, and electrophoretic gel retardation assays were carried out on recombinant 80-aa residue proteins that encompass the wild-type and mutant homeodomains. The mutant A35T vnd/NK-2 homeodomain is unable to adopt a folded conformation free in solution at temperatures down to -5 degreesC in contrast to the behavior of the corresponding wild-type vnd/NK-2 homeodomain, which is folded into a functional three-dimensional structure below 25 degreesC. The A35T vnd/NK-2 binds specifically to the vnd/NK-2 target DNA sequence, but with an affinity that is 50-fold lower than that of the wild-type homeodomain. Although the three-dimensional structure of the mutant A35T vnd/NK-2 in the DNA bound state shows characteristic helix-turn-helix behavior similar to that of the wild-type homeodomain, a notable structural deviation in the mutant A35T analog is observed for the amide proton of leucine-40. The wild-type homeodomain forms an unusual i,i-5 hydrogen bond with the backbone amide oxygen of residue 35. In the A35T mutant this amide proton resonance is shifted upfield by 1.27 ppm relative to the resonance frequency for the wild-type analog, thereby indicating a significant alteration of this i,i-5 hydrogen bond.

Secondary structures, DNA binding properties, and thermal denaturation behavior of six site-directed mutant homeodomains encoded by the vnd/NK-2 gene from Drosophila melanogaster are described. Three single site H52R, Y54M, and T56W mutations, two double site H52R/T56W and Y54M/T56W mutations, and one triple site H52R/Y54M/T56W mutation were investigated. These positions were chosen based on their variability across homeodomains displaying differences in secondary structure and DNA binding specificity. Multidimensional NMR, electrophoretic mobility shift assays, and circular dichroism spectropolarimetry studies were carried out on recombinant 80-amino acid residue proteins containing the homeodomain. Position 56, but more importantly position 56 in combination with position 52, plays an important role in determining the length of the recognition helix. The H52R mutation alone does not affect the length of this helix but does increase the thermal stability. Introduction of site mutations at positions 52 and 56 in vnd/NK-2 does not modify their high affinity binding to the 18-base pair DNA fragment containing the vnd/NK-2 consensus binding sequence, CAAGTG. Site mutations involving position 54 (Y54M, Y54M/T56W, and H52R/Y54M/T56W) all show a decrease of 1 order of magnitude in their binding affinity. The roles in structure and sequence specificity of individual atom-atom interactions are described.

The interactions responsible for the nucleotide sequence-specific binding of the vnd/NK-2 homeodomain of Drosophila melanogaster to its consensus DNA binding site have been identified. A three-dimensional structure of the vnd/NK-2 homeodomain-DNA complex is presented, with emphasis on the structure of regions of observed protein-DNA contacts. This structure is based on protein-DNA distance restraints derived from NMR data, along with homology modeling, solvated molecular dynamics, and results from methylation and ethylation interference experiments. Helix III of the homeodomain binds in the major groove of the DNA and the N-terminal arm binds in the minor groove, in analogy with other homeodomain-DNA complexes whose structures have been reported. The vnd/NK-2 homeodomain recognizes the unusual DNA consensus sequence 5'-CAAGTG-3'. The roles in sequence specificity and strength of binding of individual amino acid residues that make contact with the DNA are described. We show, based primarily on the observed protein-DNA contacts, that the interaction of Y54 with the DNA is the major determinant of this uncommon nucleotide binding specificity in the vnd/NK-2 homeodomain-DNA complex.

We describe the NMR determination of the three-dimensional structure of a 77 amino acid residue protein, which consists of the 60 residue NK-2 homeodomain from Drosophila melanogaster and adjacent amino acid residues. The NK-2 homeodomain protein is part of a 723 amino acid residue protein which is expressed early in embryonic development in part of the central nervous system. NK-2 was characterized using both a natural abundance and a uniformly 15N enriched sample by two-dimensional and three-dimensional NMR experiments. The average root-mean-square deviation for 30 structures for residues 8 to 53 is 0.40 A for the backbone heavy-atoms and 0.72 A for the backbone and side-chain heavy-atoms. These structures were obtained from 986 NOE-derived upper and lower bound restraints. The three-dimensional structure contains three helices which consist of homeodomain amino acid residues 10 to 22, 28 to 38 and 42 to 52, as well as a turn between helix II and III, characteristic of homeodomains. Residues 53 to 60 of the DNA recognition helix are not fully ordered in the absence of DNA. In the free state this segment adopts a flexible but helix-like structure between residues 53 and 56 and is disordered from residues 57 to 60 although, as shown previously, the helix elongates by eight residues upon binding to DNA. The role of variable residues 52, 54 and 56 in determining the structure and flexibility of the recognition helix, as well as the stability of the NK-2 homeodomain as manifested by its thermal denaturation, are discussed.

The secondary structure of the homeodomain encoded by the NK-2 gene from Drosophila melanogaster, in both the free and DNA-bound states, was determined in solution using two- and three-dimensional (2D and 3D) NMR spectroscopy. Proton and 15N studies were carried out on a 77 amino acid residue protein that contains the homeodomain, which was synthesized in Escherichia coli. On the basis of NOE connectivities, vicinal coupling constants, and proton-deuterium exchange behavior, three helical segments were found that consist of homeodomain amino acid residues 10-22, 28-38, and 42-52 for the protein in the absence of DNA. The major structural differences between free NK-2 and other homeodomains are the increased internal mobility of the second helix and the shorter length of the third helix, also termed the recognition helix. Despite this shorter helix, NK-2 exhibits high-affinity binding to DNA compared to other homeodomains (kD = 2.0 x 10(-10) M; L.-H. Wang and M. Nirenberg, unpublished results). The formation of the complex of NK-2 with the duplex DNA (TGTGTCAAGTG-GCTGT) significantly increases the thermal stability of the protein. The Tm increases from 25 degrees C (free NK-2) to > 47 degrees C (DNA-bound NK-2). Also, a dramatic increase in the length of helix III is observed. In the absence of DNA, the DNA recognition helix is 11 amino acid residues long (residues 42-52), whereas in the presence of DNA, the length of this helix extends to 19 amino acids (residues 42-60).(ABSTRACT TRUNCATED AT 250 WORDS)