Molecular dynamics simulations reveal the effect of mutations in the RING domains of BRCA1-BARD1 complex and its relevance to the prognosis of breast cancer

Abstract

The N-terminal RING-RING domain of BRCA1-BARD1 is an E3 ubiquitin ligase complex that plays a critical role in tumor suppression through DNA double stranded repair mechanism. Mutations in the BRCA1-BARD1 heterodimer RING domains were found to have an association with breast and ovarian cancer by a way of hampering the E3 ubiquitin ligase activity. Herein, the molecular mechanism of interaction, conformational change due to the specific mutations on the BRCA1-BARD1 complex at atomic level has been examined by employing molecular modeling techniques. Sixteen mutations have been selected for the study. Molecular dynamics simulation results reveal that the mutant complexes have more local perturbation with a high residual fluctuation in the zinc binding sites and central helix. A few of the BRCA1 (V11A, I21V, I42V, R71G, I31M and L51W) mutants have been experimentally identified that do not impair E3 ligase activity, display an enhanced number of H-bonds and non-bonded contacts at the interacting interface as revealed by MD simulation. The mutation of BRCA1 (C61G, C64Y, C39Y and C24R) and BARD1 (C53W, C71Y and C83R) zinc binding residues displayed a smaller number of significant H-bonds, other interactions and also loss of some of the hotspot residues. Additionally, most of the mutant complexes display relatively lower electrostatic energy, H-bonding and total stabilizing energy as compared to wild-type. The current study attempts to unravel the role of BRCA1-BARD1 mutations and delineates the structural and conformational dynamics in the progression of breast cancer. Communicated by Ramaswamy H. Sarma