Modeling and simulation of the mechanical response from nanoindentation test of DNA-filled viral capsids

Viruses can be described as biological objects composed mainly of two parts: a stiff protein shell called a capsid, and a core inside the capsid containing the nucleic acid and liquid. In many double-stranded DNA bacterial viruses (aka phage), the volume ratio between the liquid and the encapsidated DNA is approximately 1:1. Due to the dominant DNA hydration force, water strongly mediates the interaction between the packaged DNA strands. Therefore, water that hydrates the DNA plays an important role in nanoindentation experiments of DNA-filled viral capsids. Nanoindentation measurements allow us to gain further insight into the nature of the hydration and electrostatic interactions between the DNA strands. With this motivation, a continuum-based numerical model for simulating the nanoindentation response of DNA-filled viral capsids is proposed here. The viral capsid is modeled as large- strain isotropic hyper-elastic material, whereas porous elasticity is adopted to capture the mechanical response of the filled viral capsid. The voids inside the viral capsid are assumed to be filled with liquid, which is modeled as a homogenous incompressible fluid. The motion of a fluid flowing through the porous medium upon capsid indentation is modeled using Darcy’s law, describing the flow of fluid through a porous medium. The nanoindentation response is simulated using three-dimensional finite element analysis and the simulations are performed using the finite element code Abaqus. Force-indentation curves for empty, partially and completely DNA-filled capsids are directly compared to the experimental data for bacteriophage λ. Material parameters such as Young’s modulus, shear modulus, and bulk modulus are determined by comparing computed force-indentation curves to the data from the atomic force microscopy (AFM) experiments. Predictions are made for pressure distribution inside the capsid, as well as the fluid volume ratio variation during the indentation test.     

Interactions of DNA with divalent metal ions. II. Proton relaxation enhancement studies

The extent and modes of binding of the divalent metal ions Mn²⁺ and Co²⁺ to DNA and the effects of salt on the binding have been studied by measurements of the effects of these paramagnetic metal ions on the longitudinal and transverse relaxation rates of the protons of the solvent water molecules, a technique that is sensitive to overall binding. The number of water molecules coordinated to the DNA–bound Mn²⁺ and Co²⁺ is found to be between five and six, and the electron spin relaxation times and the electron-nuclear hyperfine constants associated with Mn²⁺ and Co²⁺ are little or not affected by the binding. These observations indicate little disturbance of the hydration sphere of Mn²⁺ and Co²⁺ upon binding to DNA. An average 2–3-fold reduction in the exchange rate of the water of hydration of the bound metal ions and an order-of-magnitude increase in their rotational correlation time are attributed to hydrogen-bond formation with the DNA. The binding constants of Mn²⁺ to DNA, at metal concentrations approaching zero, are found to be inversely proportional to the second power of the salt concentration, in agreement with the predictions of Manning's polyelectrolyte theory. A remarkable quantitative agreement with the polyelectrolyte theory is also obtained for the anticooperativity in the binding of Mn²⁺ to DNA, although the experimental results can be well accounted for by another simple electrostatic model. The various modes of binding of divalent metal ions to DNA are discussed.     

Cations Form Sequence Selective Motifs within DNA Grooves via a Combination of Cation-Pi and Ion-Dipole/Hydrogen Bond Interactions

The fine conformational subtleties of DNA structure modulate many fundamental cellular processes including gene activation/repression, cellular division, and DNA repair. Most of these cellular processes rely on the conformational heterogeneity of specific DNA sequences. Factors including those structural characteristics inherent in the particular base sequence as well as those induced through interaction with solvent components combine to produce fine DNA structural variation including helical flexibility and conformation. Cation-pi interactions between solvent cations or their first hydration shell waters and the faces of DNA bases form sequence selectively and contribute to DNA structural heterogeneity. In this paper, we detect and characterize the binding patterns found in cation-pi interactions between solvent cations and DNA bases in a set of high resolution x-ray crystal structures. Specifically, we found that monovalent cations (Tl⁺) and the polarized first hydration shell waters of divalent cations (Mg²⁺, Ca²⁺) form cation-pi interactions with DNA bases stabilizing unstacked conformations. When these cation-pi interactions are combined with electrostatic interactions a pattern of specific binding motifs is formed within the grooves.     

