Human PD-1 binds differently to its human ligands: a comprehensive modeling study.

Programmed death-1 (PD-1) is a potent inhibitory receptor of T cells which binds to two different ligands, namely PD-L1 and PD-L2, and upon binding, it inhibits T cell activation, differentiation, and proliferation, leading to a state of immune tolerance. Blocking these interactions recently emerged as a 'game changer' approach in immunotherapy. Despite the significant therapeutic potential of targeting the PD-1 pathway, the interaction between human PD-1 and its two human ligands is not fully understood. Current crystal structures describe the interactions of mouse PD-1 with human PD-L1 or mouse PD-L2. However, recent mutational and nuclear magnetic resonance (NMR) analyses suggest that human PD-1 binds its human ligands differently compared to their mouse counterparts. No detailed model is currently available to consistently fit these data. The lack of these accurate structures constitutes a high barrier against rationally developing more effective and safer agents targeting these interactions. Here we describe for the first time two accurate models for human PD-1 bound to its two human ligands. Our methodology involved combining molecular dynamics (MD) simulations with protein-protein docking and binding energy analysis to predict the most probable binding conformations for PD1 to its ligands. Our results confirm the available experimental NMR and mutational data and reveal the most accurate atomistic details so far of how human PD-1 binds to human PD-Ls and why the two ligands bind with different affinities to the same receptor.

[1]  Jürgen Bajorath,et al.  Molecular Modeling and Functional Mapping of B7-H1 and B7-DC Uncouple Costimulatory Function from PD-1 Interaction , 2003, The Journal of experimental medicine.

[2]  G. Freeman,et al.  Programmed death-1 ligand 1 interacts specifically with the B7-1 costimulatory molecule to inhibit T cell responses. , 2007, Immunity.

[3]  Jack A. Tuszynski,et al.  Rational Drug Design: One Target, Many Paths to It , 2014, Int. J. Comput. Model. Algorithms Medicine.

[4]  B. Creelan Update on immune checkpoint inhibitors in lung cancer. , 2014, Cancer control : journal of the Moffitt Cancer Center.

[5]  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..

[6]  A. Eggermont,et al.  Drug of the year: programmed death-1 receptor/programmed death-1 ligand-1 receptor monoclonal antibodies. , 2013, European journal of cancer.

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

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

[9]  Z. Weng,et al.  A novel shape complementarity scoring function for protein‐protein docking , 2003, Proteins.

[10]  T. Okazaki,et al.  The PD-1-PD-L pathway in immunological tolerance. , 2006, Trends in immunology.

[11]  John E. Ladbury,et al.  Structure and Interactions of the Human Programmed Cell Death 1 Receptor* , 2013, The Journal of Biological Chemistry.

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

[13]  Donald W. Bouldin,et al.  A Cluster Separation Measure , 1979, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[14]  Samuel Genheden,et al.  How to obtain statistically converged MM/GBSA results , 2009, J. Comput. Chem..

[15]  M. Mandalà,et al.  Targeting the PD1/PD-L1 axis in melanoma: biological rationale, clinical challenges and opportunities. , 2014, Critical reviews in oncology/hematology.

[16]  Tingjun Hou,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. , 2014, Physical chemistry chemical physics : PCCP.

[17]  Kui Wang,et al.  Human structural proteome-wide characterization of Cyclosporine A targets , 2014, Bioinform..

[18]  M. Sternberg,et al.  Modelling protein docking using shape complementarity, electrostatics and biochemical information. , 1997, Journal of molecular biology.

[19]  C. Thompson,et al.  At the Bench: Preclinical rationale for CTLA‐4 and PD‐1 blockade as cancer immunotherapy , 2013, Journal of leukocyte biology.

[20]  Dan Li,et al.  Structural basis of the interactions between CXCR4 and CXCL12/SDF-1 revealed by theoretical approaches. , 2013, Molecular bioSystems.

[21]  Sheng-You Huang,et al.  Search strategies and evaluation in protein-protein docking: principles, advances and challenges. , 2014, Drug discovery today.

[22]  Jonathan Y. Mane,et al.  Discovery of Small Molecule Inhibitors that Interact with γ‐Tubulin , 2012, Chemical biology & drug design.

[23]  Lijuan Huang,et al.  The PD-1/PD-Ls pathway and autoimmune diseases. , 2014, Cellular immunology.

[24]  E. Katchalski‐Katzir,et al.  Molecular surface recognition: determination of geometric fit between proteins and their ligands by correlation techniques. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[25]  J. Tuszynski,et al.  Relaxed complex scheme suggests novel inhibitors for the lyase activity of DNA polymerase beta. , 2011, Journal of molecular graphics & modelling.

[26]  R. Madan,et al.  Nivolumab: promising survival signal coupled with limited toxicity raises expectations. , 2014, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[27]  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.

[28]  Jonathan Y. Mane,et al.  Modeling the Yew Tree Tubulin and a Comparison of its Interaction with Paclitaxel to Human Tubulin , 2012, Pharmaceutical Research.

[29]  K. E. Jordan,et al.  A computational model for overcoming drug resistance using selective dual-inhibitors for aurora kinase A and its T217D variant. , 2013, Molecular pharmaceutics.

[30]  Jack A. Tuszynski,et al.  Detailed Computational Study of the Active Site of the Hepatitis C Viral RNA Polymerase to Aid Novel Drug Design , 2013, J. Chem. Inf. Model..

[31]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[32]  Lieping Chen Co-inhibitory molecules of the B7–CD28 family in the control of T-cell immunity , 2004, Nature Reviews Immunology.

[33]  Jonathan Y. Mane,et al.  Ensemble-based virtual screening reveals dual-inhibitors for the p53-MDM2/MDMX interactions. , 2010, Journal of molecular graphics & modelling.

[34]  J. Tuszynski,et al.  Interactions of laulimalide, peloruside, and their derivatives with the isoforms of β-tubulin , 2013 .

[35]  Zhiping Weng,et al.  ZDOCK and RDOCK performance in CAPRI rounds 3, 4, and 5 , 2005, Proteins.

[36]  K. Barakat Do We Need Small Molecule Inhibitors for the Immune Checkpoints , 2014 .

[37]  Zhiping Weng,et al.  Docking unbound proteins using shape complementarity, desolvation, and electrostatics , 2002, Proteins.

[38]  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..

[39]  Rolando Perez-Pineiro,et al.  Small Molecule Inhibitors of ERCC1-XPF Protein-Protein Interaction Synergize Alkylating Agents in Cancer Cells , 2013, Molecular Pharmacology.

[40]  C. Langer Emerging Immunotherapies in the Treatment of Non–small Cell Lung Cancer (NSCLC): The Role of Immune Checkpoint Inhibitors , 2015, American journal of clinical oncology.

[41]  K. Lafferty,et al.  A new analysis of allogeneic interactions. , 1975, The Australian journal of experimental biology and medical science.

[42]  S. Almo,et al.  Crystal structure of the complex between programmed death-1 (PD-1) and its ligand PD-L2 , 2008, Proceedings of the National Academy of Sciences.

[43]  Gerhard Klebe,et al.  PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations , 2007, Nucleic Acids Res..

[44]  Marawan Ahmed,et al.  Assessment of new anti-HER2 ligands using combined docking, QM/MM scoring and MD simulation. , 2013, Journal of molecular graphics & modelling.

[45]  Antoni Ribas,et al.  Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma. , 2013, The New England journal of medicine.

[46]  A. Dömling,et al.  Programmed death-1: therapeutic success after more than 100 years of cancer immunotherapy. , 2014, Angewandte Chemie.

[47]  P. Hwu,et al.  A New Approach to Simultaneously Quantify Both TCR α- and β-Chain Diversity after Adoptive Immunotherapy , 2012, Clinical Cancer Research.

[48]  Z. Weng,et al.  ZDOCK: An initial‐stage protein‐docking algorithm , 2003, Proteins.

[49]  Zhiping Weng,et al.  Accelerating Protein Docking in ZDOCK Using an Advanced 3D Convolution Library , 2011, PloS one.

[50]  Michael Houghton,et al.  A human ether-á-go-go-related (hERG) ion channel atomistic model generated by long supercomputer molecular dynamics simulations and its use in predicting drug cardiotoxicity. , 2014, Toxicology letters.

[51]  D. Y. Lin,et al.  The PD-1/PD-L1 complex resembles the antigen-binding Fv domains of antibodies and T cell receptors , 2008, Proceedings of the National Academy of Sciences.

[52]  Marawan Ahmed,et al.  In silico design: extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2. , 2013, Journal of molecular graphics & modelling.

[53]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[54]  Jack A. Tuszynski,et al.  A Refined Model of the HCV NS5A Protein Bound to Daclatasvir Explains Drug-Resistant Mutations and Activity against Divergent Genotypes , 2015, J. Chem. Inf. Model..

[55]  Jianyin Shao,et al.  Clustering Molecular Dynamics Trajectories: 1. Characterizing the Performance of Different Clustering Algorithms. , 2007, Journal of chemical theory and computation.

[56]  G. Freeman,et al.  PD-1 and its ligands in tolerance and immunity. , 2008, Annual review of immunology.

[57]  J. Tuszynski,et al.  DNA polymerase beta (pol β) inhibitors: a comprehensive overview. , 2012, Drug discovery today.

[58]  Marawan Ahmed,et al.  In silico investigation of lactone and thiolactone inhibitors in bacterial quorum sensing using molecular modeling , 2013, 1305.3691.

[59]  G. Freeman,et al.  The B7 family revisited. , 2005, Annual review of immunology.

[60]  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.