Importance of Incorporating Protein Flexibility in Molecule Modeling: A Theoretical Study on Type I1/2 NIK Inhibitors

NF-κB inducing kinase (NIK), which is considered as the central component of the non-canonical NF-κB pathway, has been proved to be an important target for the regulation of the immune system. In the past few years, NIK inhibitors with various scaffolds have been successively reported, among which type I1/2 inhibitors that can not only bind in the ATP-binding pocket at the DFG-in state but also extend into an additional back pocket, make up the largest proportion of the NIK inhibitors, and are worthy of more attention. In this study, an integration protocol that combines molecule docking, MD simulations, ensemble docking, MM/GB(PB)SA binding free energy calculations, and decomposition was employed to understand the binding mechanism of 21 tricyclic type I1/2 NIK inhibitors. It is found that the docking accuracy is largely dependent on the selection of docking protocols as well as the crystal structures. The predictions given by the ensemble docking based on multiple receptor conformations (MRCs) and the MM/GB(PB)SA calculations based on MD simulations showed higher linear correlations with the experimental data than those given by conventional rigid receptor docking (RRD) methods (Glide, GOLD, and Autodock Vina), highlighting the importance of incorporating protein flexibility in predicting protein–ligand interactions. Further analysis based on MM/GBSA demonstrates that the hydrophobic interactions play the most essential role in the ligand binding to NIK, and the polar interactions also make an important contribution to the NIK-ligand recognition. A deeper comparison of several pairs of representative derivatives reveals that the hydrophobic interactions are vitally important in the structural optimization of analogs as well. Besides, the H-bond interactions with some key residues and the large desolvation effect in the back pocket devote to the affinity distinction. It is expected that our study could provide valuable insights into the design of novel and potent type I1/2 NIK inhibitors.

[1]  J. Nagel,et al.  Assessing Performance , 2019, The United States and Canada.

[2]  Jianwen A. Feng,et al.  Scaffold-Hopping Approach To Discover Potent, Selective, and Efficacious Inhibitors of NF-κB Inducing Kinase. , 2018, Journal of medicinal chemistry.

[3]  M. Kubbutat,et al.  N-Acetyl-3-aminopyrazoles block the non-canonical NF-kB cascade by selectively inhibiting NIK. , 2018, MedChemComm.

[4]  A. Baldwin,et al.  Noncanonical NF-κB in Cancer , 2018, Biomedicines.

[5]  J. H. Zhang,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 7. Entropy effects on the performance of end-point binding free energy calculation approaches. , 2018, Physical chemistry chemical physics : PCCP.

[6]  Danfeng Shi,et al.  Understanding the structural and energetic basis of PD-1 and monoclonal antibodies bound to PD-L1: A molecular modeling perspective. , 2018, Biochimica et biophysica acta. General subjects.

[7]  Yuzong Chen,et al.  Computational identification of the binding mechanism of a triple reuptake inhibitor amitifadine for the treatment of major depressive disorder. , 2018, Physical chemistry chemical physics : PCCP.

[8]  Dan Li,et al.  Importance of protein flexibility on molecular recognition: modeling binding mechanisms of aminopyrazine inhibitors to Nek2. , 2018, Physical chemistry chemical physics : PCCP.

[9]  Tingjun Hou,et al.  Importance of protein flexibility in molecular recognition: a case study on Type-I1/2 inhibitors of ALK. , 2018, Physical chemistry chemical physics : PCCP.

[10]  M. Karin,et al.  NF-κB, inflammation, immunity and cancer: coming of age , 2018, Nature Reviews Immunology.

[11]  Adam R. Johnson,et al.  NF-κB inducing kinase is a therapeutic target for systemic lupus erythematosus , 2018, Nature Communications.

[12]  Qingchuan Zheng,et al.  Theoretical research in structure characteristics of different inhibitors and differences of binding modes with CBP bromodomain. , 2017, Bioorganic & medicinal chemistry.

[13]  Dan Li,et al.  Combating Drug-Resistant Mutants of Anaplastic Lymphoma Kinase with Potent and Selective Type-I1/2 Inhibitors by Stabilizing Unique DFG-Shifted Loop Conformation , 2017, ACS central science.

[14]  Feng Zhu,et al.  Revealing vilazodone's binding mechanism underlying its partial agonism to the 5-HT1A receptor in the treatment of major depressive disorder. , 2017, Physical chemistry chemical physics : PCCP.

[15]  Youyong Li,et al.  How Does the L884P Mutation Confer Resistance to Type-II Inhibitors of JAK2 Kinase: A Comprehensive Molecular Modeling Study , 2017, Scientific Reports.

[16]  Robert B. Kargbo.  New Substituted Cyanoindoline Derivatives as MAP3K14 Kinase Inhibitors for the Treatment of Cancer and Autoimmune Disorders. , 2017, ACS medicinal chemistry letters.

[17]  Shao-Cong Sun.  The non-canonical NF-κB pathway in immunity and inflammation , 2017, Nature Reviews Immunology.

[18]  Youyong Li,et al.  Discovery of Novel and Selective Adenosine A2A Receptor Antagonists for Treating Parkinson's Disease through Comparative Structure-Based Virtual Screening , 2017, J. Chem. Inf. Model..

[19]  Shan Chang,et al.  Molecular Simulation Studies on the Binding Selectivity of Type-I Inhibitors in the Complexes with ROS1 versus ALK , 2017, J. Chem. Inf. Model..

[20]  Michelle L Coote,et al.  Molecular dynamics-driven drug discovery: leaping forward with confidence. , 2017, Drug discovery today.

[21]  Adam R. Johnson,et al.  Structure-Based Design of Tricyclic NF-κB Inducing Kinase (NIK) Inhibitors That Have High Selectivity over Phosphoinositide-3-kinase (PI3K). , 2017, Journal of medicinal chemistry.

[22]  Qing-Chuan Zheng,et al.  Exploring the mechanism how AF9 recognizes and binds H3K9ac by molecular dynamics simulations and free energy calculations , 2016, Biopolymers.

[23]  Youyong Li,et al.  Comprehensive evaluation of ten docking programs on a diverse set of protein-ligand complexes: the prediction accuracy of sampling power and scoring power. , 2016, Physical chemistry chemical physics : PCCP.

[24]  Vinay Tergaonkar,et al.  Noncanonical NF-κB Signaling in Health and Disease. , 2016, Trends in molecular medicine.

[25]  A. Cavalli,et al.  Role of Molecular Dynamics and Related Methods in Drug Discovery. , 2016, Journal of medicinal chemistry.

[26]  Dan Li,et al.  Binding mechanisms of 1,4-dihydropyridine derivatives to L-type calcium channel Cav1.2: a molecular modeling study. , 2016, Molecular bioSystems.

[27]  Dan Li,et al.  Molecular principle of the cyclin-dependent kinase selectivity of 4-(thiazol-5-yl)-2-(phenylamino) pyrimidine-5-carbonitrile derivatives revealed by molecular modeling studies. , 2016, Physical chemistry chemical physics : PCCP.

[28]  Dan Li,et al.  Discovery of Novel ROCK1 Inhibitors via Integrated Virtual Screening Strategy and Bioassays , 2015, Scientific Reports.

[29]  S. Knapp,et al.  The ins and outs of selective kinase inhibitor development. , 2015, Nature chemical biology.

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

[31]  Dan Li,et al.  Importance of protein flexibility in ranking inhibitor affinities: modeling the binding mechanisms of piperidine carboxamides as Type I1/2 ALK inhibitors. , 2015, Physical chemistry chemical physics : PCCP.

[32]  Tingjun Hou,et al.  Assessing an Ensemble Docking-Based Virtual Screening Strategy for Kinase Targets by Considering Protein Flexibility , 2014, J. Chem. Inf. Model..

