DNA Scrunching in the Packaging of Viral Genomes

The motors that drive double-stranded DNA (dsDNA) genomes into viral capsids are among the strongest of all biological motors for which forces have been measured, but it is not known how they generate force. We previously proposed that the DNA is not a passive substrate but that it plays an active role in force generation. This "scrunchworm hypothesis" holds that the motor proteins repeatedly dehydrate and rehydrate the DNA, which then undergoes cyclic shortening and lengthening motions. These are captured by a coupled protein-DNA grip-and-release cycle to rectify the motion and translocate the DNA into the capsid. In this study, we examined the interactions of dsDNA with the dodecameric connector protein of bacteriophage ϕ29, using molecular dynamics simulations on four different DNA sequences, starting from two different conformations (A-DNA and B-DNA). In all four simulations starting with the protein equilibrated with A-DNA in the channel, we observed transitions to a common, metastable, highly scrunched conformation, designated A*. This conformation is very similar to one recently reported by Kumar and Grubmüller in much longer MD simulations on B-DNA docked into the ϕ29 connector. These results are significant for four reasons. First, the scrunched conformations occur spontaneously, without requiring lever-like protein motions often believed to be necessary for DNA translocation. Second, the transition takes place within the connector, providing the location of the putative "dehydrator". Third, the protein has more contacts with one strand of the DNA than with the other; the former was identified in single-molecule laser tweezer experiments as the "load-bearing strand". Finally, the spontaneity of the DNA-protein interaction suggests that it may play a role in the initial docking of DNA in motors like that of T4 that can load and package any sequence.

Force distribution in a semiflexible loop

Loops undergoing thermal fluctuations are prevalent in nature. Ringlike or cross-linked polymers, cyclic macromolecules, and protein-mediated DNA loops all belong to this category. Stability of these molecules are generally described in terms of free energy, an average quantity, but it may also be impacted by local fluctuating forces acting within these systems. The full distribution of these forces can thus give us insights into mechanochemistry beyond the predictive capability of thermodynamics. In this paper, we study the force exerted by an inextensible semiflexible polymer constrained in a looped state. By using a simulation method termed "phase-space sampling," we generate the equilibrium distribution of chain conformations in both position and momentum space. We compute the constraint forces between the two ends of the loop in this chain ensemble using Lagrangian mechanics, and show that the mean of these forces is equal to the thermodynamic force. By analyzing kinetic and potential contributions to the forces, we find that the mean force acts in the direction of increasing extension not because of bending stress, but in spite of it. Furthermore, we obtain a distribution of constraint forces as a function of chain length, extension, and stiffness. Notably, increasing contour length decreases the average force, but the additional freedom allows fluctuations in the constraint force to increase. The force distribution is asymmetric and falls off less sharply than a Gaussian distribution. Our work exemplifies a system where large-amplitude fluctuations occur in a way unforeseen by a purely thermodynamic framework, and offers computational tools useful for efficient, unbiased simulation of a constrained system.

Transitions of Double-Stranded DNA Between the A- and B-Forms

The structure of double-stranded DNA (dsDNA) is sensitive to solvent conditions. In solution, B-DNA is the favored conformation under physiological conditions, while A-DNA is the form found under low water activity. The A-form is induced locally in some protein-DNA complexes, and repeated transitions between the B- and A-forms have been proposed to generate the forces used to drive dsDNA into viral capsids during genome packaging. Here, we report analyses on previous molecular dynamics (MD) simulations on B-DNA, along with new MD simulations on the transition from A-DNA to B-DNA in solution. We introduce the A-B Index (ABI), a new metric along the A-B continuum, to quantify our results. When A-DNA is placed in an equilibrated solution at physiological ionic strength, there is no energy barrier to the transition to the B-form, which begins within about 1 ns. The transition is essentially complete within 5 ns, although occasionally a stretch of a few base pairs will remain A-like for up to ∼10 ns. A comparison of four sequences with a range of predicted A-phobicities shows that more A-phobic sequences make the transition more rapidly than less A-phobic sequences. Simulations on dsDNA with a region of roughly one turn locked in the A-form allow us to characterize the A/B junction, which has an average bend angle of 20-30°. Fluctuations in this angle occur with characteristic times of about 10 ns.

The Effect of Basepair Mismatch on DNA Strand Displacement

