JBC

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Originally published In Press as doi:10.1074/jbc.M200713200 on April 8, 2002

J. Biol. Chem., Vol. 277, Issue 24, 21514-21521, June 14, 2002
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
277/24/21514    most recent
M200713200v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mezgueldi, M. E.
Right arrow Articles by Ostap, E. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mezgueldi, M. E.
Right arrow Articles by Ostap, E. M.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

The Kinetic Mechanism of Myo1e (Human Myosin-IC)*

Mohammed El MezgueldiDagger §, Nanyun TangDagger , Steven S. Rosenfeld, and E. Michael OstapDagger ||

From the Dagger  Department of Physiology and The Pennsylvania Muscle Institute, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania 19104-6085 and the  Department of Neurology, University of Alabama, Birmingham, Alabama 35294

Received for publication, January 23, 2002, and in revised form, April 5, 2002

    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Myo1e is the widely expressed subclass-1 member of the myosin-I family. We performed a kinetic analysis of a truncated myo1e that consists of the motor and the single IQ motif with a bound calmodulin. We determined the rates and equilibrium constants for the key steps in the ATPase cycle. The maximum actin activated ATPase rate (Vmax) and the actin concentration at half-maximum of Vmax (KATPase) of myo1e are similar to those of the native protein. The KATPase is low (~1 µM), however the affinity of myo1e for actin in the presence of ATP is very weak. A weak actin affinity and a rapid rate of phosphate release result in a pathway under in vitro assay conditions in which phosphate is released while myo1e is dissociated from actin. Actin activation of the ATPase activity and the low KATPase are the result of actin activation of ADP release. We propose that myo1e is tuned to function in regions of high concentrations of cross-linked actin filaments. Additionally, we found that ADP release from actomyo1e is > 10-fold faster than other vertebrate myosin-I isoforms. We propose that subclass-1 myosin-Is are tuned for rapid sliding, whereas subclass-2 isoforms are tuned for tension maintenance or stress sensing.

    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Myosin-I motors are the single-headed, low molecular weight, membrane-associated members of the myosin superfamily (1). The molecular details of myosin-I function in the cell are not known, but myosin-Is are thought to play important roles in cortical membrane tension (2), endocytosis and endocytic trafficking (3-6), signal transduction (7), and membrane ruffling (6, 8).

Myosin-Is are the most diverse of the unconventional myosins and are represented by at least two phylogenetically distinct subclasses based on sequence comparison of motor domains (9). Subclass-1 myosin-I isoforms have long tails that contain lipid binding (TH1), proline-rich (TH2), and Src homology-3 (TH3) domains. Subclass-1 isoforms were first discovered in lower eukaryotes (10) and are sometimes termed "amoeboid" isoforms. However, vertebrate members of this subclass have been identified (11, 12). Subclass-2 myosin-I isoforms have short tails that contain only TH1 domains and are also widely expressed (9).

The kinetic mechanisms of all characterized myosin-Is follow the same pathway with the same biochemical intermediates (13-16). However, considerable kinetic variability exists within the myosin-I family. Key rate constants of isoforms from lower eukaryotes are significantly faster (3-10-fold) than those of the vertebrate isoforms. The characterized lower eukaryote isoforms belong to subclass-1 (13), and the characterized vertebrate isoforms are members of subclass-2 (14-16), thus suggesting the possibility that kinetic properties correlate with the isoform subclass rather than with the specific organism. Measurement of the rate constants of a vertebrate subclass-1 isoform will help determine whether kinetic properties are conserved within the subclasses. Additionally, knowing the rate constants of the different isoforms allows us to understand the physical limitations and properties of myosin-I isoforms and thus provides insight into the cellular functions of the subclasses. Comparisons of the sequences of myosin-I subclasses may also help unravel the structural modifications required for kinetic tuning of all myosin motors.

Myo1e, also known as myr-3 (12) and human myosin-IC (11), is the subclass-1 myosin-I expressed in nearly all vertebrate cells. Localization studies show that myo1e is found in regions of high actin concentration (17, 18), and steady-state biochemical investigations suggest that myo1e is capable of cross-linking actin filaments via an ATP-insensitive actin binding site in its tail (19). This cross-linking results in triphasic actin activated ATPase activity. Unlike myosin-Is from lower eukaryotes, the ATPase activity of myo1e appears to be regulated allosterically by its tail domain (19), and the ATPase activity is not predicted to be regulated by heavy chain phosphorylation by p21-activated kinase homologs (20). We measured all of the key rate constants in the myo1e ATPase pathway and found that, like subclass-1 isoforms from lower eukaryotes, the rate constants are significantly faster than the characterized vertebrate subclass-2 isoforms, thus suggesting conservation of kinetics and function within a single subclass across phylogenetic boundaries. However, we also found differences in the myo1e kinetics which may tune the motor for function in regions of high actin concentration.

    EXPERIMENTAL PROCEDURES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Reagents, Proteins, and Buffers-- N-Methyanthraniloyl (mant)1 nucleotide derivatives were purchased from Molecular Probes, and 2'-deoxymantATP (dmantATP) was synthesized as described by Hiratsuka (21). ADP and ATP concentrations were determined spectrophotometrically before each experiment by absorbance at 259 nm, epsilon 259 = 15,400 M-1 cm-1. MantATP and mantADP concentrations were determined by absorbance at 255 nm, epsilon 255 = 23,300 M-1 cm-1 (21).

Rabbit skeletal muscle actin was prepared and gel filtered (22). Actin concentrations were determined by absorbance at 290 nm, epsilon 259 = 26,600 M-1 cm-1 (23). Actin was labeled with pyrenyl iodoacetamide (pyrene-actin) and gel filtered (24). All actin was stabilized with a molar equivalent of phalloidin (Sigma). Calmodulin (CaM) was expressed in bacteria and purified as described (25).

Steady-state and transient experiments were performed at 25 ± 0.1 °C in KMg50 buffer (10 mM imidazole, pH 7.0, 50 mM KCl, 1 mM MgCl2, 1 mM EGTA, 1 mM DTT (26-28)) or KMg0 buffer (10 mM imidazole, pH 7.0, 1 mM MgCl2, 1 mM EGTA, 1 mM DTT). Free CaM (1-3 µM) was included in all solutions.

