Structural and thermodynamic principles of viral packaging

Packaging of genetic material inside a capsid is one of the major processes in the lifecycle of bacteriophages. To establish the basic principles of packing double-stranded DNA into a phage, we present a low-resolution model of bacteriophage varphi29 and report simulations of DNA packaging. The simulations show excellent agreement with available experimental data, including the forces of packaging and the average structures seen in cryo-electron microscopy. The conformation of DNA inside the bacteriophage is primarily determined by the shape of the capsid and the elastic properties of DNA, but the energetics of packaging are dominated by electrostatic repulsions and the large entropic penalty associated with DNA confinement. In this slightly elongated capsid, the DNA assumes a folded toroidal conformation, rather than a coaxial spool. The model can be used to study packaging of other bacteriophages with different shapes under a range of environmental conditions.     

Double-stranded DNA organization in bacteriophage heads: an alternative toroid-based model.

Studies of the organization of double-stranded DNA within bacteriophage heads during the past four decades have produced a wealth of data. However, despite the presentation of numerous models, the true organization of DNA within phage heads remains unresolved. The observations of toroidal DNA structures in electron micrographs of phage lysates have long been cited as support for the organization of DNA in a spool-like fashion. This particular model, like all other models, has not been found to be consistent will all available data. Recently we proposed that DNA within toroidal condensates produced in vitro is organized in a manner significantly different from that suggested by the spool model. This new toroid model has allowed the development of an alternative model for DNA organization within bacteriophage heads that is consistent with a wide range of biophysical data. Here we propose that bacteriophage DNA is packaged in a toroid that is folded into a highly compact structure.     

Modulating DNA configuration by interfacial traction: an elastic rod model to characterize DNA folding and unfolding

As a continuum model of DNA, a thin elastic rod subjected to interfacial interactions is used to investigate the equilibrium configuration of DNA in intracellular solution. The interfacial traction between the rod and the solution environment is derived in detail. Kirchhoff’s theory of elastic rods is used to analyze the equilibrium configuration of a DNA segment under the action of the interfacial traction. The influences of the interfacial energy factor and bending stiffness on the toroidal spool formation of the DNA segment are discussed. The results show that the equilibrium configuration of DNA is mainly determined by competition between the interfacial energy and elastic strain energy of the DNA itself, and the interfacial traction is one of the forces that drives DNA folding and unfolding.     

Evaluation of elastic properties of atomistic DNA models

A number of intriguing aspects in dynamics of double-helical DNA is related to the coupling between its macroscopic and microscopic states. A link between the elastic properties of long DNA chains and their atom-level dynamics can be established by comparing the worm-like chain model of polymer DNA with the conformational ensembles produced by molecular dynamics simulations. This problem is complicated by the complexity of the DNA structure, the small size of DNA fragments, and relatively short trajectory durations accessible in computer simulations of microscopic DNA dynamics. A careful study of all these aspects has been performed by using longer DNA fragments and increased durations of MD trajectories as compared to earlier such investigations. Special attention is paid to the necessary conditions and criteria of time convergence, and the possibility to increase the sampling by using constrained DNA models and simplified simulation conditions. It is found that dynamics of 25-mer duplexes with regular sequences agrees well with the worm-like chain theory and that accurate evaluation of DNA elastic parameters requires at least two turns of the double helix and approximately 20-ns duration of trajectories. Bond length and bond-angle constraints affect the estimates within numerical errors. In contrast, simplified treatment of solvation can strongly change the observed elastic parameters of DNA. The elastic parameters evaluated for AT- and GC-alternating duplexes reasonably agree with experimental data and suggest that, in different basepair sequences, the torsional and stretching elasticities vary stronger than the bending stiffness.     

Forces and pressures in DNA packaging and release from viral capsids

