Crenarchaeal CdvA forms double-helical filaments containing DNA and interacts with ESCRT-III-like CdvB

Background

The phylum Crenarchaeota lacks the FtsZ cell division hallmark of bacteria and employs instead Cdv proteins. While CdvB and CdvC are homologues of the eukaryotic ESCRT-III and Vps4 proteins, implicated in membrane fission processes during multivesicular body biogenesis, cytokinesis and budding of some enveloped viruses, little is known about the structure and function of CdvA. Here, we report the biochemical and biophysical characterization of the three Cdv proteins from the hyperthermophilic archaeon Metallospherae sedula.

Methodology/Principal Findings

Using sucrose density gradient ultracentrifugation and negative staining electron microscopy, we evidenced for the first time that CdvA forms polymers in association with DNA, similar to known bacterial DNA partitioning proteins. We also observed that, in contrast to full-lengh CdvB that was purified as a monodisperse protein, the C-terminally deleted CdvB construct forms filamentous polymers, a phenomenon previously observed with eukaryotic ESCRT-III proteins. Based on size exclusion chromatography data combined with detection by multi-angle laser light scattering analysis, we demonstrated that CdvC assembles, in a nucleotide-independent way, as homopolymers resembling dodecamers and endowed with ATPase activity in vitro. The interactions between these putative cell division partners were further explored. Thus, besides confirming the previous observations that CdvB interacts with both CdvA and CdvC, our data demonstrate that CdvA/CdvB and CdvC/CdvB interactions are not mutually exclusive.

Conclusions/Significance

Our data reinforce the concept that Cdv proteins are closely related to the eukaryotic ESCRT-III counterparts and suggest that the organization of the ESCRT-III machinery at the Crenarchaeal cell division septum is organized by CdvA an ancient cytoskeleton protein that might help to coordinate genome segregation.     

Interaction of polyamines with chromatin and DNA: formation of compact structures

We have used flow linear dichroism (LD) and light scattering at 90 degrees to study the condensation of both DNA and calf thymus chromatin induced by spermine, triamines NH3+(CH2)iNH+(CH2)jNH3+, designated as much less than i, j much greater than: much less than 3, 4 much greater than (spermidine), much less than 3, 3 much greater than, much less than 2, 3 much greater than, much less than 2, 2 much greater than; the diamines putrescine and cadaverine and MgCl2. It is found that the different polyamines affected DNA and chromatin in a similar way. The degree of compaction of the chromatin fibers induced by spermine, triamines except much less than 2, 2 much greater than and Mg2+ has been found to be identical. The triamine much less than 2, 2 much greater than and the diamines studied do not condense either chromatin of DNA. Such a big difference in the action of the triamines indicates that not only the charge, but also the structure of the polycations are important for their interactions with DNA and chromatin. The stoichiometry of polyamine binding to chromatin at which condensation occurred is found to be 2 polyamine molecules per DNA helical turn. Polyamines are supposed to bind to the exposed sites of core DNA every 10 b.p. The extent of DNA phosphate neutralization by the histones is estimated to be about 55%. It has been shown that a mixture of mono- and multivalent cations affected DNA and chromatin condensation competitively and not synergistically, as claimed in a recent report by Sen and Crothers.     

Hydration of the phosphate group in double-helical DNA.

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.     

Binding of alkaline cations to the double-helical form of gramicidin.

Gramicidin is a polypeptide antibiotic that forms monovalent cation-specific channels in membrane environments. In organic solvents and in lipids containing unsaturated fatty acid chains, it forms a double-helical "pore" structure, in which two monomers are intertwined. This form of gramicidin can bind two cations inside its lumen, and the crystal structures of both an ion complex and an ion-free form have been determined. In this study, we have used circular dichroism (CD) spectroscopy to examine the binding mechanism and the binding constants (K1 and K2) of cations to gramicidin in the double helical form in methanol solution. The dramatic change in optical rotation in the far-ultraviolet CD spectrum of gramicidin provides a useful tool for monitoring the binding. The binding mechanism appears to involve a large conformation change associated with the binding of ions to the first of the two sites. The calculated values for the K1 binding constants for alkaline cations are considerably smaller than the K2 binding constants. The order of binding affinity for alkaline cations is similar to that for the helical dimer "channel" form of gramicidin, i.e., Cs+ approximately Rb+ > > K+ > Li+, but in comparison to the helical dimer form, the binding to double-helical dimers is dominated by a cation size-dependent conformational change in the gramicidin structure.     

Natively unfolded regions of the vertebrate fibrinogen molecule

