The role for glutamic acid at position 196 in human hypoxanthine phosphoribosyltransferase (HPRT) as investigated using site-directed mutagenesis

The crystal structure of human HPRT reveals the involvement of E196 side chain at the A-B dimer interface. Interference by valine substitution at this position (E196V), as identified in patients with Lesch-Nyhan disease, nearly abolishes enzymatic activity. Kinetic analysis of the active mutants (E196A, E196D, E196Q, and E196R) suggests that interaction between K68 and E196 side chains contributes to stabilization of cis-configuration during the catalytic cycle. The study also provides further insight into the role of A-B dimer interactions relating to K68 in the regulation of cis-trans isomerization that potentially governs the rate-limiting steps in the HPRT reaction.     

Interactions at the dimer interface influence the relative efficiencies for purine nucleotide synthesis and pyrophosphorolysis in a phosphoribosyltransferase

Enzymes that salvage 6-oxopurines, including hypoxanthine phosphoribosyltransferases (HPRTs), are potential targets for drugs in the treatment of diseases caused by protozoan parasites. For this reason, a number of high-resolution X-ray crystal structures of the HPRTs from protozoa have been reported. Although these structures did not reveal why HPRTs need to form dimers for catalysis, they revealed the existence of potentially relevant interactions involving residues in a loop of amino acid residues adjacent to the dimer interface, but the contributions of these interactions to catalysis remained poorly understood. The loop, referred to as active-site loop I, contains an unusual non-proline cis-peptide and is composed of residues that are structurally analogous with Leu67, Lys68, and Gly69 in the human HPRT. Functional analyses of site-directed mutations (K68D, K68E, K68N, K68P, and K68R) in the HPRT from Trypanosoma cruzi, etiologic agent of Chagas' disease, show that the side-chain at position 68 can differentially influence the K(m) values for all four substrates as well as the k(cat) values for both IMP formation and pyrophosphorolysis. Also, the results for the K68P mutant are inconsistent with a cis-trans peptide isomerization-assisted catalytic mechanism. These data, together with the results of structural studies of the K68R mutant, reveal that the side-chain of residue 68 does not participate directly in reaction chemistry, but it strongly influences the relative efficiencies for IMP formation and pyrophosphorolysis, and the prevalence of lysine at position 68 in the HPRT of the majority of eukaryotes is consistent with there being a biological role for nucleotide pyrophosphorolysis.     

Pressure dissociation studies provide insight into oligomerization competence of temperature-sensitive folding mutants of P22 tailspike

Several temperature-sensitive folding (tsf) mutants of the tailspike protein from bacteriophage P22 have been found to fold with lower efficiency than the wild-type sequence, even at lowered temperatures. Previous refolding studies initiated from the unfolded monomer have indicated that the tsf mutations decrease the rate of structured monomer formation. We demonstrate that pressure treatment of the tailspike aggregates provides a useful tool to explore the effects of tsf mutants on the assembly pathway of the P22 tailspike trimer. The effects of pressure on two different tsf mutants, G244R and E196K, were explored. Pressure treatment of both G244R and E196K aggregates produced a folded trimer. E196K forms almost no native trimer in in vitro refolding experiments, yet it forms a trimer following pressure in a manner similar to the native tailspike protein. In contrast, trimer formation from pressure-treated G244R aggregates was not rapid, despite the presence of a G244R dimer after pressure treatment. The center-of-mass shifts of the fluorescence spectra under pressure are nearly identical for both tsf aggregates, indicating that pressure generates similar intermediates. Taken together, these results suggest that E196K has a primary defect in formation of the beta-helix during monomer collapse, while G244R is primarily an assembly defect.     

Different misfolding mechanisms converge on common conformational changes

