Conformational and helicoidal analysis of 30 PS of molecular dynamics on the d(CGCGAATTCGCG) double helix: "curves", dials and windows

Interactions of the DNA Repair Enzyme Human Thymine DNA Glycosylase with Cognate and Noncognate DNA

Glycosylases specifically recognize and flip their target base out of the DNA helix into the enzyme's active site. Our simulations show that a partially flipped state, already present in free DNA carrying a T:G mispair, becomes the more probable state compared to the closed state after binding of thymine DNA glycosylase (TDG). Paired thymine (T:A) or methyl-cytosine (mC:G) does not exhibit a partially flipped state in free or complexed DNA. Important enzyme-DNA interactions exhibit significant strength in the intrahelical and extrahelical TDG-DNA complexes. The computed binding free energy differences suggest these interactions account for the stabilization of the partially flipped state, thereby driving the T:G mispair toward base flip. In the fully flipped state, the cognate base thymine is significantly better accommodated in the enzyme's active site than noncognate bases are, suggesting the hydrolysis step as the last of several stages at which base recognition can be achieved.

Xylonucleic acid: synthesis, structure, and orthogonal pairing properties

There is a common interest for studying xeno-nucleic acid systems in the fields of synthetic biology and the origin of life, in particular, those with an engineered backbone and possessing novel properties. Along this line, we have investigated xylonucleic acid (XyloNA) containing a potentially prebiotic xylose sugar (a 3'-epimer of ribose) in its backbone. Herein, we report for the first time the synthesis of four XyloNA nucleotide building blocks and the assembly of XyloNA oligonucleotides containing all the natural nucleobases. A detailed investigation of pairing and structural properties of XyloNAs in comparison to DNA/RNA has been performed by thermal UV-melting, CD, and solution state NMR spectroscopic studies. XyloNA has been shown to be an orthogonal self-pairing system which adopts a slightly right-handed extended helical geometry. Our study on one hand, provides understanding for superior structure-function (-pairing) properties of DNA/RNA over XyloNA for selection as an informational polymer in the prebiotic context, while on the other hand, finds potential of XyloNA as an orthogonal genetic system for application in synthetic biology.

All-atom polarizable force field for DNA based on the classical Drude oscillator model

Presented is a first generation atomistic force field (FF) for DNA in which electronic polarization is modeled based on the classical Drude oscillator formalism. The DNA model is based on parameters for small molecules representative of nucleic acids, including alkanes, ethers, dimethylphosphate, and the nucleic acid bases and empirical adjustment of key dihedral parameters associated with the phosphodiester backbone, glycosidic linkages, and sugar moiety of DNA. Our optimization strategy is based on achieving a compromise between satisfying the properties of the underlying model compounds in the gas phase targeting quantum mechanical (QM) data and reproducing a number of experimental properties of DNA duplexes in the condensed phase. The resulting Drude FF yields stable DNA duplexes on the 100-ns time scale and satisfactorily reproduce (1) the equilibrium between A and B forms of DNA and (2) transitions between the BI and BII substates of B form DNA. Consistency with the gas phase QM data for the model compounds is significantly better for the Drude model as compared to the CHARMM36 additive FF, which is suggested to be due to the improved response of the model to changes in the environment associated with the explicit inclusion of polarizability. Analysis of dipole moments associated with the nucleic acid bases shows the Drude model to have significantly larger values than those present in CHARMM36, with the dipoles of individual bases undergoing significant variations during the MD simulations. Additionally, the dipole moment of water was observed to be perturbed in the grooves of DNA.

Structural basis for the recognition of diastereomeric 5',8-cyclo-2'-deoxypurine lesions by the human nucleotide excision repair system

The hydroxyl radical is a powerful oxidant that generates DNA lesions including the stereoisomeric R and S 5′,8-cyclo-2′-deoxyadenosine (cdA) and 5′,8-cyclo-2′-deoxyguanosine (cdG) pairs that have been detected in cellular DNA. Unlike some other oxidatively generated DNA lesions, cdG and cdA are repaired by the human nucleotide excision repair (NER) apparatus. The relative NER efficiencies of all four cyclopurines were measured and compared in identical human HeLa cell extracts for the first time under identical conditions, using identical sequence contexts. The cdA and cdG lesions were excised with similar efficiencies, but the efficiencies for both 5′R cyclopurines were greater by a factor of ∼2 than for the 5′S lesions. Molecular modeling and dynamics simulations have revealed structural and energetic origins of this difference in NER-incision efficiencies. These lesions cause greater DNA backbone distortions and dynamics relative to unmodified DNA in 5′R than in 5′S stereoisomers, producing greater impairment in van der Waals stacking interaction energies in the 5′R cases. The locally impaired stacking interaction energies correlate with relative NER incision efficiencies, and explain these results on a structural basis in terms of differences in dynamic perturbations of the DNA backbone imposed by the R and S covalent 5′,8 bonds.

Free energy profiles of base flipping in intercalative polycyclic aromatic hydrocarbon-damaged DNA duplexes: energetic and structural relationships to nucleotide excision repair susceptibility

The crystal structure of Rad4/Rad23, the yeast homolog of the human nucleotide excision repair (NER) lesion recognition factor XPC-RAD23B ( Min , J. H. and Pavletich , N. P. ( 2007 ) Nature 449 , 570 - 575 ) reveals that the lesion-partner base is flipped out of the helix and binds to amino acids of the protein. This suggests the hypothesis that the flipping of this partner base must overcome a free energy barrier, which constitutes one element contributing to changes in the thermodynamic properties induced by the DNA damage and sensed by the recognition protein. We explored this hypothesis by computing complete flipping free energy profiles for two lesions derived from the procarcinogenic polycyclic aromatic hydrocarbons (PAHs), dibenzo[a,l]pyrene (DB[a,l]P) and benzo[a]pyrene (B[a]P), R-trans-anti-DB[a,l]P-N(6)-dA (R-DB[a,l]P-dA) and R-trans-anti-B[a]P-N(6)-dA (R-B[a]P-dA), and the corresponding unmodified duplex. The DB[a,l]P and B[a]P adducts differ in number and organization of their aromatic rings. We integrate these results with prior profiles for the R-trans-anti-DB[a,l]P-dG adduct ( Zheng , H. et al. ( 2010 ) Chem. Res. Toxicol. 23 , 1868 - 1870 ). All adopt conformational themes involving intercalation of the PAH aromatic ring system into the DNA duplex; however, R-DB[a,l]P-dA and R-B[a]P-dA intercalate from the major groove, while R-DB[a,l]P-dG intercalates from the minor groove. These structural differences produce different computed van der Waals stacking interaction energies between the flipping partner base with the lesion aromatic ring system and adjacent bases; we find that the better the stacking, the higher the relative flipping free energy barrier and hence lower flipping probability. The better relative NER susceptibilities correlate with greater ease of flipping in these three differently intercalated lesions. In addition to partner base flipping, the Rad4/Rad23 crystal structure shows that a protein-β-hairpin, BHD3, intrudes from the major groove side between the DNA strands at the lesion site. We present a molecular modeling study for the R-DB[a,l]P-dG lesion in Rad4/Rad23 showing BHD3 β-hairpin intrusion with lesion eviction, and we hypothesize that lesion steric effects play a role in the recognition of intercalated adducts.

Impact of geometry optimization on base-base stacking interaction energies in the canonical A- and B-forms of DNA

Base stacking is known to make an important contribution to the stability of DNA and RNA, and accordingly, significant efforts are ongoing to calculate stacking energies using ab initio quantum mechanical methods. To date, impressive improvements have been made in the model chemistries used to perform stacking energy calculations, including extensions that include robust treatments of electron correlation with extended basis sets, as required to treat interactions where dispersion makes a significant contribution. However, those efforts typically use rigid monomer geometries when calculating the interaction energies. To overcome this, in the present work, we describe a novel internal coordinate definition that allows the relative, intermolecular orientation of stacked base monomers to be constrained during geometry optimizations while allowing full optimization of the intramolecular degrees of freedom. Use of the novel reference frame to calculate the impact of full geometry optimization versus constraining the bases to be planar on base monomer stacking energies, combined with density-fitted, spin-component scaling MP2 treatment of electron correlation, shows that full optimization makes the average stacking energy more favorable by -3.4 and -1.5 kcal/mol for the canonical A and B conformations of the 16 5' to 3' base stacked monomers. Thus, treatment of geometry optimization impacts the stacking energies to an extent similar to or greater than the impact of current state of the art increases in the rigor of the model chemistry itself used to treat base stacking. Results also indicate that stacking favors the B-form of DNA, though the average difference versus the A-form decreases from -2.6 to -0.6 kcal/mol when the intramolecular geometry is allowed to fully relax. However, stacking involving cytosine is shown to favor the A-form of DNA, with that contribution generally larger in the fully optimized bases. The present results show the importance of allowing geometry optimization, as well as properly treating the appropriate model chemistry, in studies of nucleic acid base stacking.

The effect of a G:T mispair on the dynamics of DNA

Distortions in the DNA sequence such as damages or mispairs are specifically recognized and processed by DNA repair enzymes. A particular challenge for the enzymatic specificity is the recognition of a wrongly-placed native nucleotide such as thymine in T:G mispairs. An important step of substrate binding which is observed in many repair proteins is the flipping of the target base out of the DNA helix into the enzyme's active site. In this work we investigate how much the intrinsic dynamics of mispaired DNA is changed compared to canonical DNA. Our molecular dynamics simulations of DNA with and without T:G mispairs show significant differences in the conformation of paired and mispaired DNA. The wobble pair T:G shows local distortions such as twist, shear and stretch which deviate from canonical B form values. Moreover, the T:G mispair is found to be kinetically less stable, exhibiting two states with respect to base opening: a closed state comparable to the canonical base pairs, and a more open state, indicating a proneness for base flip. In addition, we observe that the thymine base in a T:G mispair is significantly more probable to be flipped than thymine in a T:A pair or cytosine in a C:G pair. Such local deformations and in particular the existence of a second, more-open state can be speculated to help the target-site recognition by repair enzymes.

Nuclear magnetic resonance solution structure of an N(2)-guanine DNA adduct derived from the potent tumorigen dibenzo[a,l]pyrene: intercalation from the minor groove with ruptured Watson-Crick base pairing

The most potent tumorigen identified among the polycyclic aromatic hydrocarbons (PAH) is the nonplanar fjord region dibenzo[a,l]pyrene (DB[a,l]P). It is metabolically activated in vivo through the widely studied diol epoxide (DE) pathway to form covalent adducts with DNA bases, predominantly guanine and adenine. The (+)-11S,12R,13R,14S DE enantiomer forms adducts via its C14 position with the exocyclic amino group of guanine. Here, we present the first nuclear magnetic resonance solution structure of a DB[a,l]P-derived adduct, the 14R-(+)-trans-anti-DB[a,l]P-N(2)-dG (DB[a,l]P-dG) lesion in double-stranded DNA. In contrast to the stereochemically identical benzo[a]pyrene-derived N(2)-dG adduct (B[a]P-dG) in which the B[a]P rings reside in the B-DNA minor groove on the 3'-side of the modifed deoxyguanosine, in the DB[a,l]P-derived adduct the DB[a,l]P rings intercalate into the duplex on the 3'-side of the modified base from the sterically crowded minor groove. Watson-Crick base pairing of the modified guanine with the partner cytosine is broken, but these bases retain some stacking with the bulky DB[a,l]P ring system. This new theme in PAH DE-DNA adduct conformation differs from (1) the classical intercalation motif in which Watson-Crick base pairing is intact at the lesion site and (2) the base-displaced intercalation motif in which the damaged base and its partner are extruded from the helix. The structural considerations that lead to the intercalated conformation of the DB[a,l]P-dG lesion in contrast to the minor groove alignment of the B[a]P-dG adduct, and the implications of the DB[a,l]P-dG conformational motif for the recognition of such DNA lesions by the human nucleotide excision repair apparatus, are discussed.

Structure and energy of non-canonical basepairs: comparison of various computational chemistry methods with crystallographic ensembles

Different types of non-canonical basepairs, in addition to the Watson-Crick ones, are observed quite frequently in RNA. Their importance in the three dimensional structure is not fully understood, but their various roles have been proposed by different groups. We have analyzed the energetics and geometry of 32 most frequently observed basepairs in the functional RNA crystal structures using different popular empirical, semi-empirical and ab initio quantum chemical methods and compared their optimized geometry with the crystal data. These basepairs are classified into three categories: polar, non-polar and sugar-mediated, depending on the types of atoms involved in hydrogen bonding. In case of polar basepairs, most of the methods give rise to optimized structures close to their initial geometry. The interaction energies also follow similar trends, with the polar ones having more attractive interaction energies. Some of the C-H...O/N hydrogen bond mediated non-polar basepairs are also found to be significantly stable in terms of their interaction energy values. Few polar basepairs, having amino or carboxyl groups not hydrogen bonded to anything, such as G:G H:W C, show large flexibility. Most of the non-polar basepairs, except A:G s:s T and A:G w:s C, are found to be stable; indicating C-H...O/N interaction also plays a prominent role in stabilizing the basepairs. The sugar mediated basepairs show variability in their structures, due to the involvement of flexible ribose sugar. These presumably indicate that the most of the polar basepairs along with few non-polar ones act as seed for RNA folding while few may act as some conformational switch in the RNA.

Nucleotide excision repair efficiencies of bulky carcinogen-DNA adducts are governed by a balance between stabilizing and destabilizing interactions

The nucleotide excision repair (NER) machinery, the primary defense against cancer-causing bulky DNA lesions, is surprisingly inefficient in recognizing certain mutagenic DNA adducts and other forms of DNA damage. However, the biochemical basis of resistance to repair remains poorly understood. To address this problem, we have investigated a series of intercalated DNA-adenine lesions derived from carcinogenic polycyclic aromatic hydrocarbon (PAH) diol epoxide metabolites that differ in their response to the mammalian NER apparatus. These stereoisomeric PAH-derived adenine lesions represent ideal model systems for elucidating the effects of structural, dynamic, and thermodynamic properties that determine the recognition of these bulky DNA lesions by NER factors. The objective of this work was to gain a systematic understanding of the relation between aromatic ring topology and adduct stereochemistry with existing experimental NER efficiencies and known thermodynamic stabilities of the damaged DNA duplexes. For this purpose, we performed 100 ns molecular dynamics studies of the lesions embedded in identical double-stranded 11-mer sequences. Our studies show that, depending on topology and stereochemistry, stabilizing PAH-DNA base van der Waals stacking interactions can compensate for destabilizing distortions caused by these lesions that can, in turn, cause resistance to NER. The results suggest that the balance between helix stabilizing and destabilizing interactions between the adduct and nearby DNA residues can account for the variability of NER efficiencies observed in this class of PAH-DNA lesions.

Optimization of the CHARMM additive force field for DNA: Improved treatment of the BI/BII conformational equilibrium

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.

DNA sequence context conceals α-anomeric lesions

DNA sequence context has long been known to modulate detection and repair of DNA damage. Recent studies using experimental and computational approaches have sought to provide a basis for this observation. We have previously shown that an α-anomeric adenosine (αA) flanked by cytosines (5'CαAC-3') resulted in a kinked DNA duplex with an enlarged minor groove. Comparison of different flanking sequences revealed that a DNA duplex containing a 5'CαAG-3' motif exhibits unique substrate properties. However, this substrate was not distinguished by unusual thermodynamic properties. To understand the structural basis of the altered recognition, we have determined the solution structure of a DNA duplex with a 5'CαAG-3' core, using an extensive set of restraints including dipolar couplings and backbone torsion angles. The NMR structure exhibits an excellent agreement with the data (total R(X) <5.3%). The αA base is intrahelical, in a reverse Watson-Crick orientation, and forms a weak base pair with a thymine of the opposite strand. In comparison to the DNA duplex with a 5'CαAC-3' core, we observe a significant reduction of the local perturbation (backbone, stacking, tilt, roll, and twist), resulting in a straighter DNA with narrower minor groove. Overall, these features result in a less perturbed DNA helix and obscure the presence of the lesion compared to the 5'CαAC-3' sequence. The improved stacking of the 5'CαAG-3' core also affects the energetics of the DNA deformation that is required to form a catalytically competent complex. These traits provide a rationale for the modulation of the recognition by endonuclease IV.

Solution structure and conformational dynamics of deoxyxylonucleic acids (dXNA): an orthogonal nucleic acid candidate

Orthogonal nucleic acids are chemically modified nucleic acid polymers that are unable to transfer information with natural nucleic acids and thus can be used in synthetic biology to store and transfer genetic information independently. Recently, it was proposed that xylose-DNA (dXNA) can be considered to be a potential candidate for an orthogonal system. Herein, we present the structure in solution and conformational analysis of two self-complementary, fully modified dXNA oligonucleotides, as determined by CD and NMR spectroscopy. These studies are the initial experimental proof of the structural orthogonality of dXNAs. In aqueous solution, dXNA duplexes predominantly form a linear ladderlike (type-1) structure. This is the first example of a furanose nucleic acid that adopts a ladderlike structure. In the presence of salt, an equilibrium exists between two types of duplex form. The corresponding nucleoside triphosphates (dXNTPs) were synthesized and evaluated for their ability to be incorporated into a growing DNA chain by using several natural and mutant DNA polymerases. Despite the structural orthogonality of dXNA, DNA polymerase β mutant is able to incorporate the dXNTPs, showing DNA-dependent dXNA polymerase activity.

Altering the electrostatic potential in the major groove: thermodynamic and structural characterization of 7-deaza-2'-deoxyadenosine:dT base pairing in DNA

As part of an ongoing effort to explore the effect of major groove electrostatics on the thermodynamic stability and structure of DNA, a 7-deaza-2′-deoxyadenosine:dT (7-deaza-dA:dT) base pair in the Dickerson–Drew dodecamer (DDD) was studied. The removal of the electronegative N7 atom on dA and the replacement with an electropositive C–H in the major groove was expected to have a significant effect on major groove electrostatics. The structure of the 7-deaza-dA:dT base pair was determined at 1.1 Å resolution in the presence of Mg²⁺. The 7-deaza-dA, which is isosteric for dA, had minimal effect on the base pairing geometry and the conformation of the DDD in the crystalline state. There was no major groove cation association with the 7-deaza-dA heterocycle. In solution, circular dichroism showed a positive Cotton effect centered at 280 nm and a negative Cotton effect centered at 250 nm that were characteristic of a right-handed helix in the B-conformation. However, temperature-dependent NMR studies showed increased exchange between the thymine N3 imino proton of the 7-deaza-dA:dT base pair and water, suggesting reduced stacking interactions and an increased rate of base pair opening. This correlated with the observed thermodynamic destabilization of the 7-deaza-dA modified duplex relative to the DDD. A combination of UV melting and differential scanning calorimetry experiments were conducted to evaluate the relative contributions of enthalpy and entropy in the thermodynamic destabilization of the DDD. The most significant contribution arose from an unfavorable enthalpy term, which probably results from less favorable stacking interactions in the modified duplex, which was accompanied by a significant reduction in the release of water and cations from the 7-deaza-dA modified DNA.

The miscoding potential of 5-hydroxycytosine arises due to template instability in the replicative polymerase active site

5-Hydroxycytosine (5-OHC) is a stable oxidation product of cytosine associated with an increased frequency of C → T transition mutations. When this lesion escapes recognition by the base excision repair pathway and persists to serve as a templating base during DNA synthesis, replicative DNA polymerases often misincorporate dAMP at the primer terminus, which can lead to fixation of mutations and subsequent disease. To characterize the dynamics of DNA synthesis opposite 5-OHC, we initiated a comparison of unmodified dCMP to 5-OHC, 5-fluorocytosine (5-FC), and 5-methylcytosine (5-MEC) in which these bases act as templates in the active site of RB69 gp43, a high-fidelity DNA polymerase sharing homology with human replicative DNA polymerases. This study presents the first crystal structure of any DNA polymerase binding this physiologically important premutagenic DNA lesion, showing that while dGMP is stabilized by 5-OHC through normal Watson-Crick base pairing, incorporation of dAMP leads to unstacking and instability in the template. Furthermore, the electronegativity of the C5 substituent appears to be important in the miscoding potential of these cytosine-like templates. While dAMP is incorporated opposite 5-OHC ~5 times more efficiently than opposite unmodified dCMP, an elevated level of incorporation is also observed opposite 5-FC but not 5-MEC. Taken together, these data imply that the nonuniform templating by 5-OHC is due to weakened stacking capabilities, which allows dAMP incorporation to proceed in a manner similar to that observed opposite abasic sites.

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.

The role of response elements organization in transcription factor selectivity: the IFN-β enhanceosome example

What is the mechanism through which transcription factors (TFs) assemble specifically along the enhancer DNA? The IFN-β enhanceosome provides a good model system: it is small; its components' crystal structures are available; and there are biochemical and cellular data. In the IFN-β enhanceosome, there are few protein-protein interactions even though consecutive DNA response elements (REs) overlap. Our molecular dynamics (MD) simulations on different motif combinations from the enhanceosome illustrate that cooperativity is achieved via unique organization of the REs: specific binding of one TF can enhance the binding of another TF to a neighboring RE and restrict others, through overlap of REs; the order of the REs can determine which complexes will form; and the alternation of consensus and non-consensus REs can regulate binding specificity by optimizing the interactions among partners. Our observations offer an explanation of how specificity and cooperativity can be attained despite the limited interactions between neighboring TFs on the enhancer DNA. To date, when addressing selective TF binding, attention has largely focused on RE sequences. Yet, the order of the REs on the DNA and the length of the spacers between them can be a key factor in specific combinatorial assembly of the TFs on the enhancer and thus in function. Our results emphasize cooperativity via RE binding sites organization.

Intercalative conformations of the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N⁶-dA adducts: molecular modeling and MD simulations

Among the polycyclic aromatic hydrocarbon class of chemical carcinogens, dibenzo[a,l]pyrene (DB[a,l]P) is the most potent tumorigen that has been identified to date. Structurally, it is bulky with six aromatic rings, and it contains the nonplanar fjord-region. The conformational properties of DB[a,l]P-derived DNA adducts responsible for its extraordinary carcinogenicity are hence of great interest. We have carried out molecular modeling and MD simulations for the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N⁶-dA adducts derived from the reactions of the DB[a,l]P diol epoxides with adenine in double-stranded DNA. The structures are based on the classically intercalated NMR solution structures of the analogous fjord-region benzo[c]phenanthrene-derived-N⁶-dA adducts. One objective was to gain insight on the impact of the more bulky DB[a,l]P ring system on the structural characteristics of the intercalative adduct conformations. A further objective was to elucidate the effect of the flexible twist associated with the sterically hindered aromatic ring in the fjord-region on the intercalated conformations, for comparison with the intercalated but planar bay-region benzo[a]pyrene-derived-N⁶-dA adducts. For the DB[a,l]P-N⁶-dA adducts, our results show that the 14R (+)-adduct is more favorably intercalated on the 5'-side of the modified adenine than the stereoisomeric 14S (-)-adduct, intercalated on its 3'-side. The 14R (+)-adduct manifests better van der Waals stacking interactions with flanking base pairs, less perturbed Watson-Crick hydrogen bonding, less local groove enlargement, less unwinding, and a lower solvent exposure than the 14S (-)-adduct. These structural findings are consistent with observed thermodynamic melting data, UV absorption properties, and fluorescence quenching studies. By contrast, the NMR solution structures for the analogous but less bulky B[c]Ph-derived adducts reveal no such stereoisomeric effect, while the planar bay-region benzo[a]pyrene-derived-N⁶-dA adducts do. Differences in nucleotide excision repair susceptibilities of the fjord and bay region adducts stem from distinctions in their intercalative conformations, produced by the intrinsic topological variations in their polycyclic aromatic ring systems.

Structural and binding study of modified siRNAs with the Argonaute 2 PAZ domain by NMR spectroscopy

By using high-resolution NMR spectroscopy, the structures of a natural short interfering RNA (siRNA) and of several altritol nucleic acid (ANA)-modified siRNAs were determined. The interaction of modified siRNAs with the PAZ domain of the Argonaute 2 protein of Drosophila melanogaster was also studied. The structures show that the modified siRNA duplexes (ANA/RNA) adopt a geometry very similar to the naturally occurring A-type siRNA duplex. All ribose residues, except for the 3' overhang, show 3'-endo conformation. The six-membered altritol sugar in ANA occurs in a chair conformation with the nucleobase in an axial position. In all siRNA duplexes, two overhanging nucleotides at the 3' end enhance the stability of the first neighboring base pair by a stacking interaction. The first overhanging nucleotide has a rather fixed position, whereas the second overhanging nucleotide shows larger flexibility. NMR binding studies of the PAZ domain with ANA-modified siRNAs demonstrate that modifications in the double-stranded region of the antisense strand have some small effects on the binding affinity as compared with the unmodified siRNA. Modification of the 3' overhang with thymidine (dTdT) residues shows a sixfold increase in the binding affinity compared with the unmodified siRNA (relative binding affinity of 17% compared with dTdT-modified overhang), whereas modification of the 3' overhang with ANA largely decreases the binding affinity.

Simulaid: a simulation facilitator and analysis program

Simulaid performs a large number of simulation-related tasks: interconversion and modification of structure and trajectory files, optimization of orientation, and a large variety of analysis functions. The program can handle structures in PDB (Berman et al., Nucleic Acids Res 2000, 28, 235), Charmm (Brooks et al., J Comput Chem 4, 187) CRD, Amber (Case et al.), Macromodel (Mohamadi et al., J Comput Chem 1990, 11, 440), Gromos/Gromacs (Hess et al.), InsightII (InsightII. Accelrys Inc.: San Diego, 2005), Grasp (Nicholls et al., Proteins: Struct Funct Genet 1991, 11, 281) .crg, Tripos (Tripos International, S. H. R., St. Louis, MO) .mol2 (input only), and in the MMC (Mezei, M.; MMC: Monte Carlo program for molecular assemblies. Available at: http://inka.mssm.edu/~mezei/mmc) formats; and trajectories in the formats of Charmm, Amber, Macromodel, and MMC. Analysis features include (but are not limited to): (1) simple distance calculations and hydrogen-bond analysis, (2) calculation of 2-D RMSD maps (produced both as text file with the data and as a color-coded matrix) and cross RMSD maps between trajectories, (3) clustering based on RMSD maps, (4) analysis of torsion angles, Ramachandran (Ramachandran and Sasiskharan, Adv Protein Chem 1968, 23, 283) angles, proline kink (Visiers et al., Protein Eng 2000, 13, 603) angles, pseudorotational (Altona and Sundaralingam, J Am Chem Soc 1972, 94, 8205; Cremer and Pople, J Am Chem Soc 1975, 97, 1354) angles, and (5) analysis based on circular variance (Mezei, J Mol Graphics Model 2003, 21, 463). Torsion angle evolutions are presented in dial plots (Ravishanker et al., J Biomol Struct Dyn 1989, 6, 669). Several of these features are unique to Simulaid.

Lysine120 interactions with p53 response elements can allosterically direct p53 organization

p53 can serve as a paradigm in studies aiming to figure out how allosteric perturbations in transcription factors (TFs) triggered by small changes in DNA response element (RE) sequences, can spell selectivity in co-factor recruitment. p53-REs are 20-base pair (bp) DNA segments specifying diverse functions. They may be located near the transcription start sites or thousands of bps away in the genome. Their number has been estimated to be in the thousands, and they all share a common motif. A key question is then how does the p53 protein recognize a particular p53-RE sequence among all the similar ones? Here, representative p53-REs regulating diverse functions including cell cycle arrest, DNA repair, and apoptosis were simulated in explicit solvent. Among the major interactions between p53 and its REs involving Lys120, Arg280 and Arg248, the bps interacting with Lys120 vary while the interacting partners of other residues are less so. We observe that each p53-RE quarter site sequence has a unique pattern of interactions with p53 Lys120. The allosteric, DNA sequence-induced conformational and dynamic changes of the altered Lys120 interactions are amplified by the perturbation of other p53-DNA interactions. The combined subtle RE sequence-specific allosteric effects propagate in the p53 and in the DNA. The resulting amplified allosteric effects far away are reflected in changes in the overall p53 organization and in the p53 surface topology and residue fluctuations which play key roles in selective co-factor recruitment. As such, these observations suggest how similar p53-RE sequences can spell the preferred co-factor binding, which is the key to the selective gene transactivation and consequently different functional effects.

Distant neighbor base sequence context effects in human nucleotide excision repair of a benzo[a]pyrene-derived DNA lesion

The effects of non-nearest base sequences, beyond the nucleotides flanking a DNA lesion on either side, on nucleotide excision repair (NER) in extracts from human cells were investigated. We constructed two duplexes containing the same minor groove-aligned 10S (+)-trans-anti-B[a]P-N(2)-dG (G*) DNA adduct, derived from the environmental carcinogen benzo[a]pyrene (B[a]P): 5'-C-C-A-T-C-G*-C-T-A-C-C-3' (CG*C-I), and 5'-C-A-C3-A4-C5-G*-C-A-C-A-C-3' (CG*C-II). We used polyacrylamide gel electrophoresis to compare the extent of DNA bending, and molecular dynamics simulations to analyze the structural characteristics of these two DNA duplexes. The NER efficiencies are 1.6(+/-0.2)-fold greater in the case of the CG*C-II than the CG*C-I sequence context in 135-mer duplexes. Gel electrophoresis and self-ligation circularization experiments revealed that the CG*C-II duplex is more bent than the CG*C-I duplex, while molecular dynamics simulations showed that the unique -C3-A4-C5- segment in the CG*C-II duplex plays a key role. The presence of a minor groove-positioned guanine amino group, the Watson-Crick partner to C3, acts as a wedge; facilitated by a highly deformable local -C3-A4- base step, this amino group allows the B[a]P ring system to produce a more enlarged minor groove in CG*C-II than in CG*C-I, as well as a local untwisting and enlarged and flexible Roll only in the CG*C-II sequence. These structural properties fit well with our earlier findings that in the case of the family of minor groove 10S (+)-trans-anti-B[a]P-N(2)-dG lesions, flexible bends and enlarged minor groove widths constitute NER recognition signals, and extend our understanding of sequence context effects on NER to the neighbors that are distant to the lesion.

Insight into F plasmid DNA segregation revealed by structures of SopB and SopB-DNA complexes

Accurate DNA segregation is essential for genome transmission. Segregation of the prototypical F plasmid requires the centromere-binding protein SopB, the NTPase SopA and the sopC centromere. SopB displays an intriguing range of DNA-binding properties essential for partition; it binds sopC to form a partition complex, which recruits SopA, and it also coats DNA to prevent non-specific SopA–DNA interactions, which inhibits SopA polymerization. To understand the myriad functions of SopB, we determined a series of SopB–DNA crystal structures. SopB does not distort its DNA site and our data suggest that SopB–sopC forms an extended rather than wrapped partition complex with the SopA-interacting domains aligned on one face. SopB is a multidomain protein, which like P1 ParB contains an all-helical DNA-binding domain that is flexibly attached to a compact (β₃–α)₂ dimer-domain. Unlike P1 ParB, the SopB dimer-domain does not bind DNA. Moreover, SopB contains a unique secondary dimerization motif that bridges between DNA duplexes. Both specific and non-specific SopB–DNA bridging structures were observed. This DNA-linking function suggests a novel mechanism for in trans DNA spreading by SopB, explaining how it might mask DNA to prevent DNA-mediated inhibition of SopA polymerization.

The energetics of the acetylation switch in p53-mediated transcriptional activation

Targeted therapeutic intervention in receptor-ligand interactions of p53-mediated tumor suppression can impact progression of disease, aging, and variation in genetic expression. Here, we conducted a number of molecular simulations, based on structures of p53 in complex with its transcriptional coactivating CBP bromodomain, determined by NMR spectroscopy, to investigate the energetics of the binding complex. Building on the observation that acetylation of K382 in p53 serves as the essential triggering switch for a specific interaction with CBP, we assessed the differential effect of acetylation on binding from simulations of an octapeptide derived from p53 with acetylated and nonacetylated K382 (residues 379-386). Cluster analysis of the simulations shows that acetylation of the free peptide does not significantly change the population of the preferred conformation of the peptide in solution for binding to CBP. Conversion of the acetylated K382 to nonacetylated form with free energy perturbation (FEP) simulations of the p53 CBP complex and the free peptide showed that the relative contribution of the acetyl group to binding is 4.8 kcal/mol. An analysis of residue contributions to the binding energy using an MM-GBSA approach agrees with the FEP results and sheds additional light on the origin of selectivity in p53 binding to the CBP bromodomain.

Binding of the human nucleotide excision repair proteins XPA and XPC/HR23B to the 5R-thymine glycol lesion and structure of the cis-(5R,6S) thymine glycol epimer in the 5'-GTgG-3' sequence: destabilization of two base pairs at the lesion site

The 5R thymine glycol (5R-Tg) DNA lesion exists as a mixture of cis-(5R,6S) and trans-(5R,6R) epimers; these modulate base excision repair. We examine the 7:3 cis-(5R,6S):trans-(5R,6R) mixture of epimers paired opposite adenine in the 5′-GTgG-3′ sequence with regard to nucleotide excision repair. Human XPA recognizes the lesion comparably to the C8-dG acetylaminoflourene (AAF) adduct, whereas XPC/HR23B recognition of Tg is superior. 5R-Tg is processed by the Escherichia coli UvrA and UvrABC proteins less efficiently than the C8-dG AAF adduct. For the cis-(5R, 6S) epimer Tg and A are inserted into the helix, remaining in the Watson–Crick alignment. The Tg N3H imine and A N⁶ amine protons undergo increased solvent exchange. Stacking between Tg and the 3′-neighbor G•C base pair is disrupted. The solvent accessible surface and T₂ relaxation of Tg increases. Molecular dynamics calculations predict that the axial conformation of the Tg CH₃ group is favored; propeller twisting of the Tg•A pair and hydrogen bonding between Tg OH6 and the N7 atom of the 3′-neighbor guanine alleviate steric clash with the 5′-neighbor base pair. Tg also destabilizes the 5′-neighbor G•C base pair. This may facilitate flipping both base pairs from the helix, enabling XPC/HR23B recognition prior to recruitment of XPA.

NMR and computational studies of stereoisomeric equine estrogen-derived DNA cytidine adducts in oligonucleotide duplexes: opposite orientations of diastereomeric forms

The equine estrogens equilin (EQ) and equilenin (EN) are the active components in the widely prescribed hormone replacement therapy formulation Premarin. Metabolic activation of EQ and EN generates the catechol 4-hydroxyequilenin (4-OHEN) that autoxidizes to the reactive o-quinone form in aerated aqueous solutions. The o-quinones react predominantly with C, and to a lesser extent with A and G, to form premutagenic cyclic covalent DNA adducts in vitro and in vivo. To obtain insights into the structural properties of these biologically important DNA lesions, we have synthesized site-specifically modified oligonucleotides containing the stereoisomeric 1'S,2'R,3'R-4-OHEN-C3 and 1'R,2'S,3'S-4-OHEN-C4 adducts derived from the reaction of 4-OHEN with the C in the oligonucleotide 5'-GGTAGCGATGG in aqueous solution. A combined NMR and computational approach was utilized to determine the conformational characteristics of the two major 4-OHEN-C3 and 4-OHEN-C4 stereoisomeric adducts formed in this oligonucleotide hybridized with its complementary strand. In both cases, the modified C adopts an anti glycosidic bond conformation; the equilenin distal ring protrudes into the minor groove while its two proximal hydroxyl groups are exposed on the major groove side of the DNA duplex. The bulky 4-OHEN-C adduct distorts the duplex within the central GC*G portion, but Watson-Crick pairing is maintained adjacent to C* in both stereoisomeric adducts. For the 4-OHEN-C3 adduct, the equilenin rings are oriented toward the 5'-end of the modified strand, while in 4-OHEN-C4 the equilenin is 3'-directed. Correspondingly, the distortions of the double-helical structures are more pronounced on the 5'- or the 3'-side of the lesion, respectively. These differences in stereoisomeric adduct conformations may play a role in the processing of these lesions in cellular environments.

The Staphylococcus aureus pSK41 plasmid-encoded ArtA protein is a master regulator of plasmid transmission genes and contains a RHH motif used in alternate DNA-binding modes

Plasmids harbored by Staphylococcus aureus are a major contributor to the spread of bacterial multi-drug resistance. Plasmid conjugation and partition are critical to the dissemination and inheritance of such plasmids. Here, we demonstrate that the ArtA protein encoded by the S. aureus multi-resistance plasmid pSK41 is a global transcriptional regulator of pSK41 genes, including those involved in conjugation and segregation. ArtA shows no sequence homology to any structurally characterized DNA-binding protein. To elucidate the mechanism by which it specifically recognizes its DNA site, we obtained the structure of ArtA bound to its cognate operator, ACATGACATG. The structure reveals that ArtA is representative of a new family of ribbon–helix–helix (RHH) DNA-binding proteins that contain extended, N-terminal basic motifs. Strikingly, unlike most well-studied RHH proteins ArtA binds its cognate operators as a dimer. However, we demonstrate that it is also able to recognize an atypical operator site by binding as a dimer-of-dimers and the extended N-terminal regions of ArtA were shown to be essential for this dimer-of-dimer binding mode. Thus, these data indicate that ArtA is a master regulator of genes critical for both horizontal and vertical transmission of pSK41 and that it can recognize DNA utilizing alternate binding modes.

Constructing atomic-resolution RNA structural ensembles using MD and motionally decoupled NMR RDCs

A broad structural landscape often needs to be characterized in order to fully understand how regulatory RNAs perform their biological functions at the atomic level. We present a protocol for visualizing thermally accessible RNA conformations at atomic-resolution and with timescales extending up to milliseconds. The protocol combines molecular dynamics (MD) simulations with experimental residual dipolar couplings (RDCs) measured in partially aligned (13)C/(15)N isotopically enriched elongated RNA samples. The structural ensembles generated in this manner provide insights into RNA dynamics and its role in functionally important transitions.

Differential nucleotide excision repair susceptibility of bulky DNA adducts in different sequence contexts: hierarchies of recognition signals

The structural origin underlying differential nucleotide excision repair (NER) susceptibilities of bulky DNA lesions remains a challenging problem. We investigated the 10S (+)-trans-anti-[BP]-N(2)-2'-deoxyguanosine (G*) adduct in double-stranded DNA. This adduct arises from the reaction, in vitro and in vivo, of a major genotoxic metabolite of benzo[a]pyrene (BP), (+)-(7R,8S,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene, with the exocyclic amino group of guanine. Removal of this lesion by the NER apparatus in cell-free extracts has been found to depend on the base sequence context in which the lesion is embedded, providing an excellent opportunity for elucidating the properties of the damaged DNA duplexes that favor NER. While the BP ring system is in the B-DNA minor groove, 5' directed along the modified strand, there are orientational distinctions that are sequence dependent and are governed by flanking amino groups [Nucleic Acids Res.35 (2007), 1555-1568]. To elucidate sequence-governed NER susceptibility, we conducted molecular dynamics simulations for the 5'-...CG*GC..., 5'-...CGG*C..., and 5'-...TCG*CT... adduct-containing duplexes. We also investigated the 5'-...CG*IC... and 5'-...CIG*C... sequences, which contain "I" (2'-deoxyinosine), with hydrogen replacing the amino group in 2'-deoxyguanosine, to further characterize the structural and dynamic roles of the flanking amino groups in the damaged duplexes. Our results pinpoint explicit roles for the amino groups in tandem GG sequences on the efficiency of NER and suggest a hierarchy of destabilizing structural features that differentially facilitate NER of the BP lesion in the sequence contexts investigated. Furthermore, combinations of several locally destabilizing features in the hierarchy, consistent with a multipartite model, may provide a relatively strong recognition signal.

Dissociation of minor groove binders from DNA: insights from metadynamics simulations

We have used metadynamics to investigate the mechanism of noncovalent dissociation from DNA by two representatives of alkylating and noncovalent minor groove (MG) binders. The compounds are anthramycin in its anhydrous form (IMI) and distamycin A (DST), which differ in mode of binding, size, flexibility and net charge. This choice enables to evaluate the influence of such factors on the mechanism of dissociation. Dissociation of IMI requires an activation free energy of approximately 12 kcal/mol and occurs via local widening of the MG and loss of contacts between the drug and one DNA strand, along with the insertion of waters in between. The detachment of DST occurs at a larger free energy cost, approximately 16.5 or approximately 18 kcal/mol depending on the binding mode. These values compare well with that of 16.6 kcal/mol extracted from stopped-flow experiments. In contrast to IMI, an intermediate is found in which the ligand is anchored to the DNA through its amidinium tail. From this conformation, binding and unbinding occur almost at the same rate. Comparison between DST and with kinetic models for the dissociation of Hoechst 33258 from DNA uncovers common characteristics across different classes of noncovalent MG ligands.

A bridging water anchors the tethered 5-(3-aminopropyl)-2'-deoxyuridine amine in the DNA major groove proximate to the N+2 C.G base pair: implications for formation of interstrand 5'-GNC-3' cross-links by nitrogen mustards

Site-specific insertion of 5-(3-aminopropyl)-2'-deoxyuridine (Z3dU) and 7-deaza-dG into the Dickerson-Drew dodecamers 5'-d(C (1)G (2)C (3)G (4)A (5)A (6)T (7)T (8)C (9) Z (10)C (11)G (12))-3'.5'-d(C (13)G (14)C (15)G (16)A (17)A (18)T (19)T (20)C (21) Z (22)C (23)G (24))-3' (named DDD (Z10)) and 5'-d(C (1)G (2)C (3)G (4)A (5)A (6)T (7) X (8)C (9) Z (10)C (11)G (12))-3'.5'-d(C (13)G (14)C (15)G (16)A (17)A (18)T (19) X (20)C (21) Z (22)C (23)G (24))-3' (named DDD (2+Z10)) (X = Z3dU; Z = 7-deaza-dG) suggests a mechanism underlying the formation of interstrand N+2 DNA cross-links by nitrogen mustards, e.g., melphalan and mechlorethamine. Analysis of the DDD (2+Z10) duplex reveals that the tethered cations at base pairs A (5).X (20) and X (8).A (17) extend within the major groove in the 3'-direction, toward conserved Mg (2+) binding sites located adjacent to N+2 base pairs C (3).Z (22) and Z (10).C (15). Bridging waters located between the tethered amines and either Z (10) or Z (22) O (6) stabilize the tethered cations and allow interactions with the N + 2 base pairs without DNA bending. Incorporation of 7-deaza-dG into the DDD (2+Z10) duplex weakens but does not eliminate electrostatic interactions between tethered amines and Z (10) O (6) and Z (22) O (6). The results suggest a mechanism by which tethered N7-dG aziridinium ions, the active species involved in formation of interstrand 5'-GNC-3' cross-links by nitrogen mustards, modify the electrostatics of the major groove and position the aziridinium ions proximate to the major groove edge of the N+2 C.G base pair, facilitating interstrand cross-linking.

p53-Induced DNA bending: the interplay between p53-DNA and p53-p53 interactions

Specific p53 binding-induced DNA bending and its underlying responsible forces are crucial for the understanding of selective transcription activation. Diverse p53-response elements exist in the genome; however, it is not known what determines the DNA bending and to what extent. In order to gain knowledge of the forces that govern the DNA bending, molecular dynamics simulations were performed on a series of p53 core domain tetramer-DNA complexes in which each p53 core domain was bound to a DNA quarter site specifically. By varying the sequence of the central 4-base pairs of each half-site, different DNA bending extents were observed. The analysis showed that the dimer-dimer interactions in p53 were similar for the complexes; on the other hand, the specific interactions between the p53 and DNA, including the interactions of Arg280, Lys120, and Arg248 with the DNA, varied more significantly. In particular, the Arg280 interactions were better maintained in the complex with the CATG-containing DNA sequence and were mostly lost in the complex with the CTAG-containing DNA sequence. Structural analysis shows that the base pairings for the CATG sequence were stable throughout the simulation trajectory, whereas those for the CTAG sequence were partially dissociated in part of the trajectory, which affected the stability of the nearby Arg280-Gua base interactions. Thus, DNA bending depends on the balance between the p53 dimer-dimer interactions and p53-DNA interactions, which is in turn related to the DNA sequence and DNA flexibility.

NMR determination of oligonucleotide structure

This unit provides an overview of the use of NMR to determine oligonucleotide structure. It covers basic NMR spectral properties, acquisition of interproton distance restraints and torsion angle restraints, structure refinement, assessment of the quality of the structure obtained. Software programs used in the process are also described.

Conformational properties of equilenin-DNA adducts: stereoisomer and base effects

Equilin and equilenin, components of the hormone replacement therapy drug Premarin, can be metabolized to the catechol 4-hydroxyequilenin (4-OHEN). The quinoids produced by 4-OHEN oxidation react with dC, dA, and dG to form unusual stable cyclic adducts, which have been found in human breast tumor tissue. Four stereoisomeric adducts have been identified for each base. These 12 Premarin-derived adducts provide a unique opportunity for analyzing effects of stereochemistry and base damage on DNA structure and consequently its function. Our computational studies have shown that these adducts, with obstructed Watson-Crick hydrogen-bond edges and near-perpendicular ring systems, have limited conformational flexibility and near-mirror-image conformations in stereoisomer pairs. The dC and dA adducts can adopt major- and minor-groove positions in the double helix, but the dG adducts are positioned only in the major groove. In all cases, opposite orientations of the equilenin rings with respect to the 5' --> 3' direction of the damaged strand are found in stereoisomer pairs derived from the same base, and no Watson-Crick pairing is possible. However, detailed structural properties in DNA duplexes are distinct for each stereoisomer of each damaged base. These differences may underlie observed differential stereoisomer and base-dependent mutagenicities and repair susceptibilities of these adducts.

Segrosome structure revealed by a complex of ParR with centromere DNA

The stable inheritance of genetic material depends on accurate DNA partition. Plasmids serve as tractable model systems to study DNA segregation because they require only a DNA centromere, a centromere-binding protein and a force-generating ATPase. The centromeres of partition (par) systems typically consist of a tandem arrangement of direct repeats. The best-characterized par system contains a centromere-binding protein called ParR and an ATPase called ParM. In the first step of segregation, multiple ParR proteins interact with the centromere repeats to form a large nucleoprotein complex of unknown structure called the segrosome, which binds ParM filaments. pSK41 ParR binds a centromere consisting of multiple 20-base-pair (bp) tandem repeats to mediate both transcription autoregulation and segregation. Here we report the structure of the pSK41 segrosome revealed in the crystal structure of a ParR-DNA complex. In the crystals, the 20-mer tandem repeats stack pseudo-continuously to generate the full-length centromere with the ribbon-helix-helix (RHH) fold of ParR binding successive DNA repeats as dimer-of-dimers. Remarkably, the dimer-of-dimers assemble in a continuous protein super-helical array, wrapping the DNA about its positive convex surface to form a large segrosome with an open, solenoid-shaped structure, suggesting a mechanism for ParM capture and subsequent plasmid segregation.

Sliding of alkylating anticancer drugs along the minor groove of DNA: new insights on sequence selectivity

Currently, little is known about the molecular recognition pathways between DNA-alkylating anticancer drugs and their targets despite their pharmacological relevance. In the framework of classical molecular dynamics simulations, here we use umbrella sampling to map the potential of mean force (PMF) associated with sliding along the DNA minor groove of two of these compounds. These are an indole derivative of duocarmycin (DSI) and the putative reactive form of anthramycin (anhydro-anthramycin, IMI). Twenty-three configurations were considered for each drug/DNA complex, corresponding to a movement along approximately 3 basepairs. The alkylation site turns out to be the most favorable for DSI, while a barrier of approximately 6 kcal/mol separates the reactive configuration of IMI.DNA from the absolute minimum. An analysis of various contributions to the PMF reveals that solvent effects play an important role for the largest and more flexible drug DSI. Instead, the PMF of IMI.DNA overall correlates with changes in the binding enthalpy. Implications of these results on the sequence selectivity of the two drugs are discussed.

Accurate prediction of protonation state as a prerequisite for reliable MM-PB(GB)SA binding free energy calculations of HIV-1 protease inhibitors

Binding free energies were calculated for the inhibitors lopinavir, ritonavir, saquinavir, indinavir, amprenavir, and nelfinavir bound to HIV-1 protease. An MMPB/SA-type analysis was applied to conformational samples from 3 ns explicit solvent molecular dynamics simulations of the enzyme-inhibitor complexes. Binding affinities and the sampled conformations of the inhibitor and enzyme were compared between different HIV-1 protease protonation states to find the most likely protonation state of the enzyme in the complex with each of the inhibitors. The resulting set of protonation states leads to good agreement between calculated and experimental binding affinities. Results from the MMPB/SA analysis are compared with an explicit/implicit hybrid scheme and with MMGB/SA methods. It is found that the inclusion of explicit water molecules may offer a slight advantage in reproducing absolute binding free energies while the use of the Generalized Born approximation significantly affects the accuracy of the calculated binding affinities.

Dynamics of a benzo[a]pyrene-derived guanine DNA lesion in TGT and CGC sequence contexts: enhanced mobility in TGT explains conformational heterogeneity, flexible bending, and greater susceptibility to nucleotide excision repair

The nucleotide excision repair (NER) machinery excises a variety of bulky DNA lesions, but with varying efficiencies. The structural features of the DNA lesions that govern these differences are not well understood. An intriguing model system for studying structure-function relationships in NER is the major adduct derived from the reaction of the highly tumorigenic metabolite of benzo[a]pyrene, (+)-anti-benzo[a]pyrene diol epoxide, with the exocyclic amino group of guanine ((+)-trans-anti-[BP]-N(2)-dG, or G*). The rates of incision of the stereochemically identical lesions catalyzed by the prokaryotic UvrABC system was shown to be greater by a factor of 2.3+/-0.3 in the TG*T than in the CG*C sequence context [Biochemistry 46 (2007) 7006-7015]. Here we employ molecular dynamics simulations to elucidate the origin of the greater excision efficiency in the TG*T case and, more broadly, to delineate structural parameters that enhance NER. Our results show that the BP aromatic ring system is 5'-directed along the modified strand in the B-DNA minor groove in both sequence contexts. However, the TG*T modified duplex is much more dynamically flexible, featuring more perturbed and mobile Watson-Crick hydrogen bonding adjacent to the lesion, a greater impairment in stacking interactions, more dynamic local roll/bending, and more minor groove flexibility. These characteristics explain a number of experimental observations concerning the (+)-trans-anti-[BP]-N(2)-dG adduct in double-stranded DNA with the TG*T sequence context: its conformational heterogeneity in NMR solution studies, its highly flexible bend, and its lower thermal stability. By contrast, the CG*C modified duplex is characterized by a single BP conformation and a rigid bend. While current recognition models of bulky lesions by NER factors have stressed the importance of impaired Watson-Crick pairing/stacking and bending, our results highlight the likelihood of an important role for the local dynamics in the vicinity of the lesion.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

Vibrational dynamics of DNA. III. Molecular dynamics simulations of DNA in water and theoretical calculations of the two-dimensional vibrational spectra

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dG(n):dC(n) and dA(n):dT(n) double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG(5)C(5) and dA(5)T(5) double helices in D(2)O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP6-31G(*) calculations for the nine basis modes. The simulated 2D IR spectra of dG(n):dC(n) double helix in D(2)O solution are directly compared with closely related experimental results. The 2D IR spectra of dG(n):dC(n) and dA(n):dT(n) are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions.

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.

4-hydroxyequilenin-adenine lesions in DNA duplexes: stereochemistry, damage site, and structure

The equine estrogens, equilin and equilenin, are major components of the drug Premarin, the most widely used formula for hormone replacement therapy. The derivative 4-hydroxyequilenin (4-OHEN), a major phase I metabolite of equilin and equilenin, autoxidizes to potent cytotoxic quinoids that can react in vitro and in vivo with cytosine and adenine in DNA. Unique cyclic adducts containing the same bicyclo[3.3.1]nonane-type connection ring are produced. Each base adduct has four stereoisomers. In order to elucidate the structural effects of A versus C modification, we have carried out molecular dynamics simulations of the stereoisomeric 4-OHEN-A adducts in DNA 11-mer duplexes and compared results with an earlier study of the C adducts (Ding, S., Shapiro, R., Geacintov, N.E., and Broyde, S. (2005) Equilenin-Derived DNA Adducts to Cytosine in DNA Duplexes: Structures and Thermodynamics, Biochemistry 44, 14565-14576). Similar stereochemical principles govern the orientations in DNA duplexes of the 4-OHEN-A adducts as for the analogous C adducts, with opposite orientations of the equilenin rings in stereoisomeric pairs of adducts characterized by near-mirror image circular dichroism (CD) spectra. However, the larger purine adducts have unique structural properties in the duplexes that distinguish their characteristics from those of the pyrimidine adducts. Significant differences are observed in terms of hydrogen bonding, stacking, bending, groove dimensions, solvent exposure, and hydrophobic interactions; also, each of the four stereoisomeric 4-OHEN-A adducts exhibit distinct structural features. Each base adduct and stereoisomer distorts the structure of the DNA duplex differently. These characteristics may manifest themselves in terms of differential nucleotide excision repair susceptibilities and mutagenic activities of the 4-OHEN-A and C adducts.

Exocyclic amino groups of flanking guanines govern sequence-dependent adduct conformations and local structural distortions for minor groove-aligned benzo[a]pyrenyl-guanine lesions in a GG mutation hotspot context

The environmental carcinogen benzo[a]pyrene (BP) is metabolized to reactive diol epoxides that bind to cellular DNA by predominantly forming N²-guanine adducts (G*). Mutation hotspots for these adducts are frequently found in 5′- ··· GG ··· dinucleotide sequences, but their origins are poorly understood. Here we used high resolution NMR and molecular dynamics simulations to investigate differences in G* adduct conformations in 5′- ··· CG*GC ··· and 5′- ··· CGG* C··· sequence contexts in otherwise identical 12-mer duplexes. The BP rings are positioned 5′ along the modified strand in the minor groove in both cases. However, subtle orientational differences cause strong distinctions in structural distortions of the DNA duplexes, because the exocyclic amino groups of flanking guanines on both strands compete for space with the BP rings in the minor groove, acting as guideposts for placement of the BP. In the 5′- ··· CGG* C ··· case, the 5′-flanking G · C base pair is severely untwisted, concomitant with a bend deduced from electrophoretic mobility. In the 5′- ··· CG*GC ··· context, there is no untwisting, but there is significant destabilization of the 5′-flanking Watson–Crick base pair. The minor groove width opens near the lesion in both cases, but more for 5′- ··· CGG*C···. Differential sequence-dependent removal rates of this lesion result and may contribute to the mutation hotspot phenomenon.

Structure of Aart, a designed six-finger zinc finger peptide, bound to DNA

Cys2-His2 zinc fingers are one of the most common types of DNA-binding domains. Modifications to zinc-finger binding specificity have recently enabled custom DNA-binding proteins to be designed to a wide array of target sequences. We present here a 1.96 A structure of Aart, a designed six-zinc finger protein, bound to a consensus DNA target site. This is the first structure of a designed protein with six fingers, and was intended to provide insights into the unusual affinity and specificity characteristics of this protein. Most protein-DNA contacts were found to be consistent with expectations, while others were unanticipated or insufficient to explain specificity. Several were unexpectedly mediated by glycerol, water molecules or amino acid-base stacking interactions. These results challenge some conventional concepts of recognition, particularly the finding that triplets containing 5'A, C, or T are typically not specified by direct interaction with the amino acid in position 6 of the recognition helix.

Flexible 5-guanidino-4-nitroimidazole DNA lesions: structures and thermodynamics

5-Guanidino-4-nitroimidazole (NI), derived from guanine oxidation by reactive oxygen and nitrogen species, contains an unusual flexible ring-opened structure, with nitro and guanidino groups which possess multiple hydrogen bonding capabilities. In vitro primer extension experiments with bacterial and mammalian polymerases show that NI incorporates C as well as A and G opposite the lesion, depending on the polymerase. To elucidate structural and thermodynamic properties of the mutagenic NI lesion, we have investigated the structure of the modified base itself and the NI-containing nucleoside with high-level quantum mechanical calculations and have employed molecular modeling and molecular dynamics simulations in solution for the lesion in B-DNA duplexes, with four partner bases opposite the NI. Our results show that NI adopts a planar structure at the damaged base level. However, in the nucleoside and in DNA duplexes, steric hindrance between the guanidino group and its linked sugar causes NI to be nonplanar. The NI lesion can adopt both syn and anti conformations on the DNA duplex level, with the guanidino group positioned in the DNA major and minor grooves, respectively; the specific preference depends on the partner base. On the basis of hydrogen bonding and stacking interactions, groove dimensions, and bending, we find that the least distorted NI-modified duplex contains partner C, consistent with observed incorporation of C opposite NI. However, hydrogen bonding interactions between NI and partner G or A are also found, which would be compatible with the observed mismatches.