Shunzhou Wan

Peptide recognition by the T cell receptor: comparison of binding free energies from thermodynamic integration, Poisson–Boltzmann and linear interaction energy approximations

The binding to the T cell receptor of wild-type and variant HTLV-1 Tax peptide complexed to the major histocompatibility complex has been investigated by means of molecular dynamics simulations. The binding free energy difference is calculated using the molecular mechanics Poisson–Boltzmann surface area and linear interaction energy methods. These methods extract useful information on the binding energetics from simulations of the physical states of the ligands, which are more computationally expedient than the commonly used thermodynamic integration method. The successful reproduction of the relative binding free energies shows that these methods can be useful for free energy calculations and the rational design of drugs and vaccines.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large‐scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA‐A*0201 bound to a decameric tumor‐specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high‐performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent‐accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected.

Molecular basis of peptide recognition by the TCR: Affinity differences calculated using large scale computing

Free energy calculations of the wild-type and the variant human T cell lymphotropic virus type 1 Tax peptide presented by the MHC to the TCR have been performed using large scale massively parallel molecular dynamics simulations. The computed free energy difference (−1.86 ± 0.44 kcal/mol) using alchemical mutation-based thermodynamic integration agrees well with experimental data (−2.9 ± 0.2 kcal/mol). Our simulations exploit state-of-the-art hardware and codes whose algorithms have been optimized for supercomputing platforms. This enables us to simulate larger, more realistic biological systems for longer durations without the imposition of artificial constraints.

Rapid, Precise, and Reproducible Prediction of Peptide-MHC Binding Affinities from Molecular Dynamics That Correlate Well with Experiment.

The presentation of potentially pathogenic peptides by major histocompatibility complex (MHC) molecules is one of the most important processes in adaptive immune defense. Prediction of peptide-MHC (pMHC) binding affinities is therefore a principal objective of theoretical immunology. Machine learning techniques achieve good results if substantial experimental training data are available. Approaches based on structural information become necessary if sufficiently similar training data are unavailable for a specific MHC allele, although they have often been deemed to lack accuracy. In this study, we use a free energy method to rank the binding affinities of 12 diverse peptides bound by a class I MHC molecule HLA-A*02:01. The method is based on enhanced sampling of molecular dynamics calculations in combination with a continuum solvent approximation and includes estimates of the configurational entropy based on either a one or a three trajectory protocol. It produces precise and reproducible free energy estimates which correlate well with experimental measurements. If the results are combined with an amino acid hydrophobicity scale, then an extremely good ranking of peptide binding affinities emerges. Our approach is rapid, robust, and applicable to a wide range of ligand-receptor interactions without further adjustment.

Clinically driven design of multi-scale cancer models: the ContraCancrum project paradigm

The challenge of modelling cancer presents a major opportunity to improve our ability to reduce mortality from malignant neoplasms, improve treatments and meet the demands associated with the individualization of care needs. This is the central motivation behind the ContraCancrum project. By developing integrated multi-scale cancer models, ContraCancrum is expected to contribute to the advancement of in silico oncology through the optimization of cancer treatment in the patient-individualized context by simulating the response to various therapeutic regimens. The aim of the present paper is to describe a novel paradigm for designing clinically driven multi-scale cancer modelling by bringing together basic science and information technology modules. In addition, the integration of the multi-scale tumour modelling components has led to novel concepts of personalized clinical decision support in the context of predictive oncology, as is also discussed in the paper. Since clinical adaptation is an inelastic prerequisite, a long-term clinical adaptation procedure of the models has been initiated for two tumour types, namely non-small cell lung cancer and glioblastoma multiforme; its current status is briefly summarized.

The Effect of Mutations on Drug Sensitivity and Kinase Activity of Fibroblast Growth Factor Receptors: A Combined Experimental and Theoretical Study

Fibroblast growth factor receptors (FGFRs) are recognized therapeutic targets in cancer. We here describe insights underpinning the impact of mutations on FGFR1 and FGFR3 kinase activity and drug efficacy, using a combination of computational calculations and experimental approaches including cellular studies, X-ray crystallography and biophysical and biochemical measurements. Our findings reveal that some of the tested compounds, in particular TKI258, could provide therapeutic opportunity not only for patients with primary alterations in FGFR but also for acquired resistance due to the gatekeeper mutation. The accuracy of the computational methodologies applied here shows a potential for their wider application in studies of drug binding and in assessments of functional and mechanistic impacts of mutations, thus assisting efforts in precision medicine.

T-cell epitope prediction and immune complex simulation using molecular dynamics: state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Rapid and accurate ranking of binding affinities of epidermal growth factor receptor sequences with selected lung cancer drugs

The epidermal growth factor receptor (EGFR) is a major target for drugs in treating lung carcinoma. Mutations in the tyrosine kinase domain of EGFR commonly arise in human cancers, which can cause drug sensitivity or resistance by influencing the relative strengths of drug and ATP-binding. In this study, we investigate the binding affinities of two tyrosine kinase inhibitors—AEE788 and Gefitinib—to EGFR using molecular dynamics simulation. The interactions between these inhibitors and the EGFR kinase domain are analysed using multiple short (ensemble) simulations and the molecular mechanics/Poisson–Boltzmann solvent area (MM/PBSA) method. Here, we show that ensemble simulations correctly rank the binding affinities for these systems: we report the successful ranking of each drug binding to a variety of EGFR sequences and of the two drugs binding to a given sequence, using petascale computing resources, within a few days.

On the calculation of equilibrium thermodynamic properties from molecular dynamics

The purpose of statistical mechanics is to provide a route to the calculation of macroscopic properties of matter from their constituent microscopic components. It is well known that the macrostates emerge as ensemble averages of microstates. However, this is more often stated than implemented in computer simulation studies. Here we consider foundational aspects of statistical mechanics which are overlooked in most textbooks and research articles that purport to compute macroscopic behaviour from microscopic descriptions based on classical mechanics and show how due attention to these issues leads in directions which have not been widely appreciated in the field of molecular dynamics simulation.

Molecular Dynamics Simulation Reveals Structural and Thermodynamic Features of Kinase Activation by Cancer Mutations Within the Epidermal Growth Factor Receptor

The epidermal growth factor receptor (EGFR) is a major target for drugs in treating lung carcinoma as it promotes cell growth and tumor progression. Structural studies have demonstrated that EGFR exists in an equilibrium between catalytically active and inactive forms, and dramatic conformational transitions occur during its activation. It is known that EGFR mutations promote such conformational changes that affect its activation and drug efficacy. The most common point mutation in lung cancer patients is a leucine to arginine substitution at amino acid 834 (L834R). In a recent article, we have studied changes in drug binding affinities due to cancer mutations of EGFR using ensemble molecular dynamics (MD) simulations. Here, we address an enhanced activation mechanism thought to be associated with this mutation. Using extended timescale MD simulations, the structural and energetic properties are studied for both active and inactive conformations of EGFR. The thermodynamic stabilities of these two conformations are characterized by free energy landscapes estimated from molecular mechanics/Poisson–Boltzmann solvent area calculations. Our study reveals that the L834R mutation introduces conformational changes in both states, adjusting the relative stabilities of active and inactive conformations and hence the activation of the EGFR kinase.