Originally published In Press as doi:10.1074/jbc.M201042200 on April 8, 2002
J. Biol. Chem., Vol. 277, Issue 25, 22863-22874, June 21, 2002
Crystal Structure of a Dual Activity IMPase/FBPase
(AF2372) from Archaeoglobus fulgidus
THE STORY OF A MOBILE LOOP*
Kimberly A.
Stieglitz
§,
Kenneth A.
Johnson¶
,
Hongying
Yang
,
Mary F.
Roberts
,
Barbara A.
Seaton**,
James F.
Head**, and
Boguslaw
Stec¶
From the
Department of Chemistry, Merkert Chemistry
Center, Boston College, Chestnut Hill, Massachusetts 02467, the
** Department of Physiology and Biophysics, Boston University
School of Medicine, Boston, Massachusetts 02118, and the
¶ Department of Biochemistry and Cell Biology, W. M. Keck Center
for Computational Biology, Rice University, Houston, Texas 77005
Received for publication, January 31, 2002, and in revised form, April 5, 2002
 |
ABSTRACT |
Several hyperthermophilic organisms contain an
unusual phosphatase that has dual activity toward inositol
monophosphates and fructose 1,6-bisphosphate. The structure of the
second member of this family, an FBPase/IMPase from
Archaeoglobus fulgidus (AF2372), has been solved. This
enzyme shares many kinetic and structural similarities with that of a
previously solved enzyme from Methanococcus jannaschii
(MJ0109). It also shows some kinetic differences in divalent metal ion
binding as well as structural variations at the dimer interface that
correlate with decreased thermal stability. The availability of
different crystal forms allowed us to investigate the effect of the
presence of ligands on the conformation of a mobile catalytic loop
independently of the crystal packing. This conformational variability
in AF2372 is compared with that observed in other members of this
structural family that are sensitive or insensitive to submillimolar
concentrations of Li+. This analysis provides support for
the previously proposed mechanism of catalysis involving three metal
ions. A direct correlation of the loop conformation with strength of
Li+ inhibition provides a useful system of classification
for this extended family of enzymes.
 |
INTRODUCTION |
In mammalian cells, inositol monophosphatases are abundant
cytosolic enzymes necessary to regenerate the supply of
myo-inositol for the synthesis of phosphatidylinositols.
However, in archaeal hyperthermophilic organisms, e.g.
Archaeoglobus fulgidus, an
IMPase1 ortholog is thought
to be used in the biosynthesis of a unique osmolyte,
di-myo-inositol 1,1-phosphate (1). The archaeal
enzymes have another twist as well. Besides catalyzing the hydrolysis of Ins-1-P (needed in response to thermal and salt stress) they also
very specifically dephosphorylate FBP at C-1 (2). Enzymes with these
characteristics exist as two distinct gene products in mammalian
systems and several bacteria (3). The dual specificity of the enzyme
from primitive hyperthermophiles strongly suggest that IMPase and
FBPase enzymes may have evolved from the same gene product (2).
Structural features of the IMPase/FBPase (MJ0109) from
Methanococcus jannaschii (2) are similar to those of
mammalian IMPase, inositol polyphosphase phosphatase (IPPase),
FBPase, and the yeast enzyme Hal2. These enzymes form an extended
family for which the general principles of chemical reactivity are the
same (3). Two or three metal ions are necessary for catalysis to occur. From the structure of human IMPase with bound Ins-1-P and
Gd3+, it was inferred that two metal ions were needed for
catalysis (4). However, crystallographic data for pig kidney
FBPase (5) and for MJ0109 (6) are consistent with a three-metal
ion-assisted catalytic mechanism.
There is strong biochemical evidence that Li+ inhibits
IMPase in mammalian cells where a dose response to this ion
correlates with a decrease in the intracellular pool of
myo-inositol (7). Li+, by inhibiting mammalian
IMPases, is thought to cause a depletion of inositol that attenuates
the synthesis of phosphatidylinositol necessary to regenerate
the phospholipase C dependent phosphatidylinositol/IP3 signaling pathway (8, 9). There are two other enzymes that are very
Li+-sensitive with known three-dimensional structures:
human IPPase (10) and yeast 3'-phosphoadenosine-5'-phosphate
phosphatase, Hal2 (11). In contrast, the dual specificity thermophilic
IMPase/FBPase enzymes from M. jannaschii and A. fulgidus and mammalian FBPase are not inhibited by Li+
in the submillimolar range (Table I).
Thus, the IMPase superfamily can be functionally subdivided into those
that are strongly inhibited by Li+ and those that are less
sensitive to Li+. Structural hints on the mechanism of
Li+ inhibition of human IMPase are indirect. In pig kidney
FBPase (13), which is inhibited by Li+ but activated by low
concentrations of K+, monovalent cations are thought to
inhibit the enzyme by distorting the active site geometry or preventing
turnover or product release. Despite extensive mutagenesis experiments
carried out on several enzymes in the family, the mechanism of
Li+ inhibition is still not entirely understood.
This work presents the results of structural studies for another dual
specificity IMPase/FBPase from the hyperthermophilic sulfate reducer
A. fulgidus. Like MJ0109, this archaeal enzyme (AF2372) has
sequence features close to human IMPase but key structural similarities more akin to pig kidney FBPase. A comparison of the structures of the two archaeal Li+-insensitive enzymes
(AF2372 and MJ0109), a phosphatase of intermediate sensitivity (pig
kidney FBPase), and three very Li+-sensitive members of
this superfamily (human IMPase, IPPase and yeast Hal2) provides a
strong correlation between the level of Li+ sensitivity and
the conformation of a critical catalytic loop. The disposition of this
catalytic loop may serve as a predictor of Li+ sensitivity
for other IMPase family members. Furthermore, changes in the loop
conformation of AF2372 induced by ligand binding shed some light into
the possible mechanism of Li+ inhibition in the human
IMPase and into how this phenomenon might have developed from the more
primitive archaeal enzymes.
 |