Myosin-I Expression and Purification-- The cDNA for human myosin-IC (accession NM 004998), kindly provided by W. M. Bement (University of Wisconsin), was truncated at Glu720, generating a construct containing the motor domain and the only IQ motif (referred to as myo1eIQ throughout the paper). A FLAG peptide sequence was inserted at the C terminus and subcloned into the baculovirus transfer vector pVL1392 (Invitrogen). Recombinant baculovirus was generated using standard procedures.

Myo1eIQ with bound calmodulin was purified from Sf9 cells that were coinfected with virus containing recombinant myo1eIQ and CaM (Fig. 1). Four liters of log phase cells (2 × 106 cells/ml) were infected and incubated at 27 °C for 60 h with shaking. Cells were harvested by centrifugation; suspended in lysis buffer (10 mM Tris, pH 7.0, 200 mM NaCl, 2 mM MgCl2, 5 mM DTT, 1 mM phenylmethylsulfonyl fluoride, 0.01 mg/ml aprotinin, 0.01 mg/ml leupeptin), 2 mM MgATP, and 0.05% Igepal at 4 °C; and homogenized with five strokes in a Dounce homogenizer. Cell extract was centrifuged at 100,000 × g for 1 h. The supernatant was loaded onto anti-FLAG antibody columns (Sigma). Columns were washed with 5 column volumes of lysis buffer + 2 mM MgATP and 5 column volumes of lysis buffer. Myo1eIQ was eluted with 10 mM Tris, pH 8.0, 100 mM NaCl, 1 mM DTT, 5 µM CaM, 0.2 mg/ml FLAG peptide (Sigma), 0.01 mg/ml aprotinin, 0.01 mg/ml leupeptin. Eluted protein was loaded directly on to an 8-ml Mono Q column (Amersham Biosciences) equilibrated in column buffer (10 mM Tris, pH 8.0, 25 mM KCl, 1 mM DTT) and eluted with a linear 25 mM-1 M KCl gradient. The Mono Q column separated myo1eIQ from FLAG peptide, ADP, ATP, and free CaM. Fractions containing myo1eIQ were dialyzed versus storage buffer (10 mM Tris, pH 7.5, 100 mM KCl, 1 mM EGTA, 1 mM DTT, 50% glycerol), which concentrated the protein and allowed for storage at -20 °C. Quantitative densitometry showed that myo1eIQ was >95% pure, and the CaM:myo1eIQ ratio was 1:1. Approximately 2 mg of pure myo1e protein was obtained from 4 liters of cells.



<UP><SC>Scheme</SC> 1</UP>



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 1.   Purification of myo1eIQ. Coomassie-stained SDS-polyacrylamide gel showing purified myo1eIQ. The associated light chain is CaM. Quantitative densitometry shows that the CaM:heavy chain stoichiometry was 1:1. The predicted molecular mass of the heavy chain from the amino acid sequence is 83 kDa.

Steady-state ATPase Activity and Sedimentation Assays-- Steady-state ATPase activities were measured in KMg50 or KMg0 buffer at 25 °C using the NADH-coupled assay as described (27). Steady-state binding of 0.2 µM myo1eIQ to 0-50 µM actin was measured in KMg50 and KMg0 at 25 °C in the presence of 2 mM MgATP by ultracentrifugation assays (350,000 × g for 20 min). The fraction of actin-bound myosin was determined by assaying the NH4/EDTA ATPase activity of the supernatant (29).

Stopped Flow, Quenched Flow, and Kinetic Modeling-- Transient kinetic measurements were made at 25 °C with an Applied Photophysics (Surrey, UK) SX.18MV stopped flow having a 1.2-ms dead time. Tryptophan fluorescence (lambda ex = 295 nm) was measured using a 320 nm WG long pass emission filter (Oriel). A 400 nm long pass filter (Oriel) was used to monitor pyrene (lambda ex = 365 nm), and mantADP and mantATP (lambda ex = 295 nm) fluorescence. Usually three to five transients were averaged before nonlinear least square fitting. The time courses presented in the figures show the average of one to four individual traces. Transients were fitted to exponential functions using the software supplied with the stopped flow. Transient Pi release was measured using the coupled assay system containing the fluorescently labeled mutant of the phosphate-binding protein (PiBP (30, 31)) with the stopped flow in sequential mixing mode using an excitation wavelength of 425 nm and a 440 nm long pass filter. The dead time of the instrument in this configuration was ~2 ms. Quenched flow measurements were performed with a KinTek (State College, PA) RQF-3 apparatus. Errors reported are standard errors in the fits.

Kinetic modeling and simulations were performed using Scheme 1 of the actomyosin ATPase, where A is actin and M is myosin. Computer simulations were performed with KSIM (Neil C. Miller) or Berkeley-Madonna (Robert I. Macey and George F. Oster; University of California, Berkeley).

    RESULTS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Steady-state ATPase and Actin Binding Parameters-- The ATPase activities of subclass-1 myosin-I isoforms from all lower eukaryotes are activated by phosphorylation of serine or threonine at the TEDS site in an actin binding, surface loop by proteins that have homology to p21-activated kinases (20). Although myo1e has serines immediately adjacent to the TEDS site, it is not predicted to be regulated by heavy chain phosphorylation (20). We were unable to phosphorylate myo1eIQ with Pak3 or Acanthamoeba myosin-I heavy chain kinase (not shown), thus confirming the TEDS rule (20).

The steady-state ATPase activity of myo1eIQ is activated ~2-fold by actin filaments. Steady-state (Fig. 2) and single turnover measurements (not shown) gave rates of ~0.6 s-1 in the absence of actin (Table I) in KMg50. The actin concentration at half-maximum of the steady-state ATPase rate (KATPase) is 1.2 ± 0.36 µM, and the maximum ATPase rate (Vmax) at saturating actin is 1.2 ± 0.03 s-1. The triphasic ATPase activity observed with native myo1e was not observed (19).


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 2.   Steady-state ATPase and actin binding activity of myo1eIQ. The actin concentration dependence of the steady-state turnover of 50 nM myo1eIQ in KMg0 (black-square) and KMg50 () was measured using the NADH-coupled assay. The data points represent acquisitions from three different preparations, and the solid lines are best fits to rectangular hyperbolas. The fraction of myo1eIQ bound to actin in KMg0 (black-diamond ) in 2 mM ATP was measured by cosedimentation.

                              
View this table:
[in this window]
[in a new window]
 
