Photosystem trap energies and spectrally-dependent energy-storage efficiencies in the Chl d-utilizing cyanobacterium, Acaryochloris marina.

Acaryochloris marina is the only species known to utilize chlorophyll (Chl) d as a principal photopigment. The peak absorption wavelength of Chl d is redshifted ≈40nm in vivo relative to Chl a, enabling this cyanobacterium to perform oxygenic phototrophy in niche environments enhanced in far-red light. We present measurements of the in vivo energy-storage (E-S) efficiency of photosynthesis in A. marina, obtained using pulsed photoacoustics (PA) over a 90-nm range of excitation wavelengths in the red and far-red. Together with modeling results, these measurements provide the first direct observation of the trap energies of PSI and PSII, and also the photosystem-specific contributions to the total E-S efficiency. We find the maximum observed efficiency in A. marina (40±1% at 735nm) is higher than in the Chl a cyanobacterium Synechococcus leopoliensis (35±1% at 690nm). The efficiency at peak absorption wavelength is also higher in A. marina (36±1% at 710nm vs. 31±1% at 670nm). In both species, the trap efficiencies are ≈40% (PSI) and ≈30% (PSII). The PSI trap in A. marina is found to lie at 740±5nm, in agreement with the value inferred from spectroscopic methods. The best fit of the model to the PA data identifies the PSII trap at 723±3nm, supporting the view that the primary electron-donor is Chl d, probably at the accessory (Chl(D1)) site. A decrease in efficiency beyond the trap wavelength, consistent with uphill energy transfer, is clearly observed and fit by the model. These results demonstrate that the E-S efficiency in A. marina is not thermodynamically limited, suggesting that oxygenic photosynthesis is viable in even redder light environments.

Persistence lengths of DNA obtained from Brownian dynamics simulations.

The persistence length of DNA has been studied for decades; however, experimentally obtained values of this quantity have not been entirely consistent. We report results from Brownian dynamics simulations that address this issue, validating and demonstrating the utility of an explicitly double-stranded model for mesoscale DNA dynamics. We find that persistence lengths calculated from rotational relaxation increase with decreasing ionic strength, corroborating experimental evidence, but contradicting results obtained from wormlike coil assumptions. Further, we find that natural curvature does not significantly affect the persistence length, corroborating cyclization efficiency measurements, but contradicting results from cryo-EM.

Brownian dynamics of double-stranded DNA in periodic systems with discrete salt.

Numerical models of mesoscale DNA dynamics relevant to in vivo scenarios require methods that incorporate important features of the intracellular environment, while maintaining computational tractability. Because the explicit inclusion of ions leads to electrostatic calculations that scale as the square of the number of charged particles, such models typically handle these calculations using low-potential, mean-field approaches, rather than by considering the discrete interactions of ions. This allows approximation of the long-range, screened self-repulsion of DNA, but is unable to capture detailed electrostatic phenomena, such as short-range attractions mediated by ion-ion correlations. Here, we develop a dynamical model of explicitly double-stranded, sequence-specific DNA in a bulk environment consisting of other polyions and explicitly represented counterions and coions. DNA is represented as two interwound chains of charged Stokes spheres, and ions as free, monovalently charged Stokes spheres. Brownian dynamics simulations performed at salt concentrations of 0.1, 1, 10, and 100 mM demonstrate this model captures anticipated behaviors of the system, including increasing compaction of the polyion by the ionic atmosphere with increasing ionic strength. The decay of the distance dependence of the ion concentrations as one moves away from the polyion approaches their equilibrium values in quantitative agreement with predictions of Poisson-Boltzmann theory. The simulation results also demonstrate quantitative agreement with experimental measurements of the persistence length of B-DNA, which increases significantly at low ionic strengths. The model also captures behaviors intimating the importance of explicitly representing ionic and polyionic structure. These include penetration of the polyion interior by both coions and counterions, and counterion-mediated accumulation of coions near the surface of the polyion. Such phenomena are likely to play an important role in the formation of alternative DNA secondary structures, suggesting the present methods will prove valuable to dynamic models of superhelical stress-induced DNA structural transitions.

Visualization of the detailed structure of plasmid DNA.