EXPERIMENTAL PROCEDURES |
Chemicals--
D,L-Ins-1-P, EDTA, sodium
chloride, magnesium chloride, zinc chloride, ammonium chloride,
malachite green oxalate, ammonium molybdate,
isopropyl-1-thio-
-D-galactopyranosideisopropyl-1-thio-
-D, SDS-PAGE molecular weight markers, gel filtration molecular weight markers Trizma base, and PEG 8000 were obtained from Sigma. Glycine, myo-inositol, and ammonium sulfate were obtained from
Aldrich. Ampicillin was obtained from Fisher. Tryptone and yeast
extract were obtained from Difco. Q-Sepharose fast flow,
phenyl-Sepharose, and Sephacryl S-200HR resins were obtained from
Pharmacia. L-Ins-1-P was synthesized enzymatically using an
Ins-1-P synthase cloned from A. fulgidus as described
previously (14). The Ins-1-P and the other small molecule fractions
were separated from the synthase with Centricon concentrators. The
filtrate was evaporated under high vacuum; a known amount of glucose
6-phosphate was used to quantify the concentration of purified
L-Ins-1-P by 31P NMR spectroscopy. Further
purification of the sugar required chromatography on an AG1-X8 anion
exchange column (elution with a gradient of 0.2-1.0 M
formic acid); fractions containing Ins-1-P were collected and the pH
value was adjusted to 7.5 with NH4OH. The purified
L-Ins-1-P was stored at 4 °C and used for enzymatic assays and crystal soaks. All other chemicals were reagent grade.
Overexpression and Purification of AF2372 from A. fulgidus--
Cloning the AF2372 gene from A. fulgidus and
overexpression of the gene product in Escherichia coli were
carried out as described elsewhere (15). Key steps of the purification
procedure included heating of crude protein at 85 °C for 30 min,
followed by chromatography on Q-Sepharose fast flow (elution with a
linear gradient of 0 to 0.5 M KCl). Column fractions were
monitored for IMPase activity using 15 mM pNPP and 30 mM MgCl2 in 50 mM Tris and 1 mM EDTA, pH 7.5. Purity of the fractions was assessed with
SDS-PAGE; fractions that were greater than 90% pure were collected,
pooled, concentrated 10-fold, and dialyzed against two changes of 50 mM Tris, pH 7.5. The concentrated solution was then loaded
onto a 1.6 × 72-cm Sephacryl S-200HR gel filtration column
pre-equilibrated with 50 mM Tris and 1 mM EDTA,
pH 7.5. Fractions containing pure Ins-1-Pase activity (as measured by
SDS-PAGE and silver staining of the gels) were concentrated to 12 mg/ml
and stored in 20 mM Tris HCl and 1.0 mM EDTA,
pH 8.0, until used for crystallization.
Activity Assays--
Production of inorganic phosphate from FBP
and Ins-1-P was monitored using a colorimetric assay (16). Substrate
concentration was typically 2 mM (except for pNPP, which
was 0.5 mM) in 50 mM Tris, pH 8.0, with varying
amounts of divalent cations (Mg2+ or Mn2+)
added. The solution with substrate and enzyme was heated to 85 or
80 °C (the temperature is indicated with each assay set) for 1 min,
then 900 µl of dye was added, and the OD660 was
measured. A standard Pi curve was used to convert
OD660 to µmol of Pi produced. Assays
typically used 1.0-2.0 µg of enzyme.
Crystallization--
AF2372 was crystallized by vapor diffusion
using hanging drops of 5 µl. The apoenzyme was crystallized in either
0.2 M ammonium nitrate and 30% PEG 3350 (P32
form) or 0.2 M dihydrogen ammonium phosphate and 30% PEG
3350 (P21 form). Enzyme solution (2.5 µl of ~10 mg/ml)
was mixed with an equal volume of mother liquor and equilibrated using
the vapor diffusion hanging drop method. To monitor the catalytic loop
conformation under various conditions, the apoenzyme crystals were
soaked in solutions containing ligands and different divalent cations.
In addition, the enzyme was co-crystallized with Ca2+ (50 mM) using 5-7% v/v solution of PEG 8000 in 50 mM Tris, pH 8.0, with 100 mM NaCl. These
crystals were then soaked with either D-Ins-1-P or FBP. The
drops were equilibrated against the reservoir with the concentration of
PEG 8000 between 5 and 7%. Crystals appeared in few days and grew to
full size within a week.
X-ray Data Collection--
Crystals having approximate
dimensions of 0.5 × 0.4 × 0.4 mm were mounted in glass
capillaries directly from the crystallization wells, or flash frozen in
a stream of nitrogen gas. Diffraction data were collected both at cryo
conditions and at room temperature. Data were collected at Boston
University School of Medicine at room temperature on a Rigaku RU-300
rotating anode generator with an R-AXIS IIC imaging plate detector.
Data were also collected at Rice University using a Rigaku H3R rotating
anode generator with an R-AXIS IV imaging plate detector with the
exception of the apo form in P21 space group that was
collected at the Brookhaven Synchrotron beamline X8C. Data collected in
Brookhaven and at Boston University was indexed and reduced using DENZO
and Scalepack (see Ref. 17), while the data collected at Rice
University was indexed and reduced using Crystal Clear software. In
addition to the apoenzyme (P32 or P21), whose
crystals diffracted to ~2 Å, several data sets were collected on
crystals soaked with ligands or co-crystallized complexes. The data
were collected to ~2.2 Å resolution for the ligand soaks.
Structural Solution and Refinement--
The phase problem was
solved by using the molecular replacement method with the MJ0109 model
(PDB code 1DK4) as a probe structure. The rotational as well
translational searches were performed using the program Amore (18).
The searches succeeded with the monomer as a model. The initial
solution was refined by rigid body refinement and simulated annealing
(crystallography NMR software) with manual rebuilding using
O (19). The apoenzyme in the crystal form characterized by the
space group P32, initially indexed as P3221,
was evaluated for merohedral twinning (20) and found to be nearly
perfectly twinned with a twinning ratio of 0.499. The reflections for
this trigonal form of the apoenzyme were then sorted and processed in
space group P32. Simulated annealing refinement in the
program Xplor (21) lowered the R-factor
significantly. Several macrocycles including manual rebuilding sessions
done on the SGI work station with the program Xtalview (22) led
to convergence. The final models were refined by the program Shelxl-97 (23). The models were evaluated using Procheck (24). Table II provides a summary of data collection
and refinement parameters for both the apoenzyme crystal forms and the
ligand-soaked ternary complexes.
View this table:
[in this window]
[in a new window]
|
Table II
Statistics for data collection and refinement of apoenzyme and ternary
complexes of AF2372 with metal ions and substrates
|
|
 |
RESULTS |
AF2372 Affinities for Substrate and Metal Ions--
At the
preparatory stage for crystallographic experiments, the substrate
specificity as well as the metal ion affinity of AF2372 were evaluated
kinetically (at 85 °C). The substrate concentration, 2 mM, was chosen to be much greater than the
Km for Ins-1-P or FBP (2), so that a comparison of
activities is likely to reflect differences in
Vmax. The substrate specificity of AF2372 is
similar to that of the M. jannaschii ortholog (Fig.
1A) with Ins-1-P, FBP,
-glycerol phosphate, 2'-AMP, and pNPP as substrates. Except for FBP,
no other -CH2OPO
phosphate of the sugar phosphates examined was removed by the enzyme.
The Km values for L-Ins-1-P and FBP have
previously been determined as 0.11 and 0.08 mM,
respectively (2). With 2 mM L-Ins-1-P or FBP
(values well above Km), the Kd values for Mg2+ and Mn2+ with these
substrates were found to be 15-30 mM (Table
III). There was a 2-fold difference in
Kd values for Mg2+ using
L-Ins-1-P or FBP as substrate. It is possible that the presence of two phosphate groups in FBP facilitate metal ion
interactions within the active site. Similar to the phosphatase from
M. jannaschii (25), Ca2+ was an
inhibitory metal ion for both Ins-1-P and FBP hydrolysis catalyzed by
AF2372.

View larger version (13K):
[in this window]
[in a new window]
|
Fig. 1.
Substrate specificity and Li+
inhibition of AF2372. In A, the activity of the enzyme
is monitored against 2.0 mM of each phosphomonoester
(except p-nitrophenyl phosphate (pNPP), which was 0.5 mM) in 20 mM Mg2+, 50 mM Tris, pH 8.0, at 85 °C. In B, the effect
of different Li+ concentrations on enzyme activity toward
2.0 mM L-Ins-1-P is shown.
|
|
The inhibition of AF2372 phosphatase activity (with
L-Ins-1-P as substrate) by Li+ was
very weak with an IC50 ~290 mM (Table I).
This very weak binding precludes an accurate determination of the
mechanism of inhibition as uncompetitive (as many other IMPases) or
noncompetitive. Regardless of which mode of inhibition is correct, the
Ki would be ~290 mM because the
substrate concentration was well above Vmax. For
comparison, 300 mM Na+ or K+ had no
effect on AF2372 activity indicating that the inhibition by
Li+ is specific to that monovalent cation. The AF2372 gene
product, like MJ0109 (25) and a similar activity from
Thermotoga maritima (26) is insensitive to
Li+. For comparison, Li+-sensitive
representatives in the IMPase superfamily like human IMPase (27) have
much lower ki values (0.3 mM).
Similarly, the Li+ ki for Hal2 was 0.07 mM (11).
General Description and Quality of Crystal Structures of AF2372, a
Comparison to MJ0109--
The A. fulgidus IMPase/FBPase
crystallized in two distinct crystal forms: P21 and
P32. Although these two crystal forms of apoAF2372 had
slight differences in orientation of one monomer with respect to the
other, the r.m.s.d. calculated for C
positions of the dimers was
~0.6 Å indicating that the two forms are very similar. The crystal
packing of the two forms in the crystal lattice was different. This
finding is consistent with the P32 form having an unusually
high solvent content (69% as measured by Ficoll gradient). In fact,
this crystal form of the apoenzyme easily dissolved under slightly
changed conditions in the well. The major crystal contacts in both
crystal forms involved a mobile loop comprised of residues 32-43. In
P21 and P32 lattices, the loop had contacts
with the same region at the back of the symmetry-related molecule
around residue 110. An additional contact was formed in the
P21 lattice to the two N-terminal helices. This contact
breaks the symmetrical contact of residues 37-110 that is most likely
responsible for crystal twinning in the P32 lattice.
The overall fold of AF2372 is highly conserved in the entire IMPase
superfamily from archaea to eukaryotes. IMPases from eukaryotes are
organized as homodimers, while FBPases are tetrameric (a dimer of
dimers). The AF2372 structure was organized in a similar fashion to
that of dimeric MJ0109 with the same domain and global secondary structure organization, although there were some changes at the dimer
interface. The dimer interface of AF2372 was large, ~1200 Å2. Hydrophobic contacts constituted 47% of the total
interface, polar contacts 14.2%, and charged contacts 33.2% of the
interface. There appeared to be 11 hydrogen bonds across the interface.
Despite a very extensive nonpolar component coming from a cluster of
hydrophobic residues (six phenylalanines), the charged and polar
contacts appeared to be key stabilizing interactions for the dimer.
Consistent with the structural observations, soaking crystals in high
Ca2+ (>100 mM) disrupted the dimer interface.
More specific interactions in the P21 apoenzyme included
hydrogen bonding at the interface as well as a possible salt bridge
between positions Arg-27 and Glu-138.
The common architecture of each monomer consisted of an
/
structure organized in two domains of alternating layers of
-helices and
-sheets connected by variable length loops (Fig.
2). The structures of proteins in the
IMPase superfamily differ mostly by the size of the connecting loops.
Both domains in AF2372 were connected by a large hinge-like loop that
divided the monomer into two structural fragments. Each domain was
organized in an alternating 



motif with one possessing a
seven-stranded sheet and the other a five-stranded sheet. The larger
sheet was flanked by two helices on both sides, while the smaller one
had only a single helix in the top layer (Fig. 2A). The
elongated loop that separated the monomer into two domains may provide
enough flexibility to cause conformational changes within the monomer
that affect ligand binding at the active site. That type of mechanism
is used in inhibiting mammalian FBPases (28).

View larger version (59K):
[in this window]
[in a new window]
|
Fig. 2.
A, ribbon representation of
AF-IMPase dimer. B, superposition of AF2372 (in
thick red line) and MJ0109 (in thin blue lines)
dimers. Substrates are in black and purple.
|
|
Overall, the shape, volume, and even to the certain degree surface
accessible area were conserved for both proteins. However, despite
having the same number of amino acids, there were significant conformational differences between AF2372 and MJ0109 (Fig.
2B). The N-terminal helix in the A. fulgidus
IMPase was slightly elongated relative to that in MJ0109. The second
N-terminal helix was also slightly elongated compared with MJ0109 (the
MJ0109 helix is from residues 45-54; in AF2372 it comprises residues
44-58). The shorter loop regions in MJ0109 might be one of the major
factors that confer greater stability for the M. jannaschii
homologue (Tm = 87 °C for AF2372 versus
95 °C for MJ0109 (15)). The remaining helices of AF2372 were very
similar to both MJ0109 and human IMPase in length and positioning
(6) with one exception. The first significant difference was at the
elongated loop 105-112 in AF2372 that provides most of the crystal
contacts in both lattices. There were some small differences at the
loop around residue 148, but an important difference was at the first
external helix. In MJ0109, this helix contributed to the dimer
interface. A significant deletion in the sequence eliminated the
interfacial helix present in MJ0109. The place of the helix was taken
by a loosely packed coil. This structural element caused a much looser
organization of the dimer that might also be associated with decreased
thermal stability of AF2372 as compared with MJ0109.
Another interesting structural variation in the AF2372 structure was
the presence of the cis-peptide bond between Tyr-155 and
Tyr-156. The cis organization of this peptide bond allowed formation of
an additional hydrogen bond between different strands of the second
-sheet. Both tyrosine residues seemed to play important roles in
dimer formation; they also contributed ligands for binding of both
substrates. When FBP was bound at the active site, it was coordinated
by the hydroxyl group of Tyr-156, while when L-Ins-1-P was
bound it was coordinated by the hydroxyl group of Tyr-155. Non prolyl
cis-peptide bonds occur in proteins in relatively low amounts (~0.15%). Their presence in an unfolded state of the protein can make a significant contribution to folding kinetics. A barrier for
cis to trans transition can be the rate-limiting step for some slow
folding molecules (29). However, it is relatively rare to observe
cis-peptide bonds in the folded state. The
cis-peptide bond between Tyr-155 and Tyr-156 may be an
additional factor in the diminished thermal stability of the protein
and in the irreversibility of the thermal denaturation (after heating
at ~95 °C AF2372 cannot refold into an active protein (15).
Structural Details of Two Apo Forms and Complexes with
Ins-1-P/2 Ca2+, FBP/3
Ca2+, and F-6-p + Pi/3
Mn2+--
There were two structural forms of the
apoenzyme. The first one was the P21 crystal form. For this
form, the data were collected, and the structure determined to 1.8 Å resolution. The electron density at the first subunit was of
exceptional quality (Fig. 3). The second
subunit seemed to have more mobility, and the resulting map was of
slightly lesser quality. Both subunits appeared to have an active site
devoid of metal ligands. The second crystal form, refined at slightly
lower resolution (~2 Å), exhibited an apparent symmetry
P3121. After a closer examination of the resulting distribution of intensities, we determined that the crystal form was an
ideal merohedrally twinned structure in the P32 space
group. This crystal form also showed no significant occupancy of metal ions at the active site (Fig.
4A). However, the packing
relations were sufficiently different in both crystal forms that it
provided an independent estimate of conformational variability of the
catalytic loop. The ordering of the catalytic loop is much higher in
the P32 crystal form; that is best expressed by the ratio
of the average temperature factors (B's) for the main chain atoms of
the mobile loop divided by the average B's for the entire
structure. This ratio drops from around 2.8 in the P21
crystal form to ~1.2 in P32.

View larger version (37K):
[in this window]
[in a new window]
|
Fig. 3.
The C-terminal residues hydrogen bonded to
the extensive water network as refined in P21 crystal form
at 1.8 Å covered by the 2 Fo Fc electron
density map contoured at 1.6 level.
|
|

View larger version (40K):
[in this window]
[in a new window]
|
Fig. 4.
The active site of AF2372 covered with 2 Fo Fc electron
densities contoured at 1.5 level with Asp-82,
Asp-85, and Asp200 as the reference points. A, an
apoenzyme empty active site in P32 crystal form.
B, P21 crystal form with D-Ins-1-P
and two Ca2+ ions; C, P32 crystal
form with FBP and three Ca2+ ions. D,
P32 crystal form with F-6-P, inorganic phosphate, and three
Mn2+ ions.
|
|
The third structure was for the enzyme with D-Ins-1-P and
two Ca2+ bound at the active site (Fig. 4B).
This structure was derived by soaking the P21 crystals with
20 mM substrate and 50 mM Ca2+. The
occupancy of metal ions as well as the substrate molecule refined
around 60%. The precise description of the active site and
coordination of the metal ions was hindered by the resolution and the
quality of the data. However, general features were easily discernible,
and they are described below. There were two Ca2+ metal
ions at each active site visible in the structure. Both calcium ions
were coordinated by five ligands arranged in a distorted octahedral
coordination with one ligand missing. The missing ligand could be
easily complemented by a water molecule. The ligands of Ca1 included
Asp-82, Asp-85, Asp-200, and two oxygens of the substrate molecule,
while Ca2 was coordinated by Asp-82, Glu-67 OE1, and a single oxygen of
the phosphate group. The inositol moiety was coordinated by the Asp-85,
backbone nitrogen of Ala-172, guanidino group of Arg-191, and Tyr-155.
The interactions of metal ions and substrates with the enzyme
are summarized in Table IV.
The fourth structure was derived from the P32 crystals
soaked with FBP and Ca2+ at concentrations of 20 mM and 50 mM, respectively. The refined occupancy for the metals as well as the substrate was 100%. Three Ca2+ ions were present in this structure (Fig.
4C). The coordination of two calcium ions was analogous to
that observed in the D-Ins-1-P-enzyme complex. The third
Ca3 metal ion was located more remotely from the active site than the
third metal ion in MJ0109. In subunit A it was bound at a classical
hepta-coordinated site formed by five water molecules and two protein
ligands originating from the mobile loop (Thr-40 O
1 and carbonyl
oxygen of Pro-41). At subunit B two of the water molecules were
missing. Ca3 was quite removed from the vicinity of other metal ions
with distances ~4.5 Å from Glu-67 and ~4.5 Å from Asp-38. Both of
these tentative ligands have a bridging water molecule coordinating
Ca3. This implies that if it participates in catalysis, additional
motion of the catalytic loop would be required to bring it to the
immediate proximity of the tentative nucleophile. Additionally, the
structure of MJ0109 suggested that activating and not inhibitory metal
ions were necessary to observe a catalytically competent third metal ion binding site.
Therefore, we have obtained a crystal structure of the product complex
with activating metal ions (3 Mn2+, fructose-6-phosphate
(F-6-P), and inorganic phosphate) in P32 crystal form. The
occupancy of the metal ions and the phosphate moiety refined to 100%,
while the product F-6-P refined at ~60%. As expected, the loop moved
toward the active site, and the direct coordination of the third metal
ion changed (Figs. 4D and 5D). The Mn3 has an
octahedral coordination with three protein ligands (side chains of
Asp-38, Thr-40, and Glu-67) as well as three water molecules directly
coordinating the metal ion. This is in marked difference with the
FBP-calcium complex in which the third metal ion was lacking a
direct protein bridge (Glu-67) to the second metal ion. Additionally,
one of the inorganic phosphate oxygens was directly coordinated to
Thr-84, a feature that clearly suggested the collapsed product complex.
This oxygen position was previously suggested to be a site for the
nucleophilic water molecule in a precatalytic complex. The detailed
description of the ternary complexes will be published separately.
Structural Differences among Two Apo Forms and Complexes with
Substrates and Metal ions, a Conformation of the Catalytic
Loop--
The first loop that directly follows the N-terminal helix
(residues 32-43) is a critical catalytic loop. This mobile loop was
relatively disordered in the apoenzyme but became more ordered upon
ligand binding. Thermal factors of the refined apoenzyme structures
consistently showed elevated values in this region; this is an
indicator of mobility and/or disorder. The structures of AF2372 in
complexes with Ins-1-P/Ca2+ or FBP indicated that, in
analogy to the behavior of FBPase and MJ0109, the catalytic loop is
stabilized by the presence of metal ions and ligands.
The refined conformation of the mobile loop (residues 32-43) was
slightly different in the apo form and when a ligand was bound (Fig.
5, A, B,
and C). Note the difference in the quality of the
electron density comparing the apoenzyme to the liganded structures.
The changes in B factors for this region when crystals were soaked in
ligand solutions did fully reflect the increased "ordering" of the
loop. There was a systematic drop in the B factor for the protein plus
ligand structures. The ratio of the average B factor for residues
32-43 to the rest of the molecule approached 1 for both of the
substrate/product complexes with three metal ions (Ca2+ or
Mn2+).

View larger version (36K):
[in this window]
[in a new window]
|
Fig. 5.
The refined atomic models of the mobile
catalytic loop covered by the 2 Fo Fc electron
density at contouring levels. A, 1 in apo-P21.
B, 1.5 in IMP + 2 Ca2+. C, 1.5 FBP + 3 Ca2+. D, 1.5 in F6P + 3 Mn2+
complexes.
|
|
Soaking in metal ions alone did not improve loop stability
significantly. In fact, crystals soaked or co-crystallized with metal
ions alone diffracted poorly (>3 Å). Under conditions that gave the
P32 form of the enzyme, the combination of both substrate and metal ions was needed to stabilize the loop in these Tris/PEG 3350 conditions. When the AF2372 apoenzyme was soaked in
D-Ins-1-P, as well as other substrate/metal complexes, a
conformational change occurred that shifted the loop on average ~2.5
Å closer to the metal ions following the behavior observed in MJ0109
structures that contained ligands at the active.
Architectural Differences with Other Members of the Family, a
Catalytic Loop Variability--
The length of the amino acid peptide
chain representing the members of the extended IMPase family
differs substantially and extends from 252 residues for MJ109 and
AF2372 to 400 for IPPase. Structural alignments carried out with
SEQUOIA identified the common catalytic core despite relatively low
levels of primary sequence homology. However, as inferred from the size
difference, extensive external loops and even small domains were formed
as additions to the core domain. The superposition in SEQUOIA indicated that AF2372 and MJ109 were the closest structural homologs and could be
superimposed with 209 positions to an r.m.s.d. of 1.5 Å (Table
V). There were fewer structurally
conserved residues and higher r.m.s.d. when AF2372 was superimposed
onto the other IMPase or FBPase structures (Table V). However, two
clusters of similarities emerged. The structures in the individual
clusters had closer similarities to each other than to proteins in the other cluster. AF2372, MJ0109, and FBPase formed one cluster, while
Hal2 and human IMPase and IPPase formed the second cluster. The extent
of the r.m.s.d. was up to 2.5 Å as measured by superposition of AF2372
with human IPPase.
View this table:
[in this window]
[in a new window]
|
Table V
The results of superposition for members of the IMPase superfamily
performed with program Sequoiaa.
Shaded cells indicate Li+-insensitive and Li+-sensitive
subfamilies.
|
|
Despite significant structural additions, the structures superimposed
well around the catalytic loop. The two N-terminal helices provided a
good anchor for judging the conformational variability of this loop.
The five structures of AF2372 determined in different crystal lattices
and with different numbers of metal ions bound provided an additional
measure of structural divergence. The two apo forms had a more open
conformation of the loop, while the metal and substrate bound forms had
the loop in a closed conformation. The overall extent of this motion
was similar to that observed for MJ0109 and amounts to a displacement
of the tip of the loop by ~2.5 Å (Fig.
6A).

View larger version (48K):
[in this window]
[in a new window]
|
Fig. 6.
Analysis of mobile loop conformation
in the IMPase superfamily. A, superposition of all five
structures of AF2372 with P21 apoenzyme (in
blue), P32 apoenzyme (in yellow), IMP + 2 Ca2+ (in purple), FBP + 3 Ca2+
(in green), and F-6-p + Pi+3 Mn2+ (in
red). B, superposition of
Li+-insensitive enzymes AF2372 (in red) with
MJ0109 (in green) and pig kidney FBPase (in
orange). C, superposition of
Li+-sensitive enzymes human IMPase (in purple),
human IPPase (in light blue) with Hal2 (in dark
blue). The mobile catalytic loop is indicated by the black
arrow.
|
|
Separate from changes in conformation of the catalytic loop in the
presence or absence of ligand, there were significant differences in
the loop conformation in the IMPase superfamily. Superimposing all six
known members of the family elucidated unexpected features. While the
common core superimposed quite well, the fragments containing the
mobile loop formed two structural clusters. The loops of AF2372, MJ0109, and pig kidney FBPase were much closer to each other but further away from the loops of human IMPase, IPPase, and Hal2 that
seemed to belong to a separate cluster (Fig 6, B and
C).
 |
DISCUSSION |
There have been a number of studies directed at elucidating the
details of Li+ effects on the IMPase superfamily. A series
of combined mutational, spectroscopic, and kinetic experiments carried
out on the human IMPase implicated two residues, Cys-218 and His-217,
in Li+ inhibition (8). Both residues are remote from the
metal ion binding site and, if confirmed as part of the Li+
binding site, they cannot be involved in catalysis. There was speculation that these residues could influence the binding of the
metal ion at site 2, which would hinder the release of phosphate and
result in inhibition of IMPase activity (8). This proposed mechanism
was termed the "bumping hypothesis." Two similar mutations were
made in one of the Li+-insensitive archaeal IMPase/FBPase
enzymes (MJ0109). Mutation of Arg-198 in the M. jannaschii
enzyme (aligned with His-217 in the mammalian enzyme) to His led to
complete loss of IMPase activity (25). However, introduction of
cysteine in place of Ala-199 (which aligned with Cys-218 of the human
IMPase) increased the sensitivity of MJ0109 to Li+
(IC50 decreased from ~150 to ~ 80 mM
(25)). Such a relatively small change in Li+ sensitivity
might be interpreted as stemming from other effects. It is more likely
that elements conferring Li+ sensitivity involve secondary
structural elements that require more than a point mutation.
Random mutagenesis studies carried out on the product of the yeast Hal2
gene provided other hypotheses. The Hal2 gene is an in vivo target of
Na+ and Li+ toxicity, and its overexpression
improves the tolerance of higher plants to salt. This enzyme is
differentially sensitive to monovalent cations: tolerant of high sodium
concentrations but inhibited by lithium. Mutating residues Val-170 and
Trp-293 or one of the residues forming a salt bridge between Glu-238
and Arg-152 caused a marked decrease in Li+ sensitivity
(11). It was speculated that the salt bridge between Arg-152 to Asp-263
might be stabilized by an ionic network involving monovalent ion
binding. The hydrophobic interaction between Val-170 and Trp-292 was
required for an effective formation of the catalytic core that, when
disrupted, changed sensitivity to Na+ and Li+
more than catalytic activity. Mutations at these positions shifted the
IC50 for wild type Hal2 from 0.07 to 4 mM in
E238K mutants and to 0.5 mM in the V70A mutant.
In the case of human IMPase, a two-metal ion-assisted catalytic
mechanism was proposed. Li+ binding was proposed to
interfere with the second metal ion binding site (4), in effect,
preventing product release. In the case of archaeal IMPases, we have
previously proposed that the catalysis should be described as a
three-metal-assisted catalysis (6). Additionally, this proposal led to
the hypothesis that the mobile loop was responsible for binding of the
third metal ion. According to this hypothesis, the catalytic cycle is
initiated when the second and the third metal ions are transiently
localized in close proximity to create the water nucleophile. The
nucleophile immediately reacts with an incoming phosphate group to
effectively hydrolyze it with a reversal of phosphorus configuration.
Subsequently, the much more mobile third metal ion facilitates release
of the phosphate.
The direct consequence of this mechanism is the following hypothesis
for Li+ inhibition. The mobile loop is important for
localization and binding of the third metal ion, and this metal ion can
be replaced by Li+. However, Li+ at this site
would not lead to the formation of the water-based nucleophile, and
enzyme activity would be inhibited. This scenario offers a direct
insight into why Li+ but not other more abundant monovalent
metal ions like Na+ or K+ are inhibitors. To
promote the nucleophile formation, the third metal ion needs to insert
deep into the active site. The Mg2+ ionic radius is 0.66 Å. Therefore, only Li+ is small enough (0.68 Å) to
substitute structurally for the Mg2+ ion. However, because
of its charge it would be unable to create a nucleophile, thus
inhibiting the enzyme. The much larger Na+ or
K+ ions (0.97 Å and 1.33 Å, respectively) would be able
to activate the enzyme by general electrostatic rearrangement but would
be unable to substitute structurally for Mg2+. As is
suggested by the difference in ionic radii they must have much lower
affinity for the enzyme site. Only at exceedingly high concentrations
would they cause significant inhibition. The existence of
conformational classes for the mobile loop as found in this paper can
be directly correlated with Li+ inhibition. The
representative enzymes with a smaller and more compact mobile loop,
like AF2372, respond to Li+ in a much more subdued manner
than enzymes with longer and more flexible catalytic loop. Therefore,
the IMPase superfamily can be subdivided into two distinct classes as
defined by a key structural feature (mobile loop conformation) and its
functional importance (Li+ inhibitory properties).
The use of "mobile loops" as sensors of monovalent ions is quite
common in proteins. There are many examples of ubiquitous sodium/potassium transporters and ion channels that utilize mobile loops and differential cation sensitivity to operate. Na+-
and Cl
-dependent transporters GLYT1 and GLYT2
remove glycine from the synaptic cleft. Extracellular loop 1 in the
GLYT2 transporter was identified as the source of "structural
heterogeneity" that is involved in the specific effect of lithium on
serotonin transport (30). Similarly, the glutamate transporter GLT-1
has a region identified as loop-like that is sodium-sensitive and
coupled to drive glutamate transport. When Ser-440 in this loop was
mutated to glycine, the GLT-1 mutate could use Li+ to drive
glutamate transport (31). The Na+/Ca2+
exchangers NCX1 and NCX3 reveal structural domains (i.e.
loops) important for differential sensitivity to external
Ni2+ or Li+ (32). Although the last of these is
a totally separate biological system, this transporter has ion
sensitivities that parallel those of Li+-sensitive IMPases.
In the extended IMPase family the mobile loop appears to have a double
role. Its mobility is associated with and controlled by the metal
ion/substrate binding. The structures of archaeal dual function
IMPase/FBPases suggest that this loop mobility is required for the
efficient release of product (i.e. phosphate). The second
role of the mobile loop is to differentiate the influence exerted by
monovalent metal ions. In Li+-sensitive enzymes, the loop
appears to be longer and in a more open (farther away from the active
site) conformation, while in Li+-insensitive enzymes, the
mobile loop is shorter and in a more closed conformation. This
structural differentiation appears to be useful in classifying
IMPase-like proteins. It remains to be seen what is the precise role of
the loop in Li+ inhibition and how universal is this
structural classification. The predictive power of this hypothesis
needs to be tested on more representatives of this family of proteins.
Recently, a new human IMPase (35, 36) as well as new human BPNTase (37)
were cloned. The new inositol monophosphatase termed A2 appears to have
significant sequence homology to already characterized
Li+-sensitive IMPase. The amino acid sequence indicates a
larger protein with N- and C-terminal extensions. However, the
significant sequence homology in terms of the length and character is
observed around the catalytic mobile loop. That homology would suggest that this enzyme also should operate by the same three-metal
ion-assisted catalysis, and it also should be
Li+-sensitive. The kinetic data (34) suggest that the same
should hold true for IPPase and Hal2 (37). Additionally, during the submission process of this paper another enzyme from this extended protein family was characterized by x-ray crystallography (38). The rat
3'-phophoadenosione-5'-phosphate phosphatase (BPNTase) with significant
homology to the corresponding human enzyme turned out to posses not
only double function but also the ability to catalyze the hydpolysis of
1 and 4 inositol polyphosphate species. This protein is also
Li+-sensitive. The structural details appear to support
fully our findings concerning the catalytic mechanism as well as
inhibition by Li+ described previously (6) as well as
observations reported in this paper.
 |
ACKNOWLEDGEMENTS |
K. A. S. thanks Angela Criswell for helpful
discussions during her visits at Rice University and Edgar Ortiz,
supervisor of the x-ray facility in the Department of Biophysics and
Physiology, Boston University School of Medicine, for help in use of
the facility.
 |
FOOTNOTES |
*
This work was supported in part by Grant 2T15LM07093 from
the W. M. Keck Center for Computational Biology (to B. S.) and by Grant MCB-9978250 from the NSF, National Institutes of Health (to M. F. R.).The costs of publication of this
article were defrayed in part by the
payment of page charges. The article
must therefore be hereby marked
"advertisement" in
accordance with 18 U.S.C. Section
1734 solely to indicate this fact.
The atomic coordinates and the structure factors (code 1LBV, 1LBW, 1LBX, 1LBY, 1LBZ) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).
§
Supported with a Hitchings-Elian fellowship from Burroughs Wellcome fund.
Present address: Dept. of Biochemical Sciences, University of
Rome "La Sapienza," I-00185 Rome, Italy.

To whom correspondence should be addressed. Tel.: 713-348-3346;
Fax: 713-348-5154; Email: stec@bioc.rice.edu.
Published, JBC Papers in Press, April 8, 2002, DOI 10.1074/jbc.M201042200
 |
ABBREVIATIONS |
The abbreviations used are:
IMPase, inositol monophosphatase;
FBP, fructose 1,6-bisphosphate;
FBPase, fructose-1,6-bisphosphatase;
IPPase, inositol polyphosphate
phosphatase;
IP3, inositol 1,4,5-trisphosphate;
Ins-1-P, inositol 1-phosphate;
PEG, polyethylene glycol;
r.m.s.d., root mean
square deviation(s);
Hal2, yeast 3'-phosphodenosine-5'-phosphate
phosphatase Hal2;
BPNTase, bisphosphate nucleotidase;
F-6-P, fructose
6-phosphate.
 |
REFERENCES |
| 1.
|
Martin, D. D.,
Ciulla, R. A.,
and Roberts, M. F.
(1999)
Appl. Environ. Microbiol.
65,
1815-1825[Free Full Text]
|
| 2.
|
Stec, B.,
Yang, H.,
Johnson, K. A.,
Chen, L.,
and Roberts, M. F.
(2000)
Nat. Struct. Biol.
7,
1046-1050[CrossRef][Medline]
[Order article via Infotrieve]
|
| 3.
|
York, J. D.,
Ponder, J. W.,
and Majerus, P. W.
(1995)
Proc. Natl. Acad. Sci. U. S. A.
92,
5149-5153[Abstract/Free Full Text]
|
| 4.
|
Bone, R.,
Frank, L.,
Springer, J. P.,
and Atack, J. R.
(1994)
Biochemistry
33,
9468-9476[CrossRef][Medline]
[Order article via Infotrieve]
|
| 5.
|
Choe, Y. Z.,
Poland, B. W.,
Fromm, H. J.,
and Honzatko, R. B.
(1998)
Biochemistry
33,
11441-11450
|
| 6.
|
Johnson, K. A.,
Chen, L.,
Yang, H.,
Roberts, M. F.,
and Stec, B.
(2001)
Biochemistry
40,
618-630[CrossRef][Medline]
[Order article via Infotrieve]
|
| 7.
|
Vadnal, R.,
and Parthasarathy, R.
(1995)
Neuropsychopharm.
12,
277-285
|
| 8.
|
Atack, J. R.,
Groughton, H. B.,
and Pollack, S. J.
(1995)
FEBS Lett.
361,
1-7[CrossRef][Medline]
[Order article via Infotrieve]
|
| 9.
|
Atack, J. R.
(1996)
Brain Res. Rev.
22,
183-190[CrossRef][Medline]
[Order article via Infotrieve]
|
| 10.
|
York, J. D.,
Ponder, J. W.,
Chen, Z. W.,
Matthews, F. S.,
and Majerus, P. W.
(1994)
Biochemistry
33,
13164-13171[CrossRef][Medline]
[Order article via Infotrieve]
|
| 11.
|
Albert, A.,
Yenush, L.,
Gil-Mascarell, M. R.,
Rodriguez, P. L.,
Patel, S.,
Martinez-Ripoll, M.,
Blundell, T. L.,
and Serrano, R.
(2000)
J. Mol. Biol.
295,
927-938[CrossRef][Medline]
[Order article via Infotrieve]
|
| 12.
|
Zhang, R.,
Villeret, V.,
Lipscomb, W. N.,
and Fromm, H. J.
(1996)
Biochemistry
35,
3038-3043[CrossRef][Medline]
[Order article via Infotrieve]
|
| 13.
|
Villeret, V.,
Huang, S.,
Fromm, H. J.,
and Lipscomb, W. N.
(1995)
Proc. Natl. Acad. Sci.
92,
8916-8920[Abstract/Free Full Text]
|
| 14.
|
Chen, L.,
Zhou, C.,
Yang, H.,
and Roberts, M. F.
(2000)
Biochemistry
39,
12415-12423[CrossRef][Medline]
[Order article via Infotrieve]
|
| 15.
| Yang, H., Mohanty, U., Stieglitz, K., and Roberts, M. F. (2002)
submitted for publication.
|
| 16.
|
Itaya, K.,
and Michio, U.
(1966)
Clin. Chim. Acta
14,
361-366[CrossRef][Medline]
[Order article via Infotrieve]
|
| 17.
|
Otwinowski, Z.
(1993)
in
Oscillation Data Reduction Program: Proceedings of the CCP4 Study Weekend Data Collection and Processing
(Sawyer, L.
, Isaacs, N.
, and Bailey, S., eds)
, pp. 56-62, SERC Daresbury Laboratory, Warrington, UK
|
| 18.
|
Navaza, J.
(1994)
Acta Crystallogr. Sect. A
50,
157-163[CrossRef]
|
| 19.
|
Jones, T. A.,
and Kjeldgaard, M.
(1994)
in
From First Map to Final Model
(Bailey, S.
, Hubbard, R.
, and Waller, D., eds)
, pp. 1-13, SERC Daresbury Laboratory, Warrington, UK
|
| 20.
|
Yeates, T. O.
(1997)
Methods Enzymol.
276,
344-358[Medline]
[Order article via Infotrieve]
|
| 21.
|
Brunger, A. T.
(1992)
X-PLOR, Version 3.1
, Yale University Press, New Haven, CT
|
| 22.
|
McRee, D. E.
(1999)
J. Struct. Biol.
125,
156-165[CrossRef][Medline]
[Order article via Infotrieve]
|
| 23.
|
Sheldrick, G. M.,
and Schneider, T. R.
(1997)
Methods Enzymol.
277,
319-343
|
| 24.
|
Laskowski, R. A.,
MacArthur, M. W.,
Moss, D. S.,
and Thornton, J. M.
(1993)
J. Appl. Cryst.
26,
283-291
|
| 25.
|
Chen, L.,
and Roberts, M. F.
(1998)
Appl. Environ. Microbiol.
64,
2609-2615[Abstract/Free Full Text]
|
| 26.
|
Chen, L.,
and Roberts, M. F.
(1999)
Appl. Environ. Microbiol.
65,
4559-4567[Abstract/Free Full Text]
|
| 27.
|
McAllister, G.,
Whiting, P.,
Hammond, E. A.,
Knowles, M. R.,
Atack, J. R.,
Bailey, F. J.,
Maigetter, R.,
and Ragan, C. I.
(1992)
Biochem. J.
284,
749-754
|
| 28.
|
Choe, J. Y.,
Fromm, H. J.,
and Honzatko, R. B.
(2000)
Biochemistry
39,
8565-8574[CrossRef][Medline]
[Order article via Infotrieve]
|
| 29.
|
Pappenberger, G.,
Aygun, H.,
Engels, J. W.,
Reimer, U.,
Fischer, G.,
and Kiefhaber, T.
(2001)
Nat. Struct. Biol.
8,
452-458[CrossRef][Medline]
[Order article via Infotrieve]
|
| 30.
|
Lopez-Corcuera, B.,
Nunez, E.,
Martinz-Maza, R.,
Geelings, A.,
and Aragon, C.
(2001)
J. Biol. Chem.
276,
43463-43470[Abstract/Free Full Text]
|
| 31.
|
Zhang, Y.,
and Kanner, B. I.
(1999)
Proc. Natl. Acad. Sci. U. S. A.
96,
1710-1715[Abstract/Free Full Text]
|
| 32.
|
Iwamoto, T.,
Uehara, A.,
Nakamura, T. Y.,
Imanaga, I.,
and Shigekawa, M.
(1999)
J. Biol. Chem.
274,
23094-23102[Abstract/Free Full Text]
|
| 33.
|
Inhorn, R. C.,
and Majerus, P. W.
(1988)
J. Biol. Chem.
263,
14559-14565[Abstract/Free Full Text]
|
| 34.
|
Lopez-Coronado, J. M.,
Belles, J. M.,
Lesage, F.,
Serrano, R.,
and Rodriguez, P. L.
(1999)
J. Biol. Chem.
274,
16034-16039[Abstract/Free Full Text]
|
| 35.
|
Sjoholt, G.,
Gulbrandsen, A. K.,
Lovlie, R.,
Berle, J. O.,
Molven, A.,
and Steen, V. M.
(2000)
Mol. Psychiatry
5,
172-180[CrossRef][Medline]
[Order article via Infotrieve]
|
| 36.
|
Yoshikawa, T.,
Padigaru, M.,
Karkera, J. D.,
Sharma, M.,
Berrettini, W. H.,
Esterling, L. E.,
and Detera-Wadleigh, S. D.
(2000)
Mol. Psychiatry
5,
165-171[CrossRef][Medline]
[Order article via Infotrieve]
|
| 37.
|
Spiegelberg, B. D.,
Xiong, J. P.,
Smith, J. J., Gu, R. F.,
and York, J. D.
(1999)
J. Biol. Chem.
274,
13619-13628[Abstract/Free Full Text]
|
| 38.
|
Patel, S.,
Yenush, L.,
Rodriguez, P. L.,
Serrano, R.,
and Blundell, T. L.
(2002)
J. Mol. Biol.
315,
677-685[CrossRef][Medline]
[Order article via Infotrieve]
|
Copyright © 2002 by The American Society for Biochemistry and Molecular Biology, Inc.

CiteULike
Complore
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:

|
 |

|
 |
 
Y. Wang, K. A. Stieglitz, M. Bubunenko, D. L. Court, B. Stec, and M. F. Roberts
The Structure of the R184A Mutant of the Inositol Monophosphatase Encoded by suhB and Implications for Its Functional Interactions in Escherichia coli
J. Biol. Chem.,
September 14, 2007;
282(37):
26989 - 26996.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. Fukuda, S. Kawai, and K. Murata
NADP(H) Phosphatase Activities of Archaeal Inositol Monophosphatase and Eubacterial 3'-Phosphoadenosine 5'-Phosphate Phosphatase
Appl. Envir. Microbiol.,
September 1, 2007;
73(17):
5447 - 5452.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. A. Rodionov, O. V. Kurnasov, B. Stec, Y. Wang, M. F. Roberts, and A. L. Osterman
Genomic identification and in vitro reconstitution of a complete biosynthetic pathway for the osmolyte di-myo-inositol-phosphate
PNAS,
March 13, 2007;
104(11):
4279 - 4284.
[Abstract]
[Full Text]
[PDF]
|
 |
|