Cancer-related mutations in BRCA1-BRCT cause long-range structural changes in protein-protein binding sites: a molecular dynamics study

Simulation and Computational Study of RING Domain Mutants of BRCA1 and Ube2k in AD/PD Pathophysiology

Lysine-based post-translational modification (PTM) such as acylation, acetylation, deamination, methylation, SUMOylation, and ubiquitination has proven to be a major regulator of gene expression, chromatin structure, protein stability, protein-protein interaction, protein degradation, and cellular localization. However, besides all the PTMs, ubiquitination stands as the second most common PTM after phosphorylation that is involved in the etiology of neurodegenerative diseases (NDDs) namely, Alzheimer's disease (AD) and Parkinson's disease (PD). NDDs are characterized by the accumulation of misfolded protein aggregates in the brain that lead to disease-related gene mutation and irregular protein homeostasis. The ubiquitin-proteasome system (UPS) is in charge of degrading these misfolded proteins, which involve an interplay of E1, E2, E3, and deubiquitinase enzymes. Impaired UPS has been commonly observed in NDDs and E3 ligases are the key members of the UPS, thus, dysfunction of the same can accelerate the neurodegeneration process. Therefore, the aim of this study is firstly, to find E3 ligases that are common in both AD and PD through data mining. Secondly, to study the impact of mutation on its structure and function. The study deciphered 74 E3 ligases that were common in both AD and PD. Later, 10 hub genes were calculated of which protein-protein interaction, pathway enrichment, lysine site prediction, domain, and motif analysis were performed. The results predicted BRCA1, PML, and TRIM33 as the top three putative lysine-modified E3 ligases involved in AD and PD pathogenesis. However, based on structural characterization, BRCA1 was taken further to study RING domain mutation that inferred K32Y, K32L, K32C, K45V, K45Y, and K45G as potential mutants that alter the structural and functional ability of BRCA1 to interact with Ube2k, E2-conjugating enzyme. The most probable mutant observed after molecular dynamics simulation of 50 ns is K32L. Therefore, our study concludes BRCA1, a potential E3 ligase common in AD and PD, and RING domain mutation at sites K32 and K45 possibly disturbs its interaction with its E2, Ube2k.

From Cognitive Agents to Cognitive Systems: Theoretical, Methodological, and Empirical Developments of van Gelder's (1998) "Dynamical Hypothesis"

Over two decades have passed since the publication of van Gelder's (1998) "dynamical hypothesis." In that paper, van Gelder proposed that cognitive agents were not digital computers-per the representational computational approach-but dynamical systems. The evolution of the dynamical hypothesis was driven by parallel advances in three areas. Theoretically, a deeper understanding of genetics, biology, neuroscience, and cognitive science inspired questions about how systems within each domain dynamically interact and extend their effects across spatiotemporal scales. Methodologically, more sophisticated and domain-general tools allowed researchers to discover, model, and quantify system dynamics, structure, and patterns across multiple scales to generate a more comprehensive system-level understanding of behaviors. Empirically, we can analyze a system's behavior while preserving its natural dynamics, revealing evidence that the reductionist approach leads to an incomplete understanding of the components and the overall system. Researchers have traditionally reduced a complex system into its component processes and assumed that the parts can be recombined to explain the whole. These three advances fundamentally altered our understanding of a "cognitive agent:" How their behaviors are driven by long-range coordination across multiple processes, how the interdependent and nested structure of interacting variables produces behaviors that are greater than the sum of its parts, and how environmental constraints shape adaptive yet stable behavioral patterns.

DNA damage response signaling pathways and targets for radiotherapy sensitization in cancer

Radiotherapy is one of the most common countermeasures for treating a wide range of tumors. However, the radioresistance of cancer cells is still a major limitation for radiotherapy applications. Efforts are continuously ongoing to explore sensitizing targets and develop radiosensitizers for improving the outcomes of radiotherapy. DNA double-strand breaks are the most lethal lesions induced by ionizing radiation and can trigger a series of cellular DNA damage responses (DDRs), including those helping cells recover from radiation injuries, such as the activation of DNA damage sensing and early transduction pathways, cell cycle arrest, and DNA repair. Obviously, these protective DDRs confer tumor radioresistance. Targeting DDR signaling pathways has become an attractive strategy for overcoming tumor radioresistance, and some important advances and breakthroughs have already been achieved in recent years. On the basis of comprehensively reviewing the DDR signal pathways, we provide an update on the novel and promising druggable targets emerging from DDR pathways that can be exploited for radiosensitization. We further discuss recent advances identified from preclinical studies, current clinical trials, and clinical application of chemical inhibitors targeting key DDR proteins, including DNA-PKcs (DNA-dependent protein kinase, catalytic subunit), ATM/ATR (ataxia–telangiectasia mutated and Rad3-related), the MRN (MRE11-RAD50-NBS1) complex, the PARP (poly[ADP-ribose] polymerase) family, MDC1, Wee1, LIG4 (ligase IV), CDK1, BRCA1 (BRCA1 C terminal), CHK1, and HIF-1 (hypoxia-inducible factor-1). Challenges for ionizing radiation-induced signal transduction and targeted therapy are also discussed based on recent achievements in the biological field of radiotherapy.

