Functional impact of missense variants in BRCA1 predicted by supervised learning

Many individuals tested for inherited cancer susceptibility at the BRCA1 gene locus are discovered to have variants of unknown clinical significance (UCVs). Most UCVs cause a single amino acid residue (missense) change in the BRCA1 protein. They can be biochemically assayed, but such evaluations are time-consuming and labor-intensive. Computational methods that classify and suggest explanations for UCV impact on protein function can complement functional tests. Here we describe a supervised learning approach to classification of BRCA1 UCVs. Using a novel combination of 16 predictive features, the algorithms were applied to retrospectively classify the impact of 36 BRCA1 C-terminal (BRCT) domain UCVs biochemically assayed to measure transactivation function and to blindly classify 54 documented UCVs. Majority vote of three supervised learning algorithms is in agreement with the assay for more than 94% of the UCVs. Two UCVs found deleterious by both the assay and the classifiers reveal a previously uncharacterized putative binding site. Clinicians may soon be able to use computational classifiers such as those described here to better inform patients. These classifiers can be adapted to other cancer susceptibility genes and systematically applied to prioritize the growing number of potential causative loci and variants found by large-scale disease association studies.     

BRCA1 and BRCA2 missense variants of high and low clinical significance influence lymphoblastoid cell line post-irradiation gene expression

The functional consequences of missense variants in disease genes are difficult to predict. We assessed if gene expression profiles could distinguish between BRCA1 or BRCA2 pathogenic truncating and missense mutation carriers and familial breast cancer cases whose disease was not attributable to BRCA1 or BRCA2 mutations (BRCAX cases). 72 cell lines from affected women in high-risk breast ovarian families were assayed after exposure to ionising irradiation, including 23 BRCA1 carriers, 22 BRCA2 carriers, and 27 BRCAX individuals. A subset of 10 BRCAX individuals carried rare BRCA1/2 sequence variants considered to be of low clinical significance (LCS). BRCA1 and BRCA2 mutation carriers had similar expression profiles, with some subclustering of missense mutation carriers. The majority of BRCAX individuals formed a distinct cluster, but BRCAX individuals with LCS variants had expression profiles similar to BRCA1/2 mutation carriers. Gaussian Process Classifier predicted BRCA1, BRCA2 and BRCAX status, with a maximum of 62% accuracy, and prediction accuracy decreased with inclusion of BRCAX samples carrying an LCS variant, and inclusion of pathogenic missense carriers. Similarly, prediction of mutation status with gene lists derived using Support Vector Machines was good for BRCAX samples without an LCS variant (82–94%), poor for BRCAX with an LCS (40–50%), and improved for pathogenic BRCA1/2 mutation carriers when the gene list used for prediction was appropriate to mutation effect being tested (71–100%). This study indicates that mutation effect, and presence of rare variants possibly associated with a low risk of cancer, must be considered in the development of array-based assays of variant pathogenicity.     

Characterization of common BRCA1 and BRCA2 variants

Many missense variants identified in BRCA1 and BRCA2, two genes responsible for the majority of hereditary breast and ovarian cancer, are of unclear clinical significance. Characterizing the significance of such variants is important for medical management of patients in whom they are identified. The aim of this study was to characterize eight of the most common reported missense mutations in BRCA1 and BRCA2 occurring in patients tested for hereditary risk of breast and ovarian cancers. The prevalence of each variant in a control population, co-segregation of the variant with cancer within families, location of the variant within the gene, the nature of the amino acid substitution and conservation of the wild-type amino acid among species were considered. In a control population, the BRCA1 variants M1652I, R1347G, and S1512I, were each observed at a frequency of 4.08%, 2.04%, and 2.04%, respectively, and the BRCA2 variants A2951T, V2728I, and D1420Y, were seen at 1.02%, 0.68%, and 0.34%, respectively. Although the BRCA2 variants T598A and R2034C were not seen in this group of controls, other clinical and published observations indicate that these variants are not deleterious. Based on epidemiological and biological criteria, we therefore conclude that the BRCA1 missense mutations R1347G, S1512I and M1652I, and the BRCA2 missense mutations T598A, D1420Y, R2034C, V2728I, and A2951T, are not deleterious mutations.     

