Molecular Modeling of the Amino-Terminal Zinc Ring Domain of BRCA1
Total Page:16
File Type:pdf, Size:1020Kb
CANCER RESEARCH56. 2539-2545. June I, 19961 Molecular Modeling of the Amino-Terminal Zinc Ring Domain of BRCA1 Rachelle J. Bienstock,1 Tom Darden, Roger Wiseman, Lee Pedersen, and J. Carl Barrett National institute of Environmental Health Sciences, Laboratories of Quantitative and Computational Biology fR. J. B.. T. D., L P.] and Molecular Carcinogenesis IR. W.. J. C. B.J, Research Triangle Park. North Carolina 27709 ABSTRACT in Fig. 1. Although a common functionality for the zinc ring domain has not been identified, it appears to interact with DNA either through The equine herpes virus zinc ring domain nuclear magnetic resonance direct binding or indirectly by mediating protein-protein interactions. structure was used for homology-based modeling of the amino-terminal The zinc ring domain of several proteins has been shown to exhibit zinc ring domain of the BRCA1 breast and ovarian cancer susceptibility gene. The zinc ring domain of BRCAJ is of particular interest because it DNA-binding activity (10) or nuclear localization (11—13). is the location of significant and frequently occurring missense (Cys6tGly, The BRCA1 zinc ring domain and its structure and function are of Cys―Gly,andCys―Tyr)andframeshift (lSSdelAG) mutations observed particular interest, because it is the location of some of the most In several high-risk kindreds. The BRCA1 zinc ring domain possesses frequently occurring mutations linked to breast and ovarian cancer. 54% sequence similarity with the equine herpes virus zinc ring domain. The l85delAG mutation (deletion oftwo nucleotide base pairs in exon The model structure undergoes little conformational variance after 140 ps 2) in the zinc ring domain has been shown to occur in 1 in 100 of solvated molecular dynamics. This model proposes BRCA1 zinc ring Ashkenazi Jewish individuals (14, 15). This is in comparison with the domain residues that may play a role in DNA binding and/or protein estimated I in 833 frequency of all BRCAI mutations in the general protein interactions. These predictions provide a point of departure for population (14, 15). In this article, a three-dimensional model for the the design of mutants to probe BRCA1 zinc ring domain functionality. zinc ring domain of BRCA1 is presented. This model proposes po tential BRCAI zinc ring domain residues that may play a role in the INTRODUCTION function of this domain, either binding directly to DNA or mediating DNA binding through protein-protein interactions. BRCAJ, the breast and ovarian cancer susceptibility gene on human chromosome l7q21, was recently identified (1, 2). Mutations in this gene account for approximately 95% of families with increased mci MATERIALS AND METHODS dence of early onset breast and ovarian cancer and 45% of families with increased incidence of breast cancer alone (3, 4). BRCAJ is Sequence Alignment. The three-dimensional structures of two zinc ring believed to encode a tumor suppressor protein, because most muta domains, the EHV2and the acute PML proto-oncoprotein, have recently been tions appear to result in loss of function. This gene is expressed in determined by multidimensional NMR (8, 9, 1 1). The sequence alignment several tissues, including breast and ovary, and encodes a I863-amino between the PML and EHV zinc ring domains is shown in Fig. 2 along with the secondary structural elements identified for each protein from the solved acid protein. NMR structures. Several criteria were used to select the EHV ring domain The amino-terminal domain of BRCA1 possesses high sequence structure rather than the PML structure as the basis for the BRCA1 model. The homology with a zinc finger topology (1). This zinc finger topology, BRCAI zinc ring domain, for example, possesses greater sequence similarity C3HC4,isreferredtoasthe“zinc-ringdomain,―“RINGfinger,―or to EHV than to PML. Using a sequence alignment of Smith and Waterman (16) A-Box (5) and is found in several viral proteins, proto-oncoproteins, with standard default parameters, the BRCA I zinc ring domain possesses 54% and regulatory and transcription factors (6). These include develop sequence similarity and 38% sequence identity with the EHV zinc ring domain, mental regulatory proteins [XNF-7, SINA, PSC, SU(z)2, and DG17}, compared with 42.8% sequence similarity and 30.6% sequence identity to the cell regulatory proteins (lAP and DROC3H), signal transduction PMLdomain.Fig. 2 shows the sequencealignmentbetween the BRCA1 and proteins (CRAF1), human oncoproteins (MDM2, MEL18, BMI-l, EHV and PML ring domains,respectively. RET, T18, CBL, and PML), proteins that regulate gene expression A BLASTP (Basic Local Alignment Search Tool for Proteins; Ref. 17) search of the Brookhaven Data Bank selected and aligned the EHV with the (RPT-l, PAS4, TLR, and HIP116), viral proteins, including some BRCA1 zinc ring domain, and yielded a P of 9.6 X l0@ (chance likelihood viral gene expression proteins, (lEl 10, VZ61, CG3O, PE-38, EPO, of the two sequences relationship). Folding search programs, such as Pro EHV, LCMV, and P11), and proteins involved in DNA repair or files-3D by Luthy et a!. (18), also indicated that the EHV and BRCA1 zinc ring recombination (RAD5, RADI 8, RAD16, RAG- 1, and UVS2). protein domains have similar folds. Hydrophilicity plots of Kyte and Doolittle The zinc ring family, first identified in 1991 (7), was found to (19), surface probability [calculated according to the method of Emini et a!. possess two finger-like domains connected by a linking region and to (20)], chain flexibility predictions (Karplus-Schulz), secondary structure pre require zinc for folding. The zinc ring C3HC4 sequence folds inde dictions of turns and a helices by Chou and Fasman (21), turns prediction by pendently and is a distinct stable domain (8, 9). This domain differs Gamier et al. (22), and the antigenicity index of Jameson and Wolf (23) all structurally from other zinc finger domains in that the first and third support greater similarity of the EHV ring domain structure to BRCAI than to pairs of cysteines form one zinc-binding site, and the second and PML(datanot shown). A multiple-sequencealignmentof the threesequences (the Pileup algorithm from the package below; Refs. 24 and 25) yields pairwise fourth pairs of cysteines form the second metal-binding site. This similarity scores ofO.44 for BRCA1 and EHV, 0.37 for BRCAI and PML, and crossed pattern is similar in configuration to a 1—3, 2—4 disulfide 0.42 for EHV and PML. Therefore, EHV was used as the basis for developing bond pattern. A schematic of the BRCA1 amino-terminal zinc ring a model for the BRCA1 zinc ring domain. domain with cysteine and histidines ligated to zinc indicated is shown Several multiple sequence alignments have been published previously, aligning the BRCAI zinc ring domain sequence to the sequences of other Received I2/28/95; accepted 4/3/96. known ring domains. The alignment by Borden et a!. (11) also shows greater The costs of publication of this article were defrayed in part by the payment of page sequence similarity between BRCA1 and EHV than between BRCA1 and charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. I To whom requests for reprints should be addressed, at National Institute of Envi 2 The abbreviations used are: EHV, equine herpes virus; PML, promyelocytic leuke ronmental Health Sciences, P.O. Box 12233, Mail Drop 10-03, Research Triangle Park, mia; NMR, nuclear magnetic resonance; CVFF, consistent valence force field; kcal, NC 27709. Phone: (919) 541-3397; Fax: (919) 541-1578; E-mail: rachelle@ kilocalorie; RMS, root-mean-square; HSV, herpes simplex virus; CRD, cysteine rich picard.niehs.nih.gov. domain of protein kinate C; WT1,Wilms' tumor gene; TFIIIA, transcription factor III A. 2539 Downloaded from cancerres.aacrjournals.org on September 30, 2021. © 1996 American Association for Cancer Research. MOLECULAR MODELING OF BRCAI AMINO-TERMINALZINC RING DOMAIN The two zinc ions were entered into the structure by superposition with the herpes virus structure, which gave a general positioning for the zinc ions in the protein. The resulting structure was energy minimized in vacuo using a distance-dependent dielectric constant, as described by Guenot and Kollman (28). The ligand-zinc ion distances were constrained in the minimization to the distances found between the corresponding ligands and zinc in the herpes virus protein. The zinc-cysteine sulfur atoms were pulled together using 500 cycles of conjugate gradients with a force of 100 kcal/mol, followed by minimization with 500 cycles of steepest descents and 1000 cycles of conjugate gradients until the maximum derivative was less than 0.1 kcal/A. All constraints were R then removed, and the structure was reminimized with 1000 cycles of conju F SQ L gate gradients until the maximum derivative was less than 0. 1 kcal/A. Charges were assigned to the zinc ions and cysteine ligands Sy and histidine Nbl . Because of the ionic and covalent nature of the zinc ligand interactions, the zinc ion charges are reduced and assigned values of + I.5 in the four cysteine-ligand environment and + 1.9 in the three-cysteine- and one-histidine ligand environment. All cysteine Sys have negative partial charges, because a hydrogen is abstracted on binding with zinc. The histidine N81 also has a negative partial charge on binding with zinc. In this model, the cysteine Sys in L the four-ligand case were assigned partial charges of —0.275,andcharges of I —0.450wereassigned in the three-ligand case. The histidine N bound to the MQK zinc was assigned a partial charge of —0.670. Distances were constrained in a several-step minimization (in vacuo) until Fig. 1. Schematic of the BRCA1 zinc ring domain with cysteine and histidines ligated the ligand-zinc distance criteria were met.