[33]  Tiziano Tuccinardi,et al.  Extensive Consensus Docking Evaluation for Ligand Pose Prediction and Virtual Screening Studies , 2014, J. Chem. Inf. Model..

[34]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. , 2014, Physical chemistry chemical physics : PCCP.

[35]  Tingjun Hou,et al.  P-loop Conformation Governed Crizotinib Resistance in G2032R-Mutated ROS1 Tyrosine Kinase: Clues from Free Energy Landscape , 2014, PLoS Comput. Biol..

[36]  Diane Joseph-McCarthy,et al.  Ensemble-Based Docking Using Biased Molecular Dynamics , 2014, J. Chem. Inf. Model..

[37]  Youyong Li,et al.  Exploring the prominent performance of CX-4945 derivatives as protein kinase CK2 inhibitors by a combined computational study. , 2014, Molecular bioSystems.

[38]  M. Shlomchik,et al.  Context-specific BAFF-R signaling by the NF-κB and PI3K pathways. , 2013, Cell reports.

[39]  Lin Li,et al.  Insights into susceptibility of antiviral drugs against the E119G mutant of 2009 influenza A (H1N1) neuraminidase by molecular dynamics simulations and free energy calculations. , 2013, Antiviral research.

[40]  Tingjun Hou,et al.  Development and Evaluation of an Integrated Virtual Screening Strategy by Combining Molecular Docking and Pharmacophore Searching Based on Multiple Protein Structures , 2013, J. Chem. Inf. Model..

[41]  Tingjun Hou,et al.  Insight into Crizotinib Resistance Mechanisms Caused by Three Mutations in ALK Tyrosine Kinase using Free Energy Calculation Approaches , 2013, J. Chem. Inf. Model..

[42]  Youyong Li,et al.  Theoretical study on the interaction of pyrrolopyrimidine derivatives as LIMK2 inhibitors: insight into structure-based inhibitor design. , 2013, Molecular bioSystems.

[43]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 3. The impact of force fields and ligand charge models. , 2013, The journal of physical chemistry. B.

[44]  Qian Zhang,et al.  Modeling Compound-Target Interaction Network of Traditional Chinese Medicines for Type II Diabetes Mellitus: Insight for Polypharmacology and Drug Design , 2013, J. Chem. Inf. Model..

[45]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[46]  Woody Sherman,et al.  Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments , 2013, Journal of Computer-Aided Molecular Design.

[47]  O. Plotnikova,et al.  Inhibiting NF-κB-inducing kinase (NIK): discovery, structure-based design, synthesis, structure-activity relationship, and co-crystal structures. , 2013, Bioorganic & medicinal chemistry letters.

[48]  Youyong Li,et al.  Discovery and optimization of triazine derivatives as ROCK1 inhibitors: molecular docking, molecular dynamics simulations and free energy calculations. , 2013, Molecular bioSystems.

[49]  Adam R. Johnson,et al.  The crystal structure of the catalytic domain of the NF-κB inducing kinase reveals a narrow but flexible active site. , 2012, Structure.

[50]  H. Zeh,et al.  NF-κB hyperactivation in tumor tissues allows tumor-selective reprogramming of the chemokine microenvironment to enhance the recruitment of cytolytic T effector cells. , 2012, Cancer research.

[51]  Huanxiang Liu,et al.  Understanding the molecular basis of MK2-p38α signaling complex assembly: insights into protein-protein interaction by molecular dynamics and free energy studies. , 2012, Molecular bioSystems.

[52]  Wei Wang,et al.  Characterization of domain-peptide interaction interface: prediction of SH3 domain-mediated protein-protein interaction network in yeast by generic structure-based models. , 2012, Journal of proteome research.

[53]  Roman A. Laskowski,et al.  LigPlot+: Multiple Ligand-Protein Interaction Diagrams for Drug Discovery , 2011, J. Chem. Inf. Model..

[54]  Tingjun Hou,et al.  Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking , 2011, J. Comput. Chem..

