Electrostatic interaction between long, rigid helical macromolecules at all interaxial angles

Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

Volumn 62, Issue 2, 2000, Pages 2576-2596

  Kornyshev, A A ^a,b   Leikin, S ^c  

^a UNIVERSITY OF CALIFORNIA   (United States)

^b FORSCHUNGSZENTRUM JÜLICH   (Germany)

^c EUNICE KENNEDY SHRIVER NATIONAL INSTITUTE OF CHILD HEALTH AND HUMAN DEVELOPMENT   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

ELECTROLYTE SOLUTION; ELECTROSTATIC INTERACTION; HELICAL MACROMOLECULES;

APPROXIMATION THEORY; BIOELECTRIC POTENTIALS; DNA; ELECTROLYTES; FOURIER TRANSFORMS; INTEGRAL EQUATIONS; MATHEMATICAL MODELS; MOLECULAR STRUCTURE; PHASE TRANSITIONS;

MOLECULAR DYNAMICS;

DNA; ELECTROLYTE; SOLVENT;

CHEMICAL STRUCTURE; CHEMISTRY; CONFORMATION; METABOLISM; MOTION; STATIC ELECTRICITY; THERMODYNAMICS; TORQUE;

DNA; ELECTROLYTES; MODELS, MOLECULAR; MOTION; NUCLEIC ACID CONFORMATION; SOLVENTS; STATIC ELECTRICITY; THERMODYNAMICS; TORQUE;

EID: 0034238007     PISSN: 1063651X     EISSN: None     Source Type: Journal    
DOI: 10.1103/PhysRevE.62.2576     Document Type: Article

References (70)