Although it has long been realized that a large portion of the fibrinogen alpha chain has little if any defined structure, the physiological significance of this flexible appendage remains mysterious.     

Populating partially unfolded forms by hydrogen exchange-directed protein engineering

The native-state hydrogen exchange of a redesigned apocytochrome b(562) suggests that at least two partially unfolded forms (PUFs) exist for this four-helix bundle protein under native conditions. The more stable PUF has the N-terminal helix unfolded. To verify the conclusion further and obtain more detailed structural information about this PUF, five hydrophobic core residues in the N-terminal helix were mutated to Gly and Asp to destabilize the native state selectively and populate the PUF for structural studies. The secondary structure and the backbone dynamics of this mutant were characterized using multidimensional NMR. Consistent with the prediction, the N-terminal region of the mutant was found to be unfolded while other parts of the proteins remained folded. These results suggest that native-state hydrogen exchange-directed protein engineering can be a useful approach to populating partially unfolded forms for detailed structural studies.     

Bend of the furanose ring of the double-helical DNA

The interpretation of conformationally dependent double-helical DNA Raman bands, which stipulates their appearance with the furanose ring-puckering vibration is suggested. Quantum-mechanical calculations of the furanose ring-puckering vibration occurring by a pseudorotation pathway are carried out. Vibrational states eigen values and eigen functions are calculated for a double-well pseudorotation potential of the furanose. The obtained value of the pseudorotation barrier 4.2-3.6 kcal/mol proves to be in a good agreement with other theoretical and experimental estimates of the barrier.     

Interdependence of conformational variables in double-helical DNA.

DNA exhibits conformational polymorphism, with the details depending on the sequence and its environment. To understand the mechanisms of conformational polymorphism and these transitions, we examine the interrelationships among the various conformational variables of DNA. In particular, we examine the stress-strain relation among conformational variables, describing base-pair morphology and their effects on the backbone conformation. For the calculation of base pairs, we use the method previously developed to calculate averages over conformational variables of DNA. Here we apply this method to calculate the Boltzmann averages of conformational variables for fixed values of one particular conformational variable, which reflects the strain in the structure responding to a particular driving stress. This averaging over all but one driving variable smooths the usual rough energy surface to permit observation of the effects of one conformational variable at a time. The stress-strain analyses of conformational variables of base pair slide, twist, and roll, which exhibit characteristic changes during the conformational transition of DNA, have shown that the conformational changes of base pairs are strongly correlated with one another. Furthermore, the stress-strain relations are not symmetrical with respect to these variables, i.e., the response of one coordinate to another is different from the reverse direction. We also examine the effect of conformational changes in base-pair variables on the sugar-backbone conformation by using the minimization method we developed. The conformational changes of base pairs affect the sugar pucker and other dihedral angles of the backbone of DNA, but each variable affects the sugar-backbone differently. In particular, twist is found to have the most influence in affecting the sugar pucker and backbone conformation. These calculated conformational changes in base pairs and backbone segments are consistent with experimental observations and serve to validate the calculation method.     

A structural similarity analysis of double-helical DNA

A database of the structural properties of all 32,896 unique DNA octamer sequences has been calculated, including information on stability, the minimum energy conformation and flexibility. The contents of the database have been analysed using a variety of Euclidean distance similarity measures. A global comparison of sequence similarity with structural similarity shows that the structural properties of DNA are much less diverse than the sequences, and that DNA sequence space is larger and more diverse than DNA structure space. Thus, there are many very different sequences that have very similar structural properties, and this may be useful for identifying DNA motifs that have similar functional properties that are not apparent from the sequences. On the other hand, there are also small numbers of almost identical sequences that have very different structural properties, and these could give rise to false-positives in methods used to identify function based on sequence alignment. A simple validation test demonstrates that structural similarity can differentiate between promoter and non-promoter DNA. Combining structural and sequence similarity improves promoter recall beyond that possible using either similarity measure alone, demonstrating that there is indeed information available in the structure of double-helical DNA that is not readily apparent from the sequence.     

Nonenymatic ligation of double-helical DNA by alternate-strand triple helix formation.

Nonenzymatic ligation of double-stranded DNA has been performed using an alternate-strand binding oligodeoxyribonucleotide template to juxtapose the duplex termini in a triple helical complex. The template associates with the duplex termini by Hoogsteen hydrogen bonding to alternate strands on opposite sides of the ligation site. Intermolecular and intramolecular ligation of linearized plasmid DNA are observed in the reaction, which depends on the template oligodeoxyribonucleotide and a condensing agent, N-cyanoimidazole. Intramolecular ligation products include those in which both strands are covalently closed in a circle. Ligation of the two strands is sequential and occurs at comparable rates for the first and second strands ligating. The covalent linkages formed in the reaction can be cleaved by the restriction endonuclease Stu I, supporting their identification as phosphodiesters.     

Calorimetric determination of base-stacking enthalpies in double-helical DNA molecules

Differential scanning calorimetry was used to directly determine the transition enthalpies accompanying the duplex-to-single-strand transition of poly[d(AT)], poly(dA)·poly(dT), poly[d(AC)]·poly[d(TG)], and d(GCGCGC). The calorimetric data allow us to define the following average base-stacking enthalpies:

Human glycine N-methyltransferase is unfolded by urea through a compact monomer state

Human recombinant glycine N-methyltransferase (GNMT) unfolding by urea was studied by enzyme activity, size-exclusion chromatography, fluorescence spectroscopy, and circular dichroism. Urea unfolding of GNMT is a two-step process. The first transition is a reversible dissociation of the GNMT tetramer to compact monomers in 1.0-3.5M urea with enzyme inactivation. The compact monomers were characterized by Stokes radius (R(s)) of 40.7A equal to that of globular proteins with the same molecular mass as GNMT monomers, absence of exposure of tryptophan residues into solvent, and presence of about 50% of secondary structure of native protein. The second step of GNMT unfolding is a reversible transition of monomers from compact to a fully unfolded state with R(s) of 50A, exposed tryptophan residues, and disrupted secondary structure in 8M urea.     

Superhelicity of nucleosomal DNA changes its double-helical repeat

Winding DNA in a superhelix can be considered a process consisting of two smooth deformations: bending and twisting. The extra twist angle introduced by winding DNA into the nucleosomal superhelix is calculated by means of the Crick formula to be -0.5 degrees per base pair (bp). This is equivalent to a change of -0.15 +/- 0.015 bp in the DNA double-helical repeat. Free DNA in solution is known to have a helical repeat of 10.55 +/- 0.1 bp. On the other hand, a weighted average of various estimates of the DNA repeat in the nucleosome is 10.38 +/- 0.02. The difference happens to be perfectly accounted for by the superhelicity of the nucleosomal DNA. This implies that the latter is essentially nonconstrained .     

Superhelicity of nucleosomal DNA changes its double-helical repeat

Winding DNA in a superhelix can be considered a process consisting of two smooth deformations: bending and twisting. The extra twist angle introduced by winding DNA into the nucleosomal superhelix is calculated by means of the Crick formula to be −0.5° per base pair (bp). This is equivalent to a change of −0.15±0.015 bp in the DNA double-helical repeat. Free DNA in solution is known to have a helical repeat of 10.55±0.1 bp. On the other hand, a weighted average of various estimates of the DNA repeat in the nucleosome is 10.38±0.02. The difference happens to be perfectly accounted for by the superhelicity of the nucleosomal DNA. This implies that the latter is essentially nonconstrained.     

Effect of 1-methyladenine on double-helical DNA structures

Methylation at the N1 site of adenine leads to the formation of cytotoxic 1-methyladenine (m1A). Since the N1 site of adenine is involved in the hydrogen bonding of T·A and A·T Watson–Crick base pairs, it is expected that the pairing interactions will be disrupted upon 1-methylation. In this study, high-resolution NMR investigations were performed to determine the effect of m1A on double-helical DNA structures. Interestingly, instead of disrupting hydrogen bonding, we found that 1-methylation altered the T·A Watson–Crick base pair to T(anti)·m1A(syn) Hoogsteen base pair, providing insights into the observed differences in AlkB-repair efficiency between dsDNA and ssDNA.     

Distance-dependent dielectric constants and their application to double-helical DNA

Detailed studies of structures of biological macromolecules, even in simplified models, involve many costly and time-consuming calculations. Any thorough methods require sampling of an extremely large conformation and momentum space. Calculations of electrostatic interactions, which depend on many physical factors, such as the details of solvent, solvent accessibility in macromolecules, and molecular polarizability, are always developed in a compromise between more rigorous, detailed models and the need for immediate application to complicated biological systems. In this paper, a middle ground is taken between the more exact theoretical models and the simplest constant values for the dielectric constant. The effects of solvent, counterions, and molecular polarizability are incorporated through a set of adjustable parameters that should be determined from experimental conditions. Several previous forms for the dielectric function are compared with the new ones. The present methods use Langevin functions to span the region of dielectric constant between bulk solvent and cavity values. Application of such dielectric models to double-helical DNA is important because base-stacking preferences were previously demonstrated [A. Sarai, J. Mazur, R. Nussinov, and R. L. Jernigan (1988) Biochemistry, vol. 27, pp. 8498-8502] to be sensitive to the electrostatic formulation. Here we find that poly(dG).poly(dC) can be A form for high screening and B form for low screening. By contrast, poly(dA).poly(dT) can only take the B form. Base stacking is more sensitive to the form of the dielectric function than are the sugar-phosphate backbone conformations. Also in B form, the backbone conformations are not so affected by the base types as in A form.     

DNA and nuclear aggregates of polyamines

Polyamines (PAs) are linear polycations that are involved in many biological functions. Putrescine, spermidine and spermine are highly represented in the nucleus of eukaryotic cells and have been the subject of decades of extensive research. Nevertheless, their capability to modulate the structure and functions of DNA has not been fully elucidated. We found that polyamines self-assemble with phosphate ions in the cell nucleus and generate three forms of compounds referred to as Nuclear Aggregates of Polyamines (NAPs), which interact with genomic DNA. In an in vitro setting that mimics the nuclear environment, the assembly of PAs occurs within well-defined ratios, independent of the presence of the DNA template. Strict structural and functional analogies exist between the in vitro NAPs (ivNAPs) and their cellular homologues. Atomic force microscopy showed that ivNAPs, as theoretically predicted, have a cyclic structure, and in the presence of DNA, they form a tube-like arrangement around the double helix. Features of the interaction between ivNAPs and genomic DNA provide evidence for the decisive role of "natural" NAPs in regulating important aspects of DNA physiology, such as conformation, protection and packaging, thus suggesting a new vision of the functions that PAs accomplish in the cell nucleus.     

The unfolded state of NTL9 is compact in the absence of denaturant

Interest in the unfolded state of proteins has grown with the realization that this state can have considerable structure in the absence of denaturants. Natively unfolded proteins, mutations that unfold proteins under native conditions, and changes in pH that induce unfolding are attractive models for the unfolded state in the absence of denaturant. The unfolded state of the N-terminal domain of ribosomal protein L9 (NTL9) was previously shown to contain significant non-native electrostatic interactions [Cho, J. H., Sato, S., and Raleigh, D. P. (2004) J. Mol. Biol. 338, 827-837]. NTL9 has a mixed alpha-beta structure and folds via a two-state mechanism. We have generated a model of the unfolded state of NTL9 in the absence of denaturant by substitution of an alanine for phenylalanine 5 located in the core of this protein. The CD spectrum of the variant, denoted as F5A, exhibits significantly less structure than the wild type; however, the mean residue ellipticity of F5A at 222 nm (-8200 deg cm(2) dmol(-)(1)) is considerably larger than expected for a fully unfolded protein, indicating that residual secondary structure is populated. F5A also has more residual structure than the urea-unfolded wild type. The stability of F5A is estimated to be at least 1 kcal/mol unfavorable, showing that the unfolded state is populated to 84% or more. NMR pulsed-field gradient measurements yield a hydrodynamic radius of 16.1 A for wild-type NTL9 and 20.8 A for the F5A variant in native buffer. The physiologically relevant unfolded state of wild-type NTL9 is likely to be even more compact than F5A since the mutation should reduce the level of hydrophobic clustering in the unfolded state in the absence of denaturant. The hydrodynamic radius of F5A increases to 25.9 A in 8 M urea, and a value of 23.5 A is obtained for the wild type under similar conditions. The results show that the unfolded state of F5A in the absence of denaturant is more compact and contains more structure than the urea-unfolded form.     

Simulation of the conformation and dynamics of a double-helical model for DNA

We propose a partially flexible, double-helical model for describing the conformational and dynamic properties of DNA. In this model, each nucleotide is represented by one element (bead), and the known geometrical features of the double helix are incorporated in the equilibrium conformation. Each bead is connected to a few neighbor beads in both strands by means of stiff springs that maintain the connectivity but still allow for some extent of flexibility and internal motion. We have used Brownian dynamics simulation to sample the conformational space and monitor the overall and internal dynamics of short DNA pieces, with up to 20 basepairs. From Brownian trajectories, we calculate the dimensions of the helix and estimate its persistence length. We obtain translational diffusion coefficient and various rotational relaxation times, including both overall rotation and internal motion. Although we have not carried out a detailed parameterization of the model, the calculated properties agree rather well with experimental data available for those oligomers.