[55]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[56]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[57]  David Calkins,et al.  Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution , 2010, J. Comput. Aided Mol. Des..

[58]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[59]  A. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[60]  David C. Thompson,et al.  Comparison of Several Molecular Docking Programs: Pose Prediction and Virtual Screening Accuracy , 2009, J. Chem. Inf. Model..

[61]  M. Karin,et al.  Regulation and function of NF-kappaB transcription factors in the immune system. , 2009, Annual review of immunology.

[62]  Wei Zhang,et al.  Predicting drug resistance of the HIV‐1 protease using molecular interaction energy components , 2009, Proteins.

[63]  J. Keats,et al.  Nonredundant and complementary functions of TRAF2 and TRAF3 in a ubiquitination cascade that activates NIK-dependent alternative NF-κB signaling , 2008, Nature Immunology.

[64]  T. Mak,et al.  Activation of noncanonical NF-κB requires coordinated assembly of a regulatory complex of the adaptors cIAP1, cIAP2, TRAF2, TRAF3 and the kinase NIK , 2008, Nature Immunology.

[65]  William A. McLaughlin,et al.  Evaluating the potency of HIV‐1 protease drugs to combat resistance , 2008, Proteins.

[66]  R. Abagyan,et al.  Flexible ligand docking to multiple receptor conformations: a practical alternative. , 2008, Current opinion in structural biology.

[67]  Heather A Carlson,et al.  Small molecule inhibitors of the MDM2-p53 interaction discovered by ensemble-based receptor models. , 2007, Journal of the American Chemical Society.

[68]  Jeremy R. Greenwood,et al.  Epik: a software program for pKa prediction and protonation state generation for drug-like molecules , 2007, J. Comput. Aided Mol. Des..

[69]  Tingjun Hou,et al.  Molecular dynamics and free energy studies on the wild-type and double mutant HIV-1 protease complexed with amprenavir and two amprenavir-related inhibitors: mechanism for binding and drug resistance. , 2007, Journal of medicinal chemistry.

[70]  J. Gready,et al.  Combining docking and molecular dynamic simulations in drug design , 2006, Medicinal research reviews.

[71]  Xiaojie Xu,et al.  Recent Advances in Free Energy Calculations with a Combination of Molecular Mechanics and Continuum Models , 2006 .

[72]  Ray Luo,et al.  How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis. , 2006, The journal of physical chemistry. B.

[73]  R. Friesner,et al.  Novel procedure for modeling ligand/receptor induced fit effects. , 2006, Journal of medicinal chemistry.

[74]  Michael Karin,et al.  NF-κB: linking inflammation and immunity to cancer development and progression , 2005, Nature Reviews Immunology.

[75]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

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

[77]  M. Karin,et al.  The two NF-κB activation pathways and their role in innate and adaptive immunity , 2004 .

[78]  Yong Duan,et al.  Distinguish protein decoys by Using a scoring function based on a new AMBER force field, short molecular dynamics simulations, and the generalized born solvent model , 2004, Proteins.

[79]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[80]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[81]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[82]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[83]  J. Caamaño,et al.  Regulation of secondary lymphoid organ development by the nuclear factor‐κB signal transduction pathway , 2003, Immunological reviews.

[84]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[85]  Michael Karin,et al.  Activation by IKKα of a Second, Evolutionary Conserved, NF-κB Signaling Pathway , 2001, Science.

[86]  E. Harhaj,et al.  NF-κB-Inducing Kinase Regulates the Processing of NF-κB2 p100 , 2001 .

[87]  T. Darden,et al.  Efficient particle-mesh Ewald based approach to fixed and induced dipolar interactions , 2000 .

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

[89]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[90]  D. Wallach,et al.  MAP3K-related kinase involved in NF-KB induction by TNF, CD95 and IL-1 , 1997, Nature.

[91]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

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

[93]  Tingjun Hou,et al.  Structural Bioinformatics Prediction of Peptides Binding to the Pka Riiα Subunit Using a Hierarchical Strategy , 2022 .