Hydration of counterions interacting with DNA double helix: a molecular dynamics study

In the present work, molecular dynamics simulations have been carried out to study the dependence of counterion distribution around the DNA double helix on the character of ion hydration. The simulated systems consisted of DNA fragment d(CGCGAATTCGCG) in water solution with the counterions Na⁺, K⁺, Cs⁺ or Mg²⁺. The characteristic binding sites of the counterions with DNA and the changes in their hydration shell have been determined. The results show that due to the interaction with DNA at least two hydration shells of the counterions undergo changes. The first hydration shell of Na⁺, K⁺, Cs⁺, and Mg²⁺ counterions in the bulk consists of six, seven, ten, and six water molecules, respectively, while the second one has several times higher values. The Mg²⁺ and Na⁺ counterions, constraining water molecules of the first hydration shell, mostly form with DNA water-mediated contacts. In this case the coordination numbers of the first hydration shell do not change, while the coordination numbers of the second one decrease about twofold. The Cs⁺ and K⁺ counterions that do not constrain surrounding water molecules may be easily dehydrated, and when interacting with DNA their first hydration shell may be decreased by three and five water molecules, respectively. Due to the dehydration effect, these counterions can squeeze through the hydration shell of DNA to the bottom of the double helix grooves. The character of ion hydration establishes the correlation between the coordination numbers of the first and the second hydration shells.

Open image in new window

Graphical Abstract Hydration of counterions interacting with DNA double helix

     

Structure and Energetics of Encapsidated DNA in Bacteriophage HK97 Studied by Scanning Calorimetry and Cryo-electron Microscopy

Encapsidation of duplex DNA by bacteriophages represents an extreme case of genome condensation, reaching near-crystalline concentrations of DNA. The HK97 system is well suited to study this phenomenon in view of the detailed knowledge of its capsid structure. To characterize the interactions involved, we combined calorimetry with cryo-electron microscopy and native gel electrophoresis. We found that, as in other phages, HK97 DNA is organized in coaxially wound nested shells. When DNA-filled capsids (heads) are scanned in buffer containing 1 mM Mg²⁺, DNA melting and capsid denaturation both contribute to the complex thermal profile between 82 °C and 96 °C. In other conditions (absence of Mg²⁺ and lower ionic strength), DNA melting shifts to lower temperatures and the two events are resolved. Heads release their DNA at temperatures well below the onset of DNA melting or capsid denaturation. We suggest that, on heating, the internal pressure increases, causing the DNA to exit—probably via the portal vertex-while the capsid, although largely intact, sustains local damage that leads to an earlier onset of thermal denaturation. Heads differ structurally from empty capsids in the curvature of their protein shell, a change attributable to outwards pressure exerted by the DNA. We propose that this transition is sensed by the portal that is embedded in the capsid wall, whereupon the structure of the portal and its interactions with terminase, the packaging enzyme, are altered, thus signaling that packaging is at or approaching completion.     

Swelling and Softening of the Cowpea Chlorotic Mottle Virus in Response to pH Shifts

Cowpea chlorotic mottle virus (CCMV) forms highly elastic icosahedral protein capsids that undergo a characteristic swelling transition when the pH is raised from 5 to 7. Here, we performed nano-indentation experiments using an atomic force microscope to track capsid swelling and measure the shells’ Young’s modulus at the same time. When we chelated Ca²⁺ ions and raised the pH, we observed a gradual swelling of the RNA-filled capsids accompanied by a softening of the shell. Control experiments with empty wild-type virus and a salt-stable mutant revealed that the softening was not strictly coupled to the swelling of the protein shells. Our data suggest that a pH increase and Ca²⁺ chelation lead primarily to a loosening of contacts within the protein shell, resulting in a softening of the capsid. This appears to render the shell metastable and make swelling possible when repulsive forces among the capsid proteins become large enough, which is known to be followed by capsid disassembly at even higher pH. Thus, softening and swelling are likely to play a role during inoculation.     

