| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |

* Department of Bioengineering, Rice University, Houston, Texas 77005; and
Graduate Program of Structural and Computational Biology and Molecular Biophysics and
Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, Texas 77030
Correspondence: Address reprint requests to Jianpeng Ma, One Baylor Plaza, BCM-125, Baylor College of Medicine, Houston, TX 77030. Tel.: 713-798-8187; Fax: 713-796-9438; E-mail: jpma{at}bcm.tmc.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
200 Å in diameter and up to 10 µm (100,000 Å) in length, which serves as a propeller (Berry and Armitage, 1999
|
11 monomers for every two turns of the one-start helix. Alternately, the filament can be depicted as 11 protofilaments arranged around a hollow central channel
30 Å in diameter (Samatey et al., 2001
51 kDa protein with 494 residues, consisting of four domains, D0, D1, D2, and D3, named in accordance with their radial positions from the inside of the filament to the outside. D0, D1, and D2 comprise the core of the filament, whereas D3 protrudes from the filament's surface (Mimori et al. 1995
|
The differential mechanical properties, such as Young's modulus E and shear modulus G, of the hook and filament are essential to the function of the flagellum and are directly linked to structural differences between the two. They have been examined using a variety of experimental methods, ranging from quasielastic scattering of laser light to optical tweezers (Table 1) (Berg, 1976
; Block et al., 1989
, 1991; Fujime et al., 1972
; Hoshikawa, 1983
; Hoshikawa and Kamiya, 1985
; Powers, 2002
; Trachtenberg and Hammel, 1992
). The measured Young's modulus E and shear modulus G of the filament are in the range of 10101012 and 10111012 dyn/cm2, respectively. The shear modulus G of the hook was also determined to be 108 dyn/cm2. These values suggested the filament to be orders of magnitude stiffer than the hook. Specifically, Block and colleagues (Block et al., 1989
, 1991
) determined that the torsional/twisting compliance of the flagellum resides mainly in the hook, leading them to speculate that this compliance could allow multiple filaments to bundle in parallel and rotate in phase (Fig. 1 C). The differential moduli of the filament and hook provide an ingenious mechanism that is the key to unified motion in these bacteria. This motion would not be possible if the entire flagellum were composed of only one protein with an identical set of mechanical properties.
|
1.3 (Vogel, 1992
In this article we employ a new computational method, the quantized elastic deformational model (QEDM) (Ming et al., 2002a
; Tama et al., 2002
), to calculate the dimensionless twist/bend ratio of both the filament and hook, thereby providing a quantitative means to compare their stiffness. QEDM is able to extract motional and mechanical data from low-resolution protein structures, such as electron density maps. The protein is modeled as an elastic object and the mass density distribution is equated to the electron density distribution of the protein. Here, QEDM is used to determine the vibrational normal modes of
1000-Å-long electron density map sections of S. typhimurium filament and hook. The eigenvalues of the first bending and twisting modes of a structure are directly related to the flexural and torsional stiffness, respectively (see Methods) (Meirovitch, 2001
).
| RESULTS |
|---|
|
|
|---|
1000 Å in length at resolutions of 10 Å and 9 Å, respectively, were discretized into 8000 Voronoi cells. QEDM analysis was then performed for the hook and filament at several cutoff distances, rc, between 35 and 55 Å (see Table 2 and Methods). These cutoffs are larger than the 13-Å cutoff that is typically employed in QEDM, due to the enormous size of the filament and hook, and the coarse-grained nature of the model (Beuron et al., 2003
|
|
|
and
where R and r represent the outer and inner radii, respectively (Roark, 1989
In contrast, an experimentally determined hook twist/bend ratio could not be obtained because Young's modulus E of the hook has never been directly determined in experiments. However, Block and colleagues did obtain an experimental value of the hook shear modulus G, 108 dyn/cm2 (Block et al., 1991
). If the maximum and minimum values of the QEDM-computed hook twist/bend ratios from Table 2 are employed in conjunction with Block's value of the shear modulus G above, it is trivial to calculate a theoretical Young's modulus E for the hook, which is found to be on the order of 106107 dyn/cm2. This calculated value is several orders of magnitude smaller than the lowest experimentally determined Young's modulus for the filament (Table 1).
Radial distributions of the hook and filament
An important issue that must be addressed is the rate of convergence of the QEDM-calculated twist/bend ratios for the hook and filament. It appears that the ratios for the filament converged slightly quicker over the range of QEDM cutoff distances than for the hook (Fig. 3). This differential convergence rate of the ratios of the two structures is believed to be the direct result of their inherently different density distributions. The radial distribution function, g(r), for each structure was calculated and plotted in Fig. 4. QEDM relies on a set of finite Voronoi cells to represent the original density map (see Methods); thus the radial distributions computed are not exactly identical to those of the initial density maps (see Fig. 4 in Mimori et al., 1995
for the filament), but the general features of the curves are similar. A comparison of the radial distribution curves of the hook and filament revealed that the density for the filament is evenly dispersed over two peaks (domains D0, D1, and D2) lying in a radial range of 1050 Å, whereas that of the hook is clustered primarily around a radial distance of 55 Å (Fig. 4). The net effect of these differences is that the filament has a more densely packed, rigid core than the hook, thus imparting a higher degree of stiffness to the filament. The more evenly distributed density of the core region of the filament also accounts for the quicker convergence of the twist/bend ratios for the filament over the range of QEDM cutoff distances as seen in Fig. 3.
| CONCLUDING DISCUSSION |
|---|
|
|
|---|
In conclusion, the theoretical nature of Young's modulus calculated for the hook in this study should be stressed. Although the QEDM-obtained hook twist/bend ratios did not converge as quickly as they did for the filament, it is likely due to the inherent difference of density distribution of the two structures, rather than a simulation artifact. The theoretical value of Young's modulus should serve as a guide for further experimental studies.
| METHODS |
|---|
|
|
|---|
atoms of the backbone. ANM is essentially a normal mode analysis (Brooks et al., 1995
![]() | (1) |
represents the strength of the potential, and is assumed to be the equivalent for all pairs, rc is the cutoff distance, and
and
are the instantaneous and equilibrium values of the pairwise distance between the ith and jth C
atoms, respectively (Tirion, 1996
ij is the Heaviside step function, which details the cutoff effect of the interaction (Hildebrand, 1976
QEDM was used to analyze both the hook and filament. A 10-Å resolution electron density map of a 960-Å-long section of S. typhimurium hook was provided by D. Thomas, T. Shaikh, and D. DeRosier, private communications). The 9-Å resolution map of a 303-Å-long section of filament, also from S. typhimurium, was obtained from Koji Yonekura and Keiichi Namba (Mimori et al., 1995
). A 1000-Å section of filament was constructed from the 303-Å section using the programs EMAN (version 1.2; ncmi.bcm.tmc.edu) and IRIS Explorer (release 5.0; www.nag.co.uk). In both cases, 8000 Voronoi cells were used to discretize the density. The number of Voronoi cells can be set to either correspond with the number of C
atoms in the structure or be much less than this number, yielding a coarse-grained model, as is the case in this work (Doruker et al., 2002
; Ming et al., 2002a
,b
). Additional calculations were performed on densities discretized into 4000 Voronoi cells (data not shown). The resulting modes were identical to those calculated for the 8000 Voronoi cell case, indicating that resolution was not a factor.
was set to 1.0 for all of the simulations because its value does not affect the distribution of the motions (Atilgan et al., 2001
; Ming et al., 2002a
,b
). Solving for the eigenvalues and eigenvectors in each simulation required
100 h of computing time on a Silicon Graphics Iris Indigo (Mountain View, CA).
Relating frequency of vibration to flexural and torsional rigidity
The equation for vibration of a beam, of length L, in bending is given by
![]() | (2) |
is the vibration frequency, and m(x) is the mass per unit length (Meirovitch, 2001
![]() | (3) |
The solution of this fourth order differential equation is
![]() | (4) |
(no bending moment at the ends of the beam) and
(no shearing force at the ends of the beam), yielding the first relationship of interest
![]() | (5) |
for the first bending mode.
Similarly, the equation for torsion in a shaft, of length L, is given by
![]() | (6) |
(x) is the twist angle, and
is the vibration frequency (Meirovitch, 2001
![]() | (7) |
The solution to this second order differential equation is
![]() | (8) |
(no torque at the ends of the shaft). Leading to the second relationship of interest,
![]() | (9) |
for the first torsion mode. For a solid rod
where m and L are the total mass and length of the rod, respectively (Feynman et al., 1964
![]() | (10) |
The last equality that must be considered is the relationship between frequencies and eigenvalues. QEDM outputs a set of displacement coordinates and eigenvalues,
, for each calculated mode, 1N. From ANM,
(Atilgan et al., 2001
), where
i is the eigenvalue of the ith mode,
is the strength of the potential,
i is the frequency of the ith mode, and mp is the point mass. Total mass, m, in Eqs. 5 and 10 is equal to the product of point mass, mp, times the number of points, np (i.e., Voronoi cells for QEDM). Substituting these quantities into Eqs. 5 and 10 yields
![]() | (11) |
Dimensionless twist/bend ratio
Taking the ratio of EI/GJ from Eq. 11 generates the desired relationship
![]() | (12) |
i and
j represent the ith and jth eigenvalues for the first bending (mode 7) and torsional (mode 9) modes, respectively. | ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This research was supported in part by grants to J.M. from the American Heart Association (AHA-TX0160107Y), the Robert A. Welch Foundation (Q-1512), the National Institutes of Health (R01-GM067801), and a National Science Foundation Career Award (MCB-0237796). J.M. is a recipient of the Award for Distinguished Young Scholars Abroad from the National Natural Science Foundation of China.
Submitted on December 16, 2003; accepted for publication January 28, 2004.
| REFERENCES |
|---|
|
|
|---|
Atilgan, A. R., S. R. Durell, R. L. Jernigan, M. C. Demirel, O. Keskin, and I. Bahar. 2001. Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys. J. 80:505515.
Berg, H. C. 1976. Does the flagellar rotary motor step? In Cell Motility. R. Goldman, T. Pollard, and J. Rosenbaum, editors. Cold Spring Harbor Laboratory, New York. 4756.
Berg, H. C., and R. A. Anderson. 1973. Bacteria swim by rotating their flagellar filaments. Nature. 245:380382.[Medline]
Berry, R. M., and J. P. Armitage. 1999. The bacterial flagella motor. Adv. Microb. Physiol. 41:291337.[Medline]
Beuron, F., T. C. Flynn, J. Ma, H. Kondo, X. Zhang, and P. S. Freemont. 2003. Motions and negative cooperativity between p97 domains revealed by cryo-electron microscopy and quantized elastic deformational model. J. Mol. Biol. 327:619629.[Medline]
Block, S. M., D. F. Blair, and H. C. Berg. 1989. Compliance of bacterial flagella measured with optical tweezers. Nature. 338:514518.[Medline]
Block, S. M., D. F. Blair, and H. C. Berg. 1991. Compliance of bacterial polyhooks measured with optical tweezers. Cytometry. 12:492496.[Medline]
Boyer, P. D. 1997. The ATP synthase: a splendid molecular machine. Annu. Rev. Biochem. 66:717749.[Medline]
Brooks, B. R., D. Janezic, and M. Karplus. 1995. Harmonic analysis of large systems. I. Methodology. J. Comput. Chem. 16:15221542.
Doruker, P., R. L. Jernigan, and I. Bahar. 2002. Dynamics of large proteins through hierarchical levels of coarse-grained structures. J. Comput. Chem. 23:119127.[Medline]
Etnier, S. 2001. Flexural and torsional stiffness in multi-jointed biological beams. Biol. Bull. 200:18.
Feynman, R. P., R. B. Leighton, and M. Sands. 1964. The Feynman Lectures on Physics. Addison-Wesley, Reading, MA.
Fujime, S., M. Maruyama, and S. Asakura. 1972. Flexural rigidity of bacterial flagella studied by quasielastic scattering of laser light. J. Mol. Biol. 68:347359.[Medline]
Glagolev, A. N., and V. P. Skulachev. 1978. The proton pump is a molecular engine of motile bacteria. Nature. 272:280282.[Medline]
Gray, R. M. 1984. Vector quantization. IEEE ASSP Mag. 1:429.
Hildebrand, F. B. 1976. Advanced Calculus for Applications. Prentice-Hall, Inc., Englewood Cliffs, NJ.
Hoshikawa, H. 1983. Elastic properties of bacterial flagellar filaments. I. Free rotation case. Biophys. Chem. 17:105109.[Medline]
Hoshikawa, H., and R. Kamiya. 1985. Elastic properties of bacterial flagellar filaments. II. Determination of the modulus of rigidity. Biophys. Chem. 22:159166.[Medline]
Kamiya, R., S. Asakura, K. Wakabayashi, and K. Namba. 1979. Transition of bacterial flagella from helical to straight forms with different subunit arrangements. J. Mol. Biol. 131:725742.[Medline]
Kong, Y., D. Ming, Y. Wu, J. K. Stoops, Z. H. Zhou, and J. Ma. 2003. Conformational flexibility of pyruvate dehydrogenase complexes: a computational analysis by quantized elastic deformational model. J. Mol. Biol. 330:129135.[Medline]
Larsen, S. H., J. Adler, J. J. Gargus, and R. W. Hogg. 1974. Chemomechanical coupling without ATP: the source of energy for motility and chemotaxis in bacteria. Proc. Natl. Acad. Sci. USA. 71:12391243.
Ma, J. 2004. New advances in normal mode analysis of supermolecular complexes and applications to structural refinement. Curr. Protein Pept. Sci. 5:119123.[Medline]
Ma, J., T. C. Flynn, Q. Cui, A. G. Leslie, J. E. Walker, and M. Karplus. 2002. A dynamic analysis of the rotation mechanism for conformational change in F1-ATPase. Struct. 10:921931.
Macnab, R. M., and M. K. Ornston. 1977. Normal-to-curly flagellar transitions and their role in bacterial tumbling. Stabilization of an alternative quaternary structure by mechanical force. J. Mol. Biol. 112:130.[Medline]
Makhoul, J., S. Roucos, and H. Gish. 1985. Vector quantization in speech coding. Proc. IEEE 73:15511588.
Manson, M. D., P. Tedesco, H. C. Berg, F. M. Harold, and C. Van der Drift. 1977. A protonmotive force drives bacterial flagella. Proc. Natl. Acad. Sci. USA. 74:30603064.
Meirovitch, L. 2001. Fundamentals of Vibrations. McGraw-Hill, New York.
Mimori, Y., I. Yamashita, K. Murata, Y. Fujiyoshi, K. Yonekura, C. Toyoshima, and K. Namba. 1995. The structure of the R-type straight flagellar filament of Salmonella at 9 Å resolution by electron cryomicroscopy. J. Mol. Biol. 249:6987.[Medline]
Ming, D., Y. Kong, M. A. Lambert, Z. Huang, and J. Ma. 2002a. How to describe protein motion without amino acid sequence and atomic coordinates. Proc. Natl. Acad. Sci. USA. 99:86208625.
Ming, D., Y. Kong, S. J. Wakil, J. Brink, and J. Ma. 2002b. Domain movements in human fatty acid synthase by quantized elastic deformational model. Proc. Natl. Acad. Sci. USA. 99:78957899.
Morgan, D. G., R. M. Macnab, N. R. Francis, and D. J. DeRosier. 1993. Domain organization of the subunit of the Salmonella typhimurium flagellar hook. J. Mol. Biol. 229:7984.[Medline]
Morgan, D. G., C. Owen, L. A. Melanson, and D. J. DeRosier. 1995. Structure of bacterial flagellar filaments at 11 Å resolution: packing of the alpha-helices. J. Mol. Biol. 249:88110.[Medline]
Niklas, K. J. 1991. The elastic moduli and mechanics of Populus tremuloides (Salicaceae) petioles in bending and torsion. Am. J. Bot. 78:989996.
O'Brien, E. J., and P. M. Bennett. 1972. Structure of straight flagella from a mutant Salmonella. J. Mol. Biol. 70:133152.[Medline]
Powers, T. R. 2002. Role of body rotation in bacterial flagellar bundling. Phys. Rev. E. 65:040903.
Roark, R. J. 1989. Formulas for Stress and Strain. McGraw-Hill, Inc., New York.
Samatey, F. A., K. Imada, S. Nagashima, F. Vonderviszt, T. Kumasaka, M. Yamamoto, and K. Namba. 2001. Structure of the bacterial flagellar protofilament and implications for a switch for supercoiling. Nature. 410:331337.[Medline]
Silverman, M., and M. Simon. 1974. Flagellar rotation and the mechanism of bacterial motility. Nature. 249:7374.[Medline]
Tama, F., W. Wriggers, and C. L. Brooks,, III. 2002. Exploring global distortions of biological macromolecules and assemblies from low-resolution structural information and elastic network theory. J. Mol. Biol. 321:297305.[Medline]
Tirion, M. M. 1996. Large amplitude elastic motions in proteins from a single-parameter, atomic analysis. Phys. Rev. Lett. 77:19051908.[Medline]
Trachtenberg, S., and I. Hammel. 1992. The rigidity of bacterial flagellar filaments and its relation to filament polymorphism. J. Struct. Biol. 109:1827.[Medline]
Vogel, S. 1992. Twist-to-bend ratios and cross-sectional shapes of petioles and stems. J. Exp. Bot. 43:15271532.
Yonekura, K., S. Maki-Yonekura, and K. Namba. 2003. Complete atomic model of the bacterial flagellar filament by electron cryomicroscopy. Nature. 424:643650.[Medline]
This article has been cited by other articles:
![]() |
N. C. Darnton and H. C. Berg Force-Extension Measurements on Bacterial Flagella: Triggering Polymorphic Transformations Biophys. J., March 15, 2007; 92(6): 2230 - 2236. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |