Diverse dynamics features of novel protein kinase C (PKC) isozymes determine the selectivity of a fluorinated balanol analogue for PKCε

Background(−)-Balanol is an ATP-mimicking inhibitor that non-selectively targets protein kinase C (PKC) isozymes and cAMP-dependent protein kinase (PKA). While PKA constantly shows tumor promoting activities, PKC isozymes can ambiguously be tumor promoters or suppressors. In particular, PKCε is frequently implicated in tumorigenesis and a potential target for anticancer drugs. We recently reported that the C5(S)-fluorinated balanol analogue (balanoid 1c) had improved binding affinity and selectivity for PKCε but not to the other novel PKC isozymes, which share a highly similar ATP site. The underlying basis for this fluorine-based selectivity is not entirely comprehended and needs to be investigated further for the development of ATP mimic inhibitors specific for PKCε.ResultsUsing molecular dynamics (MD) simulations assisted by homology modelling and sequence analysis, we have studied the fluorine-based selectivity in the highly similar ATP sites of novel PKC (nPKC) isozymes. The study suggests that every nPKC isozyme has different dynamics behaviour in both apo and 1c-bound forms. Interestingly, the apo form of PKCε, where 1c binds strongly, shows the highest degree of flexibility which dramatically decreases after binding 1c.ConclusionsFor the first time to the best of our knowledge, we found that the origin of 1c selectivity for PKCε comes from the unique dynamics feature of each PKC isozyme. Fluorine conformational control in 1c can synergize with and lock down the dynamics of PKCε, which optimize binding interactions with the ATP site residues of the enzyme, particularly the invariant Lys437. This finding has implications for further rational design of balanol-based PKCε inhibitors for cancer drug development.

[1]  K. Das,et al.  Protein kinase C, an elusive therapeutic target? , 2012, Nature Reviews Drug Discovery.

[2]  G. E. Jagdmann,et al.  Increasing the cellular PKC inhibitory activity of balanol: A study of ester analogs , 1995 .

[3]  L. Hunter,et al.  Stereoselectively fluorinated N-heterocycles: a brief survey , 2013, Beilstein journal of organic chemistry.

[4]  Ledion Bitincka,et al.  MatGAT: An application that generates similarity/identity matrices using protein or DNA sequences , 2003, BMC Bioinformatics.

[5]  Ari Hardianto,et al.  Divergent response of homologous ATP sites to stereospecific ligand fluorination for selectivity enhancement. , 2017, Organic & biomolecular chemistry.

[6]  Luigi Atzori,et al.  Structural and Dynamical Insights on HLA-DR2 Complexes That Confer Susceptibility to Multiple Sclerosis in Sardinia: A Molecular Dynamics Simulation Study , 2013, PloS one.

[7]  Ying Liu,et al.  Evol and ProDy for bridging protein sequence evolution and structural dynamics , 2014, Bioinform..

[8]  Jörg Weiser,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999, J. Comput. Chem..

[9]  Cedric E. Ginestet.  ggplot2: Elegant Graphics for Data Analysis , 2011 .

[10]  Rodrigo Lopez,et al.  Analysis Tool Web Services from the EMBL-EBI , 2013, Nucleic Acids Res..

[11]  Fei Liu,et al.  Exploration of charge states of balanol analogues acting as ATP-competitive inhibitors in kinases , 2017, BMC Bioinformatics.

[12]  Anselm H. C. Horn,et al.  AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine , 2006, Journal of molecular modeling.

[13]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[14]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[15]  M. Bunnage,et al.  Total Synthesis of Balanol and Designed Analogues. , 1996 .

[16]  Ari Hardianto,et al.  Molecular Dynamics Pinpoint the Global Fluorine Effect in Balanoid Binding to PKCε and PKA , 2018, J. Chem. Inf. Model..

[17]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[18]  W. Janzen,et al.  Synthesis and protein kinase C inhibitory activities of balanol analogues with modification of 4-hydroxybenzamido moiety. , 1997, Bioorganic & medicinal chemistry.

