Originally published In Press as doi:10.1074/jbc.M204025200 on August 20, 2002
J. Biol. Chem., Vol. 277, Issue 43, 41101-41109, October 25, 2002
GTP-Yeast Actin*
Kuo-Kuang
Wen,
Xiaoyi
Yao, and
Peter A.
Rubenstein
From the Department of Biochemistry, University of Iowa College of
Medicine, Iowa City, Iowa 52242
Received for publication, April 25, 2002, and in revised form, August 7, 2002
 |
ABSTRACT |
Because of the apparently greater
conformational flexibility of yeast versus muscle actin and
the ability of other members in the actin protein superfamily to
efficiently use both ATP and GTP, we assessed the ability of yeast
actin to function with GTP. Etheno-ATP exchange studies showed that the
binding of GTP to yeast actin is about 1/9 as tight as that of ATP in
contrast to the 1/1,240 ratio for muscle actin. Proteolysis of
GTP-bound G-yeast actin suggests that the conformation of subdomain 2 is very much like that of ATP-bound actin, but CD studies show that
GTP-bound actin is less thermostable than ATP-bound actin. GTP-actin
polymerizes with an apparent critical concentration of 1.5 µM, higher than that of ATP-actin (0.3 µM)
although filament structures observed by electron microscopy
were similar. Yeast actin hydrolyzes GTP in a
polymerization-dependent manner, and GTP-bound F-actin
decorates with the myosin S1. Conversion of Phe306 in the
nucleotide binding site to the Tyr found in muscle actin raised the
nucleotide discrimination ratio from the 1/9 of wild-type actin
to 1/125. This result agrees with modeling that predicts that removal
of the Tyr hydroxyl will create a space for the C2 amino group of the
GTP guanine.
 |
INTRODUCTION |
Two cytoskeletal filament systems in the cell, actin
microfilaments and tubulin-containing microtubules, depend on the
ability to bind and hydrolyze nucleoside triphosphates for proper
function. For microtubules, the nucleotide is GTP, and for the actin
filaments, the nucleotide is assumed to be ATP (1). ATP F-actin is more stable than ADP F-actin, and the instability brought about by hydrolysis of the nucleotide is believed to be an important factor in
the regulation of the dynamics of the actin cytoskeleton (2).
The ATP requirement for actin is based on two criteria: the isolation
of actin from a number of sources with a bound ATP, and the
demonstration that there is a 500-fold higher affinity of muscle actin
for ATP than for GTP (3). This difference in nucleotide affinity was
also observed in polymerization experiments. Martonosi and Gouvea (4)
first reported that actin, in the presence of 0.1 M KCl,
would not polymerize in the presence of a stoichiometric amount of GTP.
Iyengar and Weber (3) also generated muscle G-actin free of
nucleoside triphosphates and examined the ability of different
nucleotides to affect polymerization when added back to the preparation
in the presence of salt. Their results showed that whereas 20 µM G-actin would completely polymerize in the presence of
a stoichiometric amount of ATP, 100 µM GTP was required
to polymerize the actin. Apparently, the free GTP needed to be well
above the nM-µM range required for ATP to
saturate the high affinity nucleotide binding site of the actin.
We have previously used yeast actin as a model system for examining
actin structure-function relationships in vitro, and our results have demonstrated significant quantitative differences in the
behavior of yeast versus muscle actin. Yeast actin is more susceptible to controlled proteolysis and exchanges nucleotide faster
(5). Kim et al. (6) showed that beginning with G-actin, the
yeast protein reaches steady state polymer levels more rapidly than
muscle actin. In another study, kinetic modeling demonstrated this
acceleration was because of faster nucleation and filament fragmentation during polymerization leading to more ends for filament extension (7). These results suggest that yeast actin might assume a
more open and flexible conformation than muscle actin, a hypothesis
that is consistent with additional results from ours and other
laboratories. Based on results from optical reconstruction of electron
micrographs of yeast actin filaments, Belmont et al. (8)
suggested that the area in the nucleotide cleft was more open in a
yeast than in a muscle actin filament. Furthermore, we have
demonstrated that His73 of yeast actin, unlike that of
higher actins, is not methylated (9) and that this lack of a methyl
group may allow the residues at the bottom of the nucleotide cleft to
more readily assume a conformation that favors the open conformation
(10) than is the case for other actins. If true, there might be a
decrease in nucleotide specificity of yeast in comparison with other
actins. Interestingly, Nyman et al. (11) have also shown
with chicken
-nonmuscle actin produced in yeast (non-methylated) or
isolated from tissue (methylated) that lack of methylation of this
residue results in an uncoupling of ATP hydrolysis from filament formation.
Actin is a member of a superfamily of proteins consisting of hsc70,
Escherichia coli DnaK protein, and sugar kinases,
such as hexokinase (12). All share a similar three-dimensional
structure around a nucleotide binding core, and all cleave the
-
bond of ATP resulting in a significant conformational change in the protein. However, in a number of cases, the nucleotide specificity is
rather low in comparison to actin. For example, DnaK hydrolyzes GTP with about one-half to one-fourth of the activity it has with ATP
(13), and hexokinase from rainbow trout is only about twice as active
with ATP as with GTP (14). Finally, it has recently been reported that
MreB, thought to be a bacterial version of actin, will also polymerize
in the presence of both ATP and GTP although relative efficiencies with
each of these nucleotides was not determined (15).
We have recently investigated the interaction of Neisseria porins with
yeast actin. These porins are nucleotide-gated channels that use either
ATP or GTP (16), and we wanted to determine the effect of nucleotide on
the interaction. As a control, we examined the effect of GTP on yeast
actin alone because of the characteristics of yeast actin and the
propensity of some actin superfamily members for GTP. The work reported
here shows that yeast actin, as opposed to muscle actin, will work
with GTP with relatively high efficiency.
 |
EXPERIMENTAL PROCEDURES |
Materials--
Affi-Gel 10-activated resin and Bio-Spin® 30 Tris columns were purchased from Bio-Rad. DNase I (grade D) was
purchased from Worthington Biochemicals. DNase I affinity columns were
made as described previously (17). ATP and GTP, purchased from Sigma, were dissolved in water to a final concentration of 50 mM.
The pH was adjusted to 7.5 with NaOH prior to storing at
20 °C.
The EnzChekTM phosphate assay kit and
1,N6-ethenoadenosine 5'-triphosphate
(
ATP)1 as a 5 mM stock solution in 50 mM Tris-HCl, pH 7.5, were purchased from Molecular Probes (Eugene, OR) and stored at
20 °C. The QuickChange site-directed mutagenesis kit was purchased
from Stratagene (La Jolla, CA), and the DNA primers used for
site-directed mutagenesis were obtained from Integrated DNA
Technologies (Coralville, IA). Yeast cakes for preparation of wild-type
actin were obtained from a local food store. The rabbit muscle acetone
powder and the bovine cardiac myosin S1 head preparation were generous
gifts from Dr. E. Reisler and Dr. L. Tobacman, respectively. All other
chemicals used were reagent grade quality.
Actin Purification--
Yeast wild-type actin was purified in
the Ca form by a DNase I affinity chromatography/DEAE ion exchange
chromatography protocol as described previously (17), and it was stored
in Ca G buffer (5 mM Tris-HCl, pH 7.5, 0.2 mM
ATP, 0.2 mM CaCl2, and 0.1 mM
dithiothreitol) at 4 °C and used within 3 days. Rabbit skeletal
muscle F-actin was prepared from acetone powder according to Spudich
and Watt (18). The Ca form of actin was converted to the Mg form by
adding EGTA to a final concentration equal to the calcium concentration and MgCl2 to a final concentration of 0.2 mM.
The actin was then placed on ice for at least 10 min and used within hours.
Oligodeoxynucleotide-directed Mutagenesis--
Using
site-directed mutagenesis, we replaced the codon for Phe306
in yeast actin with a Tyr codon. The actin coding sequence was contained in the centromeric plasmid pRS314 adjacent to the ACT1 promoter (19). The set of DNA primers,
5'-GGTACCACCATGTACCCAGGTATTGCC-3' and 5'-GGCAATACCTGGGTACATGGTGGTACC,
was used in which the codon for Phe306 (TTC) was
changed to Tyr (TAC). A haploid yeast strain expressing the mutant
actin as the only actin in the cell was generated as described
previously using a plasmid shuffling procedure (19). The mutant plasmid
was rescued from the cell and sequenced to ensure that the mutation was intact.
ATP- and GTP-bound Actin--
The free nucleotide was
depleted from G-actin by passing Ca-actin through a Bio-Spin® 30 Tris
column pre-equilibrated with Ca G buffer without nucleotide. Only 1/50
to 1/100 of the free nucleotide from the sample remains in the
flow-through based on results with protein-free samples (data not
shown).
ATP-bound Ca-actin was prepared by incubating 500 µM
ATP with 20 µM ATP-bound Ca-actin in
the absence of free ATP for at least 1 h on ice (overnight at
4 °C for the muscle actin). GTP-bound Ca-actin was prepared by
adding 500 µM GTP to 20 µM Ca
ATP-bound
actin from which unbound nucleotide had been removed as above and
monitoring the decrease in fluorescence until it had reached background
levels. Interestingly, the t1/2 values that
characterized the decrease in
ATP fluorescence were different for
ATP (35-40 s) and GTP (11-15 s), perhaps reflecting differences in on
rates for the two nucleotides. The GTP-actin was then depleted of
unbound nucleotide as above, and 200 µM fresh GTP was
added to the sample to yield the final GTP-actin preparation. The
ATP-bound actin and GTP-bound actin were used within 4 h.
Fluorescence Assays--
Fluorescence measurements were obtained
using a Fluorolog-3 (Jobin Yvon Inc.) fluorescence spectrometer
with a thermostatted water bath attached to the cuvette chamber. Actin
polymerization was monitored by a change in light scattering with
excitation and emission wavelengths set at 360 nm. The
ATP
fluorescence intensity change was monitored with excitation and
emission wavelengths of 340 and 410 nm, respectively.
Nucleotide-binding Titration--
Aliquots of 2.3 or 4.6 µM Ca
ATP-bound actin in the absence of free
nucleotide were added to various concentrations of ATP or GTP at
25 °C. Alternatively, ATP or GTP in small aliquots from a 50 mM stock solution were added progressively with constant stirring at 25 °C to 2.3 or 4.6 µM
ATP-actin to
achieve a final concentration of 700 µM. In both methods,
the fluorescence decrease was monitored after the nucleotide was added
until the fluorescence signal reached equilibrium. The results obtained
from both methods were identical. The observed fluorescence intensity
(Fobs) obtained by the experiments was the sum
of the fluorescence contributed by
ATP bound by actin and actin-free
ATP, and is shown in Equation 1,
|
(Eq. 1)
|
where f
A and f
are the
specific fluorescence yields derived from experiments for the bound and
free
ATP forms, respectively. [
A] is the concentration for
ATP-bound actin, [
] is that for free
ATP, and the sum of
[
A] and [
] is either 2.3 or 4.6 µM depending on
the experiment. The values of [
A] were calculated from Equation 1
and were plotted as a function of total concentration of the competing
ATP or GTP with the actin concentration constant. The solution to the
quadratic expression modified from Ref. 20 was applied to the data with
the dissociation constant ratio (KP) as a
fitting parameter using Microsoft Excel (Microsoft) and is shown
in Equation 2,
|
(Eq. 2)
|
where KP is equal to
K
/KN, KN
is the equilibrium dissociation constant for nucleotide (either ATP or
GTP) binding to actin, K
is the equilibrium
dissociation constant for
ATP, [
A] is the concentration of
ATP-bound actin, AT is the concentration of
the total
ATP, NT is the concentration of the total ATP/GTP, and ET is the concentration of
the total concentration of actin.
Thermal Denaturation--
The apparent melting temperature of
G-actin was determined by circular dichroism according to Yao et
al. (9). 1 µM G-actin in G buffer was heated at a
constant rate of 1 °C/min from 20 to 80 °C with constant
stirring. Changes in the actin ellipticity were monitored at 222 nm
using an AVIV 62 DS spectropolarimeter. Data were normalized as the
fraction of native protein based on the net change in ellipticity then
fitted to a two-state model with a single transition between a native
and a denatured form of the protein. The Tm
value was defined as the temperature at which 50% G-actin was denatured.
Limited Proteolysis of G-actin--
The digestion protocol was
essentially the same as described by Kuang and Rubenstein (21) with
modification. Briefly, 9.2 µM G-actin was incubated with
trypsin (5 µg/ml), subtilisin (1 µg/ml), or
-chymotrypsin (22.5 µg/ml) at room temperature for the desired time. All the digestion
reactions were stopped by adding SDS-sample loading buffer and then
heating at 95 °C for 3 min. The samples were analyzed by SDS-PAGE on
a 12.5% polyacrylamide gel, and the bands were visualized by staining
with Coomassie Blue.
Actin Polymerization--
Actin polymerization was triggered by
the addition of MgCl2 and KCl to final concentrations of 2 and 50 mM, respectively and was monitored by the increase
in light scattering as a function of time at 25 °C. The extent of
polymerization of GTP-actin was normalized against the extent of
polymerization of an equivalent concentration of ATP-actin for
comparison. To measure the critical concentration of ATP/GTP-bound
actin, increasing concentrations of actin were induced to polymerize
until the steady state was reached, and the net change in light
scattering in each experiment was recorded. The net changes were
plotted against the concentration of actin, and the critical
concentration of actin was obtained by drawing a line through the
points and determining its intersection on the x axis. To
assess the ability of the actin filament to withstand cold temperature,
9.2 µM actin was polymerized at 25 °C until the steady
state was reached. The temperature of the sample was then lowered from
25 to 5 °C using a digitally controlled thermostatted waterbath
connected to the cuvette chamber.
Phosphate Release Assay--
The inorganic phosphate released
from F-actin during polymerization was measured in a coupled assay
using the EnzChekTM assay developed by Webb (22). The assay
is based on a Pi-dependent enzymatic
conversion of 2-amino-6-mercapto-7-methylpurine riboside to
2-amino-6-mercapto-7-methylpurine and ribose 1-phospate by a
purine nucleoside phosphorylase. The production of
2-amino-6-mercapto-7-methylpurine in the reaction results in an
increased absorbance at 360 nm. Each reaction contained actin at the
desired concentration, 0.5 unit/ml of purine nucleoside phosphorylase,
and 200 µM 2-amino-6-mercapto-7-methylpurine riboside.
Pi released from actin in F buffer (5 mM
Tris-HCl, pH 7.5, 2 mM MgCl2, 50 mM
KCl, 0.2 mM CaCl2, 0.2 mM
dithiothreitol, and 0.2 mM nucleotide triphosphate) upon
actin polymerization was monitored continuously by the change in
absorbance at 360 nm. A standard curve allowed conversion of the
absorbance change to Pi concentration.
Myosin Decoration and Electron Microscopy--
To decorate actin
filaments with myosin S1, GTP-bound actin (9.2 µM) was
polymerized in G buffer with 100 µM GTP at 25 °C for
30 min and then mixed with 4.6 µM myosin S1 for another
30 min. To prepare the samples for electron microscopy, 4.6 µM F-actin or myosin-decorated actin filaments were
diluted 4-fold with F-buffer (5 mM Tris-HCl, pH 7.5, 2 mM MgCl2, 50 mM KCl, 0.1 mM GTP, and 0.1 mM dithiothreitol) and were
applied to carbon-coated Formvar grids and visualized following
negative staining with 1.5% (w/v) uranyl acetate using a Hitachi 7000 electron microscope (University of Iowa Electron Microscope Facility).
 |
RESULTS |
Relative Binding Affinity Assays--
To determine the capacity of
yeast actin to bind GTP, we measured the dissociation constant ratio of
ATP and GTP versus
ATP for the nucleotide tight binding
site in yeast actin by adding free ATP or GTP to
ATP-bound Ca-actin
in nucleotide-free G buffer. Parallel ATP and GTP titrations were
performed on the same preparation of
ATP-actin to eliminate
inconsistencies arising from differences between actin preparations.
Fig. 1, A and B
demonstrates the relative ability of GTP or ATP to compete with
ATP
for binding to Ca yeast actin. These data were then fit to Equation 2.
The results demonstrate that the dissociation constant ratio of
ATP
versus ATP,
K
ATP/KATP, is 2.3 whereas the ratio of
ATP versus GTP,
K
ATP/KGTP, is 0.25. Using these two numbers, we calculate a dissociation constant ratio of
GTP versus ATP,
KGTP/KATP, of 9.2. This
value is ~50-fold less then than the ratio, 500, previously reported
by Iyengar and Weber (3) for muscle actin.

View larger version (27K):
[in this window]
[in a new window]
|
Fig. 1.
Ability of ATP versus GTP to
compete with ATP for binding to Ca yeast and
skeletal muscle actins. Different amounts of ATP or GTP were added
to 4.6 µM yeast actin or 2.3 µM muscle
actin containing ATP, and the decrease in fluorescence was observed
until equilibrium had been reached. ATP-Actin concentrations were
calculated from the final fluorescence values as described under
"Experimental Procedures." Open circles represent the
experimental data, and the lines through the points
represent those predicted with Equation 2 under "Experimental
Procedures." Panel A, yeast actin with ATP; panel
B, yeast actin with GTP; panel C, muscle actin with
ATP; and panel D, muscle actin with GTP.
|
|
To ensure that the difference observed between yeast and muscle actins
was not simply due to differences in experimental methodologies employed in the two studies, we repeated the identical experiments with
Ca muscle actin. From the results shown in Fig. 1, C and D, we calculated the ratios for Ca muscle actin to be 6.2 for K
ATP/KATP and
0.005 for K
ATP/KGTP.
Thus, the dissociation constant ratio of GTP versus ATP,
KGTP/KATP, is 1.24 × 103 for our Ca muscle actin experiments. This number is
within about a factor of 2 of what was reported previously (3),
indicating that yeast actin is not as highly selective as muscle
actin in its ability to interact with GTP at the tight binding
pocket. The 6.2-fold stronger affinity of muscle actin for ATP compared with
ATP reported here is slightly higher than the 4-fold difference reported by others (23, 24), and the small difference in results may
reflect different methods used in the different studies.
Effect of Divalent Cation on Nucleotide Binding--
When
Mg2+ is bound in the high affinity divalent cation binding
site of actin, it changes the flexibility of the actin molecule (25).
We, thus, examined the dissociation constant ratio of GTP
versus ATP for yeast Mg-actin. Fig.
2 shows that
K
ATP/KATP is 9.44, and
K
ATP/KGTP is 0.8 for
Mg yeast actin, demonstrating an ~11-fold stronger interaction for Mg
yeast actin with ATP than GTP. The slight difference in the values of
KGTP/KATP for Ca-actin, 9.2, and Mg actin, 11, indicates that the divalent cation bound in the
cation tight binding site does not significantly alter the ability of
actin to interact with GTP.

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 2.
Ability of ATP versus GTP to
compete with ATP for binding to Mg yeast
actin. The experiment was performed as described above except that
the Ca-actin was first converted to Mg-actin as described under
"Experimental Procedures." , ATP; , GTP.
|
|
The Conformation of Subdomain 2 of GTP-bound Yeast
Ca-actin--
Limited proteolysis has been used widely to probe the
conformation of actin subdomain 2. Strzelecka-Golaszewska
et al. (26) showed that ADP-bound muscle actin is more
susceptible to trypsin digestion and less susceptible to subtilisin
digestion than ATP-bound muscle actin, indicating distinct
conformations of subdomain 2 when muscle actin is bound to ATP or ADP.
We thus used limited proteolysis to assess the effect of GTP binding on
the structure of subdomain 2. We subjected G-actin to limited
proteolysis with trypsin, subtilisin, or
-chymotrypsin, which
specifically cleave G-actin at
Arg62-Gly63/Lys68-Tyr69,
Met47-Gly48, and
Met44-Val45/Leu67-Lys68, respectively.
SDS-PAGE analysis (Fig. 3) shows no
significant difference in the digestion pattern between the GTP-bound
and ATP-bound yeast Ca-actins suggesting that the exposure of the
subdomain 2 sites is not significantly different for the two actins.
Actin crystallographic data (10, 27, 28) all show that the nucleotide
locates at the interface of actin's two major domains near the base of
the interdomain cleft; the polyphosphate part of the molecule interacts with both domains, and the base portion interacts with subdomains 3 and
4. Thus, the bound nucleotide functions as a connecting bridge between
the two domains. The amino acids rearrange when they interact with
either the tri- or diphosphate and cause a subdomain 2 conformational
change reported by Strzelecka-Golaszewska et al. (26).
Therefore, our results imply that the triphosphate of GTP locates at or
near the same position as the triphosphate of ATP.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 3.
Susceptibility of ATP versus
GTP Ca yeast actin to controlled proteolysis. The actin, at
a concentration of 9 µM, was subjected to proteolysis by
trypsin, subtilis, and chymotrypsin as described under "Experimental
Procedures," and the digestion products were resolved by SDS-PAGE
followed by visualization with Coomassie Blue. The electrophoretograms
are shown. A, ATP-actin; G, GTP-actin.
|
|
Thermal Stability of GTP-bound G-actin--
A structural change in
yeast actin subdomains 3 and 4 needed to accommodate GTP might not be
detected by limited proteolysis. Thus, we examined the thermal
stability of GTP-bound Ca yeast actin by circular dichroism. The molar
ellipticity at 222 nm, which reflects the content of
-helical
structure in G-actin, was monitored as a function of temperature. The
thermal stability of a protein is determined by the apparent melting
temperature (Tm), which is defined as the
midpoint of the thermal transition curve. Although the transition is
not reversible for actin, a comparison of the curves generated by two
different actins gives a relative indication of their thermal
stabilities. This approach was first used for actin by
Strzelecka-Golaszewska et al. (29) and Bertazzon et
al. (30) for muscle actin and validated for yeast actin by us
previously (19). The results in Fig. 4
show that the Tm is 49.5 °C for GTP-bound Ca
yeast G-actin, whereas the Tm is 59 °C for
ATP-bound G-actin. This substantial decrease in the "melting
temperature" of the actin in the presence of GTP may result either
from a weaker binding of GTP versus ATP to actin giving rise
to a less stable actin conformation at higher temperature or from a
GTP-induced conformational change in the protein itself.

View larger version (27K):
[in this window]
[in a new window]
|
Fig. 4.
Relative thermal stability of ATP
versus GTP Ca G-actin. Thermal stabilities were
measured by subjecting 1 µM of the appropriate actin to
increasing temperatures between 20 and 80 °C and monitoring the
change in ellipticity at 222 nm in a CD spectropolarimeter as described
under "Experimental Procedures." The data were normalized to
reflect the fraction of protein remaining in the native conformation at
each temperature. , ATP; , GTP.
|
|
Polymerizibility of GTP-bound Yeast Ca-actin--
To determine the
effect of the replacement of ATP with GTP on actin polymerization, we
examined the polymerizability of GTP-bound yeast Ca-actin at different
GTP concentrations at 25 °C. The starting material was a solution of
a 1:1 complex of actin and GTP from which unbound nucleotide had been
removed. Fig. 5 shows that GTP-bound yeast Ca-actin, contrary to the case for muscle actin, polymerizes in
the absence of exogenous nucleotide. Samples at the end of each
polymerization assay were observed by electron microscopy, and the
structure of GTP-bound yeast Ca-actin filaments appears similar to
ATP-bound yeast Ca-actin filaments (Fig.
6).

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 5.
Polymerization of GTP Ca yeast actin. Ca
yeast ATP- or GTP G-actin, 4.6 µM, in the presence of
either ATP or different concentrations of GTP, respectively, was
polymerized by the addition of MgCl2 and KCl to 2 and 50 mM, respectively, at 25 °C, and the change in light
scattering observed was as described under "Experimental
Procedures." 100 µM ATP (X); GTP
concentrations were as follows: 0 µM, ; 10 µM, ; and 100 µM, .
|
|

View larger version (90K):
[in this window]
[in a new window]
|
Fig. 6.
Electron micrographs of GTP and ATP yeast Ca
F-actin. At the conclusion of the experiments in Fig. 5, actin was
removed, stained with uranyl acetate, and viewed under the electron
microscope as described under "Experimental Procedures."
Panel A, ATP-actin; panel B, GTP-actin.
Bar = 50 nm.
|
|
The purine binding region of actin lies near the intersection of
subdomains 3 and 4. Because thermal denaturation results suggested that
the binding of GTP altered the conformation of actin, it is possible
that the integrity of the 3/4 boundary was affected. Previous results
from our laboratory demonstrated that mutations in this region might
produce a cold-sensitive polymerization defect (19, 21). To determine
whether the binding of GTP to actin produced such an effect, we
polymerized 9.2 µM GTP-bound yeast Ca-actin at 25 °C
until a steady state was reached. We then lowered the temperature to
5 °C and monitored light scattering as a function of temperature. We
observed no cold sensitivity of the preformed filaments (data not shown).
Critical Concentration of GTP-bound Actin--
At a concentration
of 4.6 µM, the extent of polymerization of GTP-bound
actin is ~35% less than that of ATP-bound actin as shown in Fig. 5
suggesting a significant difference in the critical concentration of
the GTP- and ATP-actin forms. We thus determined the critical
concentrations of the two actin preparations by measuring the increase
in light scattering when solutions of different actin concentrations
were allowed to polymerize and graphing the magnitude of these changes
against the actin concentration. Fig. 7
shows that although the critical concentration of ATP-actin is 0.3 µM, that of the GTP-actin is 1.5 µM. To
determine whether the GTP- and ATP-actin conformations were
interconvertible, we added ATP to a sample of GTP F-actin. Fig.
8 demonstrates that the extent of
polymerization increased to that of a parallel sample of ATP F-actin
showing the GTP did not cause irreversible changes in actin
conformation.

View larger version (14K):
[in this window]
[in a new window]
|
Fig. 7.
Critical concentration determined for ATP
versus GTP Ca yeast actin. As described under
"Experimental Procedures," different concentrations of the actins
were allowed to polymerize in the presence of 100 µM of
their respective nucleotides, and the final light scattering values
were plotted against the total actin concentration. The critical
concentration was calculated at the point of intersection of each line
with the x axis. ATP-actin, ; GTP-actin, .
|
|

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 8.
Stimulation of additional actin
polymerization by the addition of ATP to Ca GTP F-actin. Both 4.6 µM Ca ATP- ( ) and GTP- ( ) Ca yeast actin in the
presence of 10 µM of their respective nucleotide were
polymerized until the steady state was approached. At 1600 s,
equivalent amounts of 100 µM ATP was then added to each
sample, and the resulting change in light scattering followed. Both
samples reached the same final light scattering value. The sharp
deflection in the middle of the plot is because of the
opening and closing of the sample changer.
|
|
Pi Release Assay--
ATP hydrolysis and
Pi release have been proposed to be important determinants
of the rate of actin filament turnover (2). The slow polymerization
rate and the higher critical concentration of GTP-bound actin described
previously suggest that polymerization-induced hydrolysis of bound
nucleotide might be affected. Therefore, we assessed the ability of
actin to hydrolyze bound GTP. Yao and Rubenstein (31) reported that
ATPase activity and Pi release are concurrent during yeast
actin polymerization; there is no appreciable retention of
Pi following nucleotide hydrolysis contrary to what is
observed with muscle actin. Therefore, as an indication of actin's
ability to hydrolyze GTP, we monitored Pi release from GTP-bound actin following the induction of polymerization. Fig. 9 shows that until polymerization reaches
a steady state, Pi release and actin polymerization curves
with either GTP-bound actin or ATP-bound actin are superimposable. The
continued release of Pi following attainment of the
polymerization steady state results from continued addition of
nucleotide triphosphate-actin monomers to and release of nucleotide
diphosphate monomers from the ends of the filaments in the absence of
net filament growth. The greater divergence of the two curves in the
ATP versus GTP plots may reflect the differences in the
relative on and off rates for GTP versus ATP monomers. These
results indicate that the mode of binding of GTP to yeast actin places
the
-phosphate in a position that allows normal hydrolysis to
occur.

View larger version (20K):
[in this window]
[in a new window]
|
Fig. 9.
Pi release during the
polymerization of GTP and ATP Ca yeast actin. Either ATP
(panel A) or GTP (panel B) Ca yeast actin (6.9 µM) was induced to polymerize. Polymerization ( ) was
followed by the increase in light scattering, and Pi
release following nucleotide hydrolysis ( ) was monitored as
described under "Experimental Procedures." The concentration of
Pi released was based on use of a standard curve, and the
percent of actin polymerized was based on the final magnitude of change
in light scattering representative of the polymerization of the entire
amount of actin in the sample.
|
|
Interaction of GTP-bound Yeast Actin with Myosin S1--
GTP-actin
appears to form filaments similar to ATP-actin filaments. If their
structure is normal, they should be able to interact with myosin S1 to
generate the appearance of arrowheads on the filament surface. Fig.
10 shows an electron micrograph of
GTP-actin filaments decorated with myosin S1 head in a typical
arrowhead pattern, further indicating the normal structure of these
filaments. Significantly, this experiment was carried out in the
presence of exogenous GTP to help stabilize the actin, and this
decoration was prevented by the presence of exogenous ATP instead of
GTP because of the ability of ATP to displace S1 from actin.

View larger version (152K):
[in this window]
[in a new window]
|
Fig. 10.
Decoration of GTP-actin with myosin S1.
Actin-myosin complexes in the presence of GTP were prepared as
described under "Experimental Procedures" and applied to a
carbon-coated Formvar grid. The preparation was examined under the
electron microscope following negative staining with 1.5% uranyl
acetate. The arrows indicate the direction in which the
myosin arrows are pointing, and the bar represents a
distance of 50 nm.
|
|
Molecular Basis for Determination of Nucleotide Binding
Specificity--
The crystal structure of muscle actin shows that
three residues in the protein provide the principal points of contact
with the adenine ring of the bound ATP or ADP: Glu214,
Lys336, and Tyr306 (32). In yeast actin,
Tyr306 is replaced with a Phe. We used the Swiss-PDB Viewer
3.7 program (www.expasy.org/spdbv/mainpage.html) to fit the
residues around the adenine binding pocket in yeast (PDB number 1YAG)
and muscle actin (PDB number 1ATN) and found no significant
differences in the nature of the contacts with the adenine. We then
used the program to calculate cavities of 1.4 Å or greater diameter in the adenine binding pockets of these two proteins. Fig.
11 shows that a cavity with a surface
area of 35 Å2 and a volume of 16 Å3 is
adjacent to the adenine in yeast actin (panel A). This
cavity is not present in either muscle actin (panel B) or in
-nonmuscle actin in either the tight or open conformation (data not
shown). This cavity appears to be at the position usually occupied by the Tyr306 hydroxyl group in a position that might allow
the C2 amino group of guanine to insert with a minimum amount of
distortion of the protein structure. To further substantiate our
hypothesis, we substituted a Tyr for Phe306 in the yeast
actin structure in the same rotamer conformation and recomputed the
presence of a possible cavity. The cavity was no longer present,
strongly suggesting that this difference between the yeast and muscle
actins may play a significant role in the relatively greater ability of
yeast actin to function with GTP instead of ATP.

View larger version (21K):
[in this window]
[in a new window]
|
Fig. 11.
Phe306 in yeast actin creates a
cavity, not present in muscle actin with Tyr306, that can
potentially accommodate a guanine amino group. The region
in the vicinity of Phe306 of yeast actin (panel
A) and Tyr306 of muscle actin (panel B)
were visualized using Swiss-PDB Viewer 3.7. The surface area of the
region was calculated for both actins using a probe of 1.4 Å, and a
cavity of 16 Å3 in yeast actin was calculated as a
difference between the two structures.
|
|
To test this hypothesis, we used site-directed mutagenesis to
substitute Tyr for Phe at position 306 of yeast actin. Cells expressing
this mutant as the sole actin were easily obtained, and they exhibited
no obvious defect in growth on complex medium, hyperosmolar medium, or
on medium containing glycerol as the carbon source. Actin was purified
from these cells and tested for its relative ability to bind GTP
versus ATP using the
ATP displacement assay employed
earlier. The results in Fig. 12 show
that the ratio of GTP to ATP dissociation constants goes from 9.2:1
observed with the Phe306 actin to about 125:1 with the
Tyr306. Thus, with our assay system, this single
substitution regains one order of magnitude of the two-order of
magnitude difference in nucleotide binding specificity observed with
muscle actin.

View larger version (13K):
[in this window]
[in a new window]
|
Fig. 12.
Ability of ATP versus GTP
to compete with ATP for binding to Ca F306Y
yeast actin. Different amounts of ATP or GTP were added to 2.3 µM actin containing ATP, and the decrease in
fluorescence was observed until equilibrium had been reached.
ATP-Actin concentrations were calculated from the final fluorescence
values as described under "Experimental Procedures." Open
circles represent the experimental data, and the lines
through the points represent those predicted with Equation 2 under
"Experimental Procedures." Panel A, yeast F306Y actin
with ATP; panel B, yeast F306Y actin with GTP.
|
|
 |
DISCUSSION |
Actin had been thought to be highly selective toward adenine
versus guanine nucleotides in terms of the nucleotide it
could use for structural stabilization and possibly as a determinant of
filament turnover within the cell. Our results with yeast actin in
terms of affinity, polymerizability, and nucleotide hydrolysis demonstrate that among actins, this selectivity is certainly not universal. There is only about a 9-fold difference in the ability of
yeast actin to use GTP versus ATP instead of the 1,000-fold difference observed with muscle actin. However, with GTP the behavior of yeast actin is somewhat affected. Overall polymerization from initiation of the process to establishment of the steady state is
slower. We did not dissect the process further to rigorously differentiate between the effects of the nucleotide on nucleation versus elongation. The critical concentration is higher, and
thermal stability of the monomer is lower for GTP- compared with
ATP-actin. These differences suggest that whereas the yeast actin
structure can more easily accommodate the guanine than can that of
muscle actin, the guanine nucleotide must still cause a distortion in the normal structure of the protein.
Based on an analysis of the difference in the structures of yeast and
muscle actin, we hypothesized that the relaxed specificity of yeast
actin might be due, in part, to the presence of Phe at residue 306 instead of the Tyr found in muscle actin. The absence of the Tyr -OH
group appeared to create a cavity in which the side chain amino group
of the guanine might insert. Conversion of this Phe to Tyr in yeast
actin resulted in a recovery of one of the two-orders of magnitude
difference in nucleotide selectivity that normally differentiate yeast
from muscle actin, in agreement with our hypothesis.
The Tyr may potentially contribute to a limitation of nucleotide
accessibility in another way by restricting the conformational flexibility of the actin monomer. Examination of the crystal structures of the two proteins shows that if bound waters are included, a single
water molecule bridges the Tyr306 hydroxyl to a number of
polar residues such as Lys213 in subdomain 4 with links to
the top of this subdomain and Thr303 in subdomain 3, along
with the N3 of adenine, part of a ligand that forms a bridge between
the two major domains of the protein (Fig.
13). In yeast actin, where Phe takes
the place of Tyr, this extensive network of interacting residues is
interrupted, possibly facilitating accommodation of the guanine
nucleotide within the cleft of the protein.

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 13.
Hydrogen bonding network involving
Tyr306 and its bound water located in the adenine binding
region of muscle actin. The figure was generated by the same
procedure as used for Fig. 11 to emphasize the hydrogen bond network
(green dashed lines) in the area of the adenine involving
Tyr306. Substitution of this residue with Phe in yeast
actin presumably would weaken subdomain 3/4 interactions and the entire
bonding network leading to a conformationally more flexible molecule.
Modified from PDB number 1J6Z (40).
|
|
Another factor that may play a major contribution in accommodating the
GTP in yeast actin is a propensity for the yeast actin monomer to
assume a more open conformation than muscle actin. Actin has been
crystallized as a monomer in two conformations, a more closed
conformation which is most often seen (33-35) and a more open
conformation relative to the size of the interdomain nucleotide binding
cleft (10). A number of studies, both biochemical and structural, have
suggested that yeast actin may be more prone to the open conformation
compared with muscle actin. Optical reconstructions of electron
micrographs of yeast actin filaments (8) reveal a decreased density in
the nucleotide pocket of monomers within the filament that can best be
modeled by the open monomer conformation described by Chik et
al. (10) whereas muscle actin filaments are best modeled using the
closed conformation. Yeast actin polymerizes more rapidly than does
muscle actin, a phenomenon that may result from more rapid filament
nucleation (7, 17), and it exchanges nucleotide 10 times faster than
does muscle actin (5, 31). Whereas muscle actin retains the
Pi generated by polymerization-dependent hydrolysis of ATP for a substantial time before releasing it into the
solution to generate ADP F-actin (36), the release of Pi from yeast actin is almost immediate following its hydrolysis of the
bound nucleotide (31). Finally, in a series of mutagenesis experiments
with yeast actin involving the conserved His73 near the
bottom of the nucleotide cleft, Yao et al. (31) demonstrated that the H73E mutation caused several polymerization and monomer structural abnormalities in yeast actin. Two acidic residues, Asp179 and Asp184, on the opposite side of a
cleft are potentially in a position to interact electrostatically with
the cationic His73. In an H73E mutant, we asked which of
the two residues when converted to Arg, 184 or 179, would be most
effective in restoring the H73E mutant to a wild-type phenotype. A
D179R mutation, as opposed to a D184R mutation, was much more effective
in restoring the wild-type properties to the H73E actin, and the
H73E/D179R interaction appeared, on the basis of molecular dynamics
modeling, to be more highly favored in the open conformation (37).
The biological significance of the ability of yeast actin to function
efficiently with GTP is unclear. The concentration of GTP in the cell
is ~0.2 mM, and that of ATP is about 1 mM
(38). Coupled with the roughly 10-fold difference in yeast actin for GTP versus ATP, on the average no more than about 2% of the
actin would be expected to be in the GTP form. Unless a particular
actin-binding protein had a particular affinity for GTP
versus ATP-actin, or unless there were relatively large
local concentrations of GTP in the cell, the impact of the GTP-actin
that did form should be minor. However, the ability of yeast to utilize
GTP, the small but significant effects it exerts on some parameters of
actin function, and the availability of different fluorescently labeled yeast actins that have been generated may lead to a unique insight into
the way in which the bound nucleotide controls domain/domain interactions in yeast actin and the processes that depend on them. Realization of the potential importance of a Phe versus Tyr
at residue 306 in controlling conformational flexibility of subdomains 3 and 4 will also provide an avenue for obtaining new insight into the
importance of this interacting network in regulating polymerization and
filament dynamics.
Although the significance of the ability of yeast actin to more easily
utilize GTP is not clear, the evolutionary segregation of the Phe and
Tyr residues at position 306 strongly suggests that the identity of the
amino acid here is important to the organism. Based on actin sequences
deposited in the NIH protein database (www.ncbi.nih.gov/entrez/),
metazoans including insects, frogs, worms, sea urchins, birds, fish,
and mammals have a Tyr at residue 306. On the other hand, actins from
plants including chlamydomonas, slime molds including physarum,
dictyostelium, and acanthamoeba, and the protist Naegleria
gruberi all carry a Phe at position 306. Fungal actins display
both alleles. Filamentous fungi including neurospora, aspergillus, and
penicillium all have Tyr. The budding yeasts Saccharomyces
cerevisiae, Candida albicans, Pichia
pastoris, and Kluyveromyces lactis all have Phe.
Interestingly, the actin from Schizosaccharomyces pombe, a
fission yeast, carries a Tyr and not the Phe found in the budding
yeast. The line ending in S. pombe, believed to have
diverged from that of S. cerevisiae, a budding yeast, around
400 million years ago, is thought to have undergone rapid evolutionary
change leading to the organism we know today (39). Whether other
fission yeasts share the same sequence at this position could not be
determined. The pervasiveness of Tyr306 among the metazoans
and Phe306 among the plants suggests that in these
evolutionarily advanced groups of organisms, the particular allele has
been fixed, because it offers some selective advantage to the cells
that reflect the particular development of the cytoskeletal machinery
in the cells in which these actins are expressed.
 |