DNA strand displacement is a key reaction in DNA homologous recombination and DNA mismatch repair and is also heavily utilized in DNA-based computation and locomotion. Despite its ubiquity in science and engineering, sequence-dependent effects of displacement kinetics have not been extensively characterized. Here, we measured toehold-mediated strand displacement kinetics using single-molecule fluorescence in the presence of a single basepair mismatch. The apparent displacement rate varied significantly when the mismatch was introduced in the invading DNA strand. The rate generally decreased as the mismatch in the invader was encountered earlier in displacement. Our data indicate that a single base pair mismatch in the invader stalls branch migration and displacement occurs via direct dissociation of the destabilized incumbent strand from the substrate strand. We combined both branch migration and direct dissociation into a model, which we term the concurrent displacement model, and used the first passage time approach to quantitatively explain the salient features of the observed relationship. We also introduce the concept of splitting probabilities to justify that the concurrent model can be simplified into a three-step sequential model in the presence of an invader mismatch. We expect our model to become a powerful tool to design DNA-based reaction schemes with broad functionality.

Calculation of a fluctuating entropic force by phase space sampling

A polymer chain pinned in space exerts a fluctuating force on the pin point in thermal equilibrium. The average of such fluctuating force is well understood from statistical mechanics as an entropic force, but little is known about the underlying force distribution. Here, we introduce two phase space sampling methods that can produce the equilibrium distribution of instantaneous forces exerted by a terminally pinned polymer. In these methods, both the positions and momenta of mass points representing a freely jointed chain are perturbed in accordance with the spatial constraints and the Boltzmann distribution of total energy. The constraint force for each conformation and momentum is calculated using Lagrangian dynamics. Using terminally pinned chains in space and on a surface, we show that the force distribution is highly asymmetric with both tensile and compressive forces. Most importantly, the mean of the distribution, which is equal to the entropic force, is not the most probable force even for long chains. Our work provides insights into the mechanistic origin of entropic forces, and an efficient computational tool for unbiased sampling of the phase space of a constrained system.

Probing the elastic limit of DNA bending

Sharp bending of double-stranded DNA (dsDNA) plays an essential role in genome structure and function. However, the elastic limit of dsDNA bending remains controversial. Here, we measured the opening rates of small dsDNA loops with contour lengths ranging between 40 and 200 bp using single-molecule Fluorescence Resonance Energy Transfer. The relationship of loop lifetime to loop size revealed a critical transition in bending stress. Above the critical loop size, the loop lifetime changed with loop size in a manner consistent with elastic bending stress, but below it, became less sensitive to loop size, indicative of softened dsDNA. The critical loop size increased from ∼60 bp to ∼100 bp with the addition of 5 mM magnesium. We show that our result is in quantitative agreement with the kinkable worm-like chain model, and furthermore, can reproduce previously reported looping probabilities of dsDNA over the range between 50 and 200 bp. Our findings shed new light on the energetics of sharply bent dsDNA.

Studying DNA looping by single-molecule FRET

Bending of double-stranded DNA (dsDNA) is associated with many important biological processes such as DNA-protein recognition and DNA packaging into nucleosomes. Thermodynamics of dsDNA bending has been studied by a method called cyclization which relies on DNA ligase to covalently join short sticky ends of a dsDNA. However, ligation efficiency can be affected by many factors that are not related to dsDNA looping such as the DNA structure surrounding the joined sticky ends, and ligase can also affect the apparent looping rate through mechanisms such as nonspecific binding. Here, we show how to measure dsDNA looping kinetics without ligase by detecting transient DNA loop formation by FRET (Fluorescence Resonance Energy Transfer). dsDNA molecules are constructed using a simple PCR-based protocol with a FRET pair and a biotin linker. The looping probability density known as the J factor is extracted from the looping rate and the annealing rate between two disconnected sticky ends. By testing two dsDNAs with different intrinsic curvatures, we show that the J factor is sensitive to the intrinsic shape of the dsDNA.

Measuring shape-dependent looping probability of DNA

Recently, several studies have shown that short doubled-stranded DNA (dsDNA) loops more readily than the wormlike chain model predicts. In most of these experiments, the intrinsic bendedness of dsDNA, which in theory can dramatically influence looping dynamics, was either avoided or unaccounted for. To investigate the effect of the shape of dsDNA on looping dynamics, we characterized the shapes of several synthetic dsDNA molecules of equal length but different sequences using gel electrophoresis. We then measured their looping rates using a FRET (Förster resonance energy transfer)-based assay and extracted the looping probability density known as the J factor (jM). We also used, for comparison, several dinucleotide angular parameter sets derived from the observed electrophoretic mobility to compute the jM predicted by the wormlike chain model. Although we found a strong correlation between curvature and jM, the measured jM was higher than most dinucleotide model predictions. This result suggests that it is difficult to reconcile the looping probability with the observed gel mobility within the wormlike chain model and underscores the importance of determining the intrinsic shape of dsDNA for proper theoretical analysis.

The effect of an intervening promoter nucleosome on gene expression

Nucleosomes, which are the basic packaging units of chromatin, are stably positioned in promoters upstream of most stress-inducible genes. These promoter nucleosomes are generally thought to repress gene expression due to exclusion; they prevent transcription factors from accessing their target sites on the DNA. However, the role of promoter nucleosomes that do not directly occlude transcription factor binding sites is not obvious. Here, we varied the stability of a non-occluding nucleosome positioned between a transcription factor binding site and the TATA box region in an inducible yeast promoter and measured downstream gene expression level. We found that gene expression level depends on the occupancy of the non-occluding nucleosome in a non-monotonic manner. We postulated that a non-occluding nucleosome can serve both as a vehicle of and a barrier to chromatin remodeling activity and built a quantitative, nonequilibrium model to explain the observed nontrivial effect of the intervening nucleosome. Our work sheds light on the dual role of nucleosome as a repressor and an activator and expands the standard model of gene expression to include irreversible promoter chromatin transitions.

Immunological impact of Taekwondo competitions

Immunological changes in elite adolescent female athletes during Taekwondo competitions were investigated on-field. 6 female athletes (16.7 ± 0.8 year-old) volunteered and performed 5 bouts of demonstration Taekwondo competitions simulating real tournaments in intensity, duration, and break-time intervals on the same day. Blood samples were taken before, after the competitions and during the recovery, respectively. Immunological changes and oxidative stress in peripheral blood mononuclear cells were evaluated by flow-cytometry. During the competitions, exercise intensity was 92.2 ± 3.8% (86.1~95.7) of the maximal heart rate. Blood lactate increased immediately after the competitions (p=0.0165) and decreased to baseline during recovery. Intracellular reactive oxygen species (ROS) in the peripheral blood increased continuously during recovery (p<0.05, respectively). Natural killer cells increased immediately after the competitions (p=0.0006), and decreased during recovery. B and T cells increased immediately after the competitions and remained elevated throughout recovery (p<0.05, respectively). CD4/CD8 ratio after the competitions was decreased (p=0.0091) and returned to baseline during recovery. These results suggest that the immunological function of the elite female adolescent athletes could be attenuated after Taekwondo competitions. Further large-scaled Taekwondo studies on immunologic and apoptotic changes related to oxidative stress should be performed for improving and protecting the health of adolescent athletes.

Transcriptional regulatory circuits: predicting numbers from alphabets

Transcriptional regulatory circuits govern how cis and trans factors transform signals into messenger RNA (mRNA) expression levels. With advances in quantitative and high-throughput technologies that allow measurement of gene expression state in different conditions, data that can be used to build and test models of transcriptional regulation is being generated at a rapid pace. Here, we review experimental and computational methods used to derive detailed quantitative circuit models on a small scale and cruder, genome-wide models on a large scale. We discuss the potential of combining small- and large-scale approaches to understand the working and wiring of transcriptional regulatory circuits.

Efficacy and prognosis of a short course of prednisolone therapy for pediatric epilepsy

PURPOSE

To evaluate the efficacy and safety of adjunctive prednisolone therapy in children with cryptogenic epileptic encephalopathy, other than infantile spasms, and to determine its prognosis.

METHODS

Prednisolone, 2mg/kg per day for 6 weeks, tapered for a further 2 weeks, was given in combination with previously prescribed antiepileptic drugs. A retrospective assessment of 41 children thus treated included measurements of seizure frequency, electroencephalographic findings, global assessments of cognitive function, and adverse drug events. Long-term patient prognoses over a mean follow-up period of 3 years and 5 months (range, 14-90 months) were also examined.

RESULTS

Of 41 patients, 32 had Lennox-Gastaut syndrome, 4 had Doose syndrome, 1 had Otahara syndrome, 2 had Landau-Kleffner syndrome, and 2 had other unspecified generalized epilepsies. After prednisolone therapy, 73% (30/41) of patients showed a reduction in seizure frequency of >50%, and 59% (24/41) became seizure free. However, only seven patients (four with Lennox-Gastaut syndrome, two with Doose syndrome, and one with unspecified generalized epilepsy) who became seizure free remained free of seizures at the time of the final follow-up. Electroencephalographic findings and global assessments of cognitive function correlated well with seizure outcomes. No significant demographic factors influenced the efficacy of prednisolone or patient prognoses after prednisolone tapering. Most adverse events were transient, or were tolerated well with conservative management, with maintenance of the medication.

CONCLUSION

Prednisolone therapy may be a safe and effective adjunct in patients with cryptogenic epileptic encephalopathies, but the high relapse rate is of concern.

The role of fluctuations in tRNA selection by the ribosome

The detailed mechanism of how the ribosome decodes protein sequence information with an abnormally high accuracy, after 40 years of study, remains elusive. A critical element in selecting correct transfer RNA (tRNA) transferring correct amino acid is "induced fit" between the ribosome and tRNA. By using single-molecule methods, the induced fit mechanism is shown to position favorably the correct tRNA after initial codon recognition. We provide evidence that this difference in positioning and thermal fluctuations constitutes the primary mechanism for the initial selection of tRNA. This work demonstrates thermal fluctuations playing a critical role in the substrate selection by an enzyme.

