Charge Profile Analysis Reveals That Activation of Pro-apoptotic Regulators Bax and Bak Relies on Charge Transfer Mediated Allosteric Regulation

The pro-apoptotic proteins Bax and Bak are essential for executing programmed cell death (apoptosis), yet the mechanism of their activation is not properly understood at the structural level. For the first time in cell death research, we calculated intra-protein charge transfer in order to study the structural alterations and their functional consequences during Bax activation. Using an electronegativity equalization model, we investigated the changes in the Bax charge profile upon activation by a functional peptide of its natural activator protein, Bim. We found that charge reorganizations upon activator binding mediate the exposure of the functional sites of Bax, rendering Bax active. The affinity of the Bax C-domain for its binding groove is decreased due to the Arg94-mediated abrogation of the Ser184-Asp98 interaction. We further identified a network of charge reorganizations that confirms previous speculations of allosteric sensing, whereby the activation information is conveyed from the activation site, through the hydrophobic core of Bax, to the well-distanced functional sites of Bax. The network was mediated by a hub of three residues on helix 5 of the hydrophobic core of Bax. Sequence and structural alignment revealed that this hub was conserved in the Bak amino acid sequence, and in the 3D structure of folded Bak. Our results suggest that allostery mediated by charge transfer is responsible for the activation of both Bax and Bak, and that this might be a prototypical mechanism for a fast activation of proteins during signal transduction. Our method can be applied to any protein or protein complex in order to map the progress of allosteric changes through the proteins' structure.

[1]  Xu Luo,et al.  A three-helix homo-oligomerization domain containing BH3 and BH1 is responsible for the apoptotic activity of Bax. , 2007, Genes & development.

[2]  Y. Hsu,et al.  Conformation of the Bax C‐terminus regulates subcellular location and cell death , 1999, The EMBO journal.

[3]  Y. Pang,et al.  Transient binding of an activator BH3 domain to the Bak BH3-binding groove initiates Bak oligomerization , 2011, The Journal of cell biology.

[4]  T. Kuwana,et al.  Bcl-2-family proteins and the role of mitochondria in apoptosis. , 2003, Current opinion in cell biology.

[5]  Yizeng Liang,et al.  A modified electronegativity equalization method for fast and accurate calculation of atomic charges in large biological molecules. , 2009, Physical chemistry chemical physics : PCCP.

[6]  C. Cavasotto,et al.  Quantum-Mechanical Molecular Dynamics of Charge Transfer , 2010 .

[7]  S. Korsmeyer,et al.  VDAC2 Inhibits BAK Activation and Mitochondrial Apoptosis , 2003, Science.

[8]  L. G. Dias,et al.  Calculation of the dipole moment for polypeptides using the generalized born-electronegativity equalization method: Results in vacuum and continuum-dielectric solvent , 2004 .

[9]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[10]  Jaroslav Koca,et al.  TRITON: a graphical tool for ligand-binding protein engineering , 2008, Bioinform..

[11]  D. Vaux,et al.  Inhibition of Bak Activation by VDAC2 Is Dependent on the Bak Transmembrane Anchor* , 2010, The Journal of Biological Chemistry.

[12]  Chang-Sheng Wang,et al.  Atom−Bond Electronegativity Equalization Method. 1. Calculation of the Charge Distribution in Large Molecules , 1997 .

[13]  Nico Tjandra,et al.  Structure of Bax Coregulation of Dimer Formation and Intracellular Localization , 2000, Cell.

[14]  Patrick Bultinck,et al.  Toward an Alternative Hardness Kernel Matrix Structure in the Electronegativity Equalization Method (EEM) , 2006, J. Chem. Inf. Model..

[15]  Y. Hsu,et al.  Cytosol-to-membrane redistribution of Bax and Bcl-X(L) during apoptosis. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Charles L. Brooks,et al.  Are Many-Body Effects Important in Protein Folding? , 2000 .

[17]  H. Scheraga,et al.  An efficient method for calculating atomic charges of peptides and proteins from electronic populations. , 2008, The journal of physical chemistry. B.

[18]  Mahlon Wilson,et al.  Comparison between the geometric and harmonic mean electronegativity equilibration techniques , 1989 .

[19]  Victor Guallar,et al.  Importance of accurate charges in molecular docking: Quantum mechanical/molecular mechanical (QM/MM) approach , 2005, J. Comput. Chem..

[20]  L. Pauling.  The Nature Of The Chemical Bond , 1939 .

[21]  D. Green,et al.  Mitochondria and cell death: outer membrane permeabilization and beyond , 2010, Nature Reviews Molecular Cell Biology.

[22]  R. Nussinov,et al.  Allostery: absence of a change in shape does not imply that allostery is not at play. , 2008, Journal of molecular biology.

[23]  W. Craigen,et al.  VDAC2 is required for truncated BID‐induced mitochondrial apoptosis by recruiting BAK to the mitochondria , 2009, EMBO reports.

[24]  Zimei Bu,et al.  Proteins move! Protein dynamics and long-range allostery in cell signaling. , 2011, Advances in protein chemistry and structural biology.

[25]  Eszter Czinki,et al.  A combined electronegativity equalization and electrostatic potential fit method for the determination of atomic point charges , 2007, J. Comput. Chem..

[26]  R. Parr,et al.  Electronegativity: The density functional viewpoint , 1978 .

[27]  D. Andrews,et al.  Membrane Binding by tBid Initiates an Ordered Series of Events Culminating in Membrane Permeabilization by Bax , 2008, Cell.

[28]  K. Smirnov.  Computer Modeling Study of Interaction of Acetonitrile with Hydroxyl Groups of HY Zeolite , 2001 .