Table I
Steady-state ATPase parameters for myo1eIQ
Conditions were 1 mM MgCl2, 1 mM EGTA, 1 mM DTT, 2 mM MgATP, 10 mM imidazole, pH 7.0, 25 °C.

Steady-state actin binding measurements show that myo1eIQ is dissociated predominantly from 0-50 µM actin in KMg0 (Fig. 2) and KMg50 (not shown) in the presence of 2 mM MgATP. Therefore, in the presence of ATP, the predominant steady-state myo1eIQ intermediates have a very weak affinity for actin with an equilibrium dissociation constant (Kd) greater than 50 µM.

Myo1eIQ Binding to Actin Filaments-- An ~80% fluorescence quenching upon strong binding of myo1eIQ to pyrene-actin allowed us to monitor the association of myo1eIQ with actin in the absence (Reaction 1) and presence (Reaction 2) of MgADP. Data were modeled as


<UP>AM</UP> <LIM><OP><ARROW>⇄</ARROW></OP><LL>k<SUB>−7</SUB></LL><UL>k<SUB>+7</SUB></UL></LIM> <UP>A</UP>*+<UP>M</UP>

<UP><SC>Reaction</SC> 1</UP>
and
<UP>AM·ADP</UP> <LIM><OP><ARROW>⇄</ARROW></OP><LL>k<SUB>−11</SUB></LL><UL>k<SUB>+11</SUB></UL></LIM> <UP>A</UP>*+<UP>M·ADP</UP>

<UP><SC>Reaction</SC> 2</UP>
where A* represents the unquenched fluorescent state of pyrene-actin. Time courses of myo1eIQ and myo1eIQ-ADP binding to pyrene-actin follow single exponentials rates that depend linearly on the actin concentration (Fig. 3A). The apparent second order rate constant for myo1eIQ binding to actin obtained from the slope is k-7 = 9.0 ± 0.7 µM-1 s-1 (Table II). The apparent second order rate constant for myo1eIQ-ADP binding to actin is k-11 = 1.2 µM-1 s-1 (Table II).


View larger version (20K):
[in this window]
[in a new window]
 
Fig. 3.   Kinetics of myo1eIQ association with actin filaments. A, concentration dependence of the observed rate (kobs) of pyrene-actin binding to 0.2 µM myo1eIQ in the absence (black-square) and presence () of 1 mM MgADP in KMg50. Solids lines are linear fits to the data. B, time course of pyrene-actin fluorescence increase after mixing 0.4 µM pyrene-actomyo1eIQ with 50 µM actin in the absence (upper curve) and presence (lower curve) of 1 mM MgADP. Smooth lines are the best fits of the data to single exponentials.

                              
View this table:
[in this window]
[in a new window]
 
