Formation of nucleoprotein RecA filament on single-stranded DNA. Analysis by stepwise increase in ligand complexity

How Enzymes, Proteins, and Antibodies Recognize Extended DNAs; General Regularities

X-ray analysis cannot provide quantitative estimates of the relative contribution of non-specific, specific, strong, and weak contacts of extended DNA molecules to their total affinity for enzymes and proteins. The interaction of different enzymes and proteins with long DNA and RNA at the quantitative molecular level can be successfully analyzed using the method of the stepwise increase in ligand complexity (SILC). The present review summarizes the data on stepwise increase in ligand complexity (SILC) analysis of nucleic acid recognition by various enzymes-replication, restriction, integration, topoisomerization, six different repair enzymes (uracil DNA glycosylase, Fpg protein from Escherichia coli, human 8-oxoguanine-DNA glycosylase, human apurinic/apyrimidinic endonuclease, RecA protein, and DNA-ligase), and five DNA-recognizing proteins (RNA helicase, human lactoferrin, alfa-lactalbumin, human blood albumin, and IgGs against DNA). The relative contributions of structural elements of DNA fragments "covered" by globules of enzymes and proteins to the total affinity of DNA have been evaluated. Thermodynamic and catalytic factors providing discrimination of unspecific and specific DNAs by these enzymes on the stages of primary complex formation following changes in enzymes and DNAs or RNAs conformations and direct processing of the catalysis of the reactions were found. General regularities of recognition of nucleic acid by DNA-dependent enzymes, proteins, and antibodies were established.

How Human H1 Histone Recognizes DNA

Linker H1 histone is one of the five main histone proteins (H1, H2A, H2B, H3, and H4), which are components of chromatin in eukaryotic cells. Here we have analyzed the patterns of DNA recognition by free H1 histone using a stepwise increase of the ligand complexity method; the affinity of H1 histone for various single- and double-stranded oligonucleotides (d(pN)n; n = 1-20) was evaluated using their competition with 12-mer [32P]labeled oligonucleotide and protein-oligonucleotide complex delaying on nitrocellulose membrane filters. It was shown that minimal ligands of H1 histone (like other DNA-dependent proteins and enzymes) are different mononucleotides (dNMPs; Kd = (1.30 ± 0.2) × 10-2 M). An increase in the length of single-stranded (ss) homo- and hetero-oligonucleotides (d(pA)n, d(pT)n, d(pC)n, and d(pN)n with different bases) by one nucleotide link regardless of their bases, leads to a monotonic increase in their affinity by a factor of f = 3.0 ± 0.2. This factor f corresponds to the Kd value = 1/f characterizing the affinity of one nucleotide of different ss d(pN)n for H1 at n = 2-6 (which are covered by this protein globule) is approximately 0.33 ± 0.02 M. The affinity of five out of six DNA nucleotide units is approximately 25 times lower than for one of the links. The affinity of duplexes of complementary homo- and hetero-d(pN)20 is only 1.3-3.3-fold higher in comparison with corresponding ss oligonucleotides. H1 histone forms mainly weak additive contacts with internucleoside phosphate groups of ssDNAs and one chain of double-stranded DNAs, but not with the bases.

How human serum albumin recognizes DNA and RNA

We show here for the first time that HSA possesses two nucleic acid-(NA) binding sites and we estimated the relative contributions of the nucleotide links of (pN)n to their total affinity for these binding sites with higher and lower affinity for NAs. The minimal ligands of these binding sites are orthophosphate (Kd=3.0 and 20.0 mm), various dNMPs (5.6-400 μm and 0.063-18 mm) and different rNMPs (4.9-30 μm and 14-250 μm). Maximal contribution to the total affinity of all NAs to the first and second sites was observed for one nucleotide and was remarkably lower for three additional nucleotide units of (pN)n (n=1-4) with a significant decrease in the contribution at n=5-6, and at n≥7-8 all dependencies reached plateaus. For d(pA)n and r(pA)n a relatively gradual decrease in the contribution to the affinity at n=1-6 was observed, while several d(pN)n, demonstrated a sharp increase in the contribution at n=2-4. Finally, all (pN)n>10 demonstrated high affinity for the first (1.4-150 nm) and the second (80-2400 nm) sites of HSA. Double-stranded NAs showed significantly lower affinity comparing with single-stranded ligands. The thermodynamic parameters characterizing the specific contribution of every nucleotide link of all (pN)1-9 (ΔG°) to their total affinity for HSA were estimated.