Effects of salt concentrations and bending energy on the extent of ejection of phage genomes

Recent work has shown that pressures inside dsDNA phage capsids can be as high as many tens of atmospheres; it is this pressure that is responsible for initiation of the delivery of phage genomes to host cells. The forces driving ejection of the genome have been shown to decrease monotonically as ejection proceeds, and hence to be strongly dependent on the genome length. Here we investigate the effects of ambient salts on the pressures inside phage-λ, for the cases of mono-, di-, and tetravalent cations, and measure how the extent of ejection against a fixed osmotic pressure (mimicking the bacterial cytoplasm) varies with cation concentration. We find, for example, that the ejection fraction is halved in 30 mM Mg²⁺ and is decreased by a factor of 10 upon addition of 1 mM spermine. These effects are calculated from a simple model of genome packaging, using DNA-DNA repulsion energies as determined independently from x-ray diffraction measurements on bulk DNA solutions. By comparing the measured ejection fractions with values implied from the bulk DNA solution data, we predict that the bending energy makes the d-spacings inside the capsid larger than those for bulk DNA at the same osmotic pressure.     

Buoyant density sedimentation of macromolecules in sodium iothalamate density gradients.

Nucleic acids, bacteriophages, phage capsids, and a DNA-capsid complex have been centrifuged to an equilibrium buoyant density in sodium iothalamate density gradients. Nucleic acids have comparatively high hydrations and are less dense than proteins in these gradients. Sodium iothalamate gradients can be used to separate DNA from RNA, single-chain DNA from double-chain DNA and to separate bacteriophage T7 and λ deletion mutants from the respective wild-type phage.The DNA packaged in bacteriophage T7 appears to be less hydrated than free DNA in sodium iothalamate gradients. There is evidence that the hydration of DNA packaged in phage T7 is restricted by the volume of the phage head. The total volume of phage T7 was estimated to be 1.32 × 10⁻¹⁶ ml. The volume available to package phage T7 DNA was estimated to be 2.2 times the volume of the B form of T7 DNA.     

Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

Calcium ions (Ca²⁺) play key roles in various fundamental biological processes such as cell signaling and brain function. Molecular dynamics (MD) simulations have been used to study such interactions, however, the accuracy of the Ca²⁺ models provided by the standard MD force fields has not been rigorously tested. Here, we assess the performance of the Ca²⁺ models from the most popular classical force fields AMBER and CHARMM by computing the osmotic pressure of model compounds and the free energy of DNA–DNA interactions. In the simulations performed using the two standard models, Ca²⁺ ions are seen to form artificial clusters with chloride, acetate, and phosphate species; the osmotic pressure of CaAc₂ and CaCl₂ solutions is a small fraction of the experimental values for both force fields. Using the standard parameterization of Ca²⁺ ions in the simulations of Ca²⁺‐mediated DNA–DNA interactions leads to qualitatively wrong outcomes: both AMBER and CHARMM simulations suggest strong inter‐DNA attraction whereas, in experiment, DNA molecules repel one another. The artificial attraction of Ca²⁺ to DNA phosphate is strong enough to affect the direction of the electric field‐driven translocation of DNA through a solid‐state nanopore. To address these shortcomings of the standard Ca²⁺ model, we introduce a custom model of a hydrated Ca²⁺ ion and show that using our model brings the results of the above MD simulations in quantitative agreement with experiment. Our improved model of Ca²⁺ can be readily applied to MD simulations of various biomolecular systems, including nucleic acids, proteins and lipid bilayer membranes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 752–763, 2016.     

Cation Selectivity in Biological Cation Channels Using Experimental Structural Information and Statistical Mechanical Simulation