Structural consequences of a cancer-causing BRCA1-BRCT missense mutation

The integrity of the carboxyl-terminal BRCT repeat region is critical for BRCA1 tumor suppressor function; however, the molecular details of how a number of clinically derived BRCT missense mutations affect BRCA1 function remain largely unknown. Here we assess the structural response of the BRCT tandem repeat domain to a well characterized, cancer-associated single amino acid substitution, Met-1775 --> Arg-1775. The structure of BRCT-M1775R reveals that the mutated side chain is extruded from the protein hydrophobic core, thereby altering the protein surface. Charge-charge repulsion, rearrangement of the hydrophobic core, and disruption of the native hydrogen bonding network at the interface between the two BRCT repeats contribute to the conformational instability of BRCT-M1775R. Destabilization and global unfolding of the mutated BRCT domain at physiological temperatures explain the pleiotropic molecular and genetic defects associated with the BRCA1-M1775R protein.     

Characterization of cancer‐linked BRCA1‐BRCT missense variants and their interaction with phosphoprotein targets

The breast cancer tumor suppressor protein BRCA1 is involved in DNA repair and cell cycle control. Mutations at the two C‐terminal tandem (BRCT) repeats of BRCA1 detected in breast tumor patients were identified either to lower the stability of the BRCT domain and/or to disrupt the interaction of BRCT with phoshpopeptides. The aim of this study was to analyze five BRCT pathogenic mutations for their effect on structural integrity and protein stability. For this purpose, the five cancer‐associated BRCT mutants: V1696L, M1775K, M1783T, V1809F, and P1812A were cloned in suitable prokaryotic protein production vectors, and the recombinant proteins were purified in soluble and stable form for further biophysical studies. The biophysical analysis of the secondary structure and the thermodynamic stability of the wild‐type, wt, and the five mutants of the BRCT domain were performed by Circular Dichroism Spectroscopy (CD) and Differential Scanning Microcalorimetry (DSC), respectively. The binding capacity of the wt and mutant BRCT with (pBACH1/BRIP1) and pCtIP were measured by Isothermal Titration Calorimetry (ITC). The experimental results demonstrated that the five mutations of the BRCT domain: (i) affected the thermal unfolding temperature as well as the unfolding enthalpy of the domain, to a varying degree depending upon the induced destabilization and (ii) altered and/or abolished their affinity to synthetic pBACH1/BRIP1 and pCtIP phosphopeptides by affecting the structural integrity of the BRCT active sites. The presented experimental results are one step towards the elucidation of the effect of various missense mutations on the structure and function of BRCA1‐BRCT. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Expression of cancer related BRCA1 missense variants decreases MMS-induced recombination in Saccharomyces cerevisiae without altering its nuclear localization

BRCA1 tumor suppressor gene is found mutated in familial breast and ovarian cancer. Most cancer related mutations were found located at the RING (Really Interesting New Gene) and at the BRCT (BRca1 C-Terminal) domain. However, 20 y after its identification, the biological role of BRCA1 and which domains are more relevant for tumor suppression are still being elucidated. We previously reported that expression of BRCA1 cancer related variants in the RING and BRCT domain increases spontaneous homologous recombination in yeast indicating that BRCA1 may interact with yeast DNA repair/recombination. To finally demonstrate whether BRCA1 interacts with yeast DNA repair, we exposed yeast cells expressing BRCA1wt, the cancer-related variants C-61G and M1775R to different doses of the alkylating agent methyl methane-sulfonate (MMS) and then evaluated the effect on survival and homologous recombination. Cells expressing BRCA1 cancer variants were more sensitive to MMS and less inducible to recombination as compared to cell expressing BRCA1wt. Moreover, BRCA1-C61G and -M1775R did not change their nuclear localization form as compared to the BRCA1wt or the neutral variant R1751Q indicating a difference in the DNA damage processing. We propose a model where BRCA1 cancer variants interact with the DNA double strand break repair pathways producing DNA recombination intermediates, that maybe less repairable and decrease MMS-induced recombination and survival. Again, this study strengthens the use of yeast as model system to characterize the mechanisms leading to cancer in humans carrying the BRCA1 missense variant.     

BRCA1 variants in a family study of African-American and Latina women

We sequenced the entire coding region of BRCA1 to improve our understanding of the frequency and nature of BRCA1 variants in African-American and Latina women identified from a multiethnic cohort in Los Angeles, California. The study included 109 African-American and 140 Latina sibships from families with two or more cases of breast or ovarian cancer among first-degree relatives. BRCA1 was sequenced in 278 breast or ovarian cancer cases and 229 unaffected sisters. The proportion of cases with known disease-causing mutations was low (0.72, 95% confidence interval: 0–1.7%). In total, 33 sequence variants were identified, including two protein truncation mutations, one deletion, and six silent and 24 missense variants. Two novel rare variants were identified that appeared to act as benign polymorphisms. Four rare variants may be unique to women of African descent based on existing literature, and three have been described exclusively in Latina women. The frequency of common variants was similar for cases and controls, but the frequency of common variants for African-American women significantly differed from those previously described for Caucasian women. We believe this to be the largest study of high-risk African-American and Latina women sequenced for variants in the BRCA1 gene to date.     

Thermal unfolding of human BRCA1 BRCT-domain variants

Missense mutations at the BRCT domain of human BRCA1 protein have been associated with an elevated risk for hereditary breast/ovarian cancer. They have been shown to affect the binding site and they have also been proposed to affect domain stability, severely hampering the protein's tumor suppressor function. In order to assess the impact of various such mutations upon the stability and the function of the BRCT domain, heat-induced denaturation has been employed to study the thermal unfolding of variants M1775R and R1699W, which have been linked with the disease, as well as of V1833M, which has been reported for patients with a family history. Calorimetric and circular dichroism results reveal that in pH 9.0, 5 mM borate buffer, 200 mM NaCl, analogously to wild type BRCT, all three variants undergo partial thermal unfolding to a denatured state, which retains most of the native's structural characteristics. With respect to wild-type BRCT, the mutation M1775R induces the most severe effects especially upon the thermostability, while R1699W also has a strong impact. On the other hand, the thermal unfolding of variant V1833M is only moderately affected relative to wild-type BRCT. Moreover, isothermal titration calorimetric measurements reveal that contrary to M1775R and R1699W variants, V1833M binds to BACH1 and CtIP phosphopeptides.     

Analysis of natural variants of the human immunodeficiency virus type 1 gag-pol frameshift stem-loop structure

Human immunodeficiency virus type 1 uses ribosomal frameshifting for translation of the Gag-Pol polyprotein. Frameshift activities are thought to be tightly regulated. Analysis of gag p1 sequences from 270 plasma virions identified in 64% of the samples the occurrence of polymorphism that could lead to changes in thermodynamic stability of the stem-loop. Expression in Saccharomyces cerevisiae of p1-beta-galactosidase fusion proteins from 10 representative natural stem-loop variants and three laboratory mutant constructs (predicted the thermodynamic stability [Delta G degrees] ranging from -23.0 to -4.3 kcal/mol) identified a reduction in frameshift activity of 13 to 67% compared with constructs with the wild-type stem-loop (Delta G degrees, -23.5 kcal/mol). Viruses carrying stem-loops associated with greater than 60% reductions in frameshift activity presented profound defects in viral replication. In contrast, viruses with stem-loop structures associated with 16 to 42% reductions in frameshift efficiency displayed no significant viral replication deficit.     

Development and application of a method for counterselectable in-frame deletion in Clostridium perfringens

