Chain flexibility and configurational dimensions of left handed Z- DNA and right handed B-DNA helices

The configurational behaviour of flexible helices of right handed B- and left handed Z-types have been analysed using statistical mechanical procedures. The configuration-dependent parameter, most importantly, the persistence length has been computed, using the heminucleotide scheme of treating polynucleotide chains under the approximation that perturbations in the backbone torsions produce sufficient flexibility in these helices. The values of persistence lengths obtained for Z-helices are very much higher than that of B-helices indicating that former is less flexible compared to the latter. These are in accordance with the results obtained recently on B- and Z-forms of poly(dG-dC). (dG-dC) using light scattering studies. Also the persistence lengths of B_II-DNA helices characterised by a skew 3'-hemiucleotide (ε-270°), and also when they coexist with B-DNA have been computed and the values lie within the range of experimentally reported values on B-helices. It is argued that the decrease in the persistence length values of B-DNA at higher salt concentration is due to additional small fluctuations in sugar residue torsions induced due to neutralisation of electrostatic repulsions between adjacent phosphates of the nucleotide. Noteworthy is that these are correlated to winding angle variations and the consequent bending of the helix. Contribution No. 659.     

Energetics of left and right handed models of DNA

It has been shown by model building studies that various right handed and left handed models are compatible with X-ray data of B-DNA and C-DNA. These models are also found to be in good agreement with infrared dichroism data. Detailed potential energy calculations have now been carried out for these models, viz., right and left handed B-DNA and right and left handed C-DNA. It is found that base sugar stacking and interactions involving the phosphate groups are the dominant forces for stabilizing a particular structure. For some sequences, viz., A-A, T-A and C-A, left handed stacking is quite favourable in both B and C structures. But intranucleotide interactions make the left B-DNA unfavourable while the left C-DNA structure is more stable, for all the sequences, than the right C-DNA structure, proposed from fibre data. For the hexanucleoside pentaphosphate fragments the same trend is observed, with the right handed B-DNA being the most stable of the four models studied. However, the left C-DNA structure is only marginally higher in energy, particularly if the shielding effect of the counter ions, on the phosphate group is taken into consideration.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Influence of flanking sequence context on the conformational flexibility of aminofluorene-modified dG adduct in dA mismatch DNA duplexes

When positioned opposite to a dA in a DNA duplex, the prototype arylamine–DNA adduct [N-(2′-deoxyguanosin-yl)-7-fluoro-2-aminofluorene (FAF)] adopts the so-called ‘wedge’ (W) conformation, in which the carcinogen resides in the minor groove of the duplex. All 16 FAF-modified 12-mer NG*N/NAN dA mismatch duplexes (G* = FAF, N = G, A, C, T) exhibited strongly positive induced circular dichroism in the 290–360 nm range (ICD_{290–360 nm}), which supports the W conformation. The ICD_{290–360 nm} intensities were the greatest for duplexes with a 3′-flanking T. The AG*N duplex series showed little adduct-induced destabilization. An exception was the AG*T duplex, which displayed two well-resolved signals in the ¹⁹F NMR spectra. This was presumably due to a strong lesion-destabilizing effect of the 3′-T. The flanking T effect was substantiated further by findings with the TG*T duplex, which exhibited greater lesion flexibility and nucleotide excision repair recognition. Adduct conformational heterogeneity decreased in order of TG*T > AG*T > CG*T > AG*A > AG*G > AG*C. The dramatic flanking T effect on W-conformeric duplexes is consistent with the strong dependence of the ICD_290-360 on both temperature and salt concentration and could be extended to the arylamine food mutagens that are biologically relevant in humans.     

Methods of computer modeling and conformational mobility of DNA duplexes

The review is focused on issues of transferability of the context-sensitive conformational characteristics of DNA estimated from crystallographic structural data on the DNA in aqueous solution. The state of the art in molecular dynamics of charged biopolymers in aqueous solution is covered. Elaboration of expedient force fields and algorithms of calculating long-range electrostatic interactions and solving combined equations of atomic motion have made it possible to generate stable nanosecond trajectories of thermal atomic motion of the biopolymer in aqueous solution in the presence of counterions and salt ions over reasonable time. Tools for analyzing the atomic statistical trajectories of DNA duplexes in aqueous solution to infer context-sensitive conformational dynamic characteristics are discussed together with advances in simulating the mechanisms of global axial bend in DNA duplexes. These techniques allow one to consecutively analyze relationships between the contextual composition of the duplex and the basic modes of essential motions, their amplitude and extent of fluctuation. Development of satisfactory methods for estimating the free energy of biopolymer conformations in solution permits qualitative assessment of the conformational thermodynamic stability of biopolymers and their complexes.     

Rigorous conformational analysis of right and left DNA helices and junction between them

On the basis of complete scanning through conformational space of dihedral angles, twelve structural genera were obtained. Subsequent energy minimization within these genera yielded a limited set of duplexes with stacking: right-handed B-form (Wilkins type), B2-form (Watson and Crick type) and left-handed Ll-form (Sasisekharan type) and the new L2-form. In the polymeric DNA only right-handed double-helices are possible, the left-handed helices are forbidden due to poor 1–5 interchain contacts. In contrast, for short fragments the left- and right-handed helicek have practically the same energies providing some physical ground for side-by-side form, which biologically is possible as recombination form and may be as replication form.     

Probing single-stranded DNA conformational flexibility using fluorescence spectroscopy

Single-stranded DNA (ssDNA) is an essential intermediate in various DNA metabolic processes and interacts with a large number of proteins. Due to its flexibility, the conformations of ssDNA in solution can only be described using statistical approaches, such as flexibly jointed or worm-like chain models. However, there is limited data available to assess such models quantitatively, especially for describing the flexibility of short ssDNA and RNA. To address this issue, we performed FRET studies of a series of oligodeoxythymidylates, (dT)N, over a wide range of salt concentrations and chain lengths (10 < or = N < or = 70 nucleotides), which provide systematic constraints for testing theoretical models. Unlike in mechanical studies where available ssDNA conformations are averaged out during the time it takes to perform measurements, fluorescence lifetimes may act here as an internal clock that influences fluorescence signals depending on how fast the ssDNA conformations fluctuate. A reasonably good agreement could be obtained between our data and the worm-like chain model provided that limited relaxations of the ssDNA conformations occur within the fluorescence lifetime of the donor. The persistence length thus estimated ranges from 1.5 nm in 2 M NaCl to 3 nm in 25 mM NaCl.     

The relative flexibility of B-DNA and A-RNA duplexes: database analysis

An extensive analysis of structural databases is carried out to investigate the relative flexibility of B-DNA and A-RNA duplexes in crystal form. Our results show that the general anisotropic concept of flexibility is not very useful to compare the deformability of B-DNA and A-RNA duplexes, since the flexibility patterns of B-DNA and A-RNA are quite different. In other words, ‘flexibility’ is a dangerous word for describing macromolecules, unless it is clearly defined. A few soft essential movements explain most of the natural flexibility of A-RNA, whereas many are necessary for B-DNA. Essential movements occurring in naked B-DNAs are identical to those necessary to deform DNA in DNA–protein complexes, which suggest that evolution has designed DNA–protein complexes so that B-DNA is deformed according to its natural tendency. DNA is generally more flexible, but for some distortions A-RNA is easier to deform. Local stiffness constants obtained for naked B-DNAs and DNA complexes are very close, demonstrating that global distortions in DNA necessary for binding to proteins are the result of the addition of small concerted deformations at the base-pair level. Finally, it is worth noting that in general the picture of the relative deformability of A-RNA and DNA derived from database analysis agrees very well with that derived from state of the art molecular dynamics (MD) simulations.     

The concentration dependence of the right to left conformational transition in natural DNA identified by Raman spectroscopy

The classical and resonance Raman spectra of DNA from Chicken Erythrocytes have been obtained for different DNA concentrations in solution with low and high ionic strengths. The classical Raman spectra of 30 mg/ml DNA solutions were measured in varying the sodium chloride concentration from 0.1 to 4.5 M NaCl. An increase in the salt content of the solution leads to spectral changes in the 600-700 cm-1 region, indicating a C2' endo/anti to C3' endo/syn conformational transition of the purine residues. Other changes around 840 cm-1, due to the antisymmetrical stretching vibration of the PO2 group, are also detected: they were characteristic for the B----Z transition in model systems such as poly(dG-dC).poly(dG-dC). The resonance Raman spectra of low (1 mg/ml) and high (30 mg/ml) concentrated DNA solutions were obtained with low (0.1 M) and high (4.5 M) NaCl contents, in using a 284 nm excitation wavelength. No change was observed in the intensities and band positions in the low and high salt solutions of low concentrated DNA. Thus it is assumed that the DNA structure remains unchanged whatever the salt concentration for low concentrated DNA. In contrast, great modifications of the intensities and positions of some lines were found in the spectra of high DNA concentration solution when the NaCl content is increased up to 4.5 M: these changes resemble to some extent those observed in the study of B----Z transition of several polynucleotide model compounds. It is assumed that the right-handed to left-handed conformational transition may occur in certain sections of natural DNA, likely containing alternating purine-pyrimidine sequences, when the DNA concentration is sufficiently important.     

Intrinsic conformational flexibility of acetylcholinesterase

Proteins have been metaphorically described - due to the introduction and extraordinary advances in biomolecular dynamics and computational biophysics over the past decades - as "kicking and screaming" molecules [G. Weber, Adv. Protein Chem. 29 (1975) 1-83]. In fact, dynamic fluctuations in protein structural conformation have been known to play an important role in protein function. However, fundamental mechanisms by which protein fluctuations couple with catalytic function of particular enzymes remain poorly understood. To understand the dynamical properties of acetylcholinesterase (AChE) in rapid termination of cationic neurotransmitter, acetylcholine at neurosynaptic junctions, multiple molecular dynamics (MD) trajectories of AChE in the presence and absence of its inhibitors [J.M. Bui, J.A. McCammon, Proc. Natl. Acad. Sci. U.S.A. 103 (2006) 15451-15456; J.M. Bui, Z. Radic, P. Taylor, J.A. McCammon, Biophys. J. 90 (2006) 3280-3287; J.M. Bui, K. Tai, J.A. McCammon, J. Am. Chem. Soc. 126 (2004) 7198-7205; J.M. Bui, R.H. Henchman, J.A. McCammon, Biophys. J. 85 (2003) 2267-2272] have been conducted and correlated with its inhibitory mechanisms. The intrinsic flexibilities of AChE, particularly of the long omega loop, are important in facilitating the ligand's inhibition of the enzyme.     

Pax-3-DNA interaction: flexibility in the DNA binding and induction of DNA conformational changes by paired domains.

The mouse Pax-3 gene encodes a protein that is a member of the Pax family of DNA binding proteins. Pax-3 contains two DNA binding domains: a paired domain (PD) and a paired type homeodomain (HD). Both domains are separated by 53 amino acids and interact synergistically with a sequence harboring an ATTA motif (binding to the HD) and a GTTCC site (binding to the PD) separated by 5 base pairs. Here we show that the interaction of Pax-3 with these two binding sites is independent of their angular orientation. In addition, the protein spacer region between the HD and the PD can be shortened without changing the spatial flexibility of the two DNA binding domains which interact with DNA. Furthermore, by using circular permutation analysis we determined that binding of Pax-3 to a DNA fragment containing a specific binding site causes conformational changes in the DNA, as indicated by the different mobilities of the Pax-3-DNA complexes. The ability to change the conformation of the DNA was found to be an intrinsic property of the Pax-3 PD and of all Pax proteins that we tested so far. These in vitro studies suggest that interaction of Pax proteins with their specific sequences in vivo may result in an altered DNA conformation.     

Long-range effects and conformational flexibility of aldolase

The conformational flexibility and long-range interactions in rabbit muscle aldolase induced by active-site ligand binding, cross-linking of the enzyme between Cys72 and Cys338, and removal of the C-terminal tyrosine residue were studied by following the changes in the microenvironments of Cys239 and Cys289 located outside the active site. It was found that substrates induced a conformational change in aldolase, which propagates from the active site to Cys239, which is located close to intersubunit contacts. The response of the enzyme is differential. Ligands having both C-1 and C-6 phosphates or C-1 phosphate only induce the enhancement of Cys239 reactivity, whereas those with C-6 phosphates only decrease Cys239 reactivity. This correlates well with a dramatic difference in kinetic parameters for a cleavage of fructose-1,6-P2 and fructose-1-P. Therefore, these changes can be interpreted as syncatalytic. Cross-linking of the aldolase subunit by an -S-S-bridge between Cys72 and Cys338 inactivates the enzyme, abolishes binding of active-site ligands, and induces a conformational change in the enzyme that can be detected far away (at Cys239 and Cys289) from the site of perturbation. Cys72 and Cys338 are not in the active site. This shows that the region of the active site and the environment of Cys72 and Cys338 are tightly coupled and that residues far away from the active site, through such coupling, can possess properties of active-site residues. Similar, although less dramatic changes are observed upon removal of the C-terminal tyrosine residue. In view of the results obtained in this paper, aldolase seems to be quite a flexible molecule, whose conformation is sensitive to the nature of a substrate bound to the enzyme and is able to transmit the information about a local perturbation over long distances within a molecule.     

Glycation of aspartate aminotransferase and conformational flexibility

Glycation of proteins alters biological function and changes cellular processes. Our study investigated the conformational changes that accompany glycation using the cardiac aspartate aminotransferase (cAAT). We examined the effects of brief and prolonged exposure of cAAT to glyceraldehyde (Glyc) and ribose 5-phosphate (R5P). When cAAT was briefly incubated (3.5 h) with Glyc (500 microM) or R5P (5 mM) at 37 degrees C, cAAT activity and 1-anilinonaphthalene 8-sulfonate (ANS) binding increased relative to control. After prolonged incubation (64 h) with Glyc (500 microM) or R5P (5 mM) at 37 degrees C, activity and ANS binding decreased relative to control. Furthermore, upon prolonged incubation of cAAT with 500 microM Glyc (14.5 h) or 2 mM R5P (64.25 h) at 37 degrees C, the denaturation curves shifted to the right relative to control. We conclude that upon brief incubation with Glyc and R5P, cAAT exhibited a more open and flexible structure and upon prolonged incubation, a more rigid structure.     

Stability of intramolecular DNA quadruplexes: comparison with DNA duplexes

We have determined the stability of intramolecular quadruplexes that are formed by a variety of G-rich sequences, using oligonucleotides containing appropriately placed fluorophores and quenchers. The stability of these quadruplexes is compared with that of the DNA duplexes that are formed on addition of complementary C-rich oligonucleotides. We find that the linkers joining the G-tracts are not essential for folding and can be replaced with nonnucleosidic moieties, though their sequence composition profoundly affects quadruplex stability. Although the human telomere repeat sequence d[G(3)(TTAG(3))(3)] folds into a quadruplex structure, this forms a duplex in the presence of the complementary C-rich strand at physiological conditions. The Tetrahymena sequence d[G(4)(T(2)G(4))(3)], the sequence d[G(3)(T(2)G(3))(3)], and sequences related to regions of the c-myc promoter d(G(4)AG(4)T)(2) and d(G(4)AG(3)T)(2) preferentially adopt the quadruplex form in potassium-containing buffers, even in the presence of a 50-fold excess of their complementary C-rich strands, though the duplex predominates in the presence of sodium. The HIV integrase inhibitor d[G(3)(TG(3))(3)] forms an extremely stable quadruplex which is not affected by addition of a 50-fold excess of the complementary C-rich strand in both potassium- and sodium-containing buffers. Replacing the TTA loops of the human telomeric repeat with AAA causes a large decrease in quadruplex stability, though a sequence with AAA in the first loop and TTT in the second and third loops is slightly more stable.     

Cryo-EM reveals conformational flexibility in apo DNA polymerase ζ

The translesion synthesis (TLS) DNA polymerases Rev1 and Polζ function together in DNA lesion bypass during DNA replication, acting as nucleotide inserter and extender polymerases, respectively. While the structural characterization of the Saccharomyces cerevisiae Polζ in its DNA-bound state has illuminated how this enzyme synthesizes DNA, a mechanistic understanding of TLS also requires probing conformational changes associated with DNA- and Rev1 binding. Here, we used single-particle cryo-electron microscopy to determine the structure of the apo Polζ holoenzyme. We show that compared with its DNA-bound state, apo Polζ displays enhanced flexibility that correlates with concerted motions associated with expansion of the Polζ DNA-binding channel upon DNA binding. We also identified a lysine residue that obstructs the DNA-binding channel in apo Polζ, suggesting a gating mechanism. The Polζ subunit Rev7 is a hub protein that directly binds Rev1 and is a component of several other protein complexes such as the shieldin DNA double-strand break repair complex. We analyzed the molecular interactions of budding yeast Rev7 in the context of Polζ and those of human Rev7 in the context of shieldin using a crystal structure of Rev7 bound to a fragment of the shieldin-3 protein. Overall, our study provides new insights into Polζ mechanism of action and the manner in which Rev7 recognizes partner proteins.     

MolMovDB: analysis and visualization of conformational change and structural flexibility

The Database of Macromolecular Movements (http://MolMovDB.org) is a collection of data and software pertaining to flexibility in protein and RNA structures. The database is organized into two parts. Firstly, a collection of 'morphs' of solved structures representing different states of a molecule provides quantitative data for flexibility and a number of graphical representations. Secondly, a classification of known motions according to type of conformational change (e.g. 'hinged domain' or 'allosteric') incorporates textual annotation and information from the literature relating to the motion, linking together many of the morphs. A variety of subsets of the morphs are being developed for use in statistical analyses. In particular, for each subset it is possible to derive distributions of various motional quantities (e.g. maximum rotation) that can be used to place a specific motion in context as being typical or atypical for a given population. Over the past year, the database has been greatly expanded and enhanced to incorporate new structures and to improve the quality of data. The 'morph server', which enables users of the database to add new morphs either from their own research or the PDB, has also been enhanced to handle nucleic acid structures and multi-chain complexes.     

Structure and conformational flexibility of Candida rugosa lipase

Anandamide, an endogenous ligand for cannabinoid receptors, loses its biological activities when it is hydrolyzed to arachidonic acid and ethanolamine by anandamide amidohydrolase. We overexpressed a recombinant rat enzyme with a hexahistidine tag in a baculovirus-insect cell expression system, and purified the enzyme with the aid of a Ni-charged resin to a specific activity as high as 5.7 micromol/min/mg protein. The purified recombinant enzyme catalyzed not only the hydrolysis of anandamide and palmitoylethanolamide, but also their reverse synthetic reactions. In order to attain an equilibrium of the anandamide hydrolysis and its reverse reaction within 10 min, we utilized a large amount of the purified enzyme. The equilibrium constant ([arachidonic acid][ethanolamine])/([anandamide][water]) was calculated as 4x10(-3) (37 degrees C, pH 9.0). These experimental results with a purified enzyme preparation quantitatively confirmed the reversibility of the enzyme reaction previously observed with crude enzyme preparations.     

HPLC photofingerprinting of conformational peculiarities and transitions in oligonucleotide duplexes.

Two self-complementary sequence-isomeric decadeoxyribonucleotides were exposed to UV light under conditions in which they assume duplex structures. After that they were analyzed in the denatured state by reversed-phase high-performance liquid chromatography (HPLC). Characterization of the separated photoproducts allowed localization of cyclobutane pyrimidine dimers in the sequences of the modified oligonucleotides. For [d(GGAAATTTCC)]2, which is known to contain in its central part a stretch of rigid B'-conformation with decreased mobility of constituent bases, lower yields of thymine dimers, as compared with that for ordinary B-form [d(CCTTTAAAGG)]2, were found. On the contrary, mixed thymine-cytosine heterodimers generated in the former oligonucleotide demonstrate the increase in photoreactivity of these residues at the B'-B junction. This is probably due to the peculiar conformation adopted by this decanucleotide. Stimulation of B'-B transition, by increasing the temperature before melting, reduced an inhibition of thymine photodimer formation. During the melting of both oligonucleotides yields of all identified photoinduced cyclobutadipyrimidines were reduced. Possible influences of some metal cations on the stability of the B'-form were also studied by this photoprobing technique. The present study demonstrates the feasibility of HPLC photofingerprinting as a new approach for structural analysis of nucleic acids.