Prion diseases are caused by misfolding and aggregation of the prion protein (PrP). Pathogenic mutations such as Y218N and E196K are known to cause Gerstmann-Sträussler-Scheinker syndrome and Creutzfeldt-Jakob disease, respectively. Here we describe molecular dynamics simulations of these mutant proteins to better characterize the detailed conformational effects of these sequence substitutions. Our results indicate that the mutations disrupt the wild-type native PrP^C structure and cause misfolding. Y218N reduced hydrophobic packing around the X-loop (residues 165–171), and E196K abolished an important wild-type salt bridge. While differences in the mutation site led PrP mutants to misfold along different pathways, we observed multiple traits of misfolding that were common to both mutants. Common traits of misfolding included: 1) detachment of the short helix (HA) from the PrP core; 2) exposure of side chain F198; and 3) formation of a nonnative strand at the N-terminus. The effect of the E196K mutation directly abolished the wild-type salt bridge E196-R156, which further destabilized the F198 hydrophobic pocket and HA. The Y218N mutation propagated its effect by increasing the HB-HC interhelical angle, which in turn disrupted the packing around F198. Furthermore, a nonnative contact formed between E221 and S132 on the S1-HA loop, which offered a direct mechanism for disrupting the hydrophobic packing between the S1-HA loop and HC. While there were common misfolding features shared between Y218N and E196K, the differences in the orientation of HB and HC and the X-loop conformation might provide a structural basis for identifying different prion strains.     

Cooperative fluctuations point to the dimerization interface of p53 core domain

Elastic network models are used for investigation of the p53 core domain functional dynamics. Global modes of motion indicate high positive correlations for residue fluctuations across the A-B interface, which are not observed at the B-C interface. Major hinge formation is observed at the A-B interface upon dimerization indicating stability of the A-B dimer. These findings imply A-B as the native dimerization interface, whereas B-C is the crystal interface. The A-B dimer exhibits an opening-closing motion about DNA, supporting the previously suggested clamp-like model of nonspecific DNA binding followed by diffusion. Monomer A has limited positive correlations with DNA, while monomer B exhibits high positive correlations with DNA in the functionally significant slow modes. Thus, monomer B might seem to maintain the stability of the dimer-DNA complex by forming the relatively fixed arm of the dimer clamp, whereas the other arm of the clamp, monomer A, might allow sliding via continuous association/dissociation mechanisms.     

Structural and functional studies of the human phosphoribosyltransferase domain containing protein 1

Human hypoxanthine-guanine phosphoribosyltransferase (HPRT) (EC 2.4.2.8) catalyzes the conversion of hypoxanthine and guanine to their respective nucleoside monophosphates. Human HPRT deficiency as a result of genetic mutations is linked to both Lesch-Nyhan disease and gout. In the present study, we have characterized phosphoribosyltransferase domain containing protein 1 (PRTFDC1), a human HPRT homolog of unknown function. The PRTFDC1 structure has been determined at 1.7 ? resolution with bound GMP. The overall structure and GMP binding mode are very similar to that observed for HPRT. Using a thermal-melt assay, a nucleotide metabolome library was screened against PRTFDC1 and revealed that hypoxanthine and guanine specifically interacted with the enzyme. It was subsequently confirmed that PRTFDC1 could convert these two bases into their corresponding nucleoside monophosphate. However, the catalytic efficiency (k(cat)/K(m)) of PRTFDC1 towards hypoxanthine and guanine was only 0.26% and 0.09%, respectively, of that of HPRT. This low activity could be explained by the fact that PRTFDC1 has a Gly in the position of the proposed catalytic Asp of HPRT. In PRTFDC1, a water molecule at the position of the aspartic acid side chain position in HPRT might be responsible for the low activity observed by acting as a weak base. The data obtained in the present study indicate that PRTFDC1 does not have a direct catalytic role in the nucleotide salvage pathway.     

Solution structure of the chromomycin-DNA complex

The structure of the chromomycin-DNA complex at the deoxyoctanucleotide duplex level has been determined from one- and two-dimensional proton NMR studies in Mg-containing aqueous solution. The NMR results demonstrate that the antitumor agent binds as a symmetrical dimer to the self-complementary d[T-T-G-G-C-C-A-A] duplex with retention of the 2-fold symmetry in the complex. A set of intermolecular nuclear Overhauser enhancements (NOEs) establishes that two chromomycin molecules in the dimer share the minor groove at the G-G-C-C.G-G-C-C segment in such a way that each hydrophilic edge of the chromophore is located next to the G-G.C-C half-site and each C-D-E trisaccharide chain extends toward the 3'-direction of the octanucleotide duplex. In addition, the A-B disaccharide segment and the hydrophilic side chain of the antitumor agent are directed toward the phosphate backbone. The observed changes in nucleic acid NOEs and coupling patterns on complex formation establish a transition to a wider and shallower minor groove at the central G-G-C-C.G-G-C-C segment required for accommodating the chromomycin dimer. The present demonstration that chromomycin binds as a dimer and switches the conformation of the DNA at its G.C-rich minor groove binding site provides new insights into antitumor agent design and the sequence specificity of antitumor agent-DNA recognition.     