Fluctuations of transfer RNAs between classical and hybrid states

Adjacent transfer RNAs (tRNAs) in the A- and P-sites of the ribosome are in dynamic equilibrium between two different conformations called classical and hybrid states before translocation. Here, we have used single-molecule fluorescence resonance energy transfer to study the effect of Mg(2+) on tRNA dynamics with and without an acetyl group on the A-site tRNA. When the A-site tRNA is not acetylated, tRNA dynamics do not depend on [Mg(2+)], indicating that the relative positions of the substrates for peptide-bond formation are not affected by Mg(2+). In sharp contrast, when the A-site tRNA is acetylated, Mg(2+) lengthens the lifetime of the classical state but does not change the lifetime of the hybrid state. Based on these findings, the classical state resembles a state with direct stabilization of tertiary structure by Mg(2+) ions whereas the hybrid state resembles a state with little Mg(2+)-assisted stabilization. The antibiotic viomycin, a translocation inhibitor, suppresses tRNA dynamics, suggesting that the enhanced fluctuations of tRNAs after peptide-bond formation drive spontaneous attempts at translocation by the ribosome.

Ion atmosphere of three-way junction nucleic acid

The ion atmosphere of three-armed symmetric Y-shaped and asymmetric y-shaped A-RNA junctions in aqueous solution containing multivalent ions is described within the framework of a polyelectrolyte model. The fraction of "screening counterions" per polyion charge that shield the residual unneutralized charges from interacting with one another and the condensed counterions per polyion charge as a function of sodium and magnesium ion concentrations are determined. The predictions for the slope of log(k(o)/k(f)) as a function of Na+ and Mg2+ concentration, where k(o) and k(f) are the opening and folding rates of the three-helix junction molecule, respectively, are compared with experimental data (Kim et al. Proc. Nat. Acad. Sci. U.S.A. 2002, 96, 9077-9082).

Peptide bond formation destabilizes Shine–Dalgarno interaction on the ribosome

The ribosome is a molecular machine that translates the genetic code contained in the messenger RNA into an amino acid sequence through repetitive cycles of transfer RNA selection, peptide bond formation and translocation. Here we demonstrate an optical tweezer assay to measure the rupture force between a single ribosome complex and mRNA. The rupture force was compared between ribosome complexes assembled on an mRNA with and without a strong Shine-Dalgarno (SD) sequence-a sequence found just upstream of the coding region of bacterial mRNAs, involved in translation initiation. The removal of the SD sequence significantly reduced the rupture force in complexes carrying an aminoacyl tRNA, Phe-tRNA(Phe), in the A site, indicating that the SD interactions contribute significantly to the stability of the ribosomal complex on the mRNA before peptide bond formation. In contrast, the presence of a peptidyl tRNA analogue, N-acetyl-Phe-tRNA(Phe), in the A site, which mimicked the post-peptidyl transfer state, weakened the rupture force as compared to the complex with Phe-tRNA(Phe), and the resultant force was the same for both the SD-containing and SD-deficient mRNAs. These results suggest that formation of the first peptide bond destabilizes the SD interaction, resulting in the weakening of the force with which the ribosome grips an mRNA. This might be an important requirement to facilitate movement of the ribosome along mRNA during the first translocation step.

Vagus nerve stimulation in pediatric intractable epilepsy: a Korean bicentric study

OBJECTIVE

To present our experience with vagus nerve stimulation (VNS) and to evaluate the long-term efficacy and safety of the procedure in pediatric intractable epilepsy.

METHODS

This study included sixteen patients, who were implanted with a vagus nerve stimulator and could be followed up for at least more than 12 months in two epilepsy centers. Data including seizure frequency, EEG, quality of life measures and adverse events were prospectively filed over a 5-year period.

RESULTS

VNS resulted in a > 50% reduction in seizure frequency in 50.0% (8/16) of children with 31.3% (5/16) of patients achieving a > 90% reduction. Additionally, enhancements in quality of life were as follows: memory in 50.0% (8/16), mood in 62.5% (10/16), behavior in 68.8% (11/16), alertness in 68.8% (11/16), achievement in 37.5% (6/16), and verbal skills in 43.8% (7/16) of the patients. Adverse events included hoarseness in two patients, dyspnea during sleep in two patients, and sialorrhea in one patient. However, these events were tolerable or could be controlled by the adjustment of output currents. In one patient, wound revision was required.

CONCLUSION

Our data supports the role of VNS as an alternative therapy for pediatric intractable epilepsy.