Germline Missense Variants in BRCA1: New Trends and Challenges for Clinical Annotation

Genetic testing allows for the identification of germline DNA variations, which are associated with a significant increase in the risk of developing breast cancer (BC) and ovarian cancer (OC). Detection of a BRCA1 or BRCA2 pathogenic variant triggers several clinical management actions, which may include increased surveillance and prophylactic surgery for healthy carriers or treatment with the PARP inhibitor therapy for carriers diagnosed with cancer. Thus, standardized validated criteria for the annotation of BRCA1 and BRCA2 variants according to their pathogenicity are necessary to support clinical decision-making and ensure improved outcomes. Upon detection, variants whose pathogenicity can be inferred by the genetic code are typically classified as pathogenic, likely pathogenic, likely benign, or benign. Variants whose impact on function cannot be directly inferred by the genetic code are labeled as variants of uncertain clinical significance (VUS) and are evaluated by multifactorial likelihood models that use personal and family history of cancer, segregation data, prediction tools, and co-occurrence with a pathogenic BRCA variant. Missense variants, coding alterations that replace a single amino acid residue with another, are a class of variants for which determination of clinical relevance is particularly challenging. Here, we discuss current issues in the missense variant classification by following a typical life cycle of a BRCA1 missense variant through detection, annotation and information dissemination. Advances in massively parallel sequencing have led to a substantial increase in VUS findings. Although the comprehensive assessment and classification of missense variants according to their pathogenicity remains the bottleneck, new developments in functional analysis, high throughput assays, data sharing, and statistical models are rapidly changing this scenario.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Promiscuous proteins are commonly observed in biological systems, such as modular domains that recognize phosphopeptides during signal transduction. The use of phosphopeptides and compounds with phosphate groups as inhibitors to protein–protein interactions have attracted increasing interest for years. By using atomistic molecular dynamics simulations, we are able to perform detailed analyses of the dihedral space to explore protein fluctuation upon ligand binding to better understand promiscuous molecular recognition. Free energy calculation can further provide insights into the mechanism of binding, including both enthalpic and entropic contributions for molecular recognition, which assist in inhibitor design. Our calculation results show that pre-rigidifying a ligand is not always advantageous, suggesting the challenge in retaining optimized intermolecular interactions in pre-rigidified ligand. Instead, certain flexible ligands with multiple binding conformations can reduce entropic penalty, and therefore improves binding affinity. According to our computations, we can introduce new intermolecular interactions to flexible ligand to strengthen attractions while maintaining small entropic penalty by retaining its plasticity in the bound conformation. The study might cast light on a new general strategy for designing inhibitors targeting promiscuous modular domains and protein–protein interactions.

Splicing analysis of 14 BRCA1 missense variants classifies nine variants as pathogenic

Pathogenic germline mutations in the BRCA1 gene predispose carriers to early onset breast and ovarian cancer. Clinical genetic screening of BRCA1 often reveals variants with uncertain clinical significance, complicating patient and family management. Therefore, functional examinations are urgently needed to classify whether these uncertain variants are pathogenic or benign. In this study, we investigated 14 BRCA1 variants by in silico splicing analysis and mini-gene splicing assay. All 14 alterations were missense variants located within the BRCT domain of BRCA1 and had previously been examined by functional analysis at the protein level. Results from a validated mini-gene splicing assay indicated that nine BRCA1 variants resulted in splicing aberrations leading to truncated transcripts and thus can be considered pathogenic (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5072C>T/p.Thr1691Ile, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala), whereas five BRCA1 variants had no effect on splicing (c.4985T>C/p.Phe1662Ser, c.5072C>A/p.Thr1691Lys, c.5153G>C/p.Trp1718Ser, c.5154G>T/p.Trp1718Cys, and c.5333A>G/p.Asp1778Gly). Eight of the variants having an effect on splicing (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala) were previously determined to have no or an uncertain effect on the protein level, whereas one variant (c.5072C>T/p.Thr1691Ile) were shown to have a strong effect on the protein level as well. In conclusion, our study emphasizes that in silico splicing prediction and mini-gene splicing analysis are important for the classification of BRCA1 missense variants located close to exon/intron boundaries.

