Research Papers

A Model for Highly Strained DNA Compressed Inside a Protein Cavity

[+] Author and Article Information
Andrew D. Hirsh

Graduate Student
Research Assistant
Department of Mechanical Engineering,
University of Michigan,
Ann Arbor, MI 48109
e-mail: adhirsh@umich.edu

Todd D. Lillian

Assistant Professor
Department of Mechanical Engineering,
Texas Tech University,
Lubbock, TX 79409
e-mail: todd.lillian@ttu.edu

Troy A. Lionberger

Postdoctoral Research Associate
Howard Hughes Medical Institute,
Department of Physics,
University of California,
Berkeley, Berkeley, CA 94720
e-mail: talionberger@gmail.com

Maryna Taranova

Graduate Student
Research Assistant
e-mail: taranova.maryna@uci.edu

Ioan Andricioaei

Associate Professor
e-mail: andricio@uci.edu
Department of Chemistry,
University of California,
Irvine, Irvine, CA 92697

N. C. Perkins

Fellow ASME
Department of Mechanical Engineering,
University of Michigan,
Ann Arbor, MI 48109
e-mail: ncp@umich.edu

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received July 31, 2011; final manuscript received March 9, 2012; published online October 30, 2012. Assoc. Editor: Aki Mikkola.

J. Comput. Nonlinear Dynam 8(3), 031001 (Oct 30, 2012) (8 pages) Paper No: CND-11-1121; doi: 10.1115/1.4007535 History: Received July 31, 2011; Revised March 09, 2012

Deoxyribonucleic acid (DNA) is a long and flexible biopolymer that contains genetic information. Building upon the discovery of the iconic double helix over 50 years ago, subsequent studies have emphasized how its biological function is related to the mechanical properties of the molecule. A remarkable system which highlights the role of DNA bending and twisting is the packing and ejection of DNA into and from viral capsids. A recent 3D reconstruction of bacteriophage φ29 reveals a novel toroidal structure of highly bent/twisted DNA contained in a small cavity below the viral capsid. Here, we extend an elastic rod model for DNA to enable simulation of the toroid as it is compacted and subsequently ejected from a small volume. We compute biologically-relevant forces required to form the toroid and predict ejection times of several nanoseconds.

Copyright © 2013 by ASME
Topics: Cavities , DNA , Proteins
Your Session has timed out. Please sign back in to continue.


