Using Molecular Dynamics Simulations to Interrogate T Cell Receptor Non-Equilibrium Kinetics

A Computational Strategy for the Rapid Identification and Ranking of Patient-Specific T Cell Receptors Bound to Neoantigens

T cell receptor (TCR) recognition of a peptide-major histocompatibility complex (pMHC) is crucial for adaptive immune response. The identification of therapeutically relevant TCR-pMHC protein pairs is a bottleneck in the implementation of TCR-based immunotherapies. The ability to computationally design TCRs to target a specific pMHC requires automated integration of next-generation sequencing, protein-protein structure prediction, molecular dynamics, and TCR ranking. A pipeline to evaluate patient-specific, sequence-based TCRs to a target pMHC is presented. Using the three most frequently expressed TCRs from 16 colorectal cancer patients, the protein-protein structure of the TCRs to the target CEA peptide-MHC is predicted using Modeller and ColabFold. TCR-pMHC structures are compared using automated equilibration and successive analysis. ColabFold generated configurations require an ≈2.5× reduction in equilibration time of TCR-pMHC structures compared to Modeller. The structural differences between Modeller and ColabFold are demonstrated by root mean square deviation (≈0.20 nm) between clusters of equilibrated configurations, which impact the number of hydrogen bonds and Lennard-Jones contacts between the TCR and pMHC. TCR ranking criteria that may prioritize TCRs for evaluation of in vitro immunogenicity are identified, and this ranking is validated by comparing to state-of-the-art machine learning-based methods trained to predict the probability of TCR-pMHC binding.

Inhibition of platelet adhesion to exposed subendothelial collagen by steric hindrance with blocking peptide nanoparticles

The inhibition of platelet adhesion to collagen in exposed vessels represents an innovative approach to the treatment of atherosclerosis and thrombosis. This study aimed to engineer peptide-based nanoparticles that prevent platelet binding to subendothelial collagen by engaging with collagen with high affinity. We examined the interactions between integrin α2/ glycoprotein VI/ von Willebrand factor A3 domain and collagen, as well as between the synthesized peptide nanoparticles and collagen, utilizing molecular dynamics simulations and empirical assays. Our findings indicated that the bond between von Willebrand factor and collagen was more robust. Specifically, the sequences SITTIDV, VDVMQRE, and YLTSEMH in von Willebrand factor were identified as essential for its attachment to collagen. Based on these sequences, three peptide nanoparticles were synthesized (BPa: Capric-GNNQQNYK-SITTIDV, BPb: Capric-GNNQQNYK-VDVMQRE, BPc: Capric-GNNQQNYK-YLTSEMH), each displaying significant affinity towards collagen. Of these, the BPa nanoparticles exhibited the most potent interaction with collagen, leading to a 75% reduction in platelet adhesion.

Repositioning of anti-infective compounds against monkeypox virus core cysteine proteinase: a molecular dynamics study

Monkeypox virus (MPXV) core cysteine proteinase (CCP) is one of the major drug targets used to examine the inhibitory action of chemical moieties. In this study, an in silico technique was applied to screen 1395 anti-infective compounds to find out the potential molecules against the MPXV-CCP. The top five hits were selected after screening and processed for exhaustive docking based on the docked score of ≤ -9.5 kcal/mol. Later, the top three hits based on the exhaustive-docking score and interaction profile were selected to perform MD simulations. The overall RMSD suggested that two compounds, SC75741 and ammonium glycyrrhizinate, showed a highly stable complex with a standard deviation of 0.18 and 0.23 nm, respectively. Later, the MM/GBSA binding free energies of complexes showed significant binding strength with ΔGTOTAL from -21.59 to -15 kcal/mol. This report reported the potential inhibitory activity of SC75741 and ammonium glycyrrhizinate against MPXV-CCP by competitively inhibiting the binding of the native substrate.

Mechanical forces amplify TCR mechanotransduction in T cell activation and function

Adoptive T cell immunotherapies, including engineered T cell receptor (eTCR) and chimeric antigen receptor (CAR) T cell immunotherapies, have shown efficacy in treating a subset of hematologic malignancies, exhibit promise in solid tumors, and have many other potential applications, such as in fibrosis, autoimmunity, and regenerative medicine. While immunoengineering has focused on designing biomaterials to present biochemical cues to manipulate T cells ex vivo and in vivo, mechanical cues that regulate their biology have been largely underappreciated. This review highlights the contributions of mechanical force to several receptor-ligand interactions critical to T cell function, with central focus on the TCR-peptide-loaded major histocompatibility complex (pMHC). We then emphasize the role of mechanical forces in (i) allosteric strengthening of the TCR-pMHC interaction in amplifying ligand discrimination during T cell antigen recognition prior to activation and (ii) T cell interactions with the extracellular matrix. We then describe approaches to design eTCRs, CARs, and biomaterials to exploit TCR mechanosensitivity in order to potentiate T cell manufacturing and function in adoptive T cell immunotherapy.

A dynamic biomimetic model of the membrane-bound CD4-CD3-TCR complex during pMHC disengagement

The coordinated (dis)engagement of the membrane-bound T cell receptor (TCR)-CD3-CD4 complex from the peptide-major histocompatibility complex (pMHC) is fundamental to TCR signal transduction and T cell effector function. As such, an atomic-scale understanding would not only enhance our basic understanding of the adaptive immune response but would also accelerate the rational design of TCRs for immunotherapy. In this study, we explore the impact of the CD4 coreceptor on the TCR-pMHC (dis)engagement by constructing a molecular-level biomimetic model of the CD3-TCR-pMHC and CD4-CD3-TCR-pMHC complexes within a lipid bilayer. After allowing the system complexes to equilibrate (engage), we use steered molecular dynamics to dissociate (disengage) the pMHC. We find that 1) the CD4 confines the pMHC closer to the T cell by 1.8 nm at equilibrium; 2) CD4 confinement shifts the TCR along the MHC binding groove engaging a different set of amino acids and enhancing the TCR-pMHC bond lifetime; 3) the CD4 translocates under load increasing the interaction strength between the CD4-pMHC, CD4-TCR, and CD4-CD3; and 4) upon dissociation, the CD3-TCR complex undergoes structural oscillation and increased energetic fluctuation between the CD3-TCR and CD3-lipids. These atomic-level simulations provide mechanistic insight on how the CD4 coreceptor impacts TCR-pMHC (dis)engagement. More specifically, our results provide further support (enhanced bond lifetime) for a force-dependent kinetic proofreading model and identify an alternate set of amino acids in the TCR that dominate the TCR-pMHC interaction and could thus impact the design of TCRs for immunotherapy.

The Energetic Landscape of Catch Bonds in TCR Interfaces

Recognition of peptide/MHC complexes by αβ TCRs has traditionally been viewed through the lens of conventional receptor-ligand theory. Recent work, however, has shown that TCR recognition and T cell signaling can be profoundly influenced and tuned by mechanical forces. One outcome of applied force is the catch bond, where TCR dissociation rates decrease (half-lives increase) when limited force is applied. Although catch bond behavior is believed to be widespread in biology, its counterintuitive nature coupled with the difficulties of describing mechanisms at the structural level have resulted in considerable mystique. In this review, we demonstrate that viewing catch bonds through the lens of energy landscapes, barriers, and the ensuing reaction rates can help demystify catch bonding and provide a foundation on which atomic-level TCR catch bond mechanisms can be built.

The Enigmatic Nature of the TCR-pMHC Interaction: Implications for CAR-T and TCR-T Engineering

The interaction of the T-cell receptor (TCR) with a peptide in the major histocompatibility complex (pMHC) plays a central role in the adaptive immunity of higher chordates. Due to the high specificity and sensitivity of this process, the immune system quickly recognizes and efficiently responds to the appearance of foreign and altered self-antigens. This is important for ensuring anti-infectious and antitumor immunity, in addition to maintaining self-tolerance. The most common parameter used for assessing the specificity of TCR-pMHC interaction is affinity. This thermodynamic characteristic is widely used not only in various theoretical aspects, but also in practice, for example, in the engineering of various T-cell products with a chimeric (CAR-T) or artificial (TCR-engineered T-cell) antigen receptor. However, increasing data reveal the fact that, in addition to the thermodynamic component, the specificity of antigen recognition is based on the kinetics and mechanics of the process, having even greater influence on the selectivity of the process and T lymphocyte activation than affinity. Therefore, the kinetic and mechanical aspects of antigen recognition should be taken into account when designing artificial antigen receptors, especially those that recognize antigens in the MHC complex. This review describes the current understanding of the nature of the TCR-pMHC interaction, in addition to the thermodynamic, kinetic, and mechanical principles underlying the specificity and high sensitivity of this interaction.

Electrostatics in Computational Biophysics and Its Implications for Disease Effects

This review outlines the role of electrostatics in computational molecular biophysics and its implication in altering wild-type characteristics of biological macromolecules, and thus the contribution of electrostatics to disease mechanisms. The work is not intended to review existing computational approaches or to propose further developments. Instead, it summarizes the outcomes of relevant studies and provides a generalized classification of major mechanisms that involve electrostatic effects in both wild-type and mutant biological macromolecules. It emphasizes the complex role of electrostatics in molecular biophysics, such that the long range of electrostatic interactions causes them to dominate all other forces at distances larger than several Angstroms, while at the same time, the alteration of short-range wild-type electrostatic pairwise interactions can have pronounced effects as well. Because of this dual nature of electrostatic interactions, being dominant at long-range and being very specific at short-range, their implications for wild-type structure and function are quite pronounced. Therefore, any disruption of the complex electrostatic network of interactions may abolish wild-type functionality and could be the dominant factor contributing to pathogenicity. However, we also outline that due to the plasticity of biological macromolecules, the effect of amino acid mutation may be reduced, and thus a charge deletion or insertion may not necessarily be deleterious.

A molecular dynamics investigation of N-glycosylation effects on T-cell receptor kinetics

The binding interaction between the T-cell receptor (TCR) and peptide-major histocompatibility complex (pMHC) is modulated by several factors (known and unknown), however, investigations into effects of glycosylation are limited. A fully glycosylated computational model of the TCR bound to the pMHC is developed to investigate the effects of glycosylation on dissociation kinetics from the pMHC. Here, we examine the effects of N-glycosylation on TCR-pMHC bond strength using steered molecular dynamic simulations. N-glycosylation is a post-translational modification that adds sugar moieties to molecules and can modulate the activity of several immune molecules. Using a TCR-pMHC pair found in melanoma as a case study, our study demonstrates that N-glycosylation of the TCR-pMHC alters the proteins' conformation; increases the bond lifetime; and increases the number of hydrogen bonds and Lennard-Jones Contacts involved in the TCR-pMHC bond. We find that weak glycan-protein or glycan-glycan interactions impact the equilibrated structure of the TCR and pMHC leading to an increase in the overall bond strength of the TCR-pMHC complex including the duration and energetic strength under constant load. These results indicate that N-glycosylation plays an important role in the TCR-pMHC bond and should be considered in future computational and experimental studies.Communicated by Ramaswamy H. Sarma.