Quantifying intracellular protein binding thermodynamics during mechanotransduction based on FRET spectroscopy

Mechanical force modulates myriad cellular functions including migration, alignment, proliferation, and gene transcription. Mechanotransduction, the transmission of mechanical forces and its translation into biochemical signals, may be mediated by force induced protein conformation changes, subsequently modulating protein signaling. For the paxillin and focal adhesion kinase interaction, we demonstrate that force-induced changes in protein complex conformation, dissociation constant, and binding Gibbs free energy can be quantified by lifetime-resolved fluorescence energy transfer microscopy combined with intensity imaging calibrated by fluorescence correlation spectroscopy. Comparison with in vitro data shows that this interaction is allosteric in vivo. Further, spatially resolved imaging and inhibitor assays show that this protein interaction and its mechano-sensitivity are equal in the cytosol and in the focal adhesions complexes indicating that the mechano-sensitivity of this interaction must be mediated by soluble factors but not based on protein tyrosine phosphorylation.

The role of engineering approaches in analysing cancer invasion and metastasis

In the past decade, novel materials, probes and tools have enabled fundamental and applied cancer researchers to take a fresh look at the complex problem of tumour invasion and metastasis. These new tools, which include imaging modalities, controlled but complex in vitro culture conditions, and the ability to model and predict complex processes in vivo, represent an integration of traditional with novel engineering approaches; and their potential effect on quantitatively understanding tumour progression and invasion looks promising.

The role of BRCA1 in DNA double-strand repair: past and present

The breast cancer type 1 susceptibility protein (BRCA1) is involved in several important cellular pathways, including DNA damage repair, chromatin remodeling and checkpoint activation. The BRCA1 tumor suppression function has been attributed to its role in homologous recombination damage repair. In this review, historical facts concerning BRCA1, together with recent research advances regarding our understanding of the BRCA1 interacting proteins that are involved in, homologous recombination (HR) double strand break (DBS) repair and how these interacting proteins maintain chromosomal integrity, are discussed. In addition, this review poses the questions as to what extent HR repair cannot be properly fulfilled when breast cancer related mutations in the BRCA1 gene occur and how the recent and excessive studied poly-ADP ribose polymerase (PARP) inhibiting therapy approach links with the proposed tumor suppression function of the different BRCA1 domains.

The BRCA1 Breast Cancer Suppressor: Regulation of Transport, Dynamics, and Function at Multiple Subcellular Locations

Inherited mutations in the BRCA1 gene predispose to a higher risk of breast/ovarian cancer. The BRCA1 tumor suppressor is a 1863 amino acid protein with multiple protein interaction domains that facilitate its roles in regulating DNA repair and maintenance, cell cycle progression, transcription, and cell survival/apoptosis. BRCA1 was first identified as a nuclear phosphoprotein, but has since been shown to contain different transport sequences including nuclear export and nuclear localization signals that enable it to shuttle between specific sites within the nucleus and cytoplasm, including DNA repair foci, centrosomes, and mitochondria. BRCA1 nuclear transport and ubiquitin E3 ligase enzymatic activity are tightly regulated by the BRCA1 dimeric binding partner BARD1 and further modulated by cancer mutations and diverse signaling pathways. This paper will focus on the transport, dynamics, and multiple intracellular destinations of BRCA1 with emphasis on how regulation of these events has impact on, and determines, a broad range of important cellular functions.

Revised AMBER parameters for bioorganic phosphates

We report AMBER force field parameters for biological simulations involving phosphorylation of serine, threonine or tyrosine. The initial parameters used RESP fitting for the atomic partial charges and standard values for all other parameters such as Lennard-Jones coefficients. These were refined with the aid of a thermodynamic cycle consisting of experimentally determined pKa values, solvation energies from molecular dynamics free energy simulations, and gas phase basicities from QM calculations. A polarization energy term was included to account for the charge density change between the gas-phase and solution, and solvation free energies were determined using thermodynamic integration. Parameter adjustment is required to obtain consistent thermodynamic results with better balanced electrostatic interactions between water and the phosphate oxygens. To achieve this we modified the phosphate oxygen radii. A thermodynamically consistent parameter set can be derived for monoanions and requires an increase of the van der Waals phosphate oxygen radii of approximately 0.09 Å. Larger, residue-specific radii appear to be needed for dianions. The revised parameters developed here should be of particular interest for environments where simulations of multiple protonation states may be of interest.

Mechanistic insights into phosphopeptide--BRCT domain association: preorganization, flexibility, and phosphate recognition

Promiscuous proteins are commonly observed in biological systems, for example, in modular domains that recognize phosphopeptides during signal transduction. This promiscuous recognition is of fundamental interest in chemistry and biology but is challenging when designing phosphopeptides in silico for cell biology studies. To investigate promiscuous recognition and binding processes of phosphopeptides and the modular domain, we selected a domain essential in breast cancer-the breast-cancer-associated protein 1 (BRCA1) C-terminal (BRCT) repeats as our model system. We performed molecular dynamics simulations and detailed analyses of the dihedral space to study protein fluctuation and conformational changes with phosphopeptide binding. We also studied the association processes of phosphorylated and unphosphorylated peptides using Brownian dynamics with a coarse-grained model. We found that the BRCT domain is preorganized for phosphopeptide binding but has a moderate arrangement of side chains to form complexes with various types of phosphopeptides. Phosphopeptide binding restricts the system motion in general, while the nonpolar phosphopeptide becomes more flexible in the bound state. Our analysis found that the BRCT domain utilizes different mechanisms, usually termed lock and key, induced-fit, and population-shift/conformational-selection models, to recognize peptides with different features. Brownian dynamics simulations revealed that the charged phosphate group may not always accelerate peptide association processes, but it helps the phosphopeptide orient into binding pockets accurately and stabilizes the complex. This work provides insights into molecular recognition in the promiscuous protein system.

Structure-Function Of The Tumor Suppressor BRCA1

BRCA1, a multi-domain protein, is mutated in a large percentage of hereditary breast and ovarian cancers. BRCA1 is most often mutated in three domains or regions: the N-terminal RING domain, exons 11-13, and the BRCT domain. The BRCA1 RING domain is responsible for the E3 ubiquitin ligase activity of BRCA1 and mediates interactions between BRCA1 and other proteins. BRCA1 ubiquitinates several proteins with various functions. The BRCA1 BRCT domain binds to phosphoproteins with specific sequences recognized by both BRCA1 and ATM/ATR kinases. Structural studies of the RING and BRCT domains have revealed the molecular basis by which cancer causing mutations impact the functions of BRCA1. While no structural data is available for the amino acids encoded by exons 11-13, multiple binding sites and functional domains exist in this region. Many mutations in exons 11-13 have deleterious effects on the function of these domains. In this mini-review, we examine the structure-function relationships of the BRCA1 protein and the relevance to cancer progression.

Multimodel assessment of BRCA1 mutations in Taiwanese (ethnic Chinese) women with early-onset, bilateral or familial breast cancer

Although evidence suggests an importance of genetic factors in the development of breast cancer in Taiwanese (ethnic Chinese) women, including a high incidence of early-onset and secondary contralateral breast cancer, a major breast cancer predisposition gene, BRCA1, has not been well studied in this population. In fact, the carcinogenic impacts of many genetic variants of BRCA1 are unknown and classified as variants of uncertain significance (VUS). It is therefore important to establish a method to characterize the BRCA1 VUSs and understand their role in Taiwanese breast cancer patients. Accordingly, we developed a multimodel assessment strategy consisting of a prescreening portion and a validated functional assay to study breast cancer patients with early-onset, bilateral or familial breast cancer. We found germ-line BRCA1 mutations in 11.1% of our cohort and identified one novel missense mutation, c.5191C>A. Two genetic variants were initially classified as VUSs (c.1155C>T and c.5191C>A). c.1155C>T is not predicted to be deleterious in the prescreening portion of our assessment strategy. c.5191C>A, on the other hand, causes p.T1691K, which is predicted to have high deleterious probability because of significant structural alteration, a high deleterious score in the predictive programs and, clinically, triple negative characteristics in breast tumors. This mutant is confirmed by transcription activation and yeast growth-inhibition assays. In conclusion, we show as high a prevalence of germ-line BRCA1 mutation in high-risk Taiwanese patients as in Caucasians and demonstrate a useful strategy for studying BRCA1 VUSs.

Stability of domain structures in multi-domain proteins

Multi-domain proteins have many advantages with respect to stability and folding inside cells. Here we attempt to understand the intricate relationship between the domain-domain interactions and the stability of domains in isolation. We provide quantitative treatment and proof for prevailing intuitive ideas on the strategies employed by nature to stabilize otherwise unstable domains. We find that domains incapable of independent stability are stabilized by favourable interactions with tethered domains in the multi-domain context. Stability of such folds to exist independently is optimized by evolution. Specific residue mutations in the sites equivalent to inter-domain interface enhance the overall solvation, thereby stabilizing these domain folds independently. A few naturally occurring variants at these sites alter communication between domains and affect stability leading to disease manifestation. Our analysis provides safe guidelines for mutagenesis which have attractive applications in obtaining stable fragments and domain constructs essential for structural studies by crystallography and NMR.

Mechanism of PhosphoThreonine/Serine Recognition and Specificity for Modular Domains from All-atom Molecular Dynamics

Background

Phosphopeptide-binding domains mediate many vital cellular processes such as signal transduction and protein recognition. We studied three well-known domains important for signal transduction: BRCT repeats, WW domain and forkhead-associated (FHA) domain. The first two recognize both phosphothreonine (pThr) and phosphoserine (pSer) residues, but FHA has high specificity for pThr residues. Here we used molecular dynamics (MD) simulations to reveal how FHA exclusively chooses pThr and how BRCT and WW recognize both pThr/pSer. The work also investigated the energies and thermodynamic information of intermolecular interactions.

Results

Simulations carried out included wide-type and mutated systems. Through analysis of MD simulations, we found that the conserved His residue defines dual loops feature of the FHA domain, which creates a small cavity reserved for only the methyl group of pThr. These well-organized loop interactions directly response to the pThr binding selectivity, while single loop (the 2nd phosphobinding site of FHA) or in combination with α-helix (BRCT repeats) or β-sheet (WW domain) fail to differentiate pThr/pSer.

Conclusions

Understanding the domain pre-organizations constructed by conserved residues and the driving force of domain-phosphopeptide recognition provides structural insight into pThr specific binding, which also helps in engineering proteins and designing peptide inhibitors.

The (non)malignancy of cancerous amino acidic substitutions

The process of natural selection acts both on individual organisms within a population and on individual cells within an organism as they develop into cancer. In this work, we have taken a first step toward understanding the differences in selection pressures exerted on the human genome under these disparate circumstances. Focusing on single amino acid substitutions, we have found that cancer-related mutations (CRMs) are frequent in evolutionarily conserved sites, whereas single amino acid polymorphisms (SAPs) tend to appear in sites having a more relaxed evolutionary pressure. Those CRMs classed as cancer driver mutations show greater enrichment for conserved sites than passenger mutations. Consistent with this, driver mutations are enriched for sites annotated as key functional residues and their neighbors, and are more likely to be located on the surface of proteins than expected by chance. Overall the pattern of CRM and polymorphism is remarkably similar, but we do see a clear signal indicative of diversifying selection for disruptive amino acid substitutions in the cancer driver mutations. The ultimate consequence of the appearance of those mutations must be advantageous for the tumor cell, leading to cell population-growth and migration events similar to those seen in natural ecosystems.

Positional analyses of BRCA1-dependent expression in Saccharomyces cerevisiae

Mutations in BRCA1 account for a significant proportion of hereditary breast and ovarian cancers, but analysis of BRCA1 function is complicated by pleiotropic effects and binding partners (Pol II holoenzyme and transcription factors, chromatin remodelers, recombination complexes and E3 ligases). In vertebrate cells, efforts to elucidate BRCA1 transcriptional effects have focused on specific genes or restricted portions of the genome-limiting analyses of BRCA1 effects on adjoining DNA sequences and along chromosome lengths. Here, we use microarray analyses on the genetically tractable yeast cell system to elucidate BRCA1-dependent genomewide positional effects on both gene induction and repression. Yeast responses may be of clinical relevance based on findings that BRCA1 severely diminishes yeast growth kinetics but that BRCA1 mutated at sites identified from breast tumors is no longer able to retard yeast cell growth kinetics. Our analysis suggests that BRCA1 acts through both transcription factors to upregulate specific loci and chromatin remodeling complexes to effect global changes in gene expression. BRCA1 also exhibits gene repression activities. Cluster-functional analysis reveals that these repressed factors are required for mitotic stability and provide a novel molecular explanation for the conditional lethality observed between BRCA1 and chromosome segregation genes.