Kinking the double helix by bending deformation

DNA bending and torsional deformations that often occur during its functioning inside the cell can cause local disruptions of the regular helical structure. The disruptions created by negative torsional stress have been studied in detail, but those caused by bending stress have only been analyzed theoretically. By probing the structure of very small DNA circles, we determined that bending stress disrupts the regular helical structure when the radius of DNA curvature is smaller than 3.5 nm. First, we developed an efficient method to obtain covalently closed DNA minicircles. To detect structural disruptions in the minicircles we treated them by single-strand-specific endonucleases. The data showed that the regular DNA structure is disrupted by bending deformation in the 64–65-bp minicircles, but not in the 85–86-bp minicircles. Our results suggest that strong DNA bending initiates kink formation while preserving base pairing.     

Investigation of flap flexibility of β-secretase using molecular dynamic simulations

Flap motif and its dynamics were extensively reported in aspartate proteases, e.g. HIV proteases and plasmepsins. Herein, we report the first account of flap dynamics amongst different conformations of β-secretase using molecular dynamics simulation. Various parameters were proposed and a selected few were picked which could appropriately describe the flap motion. Three systems were studied, namely Free (BACE_Free) and two ligand-bound conformations, which belonged to space groups P6₁22 (BACE_Bound1) and C222₁ (BACE_Bound2), respectively and four parameters (distance between the flaps tip residue, Thr72 and Ser325, d₁; dihedral angle, ? (Thr72-Asp32-Asp228-Ser325); TriCα angles, θ₁ (Thr72-Asp32-Ser325), and θ₂ (Thr72-Asp228-Ser325)) were proposed to understand the change in dynamics of flap domain and the extent of flap opening and closing. Analysis of, θ₂, d₁, θ₁ and ? confirmed that the BACE_Free adopted semi-open, open and closed conformations with slight twisting during flap opening. However, BACE_Bound1 (P6122) showed an adaptation to open conformation due to lack of hydrogen bond interaction between the ligand and flap tip residue. A slight flap twisting, ? (lateral twisting) was observed for BACE_Bound1 during flap opening which correlates with the opening of BACE_Free. Contradictory to the BACE_Bound1, the BACE_Bound2 locked the flap in a closed conformation throughout the simulation due to formation of a stable hydrogen bond interaction between the flap tip residue and ligand. Analyses of all three systems highlight that d₁, θ₂ and ? can be precisely used to describe the extent of flap opening and closing concurrently with snapshots along the molecular dynamics trajectory across several conformations of β-secretase.     

Nucleotide excision repair of 2-acetylaminofluorene- and 2-aminofluorene-(C8)-guanine adducts: molecular dynamics simulations elucidate how lesion structure and base sequence context impact repair efficiencies

Nucleotide excision repair (NER) efficiencies of DNA lesions can vary by orders of magnitude, for reasons that remain unclear. An example is the pair of N-(2′-deoxyguanosin-8-yl)-2-aminofluorene (dG-C8-AF) and N-(2′-deoxyguanosin-8-yl)-2-acetylaminofluorene (dG-C8-AAF) adducts that differ by a single acetyl group. The NER efficiencies in human HeLa cell extracts of these lesions are significantly different when placed at G₁, G₂ or G₃ in the duplex sequence (5′-CTCG₁G₂CG₃CCATC-3′) containing the NarI mutational hot spot. Furthermore, the dG-C8-AAF adduct is a better substrate of NER than dG-C8-AF in all three NarI sequence contexts. The conformations of each of these adducts were investigated by Molecular dynamics (MD) simulation methods. In the base-displaced conformational family, the greater repair susceptibility of dG-C8-AAF in all sequences stems from steric hindrance effects of the acetyl group which significantly diminish the adduct-base stabilizing van der Waals stacking interactions relative to the dG-C8-AF case. Base sequence context effects for each adduct are caused by differences in helix untwisting and minor groove opening that are derived from the differences in stacking patterns. Overall, the greater NER efficiencies are correlated with greater extents of base sequence-dependent local untwisting and minor groove opening together with weaker stacking interactions.     

Molecular Dynamics Simulations of Chlorambucil/DNA Adducts. A Structural Basis for the 5′ -GNC Interstrand DNA Crosslink Formed by Nitrogen Mustards

Abstract