Mechanism of the formation of the RecA-ssDNA nucleoprotein filament structure: a coarse-grained approach

In prokaryotes, the RecA protein catalyzes the repair and strand exchange of double-stranded DNA. RecA binds to single-stranded DNA (ssDNA) and forms a presynaptic complex in which the protein polymerizes around the ssDNA to form a right-handed helical nucleoprotein filament structure. In the present work, the mechanism for the formation of the RecA-ssDNA filament structure is modeled using coarse-grained molecular dynamics simulations. Information from the X-ray structure was used to model the protein itself but not its interactions; the interactions between the protein and the ssDNA were modeled solely by electrostatic, aromatic, and repulsive energies. For the present study, the monomeric, dimeric, and trimeric units of RecA and 4, 8, and 11 NT-long ssDNA, respectively, were studied. Our results indicate that monomeric RecA is not sufficient for nucleoprotein filament formation; rather, dimeric RecA is the elementary binding unit, with higher multimeric units of RecA facilitating filament formation. Our results reveal that loop region flexibility at the primary binding site of RecA is essential for it to bind the incoming ssDNA, that the aromatic residues present in the loop region play an important role in ssDNA binding, and that ATP may play a role in guiding the ssDNA by changing the electrostatic potential of the RecA protein.

How human IgGs against myelin basic protein (MBP) recognize oligopeptides and MBP

Myelin basic protein (MBP) is a major protein of myelin-proteolipid shell of axons, and it plays an important role in pathogenesis of multiple sclerosis. In the literature, there are no data on how antibodies recognize different protein antigens including MBP. A stepwise increase in ligand complexity was used to estimate the relative contributions of virtually every amino acid residue (AA) of a specific 12-mer LSRFSWGAEGQK oligopeptide corresponding to immunodominant sequence of MBP to the light chains and to intact anti-MBP IgGs from sera of patients with multiple sclerosis. It was shown that the minimal ligands of the light chains of IgGs are many different free AAs (K_d  = 0.51-0.016 M), and each free AA interacts with the specific subsite of the light chain intended for recognition of this AA in specific LSRFSW oligopeptide. A gradual transition from Leu to LSRFSWGAEGQK leads to an increase in the affinity from 10^-1 to 2.3 × 10^-4  M because of additive interactions of the light chain with 6 AAs of this oligopeptide and then the affinity reaches plateau. The contributions of 6 various AAs to the affinity of the oligopeptide are different (K_d , M): 0.71 (S), 0.44 (R), 0.14 (F), 0.17 (S), and 0.62 (W). Affinity of nonspecific oligopeptides to the light chains of IgGs is significantly lower. Intact MBP interacts with both light and heavy chains of IgGs demonstrating 192-fold higher affinity than the specific oligopeptide. It is a first quantitative analysis of the mechanism of proteins recognition by antibodies. The thermodynamic model was constructed to describe the interactions of IgGs with MBP. The data obtained can be very useful for understanding how antibodies against many different proteins can recognize these proteins.

How human IgGs against DNA recognize oligonucleotides and DNA

In the literature, there are no available data on how anti-DNA antibodies recognize DNA. In the present work, to study the molecular mechanism of DNA recognition by antibodies, we have used anti-DNA IgGs from blood sera of patients with multiple sclerosis. A stepwise increase in ligand complexity approach was used to estimate the relative contributions of virtually every nucleotide unit of different single- (ss) and double-stranded (ds) oligonucleotides to their affinity for IgG fraction having high affinity to DNA-cellulose. DNA-binding site disposed on the heavy chain demonstrates higher affinity to different dNMPs (K_d  = 0.63μM-3.8μM) than the site located on the light chain (28μM-170μM). The heavy and light chains interact independently forming relatively strong contacts with 2 to 4 nucleotides of short homo- and hetero-d(pN)_2-9 . Then the increase in the affinity of different d(pN)_n became minimal, and at n ≥ 8 to 9, all dependencies reached plateaus: approximately 3.2nM to 20nM and approximately 200nM to 460nM for the heavy and light chains, respectively. A similar situation was observed for different ribooligonucleotides, in which their affinity is 6-fold to 100-fold lower than that for d(pN)_n . Transition from ss to ds d(pN)_n leads to a moderate increase in affinity of ligands to DNA-binding site of heavy chains, while light chains demonstrate the same affinity for ss and ds d(pN)_n . Long supercoiled DNA interacts with both heavy and light chains with affinity of approximately 10-fold higher than that for short oligonucleotides. The thermodynamic models were constructed to describe the interactions of IgGs light and heavy chains with DNA.

Recombinational branch migration by the RadA/Sms paralog of RecA in Escherichia coli

RadA (also known as 'Sms') is a highly conserved protein, found in almost all eubacteria and plants, with sequence similarity to the RecA strand exchange protein and a role in homologous recombination. We investigate here the biochemical properties of the E. coli RadA protein and several mutant forms. RadA is a DNA-dependent ATPase, a DNA-binding protein and can stimulate the branch migration phase of RecA-mediated strand transfer reactions. RadA cannot mediate synaptic pairing between homologous DNA molecules but can drive branch migration to extend the region of heteroduplex DNA, even without RecA. Unlike other branch migration factors RecG and RuvAB, RadA stimulates branch migration within the context of the RecA filament, in the direction of RecA-mediated strand exchange. We propose that RadA-mediated branch migration aids recombination by allowing the 3’ invading strand to be incorporated into heteroduplex DNA and to be extended by DNA polymerases.

DOI:http://dx.doi.org/10.7554/eLife.10807.001

Damage to the DNA of a cell can cause serious harm, and so cells have several ways in which they can repair DNA. Most of these processes rely on the fact that each of the two strands that make up a DNA molecule can be used as a template to build the other strand. However, this is not possible if both strands of the DNA break in the same place. This form of damage can be repaired in a process called homologous recombination, which uses an identical copy of the broken DNA molecule to repair the broken strands. As a result, this process can only occur during cell division shortly after a cell has duplicated its DNA.

One important step of homologous recombination is called strand exchange. This involves one of the broken strands swapping places with part of the equivalent strand in the intact DNA molecule. To do so, the strands of the intact DNA molecule separate in the region that will be used for the repair, and the broken strand can then use the other non-broken DNA strand as a template to replace any missing sections of DNA. The region of the intact DNA molecule where the strands need to separate often grows during this process: this is known as branch migration. In bacteria, a protein called RecA plays a fundamental role in controlling strand exchange, but there are other, similar proteins whose roles in homologous recombination are less well known.

Cooper and Lovett have now purified one of these proteins, called RadA, from the Escherichia coli species of bacteriato study how it affects homologous recombination. This revealed that RadA can bind to single-stranded DNA and stimulate branch migration to increase the rate of homologous recombination. Further investigation revealed that RadA allows branch migration to occur even when RecA is missing, but that RadA is unable to begin strand exchange if RecA is not present. The process of branch migration stabilizes the DNA molecules during homologous recombination and may also allow the repaired DNA strand to engage the machinery that copies DNA.

Cooper and Lovett also used genetic techniques to alter the structure of specific regions of RadA and found out which parts of the protein affect the ability of RadA to stimulate branch migration. Future challenges are to find out what effect RadA has on the structure of RecA and how RadA promotes branch migration.

DOI:http://dx.doi.org/10.7554/eLife.10807.002

AFAL: a web service for profiling amino acids surrounding ligands in proteins

With advancements in crystallographic technology and the increasing wealth of information populating structural databases, there is an increasing need for prediction tools based on spatial information that will support the characterization of proteins and protein-ligand interactions. Herein, a new web service is presented termed amino acid frequency around ligand (AFAL) for determining amino acids type and frequencies surrounding ligands within proteins deposited in the Protein Data Bank and for assessing the atoms and atom-ligand distances involved in each interaction (availability: http://structuralbio.utalca.cl/AFAL/index.html ). AFAL allows the user to define a wide variety of filtering criteria (protein family, source organism, resolution, sequence redundancy and distance) in order to uncover trends and evolutionary differences in amino acid preferences that define interactions with particular ligands. Results obtained from AFAL provide valuable statistical information about amino acids that may be responsible for establishing particular ligand-protein interactions. The analysis will enable investigators to compare ligand-binding sites of different proteins and to uncover general as well as specific interaction patterns from existing data. Such patterns can be used subsequently to predict ligand binding in proteins that currently have no structural information and to refine the interpretation of existing protein models. The application of AFAL is illustrated by the analysis of proteins interacting with adenosine-5'-triphosphate.

Recognition of specific and nonspecific DNA by human lactoferrin

The general principles of recognition of nucleic acids by proteins are among the most exciting problems of molecular biology. Human lactoferrin (LF) is a remarkable protein possessing many independent biological functions, including interaction with DNA. In human milk, LF is a major DNase featuring two DNA-binding sites with different affinities for DNA. The mechanism of DNA recognition by LF was studied here for the first time. Electrophoretic mobility shift assay and fluorescence measurements were used to probe for interactions of the high-affinity DNA-binding site of LF with a series of model-specific and nonspecific DNA ligands, and the structural determinants of DNA recognition by LF were characterized quantitatively. The minimal ligands for this binding site were orthophosphate (K(i)  = 5 mM), deoxyribose 5'-phosphate (K(i)  = 3 mM), and different dNMPs (K(i)  = 0.56-1.6 mM). LF interacted additionally with 9-12 nucleotides or nucleotide pairs of single- and double-stranded ribo- and deoxyribooligonucleotides of different lengths and sequences, mainly through weak additive contacts with internucleoside phosphate groups. Such nonspecific interactions of LF with noncognate single- and double-stranded d(pN)(10) provided ~6 to ~7.5 orders of magnitude of the enzyme affinity for any DNA. This corresponds to the Gibbs free energy of binding (ΔG(0)) of -8.5 to -10.0 kcal/mol. Formation of specific contacts between the LF and its cognate DNA results in an increase of the DNA affinity for the enzyme by approximately 1 order of magnitude (K(d)  = 10 nM; ΔG(0)  ≈ -11.1 kcal/mol). A general function for the LF affinity for nonspecific d(pN)(n) of different sequences and lengths was obtained, giving the K(d) values comparable with the experimentally measured ones. A thermodynamic model was constructed to describe the interactions of LF with DNA.

Structural, thermodynamic, and kinetic basis for the activities of some nucleic acid repair enzymes

X-ray structural analysis provides no quantitative estimate of the relative contribution of specific and nonspecific or strong and weak interactions to the total affinity of enzymes for nucleic acids. We have shown that the interaction between enzymes and long nucleic acids at the molecular level can be successfully analyzed by the method of stepwise increase in ligand complexity (SILC). In the present review we summarize our studies of human uracil DNA glycosylase and apurinic/apyrimidinic endonuclease, E. coli 8-oxoguanine DNA glycosylase and RecA protein using the SILC approach. The relative contribution of structural (X-ray analysis data), thermodynamic, and catalytic factors to the discrimination of specific and nonspecific DNA by these enzymes at the stages of complex formation, the following changes in DNA and enzyme conformations and especially the catalysis of the reactions is discussed.

Main factors providing specificity of repair enzymes

Specific and nonspecific DNA complex formation with human uracil-DNA glycosylase, 8-oxoguanine-DNA glycosylase, and apurine/apyrimidine endonuclease, as well as with E. coli 8-oxoguanine-DNA glycosylase and RecA protein was analyzed using the method of stepwise increase in DNA-ligand complexity. It is shown that high affinity of these enzymes to any DNA (10(-4)-10(-8) M) is provided by a large number of weak additive contacts mainly with DNA internucleoside phosphate groups and in a less degree with bases of nucleotide links "covered" by protein globules. Enzyme interactions with specific DNA links are comparable in efficiency with weak unspecific contacts and provide only for one-two orders of affinity (10(-1)-10(-2) M), but these contacts are extremely important at stages of DNA and enzyme structural adaptation and catalysis proper. Only in the case of specific DNA individual for each enzyme alterations in DNA structure provide for efficient adjustment of reacting enzyme atoms and DNA orbitals with accuracy up to 10-15° and, as a result, for high reaction rate. Upon transition from nonspecific to specific DNA, reaction rate (k(cat)) increases by 4-8 orders of magnitude. Thus, stages of DNA and enzyme structural adaptation as well as catalysis proper are the basis of specificity of repair enzymes.

Thermodynamic and kinetic basis for recognition and repair of 8-oxoguanine in DNA by human 8-oxoguanine-DNA glycosylase

We have used a stepwise increase in ligand complexity approach to estimate the relative contributions of the nucleotide units of DNA containing 7,8-dihydro-8-oxoguanine (oxoG) to its total affinity for human 8-oxoguanine DNA glycosylase (OGG1) and construct thermodynamic models of the enzyme interaction with cognate and non-cognate DNA. Non-specific OGG1 interactions with 10–13 nt pairs within its DNA-binding cleft provides approximately 5 orders of magnitude of its affinity for DNA (ΔG° approximately −6.7 kcal/mol). The relative contribution of the oxoG unit of DNA (ΔG° approximately −3.3 kcal/mol) together with other specific interactions (ΔG° approximately −0.7 kcal/mol) provide approximately 3 orders of magnitude of the affinity. Formation of the Michaelis complex of OGG1 with the cognate DNA cannot account for the major part of the enzyme specificity, which lies in the k_cat term instead; the rate increases by 6–7 orders of magnitude for cognate DNA as compared with non-cognate one. The k_cat values for substrates of different sequences correlate with the DNA twist, while the K_M values correlate with ΔG° of the DNA fragments surrounding the lesion (position from −6 to +6). The functions for predicting the K_M and k_cat values for different sequences containing oxoG were found.

Binding selectivity of RecA to a single stranded DNA, a computational approach

Homologous recombination (HR) is the major DNA double strand break repair pathway which maintains the genomic integrity. It is fundamental for the survivability and functionality of all organisms. One of the initial steps in HR is the formation of the nucleoprotein filament composed by a single stranded DNA chain surrounded by the recombinases protein. The filament orchestrates the search for an undamaged homologue, as a template for the repair process. Our theoretical study was aimed at elucidating the selectivity of the interaction between a monomer of the recombinases enzyme in the Escherichia coli, EcRecA, the bacterial homologue of human Rad51, with a series of oligonucleotides of nine bases length. The complex, equilibrated for 20 ns with Langevian dynamics, was inserted in a periodic box with a 8 Å buffer of water molecules explicitly described by the TIP3P model. The absolute binding free energies are calculated in an implicit solvent using the Poisson-Boltzmann (PB) and the generalized Born (GB) solvent accessible surface area, using the MM-PB(GB)SA model. The solute entropic contribution is also calculated by normal mode analysis. The results underline how a significant contribution of the binding free energy is due to the interaction with the Arg196, a critical amino acid for the activity of the enzyme. The study revealed how the binding affinity of EcRecA is significantly higher toward dT₉ rather than dA₉, as expected from the experimental results.

Uracil-DNA glycosylase: Structural, thermodynamic and kinetic aspects of lesion search and recognition

Uracil appears in DNA as a result of cytosine deamination and by incorporation from the dUTP pool. As potentially mutagenic and deleterious for cell regulation, uracil must be removed from DNA. The major pathway of its repair is initiated by uracil-DNA glycosylases (UNG), ubiquitously found enzymes that hydrolyze the N-glycosidic bond of deoxyuridine in DNA. This review describes the current understanding of the mechanism of uracil search and recognition by UNG. The structure of UNG proteins from several species has been solved, revealing a specific uracil-binding pocket located in a DNA-binding groove. DNA in the complex with UNG is highly distorted to allow the extrahelical recognition of uracil. Thermodynamic studies suggest that UNG binds with appreciable affinity to any DNA, mainly due to the interactions with the charged backbone. The increase in the affinity for damaged DNA is insufficient to account for the exquisite specificity of UNG for uracil. This specificity is likely to result from multistep lesion recognition process, in which normal bases are rejected at one or several pre-excision stages of enzyme-substrate complex isomerization, and only uracil can proceed to enter the active site in a catalytically competent conformation. Search for the lesion by UNG involves random sliding along DNA alternating with dissociation-association events and partial eversion of undamaged bases for initial sampling.

Interaction of single-stranded DNA with the second DNA-binding site of the RecA nucleoprotein filament

Bacterial RecA protein is a prototype of ATP-dependent homologous recombinases found ubiquitously from bacteriophages up to human beings. When RecA filament is forming on single-stranded DNA in the presence of ATP, it initiates the strand exchange reaction with homologous double-stranded DNA. Among three phases of the reaction (the search for homology, the three-stranded structure annealing in conjunction with the switch of pairing, and the strand displacement) the first one is the most enigmatic and least studied. As commonly recognized, this phase is directed by a special (stretched) filament structure and does not required any additional consumption of energy in ATP hydrolysis. The novel approaches in the study of strand exchange reaction, using short oligonucleotides as DNA substrates and sensitive methods for a real-time monitoring of the reaction suggest that all three phases of the reaction depend on the ATP hydrolysis.

Structural modeling of sequence specificity by an autoantibody against single-stranded DNA

11F8 is a sequence-specific pathogenic anti-single-stranded (ss)DNA autoantibody isolated from a lupus prone mouse. Site-directed mutagenesis of 11F8 has shown that six binding site residues (R31VH, W33VH, L97VH, R98VH, Y100VH, and Y32VL) contribute 80% of the free energy for complex formation. Mutagenesis results along with intermolecular distances obtained from fluorescence resonance energy transfer were implemented here as restraints to model docking between 11F8 and the sequence-specific ssDNA. The model of the complex suggests that aromatic stacking and two sets of bidentate hydrogen bonds between binding site arginine residues (R31VH and R96VH) and loop nucleotides provide the molecular basis for high affinity and specificity. In part, 11F8 utilizes the same ssDNA binding motif of Y32VL, H91VL, and an aromatic residue in the third complementarity-determining region to recognize thymine-rich sequences as do two anti-ssDNA autoantibodies crystallized in complex with thymine. R31SVH is a dominant somatic mutation found in the J558 germline sequence that is implicated in 11F8 sequence specificity. A model of the mutant R31S11F8.ssDNA complex suggests that different interface contacts occur when serine replaces arginine 31 at the binding site. The modeled contacts between the R31S11F8 mutant and thymine are closely related to those observed in other anti-ssDNA binding antibodies, while we find additional contacts between 11F8 and ssDNA that involve amino acids not utilized by the other antibodies. These data-driven 11F8.ssDNA models provide testable hypotheses concerning interactions that mediate sequence specificity in 11F8 and the effects of somatic mutation on ssDNA recognition.

Mechanism of RecA-mediated homologous recombination revisited by single molecule nanomanipulation

The mechanisms of RecA-mediated three-strand homologous recombination are investigated at the single-molecule level, using magnetic tweezers. Probing the mechanical response of DNA molecules and nucleoprotein filaments in tension and in torsion allows a monitoring of the progression of the exchange in real time, both from the point of view of the RecA-bound single-stranded DNA and from that of the naked double-stranded DNA (dsDNA). We show that strand exchange is able to generate torsion even along a molecule with freely rotating ends. RecA readily depolymerizes during the reaction, a process presenting numerous advantages for the cell's 'protein economy' and for the management of topological constraints. Invasion of an untwisted dsDNA by a nucleoprotein filament leads to an exchanged duplex that remains topologically linked to the exchanged single strand, suggesting multiple initiations of strand exchange on the same molecule. Overall, our results seem to support several important assumptions of the monomer redistribution model.