The effect of valine substitution for glycine in the dimer interface of citrate synthase from Thermoplasma acidophilum on stability and activity

To determine the role of hydrophobic interactions in the dimer interface of citrate synthase (CS) from Thermoplasma (Tp) acidophilum in thermostabilization, we have used site-directed mutagenesis to replace Gly 196 by Val on the helix L of the subunit interface. Recombinant wild-type and Gly 196 mutant TpCS enzymes were largely identical in terms of substrate specificities (K(m) for oxaloacetate and acetyl CoA). However, the mutation not only reduced catalytic activity (about 10-fold) (i.e., V(max), k(cat) and specific activity) of the TpCS, but also decreased its thermal and chemical stability. Archaeal citrate synthase is active as a dimer, since residues from both monomers participate in the active site. Our results suggest that Gly196 --> Val mutation interferes with dimerization, so that improper dimerization or dissociation of the dimer would have a profound affect on the activity as well as the conformational stability of TpCS.     

Crystal structures of free,IMP-, and GMP-bound Escherichia coli hypoxanthine phosphoribosyltransferase

Crystal structures have been determined for free Escherichia coli hypoxanthine phosphoribosyltransferase (HPRT) (2.9 A resolution) and for the enzyme in complex with the reaction products, inosine 5'-monophosphate (IMP) and guanosine 5'-monophosphate (GMP) (2.8 A resolution). Of the known 6-oxopurine phosphoribosyltransferase (PRTase) structures, E. coli HPRT is most similar in structure to that of Tritrichomonas foetus HGXPRT, with a rmsd for 150 Calpha atoms of 1.0 A. Comparison of the free and product bound structures shows that the side chain of Phe156 and the polypeptide backbone in this vicinity move to bind IMP or GMP. A nonproline cis peptide bond, also found in some other 6-oxopurine PRTases, is observed between Leu46 and Arg47 in both the free and complexed structures. For catalysis to occur, the 6-oxopurine PRTases have a requirement for divalent metal ion, usually Mg(2+) in vivo. In the free structure, a Mg(2+) is coordinated to the side chains of Glu103 and Asp104. This interaction may be important for stabilization of the enzyme before catalysis. E. coli HPRT is unique among the known 6-oxopurine PRTases in that it exhibits a marked preference for hypoxanthine as substrate over both xanthine and guanine. The structures suggest that its substrate specificity is due to the modes of binding of the bases. In E. coli HPRT, the carbonyl oxygen of Asp163 would likely form a hydrogen bond with the 2-exocyclic nitrogen of guanine (in the HPRT-guanine-PRib-PP-Mg(2+) complex). However, hypoxanthine does not have a 2-exocyclic atom and the HPRT-IMP structure suggests that hypoxanthine is likely to occupy a different position in the purine-binding pocket.     

Rare E196A mutation in PRNP gene of 3 Chinese patients with Creutzfeldt-Jacob disease

Inherited prion diseases are characterized by mutations in the PRNP gene, which account for 5–15% of human prion diseases. Here we reported 3 Chinese genetic Creutzfeldt-Jacob disease cases (gCJD) with a rare mutation in PRNP leading to an exchange of amino acid from glutamic acid (E) to alanine (A) at codon 196 (E196A). All three patients were Han Chinese without any sibship among them. They showed various unspecific symptoms at onset and displayed typical clinical manifestations of sporadic CJD with progress of disease. The same time, 2 cases showed psychotic symptoms during the clinical courses. 14-3-3 proteins were positive in cerebrospinal fluid (CSF) and special abnormality were detected in MRI of all the cases. The polymorphism of codon 129 was methionin homozygote and that of codon 219 was glutamate homozygote in all 3 patients. The disease durations of the 3 cases varied from 10 to 22 months and no disease associated family history was figured out in all the cases.     

Synthesis of a hexasaccharide fragment of group E streptococci polysaccharide and the tetrasaccharide repeating unit of E. coli O7:K98:H6

Syntheses of a hexasaccharide, the dimer of the repeating unit of the group E streptococci polysaccharide, and a tetrasaccharide, the repeating unit of the E. coli O7:K98:H6, were achieved by constructing alternate alpha-L-(1-->2)- and alpha-L-(1-->3)-linked L-rhamnopyranose backbones and substituting with beta-linked D-glucopyranose side chains for the former, and a D-glucopyranosyluronate branch for the latter, respectively, at O-2 of the L-rhamnose ring.     

Ionization of His 55 at the dimer interface of dynein light-chain LC8 is coupled to dimer dissociation

LC8 is a highly conserved light-chain subunit of cytoplasmic dynein that interacts with a wide variety of cellular proteins and is presumed to play a fundamental role in dynein assembly and cargo recruitment and in the assembly of protein complexes unrelated to dynein. LC8 is a dimer at physiological pH but dissociates to a folded monomer at pH < 4.8. We have suggested that acid-induced dimer dissociation is due to protonation of His 55, which is stacked against His 55' and completely buried in the dimer interface. In this work, we show that the pH-induced dissociation is reversible and indeed governed by the ionization state of His 55. Mutagenesis of His 55 to Lys results in a monomer in the pH range of 3-8, while the mutation to Ala results in a dimer in the same pH range. Mutations that disrupt intermolecular hydrogen bonds between Tyr 65 and Lys 44' and His 55 and Thr 67' do not change the association state of the dimer. Titration curves for His 55 and the two other histidines, His 72 and 68, were determined by (13)C-(1)H NMR for H55K and for WT-LC8 in the monomeric and dimeric states. The pK(a) values of His 72 and His 68 are 6 in the WT dimer and 6.2-6.5 in monomeric H55K, while the pK(a) of His 55 is about 4.5 in the WT dimer. These results indicate that deprotonation of His 55 is linked to dimer formation and that mutation of His 55 to a small neutral residue or to a positively charged residue uncouples the protonation and dissociation processes.     

Dominant Intermediate Charcot-Marie-Tooth disorder is not due to a catalytic defect in tyrosyl-tRNA synthetase

Charcot-Marie-Tooth disorder (CMT) is the most common inherited peripheral neuropathy, afflicting 1 in every 2500 Americans. One form of this disease, Dominant Intermediate Charcot-Marie-Tooth disorder type C (DI-CMTC), is due to mutation of the gene encoding the cytoplasmic tyrosyl-tRNA synthetase (TyrRS). Three different TyrRS variants have been found to give rise to DI-CMTC: replacing glycine at position 41 by arginine (G41R), replacing glutamic acid at position 196 by lysine (E196K), and deleting amino acids 153-156 (Δ(153-156)). To test the hypothesis that DI-CMTC is due to a defect in the ability of tyrosyl-tRNA synthetase to catalyze the aminoacylation of tRNA(Tyr), we have expressed each of these variants as recombinant proteins and used single turnover kinetics to characterize their abilities to catalyze the activation of tyrosine and its subsequent transfer to the 3' end of tRNA(Tyr). Two of the variants, G41R and Δ(153-156), display a substantial decrease in their ability to bind tyrosine (>100-fold). In contrast, the E196K substitution does not significantly affect the kinetics for formation of the tyrosyl-adenylate intermediate and actually increases the rate at which the tyrosyl moiety is transferred to tRNA(Tyr). The observation that the E196K substitution does not decrease the rate of catalysis indicates that DI-CMTC is not due to a catalytic defect in tyrosyl-tRNA synthetase.     

The role of conserved arginine residue in loop 4 of glycoside hydrolase family 10 xylanases

An arginine residue in loop 4 connecting beta strand 4 and alpha-helix 4 is conserved in glycoside hydrolase family 10 (GH10) xylanases. The arginine residues, Arg(204) in xylanase A from Bacillus halodurans C-125 (XynA) and Arg(196) in xylanase B from Clostridium stercorarium F9 (XynB), were replaced by glutamic acid, lysine, or glutamine residues (XynA R204E, K and Q, and XynB R196E, K and Q). The pH-k(cat)/K(m) and the pH-k(cat) relationships of these mutant enzymes were measured. The pK(e2) and pK(es2) values calculated from these curves were 8.59 and 8.29 (R204E), 8.59 and 8.10 (R204K), 8.61 and 8.19 (R204Q), 7.42 and 7.19 (R196E), 7.49 and 7.18 (R196K), and 7.86 and 7.38 (R196Q) respectively. Only the pK(es2) value of arginine derivatives was less than those of the wild types (8.49 and 9.39 [XynA] and 7.62 and 7.82 [XynB]). These results suggest that the conserved arginine residue in GH10 xylanases increases the pK(a) value of the proton donor Glu during substrate binding. The arginine residue is considered to clamp the proton donor and subsite +1 to prevent structural change during substrate binding.     