Cation selective channels constitute the gate for ion currents through the cell membrane. Here we present an improved statistical mechanical model based on atomistic structural information, cation hydration state and without tuned parameters that reproduces the selectivity of biological Na⁺ and Ca²⁺ ion channels. The importance of the inclusion of step-wise cation hydration in these results confirms the essential role partial dehydration plays in the bacterial Na⁺ channels. The model, proven reliable against experimental data, could be straightforwardly used for designing Na⁺ and Ca²⁺ selective nanopores.     

Replication of herpesvirus DNA. V. Maturation of concatemeric DNA of pseudorabies virus to genome length is related to capsid formation.

The maturation of pseudorabies virus DNA from the replicative concatemeric form to molecules of genome length was examined using nine DNA+ temperature-sensitive mutants of pseudorabies virus, each belonging to a different complementation group. At the nonpermissive temperature, cells infected with each of the mutants synthesized concatemeric DNA. Cleavage of the concatemeric DNA to genome-length viral DNA was defective in all the DNA+ ts mutants tested, indicating that several viral gene products are involved in the DNA maturation process. In none of the ts mutant-infected cells were capsids with electron-dense cores (containing DNA) formed. Empty capsids with electron-translucent cores were, however, formed in cells infected with six of the nine temperature-sensitive mutants; in cells infected with three of the mutants, no capsid assembly occurred. Because these three mutants are deficient both in maturation of DNA and in the assembly of viral capsids, we conclude that maturation of viral DNA is dependent upon the assembly of capsids. In cells infected with two of the mutants (tsN and tsIE13), normal maturation of viral DNA occurred after shiftdown of the cells to the permissive temperature in the presence of cycloheximide, indicating that the temperature-sensitive proteins involved in DNA maturation became functional after shiftdown. Furthermore, because cycloheximide reduces maturation of DNA in wild-type-infected cells but not in cells infected with these two mutants, we conclude that a protein(s) necessary for the maturation of concatemeric DNA, which is present in limiting amounts during the normal course of infection, accumulated in the mutant-infected cells at the nonpermissive temperature. Concomitant with cleavage of concatemeric DNA, full capsids with electron-dense cores appeared after shiftdown of tsN-infected cells to the permissive temperature, indicating that there may be a correlation between maturation of DNA and formation of full capsids. The number of empty and full capsids (containing electron-dense cores) present in tsN-infected cells incubated at the nonpermissive temperature, as well as after shiftdown to the permissive temperature in the presence of cycloheximide, was determined by electron microscopy and by sedimentation analysis in sucrose gradients. After shiftdown to the permissive temperature in the presence of cycloheximide, the number of empty capsids present in tsN-infected cells decreased with a concomitant accumulation of full capsids, indicating that empty capsids are precursors to full capsids.     

Molecular Dynamics Simulation Studies of the Limiting Conductances of MgCl2 and CaCl2 in Supercritical Water Using SPC/E Model for Water

We report results of molecular dynamics (MD) simulations of the limiting conductances of MgCl² and CaCl² in supercritical water as a function of water density using the SPC/E model for water. The limiting conductances of Mg²⁺, Ca²⁺ and Cl^- over the whole range of water density considered exhibits a linear dependence of the limiting conductance on the water density. In the cases of Mg²⁺ and Ca²⁺, a solventberg picture for the behavior of small divalent cation emerges from our studies. From the view of the solventberg picture, the ion and its shell moving together as an entity interacts with the second hydration shell water molecules, and its mobility is restricted mostly by the number of the second hydration shell water which is proportional to the water density of the whole system. In the case of Cl^-, the range of water density considered in this study belongs to the higher-density region (above 0.45?g/cm³) in which the effect of the number of hydration water molecules around ions dominated. As the water density increases, the water molecules of the first hydration shell restrict the mobility of Cl^- and the limiting conductance of Cl^- decreases nearly linearly. Significant different dependence on the water density is observed between the calculated limiting conductances of MgCl² and CaCl² at 673?K and the experimental results over the water density of 0.60–0.90?g/cm³. Possible limitation of the extended simple point charge (SPC/E) model with regard to this difference should be pointed out and the use of a more precise model like the revised polarizable (RPOL) model is indispensable for a further MD study to gain a complete picture of the chemical circumstance around the ions.     