[19]  Zaida Luthey-Schulten,et al.  MultiSeq: unifying sequence and structure data for evolutionary analysis , 2006, BMC Bioinformatics.

[20]  M. Kazanietz,et al.  Protein kinase C and cancer: what we know and what we do not , 2014, Oncogene.

[21]  Jacob Kongsted,et al.  Binding free energy based structural dynamics analysis of HIV-1 RT RNase H-inhibitor complexes. , 2014, Integrative biology : quantitative biosciences from nano to macro.

[22]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[23]  Lars V. Schäfer,et al.  Molecular mechanism of peptide editing in the tapasin–MHC I complex , 2016, Scientific Reports.

[24]  Josep Ramón Goñi,et al.  Molecular dynamics simulations: advances and applications , 2015, Advances and applications in bioinformatics and chemistry : AABC.

[25]  M. G. Johnson,et al.  Synthesis and protein kinase inhibitory activity of balanol analogues with modified benzophenone subunits. , 2002, Journal of medicinal chemistry.

[26]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[27]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[28]  M. D. Hill,et al.  Applications of Fluorine in Medicinal Chemistry. , 2015, Journal of medicinal chemistry.

[29]  Jennifer L. Knight,et al.  OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. , 2016, Journal of chemical theory and computation.

[30]  The Uniprot Consortium,et al.  UniProt: a hub for protein information , 2014, Nucleic Acids Res..

[31]  Narayanan Eswar,et al.  Protein structure modeling with MODELLER. , 2008, Methods in molecular biology.

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

[33]  B. Hess,et al.  Similarities between principal components of protein dynamics and random diffusion , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[35]  R. Mannella,et al.  Langevin stabilization of molecular-dynamics simulations of polymers by means of quasisymplectic algorithms. , 2007, The Journal of chemical physics.

[36]  L. Hunter,et al.  The C–F bond as a conformational tool in organic and biological chemistry , 2010, Beilstein journal of organic chemistry.

[37]  S. Genheden,et al.  The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities , 2015, Expert opinion on drug discovery.

[38]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[39]  Susan S. Taylor,et al.  Protein kinases: evolution of dynamic regulatory proteins. , 2011, Trends in biochemical sciences.

[40]  Matthias Stein,et al.  webPIPSA: a web server for the comparison of protein interaction properties , 2008, Nucleic Acids Res..

[41]  Janet M. Thornton,et al.  PROCHECK: validation of protein‐structure coordinates , 2012 .

[42]  Gerhard Wolber,et al.  The impact of molecular dynamics on drug design: applications for the characterization of ligand-macromolecule complexes. , 2015, Drug discovery today.

[43]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[44]  A. Sali,et al.  Statistical potential for assessment and prediction of protein structures , 2006, Protein science : a publication of the Protein Society.

[45]  W. Janzen,et al.  Balanol: A Novel and Potent Inhibitor of Protein Kinase C from the Fungus Verticillium balanoides. , 1993 .

[46]  María Martín,et al.  UniProt: A hub for protein information , 2015 .

[47]  Y. Cho‐Chung,et al.  Biochemical characterization of extracellular cAMP-dependent protein kinase as a tumor marker. , 2000, Biochemical and biophysical research communications.

[48]  J. Lampe,et al.  Protein kinase C inhibitory activities of balanol analogs bearing carboxylic acid replacements , 1995 .

[49]  Susan S. Taylor,et al.  PKA: a portrait of protein kinase dynamics. , 2004, Biochimica et biophysica acta.

[50]  K C Nicolaou,et al.  Molecular design and biological activity of potent and selective protein kinase inhibitors related to balanol. , 1995, Chemistry & biology.

[51]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[52]  Holger Gohlke,et al.  Converging free energy estimates: MM‐PB(GB)SA studies on the protein–protein complex Ras–Raf , 2004, J. Comput. Chem..