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



* International School for Advanced Studies, Trieste, Italy;
Department of Physics "G. Galilei" of the University of Padua, Padua, Italy; and
Laboratory of Ultrastructural Analysis, Department of Biology and Medicine, University of Lausanne, 1015 Lausanne-Dorigny, Switzerland
Correspondence: Address reprint requests to Andrzej Stasiak, Tel.: 41-21-692-42-82; Fax: 41-21-692-41-05; E-mail: andrzej.stasiak{at}lau.unil.ch.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In recent years several independent studies investigated how ATP-hydrolysis-driven type II topoisomerases can selectively lower the frequency of DNA knotting (Rybenkov et al., 1997
; Yan et al., 1999
; Vologodskii et al., 2001
). However, there were no systematic studies that investigated the most likely relaxation path of a given DNA knot by a hypothetical topoisomerase that does not show a chirality bias and is just driven by the free energy gradient. Characterization of this relatively simple situation provides a necessary reference required for the understanding of such subtleties of topoisomerase action as the chirality bias (Roca, 2001
; Charvin et al., 2003
; Stone et al., 2003
; Schvartzman and Stasiak, 2004
; Trigueros et al., 2004
) or differential relaxation of various types of crossings in DNA twist knots (Mann et al., 2004
).
We investigate here the case of knotted polymers with 33 statistical Kuhn segments that form freely jointed isolateral polygons with 33 edges. We have selected this size since it corresponds to double-stranded DNA of bacteriophage P4 that provides a convenient experimental system of DNA knotting (Arsuaga et al., 2002a
,b
). However, our study has a more general perspective, as we aim to investigate generic knot spaces accessible to all polymers that are conveniently modeled as freely jointed, isosegemental chains (Vologodskii et al., 1974
; Smith et al., 1992
; Deutsch, 1999
; Katritch et al., 2000
; Arsuaga et al., 2002b
; Millett and Rawdon, 2003
).
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
| RESULTS |
|---|
|
|
|---|
Table 1 lists the probabilities with which a random passage in a given knot type (the starting knot types are indicated in the head of the table) leads to emergence of a new trajectory of the same or a different knot type (the types of resulting knots are indicated on the left side of the table). There are three sorts of data in the Table 1. Upper entries in corresponding cells are the transition probabilities with which a given starting knot can change its knot type to the respective resulting knots in an experiment in which passages resulting in no change of original topology are also taken into account. The middle entries (underlined) in diagonal fields of the table list probabilities that a given knot does not change its knot type after one segment-segment passage. The lower entries in respective cells result from a renormalization of the data by taking into account only these passages that have led to the change of the knot type. Looking at Table 1 one can see, for example, that random intersegmental passages in unknots (first column) led in almost 95% of cases back to the same knot type (underlined). Fig. 1 a presents an example of such a passage that does not change the knot type. Table 1 also shows that if only the passages that resulted in a change of the knot type were analyzed (lower entries), then right- and left-handed trefoil knots were both produced with almost 42% efficiency from unknots (see Methods section for the information how left-right symmetry was obtained from the data). Figure-of-eight knots, that are achiral, were produced in
13% of knot type-changing passages from unknots, whereas 52R and 52L knots were both produced with
1% probability. More complex knots emerged from unknots with still lower probabilities.
|
The observed one-passage connections within knots with up to six crossings agree with earlier analytical treatment of this subject (Darcy and Sumners, 1997
, 2000
). We observed transitions between all these knot spaces that are separated by one passage and didn't observe any when more than one passage is needed for such a change (Darcy and Sumners, 1997
, 2000
). We believe, therefore, that our simulation procedure probes adequately the connectivity between different knot spaces. Our studies go however further than revealing one-passage connectivity between different knot spaces since we obtain quantitative information about the extent of contact between the neighboring knot spaces.
Interfaces and shapes of individual knot spaces
By monitoring ratios of probabilities with which random configurations of studied knots exit by random passages to different neighboring knot spaces we can probe the effective interface area between corresponding knot spaces. The effective interface areas between different knot spaces are defined operationally as proportional to the number of observed passages between the knots in question. However, if different interfaces vary in their curvature and corrugation, then the actual interface areas may not be proportional to the corresponding effective interface areas. The lower entries in corresponding cells in Table 1 reflect the relative effective interface area of individual knot spaces with other knot spaces. Notice that the global interface of the starting knot space with all neighboring spaces is normalized to one for every individual knot type (sum of lower entries in every column is 1, taking into account that there are always two different chiral knots of the corresponding type). One can see, for example, that >41% of the surface enclosing the space of unknots forms an interface with the space of right-handed trefoils (first column, second row in Table 1). On the other side, >90% of the surface enclosing the configuration space occupied by right-handed trefoils is shared with the surface that encloses unknots (second column, first row). It should be obvious that the interface of the space of unknots with the space of right-handed trefoils has to be equal to the interface of the space of right-handed trefoils with the space of unknots. Therefore, one can simply conclude that the total surface enclosing the space of unknots is
2.16 times bigger than the corresponding surface enclosing right-handed trefoils.
The same type of calculation can be performed for any pair of knots with one-passage connectivity. Using this method we can express the delimiting surface of various knot spaces as a fraction of the total surface of all knot spaces summed together and taken as 1. For 33-segments long freely jointed polygons the knots 01, 31R or 31L, 41, 51R or 51L, 52R or 52L, 61R or 61L, 62R or 62L, 63, 31R#31R or 31L#31L, and 31R#31L have the following delimiting surfaces: 0.4518, 0.2090, 0.0655, 0.0091, 0.0128, 0.0022, 0.0031, 0.0035, 0.0015, and 0.0036, respectively, with a
10% error.
Diagonal entries in Table 1 (underlined) tell us how frequently configurations representing a given knot type remain in the same knot space after a random strand passage. It is visible that, as the knots get more complex their tendency to remain in the same knot type decreases. The probability of remaining in the same knot space after a random passage reflects the ratio between the area of internal invaginations (passage through which does not change the knot type) and the area of the external surfaces which separate different knot types. If the density of internal invaginations in various knot spaces is similar, the probability of remaining in the same knot space after a random passage reflects in an indirect way the volume-to-surface ratio of individual knot spaces. Thus, for example, a space of unknots seems to have the highest volume to surface ratio, whereas each of the spaces of various six-crossing knots has a much smaller volume-to-surface ratio (Table 1).
It was shown earlier that the so-called ideal geometric representations of knots decrease their volume-to-surface ratios with increasing complexity of knots (Katritch et al., 1996
). More recently it was proposed that the relative probability of occurrence of various knots at topological equilibrium should decrease exponentially with the length/diameter ratio (L/D) of their corresponding ideal geometric representations (Grosberg, 1998
).
The underlying hypothesis was that the diameter of tubes formed by ideal knots approximates the 3-D space available to a given knot type (Grosberg et al., 1996
; Grosberg, 1998
) and that this, in turn, is directly correlated to the probability of occurrence of the knot, or, in other words, to the volume of the high-dimensional configuration space available to that knot type.
This proposal, however, failed to explain why the probability of occurrence of 52 knots is significantly higher than that of 51 knots, despite the fact that 52 knots show higher L/D ratio of their ideal configurations (Deguchi and Tsurusaki, 1994
; Katritch et al., 1996
).
The core of the issue is, therefore, whether there is a relation between the physical 3-D space available for configurations of a given knot and the high dimensional configuration space occupied by the same knot type.
In the studied system the equilibrium probability of occurrence of a given knot type (volume of the high-dimensional configuration space) is the result of the complex dynamical process that, via intersegmental passages, allow the chain to maintain or change its knot type with different probabilities. The probability to stay in a given knot type after a pivotal type of move depends, however, on the extent of the 3-D space within which the configuration can move without self intersections leading to the change of the knot type. We suggest, therefore, that only the probability to maintain the knot topology is directly related to L/D ratio of ideal representations of knot. Fig. 2 shows that the probability of remaining in the same knot type upon a random segment-segment passage exhibits in fact an exponential decrease with L/D ratio of ideal knots and is consistent with the fact that 51 knots, despite their lower probability of formation, are more "stable" than 52 knots.
|
Intensity of exchanges between all investigated knot spaces
Fig. 3 illustrates the connectivity diagram of the analyzed knot types together with their expected probability. Fig. 3 presents only achiral knots and right-handed isoforms of chiral knots. Left-handed isoforms would be placed symmetrically on the left side of the "reflecting plane" occupied by achiral knots 01, 41, 63, and 31R#31L. Fig. 3 allows us to quickly find the minimal number of passages required for passing from one knot type to another. For example, to pass from the space of 51R knots to that of unknots, two passages are needed and two more would be needed to reach the space of 51L knots. In Fig. 3 we also indicated the intensity of exchanges or fluxes between neighboring knot spaces at thermodynamic equilibrium. The fluxes resulting from random passages between neighboring knot spaces are normalized with respect to the flux between unknots and each of trefoil-knots (the intensity of which was taken as 1).
|
|
Knot spaces are usually sampled using a phantom chain evolution; however, such a sampling can only provide information about expected stationary probabilities of different knots without revealing the overall organization of the entire knot space (Rybenkov et al., 1993
; Deguchi and Tsurusaki, 1994
; Katritch et al., 2000
). It is important to mention though, that the expected stationary probabilities obtained with phantom sampling closely coincided with these obtained from nonphantom exploration (data not shown) indicating that the structure of the foamy knot space is not affected by the particular polymer dynamics applied in our simulation studies.
Using the topoisomerase-inspired approach presented here, we were able to verify earlier theoretical studies that investigated which knots are neighbors in the knots configuration space (Darcy and Sumners, 1997
, 2000
). In addition, we were able to characterize the extent of contacts, and thus exchanges, between configuration spaces of different knots. The probabilities of staying in the same knot type after a random passage provided us with an estimation of the surface-to-volume ratio of configuration spaces of various knots.
Biological applications
The presented methodology can be applied to quantify biases manifested by some DNA topoisomerases during relaxation of particular types of DNA knots. It was reported recently that type II topoisomerases act preferentially on clasp crossings in twist 52 knots (Mann et al., 2004
). However, the actual extent of the bias can be only quantified when one knows what the probability of acting on the clasp crossing in the absence of any bias is. Table I shows that without any bias the action on the clasp crossing on 52 knot tied on 10-kb-long DNA should happen in
36% of strand passages that result in the change of the knot type. Therefore, if the probability of action of the investigated topoisomerase II on the clasp crossing is significantly >36% the bias can be confirmed.
Another example where the probabilities of passages between different knot types can be used to obtain biologically significant information concerns the geometry of DNA packing inside phage heads. The DNA knots formed within phage heads can provide hints about the nature of DNA packing (Arsuaga et al., 2002b
) but the knots are too complex for their knot type to be determined by electron-microscopy (Krasnow et al., 1983
) or by their electrophoretic migration (Vologodskii et al., 1998
). However, upon partial relaxation knots become simpler and can be correctly recognized by their electrophoretic migration (Arsuaga et al., 2002b
). The methodology presented here should allow finding out what were the knots that produced the partially relaxed knots of a known knot type. The probability of a passage from one knot to another as a result of individual topoisomerase reactions can not be obtained from the classical exploration of the knot space where the equilibrium probabilities of different knots are calculated (Rybenkov et al., 1993
).
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This research was funded in part by Swiss National Science Foundation grants 3152-068151 and 3100A0-103962.
| FOOTNOTES |
|---|
Submitted on May 12, 2004; accepted for publication August 6, 2004.
| REFERENCES |
|---|
|
|
|---|
Arsuaga, J., R. K. Tan, M. Vazquez, D. W. Sumners, and S. C. Harvey. 2002a. Investigation of viral DNA packaging using molecular mechanics models. Biophys. Chem. 101102:475484.[CrossRef]
Arsuaga, J., M. Vazquez, S. Trigueros, D. W. Sumners, and J. Roca. 2002b. Knotting probability of DNA molecules confined in restricted volumes: DNA knotting in phage capsids. Proc. Natl. Acad. Sci. USA. 99:53735377.
Buck, G. 1998. Four-thirds power law for knots and links. Nature. 392:238239.[CrossRef]
Champoux, J. J. 2001. DNA topoisomerases: structure, function, and mechanism. Annu. Rev. Biochem. 70:369413.[CrossRef][Medline]
Charvin, G., D. Bensimon, and V. Croquette. 2003. Single-molecule study of DNA unlinking by eukaryotic and prokaryotic type-II topoisomerases. Proc. Natl. Acad. Sci. USA. 100:98209825.
Darcy, I. K., and D. W. Sumners. 1997. A strand passage metric for topoisomerase action. In Knots '96: Proceedings of the Fifth MSJ International Research Institute of Mathematical Society of Japan. S. Suzuki, editor. World Scientific, Singapore.
Darcy, I. K., and D. W. Sumners. 2000. Rational tangle distance on knots and links. Math. Proc. Cambridge Philos. Soc. 128:497510.[CrossRef]
Dean, F. B., A. Stasiak, T. Koller, and N. R. Cozzarelli. 1985. Duplex DNA knots produced by Escherichia coli topoisomerase I, structure and requirements for formation. J. Biol. Chem. 260:49754983.
Deguchi, T., and K. Tsurusaki. 1994. A statistical study of random knotting using the Vassiliev invariants. Journal of Knot Theory and Its Ramifications. 3:321353.[CrossRef]
Deutsch, J. M. 1999. Equilibrium size of large ring molecules. Phys. Rev. E. 59:R2539R2541.[CrossRef]
Dobay, A., J. Dubochet, K. Millett, P.-E. Sottas, and A. Stasiak. 2003. Scaling behavior of random knots. Proc. Natl. Acad. Sci. USA. 100:56115615.
Ewing, B., and K. C. Millett. 1996. Computational algorithms and the complexity of link polynomials. In Progress in Knot Theory and Related Topics, Trauaux en Cours. 56. Hermann. Paris. 5168.
Freyd, P., D. Yetter, J. Hoste, W. Lickorish, K. Millett, and A. Ocneau. 1985. A new polynomial invariant of knots and links. 12:239246.
Grimmet, G., and D. Stirzaker. 1992. Probability and Random Processes. Oxford University Press, Oxford.
Grosberg, A. Y. 1998. Entropy of a knot: simple arguments about difficult problem. In Ideal Knots. A. Stasiak, S. Katritch, and L. H. Kauffman, editors. World Scientific, Singapore. 129142.
Grosberg, A. Y., A. Feigel, and A. Rabin. 1996. Flory-type theory of a knotted ring polymer. Phys. Rev. E. 54:66186622.[CrossRef]
Janse Van Rensburg, E. J., D. W. Sumners, and S. G. Whittington. 1998. The writhe of knots and links. In Ideal Knots. A. Stasiak, S. Katritch, and L. H. Kauffman, editors. World Scientific, Singapore. 7087.
Katritch, V., J. Bednar, D. Michoud, R. G. Scharein, J. Dubochet, and A. Stasiak. 1996. Geometry and physics of knots. Nature. 384:142145.[CrossRef]
Katritch, V., W. K. Olson, A. V. Vologodskii, J. Dubochet, and A. Stasiak. 2000. Tightness of random knotting. Phys. Rev. E. 61:55455549.[CrossRef]
Krasnow, M. A., A. Stasiak, S. J. Spengler, F. Dean, T. Koller, and N. R. Cozzarelli. 1983. Determination of the absolute handedness of knots and catenanes of DNA. Nature. 304:559560.[CrossRef][Medline]
Le Bret, M. 1980. Monte Carlo computation of the supercoiling energy, the sedimentation constant and the radius of gyration of unknotted and knotted circular DNA. Biopolymers. 19:619637.[CrossRef][Medline]
Mann, J. K., R. W. Deibler, D. W. Sumners, and E. L. Zechiedrich. 2004. Unknotting by type II topoisomerases. Abstracts of papers presented to the American Mathematical Society. 25:99492187.
Millett, K. C., and E. J. Rawdon. 2003. Energy, rope length, and other physical aspects of equilateral knots. J. Comp. Phys. 186:426456.[CrossRef]
Roca, J. 2001. Varying levels of positive and negative supercoiling differently affect the efficiency with which topoisomerase II catenates and decatenates DNA1. J. Mol. Biol. 305:441450.[CrossRef][Medline]
Rolfsen, D. 1976. Knots and Links. Publish or Perish Books, Berkeley, CA.
Rybenkov, V. V., N. R. Cozzarelli, and A. V. Vologodskii. 1993. The probability of DNA knotting and the effective diameter of the DNA double helix. Proc. Natl. Acad. Sci. USA. 90:53075311.
Rybenkov, V. V., C. Ulsperger, A. V. Vologodskii, and N. R. Cozzarelli. 1997. Simplification of DNA topology below equilibrium values by type II topoisomerases. Science. 277:690693.
Schvartzman, J. B., and A. Stasiak. 2004. A topological view of the replicon. EMBO Rep. 5:256261.[CrossRef][Medline]
Shaw, S. Y., and J. C. Wang. 1993. Knotted DNA rings: probability of formation and resolution of the two chiral trefoils. Science. 260:533536.
Smith, S. B., L. Finzi, and C. Bustamante. 1992. Direct mechanical measurements of the elasticity of single DNA molecules by using magnetic beads. Science. 258:11221126.
Stasiak, A. 2003. Topoisomerases. In Encyclopedia of the Human Genome. Nature Publishing Group, London.
Stasiak, A., J. Dubochet, V. Katritch, and P. Pieranski. 1998. Ideal knots and their relation to the physics of real knots. In Ideal Knots. A. Stasiak, S. Katritch, and L. H. Kauffman, editors. World Scientific, Singapore. 119.
Stone, M. D., Z. Bryant, N. J. Crisona, S. B. Smith, A. Vologodskii, C. Bustamante, and N. R. Cozzarelli. 2003. Chirality sensing by Escherichia coli topoisomerase IV and the mechanism of type II topoisomerases. Proc. Natl. Acad. Sci. USA. 100:86548659.
Trigueros, S., J. Salceda, I. Bermudez, X. Fernandez, and J. Roca. 2004. Asymmetric removal of supercoils suggests how topoisomerase II simplifies DNA topology. J. Mol. Biol. 335:723731.[CrossRef][Medline]
Vologodskii, A. V., N. J. Crisona, B. Laurie, P. Pieranski, V. Katritch, J. Dubochet, and A. Stasiak. 1998. Sedimentation and electrophoretic migration of DNA knots and catenanes. J. Mol. Biol. 278:13.[CrossRef][Medline]
Vologodskii, A. V., S. D. Levene, K. V. Klenin, M. Frank-Kamenetskii, and N. R. Cozzarelli. 1992. Conformational and thermodynamic properties of supercoiled DNA. J. Mol. Biol. 227:12241243.[CrossRef][Medline]
Vologodskii, A. V., A. V. Lukashin, M. D. Frank-Kamenetskii, and V. V. Anshelevitch. 1974. The knot problem in statistical mechanics of polymer chains. Sov. Phys. JETP (Engl. Transl.). 39:10591063.
Vologodskii, A. V., W. Zhang, V. V. Rybenkov, A. A. Podtelezhnikov, D. Subramanian, J. D. Griffith, and N. R. Cozzarelli. 2001. Mechanism of topology simplification by type II DNA topoisomerases. Proc. Natl. Acad. Sci. USA. 98:30453049.
Yan, J., M. O. Magnasco, and J. F. Marko. 1999. A kinetic proofreading mechanism for disentanglement of DNA by topoisomerases. Nature. 401:932935.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
I. K. Darcy, R. G. Scharein, and A. Stasiak 3D visualization software to analyze topological outcomes of topoisomerase reactions Nucleic Acids Res., June 1, 2008; 36(11): 3515 - 3521. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Liu, E. L. Zechiedrich, and H. S. Chan Inferring Global Topology from Local Juxtaposition Geometry: Interlinking Polymer Rings and Ramifications for Topoisomerase Action Biophys. J., April 1, 2006; 90(7): 2344 - 2355. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |