Mechanism for verification of mismatched and homoduplex DNAs by nucleotides‐bound MutS analyzed by molecular dynamics simulations

In order to understand how MutS recognizes mismatched DNA and induces the reaction of DNA repair using ATP, the dynamics of the complexes of MutS (bound to the ADP and ATP nucleotides, or not) and DNA (with mismatched and matched base‐pairs) were investigated using molecular dynamics simulations. As for DNA, the structure of the base‐pairs of the homoduplex DNA which interacted with the DNA recognition site of MutS was intermittently disturbed, indicating that the homoduplex DNA was unstable. As for MutS, the disordered loops in the ATPase domains, which are considered to be necessary for the induction of DNA repair, were close to (away from) the nucleotide‐binding sites in the ATPase domains when the nucleotides were (not) bound to MutS. This indicates that the ATPase domains changed their structural stability upon ATP binding using the disordered loop. Conformational analysis by principal component analysis showed that the nucleotide binding changed modes which have structurally solid ATPase domains and the large bending motion of the DNA from higher to lower frequencies. In the MutS–mismatched DNA complex bound to two nucleotides, the bending motion of the DNA at low frequency modes may play a role in triggering the formation of the sliding clamp for the following DNA‐repair reaction step. Moreover, MM‐PBSA/GBSA showed that the MutS–homoduplex DNA complex bound to two nucleotides was unstable because of the unfavorable interactions between MutS and DNA. This would trigger the ATP hydrolysis or separation of MutS and DNA to continue searching for mismatch base‐pairs. Proteins 2016; 84:1287–1303. © 2016 Wiley Periodicals, Inc.     

Related base sequences in the DNA of simple and complex organisms. V. The specificity of interactions between oligonucleotides and denatured DNA

A mixture of differentially labeled mouse and Bacillus subtilis DNA was used as a source of oligodeoxynucleotides of chain lengths from 15 to 40 nucleotides. The extent of interaction of these oligonucleotides with homologous or heterologous DNA bound to membrane filters was measured. The specificity of such interactions increases with chain length and with the incubation temperature. The thermal stability of the complexes is a function of chain length. Homologous oligonucleotide/DNA duplexes of B. subtilis are more stable than those of mouse of corresponding size, consistent with the incidence of partially related base sequences in mouse DNA. Oligonucleotides in this size range are also able to recognize partially complementary base sequences in the DNA of closely related organisms. This approach shows promise as a means of obtaining quantitative estimates of base sequence divergence between the DNAs of related organisms.This research was supported by a grant from the National Science Foundation (GB 6099).     

The solution structure of a d[C(TTCG)G] DNA hairpin and comparison to the unusually stable RNA analogue.

The solution structure of the DNA analogue of the unusually stable r[C(UUCG)G] RNA hairpin, 5'-d[GGA-C(TTCG)GTCC]-3', has been determined by NMR spectroscopy, and its structure has been compared to that of the RNA molecule. The RNA molecule is compact and rigid with a highly structured loop. However, the DNA molecule is much less structured. The DNA hairpin contains a B-form stem of four base pairs. The terminal base pair frays, and the 3'-terminal nucleotides, C11 and C12, are in equilibrium between 2'-endo and 3'-endo conformations. Unlike the RNA loop, the DNA loop contains no syn nucleotides, and there is no evidence for base-base or base-phosphate hydrogen bonding in the loop. The loop is flexible, and reveals no specific internucleotide interactions.     

Optimization of nucleic acid sequences

Base sequence influences the structure, mechanics, dynamics, and interactions of nucleic acids. However, studying all possible sequences for a given fragment leads to a number of base combinations that increases exponentially with length. We present here a novel methodology based on a multi-copy approach enabling us to determine which base sequence favors a given structural change or interaction via a single energy minimization. This methodology, termed ADAPT, has been implemented starting from the JUMNA molecular mechanics program by adding special nucleotides, "lexides," containing all four bases, whose contribution to the energy of the system is weighted by continuously variable coefficients. We illustrate the application of this approach in the case of double-stranded DNA by determining the optimal sequences satisfying structural (B-Z transition), mechanical (intrinsic curvature), and interaction (ligand-binding) properties.     

