The flexibility of the nucleic acids: (II). The calculation of internal energy and applications to mononucleotide repeat DNA

Conformation and Rheological Properties of Calf-Thymus DNA in Solution

Studies of DNA molecule behavior in aqueous solutions performed through different approaches allow assessment of the solute-solvent interactions and examination of the strong influence of conformation on its physicochemical properties, in the presence of different ionic species and ionic concentrations. Firstly, the conformational behavior of calf-thymus DNA molecules in TE buffer solution is presented as a function of temperature. Secondly, their rheological behavior is discussed, as well as the evidence of the critical concentrations, i.e., the overlap and the entanglement concentrations (C* and C_e, respectively) from steady state flow and oscillatory dynamic shear experiments. The determination of the viscosity in the Newtonian plateau obtained from flow curves η (γ˙) allows estimation of the intrinsic viscosity and the specific viscosities at zero shear when C[η] < 40. At end, a generalized master curve is obtained from the variation of the specific viscosity as a function of the overlap parameter C[η]. The variation of the exponent s obtained from the power law η~γ˙^−s for both flow and dynamic results is discussed in terms of Graessley’s analysis. In the semi-dilute regime with entanglements, a dynamic master curve is obtained as a function of DNA concentration (C_DNA > 2.0 mg/mL) and temperature (10 °C < T < 40 °C).

Atomistic Monte Carlo simulation of lipid membranes

Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.

THZ-FREQUENCY SPECTROSCOPIC SENSING OF DNA AND RELATED BIOLOGICAL MATERIALS

The terahertz frequency absorption spectra of DNA molecules reflect low-frequency internal helical vibrations involving rigidly bound subgroups that are connected by the weakest bonds, including the hydrogen bonds of the DNA base pairs, and/or non-bonded interactions. Although numerous difficulties make the direct identification of terahertz phonon modes in biological materials very challenging, recent studies have shown that such measurements are both possible and useful. Spectra of different DNA samples reveal a large number of modes and a reasonable level of sequence-specific uniqueness. This chapter utilizes computational methods for normal mode analysis and theoretical spectroscopy to predict the low-frequency vibrational absorption spectra of short artificial DNA and RNA. Here the experimental technique is described in detail, including the procedure for sample preparation. Careful attention was paid to the possibility of interference or etalon effects in the samples, and phenomena were clearly differentiated from the actual phonon modes. The results from Fourier-transform infrared spectroscopy of DNA macromolecules and related biological materials in the terahertz frequency range are presented. In addition, a strong anisotropy of terahertz characteristics is demonstrated. Detailed tests of the ability of normal mode analysis to reproduce RNA vibrational spectra are also conducted. A direct comparison demonstrates a correlation between calculated and experimentally observed spectra of the RNA polymers, thus confirming that the fundamental physical nature of the observed resonance structure is caused by the internal vibration modes in the macromolecules. Application of artificial neural network analysis for recognition and discrimination between different DNA molecules is discussed.

Mechanism of reactant and product dissociation from the anthrax edema factor: a locally enhanced sampling and steered molecular dynamics study

The anthrax edema factor is a toxin overproducing damaging levels of cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi) from ATP. Here, mechanisms of dissociation of ATP and products (cAMP, PPi) from the active site are studied using locally enhanced sampling (LES) and steered molecular dynamics simulations. Various substrate conformations and ionic binding modes found in crystallographic structures are considered. LES simulations show that PPi and cAMP dissociate through different solvent accessible channels, while ATP dissociation requires significant active site exposure to solvent. The ionic content of the active site directly affects the dissociation of ATP and products. Only one ion dissociates along with ATP in the two-Mg(2+) binding site, suggesting that the other ion binds EF prior to ATP association. Dissociation of reaction products cAMP and PPi is impaired by direct electrostatic interactions between products and Mg(2+) ions. This provides an explanation for the inhibitory effect of high Mg(2+) concentrations on EF enzymatic activity. Breaking of electrostatic interactions is dependent on a competitive binding of water molecules to the ions, and thus on the solvent accessibility of the active site. Consequently, product dissociation seems to be a two-step process. First, ligands are progressively solvated while preserving the most important electrostatic interactions, in a process that is dependent on the flexibility of the active site. Second, breakage of the electrostatic bonds follows, and ligands diffuse into solvent. In agreement with this mechanism, product protonation facilitates dissociation.

Activation of the edema factor of Bacillus anthracis by calmodulin: evidence of an interplay between the EF-calmodulin interaction and calcium binding

Calmodulin (CaM) is a remarkably flexible protein which can bind multiple targets in response to changes in intracellular calcium concentration. It contains four calcium-binding sites, arranged in two globular domains. The calcium affinity of CaM N-terminal domain (N-CaM) is dramatically reduced when the complex with the edema factor (EF) of Bacillus anthracis is formed. Here, an atomic explanation for this reduced affinity is proposed through molecular dynamics simulations and free energy perturbation calculations of the EF-CaM complex starting from different crystallographic models. The simulations show that electrostatic interactions between CaM and EF disfavor the opening of N-CaM domains usually induced by calcium binding. Relative calcium affinities of the N-CaM binding sites are probed by free energy perturbation, and dissociation probabilities are evaluated with locally enhanced sampling simulations. We show that EF impairs calcium binding on N-CaM through a direct conformational restraint on Site 1, by an indirect destabilization of Site 2, and by reducing the cooperativity between the two sites.

Toward a consensus view of duplex RNA flexibility

The structure and flexibility of the RNA duplex has been studied using extended molecular dynamics simulations on four diverse 18-mer oligonucleotides designed to contain many copies of the 10 unique dinucleotide steps in different sequence environments. Simulations were performed using the two most popular force fields for nucleic acids simulations (AMBER and CHARMM) in their latest versions, trying to arrive to a consensus picture of the RNA flexibility. Contrary to what was found for DNA duplex (DNA(2)), no clear convergence is found for the RNA duplex (RNA(2)), but one of the force field seems to agree better with experimental data. MD simulations performed with this force field were used to fully characterize, for the first time to our knowledge, the sequence-dependent elastic properties of RNA duplexes at different levels of resolutions. The flexibility pattern of RNA(2) shows similarities with DNA(2), but also surprising differences, which help us to understand the different biological functions of both molecules. A full mesoscopic model of RNA duplex at different resolution levels is derived to be used for genome-wide description of the flexibility of double-helical fragments of RNA.

Use of allostery to identify inhibitors of calmodulin-induced activation of Bacillus anthracis edema factor

Allostery plays a key role in the regulation of the activity and function of many biomolecules. And although many ligands act through allostery, no systematic use is made of it in drug design strategies. Here we describe a procedure for identifying the regions of a protein that can be used to control its activity through allostery. This procedure is based on the construction of a plausible conformational path, which describes protein transition between known active and inactive conformations. The path is calculated by using a framework approach that steers and markedly improves the conjugate peak refinement method. The evolution of conformations along this path was used to identify a putative allosteric site that could regulate activation of Bacillus anthracis adenylyl cyclase toxin (EF) by calmodulin. Conformations of the allosteric site at different steps along the path from the inactive (free) to the active (bound to calmodulin) forms of EF were used to perform virtual screenings and propose candidate EF inhibitors. Several candidates then proved to inhibit calmodulin-induced activation in an in vitro assay. The most potent compound fully inhibited EF at a concentration of 10 microM. The compounds also inhibited the related adenylyl cyclase toxin from Bordetella pertussis (CyaA). The specific homology between the putative allosteric sites in both toxins supports that these pockets are the actual binding sites of the selected inhibitors.

Solvation oscillations and excited-state dynamics of 2-amino- and 2-hydroxy-7-nitrofluorene and its 2'-deoxyriboside

Push-pull substituted fluorenes are considered for use as dynamic solvation probes in polynucleotides. Their fluorescence band is predicted (by simulations) to show weak spectral oscillations on the subpicosecond time scale depending on the nucleotide sequence. The oscillations reflect the local far-infrared spectrum of the environment around the probe molecule. A connection is provided by the continuum theory of polar solvation which, however, neglects molecular aspects. We examine the latter using acetonitrile solution as a test case. A collective librational solvent mode at 100 cm(-1) is observed with 2-amino-7-nitrofluorene, 2-dimethylamino-7-nitrofluorene, 2-hydroxy-7-nitrofluorene, and its 2'-deoxyriboside. Different strengths of the oscillation indicate that rotational friction of nearby acetonitrile molecules depends on the solute structure or that H bonding is involved in launching the librational coherence. Polar solvation in methanol is used for comparison. With hydroxynitrofluorenes, the observation window is limited by intersystem crossing for which rates are reported. A prominent excited-state absorption band of nitrofluorenes at 430 nm can be used to monitor polar solvation. Structural and electronic relaxation pathways are discussed with the help of quantum chemical calculations.

Design of sequence-specific bifunctional nucleic acid ligands

Homopyrimidine oligodeoxynucleotides have been covalently linked to intercalating agents. These bifunctional nucleic acid ligands bind to the major groove of DNA at homopurine.homopyrimidine sequences, where they form triple helices. The homopyrimidine oligonucleotide binds parallel to the purine strand of the double helix. Two hydrogen bonds are formed between bases of the oligonucleotide and the purines engaged in Watson-Crick base pairs. The intercalating agent inserts its aromatic ring at the triplex-duplex junction, resulting in a strong stabilization of the triple helical structure. Bifunctional oligonucleotide-intercalator conjugates provide new tools for a selective control of gene expression. In addition, irreversible reactions can be targeted to the oligonucleotide recognition sequence. Cleavage reactions can be induced by a copper-phenanthroline chelate or an ellipticine derivative covalently linked to the triple helix-forming oligonucleotide.

Molecular flexibility in ab initio drug docking to DNA: binding-site and binding-mode transitions in all-atom Monte Carlo simulations

The dynamics of biological processes depend on the structure and flexibility of the interacting molecules. In particular, the conformational diversity of DNA allows for large deformations upon binding. Drug–DNA interactions are of high pharmaceutical interest since the mode of action of anticancer, antiviral, antibacterial and other drugs is directly associated with their binding to DNA. A reliable prediction of drug–DNA binding at the atomic level by molecular docking methods provides the basis for the design of new drug compounds. Here, we propose a novel Monte Carlo (MC) algorithm for drug–DNA docking that accounts for the molecular flexibility of both constituents and samples the docking geometry without any prior binding-site selection. The binding of the antimalarial drug methylene blue at the DNA minor groove with a preference of binding to AT-rich over GC-rich base sequences is obtained in MC simulations in accordance with experimental data. In addition, the transition between two drug–DNA-binding modes, intercalation and minor-groove binding, has been achieved in dependence on the DNA base sequence. The reliable ab initio prediction of drug–DNA binding achieved by our new MC docking algorithm is an important step towards a realistic description of the structure and dynamics of molecular recognition in biological systems.

Position- and orientation-specific enhancement of topoisomerase I cleavage complexes by triplex DNA structures

Topoisomerase I (Top1) activities are sensitive to various endogenous base modifications, and anticancer drugs including the natural alkaloid camptothecin. Here, we show that triple helix-forming oligonucleotides (TFOs) can enhance Top1-mediated DNA cleavage by affecting either or both the nicking and the closing activities of Top1 depending on the position and the orientation of the triplex DNA structure relative to the Top1 site. TFO binding 1 bp downstream from the Top1 site enhances cleavage by inhibiting religation and to a lesser extent DNA nicking. In contrast, TFO binding 4 bp downstream from the Top1 site enhances DNA nicking especially when the 3' end of the TFO is proximal to the Top1 site. However, when the orientation of the triplex is inverted, with its 5' terminus 4 bp downstream from the Top1 site, religation is also inhibited. These position- and orientation-dependent effects of triplex structures on the Top1-mediated DNA cleavage and religation are discussed in the context of molecular modeling and effects of TFO on DNA twist and mobility at the duplex/triplex junction.

Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5'-YpR-3' intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.

Mechanism of HIV-1 integrase inhibition by styrylquinoline derivatives in vitro

Styrylquinoline derivatives (SQ) efficiently inhibit the 3'-processing activity of integrase (IN) with IC50 values of between 0.5 and 5 microM. We studied the mechanism of action of these compounds in vitro. First, we used steady-state fluorescence anisotropy to assay the effects of the SQ derivatives on the formation of IN-viral DNA complexes independently of the catalytic process. The IC50 values obtained in activity and DNA-binding tests were similar, suggesting that the inhibition of 3'-processing can be fully explained by the prevention of IN-DNA recognition. SQ compounds act in a competitive manner, with Ki values of between 400 and 900 nM. In contrast, SQs did not inhibit 3'-processing when IN-DNA complexes were preassembled. Computational docking followed or not by molecular dynamics using the catalytic core of HIV-1 IN suggested a competitive inhibition mechanism, which is consistent with our previous data obtained with the corresponding Rous sarcoma virus domain. Second, we used preassembled IN-preprocessed DNA complexes to assay the potency of SQs against the strand transfer reaction, independently of 3'-processing. Inhibition occurred even if the efficiency was decreased by about 5- to 10-fold. Our results suggest that two inhibitor-binding modes exist: the first one prevents the binding of the viral DNA and then the two subsequent reactions (i.e., 3'-processing and strand transfer), whereas the second one prevents the binding of target DNA, thus inhibiting strand transfer. SQ derivatives have a higher affinity for the first site, in contrast to that observed for the diketo acids, which preferentially bind to the second one.

Critical assessment of side-chain conformational space sampling procedures designed for quantifying the effect of side-chain environment

We introduce a family of procedures designed to sample side-chain conformational space at particular locations in protein structures. These procedures (CRSP) use intensive cycles of random assignment of side-chain conformations followed by minimization to determine all the conformations that a group of side-chains can adopt simultaneously. First, we consider a procedure evolving in the dihedral space (dCRSP). Our results suggest that it can accurately map low-energy conformations adopted by clusters of side-chains of a protein. dCRSP is relatively insensitive to various important parameters, and it is sufficiently accurate to capture efficiently the constraint induced by the environment on the conformations a particular side-chain can adopt. Our results show that dCRSP, compared with molecular dynamics (MD), can overcome the problem of the limited set of conformations reached in a reasonable amount of simulations. Next, we introduce procedures (vCRSP) in which valence angles are relaxed, and we assess how efficiently they quantify the conformational entropy of side-chains in the protein native state. For simple peptides, entropies obtained with vCRSP are fully compatible with those obtained with a Monte Carlo procedure. For side-chains in a protein environment, however, vCRSP appears of limited use. Finally, we consider a two-step procedure that combines dCRSP and vCRSP. Our tests suggest that it is able to overcome the limitations of vCRSP. We also note that dCRSP provides a reasonable initial approximation. This family of procedures offers promise in quantifying the contribution of conformational entropy to the energetics of protein structures.

Docking macromolecules with flexible segments

We address a major obstacle to macromolecular docking algorithms by presenting a new method that takes into account the induced conformational adjustment of flexible loops situated at a protein/macromolecule interface. The method, MC2, is based on a multiple copy representation of the loops, coupled with a Monte Carlo conformational search of the relative position of the macromolecules and their side chain conformations. The selection of optimal loop conformations takes place during Monte Carlo cycling by the iterative adjustment of the weight of each copy. We describe here the parameterization of the method and trials on a protein-DNA complex of known 3-D structure, involving the Drosophila prd paired domain protein and its target oligonucleotide Wenqing, X. et al., Cell 1995, 80, 639. We demonstrate that our algorithm can correctly configure and position this protein, despite its relatively complex interactions with both grooves of DNA.

Optimization of triple-helix-directed DNA cleavage by benzoquinoquinoxaline-ethylenediaminetetraacetic acid conjugates

The formation of triple-helical structures of DNA is based on sequence-specific recognition of oligopyrimidine.oligopurine stretches of double-helical DNA. Triple-helical structures can be stabilized by DNA-binding ligands. Benzoquinoquinoxaline (BQQ) derivatives are among the most potent intercalating-type agents known to stabilize DNA triple-helical structures. We previously reported the conversion of BQQ into a triplex-directed DNA cleaving agent, namely BQQ-ethylenediaminetetraacetic acid (EDTA), by coupling of 6-(3-aminopropylamino)BQQ to a suitable ethylenediaminetetraacetic acid derivative, and we demonstrated the ability of this conjugate to cause double-stranded cleavage of DNA at the triplex site. However, this prototype derivative BQQ-EDTA conjugate showed lower affinity towards triplex DNA than BQQ itself. In the light of this observation, and guided by molecular modeling studies, we synthesized a second generation of BQQ-EDTA conjugates based on 6-[bis(2-aminoethyl)amino]- and 6-(3,3'-diamino-N-methyldipropylamino)-BQQ derivatives. We confirmed by DNA melting experiments that the new conjugates displayed an increased specific affinity towards triple helices when compared to the previously synthesized BQQ-EDTA. In addition, the efficiency of these new agents in triplex-specific binding and cleavage was demonstrated by triplex-directed double-stranded cleavage of plasmid DNA.

Molecular dynamics and mechanics calculations on a DNA duplex with A(+)-C, G-T and T-C mispairs

Despite major advances in characterizing purine(R)-purine(R), purine(R)-pyrimidine(Y) and pyrimidine(Y)-pyrimidine(Y) mismatches in DNA, there have not been any structural studies on a synthetic DNA duplex containing several different mispairs. Here, using NMR restrained molecular mechanics and dynamics simulations we have structurally characterized a 12 nucleotide long antiparallel DNA duplex with three different mispairs, namely A+-C, G-T and T-C. Our results show that the overall conformation of the antiparallel DNA duplex is B-DNA-like with slight structural distortions at or near the mispairs' sites. All these mispairs are properly stacked with their flanking base pairs. Each mispair is stabilized by two hydrogen bonds and the decreasing order of the hydrogen-bonding interactions is G-T>T-C>A+-C. G-T mispair has smaller configurational space while the structure is slightly bent at A+-C mispair's site. Overall, this study is the first ever structural characterization of a DNA duplex with three different mismatched base pairs and throws light upon the local conformations of the three mispairs present in the DNA duplex.

Solvation model dependency of helix-coil transition in polyalanine

Helix-coil transitions in polyalanine molecules of length 10 are studied by multi-canonical Monte Carlo simulations. The solvation effects are included by either a distance-dependent dielectric permittivity or by a term that is proportional to the solvent-accessible surface area of the peptide. We found a strong dependence of the characteristics of the helix-coil transition from the details of the solvation model.

Theoretical study of fast repair of DNA damage by cistanoside C and analogs: mechanism and docking

Experiments show that the natural substances phenylpropanoid glycosides (PPGs) extracted from pelicularis spicata are capable of repairing DNA damaged by oxygen radicals. Based on kinetic measurements and experiments on tumor cells, a theoretical study of the interaction between PPG molecules and isolated DNA bases, as well as a DNA fragment has been performed. An interaction mechanism reported early has been refined. The docking calculations performed using junction minimization of nucleic acids (JUMNA) software showed that the PPG molecules can be docked into the minor groove of DNA and form complexes with the geometry suitable for an electron transfer between guanine radical and the ligand. Such complexes can be formed without major distortions of DNA structure and are further stabilized by the interaction with the rhamnosyl side-groups.

Bending and adaptability to proteins of the cAMP DNA-responsive element: molecular dynamics contrasted with NMR

DNA bending is assumed to play a crucial role during recognition of the cAMP-responsive element (CRE) by transcription factors. However, diverging results have been obtained for the bending direction of the unbound double helix. The refined NMR structures present a bend directed toward the minor groove, while biochemical methods conclude that there is a bend toward the major groove. The present 10-ns molecular dynamics (MD) simulation of d(GAGATGACGTCATCTC)(2), which contains the octamer CRE in its center, was carried out with AMBER in explicit water and counterions. It shows that CRE is a flexible segment, although it is bent slightly toward the major groove (7 degrees -8 degrees ) on the average. The MD structure agrees with both the biochemical results and unrefined NMR data. The divergence with the NMR refined structures suggests an improper electrostatic parameterization in the refinement software. The malleability of the central CpG is certainly the major contribution to the curving of the whole CRE segment in both the unbound and bound states. Comparison with the crystal structure of CRE bound to GCN4 shows that the deformation induced by the protein is concentrated mainly on the CpG step, rendering the bound structure of CRE closer to the structure of the 12-0 tetradecanoylphorbol-beta-acetate-responsive element.

Photophysical consequences of coupling bacteriochlorophyll a with serine and its resulting solubility in water

We investigated the dependence on solvents of optical absorption and emission of the bacteriochlorophyll a-serine (BChl-ser), a water soluble bacteriochlorophyll (BChl) derivative. Comparison between the experimental data and those collected for BChl in nonaqueous solvents shows that only a minor interaction takes place between serine and the macrocycle's pi-electron system. Nevertheless, the coupling with serine results in a small enhancement of the nonradiative relaxation rate from the first excited singlet state S1. In buffered aqueous solution (pH = 7.4), the Stokes shift of the BChl-ser fluorescence and its nonradiative relaxation rate are enhanced compared with those in nonaqueous solutions (Scherz, A., S. Katz, Y. Vakrat, V. Brumfeld, E. Gabelmann, D. Leupold, J. R. Norris, H. Scheer and Y. Salomon (1998) Photosynthesis: Mechanisms and Effects, Vol. V (Edited by G. Garab), pp. 4207-4212. Kluwer Academic, Dordrecht.), probably as a result of a hydrogen bonding between the BChl macrocycle and the water molecules. In aprotic solvents, without hydrogen bonds, the permanent dipole moment of the first excited singlet state in both BChl and BChl-ser is increased compared with the ground state by at least 2.5 Debye.

Impact of C5-cytosine methylation on the solution structure of d(GAAAACGTTTTC)2. An NMR and molecular modelling investigation

The solution structures of d(GAAAACGTTTTC)2 and of its methylated derivative d(GAAAAMe5CGTTTTC)2 have been determined by NMR and molecular modelling in order to examine the impact of cytosine methylation on the central CpG conformation. Detailed 1H NMR and 31P NMR investigation of the two oligomers includes quantitative NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Back-calculations of NOESY spectra and simulations of double-quantum-filtered COSY patterns were performed to gain accurate information on interproton distances and sugar phase angles. Molecular models under experimental constraints were generated by energy minimization by means of the molecular mechanics program JUMNA. The MORASS software was used to iteratively refine the structures obtained. After methylation, the oligomer still has a B-DNA conformation. However, there are differences in the structural parameters and the thermal stability as compared to the unmethylated molecule. Careful structural analysis shows that after methylation CpG departs from the usual conformation observed in other ACGT tetramers with different surroundings. Subtle displacements of bases, sugars and backbone imposed by the steric interaction of the two methyl groups inside the major groove are accompanied by severe pinching of the minor groove at the C-G residues.

How NF-kappaB can be attracted by its cognate DNA

NF-kappaB is involved in the transcriptional regulation of a large number of genes, in particular those of human immunodeficiency virus (HIV). Recently, we used NMR spectroscopy and molecular modelling to study the solution structure of a native duplex related to the HIV-1 kappaB site, together with a mutated duplex for which a three base-pair change abolishes NF-kappaB binding. The native duplex shows unusual dynamics of the four steps surrounding the kappaB site. Here, we explore the intrinsic properties of the NMR-refined structures of both duplexes in order to understand why the native sequence is recognised by NF-kappaB among other DNA sequences. We establish that only the native kappaB site can adopt a conformation where its structure (curvature and base displacement), the accessibility and the electrostatic potentials of key atoms become very favourable for binding the large loops of NF-kappaB, in contrast to the mutated duplex. Finally, we show that the neutralisation of phosphate groups contacted by NF-kappaB favours a more canonical DNA structure. These findings lead to a new hypothesis for specific recognition through the phosphodiester backbone dynamics of the sequences flanking a binding site. Such unusual behaviour confers upon the overall duplex properties that can be used by NF-kappaB to select its binding site. Thus, the selectivity determinants for NF-kappaB binding appear to depend on deformability of an "extended" consensus sequence.

A molecular model for RecA-promoted strand exchange via parallel triple-stranded helices

A number of studies have concluded that strand exchange between a RecA-complexed DNA single strand and a homologous DNA duplex occurs via a single-strand invasion of the minor groove of the duplex. Using molecular modeling, we have previously demonstrated the possibility of forming a parallel triple helix in which the single strand interacts with the intact duplex in the minor groove, via novel base interactions (Bertucat et al., J. Biomol. Struct. Dynam. 16:535-546). This triplex is stabilized by the stretching and unwinding imposed by RecA. In the present study, we show that the bases within this triplex are appropriately placed to undergo strand exchange. Strand exchange is found to be exothermic and to result in a triple helix in which the new single strand occupies the major groove. This structure, which can be equated to so-called R-form DNA, can be further stabilized by compression and rewinding. We are consequently able to propose a detailed, atomic-scale model of RecA-promoted strand exchange. This model, which is supported by a variety of experimental data, suggests that the role of RecA is principally to prepare the single strand for its future interactions, to guide a minor groove attack on duplex DNA, and to stabilize the resulting, stretched triplex, which intrinsically favors strand exchange. We also discuss how this mechanism can incorporate homologous recognition.

Stability of G,A triple helices

In this work we selected double-stranded DNA sequences capable of forming stable triplexes at 20 or 50 degrees C with corresponding 13mer purine oligonucleotides. This selection was obtained by a double aptamer approach where both the starting sequences of the oligonucleotides and the target DNA duplex were random. The results of selection were confirmed by a cold exchange method and the influence of the position of a 'mismatch' on the stability of the triplex was documented in several cases. The selected sequences obey two rules: (i) they have a high G content; (ii) for a given G content the stability of the resulting triplex is higher if the G residues lie in stretches. The computer simulation of the Mg2+, Na+and Cl-environment around three triplexes by a density scaled Monte Carlo method provides an interpretation of the experimental observations. The Mg2+cations are statistically close to the G N7 and relatively far from the A N7. The presence of an A repels the Mg2+from adjacent G residues. Therefore, the triplexes are stabilized when the Mg2+can form a continuous spine on G N7.

Conformational variation of the central CG site in d(ATGACGTCAT)2 and d(GAAAACGTTTTC)2. An NMR, molecular modelling and 3D-homology investigation

The determination of the solution structure of two self-complementary oligomers d(ATGACGTCAT)2 (CG10) and d(GAAAACGTTTTC)2 (CG12), both containing the 5'-pur-ACGT-pyr-3' sequence, is reported. The impact of the base context on the conformation of the central CpG site has been examined by a combined approach of: (a) 2D 1H-NMR and 31P-NMR; (b) molecular mechanics under experimental constraints; (c) back-calculations of NOESY spectra and iterative refinements of distances; and (d) 3D-homology search of the central tetrad ACGT within the complete oligonucleotides. A full NMR study of each fragment is achieved by means of standard 2D experiments: NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Sugar phase angle, epsilon-zeta difference angle and NOE-derived distances are input as experimental constraints to generate molecular models by energy minimization with the help of jumna. The morass program is used to iteratively refine the structures obtained. The similarity of the two ACGTs within the whole oligonucleotides is investigated. Both the decamer and the dodecamer adopt a B-like DNA conformation. However, the helical parameters within this conformational type are significantly different in CG12 and CG10. The central CpG step conformation is not locked by its nearest environment (5'A and 3'T) as seen from the structural analysis of ACGT in the two molecules. In CG12, despite the presence of runs of A-T pairs, CpG presents a high twist of 43 degrees and a sugar phase at the guanine of about 180 degrees, previously observed in other ACGT-containing-oligomers. Conversely, ACGT in CG10 exhibits strong inclinations, positive rolls, a flat profile of sugar phase, twist and glycosidic angles, as a result of the nucleotide sequence extending beyond the tetrad. The structural specificity of CG10 and its flexibility (as reflected by its energy) are tentatively related to the process of recognition of the cyclic AMP response element by its cognate protein.

Characterization of flexible molecules in solution: the RGDW peptide

Molecular dynamics simulations with adaptive umbrella sampling of the potential energy are used to study conformations of the adhesion peptide RGDW. The peptide is simulated in a box of explicit water. It results in a combination of room temperature (300 K) simulations, in which conformations dominating the average properties of the system are sampled, with high temperature ( approximately 1000 K) simulations in which free energy barriers separating different local minima are crossed efficiently. The simulations with explicit water are compared to simulations of the isolated peptide using different treatments of the electrostatics, and to published experimental data. There is good agreement for data related to the backbone conformation of the peptide. Some discrepancies are evident for data related to side-chain conformations. Together the simulations and experiments provide a description of the RGDW system that is more detailed and reliable than what can be obtained by either simulations or experiments alone.

A model for parallel triple helix formation by RecA: single-single association with a homologous duplex via the minor groove

The nucleoproteic filaments of RecA polymerized on single stranded DNA are able to integrate double stranded DNA in a coaxial arrangement (with DNA stretched by a factor 1.5), to recognize homologous sequences in the duplex and to perform strand exchange between the single stranded and double stranded molecules. While experimental results favor the hypothesis of an invasion of the minor groove of the duplex by the single strand, parallel minor groove triple helices have never been isolated or even modeled, the minor groove offering little space for a third strand to interact. Based on an internal coordinate modeling study, we show here that such a structure is perfectly conceivable when the two interacting oligomers are stretched by a factor 1.5, in order to open the minor groove of the duplex. The model helix presents characteristics that coincide with known experimental data on unwinding, base pair inclination and inter-proton distances. Moreover, we show that extension and unwinding stabilize the triple helix. New patterns of triplet interaction via the minor groove are presented.

New junction models for alternate-strand triple-helix formation

BACKGROUND

[corrected] Oligonucleotide-directed triple-helix (triplex) formation can interfere with gene expression but only long tracts of oligopyrimidine*oligopurine sequences can be targeted. Attempts have been made to recognize short oligopurine sequences alternating on the two strands of double-stranded DNA by the covalent linkage of two triplex-forming oligonucleotides. Here we focus on the rational optimization of such an alternate-strand triplex formation on a DNA duplex containing a 5'-GpT-3'/3'-CpA-5' or a 5'-TpG-3'/3'-ApC-5' step by combination of (G,T)- and (G,A)-containing oligonucleotides that bind to the oligopurine strands in opposite orientations.

RESULTS

The deletion of one nucleotide in the reverse Hoogsteen region of the oligonucleotide provides the best binding at the 5'GpT-3'/3'-CpA-5' step, whereas the addition of two cytosines as a linker between the two oligonucleotides is the best strategy to cross a 5'-TpG-3'/3'-ApC-5' step. Energy minimization and experimental data suggest that these two cytosines are involved in the formation of two novel base quadruplets.

CONCLUSIONS

These data provide a rational basis for the design of oligonucleotides capable of binding to oligopurine sequences that alternate on the two strands of double-stranded DNA with a 5'-GpT-3'/3'-CpA-5' or a 5'-TpG-3'/3'-ApC-5' step at the junction.

Molecular mechanics calculations on a triple stranded DNA involving C+.G-T and T.A+-C mismatched base triples

We have carried out molecular modeling of a triple stranded pyrimidine(Y). purine(R): pyrimidine(Y) (where ':' refers to Watson-Crick and '.' to Hoogsteen bonding) DNA, formed by a homopurine (d-TGAGGAAAGAAGGT) and homo-pyrimidine (d-CTCCTTTCTTCC). Molecular mechanics calculations using NMR constraints have provided a detailed three dimensional structure of the triplex. The entire stretches of purine and the pyrimidine nucleotides have a conformation close to B-DNA. The three strands are held by the canonical C+.G:C and T.A:T hydrogen bonds. The structure also contains two mismatch C+.G-T and T.A+-C base triples which have been characterized for the first time. In the A+-C base-pair of the T.A+-C triple, both hydrogen donors are situated on the purine (A+(1N) and A+(6N)). We observe a unique hydrogen bonding interaction scheme in case of C+.G-T where one acceptor, G(60), is bonded to three donors (C+(3NH), C+(4NH2) and T(3NH)). Though the C+.G-T base triple is less stable than C+.G:C, it is significantly more stable than T.A:T. On the other hand, T.A+-C is as stable as the T.A:T base triad.

A novel palindromic triple-stranded structure formed by homopyrimidine dodecamer d-CTTCTCCTCTTC and homopurine hexamer d-GAAGAG

We have carried out NMR and molecular mechanics studies on a complex formed when a palindromic homopyrimidine dodecamer (d-CTTCTCCTCTTC) and a homopurine hexamer (d-GAAGAG) are mixed in 1:1 molar ratio in aqueous solutions. Such studies unequivocally establish that two strands of each oligomer combine to form a triple-stranded DNA structure with a palindromic symmetry and with six T.A:T and six C+. G:C hydrogen-bonded base triads. The two purine strands are placed head to head, with their 3' ends facing each other in the center of the structure. One-half of each pyrimidine strand contains protonated and the other half contains non-protonated cytosines. The two half segments containing protonated cytosines are hydrogen bonded to each of the two purine hexamers through Hoogsteen T.A and C+.G base pairing. The segments containing non-protonated cytosines are involved in Watson-Crick (A:T and G:C) base pairing. This leads to a palindromic triplex with a C2-dyad symmetry with respect to the center of the structure. The complex is less stable at neutral pH, but the cytosines involved in Hoogsteen base pairing remain protonated even under these conditions. Molecular mechanics calculations using NMR constraints have provided a detailed three-dimensional structure of the complex. The entire stretches of purine, and the pyrimidine nucleotides have a conformation close to B-DNA.

DNA bending by small, mobile multivalent cations

We propose a purely electrostatic mechanism by which small, mobile, multivalent cations can induce DNA bending. A multivalent cation binds at the entrance to the B-DNA major groove, between the two phosphate strands, electrostatically repelling sodium counterions from the neighboring phosphates. The unscreened phosphates on both strands are strongly attracted to the groove-bound cation. This leads to groove closure, accompanied by DNA bending toward the cationic ligand. We explicitly treat the dynamic character of the cation-DNA interaction using an adiabatic approximation, noting that DNA bending is much slower than the diffusion of nonspecifically bound, mobile cations. We make semiquantitative estimates of the free energy components of bending-electrostatic (with a sigmoidal distance-dependent dielectric function), elastic, and entropic cation localization-and find that the equilibrium state is bent B-DNA stabilized with a self-localized cation. This is a bending polaron, formation of which should be critically dependent on the strength of electrostatic interaction and the concentration of highly mobile cations available for self-localization. We predict that the resultant bend will be large (approximately 20-40 degrees), smooth (because it is spread over 6 bp), and infrequent. The stability of such a bend can be variable, from transient to highly stable (static) bending, observable with standard curvature-measuring techniques. We further predict that this bending mechanism will have an unusual sequence dependence: sequences with less binding specificity will be more bent, unless the specific binding site is in the major groove.

Solution structure of a non-palindromic 16 base-pair DNA related to the HIV-1 kappa B site: evidence for BI-BII equilibrium inducing a global dynamic curvature of the duplex

1H and 31P NMR spectroscopy have been used together with molecular modelling to determine the fine structure of a non-palindromic 16 bp DNA containing the NF-kappa B binding site. Much emphasis has been placed upon NMR optimization of both two-dimensional 31P NMR techniques to extract structural information defining the phosphodiester backbone conformation and selective homonuclear 2D COSY experiments to determine sugar conformations. NMR data show evidence for a dynamic behaviour of steps flanking the ten base-pairs of the NF-kappa B binding site. A BI-BII equilibrium at these steps is demonstrated and two models for each extreme conformation are proposed in agreement with NMR data. In the refined BII structures, the NF-kappa B binding site exhibits an intrinsic curvature towards the major groove that is magnified by the four flanking steps in the BII conformation. Furthermore, the base-pairs are translated into the major groove. Thus, we present a novel mode of dynamic intrinsic curvature compatible with the DNA curvature observed in the X-ray structure of the p50-DNA complex.

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.

Hydration of the dTn.dAn x dTn parallel triple helix: a Fourier transform infrared and gravimetric study correlated with molecular dynamics simulations

We present a comparative analysis of the water organization around the dTn.dAn x dTn triple helix and the Watson-Crick double helix dTn.dAn respectively by means of gravimetric measurements, infrared spectroscopy and molecular dynamics simulations. The hydration per nucleotide determined by gravimetric and spectroscopic methods correlated with the molecular dynamics simulations shows that at high relative humidity (98% RH) the triple helix is less solvated than the duplex (17 +/- 2 water molecules per nucleotide instead of 21 +/-1). The experimental desorption curves are different for both structures and indicate that below 81% RH the triplex becomes more hydrated than the duplex. At this RH the FTIR spectra show the emergence of N-type sugars in the adenosine strand of the triplex. When the third strand is bound in the major groove of the Watson-Crick duplex molecular dynamics simulations show the formation of a spine of water molecules between the two thymidine strands.

Base-base and deoxyribose-base stacking interactions in B-DNA and Z-DNA: a quantum-chemical study

Base-stacking interactions in canonical and crystal B-DNA and in Z-DNA steps are studied using the ab initio quantum-chemical method with inclusion of electron correlation. The stacking energies in canonical B-DNA base-pair steps vary from -9.5 kcal/mol (GG) to -13.2 kcal/mol (GC). The many-body nonadditivity term, although rather small in absolute value, influences the sequence dependence of stacking energy. The base-stacking energies calculated for CGC and a hypothetical TAT sequence in Z-configuration are similar to those in B-DNA. Comparison with older quantum-chemical studies shows that they do not provide even a qualitatively correct description of base stacking. We also evaluate the base-(deoxy)ribose stacking geometry that occurs in Z-DNA and in nucleotides linked by 2',5'-phosphodiester bonds. Although the molecular orbital analysis does not rule out the charge-transfer n-pi* interaction of the sugar 04' with the aromatic base, the base-sugar contact is stabilized by dispersion energy similar to that of stacked bases. The stabilization amounts to almost 4 kcal/mol and is thus comparable to that afforded by normal base-base stacking. This enhancement of the total stacking interaction could contribute to the propensity of short d(CG)n sequences to adopt the Z-conformation.