1. F. Livolant, J. Phys. (Paris) 47, 1605 (1986).
2. T. E. Strzelecka and R. L. Rill, J. Am. Chem. Soc. 109, 4513 (1987).
3. T. E. Strzelecka, M. W. Davidson, and R. L. Rill, Nature (London) 331, 457 (1988).
4. Yu. M. Yevdokimov, S. G. Skurdin, and V. I. Salyanov, Liq. Cryst. 3, 1443 (1988).
5. F. Livolant, Physica A 176, 117 (1991).
6. A. Leforestier and F. Livolant, Biophys. J. 65, 56 (1993).
7. C. Robinson, Tetrahedron 13, 219 (1961).
8. D. B. Dupre and R. W. Duke, J. Chem. Phys. 63, 143 (1975).
9. M. M. Giraud-Guille, Int. Rev. Cytol. 166, 59 (1996).
10. D. H. Van Winkle, M. W. Davidson, W.-H. Chen, and R. L. Rill, Macromolecules 23, 4140 (1990).
11. R. Podgornik, H. Strey, and V. A. Parsegian, Curr. Opin. Colloid Interface Sci. 3, 534 (1998).
12. H. H. Strey, J. Wang, R. Podgornik, A. Rupprecht, L. Yu, V. A. Parsegian, and E. B. Sirota, Phys. Rev. Lett. 84, 3105 (2000).
13. R. D. Kamien, Mol. Cryst. Liq. Cryst. Sci. Technol., Sect. A 299, 265 (1997).
14. A. B. Harris, R. D. Kamien, and T. C. Lubensky, Phys. Rev. Lett. 78, 1476 (1997).
15. A. B. Harris, R. D. Kamien, and T. C. Lubensky, Rev. Mod. Phys. 71, 1745 (1999).
16. T. C. Lubensky, A. B. Harris, R. D. Kamien, and G. Yan, Ferroelectrics 212, 1 (1998).
17. L. Onsager, Ann. (N.Y.) Acad. Sci. 51, 627 (1949).
18. S. L. Brenner and V. A. Parsegian, Biophys. J. 14, 327 (1974).
19. J. L. Barrat and J. F. Joanny, Adv. Chem. Phys. 94, 1 (1996).
20. R. Podgornik and V. A. Parsegian, Phys. Rev. Lett. 80, 1560 (1998).
21. B.-Y. Ha and A. J. Liu, Europhys. Lett. 46, 624 (1999).
22. J. M. Bailey, Biopolymers 12, 559 (1973).
23. D. Soumpasis, J. Chem. Phys. 69, 3190 (1978).
24. B. Jayaram, K. A. Sharp, and B. Honig, Biopolymers 28, 975 (1989).
25. K. Wagner, E. Keyes, T. W. Kephart, and G. Edwards, Biophys. J. 73, 21 (1997).
26. B. Jayaram and D. L. Beveridge, J. Phys. Chem. 94, 4666 (1990).
27. A. A. Kornyshev and S. Leikin, J. Chem. Phys. 107, 3656 (1997)
28. J. Chem. Phys.A. A. KornyshevS. Leikin108, 7035(E) (1998).
29. A. A. Kornyshev and S. Leikin, Biophys. J. 75, 2513 (1998).
30. A. A. Kornyshev and S. Leikin, Proc. Natl. Acad. Sci. USA 95, 13579 (1998).
31. A. A. Kornyshev and S. Leikin, Phys. Rev. Lett. 82, 4138 (1999).
32. Such Hamiltonian includes both electrostatic and nonelectrostatic terms. Generally, the latter are not well known.
33. F. Oosawa, Polyelectrolytes (Marcel Dekker, New York, 1971).
34. B.-Y. Ha and A. J. Liu, Phys. Rev. Lett. 79, 1289 (1997).
35. B.-Y. Ha and A. J. Liu, Phys. Rev. E 60, 803 (1999).
36. P. Mariani and L. Saturni, Biophys. J. 70, 2867 (1996).
37. P. Mariani, F. Ciuchi, and L. Saturni, Biophys. J. 74, 430 (1998).
38. V. A. Bloomfield and I. Rouzina, Methods Enzymol. 295, 364 (1999).
39. We neglect dielectric saturation of the solvent in immediate vicinity of ions since we do not apply our results at such short distances. Furthermore, we do not account for nonlocal dielectric response of water and neglect hydration forces that may be caused by specific water ordering near molecular groups.
40. M. E. Fisher and Y. Levin, Phys. Rev. Lett. 71, 3826 (1993).
41. B. P. Lee and M. E. Fisher, Phys. Rev. Lett. 76, 2906 (1996).
42. M. E. Fisher, J. Stat. Phys. 75, 1 (1994).
43. D. M. Zuckerman, M. E. Fisher, and B. P. Lee, Phys. Rev. E 56, 6569 (1997).
44. DNA contains fixed phosphate charges and bound counterions so that (Formula presented) is the surface charge density of phosphates. Protein helices typically contain both negative and positive charges. If fixed charges of one sign dominate, we define (Formula presented) as their mean surface charge density. If a helix contains equal number of positive and negative charges, either one can be used to define (Formula presented)
45. See, e.g., W. Saenger, Principles of Nucleic Acid Structure (Springer, New York, 1984).
46. C. G. Baumann, S. B. Smith, V. A. Bloomfield, and C. Bustamante, Proc. Natl. Acad. Sci. USA 94, 6185 (1997).
47. Based on the crystallographic data 44 the deviation of the axial rise per base pair (bp) from its average value 3.34 Å is about 0.1 Å. For 150-bp -long DNA with random base pair sequence, this leads to axial deviation from the ideal helix (Formula presented) that is negligibly small compared to the helical pitch (Formula presented) Similarly, the deviation of the angle between adjacent base pairs from its average value 36° is ∼4°. The deviation over the whole length of the molecule is, therefore, ∼50°, which is still small enough compared to the full 360° turn that would disrupt orientational correlations between strands on opposing molecules.
48. The average distance between neighboring charges on DNA surface is ∼7 Å, much smaller than the helical pitch and than the separation between molecules in the cholesteric phase. Even when the molecules are in direct contact with each other, the effects from charge discreteness contribute only ∼10% to the interaction energy 29. In the cholesteric phase, their contribution should be much smaller.
49. Equation (17) differs from the result reported in Ref. 18 by a factor of 2 due to a misprint there [V. A. Parsegian, (private communication)].
50. The radius of the cylinder formed by centers of phosphate groups on B-DNA is (Formula presented)
51. The distance range where this approximation is valid depends on the symmetry of molecular charge pattern that determines the “weight factors” (Formula presented) For example, when (Formula presented) at all (Formula presented) this approximation is valid at all separations. In contrast, when (Formula presented) this approximation cannot be used at all.
52. In Eq. (43) we used that zero surface-to-surface separation corresponds to (Formula presented)
53. D. Durand, J. Doucet, and F. Livolant, J. Phys. II 2, 1769 (1992).
54. F. Livolant and A. Leforestier, Mol. Cryst. Liq. Cryst. 215, 47 (1992).
55. G. S. Manning, Q. Rev. Biophys. 11, 179 (1978).
56. The exact nature of this phase is not entirely clear. Originally, it was reported as a hexagonal columnar phase (see, e.g., 56 and references therein). Later studies of phase behavior of long DNA raised the possibility that this may be a line hexatic pase 1256. It is clear, however, that this phase is not chiral, the averaged local packing of molecules is hexagonal, and the averaged orientation of all molecules is the same.
57. R. Podgornik, H. H. Strey, K. Gawrisch, D. C. Rau, A. Rupprecht, and V. A. Parsegian, Proc. Natl. Acad. Sci. USA 93, 4261 (1996).
58. H. H. Strey, V. A. Parsegian, and R. P. Podgornik, Phys. Rev. E 59, 999 (1999).
59. This is true only in the presence of at least some mutual alignment between nearest neighbors. At completely uncorrelated rotations of molecules about their principal axes, only quantum interactions between fluctuating charge moments contribute to the net torque [S. A. Issaenko, A. B. Harris, and T. C. Lubensky, Phys. Rev. E 60, 578 (1999)]. In the limit of strong correlations in the rotations, the quantum torque is several orders of magnitude weaker than the electrostatic one.
60. P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Clarendon, Oxford, 1993).
61. B. W. van der Meer, G. Vetrogen, A. J. Dekker, and J. G. J. Yptma, J. Chem. Phys. 65, 3935 (1976).
62. Based on Eqs. (23242526), such screw motion has zero energetic cost. However, to derive Eqs. (23242526) we neglected edge effects for finite length helices. The edge effects will create some energetic cost, but to the extent that they are weak, the cost should be small.
63. The persistence length of DNA (500–1000 Å 45) is larger than the helical pitch by more than an order of magnitude and it is even larger than the total length of 150-bp fragments. In the first order approximation, undulations are essentially equivalent to long-wavelength fluctuations of the interaxial distance that do not change the (Formula presented) alignment of neighboring molecules.
64. R. Podgornik, D. C. Rau, and V. A. Parsegian, Macromolecules 22, 1780 (1989).
65. Since molecular packing in the cholesteric phase is nearly hexagonal, the distance between neighboring layers is (Formula presented) The rotation per layer is ψ. The pitch corresponds to the rotation by 2π.
66. L. Salem, X. Chapuisat, G. Segal, P. C. Hiberty, C. Minot, C. Leforestier, and P. Sautet, J. Am. Chem. Soc. 109, 2887 (1987).
67. F. Livolant and M. F. Maestre, Biochemistry 27, 3056 (1988).
68. Circular dichroism of DNA in the base pair absorption band (∼260 nm) is due to electronic transitions caused by interaction of the electric vector of light with the conjugated π electron system of DNA bases. The corresponding transition dipoles lie in the plane of the bases which is normal 44 to the principal axis of B-DNA. This plane is parallel to the cholesteric twist axis and, correspondingly, to the direction in which the anomalous CD signal is supposed to be observed. Then, only one component of the electric vector of light is parallel to transition dipoles while the second component is normal to them and it cannot cause electronic transitions. Only the first component can be absorbed so that B-DNA should not be able to distinguish left from right circularly polarized light in the direction of the cholesteric twist axis. Of course, real B-DNA molecules may have base pairs that are not perfectly parallel to the cholesteric twist axis so that some circular dichroism may be observed in this direction. However, the corresponding CD signal should be small rather than anomalously large.
69. It may be useful to recall a similar measurement of the twist sense in the cholesteric phase of xanthan (a charged, double-stranded, right-handed helical polysaccharide) [G. Maret, M. Milas, and M. Rinaudo, Polym. Bull. 4, 291 (1981)]. The anomalous optical rotary dispersion in the absorption band of xanthan was attributed to the left-handed cholesteric twist. However, later direct measurements revealed a right-handed cholesteric twist 1.
70. Note that for biological helices immersed into an aqueous electrolyte solution (Formula presented) because (Formula presented) In this case, Eq. (A15) can be simplified to yield the expression (Formula presented) used in our previous work 27. However, here we will also discuss interactions between macromolecules in nonpolar media when (Formula presented) and (Formula presented) Then, such approximation cannot be used and, thus, we avoid it.