Cooperative binding of DNA and CBFβ to the Runt domain of the CBFα studied via MD simulations

The Runt domain (RD) is the DNA-binding region of the Runx genes. A related protein, known as core binding factor β (CBFβ) also binds to the RD to enhance RD–DNA interaction by 6- to 10-fold. Here, we report results from molecular dynamics (MD) simulations of RD alone, as a dimer in complexes with DNA and CBFβ and in a ternary complex with DNA and CBFβ. Consistent with the experimental findings, in the presence of CBFβ the estimated free energy of binding of RD to the DNA is more favorable, which is shown to be due to more favorable intermolecular interactions and desolvation contributions. Also contributing to the enhanced binding are favorable intramolecular interactions between the ‘wing’ residues (RD residues 139–145) and the ‘wing1’ residues (RD residues 104–116). The simulation studies also indicate that the RD–CBFβ binding is more favorable in the presence of DNA due to a more favorable RD–CBFβ interaction energy. In addition, it is predicted that long-range interactions involving ionic residues contribute to binding cooperativity. Results from the MD calculations are used to interpret a variety of experimental mutagenesis data. A novel role for RD Glu116 to the RD–CBFβ interaction is predicted.

[1]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[2]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[3]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[4]  D W Hukins,et al.  Refinement of the structure of B-DNA and implications for the analysis of x-ray diffraction data from fibers of biopolymers. , 1973, Journal of molecular biology.

[5]  R. Nussinov,et al.  Protein–protein interactions: Structurally conserved residues distinguish between binding sites and exposed protein surfaces , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Bernard R. Brooks,et al.  New spherical‐cutoff methods for long‐range forces in macromolecular simulation , 1994, J. Comput. Chem..

[7]  S. Itohara,et al.  Causal Relationship between the Loss of RUNX3 Expression and Gastric Cancer , 2002, Cell.

[8]  J. Bravo,et al.  Energetic Contribution of Residues in the Runx1 Runt Domain to DNA Binding* , 2003, Journal of Biological Chemistry.

[9]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[10]  J. Kelley,et al.  Biophysical characterization of interactions between the core binding factor alpha and beta subunits and DNA. , 2000, FEBS letters.

[11]  J. Bushweller,et al.  Biochemical and biophysical properties of the core-binding factor alpha2 (AML1) DNA-binding domain. , 1996, Journal of Biological Chemistry.

[12]  J. Bravo,et al.  Energetic and Functional Contribution of Residues in the Core Binding Factor β (CBFβ) Subunit to Heterodimerization with CBFα* , 2000, The Journal of Biological Chemistry.

[13]  S. Mundlos,et al.  Mutations in the RUNX2 gene in patients with cleidocranial dysplasia , 2002, Human mutation.

[14]  T. Tahirov,et al.  Crystallization and preliminary X-ray analyses of quaternary, ternary and binary protein–DNA complexes with involvement of AML1/Runx-1/CBFα Runt domain, CBFβ and the C/EBPβ bZip region , 2001 .

[15]  J. Bravo,et al.  Mutagenesis of the Runt Domain Defines Two Energetic Hot Spots for Heterodimerization with the Core Binding Factor β Subunit* , 2003, Journal of Biological Chemistry.

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

[17]  E A Merritt,et al.  Raster3D: photorealistic molecular graphics. , 1997, Methods in enzymology.

[18]  J. Downing,et al.  AML1, the Target of Multiple Chromosomal Translocations in Human Leukemia, Is Essential for Normal Fetal Liver Hematopoiesis , 1996, Cell.

[19]  Jerónimo Bravo,et al.  The leukemia-associated AML1 (Runx1)–CBFβ complex functions as a DNA-induced molecular clamp , 2001, Nature Structural Biology.

[20]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data , 2000 .

[21]  J. Kelley,et al.  Biophysical characterization of interactions between the core binding factor α and β subunits and DNA , 2000 .

[22]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[23]  Chaohong Sun,et al.  The Ig fold of the core binding factor alpha Runt domain is a member of a family of structurally and functionally related Ig-fold DNA-binding domains. , 1999, Structure.

[24]  D. Rabinovich,et al.  DNA recognition by the RUNX1 transcription factor is mediated by an allosteric transition in the RUNT domain and by DNA bending. , 2002, Structure.

[25]  T. Shimada,et al.  Identification of a novel frameshift mutation (383insT) in the RUNX2 (PEBP2 α/CBFA1/AML3) gene in a Japanese patient with cleidocranial dysplasia , 2001, Journal of Bone and Mineral Metabolism.

[26]  T. Motoyama,et al.  Hypermethylation of the RUNX3 gene promoter in testicular yolk sac tumor of infants. , 2003, The American journal of pathology.

[27]  U. Sauer,et al.  The RUNX1 Runt domain at 1.25A resolution: a structural switch and specifically bound chloride ions modulate DNA binding. , 2002, Journal of molecular biology.

[28]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution , 2000 .

[29]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[30]  Alexander D. MacKerell,et al.  Computational Biochemistry and Biophysics , 2001 .

[31]  A. Das,et al.  Free‐energy component analysis of 40 protein–DNA complexes: A consensus view on the thermodynamics of binding at the molecular level , 2002, J. Comput. Chem..

[32]  A. Lund,et al.  RUNX: a trilogy of cancer genes. , 2002, Cancer cell.

[33]  A. Passaniti,et al.  Runt-related gene 2 in endothelial cells: inducible expression and specific regulation of cell migration and invasion. , 2001, Cancer research.

[34]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[35]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[36]  Roger L. Williams,et al.  Structural basis for the heterodimeric interaction between the acute leukaemia‐associated transcription factors AML1 and CBFβ , 2000, The EMBO journal.

[37]  M. Abe,et al.  Indispensable role of the transcription factor PEBP2/CBF in angiogenic activity of a murine endothelial cell MSS31 , 2000, Oncogene.

[38]  Alexander D. MacKerell,et al.  All-atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution , 2000, J. Comput. Chem..

[39]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[40]  D. Foguel,et al.  DNA tightens the dimeric DNA-binding domain of human papillomavirus E2 protein without changes in volume. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[41]  C. Brooks,et al.  Novel generalized Born methods , 2002 .

[42]  J. Bushweller,et al.  CBFβ allosterically regulates the Runx1 Runt domain via a dynamic conformational equilibrium , 2004, Nature Structural &Molecular Biology.

[43]  J. Bushweller,et al.  Biochemical and Biophysical Properties of the Core-binding Factor α2 (AML1) DNA-binding Domain* , 1996, The Journal of Biological Chemistry.

[44]  B. Brooks,et al.  Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamide , 1992, Biopolymers.

[45]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[46]  Y. Ito,et al.  PEBP2/PEA2 represents a family of transcription factors homologous to the products of the Drosophila runt gene and the human AML1 gene. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[47]  Takashi Kumasaka,et al.  Structural Analyses of DNA Recognition by the AML1/Runx-1 Runt Domain and Its Allosteric Control by CBFβ , 2001, Cell.

[48]  William L. Jorgensen,et al.  Theoretical studies of medium effects on conformational equilibria , 1983 .

[49]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.