In a previous communication (Kindt et al., 2001) we reported preliminary results of Brownian dynamics simulation and analytical theory which address the packaging and ejection forces involving DNA in bacteriophage capsids. In the present work we provide a systematic formulation of the underlying theory, featuring the energetic and structural aspects of the strongly confined DNA. The free energy of the DNA chain is expressed as a sum of contributions from its encapsidated and released portions, each expressed as a sum of bending and interstrand energies but subjected to different boundary conditions. The equilibrium structure and energy of the capsid-confined and free chain portions are determined, for each ejected length, by variational minimization of the free energy with respect to their shape profiles and interaxial spacings. Numerical results are derived for a model system mimicking the lambda-phage. We find that the fully encapsidated genome is highly compressed and strongly bent, forming a spool-like condensate, storing enormous elastic energy. The elastic stress is rapidly released during the first stage of DNA injection, indicating the large force (tens of pico Newtons) needed to complete the (inverse) loading process. The second injection stage sets in when approximately 1/3 of the genome has been released, and the interaxial distance has nearly reached its equilibrium value (corresponding to that of a relaxed torus in solution); concomitantly the encapsidated genome begins a gradual morphological transformation from a spool to a torus. We also calculate the loading force, the average pressure on the capsid's walls, and the anisotropic pressure profile within the capsid. The results are interpreted in terms of the (competing) bending and interaction components of the packing energy, and are shown to be in good agreement with available experimental data.     

Zinc affects the conformation of nucleoprotein filaments formed by replication protein A (RPA) and long natural DNA molecules

Replication protein A is the major single strand DNA binding protein of human cells, composed of three subunits with molecular weights of 70, 32, and 14 kDa. Most of the DNA binding activity of RPA has been mapped to the largest subunit that contains two OB-fold DNA binding domains and a third, OB-like structure in the carboxyterminal domain (CTD). This third domain resembles an OB-fold with a zinc binding domain inserted in the middle of the structure, and has recently been shown to carry a coordinated Zn(II) ion. The bound metal ion is essential for the tertiary structure of the RPA70-CTD, and appears to modulate its DNA binding activity when tested with synthetic oligonucleotides. We show here that zinc strongly affects the conformation of nucleoprotein filaments formed between RPA and long natural DNA molecules. In these experiments, the CTD is dispensable for DNA binding and the unwinding of long double stranded DNA molecules. However, using band shift assays and electron microscopy, we found that RPA-DNA complexes contract at zinc concentrations that do not affect the conformations of complexes formed between DNA and a RPA70 deletion construct lacking the CTD. Our data suggest that nucleoprotein complexes with RPA in its natural, zinc-bearing form may have a compact rather than an extended conformation.     

Tests of spool models for DNA packaging in phage lambda

Experiments are reported which bear on two spool models proposed for packaging the DNA of phage lambda. Both spool models fill an assumed spherical cavity with DNA wrapped in cylindrical or quasi-cylindrical layers composed of adjacent circular turns. In the curved-spool model, a single continuous segment of DNA, about 20% of the DNA length and probably located near the left end of the DNA, is in contact with the coat protein of the phage capsid. In the straight spool model, there are several DNA segments in contact with the capsid; they are concentrated in one half (probably the left half) of lambda DNA. We have identified the loci on the DNA which are in contact with the capsid by chemical crosslinking, induced by ultraviolet-irradiation of phage containing 5-bromodeoxyuridine in place of thymine.In an electron microscope experiment, phage are first lysed with EDTA, and then spread in a cytochrome c film by the formamide method. The disrupted capsid, which has the appearance of a phage ghost, serves as a marker showing where the DNA is crosslinked to the coat. The left end of the DNA is not distinguished from the right end, and so the map of DNA-capsid contacts is folded over on itself. Contacts are found nearly randomly over the entire map.In a second experiment, DNA from lysed, crosslinked phage is cut either with EcoRI or HindIII restriction endonucleases and the cut restriction fragments are labeled at their ends with ³²P. Density centrifugation in a CsCl gradient separates free DNA from restriction fragments crosslinked to protein. After digestion with proteinase k, the DNA fragments previously crosslinked to protein are identified by size after agarose gel electrophoresis. DNA fragments from all parts of the genome are found.These two experiments show that, if the DNA of each phage is packaged identically, then the curved-spool model is ruled out and the straight spool model is unlikely. Alternatively, the manner of packaging the DNA may vary from one phage to the next. These results agree with other recent experiments on λ DNA packaging by Hall & Schellman (1982a,b), and by Haas et al. (1982).A different experiment is also reported. The psoralen derivative aminomethyltrioxalen (AMT) is allowed to intercalate into λ phage and then the DNA strands are crosslinked by ultraviolet-irradiation after the rapid phase of AMT intercalation is complete. The DNA is subsequently denatured by glyoxal modification and spread for electron microscopy in a cytochrome c film by the formamide method. Sites of AMT crosslinking appear duplex; uncrosslinked regions appear as single-stranded loops. AMT is found to intercalate throughout the λ DNA. Patterns of reacted sites appear different from one DNA molecule to the next, and no consistent pattern can be found. More extensive intercalation occurs with the deletion mutant λb221 than with phage of wild-type DNA length, and free DNA shows much more reaction than the DNA inside either phage type. In order for intercalation to occur, the DNA helix must unwind and become further extended. This experiment shows that regions throughout the entire DNA molecule can unwind and be extended by intercalation, which is not confined to a single DNA segment or to segments in one half of the DNA molecule, as would be expected for the two spool models if only the DNA in contact with the capsid were accessible to the dye.     