Cyclic AMP-dependent signaling system is a primary metabolic target for non-thermal effect of microwaves on heart muscle hydration

Previously, we have suggested that cell hydration is a universal and extra-sensitive sensor for the structural changes of cell aqua medium caused by the impact of weak chemical and physical factors. The aim of present work is to elucidate the nature of the metabolic messenger through which physiological solution (PS) treated by non-thermal (NT) microwaves (MW) could modulate heart muscle hydration of rats. For this purpose, the effects of NT MW–treated PS on heart muscle hydration, [³H]-ouabain binding with cell membrane, ⁴⁵Ca²⁺ uptake and intracellular cyclic nucleotides contents in vivo and in vitro experiments were studied. It is shown that intraperitoneal injections of both Sham-treated PS and NT MW–treated PS elevate heart muscle hydration. However, the effect of NT MW–treated PS on muscle hydration is more pronounced than the effect of Sham-treated PS. In vitro experiments NT MW–treated PS has dehydration effect on muscle, which is not changed by decreasing Na⁺ gradients on membrane. Intraperitoneal injection of Sham- and NT MW–treated PS containing ⁴⁵Ca²⁺ have similar dehydration effect on muscle, while NT MW–treated PS has activation effect on Na⁺/Ca²⁺ exchange in reverse mode. The intraperitoneal injection of NT MW–treated PS depresses [³H]-ouabain binding with its high-affinity membrane receptors, elevates intracellular cAMP and decreases cGMP contents. Based on the obtained data, it is suggested that cAMP-dependent signaling system serves as a primary metabolic target for NT MW effect on heart muscle hydration.     

Simulate the diffusion of hydrated ions by nanofiltration membrane process with random walk

We propose a new diffusion model describing the diffusion behaviours of hydrated ions in the process of nanofiltration (NF) based on the random walk (RW) theory when the NF membrane is uncharged or low charged. In this model, the hydration of ions and their deformation capacity are considered. The structure of the membrane is idealised into a lozenge shape and the diameter of membrane pore is defined as gapsize. A computer program named RW system in chemistry is developed to simulate based on this model. Six familiar ions Li⁺, Na⁺, Mg²⁺, Al³⁺, K⁺ and Ca²⁺ are investigated. Their characteristics are calculated by Gaussian 03, Pople, Inc., Wallingford, CT. The diffusivities of hydrated ions are calculated and discussed. The results show that the hydration of ions cannot be ignored in NF process when the membrane pore size is near the dimensions of the hydrated ions.     

Structure of T4 polyheads: II. A pathway of polyhead transformations as a model for T4 capsid maturation

We have studied the aberrant tubular polyheads of bacteriophages T4D and T2L as a model system for capsid maturation. Six different types of polyhead surface lattice morphology, and the corresponding protein compositions are reported and discussed. Using in vitro systems to induce transformations between particular polyhead types, we have deduced that the structural classes represent successive points in a transitional pathway. In the first step, coarse polyheads (analogous to the prohead τ-particle) are proteolytically cleaved by a phagecoded protease, a fragment of the gene 21 product. This cleavage of P23 to P23¹ induces a co-operative lattice transformation in the protein of the surface shell, to a conformation equivalent to that of T2L giant phage capsids. These polyheads (derived either from T4 or T2L lysates) can accept further T4-coded proteins. In doing so, they pass through intermediate structural states, eventually reaching an end point whose unit cell morphology is indistinguishable from that of the giant T4 capsids. At least one protein (called soc (Ishii & Yanagida, 1975)) is bound stoichiometrically to P23¹ in the end-state conformation. The simulation of several aspects of capsid maturation (cleavage of P23 to P23¹, stabilization, and lattice expansion) in the polyhead pathway suggest that it parallels the major events of phage T-even capsid maturation, decoupled from any involvement of DNA packaging.     

Ion Competition in Condensed DNA Arrays in the Attractive Regime