Cysteine pKa depression by a protonated glutamic acid in human DJ-1

Human DJ-1, a disease-associated protein that protects cells from oxidative stress, contains an oxidation-sensitive cysteine (C106) that is essential for its cytoprotective activity. The origin of C106 reactivity is obscure, due in part to the absence of an experimentally determined p K a value for this residue. We have used atomic-resolution X-ray crystallography and UV spectroscopy to show that C106 has a depressed p K a of 5.4 +/- 0.1 and that the C106 thiolate accepts a hydrogen bond from a protonated glutamic acid side chain (E18). X-ray crystal structures and cysteine p K a analysis of several site-directed substitutions at residue 18 demonstrate that the protonated carboxylic acid side chain of E18 is required for the maximal stabilization of the C106 thiolate. A nearby arginine residue (R48) participates in a guanidinium stacking interaction with R28 from the other monomer in the DJ-1 dimer and elevates the p K a of C106 by binding an anion that electrostatically suppresses thiol ionization. Our results show that the ionizable residues (E18, R48, and R28) surrounding C106 affect its p K a in a way that is contrary to expectations based on the typical ionization behavior of glutamic acid and arginine. Lastly, a search of the Protein Data Bank (PDB) produces several candidate hydrogen-bonded aspartic/glutamic acid-cysteine interactions, which we propose are particularly common in the DJ-1 superfamily.     

Frequency of reactivation and variability in expression of X-linked enzyme loci.

A mouse-human cell hybrid clone retaining an inactive human X chromosome was treated with 5-azacytidine. Following treatment, expression of the X-linked enzyme markers, hypoxanthine-guanine phosphoribosyltransferase (HPRT), glucose-6-phosphate dehydrogenase (G6PD), phosphoglycerate kinase (PGK), and alpha-galactosidase A (GLA) was examined. Results presented here show that 45 of the 62 clones positive for human HPRT expressed human GLA, while only four of 68 clones negative for human HPRT expressed human GLA. These results strongly suggest that there is coordinate reactivation of GLA and HPRT. Reactivated expression of G6PD was studied in detail. The studies show that 5-azacytidine can induce heritable changes in the inactive human X chromosome resulting in the expression of G6PD activity at a level lower than that from an active human X chromosome.     

Effect of Asp69 and Arg310 on the pK of His68, a key catalytic residue of adenylosuccinate lyase

Adenylosuccinate lyase (ASL) of Bacillus subtilis contains three conserved histidines, His(68), His(89), and His(141), identified by affinity labeling and site-directed mutagenesis as critical to the intersubunit catalytic site. The pH-V(max) profile for wild-type ASL is bell-shaped (pK (1) = 6.74 and pK (2) = 8.28). Only the alkaline side changes with temperature, characteristic of histidine pKs. To identify determinants of pK (2) in the enzyme-substrate complex, we replaced residues at two positions close to His(68) (but not to His(89) or His(141)) in the structure. Compared with the specific activity of 1.75 mumol adenylosuccinate reacting/min/mg of wild-type enzyme at pH 7.0, mutant enzymes D69E, D69N, R310Q, and R310K exhibit specific activities of 0.40, 0.04, 0.00083, and 0.10, respectively. While D69E has a K (m) for adenylosuccinate similar to that of wild-type ASL, D69N and R310K exhibit modest increases in K (m), and R310Q has an 11-fold increase in K (m). The mutant enzymes show no significant change in molecular weight or secondary structure. The major change is in the pH-V(max) profile: pK (2) is 8.48 for the D69E mutant and is decreased to 7.83 in D69N, suggesting a proximal negative charge is needed to maintain the high pK of 8.28 observed for wild-type enzyme and attributed to His(68). Similarly, R310Q exhibits a decrease in its pK (2) (7.33), whereas R310K shows little change in pK (2) (8.24). These results suggest that Asp(69) interacts with His(68), that Arg(310) interacts with and orients the beta-carboxylate of Asp(69), and that His(68) must be protonated for ASL to be active.     

