DNA persistence length revisited

DNA restriction fragments ranging from 79 to 789 base pairs in length have been characterized by transient electric birefringence (TEB) measurements at various temperatures between 4 and 43 degrees C. The DNA fragments do not contain runs of four or more adenine residues in a row and migrate with normal electrophoretic mobilities in polyacrylamide gels, indicating that they are not intrinsically curved or bent. The low ionic strength buffers used for the measurements contained 1 mM Tris Cl, pH 8.0, EDTA, and variable concentrations of Na(+) or Mg(2+) ions. The rotational relaxation times were obtained by fitting the TEB field-free decay signals with a nonlinear least-squared fitting program; the decay of the birefringence was monoexponential for fragments < or = 241 base pair (bp) in length and multiexponential for larger fragments. The terminal relaxation times, characteristic of the end-over-end rotation of the DNA molecules, were then used to determine the persistence length (p) and hydrodynamic radius (r) of DNA as a function of temperature and ionic strength, using several different hydrodynamic models. The specific values obtained for p and r are model dependent. The wormlike chain model of P. J. Hagerman and B. H. Zimm (Biopolymers 1981, Vol. 20, pp. 1481-1502) combined with the revised Broersma equation (J. Newman et al., Journal of Mol Biol 1997, Vol. 116, pp. 593-606) appears to be the most suitable for describing the flexibility of DNA in low ionic strength solutions. The values of p and r obtained from the global least squares fitting of this equation are independent of DNA length, and the deviations of the individual values from the average are reasonably small. The consensus r value calculated for DNA in various low ionic strength solutions containing 1 mM Tris buffer is 14.7 +/- 0.4 A at 20 degrees C. The consensus p values decrease from 814 approximately 564 A in solutions containing 1 mM Tris buffer plus 0.2-1 mM NaCl and decrease still further to 440 A in solutions containing 0.2 mM Mg(2+) ions. The persistence length exhibits a shallow maximum at 20 degrees C and decreases slowly upon either increasing or decreasing the temperature, regardless of the model used to fit the data. By contrast, the consensus values of the hydrodynamic radius are independent of temperature. The calculated persistence lengths and hydrodynamic radii are compared with other data in the literature.     

Modulation of intramolecular interactions in superhelical DNA by curved sequences: a Monte Carlo simulation study.

A Monte Carlo model for the generation of superhelical DNA structures at thermodynamic equilibrium (Klenin et al., 1991; Vologodskii et al., 1992) was modified to account for the presence of local curvature. Equilibrium ensembles of a 2700-bp DNA chain at linking number difference delta Lk = -15 were generated, with one or two permanent bends up to 120 degrees inserted at different positions. The computed structures were then analyzed with respect to the number and positions of the end loops of the interwound superhelix, and the intramolecular interaction probability of different segments of the DNA. We find that the superhelix structure is strongly organized by permanent bends. A DNA segment with a 30 degrees bend already has a significantly higher probability of being at the apex of a superhelix than the control, and for a 120 degrees bend the majority of DNAs have one end loop at the position of the bend. The entropy change due to the localization of a 120 permanent bend in the end loop is estimated to be -17 kJ mol-1 K-1. When two bends are inserted, the conformation of the superhelix is found to be strongly dependent on their relative positions: the straight interwound form dominates when the two bends are separated by 50% of the total DNA length, whereas the majority of the superhelices are in a branched conformation when the bends are separated by 33%. DNA segments in the vicinity of the permanent bend are strongly oriented with respect to each other.     

The gamma delta resolvase bends the res site into a recombinogenic complex.

We have characterized complexes between the gamma delta resolvase and its recombination site, res, using both a gel retardation assay and DNase I cleavage. The mobility of resolvase-res complexes in polyacrylamide gels is sensitive to the location of res within the DNA fragment and is at a minimum when res is at its center. This behavior is characteristic of a protein-dependent bend. By the same assay we have found that bends are induced upon the binding of resolvase to each of the three individual binding sites that constitute res. In the wild-type res, the centers of binding sites I and II are 53 bp apart and the central section of the intersite DNA is sensitive to DNase I cleavage. We find that insertions of 10 or 21 bp (one or two turns of the DNA helix) have no discernible effect on the ability of res to recombine or to form complexes with resolvase. However, insertions of short segment (e.g. 6 or 17 bp) equivalent to nonintegral numbers of helical turns, inhibit recombination and prevent the formation of the normally compact resolvase-res complex. Complexes of resolvase with res containing 10 or 21 bp insertions exhibit a pattern of enhanced and suppressed DNase I cleavages that suggest that the intersite segment is curved. This curvature requires both that site I and II are appropriately spaced, and that site III is also present and occupied.     

Chromatin reconstitution on small DNA rings. II. DNA supercoiling on the nucleosome

DNA supercoiling on the nucleosome was investigated by relaxing with topoisomerase I mono- and dinucleosomes reconstituted on small DNA rings. Besides 359 base-pair (bp) rings whose linking differences were integers, two additional series of rings with fractional differences, 341 and 354 bp in size, were used. Mononucleosomes reconstituted on 359 bp rings were found to relax into a single mononucleosome form. In contrast, 341 and 354 bp mononucleosomes relaxed into a mixture of two forms, corresponding to two adjacent topoisomers. The observation that the ratio between these two forms was, within each ring series, virtually independent of the initial linking number of the topoisomer used for the reconstitution suggested that each partition reflected an equilibrium. Comparison with the equilibria observed for the same rings in the absence of histones showed that the formation of a single nucleosome is associated with a linking number change of -1.1(+/-0.1) turn. Dinucleosomes, in contrast, were not relaxed to completion and do not reach equilibria. The corresponding linking number change per nucleosome was, however, estimated to be similar to the above figure, in agreement with previous data from the literature obtained with circular chromatins containing larger numbers of nucleosomes. DNA structure in mononucleosomes was subsequently investigated by means of high-resolution electron microscopy and gel electrophoresis. It was found that the above linking number reduction could be ascribed to a particle with a large open extranucleosomal DNA loop and with no more than 1.5 turns of a superhelix around the histone core. A theoretical model of a nucleosome on a small ring was constructed in which one part of the DNA was wrapped around a cylinder and the other part was free to vary both in torsion and flexion. The linking number reduction predicted was found to be most consistent with experimental data when the twist of the DNA in the superhelix was between 10.5 and 10.65 pb per turn, suggesting that wrapping on the nucleosome does not alter the twist of the DNA significantly. A lower estimate of the linking number reduction associated with a two-turn nucleosome was also derived, based on an analysis of recent data obtained upon treatment of reconstituted minichromosomes with gyrase. The value, 1.6 turns, set a lower limit of 10.44 bp per turn for the twist of nucleosomal DNA, in agreement with the above estimate.(ABSTRACT TRUNCATED AT 400 WORDS)     

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816–3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed ‘quasi-regularly’ throughout the genome, including both the coding and intergenic sequences. Second, there is a 10–12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in ~100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the ~100 bp long clusters of the phased A-tracts constitute the ‘structural code’ for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins.     

Molecular modeling of closed circular DNA thermodynamic ensembles

Many modeling studies of supercoiled DNA are based on equilibrium structures from theoretical calculations or energy minimization. Since closed circular DNAs are flexible, it is possible that errors are introduced by calculating properties from a single minimum energy structure, rather than from a complete thermodynamic ensemble. We have investigated this question using molecular dynamics simulations on a low resolution molecular mechanics model in which each base pair is represented by three points (a plane). This allows the inclusion of sequence-dependent variations of tip, inclination, and twist. Three kinds of sequences were tested: (1) homogeneous DNA, in which all base pairs have the helicoidal parameters of an ideal, average B-DNA; (2) random sequence DNA; and (3) curved DNA. We examined the rate of convergence of various structural parameters. Convergence for most of these is slowest for homogeneous sequences, more rapid for random sequences, and most rapid for curved sequences. The most slowly converging parameter is the antipodes profile. In a plasmid with N base pairs (bp), the antipodes distance is the distance d _ij from base pair i to base pair j halfway around the plasmid, j = i + N/2. The antipodes profile at time t is a plot of d _ij over the range i = 1, N/2. In a homogeneous plasmid, convergence requires that the antipodes profile averaged over time must be flat. Even in the small plasmids examined here, the average properties of the ensembles were found to differ from those of static equilibrium structures. These effects will be even more dramatic for larger plasmids. Further, average and dynamic properties are affected by both plasmid size and sequence. © 1996 John Wiley & Sons, Inc.     

Topological measurement of an A-tract bend angle: variation of duplex winding