Watson, J. D. and Crick, F. H. C., 1953, “Molecular Structure of Nucleic Acids – A Structure for Deoxyribose Nucleic Acid,” Nature, 171(4356), pp. 737–738. [CrossRef] [PubMed]
Calladine, C. R., Drew, H. R., Luisi, B. F., and Travers, A. A., 2004, Understanding DNA: The Molecule and How It Works, 3rd ed., Elsevier, New York.
Gonzalez-Huici, V., Salas, M., and Hermoso, J. M., 2004, “The Push-Pull Mechanism of Bacteriophage Phi 29 DNA Injection,” Mol. Microbiol., 52(2), pp. 529–540. [CrossRef] [PubMed]
Phillips, R., Kondev, J., and Theriot, J., 2008. Physical Biology of the Cell, Garland Science, New York.
Smith, D. E., Tans, S. J., Smith, S. B., Grimes, S., Anderson, D. L., and Bustamante, C., 2001, “The Bacteriophage Phi 29 Portal Motor Can Package DNA Against a Large Internal Force,” Nature, 413(6857), pp. 748–752. [CrossRef] [PubMed]
Tang, J., Olson, N., Jardine, P. J., Girimes, S., Anderson, D. L. and Baker, T. S., 2008. “DNA Poised for Release in Bacteriophage Phi 29,” Structure, 16(6), pp. 935–943. [CrossRef] [PubMed]
Purohit, P. K., Inamdar, M. M., Grayson, P. D., Squires, T. M., Kondev, J., and Phillips, R., 2005, “Forces During Bacteriophage DNA Packaging and Ejection,” Biophys. J., 88(2), pp. 851–866. [CrossRef] [PubMed]
Spakowitz, A. J., and Wang, Z. G., 2005, “DNA Packaging in Bacteriophage: Is Twist Important?” Biophys. J., 88(6), pp. 3912–3923. [CrossRef] [PubMed]
Petrov, A. S., and Harvey, S. C., 2008. “Packaging Double-Helical DNA Into Viral Capsids: Structures, Forces, and Energetics,” Biophys. J., 95(2), pp. 497–502. [CrossRef] [PubMed]
Goyal, S., Perkins, N. C., and Lee, C. L., 2005, “Nonlinear Dynamics and Loop Formation in Kirchhoff Rods With Implications to the Mechanics of DNA and Cables,” J. Comput. Phys., 209(1), pp. 371–389. [CrossRef]
Goyal, S., Lillian, T., Blumberg, S., Meiners, J.-C., Meyhöfer, E., and Perkins, N. C., 2007. “Intrinsic Curvature of DNA Influences LacR-Mediated Looping,” Biophys. J., 93(12), pp. 4342–4359. [CrossRef] [PubMed]
Lillian, T. D., Goyal, S., Kahn, J. D., Meyhöfer, E., and Perkins, N. C., 2008. “Computational Analysis of Looping of a Large Family of Highly Bent DNA by LacI,” Biophys. J., 95(12), pp. 5832–5842. [CrossRef] [PubMed]
Goyal, S., Perkins, N. C., and Lee, C. L., 2008. “Non-linear Dynamic Intertwining of Rods With Self-Contact,” Int. J. Nonlinear Mech., 43(1), pp. 65–73. [CrossRef]
Lillian, T. D., and Perkins, N. C., 2011, “Electrostatics and Self-Contact in an Elastic Rod Approximation for DNA,” J. Comput. Nonlinear Dyn., 6(1), p. 011008. [CrossRef]
Lillian, T. D., Taranova, M., Wereszczynski, J., Andricioaei, I., and Perkins, N.C., 2011, “A Multiscale Dynamic Model of DNA Supercoil Relaxation by Topoisomerase IB,” Biophys. J., 100(8), pp. 2016–2023. [CrossRef] [PubMed]
Goyal, S., 2006, “A Dynamic Rod Model to Simulate Mechanics of Cables and DNA,” PhD thesis, University of Michigan, Ann Arbor, MI.
Schlick, T., 1995, “Modeling Superhelical DNA—Recent Analytical and Dynamic Approaches,” Curr. Opin. Struct. Biol., 5(2), pp. 245–262. [CrossRef] [PubMed]
Strick, T., Allemand, J. F., Croquette, V., and Bensimon, D., 2000, “Twisting and Stretching Single DNA Molecules,” Prog. Biophys. Mol. Biol., 74(1–2), pp. 115–140. [CrossRef] [PubMed]
Bustamante, C., Bryant, Z., and Smith, S. B., 2003. “Ten Years of Tension: Single-Molecule DNA Mechanics,” Nature, 421(6921), pp. 423–427. [CrossRef] [PubMed]
Beveridge, D. L., Young, M. A., and Sprous, D., 1998, “Modeling of DNA Via Molecular Dynamics Simulation: Structure, Bending, and Conformational Transitions,” Molecular Modeling of Nucleic Acids, American Chemical Society, Washington, DC, pp. 260–284.
Wereszczynski, J., and Andricioaei, I., 2006, “On Structural Transitions, Thermodynamic Equilibrium, and the Phase Diagram of DNA and RNA Duplexes Under Torque and Tension,” Proc. Natl. Acad. Sci. U.S.A., 103(44), pp. 16200–16205. [CrossRef] [PubMed]
Hagerman, P. J., 1988, “Flexibility of DNA,” Annu. Rev. Biophys. Chem., 17, pp. 265–286. [CrossRef]
Manning, R. S., Maddocks, J. H., and Kahn, J. D., 1996. “A Continuum Rod Model of Sequence-Dependent DNA Structure,” J. Chem. Phys., 105(13), pp. 5626–5646. [CrossRef]
Balaeff, A., Mahadevan, L., and Schulten, K., 2006. “Modeling DNA Loops Using the Theory of Elasticity,” Phys. Rev. E, 73(3), p. 031919. [CrossRef]
Chung, J., and Hulbert, G. M., 1993, “A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation—The Generalized-Alpha Method,” J. Appl. Mech., 60(2), pp. 371–375. [CrossRef]
Howard, J., 2001, Mechanics of Motor Proteins and the Cytoskeleton, Sinauer Associates, Sunderland, MA.
Rickgauer, J. P., Fuller, D. N., Grimes, S., Jardine, P. J., Anderson, D. L., and Smith, D. E., 2008, “Portal Motor Velocity and Internal Force Resisting Viral DNA Packaging in Bacteriophage Phi 29,” Biophys. J., 94(1), pp. 159–167. [CrossRef] [PubMed]
Duda, R. L., and Conway, J. F., 2008, “Asymmetric EM Reveals New Twists in Phage Phi 29 Biology,” Structure, 16(6), pp. 831–832. [CrossRef] [PubMed]
Lillian, T. D., Perkins, N. C., and Goyal, S., 2008, “Computational Elastic Rod Model Applied to DNA Looping,” Proceedings of ASME IDETC/CIE 2007, Las Vegas, NV, Sep. 4–7, Vol. 5, pp. 1449–1456.