[29]  P. Cramer,et al.  CPD Damage Recognition by Transcribing RNA Polymerase II , 2007, Science.

[30]  S. Korsmeyer,et al.  A stapled BID BH3 helix directly binds and activates BAX. , 2006, Molecular cell.

[31]  Yi-Te Hsu,et al.  Movement of Bax from the Cytosol to Mitochondria during Apoptosis , 1997, The Journal of cell biology.

[32]  D. Green,et al.  How do BCL-2 proteins induce mitochondrial outer membrane permeabilization? , 2008, Trends in cell biology.

[33]  P. Colman,et al.  To trigger apoptosis, Bak exposes its BH3 domain and homodimerizes via BH3:groove interactions. , 2008, Molecular cell.

[34]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[35]  X. Bao,et al.  Investigation of the regio- and stereoselectivity of Diels–Alder reactions by newly developed ABEEM σπ model on the basis of local HSAB principle and maximum hardness principle , 2002 .

[36]  Zhong-Zhi Yang,et al.  General atom-bond electronegativity equalization method and its application in prediction of charge distributions in polypeptide , 2000 .

[37]  Simon Suchomel,et al.  Electronegativity Equalization Method: Parameterization and Validation for Large Sets of Organic, Organohalogene and Organometal Molecule , 2007, International Journal of Molecular Sciences.

[38]  Erinna F. Lee,et al.  Apoptosis Initiated When BH3 Ligands Engage Multiple Bcl-2 Homologs, Not Bax or Bak , 2007, Science.

[39]  S. Korsmeyer,et al.  Distinct BH3 domains either sensitize or activate mitochondrial apoptosis, serving as prototype cancer therapeutics. , 2002, Cancer cell.

[40]  P. Colman,et al.  Bax activation by Bim? , 2009, Cell Death and Differentiation.

[41]  G. Dewson,et al.  Molecular biology of Bax and Bak activation and action. , 2011, Biochimica et biophysica acta.

[42]  Rodrigo Lopez,et al.  Clustal W and Clustal X version 2.0 , 2007, Bioinform..

[43]  Ursula Rothlisberger,et al.  Polarization effects and charge transfer in the KcsA potassium channel. , 2006, Biophysical chemistry.

[44]  Mark S. Johnson,et al.  Accurate conformation‐dependent molecular electrostatic potentials for high‐throughput in silico drug discovery , 2009, J. Comput. Chem..

[45]  R. T. Sanderson,et al.  An Interpretation of Bond Lengths and a Classification of Bonds. , 1951, Science.

[46]  T. Kuwana,et al.  BH3 domains of BH3-only proteins differentially regulate Bax-mediated mitochondrial membrane permeabilization both directly and indirectly. , 2005, Molecular cell.

[47]  Jens Carlsson,et al.  Charges for Large Scale Binding Free Energy Calculations with the Linear Interaction Energy Method. , 2009, Journal of chemical theory and computation.

[48]  J. Martinou,et al.  Bid Induces the Oligomerization and Insertion of Bax into the Outer Mitochondrial Membrane , 2000, Molecular and Cellular Biology.

[49]  Ralph G. Pearson,et al.  Absolute hardness: companion parameter to absolute electronegativity , 1983 .

[50]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[51]  N. Tjandra,et al.  BAX Activation is Initiated at a Novel Interaction Site , 2008, Nature.

[52]  Linus Pauling,et al.  THE NATURE OF THE CHEMICAL BOND. IV. THE ENERGY OF SINGLE BONDS AND THE RELATIVE ELECTRONEGATIVITY OF ATOMS , 1932 .

[53]  M. Rehm,et al.  Single-cell quantification of Bax activation and mathematical modelling suggest pore formation on minimal mitochondrial Bax accumulation , 2010, Cell Death and Differentiation.

[54]  L. Walensky,et al.  BH3-triggered structural reorganization drives the activation of proapoptotic BAX. , 2010, Molecular cell.

[55]  Kevin A. Roth,et al.  Dysfunction and Death Proapoptotic BAX and BAK : A Requisite Gateway to Mitochondrial , 2014 .

[56]  Michel Waroquier,et al.  The Electronegativity Equalization Method I: Parametrization and Validation for Atomic Charge Calculations , 2002 .

[57]  A. Strasser,et al.  Control of apoptosis in the immune system: Bcl-2, BH3-only proteins and more. , 2003, Annual review of immunology.

[58]  L. G. Dias,et al.  Parameterization of the electronegativity equalization method based on the charge model 1 , 2002 .

[59]  Toon Verstraelen,et al.  The electronegativity equalization method and the split charge equilibration applied to organic systems: parametrization, validation, and comparison. , 2009, The Journal of chemical physics.

[60]  Lorenzo Galluzzi,et al.  Mitochondrial membrane permeabilization in cell death. , 2007, Physiological reviews.

[61]  David W. Andrews,et al.  Embedded together: The life and death consequences of interaction of the Bcl-2 family with membranes , 2007, Apoptosis.

[62]  S. Korsmeyer,et al.  Bcl-2 heterodimerizes in vivo with a conserved homolog, Bax, that accelerates programed cell death , 1993, Cell.

[63]  Radka Svobodová Vareková,et al.  Optimized and parallelized implementation of the electronegativity equalization method and the atom‐bond electronegativity equalization method , 2006, J. Comput. Chem..

[64]  A. Allred,et al.  Electronegativity values from thermochemical data , 1961 .

[65]  A. Tocilj,et al.  The X-ray structure of a BAK homodimer reveals an inhibitory zinc binding site. , 2006, Molecular cell.

[66]  R. Landick,et al.  Allosteric Control of RNA Polymerase by a Site That Contacts Nascent RNA Hairpins , 2001, Science.