Lankaš, F*
Dotaz
Zobrazit nápovědu
Oxidatively-generated interstrand cross-links rank among the most deleterious DNA lesions. They originate from abasic sites, whose aldehyde group can form a covalent adduct after condensation with the exocyclic amino group of purines, sometimes with remarkably high yields. We use explicit solvent molecular dynamics simulations to unravel the structures and mechanical properties of two DNA sequences containing an interstrand cross-link. Our simulations palliate the absence of experimental structural and stiffness information for such DNA lesions and provide an unprecedented insight into the DNA embedding of lesions that represent a major challenge for DNA replication, transcription and gene regulation by preventing strand separation. Our results based on quantum chemical calculations also suggest that the embedding of the ICL within the duplex can tune the reaction profile, and hence can be responsible for the high difference in yields of formation.
A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ∼10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simulation of the same DNA sequence, the periodicity of the correlation coefficient is not observed. To study the problem further, we extend an earlier mechanical model of DNA allostery based on constrained minimization of effective quadratic deformation energy of the DNA. We demonstrate that, if the constraints mimicking the bound proteins are properly applied, the periodicity in the binding energy is indeed recovered.
A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.
Reliable representation of the B-DNA base-pair step twist is one of the crucial requirements for theoretical modeling of DNA supercoiling and other biologically relevant phenomena in B-DNA. It has long been suspected that the twist is inaccurately described by current empirical force fields. Unfortunately, comparison of simulation results with experiments is not straightforward because of the presence of BII backbone substates, whose populations may differ in experimental and simulation ensembles. In this work, we provide a comprehensive view of the effect of BII substates on the overall B-DNA helix twist and show how to reliably compare twist values from experiment and simulation in two scenarios. First, for longer DNA segments freely moving in solution, we show that sequence-averaged twists of different BI/BII ensembles can be compared directly because of approximate cancellation of the opposing BII effects. Second, for sequence-specific data, such as a particular base-pair step or tetranucleotide twist, can be compared only for a clearly defined BI/BII backbone conformation. For the purpose of force field testing, we designed a compact set of fourteen 22-base-pair B-DNA duplexes (Set 14) containing all 136 distinct tetranucleotide sequences and carried out a total of 84 μs of molecular dynamics simulations, primarily with the OL15 force field. Our results show that the ff99bsc0εζOL1χOL4, parmbsc1, and OL15 force fields model the B-DNA helical twist in good agreement with X-ray and minicircle ligation experiments. The comprehensive understanding obtained regarding the effect of BII substates on the base-pair step geometry should aid meaningful comparisons of various conformational ensembles in future research.
DNA is the carrier of all cellular genetic information and increasingly used in nanotechnology. Quantitative understanding and optimization of its functions requires precise experimental characterization and accurate modeling of DNA properties. A defining feature of DNA is its helicity. DNA unwinds with increasing temperature, even for temperatures well below the melting temperature. However, accurate quantitation of DNA unwinding under external forces and a microscopic understanding of the corresponding structural changes are currently lacking. Here we combine single-molecule magnetic tweezers measurements with atomistic molecular dynamics and coarse-grained simulations to obtain a comprehensive view of the temperature dependence of DNA twist. Experimentally, we find that DNA twist changes by ΔTw(T) = (-11.0 ± 1.2)°/(°C·kbp), independent of applied force, in the range of forces where torque-induced melting is negligible. Our atomistic simulations predict ΔTw(T) = (-11.1 ± 0.3)°/(°C·kbp), in quantitative agreement with experiments, and suggest that the untwisting of DNA with temperature is predominantly due to changes in DNA structure for defined backbone substates, while the effects of changes in substate populations are minor. Coarse-grained simulations using the oxDNA framework yield a value of ΔTw(T) = (-6.4 ± 0.2)°/(°C·kbp) in semi-quantitative agreement with experiments.
The pyrimidine-pyrimidone (6-4) photoproduct (64-PP) is an important photoinduced DNA lesion constituting a mutational signature for melanoma. The structural impact of 64-PP on DNA complexed with histones affects the lesion mutagenicity and repair but remains poorly understood. Here we investigate the conformational dynamics of DNA-containing 64-PP within the nucleosome core particle by atomic-resolution molecular dynamics simulations and multiscale data analysis. We demonstrate that the histone core exerts important mechanical restraints that largely decrease global DNA structural fluctuations. However, the local DNA flexibility at the damaged site is enhanced due to imperfect structural adaptation to restraints imposed by the histone core. If 64-PP faces the histone core and is therefore not directly accessible by the repair protein, the complementary strand facing the solvent is deformed and exhibits higher flexibility than the corresponding strand in a naked, undamaged DNA. This may serve as an initial recognition signal for repair. Our simulations also pinpoint the structural role of proximal residues from the truncated histone tails.
We present a systematic study of the long-timescale dynamics of the Drew-Dickerson dodecamer (DDD: d(CGCGAATTGCGC)2) a prototypical B-DNA duplex. Using our newly parameterized PARMBSC1 force field, we describe the conformational landscape of DDD in a variety of ionic environments from minimal salt to 2 M Na(+)Cl(-) or K(+)Cl(-) The sensitivity of the simulations to the use of different solvent and ion models is analyzed in detail using multi-microsecond simulations. Finally, an extended (10 μs) simulation is used to characterize slow and infrequent conformational changes in DDD, leading to the identification of previously uncharacterized conformational states of this duplex which can explain biologically relevant conformational transitions. With a total of more than 43 μs of unrestrained molecular dynamics simulation, this study is the most extensive investigation of the dynamics of the most prototypical DNA duplex.
The L1 stalk is a key mobile element of the large ribosomal subunit which interacts with tRNA during translocation. Here, we investigate the structure and mechanical properties of the rRNA H76/H75/H79 three-way junction at the base of the L1 stalk from four different prokaryotic organisms. We propose a coarse-grained elastic model and parameterize it using large-scale atomistic molecular dynamics simulations. Global properties of the junction are well described by a model in which the H76 helix is represented by a straight, isotropically flexible elastic rod, while the junction core is represented by an isotropically flexible spherical hinge. Both the core and the helix contribute substantially to the overall H76 bending fluctuations. The presence of wobble pairs in H76 does not induce any increased flexibility or anisotropy to the helix. The half-closed conformation of the L1 stalk seems to be accessible by thermal fluctuations of the junction itself, without any long-range allosteric effects. Bending fluctuations of H76 with a bulge introduced in it suggest a rationale for the precise position of the bulge in eukaryotes. Our elastic model can be generalized to other RNA junctions found in biological systems or in nanotechnology.
- MeSH
- biomechanika MeSH
- konformace nukleové kyseliny MeSH
- ribozomální proteiny chemie MeSH
- RNA ribozomální 23S chemie MeSH
- simulace molekulární dynamiky MeSH
- velké podjednotky ribozomu archebakteriální chemie MeSH
- velké podjednotky ribozomu bakteriální chemie MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
