In Silico and Structure-Based Assessment of Similar Variants Discovered in Tandem Repeats of BRCT Domains of BRCA1 and BARD1 To Characterize the Folding Pattern

Structural dynamics of clinically-reported VUS in the BARD1 ARD-BRCT region to predict the molecular basis of alterations

The functional domains of BARD1, comprise the Ankyrin Repeat Domain (ARD), C-Terminal domains (BRCTs), and a linker between ARD and the BRCTs, which are known to bind to Cleavage stimulation Factor complex-subunit of 50 kDa (CstF-50). The pathogenic mutation Q564H in the BARD1 ARD-linker-BRCT region has been reported to abrogate the binding between BARD1 and CstF-50. Intermediate penetrance variants of BARD1 are associated with the occurrence of breast cancer. Therefore, seven missense variants of unknown significance (VUS), L447V, P454L, N470S, V507M, I509T, C557S, and Q564H of BARD1, reported in the ARD domain and the linker region were evaluated via molecular dynamics (MD) simulations. The mutants revealed statistically significantly different distributions of RMSD (root mean square deviation), residuewise RMSF (root mean square fluctuation), Rg (radius of gyration), SASA (solvent accessible surface area), and COM (centre of mass)-to-COM distance between the ARD and the BRCT repeat, between the wild type and each mutant. The secondary structural composition of the mutants was slightly altered relative to that of the wild type. However, the reported in-silico based prediction require further validation using in-vitro, biophysical and structure-based approachCommunicated by Ramaswamy H. Sarma.

In silico investigation on the mutational analysis of BRCA1-BARD1 RING domains and its effect on nucleosome recognition and ubiquitination

The BRCA1-BARD1 complex is a crucial tumor suppressor E3 ubiquitin ligase involved in DNA double-stranded break repair. The BRCA1-BARD1 RING domains interact with UBE2D3 through the BRCA1 interface and this complex flexibly tether to the nucleosome core particle (NCP), where BRCA1 and BARD1 interacts with histone H2A and H2B of NCP. Mutations in the BRCA1-BARD1 RING domains have been linked to familial breast and ovarian cancer. Seven mutations were analyzed to understand their effect on the binding interface of protein partners and changes in conformational dynamics. Molecular dynamics simulations revealed that mutant complexes were less conformationally flexible than the wildtype complex. Protein-protein interaction profiling showed the importance of specific molecular interactions, hotspot and hub residues, and some of these were lost in the mutant complexes. Two mutations (BRCA1L51W-K65R and BARD1C53W) hindered significant interaction between protein partners and may prevent signaling for ubiquitination of histones in NCP and other cellular targets. The structural compactness and reduced significant interaction in mutant complexes may be the possible reason of preventing ubiquitination and hinder DNA repair, resulting cancer.

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

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.