Grahic Jump Location
Fig. 1

A schematic of mature bacteriophage φ29 based on the cryo-electron micoscopy reconstruction [6,28]. The toroidal DNA is contained within a cavity just below the capsid and suspected to contain 30–40 bp of DNA. The cavity is formed from the void between the connector and the lower collar. Shown inside the capsid is a cross sectional view of concentric hoops of DNA which ultimately fill the entire capsid.

Grahic Jump Location
Fig. 2

The atomic structure of DNA superimposed with a an elastic rod with equivalent elastic properties. R→(s,t) tracks the position of the helical axis as a function of contour length s and time t with respect to the inertial frame e. We also define a body-fixed reference frame a(s,t) which is also a function of s and t. Figure adaped from [29].

Grahic Jump Location
Fig. 3

Side view of the 3-D cavity structure estimated from the connector and lower collar geometry in the cryo-EM images [6] with relevant dimensions labeled. We assume the cavity is symmetric about the vertical axis. Cavity grid points are spaced in stacked rings of points separated by dc above and below one another and dc along the circumference of each ring (see red arrows). Here, we have set dc to 4 Å to reduce the number of grid points by an order of magnitude (over 1 Å separation) and thereby gain computational efficiency. The shaded blue box indicates the approximate area A=dc2 surrounding each grid point.

Grahic Jump Location
Fig. 4

The interaction forces are dependent on all pairwise vectors between rod grid points p and points q representing the cavity surface which are fixed in space

Grahic Jump Location
Fig. 5

(a) Computational snapshots (aligned to the dots on the plot directly below) from simulating 90 bp of DNA compressed within the cavity. (b) Internal (compressive) force (pN) and elastic energy (kT) following Eq. (6) as functions of the shortening of the rod δ = L – d where L is the length of DNA (306 Å) and d is the distance between the upper and lower rod boundaries. The green box illustrates where we report internal force along the contour length. (c) Top and side view of the final toroidal structure.

Grahic Jump Location
Fig. 6

Equilibrium conformation from Fig. 5(c) upon altering the cavity grid spacing parameter dc to (a) 8 Å, (b) 4 Å, (c) 2 Å, (d) 1 Å

Grahic Jump Location
Fig. 7

(a) Dynamic ejection and toroid collapse. Snapshots illustrate DNA conformation at each of the 5 ns increments denoted in the figure below. (b) DNA reaction force (solid green line and circles) and torque (solid red line and circles) on the remaining capsid/DNA (upper end) as a function of time. The solid blue line represents the component of the viscous drag force (integrated over the length of the rod) acting along the z-axis (vertical axis).



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In