A backbone‐dependent rotamer library with high (ϕ, ψ) coverage using metadynamics simulations

Backbone‐dependent rotamer libraries are commonly used to assign the side chain dihedral angles of amino acids when modeling protein structures. Most rotamer libraries are created by curating protein crystal structure data and using various methods to extrapolate the existing data to cover all possible backbone conformations. However, these rotamer libraries may not be suitable for modeling the structures of cyclic peptides and other constrained peptides because these molecules frequently sample backbone conformations rarely seen in the crystal structures of linear proteins. To provide backbone‐dependent side chain information beyond the α‐helix, β‐sheet, and PPII regions, we used explicit‐solvent metadynamics simulations of model dipeptides to create a new rotamer library that has high coverage in the (ϕ, ψ) space. Furthermore, this approach can be applied to build high‐coverage rotamer libraries for noncanonical amino acids. The resulting Metadynamics of Dipeptides for Rotamer Distribution (MEDFORD) rotamer library predicts the side chain conformations of high‐resolution protein crystal structures with similar accuracy (~80%) to a state‐of‐the‐art rotamer library. Our ability to test the accuracy of MEDFORD at predicting the side chain dihedral angles of amino acids in noncanonical backbone conformation is restricted by the limited structural data available for cyclic peptides. For the cyclic peptide data that are currently available, MEDFORD and the state‐of‐the‐art rotamer library perform comparably. However, the two rotamer libraries indeed make different rotamer predictions in noncanonical (ϕ, ψ) regions. For noncanonical amino acids, the MEDFORD rotamer library predicts the χ ₁ values with approximately 75% accuracy.     

Alpha/gamma transitions in the B-DNA backbone

In the crystal structures of protein complexes with B-DNA, α and γ DNA backbone torsion angles often exhibit non-canonical values. It is not known if these alternative backbone conformations are easily accessible in solution and can contribute to the specific recognition of DNA by proteins. We have analysed the coupled transition of the α and γ torsion angles within the central GpC step of a B-DNA dodecamer by computer simulations. Five stable or metastable non-canonical α/γ sub-states are found. The most favourable pathway from the canonical α/γ structure to any unusual form involves a counter-rotation of α and γ, via the trans conformation. However, the corresponding free energy indicates that spontaneous flipping of the torsions is improbable in free B-DNA. This is supported by an analysis of the available high resolution crystallographic structures showing that unusual α/γ states are only encountered in B-DNA complexed to proteins. An analysis of the structural consequences of α/γ transitions shows that the non-canonical backbone geometry influences essentially the roll and twist values and reduces the equilibrium dispersion of structural parameters. Our results support the hypothesis that unusual α/γ backbones arise during protein–DNA complexation, assisting the fine structural adjustments between the two partners and playing a role in the overall complexation free energy.     

Highly Mutagenic Exocyclic DNA Adducts Are Substrates for the Human Nucleotide Incision Repair Pathway

Background

Oxygen free radicals induce lipid peroxidation (LPO) that damages and breaks polyunsaturated fatty acids in cell membranes. LPO-derived aldehydes and hydroxyalkenals react with DNA leading to the formation of etheno(ε)-bases including 1,N ⁶-ethenoadenine (εA) and 3,N ⁴-ethenocytosine (εC). The εA and εC residues are highly mutagenic in mammalian cells and eliminated in the base excision repair (BER) pathway and/or by AlkB family proteins in the direct damage reversal process. BER initiated by DNA glycosylases is thought to be the major pathway for the removal of non-bulky endogenous base damage. Alternatively, in the nucleotide incision repair (NIR) pathway, the apurinic/apyrimidinic (AP) endonucleases can directly incise DNA duplex 5′ to a damaged base in a DNA glycosylase-independent manner.

Methodology/Principal Findings

Here we have characterized the substrate specificity of human major AP endonuclease 1, APE1, towards εA, εC, thymine glycol (Tg) and 7,8-dihydro-8-oxoguanine (8oxoG) residues when present in duplex DNA. APE1 cleaves oligonucleotide duplexes containing εA, εC and Tg, but not those containing 8oxoG. Activity depends strongly on sequence context. The apparent kinetic parameters of the reactions suggest that APE1 has a high affinity for DNA containing ε-bases but cleaves DNA duplexes at an extremely slow rate. Consistent with this observation, oligonucleotide duplexes containing an ε-base strongly inhibit AP site nicking activity of APE1 with IC₅₀ values in the range of 5–10 nM. MALDI-TOF MS analysis of the reaction products demonstrated that APE1-catalyzed cleavage of εA•T and εC•G duplexes generates, as expected, DNA fragments containing 5′-terminal ε-base residue.