Mechanics of DNA flexibility visualized by selective 2'-amine acylation at nucleotide bulges

We used selective acylation of 2'-amine-substituted nucleotides to visualize local backbone conformations that occur preferentially at bulged sites in DNA duplexes. 2'-Amine acylation reports local nucleotide flexibility because unconstrained 2'-amino nucleotides more readily reach a reactive conformation in which the amide-forming transition state is stabilized by interactions between the amine nucleophile and the adjacent 3'-phosphodiester group. Bulged 2'-amine-substituted cytidine nucleotides react approximately 20-fold more rapidly than nucleotides constrained by base-pairing at 35 degrees C. In contrast, base-paired 2'-amine-substituted nucleotides flanked by a 5' or 3' bulge react two- or six-fold more rapidly, respectively, than the perfectly paired duplex. The relative lack of 2'-amine reactivity for nucleotides adjacent to a DNA bulge emphasizes, first, that structural perturbations do not extend significantly into the flanking duplex structure. Second, the exquisite sensitivity towards very local perturbations in nucleic acid structure suggests that 2'-amine acylation can be used to chemically interrogate deletion mutations in DNA. Finally, these data support the mechanical interpretation that the reactive ribose conformation for 2'-amine acylation requires that the base lies out of the helix and in the major groove, a mechanistic insight useful for designing 2'-amine-based sensors.     

Helicase from hepatitis C virus,energetics of DNA binding

The ability of a helicase to bind single-stranded nucleic acid is critical for nucleic acid unwinding. The helicase from the hepatitis C virus, NS3 protein, binds to the 3'-DNA or the RNA strand during unwinding. As a step to understand the mechanism of unwinding, DNA binding properties of the helicase domain of NS3 (NS3h) were investigated by fluorimetric binding equilibrium titrations. The global analysis of the binding data by a combinatorial approach was done using MATLAB. NS3h interactions with single-stranded DNA (ssDNA) are 300-1000-fold tighter relative to duplex DNA. The NS3h protein binds to ssDNA less than 15 nt in length with a stoichiometry of one protein per DNA. The minimal ssDNA binding site of NS3h helicase was determined to be 8 nucleotides with the microscopic K(d) of 2-4 nm or an observed free energy of -50 kJ/mol. These NS3h-DNA interactions are highly sensitive to salt, and the K(d) increases 4 times when the NaCl concentration is doubled. Multiple HCV helicase proteins bind to ssDNA >15 nucleotides in length, with an apparent occluded site of 8-11 nucleotides. The DNA binding data indicate that the interactions of multiple NS3h protein molecules with long ssDNA are both noncooperative and sequence-independent. We discuss the DNA binding properties of HCV helicase in relation to other superfamily 1 and 2 helicases. These studies provide the basis to investigate the DNA binding interactions with the unwinding substrate and their modulation by the ATPase activity of HCV helicase.     

Structural factors that determine selectivity of a high fidelity DNA polymerase for deoxy-, dideoxy-, and ribonucleotides

In addition to discriminating against base pair mismatches, DNA polymerases exhibit a high degree of selectivity for deoxyribonucleotides over ribo- or dideoxynucleotides. It has been proposed that a single active site residue (steric gate) blocks productive binding of nucleotides containing 2'-hydroxyls. Although this steric gate plays a role in sugar moiety discrimination, its interactions do not account fully for the observed behavior of mutants. Here we present 10 high resolution crystal structures and enzyme kinetic analyses of Bacillus DNA polymerase I large fragment variants complexed with deoxy-, ribo-, and dideoxynucleotides and a DNA substrate. Taken together, these data present a more nuanced and general mechanism for nucleotide discrimination in which ensembles of intermediate conformations in the active site trap non-cognate substrates. It is known that the active site O-helix transitions from an open state in the absence of nucleotide substrates to a ternary complex closed state in which the reactive groups are aligned for catalysis. Substrate misalignment in the closed state plays a fundamental part in preventing non-cognate nucleotide misincorpation. The structures presented here show that additional O-helix conformations intermediate between the open and closed state extremes create an ensemble of binding sites that trap and misalign non-cognate nucleotides. Water-mediated interactions, absent in the fully closed state, play an important role in formation of these binding sites and can be remodeled to accommodate different non-cognate substrates. This mechanism may extend also to base pair discrimination.     

