Comparative analysis of structural dynamics and allosteric mechanisms of RecA/Rad51 family proteins: Integrated atomistic MD simulation and network-based analysis

Homologous recombination plays a key role in double-strand break repair, stalled replication fork repair, and meiosis. The RecA/Rad51 family recombinases catalyze the DNA strand invasion reaction that occurs during homologous recombination. However, the high sequence differences between homologous groups have hindered the thoroughly studies of this ancient protein family. The dynamic mechanisms of the family, particularly at the residual level, remain poorly understood. In this work, five representative RecA/Rad51 recombinase family members from all major kingdoms of living organisms: prokaryotes, eukaryotes, archaea, and viruses, were selected to explore the molecular mechanisms behind their conserved biological significance. A variety of techniques, including all-atom molecular dynamics simulation, perturbation response scanning, and protein structure network analysis, were used to examine the flexibility and correlation of protein domains, distribution of sensors and effectors and conserved hub residues. Furthermore, the potential communication routes between the ATP-binding region and the DNA-binding region of each recombinase were identified. Our results demonstrate the conserved molecular dynamics of these recombinases in the early stage of homologous recombination, including cooperative motions between regions, conserved sensing and effecting functional residue distribution, and conserved hub residues. Meanwhile, the unique ATP-DNA communication routes of each recombinase was also revealed. These results provide new insights into the mechanism of RecA/Rad51 family proteins, and provide new theoretical guidance for the development of allosteric inhibitors and the application of RecA/Rad51 family proteins.