Physical origin of DNA condensation by multivalent cations remains unsettled. Here, we report quantitative studies of how one DNA-condensing ion (Cobalt³⁺ Hexammine, or Co³⁺Hex) and one nonDNA-condensing ion (Mg²⁺) compete within the interstitial space in spontaneously condensed DNA arrays. As the ion concentrations in the bath solution are systematically varied, the ion contents and DNA-DNA spacings of the DNA arrays are determined by atomic emission spectroscopy and x-ray diffraction, respectively. To gain quantitative insights, we first compare the experimentally determined ion contents with predictions from exact numerical calculations based on nonlinear Poisson-Boltzmann equations. Such calculations are shown to significantly underestimate the number of Co³⁺Hex ions, consistent with the deficiencies of nonlinear Poisson-Boltzmann approaches in describing multivalent cations. Upon increasing the concentration of Mg²⁺, the Co³⁺Hex-condensed DNA array expands and eventually redissolves as a result of ion competition weakening DNA-DNA attraction. Although the DNA-DNA spacing depends on both Mg²⁺ and Co³⁺Hex concentrations in the bath solution, it is observed that the spacing is largely determined by a single parameter of the DNA array, the fraction of DNA charges neutralized by Co³⁺Hex. It is also observed that only ∼20% DNA charge neutralization by Co³⁺Hex is necessary for spontaneous DNA condensation. We then show that the bath ion conditions can be reduced to one variable with a simplistic ion binding model, which is able to describe the variations of both ion contents and DNA-DNA spacings reasonably well. Finally, we discuss the implications on the nature of interstitial ions and cation-mediated DNA-DNA interactions.     

Frequency of mispackaging of Prochlorococcus DNA by cyanophage

Prochlorococcus cells are the numerically dominant phototrophs in the open ocean. Cyanophages that infect them are a notable fraction of the total viral population in the euphotic zone, and, as vehicles of horizontal gene transfer, appear to drive their evolution. Here we examine the propensity of three cyanophages—a podovirus, a siphovirus, and a myovirus—to mispackage host DNA in their capsids while infecting Prochlorococcus, the first step in phage-mediated horizontal gene transfer. We find the mispackaging frequencies are distinctly different among the three phages. Myoviruses mispackage host DNA at low and seemingly fixed frequencies, while podo- and siphoviruses vary in their mispackaging frequencies by orders of magnitude depending on growth light intensity. We link this difference to the concentration of intracellular reactive oxygen species and protein synthesis rates, both parameters increasing in response to higher light intensity. Based on our findings, we propose a model of mispackaging frequency determined by the imbalance between the production of capsids and the number of phage genome copies during infection: when protein synthesis rate increase to levels that the phage cannot regulate, they lead to an accumulation of empty capsids, in turn triggering more frequent host DNA mispackaging errors.Subject terms: Environmental microbiology, Ecology     

Gating-related Structural Dynamics of the MgtE Magnesium Channel in Membrane-Mimetics Utilizing Site-Directed Tryptophan Fluorescence

Magnesium is the most abundant divalent cation present in the cell, and an abnormal Mg²⁺ homeostasis is associated with several diseases in humans. However, among ion channels, the mechanisms of intracellular regulation and transport of Mg²⁺ are poorly understood. MgtE is a homodimeric Mg²⁺-selective channel and is negatively regulated by high intracellular Mg²⁺ concentration where the cytoplasmic domain of MgtE acts as a Mg²⁺ sensor. Most of the previous biophysical studies on MgtE have been carried out in detergent micelles and the information regarding gating-related structural dynamics of MgtE in physiologically-relevant membrane environment is scarce. In this work, we monitored the changes in gating-related structural dynamics, hydration dynamics and conformational heterogeneity of MgtE in micelles and membranes using the intrinsic site-directed Trp fluorescence. For this purpose, we have engineered six single-Trp mutants in the functional Trp-less background of MgtE to obtain site-specific information on the gating-related structural dynamics of MgtE in membrane-mimetic systems. Our results indicate that Mg²⁺-induced gating might involve the possibility of a ‘conformational wave’ from the cytosolic N-domain to transmembrane domain of MgtE. Although MgtE is responsive to Mg²⁺-induced gating in both micelles and membranes, the organization and dynamics of MgtE is substantially altered in physiologically important phospholipid membranes compared to micelles. This is accompanied by significant changes in hydration dynamics and conformational heterogeneity. Overall, our results highlight the importance of lipid-protein interactions and are relevant for understanding gating mechanism of magnesium channels in general, and MgtE in particular.     