The alkylation of DNA by chlorambucil has been studied using a computational approach. Molecular dynamics simulations were performed on the fully solvated non-covalent complex, two monoadducts and a crosslinked diadduct of chlorambucil with the d(CGG₃G₂CGC).- d(GCG₁CCCG) duplex, in which the N7 atoms of G₁, G₂ and G₃ are potential alkylation sites. The results provide a structural basis for the preference of nitrogen mustards to crosslink DNA duplexes at a 5′-GNC site (a 1,3 crosslink, G₁ -G₃) rather than at a 5′-GC sites (a 1,2 crosslink, G₁ -G₂).

In the non-covalent complex simulation the drug reoriented from a non-interstrand crosslinking location to a position favorable for G₁ -G₃ diadduct formation. It proved possible to construct a G₁ -G₃ diadduct from a structure from the non-covalent simulation, and continue the molecular dynamics calculation without further disruption of the DNA structure. A crosslinked diadduct developed with four B_II conformations on the 3′ side of each alkylated guanine and of their respective complementary cytosine. In the first monoadduct simulation the starting point was the same DNA conformation used in the crosslinked diadduct simulation with alkylation at G₁. In this simulation the DNA deformation was reduced, with the helix returning to a more canonical form. A second monoadduct simulation was started from a canonical DNA conformation alkylated at G₃. Here, no significant motion towards a potential crosslinking conformation occurred. Collectively, the results suggest that crosslink formation is dependent upon the drug orientation prior to alkylation and the required deformation of the DNA to permit 1,3 crosslinking can largely be achieved in the non-covalent complex.     

Cooperative kinking at distant sites in mechanically stressed DNA

In cells, DNA is routinely subjected to significant levels of bending and twisting. In some cases, such as under physiological levels of supercoiling, DNA can be so highly strained, that it transitions into non-canonical structural conformations that are capable of relieving mechanical stress within the template. DNA minicircles offer a robust model system to study stress-induced DNA structures. Using DNA minicircles on the order of 100 bp in size, we have been able to control the bending and torsional stresses within a looped DNA construct. Through a combination of cryo-EM image reconstructions, Bal31 sensitivity assays and Brownian dynamics simulations, we have been able to analyze the effects of biologically relevant underwinding-induced kinks in DNA on the overall shape of DNA minicircles. Our results indicate that strongly underwound DNA minicircles, which mimic the physical behavior of small regulatory DNA loops, minimize their free energy by undergoing sequential, cooperative kinking at two sites that are located about 180° apart along the periphery of the minicircle. This novel form of structural cooperativity in DNA demonstrates that bending strain can localize hyperflexible kinks within the DNA template, which in turn reduces the energetic cost to tightly loop DNA.     

Temperature Dependence of the DNA Double Helix at the Nanoscale: Structure,Elasticity, and Fluctuations

Biological organisms exist over a broad temperature range of −15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ ?  + 0.03) in the interval 0–100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.     

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.     

Topoisomer gel retardation: detection of anti-Z-DNA antibodies bound to Z-DNA within supercoiled DNA minicircles.

Small DNA fragments of approximately 350 bp in length, either with or without d(CG)n tracts, are ligated into underwound DNA minicircles to generate topoisomeric rings with different topological linking numbers, Lk. These minicircles, differing by an Lk of one, can be separated by acrylamide gel electrophoresis. Furthermore, electrophoresis can be used to reveal DNA double helix conformational changes that are induced by supercoiling, such as left-handed Z-DNA. When anti-Z-DNA antibodies are added to such minicircles, their binding leads to a selective retardation of the electrophoretic migration of the Z-DNA containing circles. This effect is not seen with relaxed minicircles and those with insufficient torsional stress to induce a conformational transition. Thus the technique of 'topoisomer gel retardation' presents a very sensitive assay for the identification of proteins that selectively bind to DNA conformations stabilized by negative DNA supercoiling.     

Free energy of imperfect nucleic acid helices. 3. Small internal loops resulting from mismatches

Physical studies of enzymioally synthesized oligoribonucleotides of defined sequence are used to evaluate quantitatively the destabilizing influence of mismatched bases in a double helix. The series (A-)₄G(-C)_n(-U)₄, N = 1 to 6, exist as imperfect dimer helices when N is equal to or less than 4, and as monomolecular hairpin helices when N is 5 and 6. Internal loops become progressively more destabilizing as their size increases from 2 to 4 to 6 nucleotides resulting from 1, 2 and 3 consecutive mismatched base pairs. However, the stability of a helix will generally be greater if a given number of mismatched pairs occur consecutively rather than in isolation from one another.These data may be used for improved calculations of stability of RNA secondary structure, to estimate the frequency of structural fluctuations in a double helix and to assess the stability of modified polynucleotide helices. An unmodified double helix of one million randomly arranged base pairs should contain on the time average approximately 10 G.C and 500 A.U pairs in non-hydrogen bonded, unstacked conformations at 25 °C. Our estimate of the effect of mismatching on T_m values of high polymers is less precise because of the long temperature extrapolation required. However, we estimate that DNA or RNA treated with mutagens which interrupt up to 20% of the nucleotide pairs will show a drop of about 1.2 deg. C in melting temperature with each unit per cent of modification.     

Computational analysis of DNA gyrase action

DNA gyrase introduces negative supercoiling into circular DNA by catalyzing the passage of one DNA segment through another. The efficiency of the reaction is many times higher than that of other topological transformations. We analyze, by a computer simulation, the reaction selectivity for a model of DNA gyrase action that assumes existence of a free loop between the G- and T- DNA segments participating in the reaction. A popular model of this type assumed that the selectivity can be provided by the conformation of the DNA segment wrapped around the enzyme into the right-handed helix turn (G-segment). We simulated the distribution of the reaction products for this model. Equilibrium sets of DNA conformations with one segment of the double helix wrapped around the enzyme were constructed. From these sets we selected conformations that had a second segment properly juxtaposed with the first one. Assuming that the juxtapositions result in the strand-passing reaction, we calculated the reaction products for all these conformations. The results show that different products have to be formed if the enzyme acts according to the model. This conclusion can be extended for any model with a free loop between the G- and T-segments. An alternative model that is consistent with the major experimental observations and the computational analysis, is suggested.     

Base Pair Opening Pathways in B-DNA

Abstract

Molecular modeling is used to study the opening pathways of bases within a B-DNA oligomer. It is demonstrated that many open states are possible for a single base pair, although a preference for opening towards the major groove of the double helix is found. In addition we show that opening is strongly influenced by the nature of the base involved and is also coupled in many cases to DNA bending.     

The α-Helix of Ribonuclease T1as an Independent Stability Unit: Direct Comparison of Peptide and Protein Stability

Measurements of the change in conformational stability, Δ(ΔG), upon mutation of two acidic residues at the C terminus of the helix of ribonuclease T₁have recently been reported. Here, we investigate peptides based on the sequence of the helix with the same mutations: Glu28 replaced with Gln, Asp29 replaced with Asn, and the double mutant. In addition, the mutant Lys25 to Gln was studied. Changes in helix content of the peptides with pH confirm the conclusion found in the intact protein, that the charged forms of the acidic residues destabilize the protein by destabilizing the helix. The pH-dependence of the change in confor mational free energy for the peptides and mutant proteins show fair correspondence for D29N and the double mutant. The mutants E28Q and K25Q, on the other hand, give striking agreement between the protein and peptide systems. This agreement suggests that the helix of ribonuclease T₁behaves as an independently stabilized structural unit of the intact protein and that stabilization of the helical form of the peptide is mirrored in the protein.     

Sequence-Dependent Anisotropic Flexibility of B-DNA A Conformational Study

Abstract

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)₂, d(TTAA) ₂, d(GGGG):d(CCCC), d(GGCC) ₂ and d(CCGG) ₂. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5–7°, are in agreement with experimental value of the DNA persistence length.

Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6–12° into this groove. The above inequality is caused by stacking interaction of the bases.

The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20°, the backbone dihedral angles vary by no more than 15°.

The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of “mini-kinks” separated by a half-pitch of the double helix, i.e. by 5–6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimentional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

Conjugation Between G1 and G2 Phase Cells in Paramecium tetraurelia1

Cells of Paramecium tetraurelia, stock hr^d, cultured in a micro-capillary containing 1 μl fresh culture medium, expressed mating activity through the whole cell cycle. Mating-reactive G₂ phase cells can conjugate with cells of other phases. The G₂ phase cells, which have double (4C) the normal micronuclear DNA content, undergo pre-meiotic DNA synthesis when conjugated with G₁ phase cells. The micronucleus of the progeny from the cross between a G₁ and a G₂ cell becomes triploid.     

Temperature Dependence of the DNA Double Helix at the Nanoscale: Structure,Elasticity, and Fluctuations

Biological organisms exist over a broad temperature range of −15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ?+0.03$\text{Δ}\sigma ?+0.03$) in the interval 0–100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.     

Isolation and characterization of kinetoplast DNA from Leishmania tarentolae

Kinetoplast DNA (? = 1.703 g/ml.) was isolated by preparative cesium chloride ultracentrifugation in a fixed-angle rotor from total cell DNA of Leishmania tarentolae and examined in terms of sedimentation properties, melting characteristics, and appearance in the electron microscope. It consisted of several molecular types, either free or bound together in associations of variable size: minicircles (molecular weight = 0.56 ± 0.03 × 10⁶), catenated minicircles, “figure 8” molecules, and long molecules. The associations seem to be held together by the long molecules threading through the smaller circles and catenanes. The large associations could be broken down by sonication, DNase II-treatment, or shear forces. Minicircles, catenated dimers, trimers, and small linear fragments were separated on preparative sucrose gradients of sonicated DNA, and S_20,w values were assigned to each molecular type by band sedimentation in the analytical ultracentrifuge.     

Models of voltage-dependent conformational changes in NaChBac channels

Models of the transmembrane region of the NaChBac channel were developed in two open/inactivated and several closed conformations. Homology models of NaChBac were developed using crystal structures of Kv1.2 and a Kv1.2/2.1 chimera as templates for open conformations, and MlotiK and KcsA channels as templates for closed conformations. Multiple molecular-dynamic simulations were performed to refine and evaluate these models. A striking difference between the S4 structures of the Kv1.2-like open models and MlotiK-like closed models is the secondary structure. In the open model, the first part of S4 forms an α-helix, and the last part forms a 3₁₀ helix, whereas in the closed model, the first part of S4 forms a 3₁₀ helix, and the last part forms an α-helix. A conformational change that involves this type of transition in secondary structure should be voltage-dependent. However, this transition alone is not sufficient to account for the large gating charge movement reported for NaChBac channels and for experimental results in other voltage-gated channels. To increase the magnitude of the motion of S4, we developed another model of an open/inactivated conformation, in which S4 is displaced farther outward, and a number of closed models in which S4 is displaced farther inward. A helical screw motion for the α-helical part of S4 and a simple axial translation for the 3₁₀ portion were used to develop models of these additional conformations. In our models, four positively charged residues of S4 moved outwardly during activation, across a transition barrier formed by highly conserved hydrophobic residues on S1, S2, and S3. The S4 movement was coupled to an opening of the activation gate formed by S6 through interactions with the segment linking S4 to S5. Consistencies of our models with experimental studies of NaChBac and K_v channels are discussed.     

Variations of intramolecular ligation rates allow the detection of protein-induced bends in DNA.

A method requiring minute amounts of DNA and protein is proposed for the detection of DNA bending in solution. Local bending, induced by the binding of a protein at its stereospecific site, must increase the probability of circularization of short DNA fragments and hence the rate of formation of the corresponding minicircles. This effect should be highly sensitive to the location of the protein-binding site on the DNA fragment. Simple controls allow other possible interpretations to be ruled out. The deformation due to the cyclic AMP receptor protein when it interacts with its stereospecific site at the lactose control region is studied by this method. It is shown in this case that untwisting is negligible and bending significant. Further measurements are required to assess the actual value of the radius of curvature of the DNA double helix in this complex.     

Effects of the nucleoid protein HU on the structure, flexibility, and ring-closure properties of DNA deduced from Monte Carlo simulations

The histone-like HU (heat unstable) protein plays a key role in the organization and regulation of the Escherichia coli genome. The nonspecific nature of HU binding to DNA complicates analysis of the mechanism by which the protein contributes to the looping of DNA. Conventional models of the looping of HU-bound duplexes attribute the changes in biophysical properties of DNA brought about by the random binding of protein to changes in the effective parameters of an ideal helical wormlike chain. Here, we introduce a novel Monte Carlo approach to study the effects of nonspecific HU binding on the configurational properties of DNA directly. We randomly decorated segments of an ideal double-helical DNA with HU molecules that induce the bends and other structural distortions of the double helix find in currently available X-ray structures. We find that the presence of HU at levels approximating those found in the cell reduces the persistence length by roughly threefold compared with that of naked DNA. The binding of protein has particularly striking effects on the cyclization properties of short duplexes, altering the dependence of ring closure on chain length in a way that cannot be mimicked by a simple wormlike model and accumulating at higher-than-expected levels on successfully closed chains. Moreover, the uptake of protein on small minicircles depends on chain length, taking advantage of the HU-induced deformations of DNA structure to facilitate ligation. Circular duplexes with bound HU show much greater propensity than protein-free DNA to exist as negatively supercoiled topoisomers, suggesting a potential role of HU in organizing the bacterial nucleoid. The local bending and undertwisting of DNA by HU, in combination with the number of bound proteins, provide a structural rationale for the condensation of DNA and the observed expression levels of reporter genes in vivo.