Many pathogenic clostridial species produce toxins and enzymes. To facilitate genome-wide identification of virulence factors and biotechnological application of their useful products, we have developed a markerless in-frame deletion method for Clostridium perfringens which allows efficient counterselection and multiple-gene disruption. The system comprises a galKT gene disruptant and a suicide galK plasmid into which two fragments of a target gene for in-frame deletion are cloned. The system was shown to be accurate and simple by using it to disrupt the alpha-toxin gene of the organism. It was also used to construct of two different virulence-attenuated strains, ΗΝ1303 and HN1314: the former is a disruptant of the virRS operon, which regulates the expression of virulence factors, and the latter is a disruptant of the six genes encoding the α, θ, and κ toxins; a clostripain-like protease; a 190-kDa secretory protein; and a putative cell wall lytic endopeptidase. Comparison of the two disruptants in terms of growth ability and the background levels of secreted proteins showed that HN1314 is more useful than ΗΝ1303 as a host for the large-scale production of recombinant proteins.     

Impact of BRCA1 BRCT domain missense substitutions on phosphopeptide recognition

The BRCA1 BRCT domain binds pSer-x-x-Phe motifs in partner proteins to regulate the cellular response to DNA damage. Approximately 120 distinct missense variants have been identified in the BRCA1 BRCT through breast cancer screening, and several of these have been linked to an increased cancer risk. Here we probe the structures and peptide-binding activities of variants that affect the BRCA1 BRCT phosphopeptide-binding groove. The results obtained from the G1656D and T1700A variants illustrate the role of Ser1655 in pSer recognition. Mutations at Arg1699 (R1699W and R1699Q) significantly reduce peptide binding through loss of contacts to the main chain of the Phe(+3) residue and, in the case of R1699W, to a destabilization of the BRCT fold. The R1835P and E1836K variants do not dramatically reduce peptide binding, in spite of the fact that these mutations significantly alter the structure of the walls of the Phe(+3) pocket.     

Identification of BRCA1-IRIS, a BRCA1 locus product

Breast cancer is the most common malignancy among women, and mutations in the BRCA genes produce increased susceptibility to these malignancies in certain families. Here we identify BRCA1-IRIS as a 1,399-amino-acid BRCA1 gene product encoded by an uninterrupted open reading frame that extends from codon 1 of the known BRCA1 open reading frame to a termination point 34 triplets into intron 11. Unlike full-length BRCA1 (p220), BRCA1-IRIS is exclusively chromatin-associated, fails to interact with BARD1 in vivo or in vitro and exhibits unique nuclear immunostaining. Unlike BRCA1FL (or p220), BRCA1-IRIS also co-immunoprecipitated with DNA-replication-licensing proteins and with known replication initiation sites. Suppression of BRCA1-IRIS expression hindered the normal departure of geminin from pre-replication complexes, and depressed the rate of cellular DNA replication and possibly initiation-related synthesis. In contrast, BRCA1-IRIS overexpression stimulated DNA replication. These data imply that endogenous BRCA1-IRIS positively influences the DNA replication initiation machinery.     

Detection of protein folding defects caused by BRCA1-BRCT truncation and missense mutations

Most cancer-associated BRCA1 mutations identified to date result in the premature translational termination of the protein, highlighting a crucial role for the C-terminal, BRCT repeat region in mediating BRCA1 tumor suppressor function. However, the molecular and genetic effects of missense mutations that map to the BRCT region remain largely unknown. Using a protease-based assay, we directly assessed the sensitivity of the folding of the BRCT domain to an extensive set of truncation and single amino acid substitutions derived from breast cancer screening programs. The protein can tolerate truncations of up to 8 amino acids, but further deletion results in drastic BRCT folding defects. This molecular phenotype can be correlated with an increased susceptibility to disease. A cross-validated computational assessment of the BRCT mutation data base suggests that as much as half of all BRCT missense mutations contribute to BRCA1 loss of function and disease through protein-destabilizing effects. The coupled use of proteolytic methods and computational predictive methods to detect mutant BRCA1 conformations at the protein level will augment the efficacy of current BRCA1 screening protocols, especially in the absence of clinical data that can be used to discriminate deleterious BRCT missense mutations from benign polymorphisms.