Nuclear egress and envelopment of herpes simplex virus capsids analyzed with dual-color fluorescence HSV1(17+)

To analyze the assembly of herpes simplex virus type 1 (HSV1) by triple-label fluorescence microscopy, we generated a bacterial artificial chromosome (BAC) and inserted eukaryotic Cre recombinase, as well as β-galactosidase expression cassettes. When the BAC pHSV1(17⁺)blueLox was transfected back into eukaryotic cells, the Cre recombinase excised the BAC sequences, which had been flanked with loxP sites, from the viral genome, leading to HSV1(17⁺)blueLox. We then tagged the capsid protein VP26 and the envelope protein glycoprotein D (gD) with fluorescent protein domains to obtain HSV1(17⁺)blueLox-GFPVP26-gDRFP and -RFPVP26-gDGFP. All HSV1 BACs had variations in the a-sequences and lost the oriL but were fully infectious. The tagged proteins behaved as their corresponding wild type, and were incorporated into virions. Fluorescent gD first accumulated in cytoplasmic membranes but was later also detected in the endoplasmic reticulum and the plasma membrane. Initially, cytoplasmic capsids did not colocalize with viral glycoproteins, indicating that they were naked, cytosolic capsids. As the infection progressed, they were enveloped and colocalized with the viral membrane proteins. We then analyzed the subcellular distribution of capsids, envelope proteins, and nuclear pores during a synchronous infection. Although the nuclear pore network had changed in ca. 20% of the cells, an HSV1-induced reorganization of the nuclear pore architecture was not required for efficient nuclear egress of capsids. Our data are consistent with an HSV1 assembly model involving primary envelopment of nuclear capsids at the inner nuclear membrane and primary fusion to transfer capsids into the cytosol, followed by their secondary envelopment on cytoplasmic membranes.     

Ion Competition in Condensed DNA Arrays in the Attractive Regime

Physical origin of DNA condensation by multivalent cations remains unsettled. Here, we report quantitative studies of how one DNA-condensing ion (Cobalt³⁺ Hexammine, or Co³⁺Hex) and one nonDNA-condensing ion (Mg²⁺) compete within the interstitial space in spontaneously condensed DNA arrays. As the ion concentrations in the bath solution are systematically varied, the ion contents and DNA-DNA spacings of the DNA arrays are determined by atomic emission spectroscopy and x-ray diffraction, respectively. To gain quantitative insights, we first compare the experimentally determined ion contents with predictions from exact numerical calculations based on nonlinear Poisson-Boltzmann equations. Such calculations are shown to significantly underestimate the number of Co³⁺Hex ions, consistent with the deficiencies of nonlinear Poisson-Boltzmann approaches in describing multivalent cations. Upon increasing the concentration of Mg²⁺, the Co³⁺Hex-condensed DNA array expands and eventually redissolves as a result of ion competition weakening DNA-DNA attraction. Although the DNA-DNA spacing depends on both Mg²⁺ and Co³⁺Hex concentrations in the bath solution, it is observed that the spacing is largely determined by a single parameter of the DNA array, the fraction of DNA charges neutralized by Co³⁺Hex. It is also observed that only ∼20% DNA charge neutralization by Co³⁺Hex is necessary for spontaneous DNA condensation. We then show that the bath ion conditions can be reduced to one variable with a simplistic ion binding model, which is able to describe the variations of both ion contents and DNA-DNA spacings reasonably well. Finally, we discuss the implications on the nature of interstitial ions and cation-mediated DNA-DNA interactions.