We report results from scanning tunneling microscopy experiments in which supercoiling of the plasmid, pTNT, leads to localized duplex denaturation at a single location. Measurements of the imaged strand-separated region suggest an extent of approximately 60 base pairs (bp's), in excellent agreement with the prediction of a statistical mechanical analysis of the pTNT base sequence that accounts for the topology-mediated global coupling of melting behaviors. Notwithstanding that the sequence of pTNT includes a 30-bp run of poly-A, where thermodynamic considerations alone predict significant instability, based on the statistical result, we propose the location of denaturation is an approximately 60-bp run, spanning positions 2050-2110, that is 82% AT-rich, terminally flanks the beta-lactamase coding region, and is known to contain two sites susceptible to cleavage by the DraI restriction enzyme.

DNA mechanics.

We review the history of DNA mechanics and its analysis. We evaluate several methods to analyze the structures of superhelical DNA molecules, each predicated on the assumption that DNA can be modeled with reasonable accuracy as an extended, linearly elastic polymer. Three main approaches are considered: mechanical equilibrium methods, which seek to compute minimum energy conformations of topologically constrained molecules; statistical mechanical methods, which seek to compute the Boltzmann distribution of equilibrium conformations that arise in a finite temperature environment; and dynamic methods, which seek to compute deterministic trajectories of the helix axis by solving equations of motion. When these methods include forces of self-contact, which prevent strand passage and preserve the topological constraint, each predicts plectonemically interwound structures. On the other hand, the extent to which these mechanical methods reliably predict energetic and thermodynamic properties of superhelical molecules is limited, in part because of their inability to account explicitly for interactions involving solvent. Monte Carlo methods predict the entropy associated with supercoiling to be negative, in conflict with a body of experimental evidence that finds it is large and positive, as would be the case if superhelical deformations significantly disrupt the ordering of ambient solvent molecules. This suggests that the large-scale conformational properties predicted by elastomechanical models are not the only ones determining the energetics and thermodynamics of supercoiling. Moreover, because all such models that preserve the topological constraint correctly predict plectonemic interwinding, despite these and other limitations, this constraint evidently dominates energetic and thermodynamic factors in determining supercoil geometry. Therefore, agreement between predicted structures and structures obtained experimentally, for example, by electron microscopy, does not in itself provide evidence for the correctness or completeness of any given model of DNA mechanics.

Brownian dynamics simulations of sequence-dependent duplex denaturation in dynamically superhelical DNA.

The topological state of DNA in vivo is dynamically regulated by a number of processes that involve interactions with bound proteins. In one such process, the tracking of RNA polymerase along the double helix during transcription, restriction of rotational motion of the polymerase and associated structures, generates waves of overtwist downstream and undertwist upstream from the site of transcription. The resulting superhelical stress is often sufficient to drive double-stranded DNA into a denatured state at locations such as promoters and origins of replication, where sequence-specific duplex opening is a prerequisite for biological function. In this way, transcription and other events that actively supercoil the DNA provide a mechanism for dynamically coupling genetic activity with regulatory and other cellular processes. Although computer modeling has provided insight into the equilibrium dynamics of DNA supercoiling, to date no model has appeared for simulating sequence-dependent DNA strand separation under the nonequilibrium conditions imposed by the dynamic introduction of torsional stress. Here, we introduce such a model and present results from an initial set of computer simulations in which the sequences of dynamically superhelical, 147 base pair DNA circles were systematically altered in order to probe the accuracy with which the model can predict location, extent, and time of stress-induced duplex denaturation. The results agree both with well-tested statistical mechanical calculations and with available experimental information. Additionally, we find that sites susceptible to denaturation show a propensity for localizing to supercoil apices, suggesting that base sequence determines locations of strand separation not only through the energetics of interstrand interactions, but also by influencing the geometry of supercoiling.

Transcription-driven twin supercoiling of a DNA loop: a Brownian dynamics study.

The torque generated by RNA polymerase as it tracks along double-stranded DNA can potentially induce long-range structural deformations integral to mechanisms of biological significance in both prokaryotes and eukaryotes. In this paper, we introduce a dynamic computer model for investigating this phenomenon. Duplex DNA is represented as a chain of hydrodynamic beads interacting through potentials of linearly elastic stretching, bending, and twisting, as well as excluded volume. The chain, linear when relaxed, is looped to form two open but topologically constrained subdomains. This permits the dynamic introduction of torsional stress via a centrally applied torque. We simulate by Brownian dynamics the 100 micros response of a 477-base pair B-DNA template to the localized torque generated by the prokaryotic transcription ensemble. Following a sharp rise at early times, the distributed twist assumes a nearly constant value in both subdomains, and a succession of supercoiling deformations occurs as superhelical stress is increasingly partitioned to writhe. The magnitude of writhe surpasses that of twist before also leveling off when the structure reaches mechanical equilibrium with the torsional load. Superhelicity is simultaneously right handed in one subdomain and left handed in the other, as predicted by the "transcription-induced twin-supercoiled-domain" model [L. F. Liu and J. C. Wang, Proc. Natl. Acad. Sci. U.S.A. 84, 7024 (1987)]. The properties of the chain at the onset of writhing agree well with predictions from theory, and the generated stress is ample for driving secondary structural transitions in physiological DNA.