FOOTNOTES |
*
This work was supported in part by National Institutes of
Health Grants AI45728 and GM33689 (to P. A. 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.
To whom correspondence should be addressed. Tel.: 319-335-7911;
Fax: 319-335-9570; E-mail: peter-rubenstein@uiowa.edu.
Published, JBC Papers in Press, August 20, 2002, DOI 10.1074/jbc.M204025200
 |
ABBREVIATIONS |
The abbreviation used is:
ATP, 1,N6-ethenoadenosine 5'- triphosphate.
 |
REFERENCES |
| 1.
|
Carlier, M. F.
(1991)
Curr. Opin. Cell Biol.
3,
12-17[CrossRef][Medline]
[Order article via Infotrieve]
|
| 2.
|
Carlier, M. F.
(1991)
J. Biol. Chem.
266,
1-4[Free Full Text]
|
| 3.
|
Iyengar, M.,
and Weber, H. H.
(1964)
Biochim. Biophys. Acta
86,
543-553[Medline]
[Order article via Infotrieve]
|
| 4.
|
Martonosi, A.,
and Gouvea, M. A.
(1961)
J. Biol. Chem.
236,
1345-1352[Free Full Text]
|
| 5.
|
Chen, X.,
Peng, J.,
Pedram, M.,
Swenson, C. A.,
and Rubenstein, P. A.
(1995)
J. Biol. Chem.
270,
11415-11423[Abstract/Free Full Text]
|
| 6.
|
Kim, E.,
Miller, C. J.,
and Reisler, E.
(1996)
Biochemistry
35,
16566-16572[CrossRef][Medline]
[Order article via Infotrieve]
|
| 7.
|
Buzan, J. M.,
and Frieden, C.
(1996)
Proc. Natl. Acad. Sci. U. S. A.
93,
91-95[Abstract/Free Full Text]
|
| 8.
|
Belmont, L. D.,
Orlova, A.,
Drubin, D. G.,
and Egelman, E. H.
(1999)
Proc. Natl. Acad. Sci. U. S. A.
96,
29-34[Abstract/Free Full Text]
|
| 9.
|
Yao, X.,
Grade, S.,
Wriggers, W.,
and Rubenstein, P. A.
(1999)
J. Biol. Chem.
274,
37443-37449[Abstract/Free Full Text]
|
| 10.
|
Chik, J. K.,
Lindberg, U.,
and Schutt, C. E.
(1996)
J. Mol. Biol.
263,
607-623[CrossRef][Medline]
[Order article via Infotrieve]
|
| 11.
|
Nyman, T.,
Schuler, H.,
Korenbaum, E.,
Schutt, C. E.,
Karlsson, R.,
and Lindberg, U.
(2002)
J. Mol. Biol.
317,
577-589[CrossRef][Medline]
[Order article via Infotrieve]
|
| 12.
|
Hurley, J. H.
(1996)
Annu. Rev. Biophys. Biomol. Struct.
25,
137-162[Medline]
[Order article via Infotrieve]
|
| 13.
|
Liberek, K.,
Skowyra, D.,
Zylicz, M.,
Johnson, C.,
and Georgopoulos, C.
(1991)
J. Biol. Chem.
266,
14491-14496[Abstract/Free Full Text]
|
| 14.
|
Nagayama, F.,
Ohshima, H.,
Suzuki, H.,
and Ohshima, T.
(1980)
Biochim. Biophys. Acta
615,
85-93[Medline]
[Order article via Infotrieve]
|
| 15.
|
van den Ent, F.,
Amos, L. A.,
and Lowe, J.
(2001)
Nature
413,
39-44[CrossRef][Medline]
[Order article via Infotrieve]
|
| 16.
|
Smith, N. H.,
Maynard Smith, J.,
and Spratt, B. G.
(1995)
Mol. Biol. Evol.
12,
363-370[Abstract]
|
| 17.
|
Cook, R. K.,
Blake, W. T.,
and Rubenstein, P. A.
(1992)
J. Biol. Chem.
267,
9430-9436[Abstract/Free Full Text]
|
| 18.
|
Spudich, J. A.,
and Watt, S.
(1971)
J. Biol. Chem.
246,
4866-4871[Abstract/Free Full Text]
|
| 19.
|
Chen, X.,
Cook, R. K.,
and Rubenstein, P. A.
(1993)
J. Cell Biol.
123,
1185-1195[Abstract/Free Full Text]
|
| 20.
|
Klumb, L. A.,
Chu, V.,
and Stayton, P. S.
(1998)
Biochemistry
37,
7657-7663[CrossRef][Medline]
[Order article via Infotrieve]
|
| 21.
|
Kuang, B.,
and Rubenstein, P. A.
(1997)
J. Biol. Chem.
272,
4412-4418[Abstract/Free Full Text]
|
| 22.
|
Webb, M. R.
(1992)
Proc. Natl. Acad. Sci. U. S. A.
89,
4884-4887[Abstract/Free Full Text]
|
| 23.
|
Kinosian, H. J.,
Selden, L. A.,
Gershman, L. C.,
and Estes, J. E.
(2000)
Biochemistry
39,
13176-13188[CrossRef][Medline]
[Order article via Infotrieve]
|
| 24.
|
Vinson, V. K., De La,
Cruz, E. M.,
Higgs, H. N.,
and Pollard, T. D.
(1998)
Biochemistry
37,
10871-10880[CrossRef][Medline]
[Order article via Infotrieve]
|
| 25.
|
Orlova, A.,
and Egelman, E. H.
(1993)
J. Mol. Biol.
232,
334-341[CrossRef][Medline]
[Order article via Infotrieve]
|
| 26.
|
Strzelecka-Golaszewska, H.,
Moraczewska, J.,
Khaitlina, S. Y.,
and Mossakowska, M.
(1993)
Eur. J. Biochem.
211,
731-742[Medline]
[Order article via Infotrieve]
|
| 27.
|
Flaherty, K. M.,
McKay, D. B.,
Kabsch, W.,
and Holmes, K. C.
(1991)
Proc. Natl. Acad. Sci. U. S. A.
88,
5041-5045[Abstract/Free Full Text]
|
| 28.
|
Cedergren-Zeppezauer, E. S.,
Goonesekere, N. C.,
Rozycki, M. D.,
Myslik, J. C.,
Dauter, Z.,
Lindberg, U.,
and Schutt, C. E.
(1994)
J. Mol. Biol.
240,
459-475[CrossRef][Medline]
[Order article via Infotrieve]
|
| 29.
|
Strzelecka-Golaszewska, H.,
Venyaminov, S.,
Zmorzynski, S.,
and Mossakowska, M.
(1985)
Eur. J. Biochem.
147,
331-342[Medline]
[Order article via Infotrieve]
|
| 30.
|
Bertazzon, A.,
Tian, G. H.,
Lamblin, A.,
and Tsong, T. Y.
(1990)
Biochemistry
29,
291-298[CrossRef][Medline]
[Order article via Infotrieve]
|
| 31.
|
Yao, X.,
and Rubenstein, P. A.
(2001)
J. Biol. Chem.
276,
25598-25604[Abstract/Free Full Text]
|
| 32.
|
Kabsch, W.,
and Holmes, K. C.
(1995)
FASEB J.
9,
167-174[Abstract]
|
| 33.
|
Kabsch, W.,
Mannherz, H. G.,
Suck, D.,
Pai, E. F.,
and Holmes, K. C.
(1990)
Nature
347,
37-44[CrossRef][Medline]
[Order article via Infotrieve]
|
| 34.
|
Mannherz, H. G.,
Gooch, J.,
Way, M.,
Weeds, A. G.,
and McLaughlin, P. J.
(1992)
J. Mol. Biol.
226,
899-901[CrossRef][Medline]
[Order article via Infotrieve]
|
| 35.
|
Schutt, C. E.,
Myslik, J. C.,
Rozycki, M. D.,
Goonesekere, N. C.,
and Lindberg, U.
(1993)
Nature
365,
810-816[CrossRef][Medline]
[Order article via Infotrieve]
|
| 36.
|
Melki, R.,
Fievez, S.,
and Carlier, M. F.
(1996)
Biochemistry
35,
12038-12045[CrossRef][Medline]
[Order article via Infotrieve]
|
| 37.
|
Yao, X.,
Nguyen, V.,
Wriggers, W.,
and Rubenstein, P. A.
(2002)
J. Biol. Chem.
277,
22875-22882[Abstract/Free Full Text]
|
| 38.
|
Rudoni, S.,
Colombo, S.,
Coccetti, P.,
and Martegani, E.
(2001)
Biochim. Biophys. Acta
1538,
181-189[Medline]
[Order article via Infotrieve]
|
| 39.
|
Wood, V.,
Gwilliam, R.,
Rajandream, M. A.,
Lyne, M.,
Lyne, R.,
Stewart, A.,
Sgouros, J.,
Peat, N.,
Hayles, J.,
Baker, S.,
Basham, D.,
Bowman, S.,
Brooks, K.,
Brown, D.,
Brown, S.,
et al..
(2002)
Nature
415,
871-880[CrossRef][Medline]
[Order article via Infotrieve]
|
| 40.
|
Otterbein, L. R.,
Graceffa, P.,
and Dominguez, R.
(2001)
Science
293,
708-711[Abstract/Free Full Text]
|
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: