DNA Recombinase Proteins, their Function and Structure in the Active Form, a Computational Study

Homologous recombination is a crucial sequence of reactions in all cells for the repair of double strand DNA (dsDNA) breaks. While it was traditionally considered as a means for generating genetic diversity, it is now known to be essential for restart of collapsed replication forks that have met a lesion on the DNA template (Cox et al., 2000). The central stage of this process requires the presence of the DNA recombinase protein, RecA in bacteria, RadA in archaea, or Rad51 in eukaryotes, which leads to an ATP-mediated DNA strand-exchange process. Despite many years of intense study, some aspects of the biochemical mechanism, and structure of the active form of recombinase proteins are not well understood. Our theoretical study is an attempt to shed light on the main structural and mechanistic issues encountered on the RecA of the e-coli, the RecA of the extremely radio resistant Deinococcus Radiodurans (promoting an inverse DNA strand-exchange repair), and the homolog human Rad51. The conformational changes are analyzed for the naked enzymes, and when they are linked to ATP and ADP. The average structures are determined over 2ns time scale of Langevian dynamics using a collision frequency of 1.0 ps(sup -1). The systems are inserted in an octahedron periodic box with a 10 Angstrom buffer of water molecules explicitly described by the TIP3P model. The corresponding binding free energies are calculated in an implicit solvent using the Poisson-Boltzmann solvent accessible surface area, MM-PBSA model. The role of the ATP is not only in stabilizing the interaction RecA-DNA, but its hydrolysis is required to allow the DNA strand-exchange to proceed. Furthermore, we extended our study, using the hybrid QM/MM method, on the mechanism of this chemical process. All the calculations were performed using the commercial code Amber 9.