The dimeric dihydroorotate dehydrogenase A from Lactococcus lactis dissociates reversibly into inactive monomers

The flavoenzyme dihydroorotate dehydrogenase A from Lactococcus lactis is a homodimeric protein of 311 residues/subunit, and the two active sites are positioned at a distance from the dimer interface. To promote formation of the monomeric form of the enzyme, we changed the residues involved in formation of two salt bridges formed between the residues Glu206 of the one polypeptide and Lys296 of the other polypeptide. The mutant enzymes formed inactive precipitates when cells were grown at 37 degrees C, but remained soluble and active when cells were grown at 25 degrees C. The salt bridges were not needed for activity, because the mutant enzymes in which one of the residues was converted to an alanine (E206A or K296A) retained almost full activity. The mutant enzymes in which the charge of one of the residues of the salt bridge was inverted (i.e., E206K or K296E) were severely impaired. The double mutant E206K-K296E, which has the possibility of forming salt bridges in the opposite orientation of the wild type, was fully active in concentrated solutions, but dissociated into inactive monomers upon dilution. The K(D) for the dimer to monomer dissociation reaction was 12 microM, and dimer formation was favored by the product, orotate, or by high ionic strength, indicating that the hydrophobic interactions are important for the subunit contacts. Wild-type dihydroorotate dehydrogenase A was similarly found to dissociate into inactive monomers, but with a K(D) for dissociation equal to 0.12 microM. These results imply that the dimeric state is necessary for activity of the enzyme.     

Dication and trication which can increase the permeability of Escherichia coli outer membrane

Recent success in the preparation of the monomer, dimer and trimer in compound 48/80 prompted us to investigate the action of these compounds on Escherichia coli cells. It was found that compound 48/80 inhibited growth of E. coli cells, while the monomer, dimer and trimer in 48/80 did not. However, the following experiments showed that the dimer and trimer disrupted the permeability barrier of the outer membrane of E. coli. First, addition of the dimer or trimer in cell suspension stimulated the uptake of tetraphenylphosphonium cation. Second, the synergistic effect of the dimer on the action of gramicidin caused the efflux of K+. In experiments using isolated cytoplasmic membrane vesicles, addition of gramicidin alone caused the efflux of K+. Thus, it was speculated that, with whole cells, the dimer formed some defect structure in the outer membrane, through which gramicidin reached the cytoplasmic membrane and increased the K+ permeability. The temperature dependence of efflux K+ showed that the dimer in 48/80 rendered the outer membrane permeable to gramicidin at temperatures above the phase transition of the outer membrane.     

Mechanism of transfer of reduced nicotinamide adenine dinucleotide among dehydrogenases. Transfer rates and equilibria with enzyme-enzyme complexes

The direct transfer of NADH between A-B pairs of dehydrogenases and also the dissociation of NADH from individual E-NADH complexes have been investigated by transient stopped-flow kinetic techniques. Such A-B transfers of NADH occur without the intermediate dissociation of coenzyme into the aqueous solvent environment [Srivastava, D.K., & Bernhard, S.A. (1985) Biochemistry 24, 623-628]. The equilibrium distributions of limiting NADH among aqueous solvent and A and B dehydrogenase sites have also been determined. At sufficiently high but realizable concentrations of dehydrogenases, both the transfer rate and the equilibrium distribution of bound NADH are virtually independent of the excessive enzyme concentrations; at excessive E2 concentration, substantial NADH is bound to the E1 site. These results further substantiate earlier kinetic arguments for the preferential formation of an EA-NADH-EB complex, within which coenzyme is directly transferred between sites. The unimolecular specific rates of coenzyme transfer from site to site are nearly invariant among different A-B dehydrogenase pairs. The equilibrium constants for the distribution of coenzyme within the EA X EB complexes are near unity. At high [E2] and for [E2] greater than [E1] greater than [NADH], E1-NADH X E2 and E1 X NADH-E2 are virtually the only coenzyme-contained species. In contrast to the nearly invariant unimolecular NADH transfer rates within EA X EB complexes, unimolecular specific rates of dissociation of NADH from E-NADH into aqueous solution are highly variable.(ABSTRACT TRUNCATED AT 250 WORDS)