The role of weak specific and nonspecific interactions in recognition and conversation by enzymes of long DNA

According to a currently accepted model, enzymes engage in high-rate sliding along DNA when searching for specific recognition sequences or structural elements (modified nucleotides, breaks, single-stranded DNA fragments, etc.). Such sliding requires these enzymes to possess sufficiently high affinity for DNA of any sequence. Thus, significant differences in the enzymes' affinity for specific and nonspecific DNA sequences cannot be expected, and formation of a complex between an enzyme and its target DNA unlikely contributes significantly in the enzyme specificity. To elucidate the factors providing the specificity we have analyzed many DNA replication, DNA repair, topoisomerization, integration, and recombination enzymes using a number of physicochemical methods, including a method of stepwise increase in ligand complexity developed in our laboratory. It was shown that high affinity of all studied enzymes for long DNA is provided by formation of many weak contacts of the enzymes with all nucleotide units covered by protein globules. Contacts of positively charged amino acid residues with internucleotide phosphate groups contribute most to such interactions; the contribution of each contact is very small and the full contact interface usually resembles interactions between oppositely charged biopolymer surfaces. In some cases significant contribution to the affinity is made through hydrophobic and/or van der Waals interactions of the enzymes with nucleobases. Overall, depending on the enzyme, such nonspecific interactions provide 5-8 orders of the enzyme affinity for DNA. Specific interactions of enzymes with long DNA, in contrast to contacts of enzymes with small ligands, are usually weak and comparable in efficiency with weak nonspecific contacts. The sum of specific interactions most often provides approximately one and rarely two orders of the affinity. According to structural data, DNA binding to any of the investigated enzymes is followed by a stage of DNA conformation adjustment including partial or complete DNA melting, deformation of its backbone, stretching, compression, bending or kinking, eversion of nucleotides from the DNA helix, etc. The full set of such changes is characteristic for each individual enzyme. The fact that all enzyme-dependent changes in DNA are effected through weak specific rather than strong interactions is very important. Enzyme-specific changes in DNA conformation are required for effective adjustment of reacting orbitals with accuracy about 10-15 degrees, which is possible only for specific DNA. A transition from nonspecific to specific DNA leads to an increase in the reaction rate (kcat) by 4-8 orders of magnitude. Thus, the stages of DNA conformation adjustment and catalysis proper provide the high specificity of enzyme action.     

Structural analysis of prokaryotic and eukaryotic 5S rRNAs by RNase H

The availabilities of single-stranded 5S rRNA regions c, d and d' for base pairing interactions were analyzed by using synthetic DNA oligomers. Hybrid formation was detected by the endonucleolytical mode of the RNA-DNA specific action of RNase H. Provided that the hybrid interaction involved 6 successive base pairs, 5S rRNA loop c nucleotides 42-47 displayed accessibility in Escherichia coli, Bacillus stearothermophilus and Thermus thermophilus 5S rRNAs as well as in eukaryotic 5S rRNAs from Saccharomyces carlsbergensis, Rattus rattus and Equisetum arvense. Investigating eubacterial 5S rRNA regions d and d' (nucleotides 71-76 and 99-105, respectively), susceptibility was observed in E. coli 5S rRNA which, however, decreases in B. stearothermophilus and even more so in T. thermophilus 5S rRNA. For additional evaluation of the data obtained by RNase H cleavage, association constants of the hexanucleotides were determined by equilibrium dialysis at 4 degrees C for B. stearothermophilus 5S rRNA. The results obtained reveal that nucleotides 36-41 of B. stearothermophilus 5S rRNA are inaccessible for Watson-Crick interaction, which suggests that this part of loop c is in a structurally constrained configuration, or buried in the tertiary structure or involved in tertiary interactions.     

Joining of DNA ends bearing non-matching 3'-overhangs

Using the Xenopus egg extract as an in vitro system for double strand break repair, the joining of DNA ends bearing non-complementary, homopolymeric, 4nt 3'-protruding single strands ("overhangs") was examined. Such 3'-overhangs can not be filled-in and cannot align and anneal by canonical base pair interactions, thus presenting a special challenge to the repair machinery. The results indicate that two such non-matching 3'-overhangs typically overlap by 2nt forming non-canonical base pairs, from which the filling-in of the remaining gaps is primed. The repair reaction is inhibited in Ku-depleted extracts. Unexpectedly, with some of the substrates the predominant repair products were joints with no nucleotide loss, suggesting that the two DNA ends aligned without overlap. However, the additional finding that an activity in the egg extract adds one or a few nucleotides to a fraction of the 3'-ends favors a model in which most or all of the zero-loss joints are the net result of 3'-overhang extension and 2bp overlap formation. The nucleotide addition reaction is stimulated by increasing the concentration of the complementary dNTP or ddNTP in the extract, suggesting a process templated by free nucleotides and the involvement of a DNA polymerase-like activity.     

Establishing broad generality of DNA catalysts for site-specific hydrolysis of single-stranded DNA

We recently reported that a DNA catalyst (deoxyribozyme) can site-specifically hydrolyze DNA on the minutes time scale. Sequence specificity is provided by Watson-Crick base pairing between the DNA substrate and two oligonucleotide binding arms that flank the 40-nt catalytic region of the deoxyribozyme. The DNA catalyst from our recent in vitro selection effort, 10MD5, can cleave a single-stranded DNA substrate sequence with the aid of Zn(2+) and Mn(2+) cofactors, as long as the substrate cleavage site encompasses the four particular nucleotides ATG^T. Thus, 10MD5 can cleave only 1 out of every 256 (4(4)) arbitrarily chosen DNA sites, which is rather poor substrate sequence tolerance. In this study, we demonstrated substantially broader generality of deoxyribozymes for site-specific DNA hydrolysis. New selection experiments were performed, revealing the optimality of presenting only one or two unpaired DNA substrate nucleotides to the N(40) DNA catalytic region. Comprehensive selections were then performed, including in some cases a key selection pressure to cleave the substrate at a predetermined site. These efforts led to identification of numerous new DNA-hydrolyzing deoxyribozymes, many of which require merely two particular nucleotide identities at the cleavage site (e.g. T^G), while retaining Watson-Crick sequence generality beyond those nucleotides along with useful cleavage rates. These findings establish experimentally that broadly sequence-tolerant and site-specific deoxyribozymes are readily identified for hydrolysis of single-stranded DNA.     

A novel double lesion in X-irradiated DNA consists of a strand break and a base modification.

A single free radical-initiating event can produce a pair of base lesions in DNA oligomers exposed to ionizing radiation. Whereas double base lesions have been observed previously, the present study shows that double lesions may sometimes consist of a base lesion and an associated strand break. The mechanism for the formation of double lesions is discussed. A redox process is postulated in which guanine is the source of the electron. It is suggested that double lesions may be formed in DNA either on adjacent nucleotides or, alternatively, on nucleotides separated by one, two or possibly more intervening nucleotides. It is hypothesized that intramolecular electron transfer facilitates the formation of double lesions on nonadjacent nucleotides.     

The intrinsic structure and stability of out-of-alternation base pairs in Z-DNA.

Alternating pyrimidine-purine sequences typically form Z-DNA, with the pyrimidines in the anti and purines in the syn conformations. The observation that dC and dT nucleotides can also adopt the syn conformation (i.e. the nucleotides are out-of-alternation) extends the range of sequences that can convert to this left-handed form of DNA. Here, we study the effects of placing two adjacent d(G*C) base pairs as opposed to a single d(G*C) base pair or two d(A*T) base pairs out-of-alternation by comparing the structure of d(m5CGGCm5CG)2with the previously published structures of d(m5CGGGm5CG)*d(m5CGCCm5CG) and d(m5CGATm5CG)2. A high buckle and loss of stacking interactions are observed as intrinsic properties of the out-of-alternation base pairs regardless of sequence and the context of the dinucleotide. From solution titrations, we find that the destabilizing effect of out-of-alternation d(G*C) base pairs are identical whether these base pairs are adjacent or isolated. We can therefore conclude that it is these intrinsic distortions in the structure of the base pairs and not neighboring effects that account for the inability of out-of-alternation base pairs to adopt the left-handed Z conformation.     

On loop folding in nucleic acid hairpin-type structures

In a series of studies, combining NMR, optical melting and T-jump experiments, it was found that DNA hairpins display a maximum stability when the loop part of the molecule comprises four or five nucleotide residues. This is in contrast with the current notion based on RNA hairpin studies, from which it had been established that a maximum hairpin stability is obtained for six or seven residues in the loop. Here we present a structural model to rationalize these observations. This model is based on the notion that to a major extent base stacking interactions determine the stability of nucleic acid conformations. The model predicts that loop folding in RNA is characterized by an extension of the base stacking at the 5'-side of the double helix by five or six bases; the remaining gap can then easily be closed by two nucleotides. Conversely, loop folding in DNA is characterized by extending base stacking at the 3'-side of the double helical stem by two or three residues; again bridging of the remaining gap can then be achieved by one or two nucleotides. As an example of loop folding in RNA the anticodon loop of yeast tRNAPhe is discussed. For the DNA hairpin formed by d(ATCCTAT4TAGGAT) it is shown that the loop structure obtained from molecular mechanics calculations obeys the above worded loop folding principles.     

Pre-steady-state kinetic studies of the fidelity of Sulfolobus solfataricus P2 DNA polymerase IV

Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4) is a thermostable archaeal enzyme and a member of the error-prone and lesion-bypass Y-family. In this paper, for the first time, the fidelity of a Y-family polymerase, Dpo4, was determined using pre-steady-state kinetic analysis of the incorporation of a single nucleotide into an undamaged DNA substrate 21/41-mer at 37 degrees C. We assessed single-turnover (with Dpo4 in molar excess over DNA) saturation kinetics for all 16 possible nucleotide incorporations. The fidelity of Dpo4 was estimated to be in the range of 10(-3)-10(-4). Interestingly, the ground-state binding affinity of correct nucleotides (70-230 microM) is 10-50-fold weaker than those of replicative DNA polymerases. Such a low affinity is consistent with the lack of interactions between Dpo4 and the bound nucleotides as revealed in the crystal structure of Dpo4, DNA, and a matched nucleotide. The affinity of incorrect nucleotides for Dpo4 is approximately 2-10-fold weaker than that of correct nucleotides. Intriguingly, the mismatched dCTP has an affinity similar to that of the matched nucleotides when it is incorporated against a pyrimidine template base flanked by a 5'-template guanine. The incoming dCTP likely skips the first available template base and base pairs with the 5'-template guanine, as observed in the crystal structure of Dpo4, DNA, and a mismatched nucleotide. The mismatch incorporation rates, regardless of the 5'-template base, were approximately 2-3 orders of magnitude slower than the incorporation rates for matched nucleotides, which is the predominant contribution to the fidelity of Dpo4.     

Precise binding of single-stranded DNA termini by human RAD52 protein

The human RAD52 protein, which exhibits a heptameric ring structure, has been shown to bind resected double strand breaks (DSBs), consistent with an early role in meiotic recombination and DSB repair. In this work, we show that RAD52 binds single-stranded and tailed duplex DNA molecules via precise interactions with the terminal base. When probed with hydroxyl radicals, ssDNA-RAD52 complexes exhibit a four-nucleotide repeat hypersensitivity pattern. This unique pattern is due to the interaction of RAD52 with either a 5' or a 3' terminus of the ssDNA, is sequence independent and is phased precisely from the terminal nucleotide. Hypersensitivity is observed over approximately 36 nucleotides, consistent with the length of DNA that is protected by RAD52 in nuclease protection assays. We propose that RAD52 binds DNA breaks via specific interactions with the terminal base, leading to the formation of a precisely organized ssDNA-RAD52 complex in which the DNA lies on an exposed surface of the protein. This protein-DNA arrangement may facilitate the DNA-DNA interactions necessary for RAD52-mediated annealing of complementary DNA strands.     

Differential Mismatch Recognition Specificities of Eukaryotic MutS Homologs,MutSα and MutSβ

In eukaryotes, the recognition of the DNA postreplication errors and initiation of the mismatch repair is carried out by two MutS homologs: MutSα and MutSβ. MutSα recognizes base mismatches and 1 to 2 unpaired nucleotides whereas MutSβ recognizes longer insertion-deletion loops (IDLs) with 1 to 15 unpaired nucleotides as well as certain mismatches. Results from molecular dynamics simulations of native MutSβ:IDL-containing DNA and MutSα:mismatch DNA complexes as well as complexes with swapped DNA substrates provide mechanistic insight into how the differential substrate specificities are achieved by MutSα and MutSβ, respectively. Our simulations results suggest more extensive interactions between MutSβ and IDL-DNA and between MutSα and mismatch-containing DNA that suggest corresponding differences in stability. Furthermore, our simulations suggest more expanded mechanistic details involving a different degree of bending when DNA is bound to either MutSα or MutSβ and a more likely opening of the clamp domains when noncognate substrates are bound. The simulation results also provide detailed information on key residues in MutSβ and MutSα that are likely involved in recognizing IDL-DNA and mismatch-containing DNA, respectively.     

Full-length Dengue virus RNA-dependent RNA polymerase-RNA/DNA complexes: stoichiometries, intrinsic affinities, cooperativities, base, and conformational specificities

Fundamental aspects of interactions of the Dengue virus type 3 full-length polymerase with the single-stranded and double-stranded RNA and DNA have been quantitatively addressed. The polymerase exists as a monomer with an elongated shape in solution. In the absence of magnesium, the total site size of the polymerase-ssRNA complex is 26 ± 2 nucleotides. In the presence of Mg(2+), the site size increases to 29 ± 2 nucleotides, indicating that magnesium affects the enzyme global conformation. The enzyme shows a preference for the homopyrimidine ssRNAs. Positive cooperativity in the binding to homopurine ssRNAs indicates that the type of nucleic acid base dramatically affects the enzyme orientation in the complex. Both the intrinsic affinity and the cooperative interactions are accompanied by a net ion release. The polymerase binds the dsDNA with an affinity comparable with the ssRNAs affinity, indicating that the binding site has an open conformation in solution. The lack of detectable dsRNA or dsRNA-DNA hybrid affinities indicates that the entry to the binding site is specific for the sugar-phosphate backbone and/or conformation of the duplex.     

Dependence of Deformability of Geometries and Characteristics of Intramolecular Hydrogen Bonds in Canonical 2′-Deoxyribonucleotides on DNA Conformations

Abstract

The molecular structure and deformability (with respect to average geometry) of methyl ethers of canonical 2′-deoxyribonucleotides thymidine-5′-phosphate (mTMP), 2-deoxycytidine-5′- phosphate (mCMP), 2-deoxyadenosine-5′-phosphate (mAMP) and 2′-deoxyguanosine-5′- phosphate (mGMP) in different types of DNA have been calculated using B3LYP/cc-pvdz method. Comparison of energy at equilibrium conformations of nucleotides and conformations with torsion angles of backbone fixed to average values for different types of DNA reveals that incorporation of nucleotides to A-DNA macromolecules requires the minimum amount of deformation energy. Therefore, this type of DNA should be the least strained from viewpoint of intramolecular deformations of monomers. Modeling of environmental effects within the PCM approach reveals that the immersion of nucleotides in polar medium results in significant decrease of energy differences between anti conformers of all DNTs and syn conformers of mGMP This also leads to reduction by almost a half nucleotides' deformation energy facilitating formation of DNA macromolecule. Change of DNTs conformation causes switch between different types of intramolecular H bonds. Every type of DNA possesses unique set of intramolar hydrogen bonds in nucleotides.