Cloning and characterization of the beta-amylase gene from Bacillus polymyxa.

The gene for beta-amylase was isolated from Bacillus polymyxa by molecular cloning in B. subtilis. B. subtilis cells containing this gene express and secrete an amylase which resembles the B. polymyxa beta-amylase and barley beta-amylase in terms of the products it generates during carbohydrate hydrolysis. Starch hydrolysis with this beta-amylase produces maltose, not glucose, whereas maltotriose and cycloheptaose are resistant to the action of this beta-amylase. The enzyme has a molecular weight of approximately 68,000. Restriction endonuclease mapping demonstrated that the DNA inserted in pBD64 and containing the gene is approximately 3 kilobases in length.     

The stereochemistry of a four-way DNA junction: a theoretical study.

The stereochemical conformation of the four-way helical junction in DNA (the Holliday junction; the postulated central intermediate of genetic recombination) has been analysed, using molecular mechanical computer modelling. A version of the AMBER program package was employed, that had been modified to include the influence of counterions and a global optimisation procedure. Starting from an extended planar structure, the conformation was varied in order to minimise the energy, and we discuss three structures obtained by this procedure. One structure is closely related to a square-planar cross, in which there is no stacking interaction between the four double helical stems. This structure is probably closely similar to that observed experimentally in the absence of cations. The remaining two structures are based on related, yet distinct, conformations, in which there is pairwise coaxial stacking of neighbouring stems. In these structures, the four DNA stems adopt the form of two quasi-continuous helices, in which base stacking is very similar to that found in standard B-DNA geometry. The two stacked helices so formed are not aligned parallel to each other, but subtend an angle of approximately 60 degrees. The strands that exchange between one stacked helix and the other are disposed about the smaller angle of the cross (i.e. 60 degrees rather than 120 degrees), generating an approximately antiparallel alignment of DNA sequences. This structure is precisely the stacked X-structure proposed on the basis of experimental data. The calculations indicate distortions from standard B-DNA conformation that are required to adopt the stacked X-structure; a widening of the minor groove at the junction, and reorientation of the central phosphate groups of the exchanging strands. An important feature of the stacked X-structure is that it presents two structurally distinct sides. These may be recognised differently by enzymes, providing a rationalisation for the points of cleavage by Holliday resolvases.     

Langevin dynamics simulations of genome packing in bacteriophage

We use Langevin dynamics simulations to study the process by which a coarse-grained DNA chain is packaged within an icosahedral container. We focus our inquiry on three areas of interest in viral packing: the evolving structure of the packaged DNA condensate; the packing velocity; and the internal buildup of energy and resultant forces. Each of these areas has been studied experimentally, and we find that we can qualitatively reproduce experimental results. However, our findings also suggest that the phage genome packing process is fundamentally different than that suggested by the inverse spool model. We suggest that packing in general does not proceed in the deterministic fashion of the inverse-spool model, but rather is stochastic in character. As the chain configuration becomes compressed within the capsid, the structure, energy, and packing velocity all become dependent upon polymer dynamics. That many observed features of the packing process are rooted in condensed-phase polymer dynamics suggests that statistical mechanics, rather than mechanics, should serve as the proper theoretical basis for genome packing. Finally we suggest that, as a result of an internal protein unique to bacteriophage T7, the T7 genome may be significantly more ordered than is true for bacteriophage in general.     