Conclusions/Significance

The fact that ε-bases and Tg in duplex DNA are recognized and cleaved by APE1 in vitro, suggests that NIR may act as a backup pathway to BER to remove a large variety of genotoxic base lesions in human cells.     

Intrinsic α‐helical and β‐sheet conformational preferences: A computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle (-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (N–C_α–C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the “bridge” region of-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between -ψ conformations and τ.     

Interactions of L-3,5,3'-Triiodothyronine,Allopregnanolone, and Ivermectin with the GABAA Receptor: Evidence for Overlapping Intersubunit Binding Modes

Structural mechanisms of modulation of γ-aminobutyric acid (GABA) type A receptors by neurosteroids and hormones remain unclear. The thyroid hormone L-3,5,3’-triiodothyronine (T3) inhibits GABA_A receptors at micromolar concentrations and has common features with neurosteroids such as allopregnanolone (ALLOP). Here we use functional experiments on α₂β₁γ₂ GABA_A receptors expressed in Xenopus oocytes to detect competitive interactions between T3 and an agonist (ivermectin, IVM) with a crystallographically determined binding site at subunit interfaces in the transmembrane domain of a homologous receptor (glutamate-gated chloride channel, GluCl). T3 and ALLOP also show competitive effects, supporting the presence of both a T3 and ALLOP binding site at one or more subunit interfaces. Molecular dynamics (MD) simulations over 200 ns are used to investigate the dynamics and energetics of T3 in the identified intersubunit sites. In these simulations, T3 molecules occupying all intersubunit sites (with the exception of the α-β interface) display numerous energetically favorable conformations with multiple hydrogen bonding partners, including previously implicated polar/acidic sidechains and a structurally conserved deformation in the M1 backbone.     

A maximum entropy approach to the study of residue-specific backbone angle distributions in α-synuclein,an intrinsically disordered protein

α-Synuclein is an intrinsically disordered protein of 140 residues that switches to an α-helical conformation upon binding phospholipid membranes. We characterize its residue-specific backbone structure in free solution with a novel maximum entropy procedure that integrates an extensive set of NMR data. These data include intraresidue and sequential H^N–H^α and H^N–H^N NOEs, values for ³J_HNHα, ¹J_HαCα, ²J_CαN, and ¹J_CαN, as well as chemical shifts of ¹⁵N, ¹³C^α, and ¹³C′ nuclei, which are sensitive to backbone torsion angles. Distributions of these torsion angles were identified that yield best agreement to the experimental data, while using an entropy term to minimize the deviation from statistical distributions seen in a large protein coil library. Results indicate that although at the individual residue level considerable deviations from the coil library distribution are seen, on average the fitted distributions agree fairly well with this library, yielding a moderate population (20–30%) of the PP_II region and a somewhat higher population of the potentially aggregation-prone β region (20–40%) than seen in the database. A generally lower population of the α_R region (10–20%) is found. Analysis of ¹H–¹H NOE data required consideration of the considerable backbone diffusion anisotropy of a disordered protein.     

Roles of Residues Arg-61 and Gln-38 of Human DNA Polymerase η in Bypass of Deoxyguanosine and 7,8-Dihydro-8-oxo-2′-deoxyguanosine

Like the other Y-family DNA polymerases, human DNA polymerase η (hpol η) has relatively low fidelity and is able to tolerate damage during DNA synthesis, including 7,8-dihydro-8-oxo-2′-deoxyguanosine (8-oxoG), one of the most abundant DNA lesions in the genome. Crystal structures show that Arg-61 and Gln-38 are located near the active site and may play important roles in the fidelity and efficiency of hpol η. Site-directed mutagenesis was used to replace these side chains either alone or together, and the wild type or mutant proteins were purified and tested by replicating DNA past deoxyguanosine (G) or 8-oxoG. The catalytic activity of hpol η was dramatically disrupted by the R61M and Q38A/R61A mutations, as opposed to the R61A and Q38A single mutants. Crystal structures of hpol η mutant ternary complexes reveal that polarized water molecules can mimic and partially compensate for the missing side chains of Arg-61 and Gln-38 in the Q38A/R61A mutant. The combined data indicate that the positioning and positive charge of Arg-61 synergistically contribute to the nucleotidyl transfer reaction, with additional influence exerted by Gln-38. In addition, gel filtration chromatography separated multimeric and monomeric forms of wild type and mutant hpol η, indicating the possibility that hpol η forms multimers in vivo.     

The efficiency and fidelity of 8-oxo-guanine bypass by DNA polymerases δ and η

A DNA lesion created by oxidative stress is 7,8-dihydro-8-oxo-guanine (8-oxoG). Because 8-oxoG can mispair with adenine during DNA synthesis, it is of interest to understand the efficiency and fidelity of 8-oxoG bypass by DNA polymerases. We quantify bypass parameters for two DNA polymerases implicated in 8-oxoG bypass, Pols δ and η. Yeast Pol δ and yeast Pol η both bypass 8-oxoG and misincorporate adenine during bypass. However, yeast Pol η is 10-fold more efficient than Pol δ, and following bypass Pol η switches to less processive synthesis, similar to that observed during bypass of a cis-syn thymine-thymine dimer. Moreover, yeast Pol η is at least 10-fold more accurate than yeast Pol δ during 8-oxoG bypass. These differences are maintained in the presence of the accessory proteins RFC, PCNA and RPA and are consistent with the established role of Pol η in suppressing ogg1-dependent mutagenesis in yeast. Surprisingly different results are obtained with human and mouse Pol η. Both mammalian enzymes bypass 8-oxoG efficiently, but they do so less processively, without a switch point and with much lower fidelity than yeast Pol η. The fact that yeast and mammalian Pol η have intrinsically different catalytic properties has potential biological implications.     

New Insights into the Interdependence between Amino Acid Stereochemistry and Protein Structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ϕ and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.     

Solution structure of an oligonucleotide containing an abasic site: evidence for an unusual deoxyribose conformation

The antitumor antibiotic bleomycin causes two major lesions in the deoxyribose backbone of DNA: formation of 4′-keto abasic sites and formation of strand breaks with 3′-phosphoglycolate and 5′-phosphate ends. As a model for the 4′-keto abasic site, we have characterized an abasic site (X) in d(CCAAAGXACTGGG)·d(CCCAGTACTTTGG) by two-dimensional NMR spectroscopy. A total of 475 NOEs and 101 dihedral angles provided the restraints for molecular modeling. Four unusual NOEs were observed between each anomer of the abasic site and the neighboring bases. In addition, four coupling constants for adjacent protons of the deoxyribose of both the α and β anomers of the abasic site were observed. The modeling suggests that for both anomers the abasic site is extrahelical, without significant distortion of the backbone opposite the lesion. The coupling constants further allowed assignment of an unusual sugar pucker for each anomer. The unique position of the abasic site in our structural model for each anomer is discussed in terms of repair of such lesions in vivo.     

Evidence for lesion bypass by yeast replicative DNA polymerases during DNA damage

The enzyme ribonucleotide reductase, responsible for the synthesis of deoxyribonucleotides (dNTP), is upregulated in response to DNA damage in all organisms. In Saccharomyces cerevisiae, dNTP concentration increases ~6- to 8-fold in response to DNA damage. This concentration increase is associated with improved tolerance of DNA damage, suggesting that translesion DNA synthesis is more efficient at elevated dNTP concentration. Here we show that in a yeast strain with all specialized translesion DNA polymerases deleted, 4-nitroquinoline oxide (4-NQO) treatment increases mutation frequency ~3-fold, and that an increase in dNTP concentration significantly improves the tolerance of this strain to 4-NQO induced damage. In vitro, under single-hit conditions, the replicative DNA polymerase ε does not bypass 7,8-dihydro-8-oxoguanine lesion (8-oxoG, one of the lesions produced by 4-NQO) at S-phase dNTP concentration, but does bypass the same lesion with 19–27% efficiency at DNA-damage-state dNTP concentration. The nucleotide inserted opposite 8-oxoG is dATP. We propose that during DNA damage in S. cerevisiae increased dNTP concentration allows replicative DNA polymerases to bypass certain DNA lesions.     

Molecular Basis of the Membrane Interaction of the β2e Subunit of Voltage-Gated Ca2+ Channels

The auxiliary β subunit plays an important role in the regulation of voltage-gated calcium (Ca_V) channels. Recently, it was revealed that β2e associates with the plasma membrane through an electrostatic interaction between N-terminal basic residues and anionic phospholipids. However, a molecular-level understanding of β-subunit membrane recruitment in structural detail has remained elusive. In this study, using a combination of site-directed mutagenesis, liposome-binding assays, and multiscale molecular-dynamics (MD) simulation, we developed a physical model of how the β2e subunit is recruited electrostatically to the plasma membrane. In a fluorescence resonance energy transfer assay with liposomes, binding of the N-terminal peptide (23 residues) to liposome was significantly increased in the presence of phosphatidylserine (PS) and phosphatidylinositol 4,5-bisphosphate (PIP₂). A mutagenesis analysis suggested that two basic residues proximal to Met-1, Lys-2 (K2) and Trp-5 (W5), are more important for membrane binding of the β2e subunit than distal residues from the N-terminus. Our MD simulations revealed that a stretched binding mode of the N-terminus to PS is required for stable membrane attachment through polar and nonpolar interactions. This mode obtained from MD simulations is consistent with experimental results showing that K2A, W5A, and K2A/W5A mutants failed to be targeted to the plasma membrane. We also investigated the effects of a mutated β2e subunit on inactivation kinetics and regulation of Ca_V channels by PIP₂. In experiments with voltage-sensing phosphatase (VSP), a double mutation in the N-terminus of β2e (K2A/W5A) increased the PIP₂ sensitivity of Ca_V2.2 and Ca_V1.3 channels by ∼3-fold compared with wild-type β2e subunit. Together, our results suggest that membrane targeting of the β2e subunit is initiated from the nonspecific electrostatic insertion of N-terminal K2 and W5 residues into the membrane. The PS-β2e interaction observed here provides a molecular insight into general principles for protein binding to the plasma membrane, as well as the regulatory roles of phospholipids in transporters and ion channels.     

High-Resolution Crystal Structures of Protein Helices Reconciled with Three-Centered Hydrogen Bonds and Multipole Electrostatics

Theoretical and experimental evidence for non-linear hydrogen bonds in protein helices is ubiquitous. In particular, amide three-centered hydrogen bonds are common features of helices in high-resolution crystal structures of proteins. These high-resolution structures (1.0 to 1.5 Å nominal crystallographic resolution) position backbone atoms without significant bias from modeling constraints and identify Φ = -62°, ψ = -43 as the consensus backbone torsional angles of protein helices. These torsional angles preserve the atomic positions of α-β carbons of the classic Pauling α-helix while allowing the amide carbonyls to form bifurcated hydrogen bonds as first suggested by Némethy et al. in 1967. Molecular dynamics simulations of a capped 12-residue oligoalanine in water with AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications), a second-generation force field that includes multipole electrostatics and polarizability, reproduces the experimentally observed high-resolution helical conformation and correctly reorients the amide-bond carbonyls into bifurcated hydrogen bonds. This simple modification of backbone torsional angles reconciles experimental and theoretical views to provide a unified view of amide three-centered hydrogen bonds as crucial components of protein helices. The reason why they have been overlooked by structural biologists depends on the small crankshaft-like changes in orientation of the amide bond that allows maintenance of the overall helical parameters (helix pitch (p) and residues per turn (n)). The Pauling 3.6₁₃ α-helix fits the high-resolution experimental data with the minor exception of the amide-carbonyl electron density, but the previously associated backbone torsional angles (Φ, Ψ) needed slight modification to be reconciled with three-atom centered H-bonds and multipole electrostatics. Thus, a new standard helix, the 3.6_13/10-, Némethy- or N-helix, is proposed. Due to the use of constraints from monopole force fields and assumed secondary structures used in low-resolution refinement of electron density of proteins, such structures in the PDB often show linear hydrogen bonding.     

In Vitro and In Vivo Oxidation of Methionine Residues in Small, Acid-Soluble Spore Proteins from Bacillus Species

Methionine residues in α/β-type small, acid-soluble spore proteins (SASP) of Bacillus species were readily oxidized to methionine sulfoxide in vitro by t-butyl hydroperoxide (tBHP) or hydrogen peroxide (H₂O₂). These oxidized α/β-type SASP no longer bound to DNA effectively, but DNA binding protected α/β-type SASP against methionine oxidation by peroxides in vitro. Incubation of an oxidized α/β-type SASP with peptidyl methionine sulfoxide reductase (MsrA), which can reduce methionine sulfoxide residues back to methionine, restored the α/β-type SASP’s ability to bind to DNA. Both tBHP and H₂O₂ caused some oxidation of the two methionine residues of an α/β-type SASP (SspC) in spores of Bacillus subtilis, although one methionine which is highly conserved in α/β-type SASP was only oxidized to a small degree. However, much more methionine sulfoxide was generated by peroxide treatment of spores carrying a mutant form of SspC which has a lower affinity for DNA. MsrA activity was present in wild-type B. subtilis spores. However, msrA mutant spores were no more sensitive to H₂O₂ than were wild-type spores. The major mechanism operating for dealing with oxidative damage to α/β-type SASP in spores is DNA binding, which protects the protein’s methionine residues from oxidation both in vitro and in vivo. This may be important in vivo since α/β-type SASP containing oxidized methionine residues no longer bind DNA well and α/β-type SASP-DNA binding is essential for long-term spore survival.     

7,8-Dihydro-8-oxo-1,N6-ethenoadenine: an exclusively Hoogsteen-paired thymine mimic in DNA that induces A→T transversions in Escherichia coli

This work investigated the structural and biological properties of DNA containing 7,8-dihydro-8-oxo-1,N⁶-ethenoadenine (oxo-ϵA), a non-natural synthetic base that combines structural features of two naturally occurring DNA lesions (7,8-dihydro-8-oxoadenine and 1,N⁶-ethenoadenine). UV-, CD-, NMR spectroscopies and molecular modeling of DNA duplexes revealed that oxo-ϵA adopts the non-canonical syn conformation (χ = 65º) and fits very well among surrounding residues without inducing major distortions in local helical architecture. The adduct remarkably mimics the natural base thymine. When considered as an adenine-derived DNA lesion, oxo-ϵA was >99% mutagenic in living cells, causing predominantly A→T transversion mutations in Escherichia coli. The adduct in a single-stranded vector was not repaired by base excision repair enzymes (MutM and MutY glycosylases) or the AlkB dioxygenase and did not detectably affect the efficacy of DNA replication in vivo. When the biological and structural data are viewed together, it is likely that the nearly exclusive syn conformation and thymine mimicry of oxo-ϵA defines the selectivity of base pairing in vitro and in vivo, resulting in lesion pairing with A during replication. The base pairing properties of oxo-ϵA, its strong fluorescence and its invisibility to enzymatic repair systems in vivo are features that are sought in novel DNA-based probes and modulators of gene expression.     

The Werner syndrome protein limits the error-prone 8-oxo-dG lesion bypass activity of human DNA polymerase kappa

Human DNA polymerase kappa (hpol κ) is the only Y-family member to preferentially insert dAMP opposite 7,8-dihydro-8-oxo-2′-deoxyguanosine (8-oxo-dG) during translesion DNA synthesis. We have studied the mechanism of action by which hpol κ activity is modulated by the Werner syndrome protein (WRN), a RecQ helicase known to influence repair of 8-oxo-dG. Here we show that WRN stimulates the 8-oxo-dG bypass activity of hpol κ in vitro by enhancing the correct base insertion opposite the lesion, as well as extension from dC:8-oxo-dG base pairs. Steady-state kinetic analysis reveals that WRN improves hpol κ-catalyzed dCMP insertion opposite 8-oxo-dG ∼10-fold and extension from dC:8-oxo-dG by 2.4-fold. Stimulation is primarily due to an increase in the rate constant for polymerization (k_pol), as assessed by pre-steady-state kinetics, and it requires the RecQ C-terminal (RQC) domain. In support of the functional data, recombinant WRN and hpol κ were found to physically interact through the exo and RQC domains of WRN, and co-localization of WRN and hpol κ was observed in human cells treated with hydrogen peroxide. Thus, WRN limits the error-prone bypass of 8-oxo-dG by hpol κ, which could influence the sensitivity to oxidative damage that has previously been observed for Werner''s syndrome cells.     

Probing the Effect of the Non-Active-Site Mutation Y229W in New Delhi Metallo-β-lactamase-1 by Site-Directed Mutagenesis,Kinetic Studies,and Molecular Dynamics Simulations

New Delhi metallo-β-lactmase-1 (NDM-1) has attracted extensive attention for its high catalytic activities of hydrolyzing almost all β-lactam antibiotics. NDM-1 shows relatively higher similarity to subclass B1 metallo-β-lactmases (MβLs), but its residue at position 229 is identical to that of B2/B3 MβLs, which is a Tyr instead of a B1-MβL-conserved Trp. To elucidate the possible role of Y229 in the bioactivity of NDM-1, we performed mutagenesis study and molecular dynamics (MD) simulations. Although residue Y229 is spatially distant from the active site and not contacting directly with the substrate or zinc ions, the Y229W mutant was found to have higher k_cat and K_m values than those of wild-type NDM-1, resulting in 1∼7 fold increases in k_cat/K_m values against tested antibiotics. In addition, our MD simulations illustrated the enhanced flexibility of Loop 2 upon Y229W mutation, which could increase the kinetics of both substrate entrance (kon) and product egress (koff). The enhanced flexibility of Loop 2 might allow the enzyme to adjust the geometry of its active site to accommodate substrates with different structures, broadening its substrate spectrum. This study indicated the possible role of the residue at position 229 in the evolution of NDM-1.     

Role of the Nfo and ExoA apurinic/apyrimidinic endonucleases in repair of DNA damage during outgrowth of Bacillus subtilis spores

Germination and outgrowth are critical steps for returning Bacillus subtilis spores to life. However, oxidative stress due to full hydration of the spore core during germination and activation of metabolism in spore outgrowth may generate oxidative DNA damage that in many species is processed by apurinic/apyrimidinic (AP) endonucleases. B. subtilis spores possess two AP endonucleases, Nfo and ExoA; the outgrowth of spores lacking both of these enzymes was slowed, and the spores had an elevated mutation frequency, suggesting that these enzymes repair DNA lesions induced by oxidative stress during spore germination and outgrowth. Addition of H₂O₂ also slowed the outgrowth of nfo exoA spores and increased the mutation frequency, and nfo and exoA mutations slowed the outgrowth of spores deficient in either RecA, nucleotide excision repair (NER), or the DNA-protective α/β-type small acid-soluble spore proteins (SASP). These results suggest that α/β-type SASP protect DNA of germinating spores against damage that can be repaired by Nfo and ExoA, which is generated either spontaneously or promoted by addition of H₂O₂. The contribution of RecA and Nfo/ExoA was similar to but greater than that of NER in repair of DNA damage generated during spore germination and outgrowth. However, nfo and exoA mutations increased the spontaneous mutation frequencies of outgrown spores lacking uvrA or recA to about the same extent, suggesting that DNA lesions generated during spore germination and outgrowth are processed by Nfo/ExoA in combination with NER and/or RecA. These results suggest that Nfo/ExoA, RecA, the NER system, and α/β-type SASP all contribute to the repair of and/or protection against oxidative damage of DNA in germinating and outgrowing spores.     

Entropy,Energy, and Bending of DNA in Viral Capsids

Inspired by novel single-molecule and bulk solution measurements, the physics underlying the forces and pressures involved in DNA packaging into bacteriophage capsids became the focus of numerous recent theoretical models. These fall into two general categories: Continuum-elastic theories (CT), and simulation studies—mostly of the molecular dynamics (MD) genre. Both types of models account for the dependence of the force, and hence the packaging free energy (ΔF), on the loaded DNA length, but differ markedly in interpreting their origin. While DNA confinement entropy is a dominant contribution to ΔF in the MD simulations, in the CT theories this role is fulfilled by interstrand repulsion, and there is no explicit entropy term. The goal of this letter is to resolve this apparent contradiction, elucidate the origin of the entropic term in the MD simulations, and point out its tacit presence in the CT treatments.The genomic double-stranded (ds) DNA inside bacteriophage heads is highly stressed, leading to internal pressures of up to ∼50 atmospheres, reflecting the tight packing and extreme bending of this highly charged and rigid molecule (1). The interaxial distance (d) between neighboring (nonbonded) dsDNA segments in the fully packaged virus is typically ≈2.5 nm (2,3), just slightly larger than the hardcore diameter of dsDNA (b = 2.0 nm) and well into the repulsive regime (d ≤ 2.8 nm) of DNA-DNA interaction in ionic solutions (4–6). Moreover, free dsDNA in (physiological) solution is a fluctuating, semiflexible, wormlike chain (WLC), with persistence length ξ ≈ 50 nm, larger than the radius of most viral capsids. Thus, on a molecular scale, packaging the long (e.g., the 330-ξ long λ-phage genome) viral DNA into its tiny capsid requires enormous mechanical work.The force needed to package the DNA is provided by an ATP-driven motor protein situated at the capsid portal. Recent single molecule measurements reveal that this force, f(L_int), increases sharply with the loaded genome length, L_int, rising to ∼30–100 pN, depending on the virus in question (7,8). These studies inspired the formulation of many theoretical models of DNA packaging in viral capsids, which fall roughly into two categories: