Structural and Molecular Biology of the eye Lens Membrane

Abstract

The DNA double helix exhibits local sequence-dependent polymorphism at the level of the single base pair and dinucleotide step. Curvature of the DNA molecule occurs in DNA regions with a specific type of nucleotide sequence periodicities. Negative supercoiling induces in vitro local nucleotide sequence-dependent DNA structures such as cruciforms, left-handed DNA, multistranded structures, etc. Techniques based on chemical probes have been proposed that make it possible to study DNA local structures in cells. Recent results suggest that the local DNA structures observed in vitro exist in the cell, but their occurrence and structural details are dependent on the DNA superhelical density in the cell and can be related to some cellular processes.     

Curved DNA

A priori considerations and the concept of the sequence-dependent local curving of the DNA axis. Experimental evidence: electric dichroism (relaxation time measurements); anomalous electrophoretic mobility and gel-filtration of some restriction fragments of DNA; one-sided binding of the nucleosomal DNA to the mica surface. Theoretical predictions concerning the nucleotide sequences of the curved DNA. Discovery of the dinucleotide periodicity in the chromatin DNA. The sequence periodicity as a tool for mapping of the nucleosomes along the sequences. Preferential binding of the histone octamers to the curved pieces of DNA--sequence analysis predictions and comparison with experiments: Theoretical and experimental estimates of the tilt and roll angles for different combinations of the neighboring base-pairs. Inherent sequence-dependent curvature and apparent persistence length of DNA.     

Simple Elastic Network Models for Exhaustive Analysis of Long Double-Stranded DNA Dynamics with Sequence Geometry Dependence

Simple elastic network models of DNA were developed to reveal the structure-dynamics relationships for several nucleotide sequences. First, we propose a simple all-atom elastic network model of DNA that can explain the profiles of temperature factors for several crystal structures of DNA. Second, we propose a coarse-grained elastic network model of DNA, where each nucleotide is described only by one node. This model could effectively reproduce the detailed dynamics obtained with the all-atom elastic network model according to the sequence-dependent geometry. Through normal-mode analysis for the coarse-grained elastic network model, we exhaustively analyzed the dynamic features of a large number of long DNA sequences, approximately ∼150 bp in length. These analyses revealed positive correlations between the nucleosome-forming abilities and the inter-strand fluctuation strength of double-stranded DNA for several DNA sequences.     

Structural bistability of repetitive DNA elements featuring CA/TG dinucleotide steps and mode of evolution of satellite DNA.

Satellite DNA sequences are known to be important components required for the construction of centromeres and are common to all higher eukaryotes. Nevertheless, their nucleotide sequences vary significantly, even in evolutionarily related species. In order to elucidate how the nucleotide sequences define the conformational character of centromeric satellite DNA, an evolutionary path toward repetitive units has been hypothesized. In that context, the DNA conformation of fish satellite DNA was evaluated in two ways: the organization of subrepeats and sequence characteristics were compared, and the differences in stacking energies between A-helix and B-helix and the sequence-dependent bendability of the helices were evaluated. Our findings suggest that the monomeric units making up currently observed repetitive sequences have evolved through stepwise amplification of shorter, ancestral sequences by increasing the length of the units. In addition, we suggest that potentially key sequences required for DNA amplification comprise highly flexible structures. Thus flexibility of the DNA structure may be a primary prerequisite for DNA amplification.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length—a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100–200 Å coherence length, which is similar to ~150 Å coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800–1200 Å, consistent with its theoretically predicted role in interactions between DNA at close separations.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Dual role of DNA intrinsic curvature and flexibility in determining nucleosome stability

A statistical mechanistic approach to evaluate the sequence-dependent thermodynamic stability of nucleosomes is proposed. The model is based on the calculation of the DNA intrinsic curvature, obtained by integrating the nucleotide step deviations from the canonical B-DNA structure, and on the evaluation of the first order elastic distortion energy to reach the nucleosomal superstructure. Literature data on the free energy of nucleosome formation as obtained by competitive nucleosome reconstitution of a significant pool of different DNA sequences were compared with the theoretical results, and a satisfactorily good correlation was found. A striking result of the comparison is the emergence of two opposite roles of the DNA intrinsic curvature and flexibility in determining nucleosome stability. Finally, the obtained results suggest that the curvature-dependent DNA hydration should play a relevant role in the sequence-dependent nucleosome stability.     

The crystal structures of psoralen cross-linked DNAs: drug-dependent formation of Holliday junctions

The single-crystal structures are presented for two DNA sequences with the thymine bases covalently cross-linked across the complementary strands by 4'-hydroxymethyl-4,5',8-trimethylpsoralen (HMT). The HMT-adduct of d(CCGCTAGCGG) forms a psoralen-induced Holliday junction, showing for the first time the effect of this important class of chemotheraputics on the structure of the recombination intermediate. In contrast, HMT-d(CCGGTACCGG) forms a sequence-dependent junction. In both structures, the DNA duplex is highly distorted at the thymine base linked to the six-member pyrone ring of the drug. The psoralen cross-link defines the intramolecular interactions of the drug-induced junction, while the sequence-dependent structure is nearly identical to the native Holliday junction of d(CCGGTACCGG) alone. The two structures contrast the effects of drug- and sequence-dependent interactions on the structure of a Holliday junction, suggesting a role for psoralen in the mechanism to initiate repair of psoralen-lesions in mammalian DNA.     

DNA shape, genetic codes, and evolution

Although the three-letter genetic code that maps nucleotide sequence to protein sequence is well known, there must exist other codes that are embedded in the human genome. Recent work points to sequence-dependent variation in DNA shape as one mechanism by which regulatory and other information could be encoded in DNA. Recent advances include the discovery of shape-dependent recognition of DNA that depends on minor groove width and electrostatics, the existence of overlapping codes in protein-coding regions of the genome, and evolutionary selection for compensatory changes in nucleotide composition that facilitate nucleosome occupancy. It is becoming clear that DNA shape is important to biological function, and therefore will be subject to evolutionary constraint.     

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments’ propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between ‘omics’ and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles’ local geometry consists of straight segments connected by sharp bends compressing the DNA’s inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles’ inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes’ nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.     

Signatures of Protein-DNA Recognition in Free DNA Binding Sites

One obstacle to achieving complete understanding of the principles underlying sequence-dependent recognition of DNA is the paucity of structural data for DNA recognition sequences in their free (unbound) state. Here, we carried out crystallization screening of 50 DNA duplexes containing cognate protein binding sites and obtained new crystal structures of free DNA binding sites for three distinct modes of DNA recognition: anti-parallel β strands (MetR), helix-turn-helix motif + hinge helices (PurR), and zinc fingers (Zif268). Structural changes between free and protein-bound DNA are manifested differently in each case. The new DNA structures reveal that distinctive sequence-dependent DNA geometry dominates recognition by MetR, protein-induced bending of DNA dictates recognition by PurR, and deformability of DNA along the A-B continuum is important in recognition by Zif268. Together, our findings show that crystal structures of free DNA binding sites provide new information about the nature of protein-DNA interactions and thus lend insights towards a structural code for DNA recognition.     

T4 endonuclease VII resolves cruciform DNA with nick and counter-nick and its activity is directed by local nucleotide sequence.

Endonuclease VII of phage T4 resolves Holliday structures in vitro by nicking pairs of strands across the junction. We report here analyses of this reaction between endonuclease VII and a Holliday structure analogue, made in vitro from synthetic oligonucleotides. The enzyme cleaves the structure in a non-concerted way and nicks each strand independently. Combinations of nicks with counter-nicks in strands across the junction resolve the construct. The specificity of the enzyme for DNA secondary structures was tested with a series of branched molecules made from oligonucleotides with the same nucleotide sequence in one strand. Results show that the number, location and relative cleavage efficiencies depend largely on the local nucleotide sequence, rather than on the branch type. In particular, endonuclease VII cleaves a complete four-armed cruciform as efficiently as a three-armed Y-junction or its derivatives, a semi-Y, a fork with two single-strand overhangs, a single-strand overhang, and a nicked DNA. However, exchange or addition of one or more nucleotides within the cleavage area flanking the structural signal for endonuclease VII strongly affects the cleavage pattern as well as their relative efficiency of usage. Examples with a single-stranded overhang are presented and show in summary that the enzyme has a fivefold preference for pyrimidines rather than purines.     

Point mutation in the polyomavirus enhancer alters local DNA conformation

A point mutation in the enhancer of polyomavirus host range mutant, PyEC F441, permits productive infection of the murine embryonal carcinoma cell line, F9. This mutation at nucleotide position 5258 introduces a local conformational change in naked viral DNA. The effect of all four possible nucleotide sequences at position 5258 on local DNA conformation was analyzed by gel electrophoresis of fragments produced by ligation of synthetic oligonucleotides having these sequences. The results indicated that both the wild-type and the F441 sequences introduced local structural polymorphism that can lead to DNA bending. The wild-type sequence had a greater effect on DNA curvature than the F441 sequence. The two other sequences at nucleotide 5258 did not appear to introduce detectable amounts of DNA curvature.     

DNA sequence-dependent deformability--insights from computer simulations

The article reviews some recent developments in studying DNA sequence-dependent deformability, with emphasis on computer modeling. After a brief outline of available experimental techniques, we proceed to computational methods and focus on atomic-resolution molecular dynamics (MD) simulations. A sequence-dependent local (base-pair step) force field inferred from MD is compared with force fields obtained by other techniques. Various methods for establishing global (flexible-rod) DNA elastic constants are reviewed, including an approach based on atomic resolution MD. The problem of defining the global deformation variables, as well as the question of anisotropy and nonlocal effects, are discussed. As an example, both local and global deformability calculations from atomic-resolution MD of EcoRI dodecamer are presented.     

The influence of the 2-amino group of guanine on DNA conformation. Uranyl and DNase I probing of inosine/diaminopurine substituted DNA.

The conformation of the DNA helix is supposed to be a critical element in site-specific recognition by ligands both large and small. Groove width is one important measure of the conformation which varies with the local nucleotide composition, perhaps because of the presence of a purine 2-amino group on G.C base pairs. We have probed DNA with G-->inosine (I) and/or A-->diaminopurine (DAP) substitutions to see whether the location of the purine 2-amino group can indeed affect the minor groove width. At acid pH, the reactivity towards uranyl nitrate is modulated in substituted DNA quite differently from natural DNA, consistent with a marked narrowing of the minor groove at sites of G-->I substitution and widening at sites of A-->DAP replacement. The latter exerts the dominant effect. The expected changes in conformation are equally evident in the patterns of susceptibility to DNase I cleavage, but not to hydroxyl radical attack. Nuclease cleavage is maximal in normal and substituted DNA at regions of inferred moderate groove width which are generally little affected by the nucleotide substitutions. Consistent with models of sequence-dependent cutting by DNase I we find that the presence of a purine 2-amino group on the base pair three places upstream of the cutting site has a profound influence on the rate of reaction.