An effective mesoscopic model of double-stranded DNA

Watson and Crick’s epochal presentation of the double helix structure in 1953 has paved the way to intense exploration of DNA’s vital functions in cells. Also, recent advances of single molecule techniques have made it possible to probe structures and mechanics of constrained DNA at length scales ranging from nanometers to microns. There have been a number of atomistic scale quantum chemical calculations or molecular level simulations, but they are too computationally demanding or analytically unfeasible to describe the DNA conformation and mechanics at mesoscopic levels. At micron scales, on the other hand, the wormlike chain model has been very instrumental in describing analytically the DNA mechanics but lacks certain molecular details that are essential in describing the hybridization, nano-scale confinement, and local denaturation. To fill this fundamental gap, we present a workable and predictive mesoscopic model of double-stranded DNA where the nucleotides beads constitute the basic degrees of freedom. With the inter-strand stacking given by an interaction between diagonally opposed monomers, the model explains with analytical simplicity the helix formation and produces a generalized wormlike chain model with the concomitant large bending modulus given in terms of the helical structure and stiffness. It also explains how the helical conformation undergoes overstretch transition to the ladder-like conformation at a force plateau, in agreement with the experiment.     

Unpaired bases in the B--Z junction of the supercoiled plasmid

The restriction analysis has been used to establish that O-beta-diethylaminoethylhydroxylamine (OHA) produces modification of unpaired cytidines in the polylinker region adjacent to the Z-insert (dG-dC)10. (dG-dC)10 in the negatively supercoiled plasmid pGC20. The length of the transition region between B- and Z-portions of DNA is not less than 36 bps. The reaction of OHA with the unpaired cytidines in the B-Z junction is a fixing one and produces no secondary despiralling of the neighboring regions. The reaction with DNA proceeds much slower than the one with monomers and single-strand polynucleotides. The structural nonuniformity has been observed, which is manifested in the alternating B and "non-B" form DNA in the B-Z junction. It is suggested that these junctions may contain nucleotide sequences which are stable to violation of the B structure during the change in superhelical density of DNA.     

Mode of DNA packing within bacteriophage heads

Electron micrographs of five different DNA bacteriophages, as prepared by drying in thin films of negative stain, frequently show their heads to be disrupted and flattened. In such cases DNA strands, no larger than 2.5 nm in diameter, become visible, either contained within partially ruptured capsids or completely ejected from severely ruptured ones. Seen in either aspect, the strands appear with circular outline; in some cases a set of concentric circles (or a tightly wound spiral) is evident.Two alternative models of DNA packing within phage heads are proposed. Both are consistent with the electron microscopic observations and, as applied specifically to T4 phage heads, they are also consistent with available data from birefringence studies. One model proposes that the DNA, in simple double-helix form, is wound into a ball. The other suggests that the DNA is wound like a spool, with a greater number of turns in the central region than at the two ends and with the spool axis perpendicular to the axis of the phage particle. The available evidence does not permit a choice to be made between the two models.     

Cleavage of a four-way DNA junction by a restriction enzyme spanning the point of strand exchange.

The four-way DNA junction is believed to fold in the presence of metal ions into an X-shaped structure, in which there is pairwise coaxial stacking of helical arms. A restriction enzyme MboII has been used to probe this structure. A junction was constructed containing a recognition site for MboII in one helical arm, positioned such that stacking of arms would result in cleavage in a neighbouring arm. Strong cleavage was observed, at the sites expected on the basis of coaxial stacking. An additional cleavage was seen corresponding to the formation of an alternative stacking isomer, suggesting that the two isomeric forms are in dynamic equilibrium in solution.     

Cloning and expression of the gene for thermostable pullulanase from Clostridium thermohydrosulfuricum in Escherichia coli

Using a pUC19-based genomic library of the anaerobic thermophilic bacterium C. thermohydrosulfuricum a DNA fragment that confers pullulanase activity to E. coli cells has been identified. Subcloning and restriction mapping procedures was carried out and the primary structure of the 5'-region of the pullulanase gene (pul) was determined. The pul enzyme was shown to be a protein with molecular weight of approximately 60,000. It was found that both pullulanase and glucoamylase activities resides in pullulanase. The intracellular distribution of pullulanase was studied. An E. coli strain that produces large amounts of thermostable pullulanase has been constructed.     

Structure prediction of a complex between the chromosomal protein HMG-D and DNA

Non-histone chromosomal proteins are an important part of nuclear structure and function due to their ability to interact with DNA to form and modulate chromatin structure and regulate gene expression. However, the understanding of the function of chromosomal proteins at the molecular level has been hampered by the lack of structures of chromosomal protein–DNA complexes. We have carried out a molecular dynamics modeling study to provide insight into the mode of DNA binding to the chromosomal HMG-domain protein, HMG-D. Three models of a complex of HMG-D bound to DNA were derived through docking the protein to two different DNA fragments of known structure. Molecular dynamics simulations of the complexes provided data indicating the most favorable model. This model was further refined by molecular dynamics simulation and extensively analyzed. The structure of the corresponding HMG-D-DNA complex exhibits many features seen in the NMR structures of the sequence-specific HMG-domain-DNA complexes, lymphoid enhancer factor 1 (LEF-1) and testis determining factor (SRY). The model reveals differences from these known structures that suggest how chromosomal proteins bind to many different DNA sequences with comparable affinity. Proteins 30:113–135, 1998. © 1998 Wiley-Liss, Inc.     

The nucleotide sequence of a HMW glutenin subunit gene located on chromosome 1A of wheat (Triticum aestivum L.).

A cloned 8.2 kb EcoRI fragment has been isolated from a genomic library of DNA derived from Triticum aestivum L. cv. Cheyenne. This fragment contains sequences related to the high molecular weight (HMW) subunits of glutenin, proteins considered to be important in determining the elastic properties of gluten. The cloned HMW subunit gene appears to be derived from chromosome 1A. The nucleotide sequence of this gene has provided new information on the structure and evolution of the HMW subunits. However, hybrid-selection translation experiments suggest that this gene is silent.     

Structure of a B-DNA decamer with a central T-A step: C-G-A-T-T-A-A-T-C-G.

The X-ray crystal structure analysis of the decamer C-G-A-T-T-A-A-T-C-G has been carried out to a resolution of 1.5 A. The crystals are space group P2(1)2(1)2(1), cell dimensions a = 38.60 A, b = 39.10 A, c = 33.07 A. The structure was solved by molecular replacement and refined with X-PLOR and NUCLSQ. The final R factor for a model with 404 DNA atoms, 108 water molecules and one magnesium hexahydrate cation is 15.7%. The double helix is essentially isostructural with C-G-A-T-C-G-A-T-C-G, with closely similar local helix parameters. The structure of the T-T-A-A center differs from that found in C-G-C-G-T-T-A-A-C-G-C-G in that the minor groove in our decamer is wide at the central T-A step rather than narrow, and the twist angle of the T-A step is small (31.1 degrees) rather than large. Whereas the tetrad model provides a convenient framework for discussing local DNA helix structure, it cannot be the entire story. The articulated helix model of DNA structure proposes that certain sequence regions of DNA show preferential twisting or bending properties, whereas other regions are less capable of deformation, in a manner that may be useful in sequence recognition by drugs and protein. Further crystal structure analyses should help to delineate the precise nature of sequence-dependent articulation in the DNA double helix.     

The complex between a four-way DNA junction and T7 endonuclease I

The junction-resolving enzyme endonuclease I is selective for the structure of the DNA four-way (Holliday) junction. The enzyme binds to a four-way junction in two possible orientations, with a 4:1 ratio, opening the DNA structure at the centre and changing the global structure into a 90 degrees cross of approximately coaxial helices. The nuclease cleaves the continuous strands of the junction in each orientation. Binding leads to pronounced regions of protection of the DNA against hydroxyl radical attack. Using all this information together with the known structure of the enzyme and the structure of the BglI-DNA complex, we have constructed a model of the complex of endonuclease I and a DNA junction. This shows how the enzyme is selective for the structure of a four-way junction, such that both continuous strands can be accommodated into the two active sites so that a productive resolution event is possible.