The rotational variant method of Lutter et al. was developed to measure the bend angle induced when a protein binds to DNA. To measure the intrinsic bend conferred by a sequence of six adenine bases (an A6 tract), the method was modified by relaxing at high temperature to remove the bend. We describe here an alternative approach that involves unwinding the duplex DNA between adjacent bends in plasmids containing tandemly repeated blocks of A-tracts. This method measures the topological difference contributed by adjacent bends when they are in two different rotational settings, and therefore does not require reference to a straight state. The interbend DNA was unwound by use of the intercalator chloroquine, or, alternatively, by raising the temperature in the relaxation reaction. The effect of this unwinding is to change the pitch of the superhelix of the tandem repeats from which the bend angle is measured. The result is a bend angle value that is consistent with that measured using the bend-straightening version of the method. This version offers several advantages that complement the conventional bent versus straight approach.     

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.     

Total weak acid concentration and effective dissociation constant of nonvolatile buffers in human plasma.

The strong ion approach provides a quantitative physicochemical method for describing the mechanism for an acid-base disturbance. The approach requires species-specific values for the total concentration of plasma nonvolatile buffers (A(tot)) and the effective dissociation constant for plasma nonvolatile buffers (K(a)), but these values have not been determined for human plasma. Accordingly, the purpose of this study was to calculate accurate A(tot) and K(a) values using data obtained from in vitro strong ion titration and CO(2) tonometry. The calculated values for A(tot) (24.1 mmol/l) and K(a) (1.05 x 10(-7)) were significantly (P < 0.05) different from the experimentally determined values for horse plasma and differed from the empirically assumed values for human plasma (A(tot) = 19.0 meq/l and K(a) = 3.0 x 10(-7)). The derivatives of pH with respect to the three independent variables [strong ion difference (SID), PCO(2), and A(tot)] of the strong ion approach were calculated as follows: dpH/dSID(+) = [1 + 10(pK(a)-pH)](2)/(2.303 x [SPCO(2)10(pH-pK'(1)[1 + 10(pK(a)-pH](2) + A(tot)10(pK(a)-PH]]; dpH/dPCO(2) = S10(-pK'(1)/[2.303[A(tot)10(pH)(10(pH + 10(pK(a))(-2) - SID(+)10(-pH)]], dpH/dA(tot) = -1/[2.303[SPCO(2)10(pH-pK'(1) + SID(+)10(pK(a)-pH)]], where S is solubility of CO(2) in plasma. The derivatives provide a useful method for calculating the effect of independent changes in SID(+), PCO(2), and A(tot) on plasma pH. The calculated values for A(tot) and K(a) should facilitate application of the strong ion approach to acid-base disturbances in humans.     

Influence of the OGG1 Ser326Cys polymorphism on oxidatively damaged DNA and repair activity

Oxidatively damaged DNA base lesions are considered to be mainly repaired by 8-oxoguanine DNA glycosylase (OGG1) mediated pathways. We investigated the effect of the OGG1 Ser326Cys polymorphism on the level and repair of oxidatively damaged DNA in mononuclear blood cells (MNBC) by means of the comet assay. We collected blood samples from 1,019 healthy subjects and genotyped for the OGG1 Ser326Cys polymorphism. We found 49 subjects homozygous for the variant genotype (Cys/Cys) and selected same numbers of age-matched subjects with the heterozygous (Ser/Cys) and homozygous wild-type genotype (Ser/Ser). Carriers of the Cys/Cys genotype had higher levels of formamidopyrimidine DNA glycosylase (FPG) sensitive sites in MNBC (0.31 ± 0.03 lesions/10(6)bp) compared to Ser/Ser (0.19 ± 0.02 lesions/10(6)bp, P<0.01). The level of hOGG1 sensitive sites in MNBC from the Ser326Cys carriers (0.19 ± 0.16 lesions/10(6) bp) was also higher compared to the Ser/Ser genotype (0.11 ± 0.09 lesions/10(6) bp, P<0.05). Still, there was no genotype-related difference in DNA repair incision activity of MNBC extracts on nucleoids with oxidatively damaged DNA induced by Ro19-8022/white light (P=0.20). In addition, there were no differences in the expression of OGG1 (P=0.69), ERCC1 (P=0.62), MUTYH (P=0.85), NEIL1 (P=0.17) or NUDT1 (P=0.48) in whole blood. Our results indicate that the OGG1 Ser326Cys polymorphism has limited influence on the DNA repair incisions by extracts of MNBC, whereas the apparent increased risk of cancer in subjects with the Cys/Cys genotype may be because of higher levels of oxidatively damaged DNA.     

Ectopic Expression of the Drosophila Homeotic Gene Proboscipedia under Antennapedia P1 Control Causes Dominant Thoracic Defects

A deletion mutation in the Antennapedia Complex of Drosophila melanogaster, Df(3R)SCBXL2, induces both dominant and recessive loss-of-function phenotypes. The deletion is associated with diminished function of proboscipedia (pb), a homeotic gene required for mouthparts formation. Df(3R)SCBXL2 also has associated dominant thoracic defects related to diminished expression of the homeotic Antennapedia (Antp) gene copy on the homologous chromosome. This is shown to be a consequence of ectopic pb expression in the thorax. Newly juxtaposed Antp sequences provide the pb gene on the deletion bearing chromosome with a second promoter, Antp P1, in addition to its own. Ectopic pb protein expression occurs under Antp P1 control, by alternate splicing, and results in diminished accumulation of Antp protein in the imaginal disc cells where Antp P1 is normally expressed. The analysis of this mutant chromosome thus demonstrates that pb protein is capable of participating in the negative regulation of a more posteriorly expressed homeotic gene, as well as serving a homeotic "selector" function in the head.     

The histone core exerts a dominant constraint on the structure of DNA in a nucleosome.

We have examined the structures of unique sequence, A/T-rich DNAs that are predicted to be relatively rigid [oligo(dA).oligo(dT)], flexible [oligo[d(A-T)]], and curved, using the hydroxyl radical as a cleavage reagent. A 50-base-pair segment containing each of these distinct DNA sequences was placed adjacent to the T7 RNA polymerase promoter, a sequence that will strongly position nucleosomes. The final length of the DNA fragments was 142 bp, enough DNA to assemble a single nucleosome. Cleavage of DNA in solution, while bound to a calcium phosphate crystal, and after incorporation into a nucleosome is examined. We find that the distinct A/T-rich DNAs have very different structural features in solution and helical periodicities when bound to a calcium phosphate. In contrast, the organization of the different DNA sequences when associated with a histone octamer is very similar. We conclude that the histone core exerts a dominant constraint on the structure of DNA in a nucleosome and that inclusion of these various unique sequences has only a very small effect on overall nucleosome stability and structure.     

Intrinsic curvature in duplex DNA inhibits Human Topoisomerase I

Human Topoisomerase I (hTopo I) have been known as a potential target for cancer therapy. A series of duplex DNA with different intrinsic curvatures have been designed as inhibitors to hTopo I. The activities of hTopo I on relaxing supercoiled plasmid pUC 19 are apparently diminished in the presence of the curved DNA. More potent inhibitions and smaller IC(50) are achieved by duplex DNA with higher curvatures. EMSA indicates that hTopo I can recognize the curved DNA through binding interactions. Our studies demonstrate that the activity of hTopo I can be modulated by the intrinsic curvature of linear DNA and provide a new avenue to design curved DNA as hTopo I inhibitors with high therapeutic efficiency and low toxicity.     

Intrinsic bends and integration host factor binding at F plasmid oriT.

F plasmid oriT DNA extending from the F kilobase coordinate 66.7 (base pair [bp] 1 on the oriT sequence map) rightward to bp 527 was analyzed for intrinsic bends (by permutation assays) and for binding of integration host factor (IHF) (by gel retardation and DNase footprinting). Intrinsic bending of the 527-bp fragment (bend center approximately at bp 240) was represented as a composite of at least two components located near bp 170 and near bp 260. IHF bound primarily to a site extending from bp 165 to 195 and with lower affinity to a site extending from bp 287 to 319. The intrinsic curvature and sequences to which IHF binds (IHF is known to bend DNA) may play a structural role in oriT function.     

Monte Carlo simulations of supercoiling free energies for unknotted and trefoil knotted DNAs.

A new Monte Carlo (MC) algorithm is proposed for simulating inextensible circular chains with finite twisting and bending rigidity. This new algorithm samples the relevant Riemann volume elements in a uniform manner, when the constraining potential vanishes. Simulations are performed for filaments comprising 170 subunits, each containing approximately 28 bp, which corresponds to a DNA length of 4770 bp. The bending rigidity is chosen to yield a persistence length, P = 500 A, and the intersubunit potential is taken to be a hard-cylinder potential with diameter d = 50 A. This value of d yields the same second virial coefficient as the electrostatic potential obtained by numerical solution of the Poisson-Boltzmann equation for 150 mM salt. Simulations are performed for unknotted circles and also for trefoil knotted circles using two different values of the torsional rigidity, C = (2.0 and 3.0) x 10(-19) dyne cm2. In the case of unknotted circles, the simulated supercoiling free energy varies practically quadratically with linking difference delta l. The simulated twist energy parameter ET, its slope dET/dT, and the mean reduced writhe <w>/delta l for C = 3 x 10(-19) dyne cm2 all agree well with recent simulations for unknotted circles using the polygon-folding algorithm with identical P, d, and C. The simulated ET vs. delta l data for C = 2.0 x 10(-19) dyne cm2 agree rather well with recent experimental data for p30 delta DNA (4752 bp), for which the torsional rigidity, C = 2.07 x 10(-19) dyne cm2, was independently measured. The experimental data for p30 delta are enormously more likely to have arisen from C = 2.0 x 10(-19) than from C = 3.0 x 10(-19) dyne cm2. Serious problems with the reported experimental assessments of ET for pBR322 and their comparison with simulated data are noted. In the case of a trefoil knotted DNA, the simulated value, (ET)tre, exceeds that of the unknotted DNA, (ET)unk, by approximately equal to 1.40-fold at magnitude of delta l = 1.0, but declines to a plateau about 1.09-fold larger than (ET)unk when magnitude of delta l > or = 15. Although the predicted ratio, (ET)tre/(ET)unk approximately equal to 1.40, agrees fairly well with recent experimental measurements on a 5600-bp DNA, the individual measured ET values, like some of those reported for pBR322, are so large that they cannot be simulated using P = 500 A, d = 50 A, and any previous experimental estimate of C.     

单纯疱疹病毒的PCR直接基因分型

本文利用ＰＣＲ技术建立一种对ＨＳＶ直接基因分型的方法。在ＨＳＶ－Ⅰ、Ⅱ两型的ＤＮＡ多聚酶基因上设计了一条两型共同的上游引物（ＨＤＰ－Ｂ）和两条型特异的下游引物（ＨＤＰ－１、ＨＤＰ－２）。三条引物共同组成一个扩增反应体系,在ＨＳＶ－Ⅰ产生５４３ｂｐ条带,ＨＳＶ－Ⅱ产生３７２ｂｐ条带,据此在基因水平上对ＨＳＶ进行分型。５株不同来源的ＨＳＶ（２株Ⅰ型,３株Ⅱ型）分型结果与病毒分离及血清学方法完全一致。该反应体系与其它来源的ＤＮＡ不产生特异反应,敏感性可达１ｆｇ。应用该法对１５１份临床可疑ＨＳＶ感染的标本进行检测并分型,结果与免疫学方法完全一致。     

Monovalent cation binding by curved DNA molecules containing variable numbers of a-tracts

Monovalent cation binding by DNA A-tracts, runs of four or more contiguous adenine or thymine residues, has been determined for two curved ∼200 basepair (bp) restriction fragments, one taken from the M13 origin of replication and the other from the VP1 gene of SV40. These two fragments have previously been shown to contain stable, centrally located bends of 44° and 46°, respectively, located within ∼60 bp “curvature modules” containing four or five irregularly spaced A-tracts. Transient electric birefringence measurements of these two fragments, sequence variants containing reduced numbers of A-tracts in the SV40 curvature module or changes in the residues flanking the A-tracts in the M13 curvature module, have been combined with the free solution electrophoretic mobilities of the same fragments using known equations to estimate the effective charge of each fragment. The effective charge is reduced, on average, by one-third charge for each A-tract in the curvature module, suggesting that each A-tract binds a monovalent cation approximately one-third of the time. Monovalent cation binding to two or more A-tracts is required to observe significant curvature of the DNA helix axis.     

Orientation of the central pair complex during flagellar bend formation in Chlamydomonas

Thin section electron micrographs of rapidly fixed Chlamydomonas cells were used to establish a relationship between flagellar bends and orientation of the central pair microtubule complex. Using conditions that preserve flagellar waveforms during both forward swimming (asymmetric bends) and backward swimming (symmetric bends), we found that central pair orientation differs in bent regions and straight regions. During forward swimming, a plane through the two central pair microtubules is parallel to the bend plane throughout principal bends, in both effective stroke and recovery stroke phases of the beat cycle. In these curved segments, the C1 microtubule always faces the outer edge of the curve. This parallel orientation twists in straight regions both proximal and distal to bends. During backward swimming episodes induced by photoshock, when Chlamydomonas flagella beat with principal and reverse bends of similar magnitude, the central pair twists by 180 degrees between successive bends. These observations support a model in which central pair orientation in Chlamydomonas is linked to doublet-specific dynein activation, and bend propagation is linked to rotation of the central pair complex.