Table II
Rate and equilibrium constants for the myo1eIQ ATPase cycle
KMg50 (50 mM KCl, 1 mM MgCl2, 1 mM EGTA, 1 mM DTT, 10 mM imidazole, (pH 7.0, 25 °C).

The rate of dissociation of myo1eIQ and myo1eI-ADP from pyrene-actin was measured by competition with 60-fold excess of unlabeled actin (Fig. 3B). Time courses fit single exponentials with rates k+7 = 0.056 s-1 and k+11 = 0.065 s-1 (Table II). Dissociation equilibrium constants (K7 and K11) calculated from the dissociation and association rates (k+7/k-7 and k+11/k-11) are K7 = 6.2 nM and K11 = 55 nM (Table II).

MantATP Binding to Myo1eIQ and Actomyo1eIQ-- We used the enhancement in fluorescence of the ATP analog mantATP to measure the rate of nucleotide binding. Unlike skeletal muscle myosin-II (32), an intrinsic tryptophan fluorescence change upon MgATP binding was not detected. The transient increase in mantATP fluorescence was best fit to two exponentials (data not shown). The two exponentials are not the result of the presence of multiple mantATP isomers (33) because a two exponential fluorescence increase was also detected when experiments were performed with dmantATP, which is a single-isomer preparation of mantATP. The amplitude of the fast exponential phase is 70-80% of the signal, and its rate is linearly related to the dmantATP concentration (Fig. 4A). We assumed a two-step ATP binding reaction,


View larger version (26K):
[in this window]
[in a new window]
 
Fig. 4.   dmantATP binding to myo1eIQ and actomyo1eIQ. A, rate of dmantATP binding to 0.5 µM myo1eIQ and B, dmantATP binding to 0.5 µM actomyo1eIQ as a function of nucleotide concentration. The observed rates (kobs) were obtained by fitting the stopped flow fluorescence data at each nucleotide concentration to the sum of two exponentials. Apparent second order association rate constants were determined from the slopes of the linear fits of the fast phase of the two exponential fits (solid lines). The insets show the nucleotide dependence of the slow phase of the two exponential fits.


<UP>M</UP>+<UP>ATP</UP> <LIM><OP><ARROW>⇄</ARROW></OP><UL>K<SUB>1</SUB></UL></LIM><UP> M</UP>(<UP>ATP</UP>) <LIM><OP><ARROW>→</ARROW></OP><UL>k<SUB>+2</SUB></UL></LIM> <UP>M·ATP*</UP>

<UP><SC>Reaction</SC> 3</UP>
where K1 is a rapid equilibrium and k+2 is a rate-limiting isomerization to the high fluorescence state of dmantATP (ATP*). The apparent second order rate constant for dmantATP binding, given by the slope of the plot of the observed rate versus nucleotide concentration, is K1k2 = 11 ± 0.4 µM-1 s-1 (Fig. 4A). Fits of the slow exponential phase yielded an apparent second order rate constant 0.65 µM-1 s-1 (Fig. 4A, inset). It was proposed for skeletal muscle subfragment-1 that this slow rate reports the conformational change that accompanies ATP hydrolysis (k3, Scheme 1 (33)).

A two-exponential fluorescence increase was also observed in the presence of actin. The amplitude of the fast exponential phase is 40-60% of the signal, and its rate is linearly related to the nucleotide concentration with an apparent second order rate constant of K1'k2' = 3.7 ± 0.3 µM-1 s-1 (Fig. 4B). The linear fit of the rates does not pass through the origin. The apparent second order rate constant determined from the slow exponential phase is 0.15 µM-1 s-1 (Fig. 4B, inset).

ATP-induced Population of the Weakly Bound States-- Pyrene-actin fluorescence was used to monitor the ATP-induced population of the weakly bound states (13, 34). Mixing ATP with pyrene-actomyo1eIQ resulted in an increase in fluorescence (Fig. 5A, upper trace in inset). Fluorescence transients were best fit to single exponentials at all ATP concentrations examined, and the rate of the fluorescence increase was hyperbolically related to the ATP concentration (Fig. 5A). The mechanism of ATP-induced fluorescence enhancement was modeled as shown in Reaction 4, 


View larger version (28K):
[in this window]
[in a new window]
 
Fig. 5.   MgATP-induced population of weakly bound actomyo1eIQ states. ATP-induced dissociation of 0.25 µM actomyo1eIQ monitored by pyrene-actin fluorescence () and light scattering (black-square). A, the ATP concentration dependence of kobs is shown for the full range of ATP concentrations. All stopped flow transients fit single exponentials. Pyrene fluorescence (solid line) and light scattering (dashed line) data are fit to rectangular hyperbolas. The inset shows pyrene fluorescence (top) and light scattering (bottom) transients obtained by mixing 0.25 µM actomyo1eIQ with 100 µM ATP fit to single exponential functions (smooth lines). B, expanded graph shows that kobs is linearly related to the ATP concentration up to 70 µM ATP and that linear fits of the data pass through the origin.


<AR><R><C><UP>AM</UP>+<UP>ATP</UP> <LIM><OP><ARROW>⇄</ARROW></OP><UL>K<SUB>1</SUB>′</UL></LIM> <UP>AM</UP>(<UP>ATP</UP>) <LIM><OP><ARROW>→</ARROW></OP><UL>k<SUB>+2</SUB>′</UL></LIM></C><C><UP>A*M·ATP
</UP></C></R><R><C></C><C><UP>k<SUB>+9</SUB> ⇵ k<SUB>−9</SUB>
</UP></C></R><R><C></C><C><UP>A</UP>*+<UP>M·ATP</UP></C></R></AR>

<UP><SC>Reaction</SC> 4</UP>
where K1' is a rapid equilibrium, k2' is a rate-limiting isomerization to the high fluorescence A*M·ATP state, and k+9 is the dissociation rate. The maximum rate is k2' = 440 ± 9.0 s-1 and K1' = 397 ± 23 µM (Table II). The association rate constant for MgATP binding to actomyo1eIQ obtained from the initial slope (K1'k2') is 0.86 ± 0.11 µM-1s-1 (Fig. 5B) with a y intercept of zero.

Light scattering measurements were used to monitor the rate of actomyo1eIQ dissociation by MgATP (Fig. 5A). This measurement is different from the pyrene-actin measurements in that pyrene-actin fluorescence reports the population of the attached and detached M·ATP states (Scheme 1). Mixing MgATP with actomyo1eIQ resulted in a rapid decrease in light scattering (Fig. 5A, lower trace in inset). The time courses follow single exponentials with rates that depend hyperbolically on the MgATP concentration (Fig. 5A). Fits to Reaction 4 yielded rate constants (K1' = 334 ± 46 µM; k+2' = 439 ± 23 s-1) nearly identical to those obtained with pyrene-actin (Table II). Therefore, the maximum rate of dissociation (k+9) upon ATP binding is limited by k+2' (Reaction 4).

ATP Hydrolysis-- The rate of ATP hydrolysis in the absence of actin (k+3 + k-3; Scheme 1) was measured directly by quenched flow (Fig. 6). There is a rapid initial burst of ATP hydrolysis which fits a single exponential with a rate of 108 ± 14 s-1 at 50 µM MgATP. The rate of ATP binding under experimental conditions is expected to be >500 s-1 (Table II), so this measurement is not limited by the rate of nucleotide association. The presence of a phosphate burst indicates that the rate-limiting step occurs after ATP hydrolysis. The amplitude of the phosphate burst is 0.55 Pi/myosin (Fig. 6). If we assume that B = K3/(1 + K3) and K3 = (k+3/k-3), then K3 = 1.2, k+3 = 59 s-1, and k-3 = 49 s-1 (Table II).


View larger version (10K):
[in this window]
[in a new window]
 
Fig. 6.   Phosphate burst of myo1eIQ. Time course of ADP·Pi formation by 2.0 µM myo1eIQ in the absence of actin measured by quenched flow after mixing with 50 µM ATP. The smooth line is the best fit of the data to a single exponential (kobs = 108 s-1).

Phosphate Release-- Fluorescently labeled PiBP was used to measure directly the rate of phosphate release (k+4 and k+4') in sequential mix, single turnover, stopped flow experiments (Fig. 7A (30)). Apyrase-treated myo1eIQ (2 µM) was mixed with 1 µM ATP, aged for 400 ms to allow for ATP binding, and mixed with 0-40 µM actin (Fig. 7B). PiBP was included with the myo1eIQ and the actin to prevent transients resulting from phosphate released during the age time or phosphate contamination in the actin.


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 7.   Rate of Pi release from myo1eIQ. A, time course of transient Pi release from myo1eIQ after mixing with ATP in a sequential mix, single turnover experiment in KMg50. The smooth curve is a single exponential fit to the data. The inset is an expanded region of the transient showing the absence of a lag phase. B, actin concentration dependence on the rate of phosphate release measured in KMg0 (black-square) and KMg50 (). Delay time after first mix was 400 ms. Final concentrations at t = 0 were 1.0 µM myo1eIQ, 0.5 µM ATP, 5.0 µM PiBP, 0-40 µM actin.

The time course of phosphate release in the absence of actin follows a single exponential with a rate of 1.5 s-1 in KMg0 and KMg50 (Fig. 7A). There is no lag in the phosphate release transients (Fig. 7A, inset), indicating that phosphate release precedes ADP release. Unlike previously characterized myosins (28, 30), actin did not increase the rate of phosphate release (Fig 7B). This result is likely caused by the weak affinity of the M·ADP·Pi state for actin (see "Discussion").

ADP Association and Dissociation-- The fluorescence increase of mantADP was used to determine the rate constants for ADP association to myo1eIQ and actomyo1eIQ (Fig. 8). In the absence of actin, time courses of fluorescence fit two exponential rates with approximately equal amplitudes. The fast phase was hyperbolically related to the nucleotide concentration (Fig. 8A) and modeled as a two-step binding reaction,


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 8.   Kinetics of mantADP binding to myo1eIQ. A, rate of mantADP binding to 0.5 µM myo1eIQ as a function of nucleotide concentration measured using fluorescence resonance energy transfer. The observed rates (kobs) were obtained by fitting the stopped flow fluorescence transients at each nucleotide concentration to the sum of two exponentials. The solid line is the best fit of the fast kobs () to a rectangular hyperbola. The inset shows the nucleotide concentration dependence of the slow kobs (black-square). B, fluorescence transient showing the time course of mantADP dissociation from myo1eIQ. An equilibrated mixture of 10 µM mantADP and 0.6 µM myo1eIQ was mixed with 2.5 mM ATP. The smooth line shows the best fit of the data to the sum of two exponentials.


<UP>M</UP>+<UP>ADP</UP> <LIM><OP><ARROW>⇄</ARROW></OP><UL>1/K<SUB>6</SUB></UL></LIM> <UP>M</UP>(<UP>ADP</UP>) <LIM><OP><ARROW>→</ARROW></OP><UL>k<SUB>−5</SUB></UL></LIM> <UP>M·ADP∗</UP>

<UP><SC>Reaction</SC> 5</UP>
where * indicates high mantADP fluorescence. The maximum rate is k-5 = 237 ± 13 s-1 and K6 = 6.0 ± 0.9 µM (Table II). The rate of the slow exponential phase shows only a slight nucleotide dependence and plateaus at ~11 s-1 at >10 µM mantADP (Fig. 8A, inset). The origin of this second exponential is not known, but it may report the population of a second M·ADP state (see "Discussion").

In the presence of actin, time courses of fluorescence were best fit to single exponential rates that were linearly related to the nucleotide concentration (Fig. 9A). The signal to noise ratio of the fluorescence signal limited the range of testable mantADP concentrations, so we were not able to measure the maximum rate of mantADP binding (k-5'). Therefore, the apparent second order rate constant for mantADP binding was obtained by a linear fit to the data assuming the two-step binding model of Reaction 5 (k-5'/K6' = 13 ± 1.8 µM-1 s-1; Fig. 9A). The non-zero y intercept reveals a dissociation rate constant of k+5' = 93 ± 5 s-1.


View larger version (22K):
[in this window]
[in a new window]
 
Fig. 9.   Kinetics of mantADP binding to actomyo1eIQ. A, rate of mantADP binding to 0.5 µM actomyo1eIQ as a function of nucleotide concentration. The observed rates (kobs) were obtained by fitting the stopped flow fluorescence data to a single exponential. The solid line is the best linear fit of kobs. The fluorescence transient in the inset shows the dissociation of mantADP from actomyo1eIQ obtained by mixing 2.5 mM ATP with an equilibrated mixture of 10 µM mantADP and 0.5 µM actomyo1eIQ. The smooth line is a single exponential fit to the data. B, dependence on the rate of 80 µM ATP binding to pyrene-actomyo1eIQ on the concentration of ADP. The observed rate (kobs) at each ADP concentration was obtained by fitting the fluorescent transient to a single exponential. The solid line is the best fit of the data to a hyperbola as described in Equation 1.

The rates of mantADP dissociation from myo1eIQ (Fig. 8B) and actomyo1eIQ (Fig. 9A, inset) were measured by competition with 250-fold excess of MgATP. Fluorescence time courses of mantADP dissociation fit two-exponential rates of 6.3 and 1.7 s-1 in the absence of actin and a single exponential rate of k+5' = 104 ± 4 s-1 in the presence of actin (Table II). The rate of MgADP dissociation from actomyo1eIQ is similar to that reported by the y intercept in the mantADP association experiment (Fig. 9A).

The affinity of MgADP for actomyo1eIQ was determined by competition experiments in which the ternary pyrene-actomyo1eIQ-ADP complex was mixed with MgATP (Fig. 9B). Fluorescence transients fit single exponential rates that had a hyperbolic dependence on the MgADP concentration, suggesting that AM·D and AM states are in rapid equilibrium. Therefore, the dependence of kobs on MgADP concentration is as follows,
k<SUB><UP>obs</UP></SUB>=<FR><NU>k<SUB>0</SUB></NU><DE>1+<FENCE><FR><NU>[<UP>ADP</UP>]</NU><DE>K<SUB>d</SUB></DE></FR></FENCE></DE></FR> (Eq. 1)
where kobs is the observed dissociation rate, k0 is the dissociation rate in the absence of ADP (K1'k2'[ATP]; Table II), and Kd is the dissociation equilibrium constant for ADP (K5'K6'). Nonlinear least square fits to the data yielded a dissociation equilibrium constant of K5'K6' = 3.9 ± 0.29 µM (Table II), which is similar to the calculated value (K5'K6' = 8 µM) using association and dissociation rate of mantADP.

    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Overview of the Myo1eIQ ATPase-- Myo1eIQ has a high basal ATPase rate (vo; Table I), yet all measured rate constant are >2-fold faster than v0 (Table II), so we cannot assign a single step as rate-limiting. The slower turnover number is the result of two sequential steps (phosphate release and ADP release) with nearly equivalent rates. Kinetic modeling of the ATPase reaction using our determined rate constants (Table II) yields an ATPase rate (v0 = 0.65 s-1 in KMg50), which is similar to the experimentally determined value (v0 = 0.58 ± 0.04 s-1; Table I). The predominant steady-state intermediates in the absence of actin are M·ADP·Pi and M·ADP (Scheme 1).

Overview of the Actomyo1eIQ ATPase-- Myo1eIQ has a low KATPase and has a ~2-fold actin activation of its ATPase rate (Table I), which is similar to that found for the native molecule. Actin (0-40 µM) does not increase the rate of Pi release in KMg0 or KMg50 (Fig. 7B). Therefore, actin activation of Pi release is not the mechanism for actin activation of the ATPase rate, as it is for other myosins. Rather, the increased ATPase rate is the result of actin activation of ADP release. The predominant pathway under actin concentrations used in our assays is one in which (a) ATP binding, actin dissociation, and ATP hydrolysis are fast and not rate-limiting; (b) Pi is released while myosin is detached from actin; (c) M·ADP binds to actin; and (d) ADP is released from actin-bound myosin (Fig. 10A, outlined pathway). The rate-limiting step at saturating actin (Vmax) is defined by the rate of Pi release (k+4), because at saturating actin ADP release is ~100 s-1 (Table II). The KATPase value is related to the affinity of the M·ADP state for actin. Because binding measurements indicate that myo1eIQ binds weakly to actin in the presence of ATP (Fig. 2), the predominant steady-state intermediate is an actin-detached M·ADP·Pi state.


View larger version (29K):
[in this window]
[in a new window]
 
Fig. 10.   Models for myo1e function. A, the predominant ATPase pathway of myo1eIQ under in vitro assay conditions is outlined. We propose that under cellular conditions, native myo1e follows the pathway proposed by Lymn and Taylor (43; gray shading). B, comparison of the loop-2 regions of myosins isoforms. Basic amino acids are boxed, and the conserved KK regions shown to be necessary for actin activation of Pi are underlined. Numbers identify the sequence positions. SmII, smooth muscle myosin-II; SkII, skeletal muscle myosin-II; DdII, Dictyostelium myosin-II; AcMIC, Acanthamoeba myosin-IC.

ATP Binding and Actomyo1eIQ Dissociation-- The binding of ATP to myo1eIQ and actomyo1eIQ, as determined by dmantATP fluorescence, is fast. The rate of binding at physiological nucleotide concentrations is >100-fold faster than the rate-limiting step(s). Therefore, the nucleotide-free state of myo1eIQ is not significantly populated. The rate of population of the weak binding states as measured by pyrene-actin fluorescence and light scattering is also fast, and ATP-induced dissociation from actin is limited only by the rate of ATP binding.

The rate of ATP binding to actomyo1eIQ as measured by mant fluorescence is ~4-fold faster than the rate measured by pyrene-actin fluorescence or light scattering. Additionally, the linear fit of the rates of dmantATP binding to actomyo1eIQ does not pass through the origin (Fig. 4B) but gives a y intercept of 25 s-1. The value of the y intercept likely reports a reverse isomerization rate (k-2') of myo1eIQ. However, because the rate of ATP-induced dissociation of actomyo1eIQ is >400 s-1 (Fig. 5) and essentially irreversible under the low actin concentration conditions of the experiment, a relatively slow reverse isomerization of the ATP binding step should not be observable. Therefore, the faster rate of mantATP binding and the observed reversal rate (k-2') are likely the result of an artifact because of the fluorescent modification of the 3' position of ATP. Although the dmantATP experiments probably do not report the true rate constants for ATP binding in the presence of actin, the experiments do reveal differences in the structural isomerization required for tight ATP binding in the presence and absence of actin. Additionally, because similar mant artifacts are not seen with other characterized myosins (e.g. 14), the experiments show that there are significant structural differences at the nucleotide binding sites of myosin-I isoforms.

Phosphate Release-- The v0 and Vmax rates and the extent of actin activation of myo1eIQ (Table I) are similar to the values measured for native myo1e purified from rat liver (19), so the biochemical rate and equilibrium constants determined for Sf9-expressed myo1eIQ are likely identical to the native molecule. Therefore, the fast rate of Pi release (k+4) is not an artifact of the expression system or of the expression of a truncated protein.

The rate of phosphate release from all previously characterized myosins is increased by actin binding (28, 30, 35). However, we did not detect actin activation of phosphate release from myo1eIQ. We propose that the lack of actin activation is caused by the low affinity of the M·ADP·Pi state for actin. Other myosins bind actin in the presence of ATP with equilibrium dissociation constants <50 µM at low ionic strength conditions (36, 37), whereas myo1eIQ clearly binds with a dissociation constant > 50 µM (Fig. 2). The ionic component of actin binding has been shown to be mediated by positive charges in surface loop-2 of myosin (36). Loop-2 of myo1e is shorter than those found in other characterized myosins (Fig. 10B) and does not contain the positively charged amino acids in the region that has been shown to be required for actin binding (Fig. 10B (36)). However, loop-2 of myo1e contains the KK region (Fig. 10B, underlined) shown to be crucial for actin activation of phosphate release (37). Therefore, it is possible that actin activates phosphate release from myo1eIQ, but we are not able to achieve high enough actin concentrations to populate the weakly bound AM·ADP·Pi state. Weak actin binding may be a kinetic adaptation required for cellular function (see below).

It is possible that phosphate release by myo1eIQ is uncoupled from actin binding, so the very high basal rate of Pi release (k+4) is also the rate when myo1eIQ is attached to actin (k+4'). Why myo1e would evolve a mechanism that seemingly wastes ATP is not clear; however, other mechanisms for the suppression of the basal ATPase rate may exist. For example, the tail domain (which is not present in myo1eIQ) is thought to regulate the ATPase activity of myo1e allosterically (19). Further experiments are required to determine whether the myo1e tail regulates phosphate release directly.

ADP Binding and Release-- Time courses of mantADP binding to myo1eIQ are best fit to two exponentials. The rate of the fast phase is hyperbolically related to the nucleotide concentration as expected for a two-step binding reaction (Reaction 6). The rate of the slow exponential does not show a significant nucleotide dependence and possibly reports the population of a second higher fluorescence state M·ADP state,


<UP>M</UP>+<UP>ADP</UP> <LIM><OP><ARROW>⇄</ARROW></OP><UL>1/K<SUB>6</SUB></UL></LIM> <UP>M</UP>(<UP>ADP</UP>) <LIM><OP><ARROW>→</ARROW></OP><UL>k<SUB>−5a</SUB></UL></LIM> <UP>M·ADP</UP>* <LIM><OP><ARROW>→</ARROW></OP><UL>k<SUB>−5<UP>b</UP></SUB></UL></LIM><UP> M·ADP**</UP>

<UP><SC>Reaction</SC> 6</UP>
where * indicates the fluorescence state of myosin, k-5a is the maximum rate of the fast phase, and k-5b is the maximum rate of the slow fluorescence change. It has been shown previously that myosins exists in more than one ADP state (14, 38), so it is not surprising that we are able to detect multiple conformations of myo1eIQ. Dissociation of mantADP in the absence of actin also shows two exponentials that are likely related to these two states. However, a more detailed analysis of ADP binding is required to characterize these states better.

The similarity of the rates of Pi release and ADP release in the absence of actin raises the possibility that the products are released simultaneously (Table II). However, if ADP and Pi were released together, steady-state measurements of v0 would be >2-fold faster (see above).

Time courses of mantADP binding to actomyo1eIQ fit single exponentials that are linearly related to the nucleotide concentration (Fig. 9A). A linear fit to the data does not pass through the origin but reveals a dissociation rate (k+5') of ~100 s-1. This dissociation rate is nearly identical to the dissociation rate determined directly by displacing myo1eIQ-bound mantADP with unlabeled ATP (Fig. 9A, inset). The dissociation equilibrium constant (K5'K6') calculated by dividing the dissociation rate (k+5') by the association rate (k-5'/K6') obtained from the mantADP experiments is within a factor of 2 of that determined from the pyrene-actin fluorescence experiments. Therefore, it is unlikely that the mantADP experiments are reporting artifactually high rates as shown for dmantATP binding to actin (see above).

Comparison with Other Myosin-I-- Myo1e is the first vertebrate subclass-1 myosin-I to be characterized. Like the subclass-1 isoforms from Acanthamoeba (13), myo1e has relatively large rate constants that are similar to those of skeletal muscle myosin-II (Table II). Most notably, the rates of ADP release from all subclass-1 isoforms are >10-fold faster than subclass-2 isoforms, i.e. the lifetimes of the AM·ADP states are >10 longer for subclass-2 isoforms (13-16). The rate of ADP release limits sliding velocity (39), so we propose that subclass-1 isoforms are better tuned for fast motility, whereas subclass-2 isoforms are better tuned for maintenance of force (14, 15).

A difference among the amoeba and vertebrate subclass-1 isoforms is that Acanthamoeba myosin-IA and -IB have low basal ATPase rates (<0.1 s-1) and a large actin activation of the steady-state ATPase (> 50-fold (40)), whereas myo1e has an high basal ATPase rate and a low actin activation (Table I). Like myo1e, vertebrate subclass-2 isoforms also have high basal ATPase rates and relatively low actin activations (41). This might be an adaptation to the cellular environment of vertebrate cells, or it might be related to the regulation of myosin-I in these cells (8). A second difference between amoeba subclass-1 isoforms and myo1e is, like subclass-2 isoforms, myo1e has a very high ADP affinity (K5'K6' < 10 µM). Therefore, it is possible that ADP release from both vertebrate subclasses is highly force-dependent (14, 15).

Relevance of Myo1e Kinetics to in Vivo Function-- Like all other characterized myosin-I isoforms (13-16), myo1e is a low duty ratio motor, so under unloaded conditions it is predominantly weakly bound or detached from actin filaments. Therefore, for myole to support motility, a high effective duty ratio must be created by bringing together locally high concentrations of myosin and actin. Such a mechanism is consistent with its observed cellular localization in regions of high F-actin concentration (17, 18).

The force-generating power stroke of myosin is thought to accompany phosphate release (42). Therefore, an ATPase pathway in which phosphate is released while myosin is detached from the actin filament (Fig. 10A, outlined pathway) would not generate force. It is possible that the power stroke occurs at a different step in the myo1e ATPase pathway (e.g. ADP release); however, this possibility is unlikely given the sequence and mechanistic conservation among myosin isoforms. We propose that the low actin affinity of myo1e in the presence of ATP is an adaptation to the high actin concentration environment in which myo1e functions (17, 18). Native myo1e cross-links actin filaments via an ATP-insensitive actin binding site in its tail, creating densely packed and cross-linked bundles (19). If the weak binding states (M·ATP and M·ADP·Pi) of myo1e had a high actin affinity, then ATP hydrolysis and the recovery stroke would occur while attached to actin (30). Therefore, an associated hydrolysis pathway would result in the reversal of myosin's power stroke. To ensure that myo1e is detached from actin during the recovery stroke, a weak affinity must be maintained in the pre-force-generating states. We therefore predict that under the high actin concentrations of the cell, the M·ADP·Pi state will bind and release phosphate while attached to actin (Fig. 10A, gray pathway) as proposed for other myosins (43). Further experiments are required to confirm this prediction and to understand better the unique kinetic adaptations of all vertebrate myosin-I isoforms.

    ACKNOWLEDGEMENTS

We thank Tianming Lin and Jade Poole for technical assistance; Li-Qiong Chen for constructing the expression virus; William M. Bement for providing the myo1e cDNA; and H. Lee Sweeney, Yale E. Goldman, and Enrique M. De La Cruz for helpful discussions and comments on the manuscript.

    FOOTNOTES

* This work was supported in part by National Institutes of Health Grant GM57247 and an American Heart Association Grant (to E. M. O.).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.

§ Supported by a Wellcome Trust Research career development fellowship.

|| To whom correspondence should be addressed: Dept. of Physiology, B400 Richards, University of Pennsylvania School of Medicine, 3700 Hamilton Walk, Philadelphia, PA 19104-6085. Tel.: 215-573-9758; Fax: 215-573-1171; E-mail: ostap@mail.med.upenn.edu.

Published, JBC Papers in Press, April 8, 2002, DOI 10.1074/jbc.M200713200

    ABBREVIATIONS

The abbreviations used are: mant, N-methyanthraniloyl; A, actin; CaM, calmodulin; dmant, 2'-deoxymant; DTT, dithiothreitol; M, myosin; PiBP, phosphate-binding protein; pyrene-actin, pyrenyl iodoacetamide-labeled actin.

    REFERENCES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

1. Berg, J. S., Powell, B. C., and Cheney, R. E. (2001) Mol. Biol. Cell 12, 780-794[Abstract/Free Full Text]
2. Dai, J., Ting-Beall, H. P., Hochmuth, R. M., Sheetz, M. P., and Titus, M. A. (1999) Biophys. J. 77, 1168-1176[Abstract/Free Full Text]
3. Geli, M. I., and Riezman, H. (1996) Science 272, 533-535[Abstract]
4. Jung, G., Wu, X., and Hammer, J. A., III (1996) J. Cell Biol. 133, 305-323[Abstract/Free Full Text]
5. Raposo, G., Cordonnier, M. N., Tenza, D., Menichi, B., Durrbach, A., Louvard, D., and Coudrier, E. (1999) Mol. Biol. Cell. 10, 1477-1494[Abstract/Free Full Text]
6. Novak, K. D., Peterson, M. D., Reedy, M. C., and Titus, M. A. (1995) J. Cell Biol. 131, 1205-1221[Abstract/Free Full Text]
7. Gillespie, P. G., Wagner, M. C., and Hudspeth, A. J. (1993) Neuron 11, 581-594[CrossRef][Medline] [Order article via Infotrieve]
8. Tang, N., and Ostap, E. M. (2001) Curr. Biol. 11, 1131-1135[CrossRef][Medline] [Order article via Infotrieve]
9. Sokac, A. M., and Bement, W. M. (2000) Int. Rev. Cytol. 200, 197-304[CrossRef][Medline] [Order article via Infotrieve]
10. Pollard, T. D., and Korn, E. D. (1973) J. Biol. Chem. 248, 4682-4690[Abstract/Free Full Text]
11. Bement, W. M., Wirth, J. A., and Mooseker, M. S. (1994) J. Mol. Biol. 243, 356-363[CrossRef][Medline] [Order article via Infotrieve]
12. Stöffler, H.-E., Ruppert, C., Reinhard, J., and Bähler, M. (1995) J. Cell Biol. 129, 819-830[Abstract/Free Full Text]
13. Ostap, E. M., and Pollard, T. D. (1996) J. Cell Biol. 32, 1053-1060
14. Jontes, J. D., Milligan, R. A., Pollard, T. D., and Ostap, E. M. (1997) Proc. Natl. Acad. Sci. U. S. A. 94, 14332-14337[Abstract/Free Full Text]
15. Coluccio, L. M., and Geeves, M. A. (1999) J. Biol. Chem. 274, 21575-21580[Abstract/Free Full Text]
16. Geeves, M. A., Perreault-Micale, C., and Coluccio, L. M. (2000) J. Biol. Chem. 275, 21624-21630[Abstract/Free Full Text]
17. Stöffler, H.-E., Honnert, U., Bauer, C. A., Hofer, D., Schwarz, H., Muller, R. T., Drenckhahn, D., and Bähler, M. (1998) J. Cell Sci. 111, 2779-2788[Abstract]
18. Skowron, J. F., Bement, W. M., and Mooseker, M. S. (1998) Cell Motil. Cytoskeleton 41, 308-324[CrossRef][Medline] [Order article via Infotrieve]
19. Stöffler, H.-E., and Bähler, M. (1998) J. Biol. Chem. 273, 14605-14611[Abstract/Free Full Text]
20. Bement, W. M., and Mooseker, M. S. (1995) Cell Motil. Cytoskeleton 31, 87-92[CrossRef][Medline] [Order article via Infotrieve]
21. Hiratsuka, T. (1983) Biochim. Biophys. Acta 742, 496-508[CrossRef][Medline] [Order article via Infotrieve]
22. Spudich, J. A., and Watt, S. (1971) J. Biol. Chem. 246, 4866-4871[Abstract/Free Full Text]
23. Houk, T. W., and Ue, K. (1974) Anal. Biochem. 62, 66-74[CrossRef][Medline] [Order article via Infotrieve]
24. Pollard, T. D. (1984) J. Cell Biol. 99, 769-777[Abstract/Free Full Text]
25. Putkey, J. A., Slaughter, G. R., and Means, A. R. (1985) J. Biol. Chem. 260, 4704-4712[Abstract/Free Full Text]
26. De La Cruz, E. M., Wells, A. L., Rosenfeld, S. S., Ostap, E. M., and Sweeney, H. L. (1999) Proc. Natl. Acad. Sci. U. S. A. 96, 13726-13731[Abstract/Free Full Text]
27. De La Cruz, E. M., Sweeney, H. L., and Ostap, E. M. (2000) Biophys. J. 79, 1524-1529[Abstract/Free Full Text]
28. De La Cruz, E. M., Wells, A. L., Sweeney, H. L., and Ostap, E. M. (2000) Biochemistry 39, 14196-14202[CrossRef][Medline] [Order article via Infotrieve]
29. Chalovich, J. M., and Eisenberg, E. (1982) J. Biol. Chem. 257, 2432-2437[Free Full Text]
30. White, H. D., Belknap, B., and Webb, M. R. (1997) Biochemistry 30, 11828-11836
31. Xing, J., Wriggers, W., Jefferson, G. M., Stein, R., Cheung, H. C., and Rosenfeld, S. S. (2000) J. Biol. Chem. 275, 35413-35423[Abstract/Free Full Text]
32. Johnson, K. A., and Taylor, E. W. (1978) Biochemistry 17, 3432-3442[CrossRef][Medline] [Order article via Infotrieve]
33. Woodward, S. K. A., Eccleston, J. F., and Geeves, M. A. (1991) Biochemistry 30, 422-430[CrossRef][Medline]