Epistasis arises from shifting the rate-limiting step during enzyme evolution of a ß-lactamase.

Epistasis, the non-additive effect of mutations, can provide combinatorial improvements to enzyme activity that substantially exceed the gains from individual mutations. Yet the molecular mechanisms of epistasis remain elusive, undermining our ability to predict pathogen evolution and engineer biocatalysts. Here we reveal how directed evolution of a ß-lactamase yielded highly epistatic activity enhancements. Evolution selected four mutations that increase antibiotic resistance 40-fold, despite their marginal individual effects (≤2-fold). Synergistic improvements coincided with the introduction of super-stochiometric burst kinetics, indicating that epistasis is rooted in the enzyme's conformational dynamics. Our analysis reveals that epistasis stemmed from distinct effects of each mutation on the catalytic cycle. The initial mutation increased protein flexibility and accelerated substrate binding, which is rate-limiting in the wild-type enzyme. Subsequent mutations predominantly boosted the chemical steps by fine-tuning substrate interactions. Our work identifies an overlooked cause for epistasis: changing the rate-limiting step can result in substantial synergy that boosts enzyme activity.

Electric Fields Are a Key Determinant of Carbapenemase Activity in Class A ß-Lactamases.

Resistance to antibiotics is a public health crisis. Although carbapenems are less susceptible to resistance than other ß-lactam antibiotics, ß-lactamases mediating resistance against these drugs are spreading. Here, we dissect the contributions of electric fields to carbapenemase activity in class A ß-lactamases. We perform QM/MM molecular dynamics simulations of meropenem acyl-enzyme hydrolysis that correctly discriminate carbapenemases. Electric field analysis shows that active-site fields in the deacylation transition state and tetrahedral intermediate are important determinants of activity. The active-site fields identify several residues, some distal, that distinguish efficient carbapenemases. Our field analysis script (www.github.com/bunzela/FieldTools) may help in understanding and combating antibiotic resistance.

emle-engine: A Flexible Electrostatic Machine Learning Embedding Package for Multiscale Molecular Dynamics Simulations.

We present in this work the emle-engine package (https://github.com/chemle/emle-engine)─the implementation of a new machine learning embedding scheme for hybrid machine learning potential/molecular-mechanics (ML/MM) dynamics simulations. The package is based on an embedding scheme that uses a physics-based model of the electronic density and induction with a handful of tunable parameters derived from in vacuo properties of the subsystem to be embedded. This scheme is completely independent of the in vacuo potential and requires only the positions of the atoms of the machine learning subsystem and the positions and partial charges of the molecular mechanics environment. These characteristics allow emle-engine to be employed in existing QM/MM software. We demonstrate that the implemented electrostatic machine learning embedding scheme (named EMLE) is stable in enhanced sampling molecular dynamics simulations. Through the calculation of free energy surfaces of alanine dipeptide in water with two different ML options for the in vacuo potential and three embedding models, we test the performance of EMLE. When compared to the reference DFT/MM surface, the EMLE embedding is clearly superior to the MM one based on fixed partial charges. The configurational dependence of the electronic density and the inclusion of the induction energy introduced by the EMLE model leads to a systematic reduction in the average error of the free energy surface when compared to MM embedding. By enabling the usage of EMLE embedding in practical ML/MM simulations, emle-engine will make it possible to accurately model systems and processes that feature significant variations in the charge distribution of the ML subsystem and/or the interacting environment.

Cooperative Conformational Transitions Underpin the Activation Heat Capacity in the Temperature Dependence of Enzyme Catalysis.

Many enzymes display non-Arrhenius behavior with curved Arrhenius plots in the absence of denaturation. There has been significant debate about the origin of this behavior and recently the role of the activation heat capacity (ΔCP⧧) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of ΔCP⧧ (arising from narrowing of the conformational space along the reaction coordinate), then curved Arrhenius plots are a consequence. To investigate these phenomena in detail, we have collected high precision temperature-rate data over a wide temperature interval for a model glycosidase enzyme MalL, and a series of mutants that change the temperature-dependence of the enzyme-catalyzed rate. We use these data to test a range of models including macromolecular rate theory (MMRT) and an equilibrium model. In addition, we have performed extensive molecular dynamics (MD) simulations to characterize the conformational landscape traversed by MalL in the enzyme-substrate complex and an enzyme-transition state complex. We have crystallized the enzyme in a transition state-like conformation in the absence of a ligand and determined an X-ray crystal structure at very high resolution (1.10 Å). We show (using simulation) that this enzyme-transition state conformation has a more restricted conformational landscape than the wildtype enzyme. We coin the term "transition state-like conformation (TLC)" to apply to this state of the enzyme. Together, these results imply a cooperative conformational transition between an enzyme-substrate conformation (ES) and a transition-state-like conformation (TLC) that precedes the chemical step. We present a two-state model as an extension of MMRT (MMRT-2S) that describes the data along with a convenient approximation with linear temperature dependence of the activation heat capacity (MMRT-1L) that can be used where fewer data points are available. Our model rationalizes disparate behavior seen for MalL and previous results for a thermophilic alcohol dehydrogenase and is consistent with a raft of data for other enzymes. Our model can be used to characterize the conformational changes required for enzyme catalysis and provides insights into the role of cooperative conformational changes in transition state stabilization that are accompanied by changes in heat capacity for the system along the reaction coordinate. TLCs are likely to be of wide importance in understanding the temperature dependence of enzyme activity and other aspects of enzyme catalysis.

4-Thiaproline accelerates the slow folding phase of proteins containing cis prolines in the native state by two orders of magnitude.

The cis/trans isomerization of peptidyl-prolyl peptide bonds is often the bottleneck of the refolding reaction for proteins containing cis proline residues in the native state. Proline (Pro) analogues, especially C4-substituted fluoroprolines, have been widely used in protein engineering to enhance the thermodynamic stability of peptides and proteins and to investigate folding kinetics. 4-thiaproline (Thp) has been shown to bias the ring pucker of Pro, to increase the cis population percentage of model peptides in comparison to Pro, and to diminish the activation energy barrier for the cis/trans isomerization reaction. Despite its intriguing properties, Thp has been seldom incorporated into proteins. Moreover, the impact of Thp on the folding kinetics of globular proteins has never been reported. In this study, we show that upon incorporation of Thp at cisPro76 into the thioredoxin variant Trx1P the half-life of the refolding reaction decreased from ~2 h to ~35 s. A dramatic acceleration of the refolding rate could be observed also for the protein pseudo wild-type barstar upon replacement of cisPro48 with Thp. Quantum chemical calculations suggested that the replacement of the Cγ H2 group by a sulfur atom in the pyrrolidine ring, might lower the barrier for cis/trans rotation due to a weakened peptide bond. The protein variants retained their thermodynamic stability upon incorporation of Thp, while the catalytic and enzymatic activities of the modified Trx1P remained unchanged. Our results show that the Pro isostere Thp might accelerate the rate of the slow refolding reaction for proteins containing cis proline residues in the native state, independent from the local structural environment.

Active Site Loop Engineering Abolishes Water Capture in Hydroxylating Sesquiterpene Synthases.

Terpene synthases (TS) catalyze complex reactions to produce a diverse array of terpene skeletons from linear isoprenyl diphosphates. Patchoulol synthase (PTS) from Pogostemon cablin converts farnesyl diphosphate into patchoulol. Using simulation-guided engineering, we obtained PTS variants that eliminate water capture. Further, we demonstrate that modifying the structurally conserved Hα-1 loop also reduces hydroxylation in PTS, as well as in germacradiene-11-ol synthase (Gd11olS), leading to cyclic neutral intermediates as products, including α-bulnesene (PTS) and isolepidozene (Gd11olS). Hα-1 loop modification could be a general strategy for engineering sesquiterpene synthases to produce complex cyclic hydrocarbons without the need for structure determination or modeling.

Free energy along drug-protein binding pathways interactively sampled in virtual reality.

We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight.

Promiscuous recognition of MR1 drives self-reactive mucosal-associated invariant T cell responses.

Mucosal-associated invariant T (MAIT) cells use canonical semi-invariant T cell receptors (TCR) to recognize microbial riboflavin precursors displayed by the antigen-presenting molecule MR1. The extent of MAIT TCR crossreactivity toward physiological, microbially unrelated antigens remains underexplored. We describe MAIT TCRs endowed with MR1-dependent reactivity to tumor and healthy cells in the absence of microbial metabolites. MAIT cells bearing TCRs crossreactive toward self are rare but commonly found within healthy donors and display T-helper-like functions in vitro. Experiments with MR1-tetramers loaded with distinct ligands revealed significant crossreactivity among MAIT TCRs both ex vivo and upon in vitro expansion. A canonical MAIT TCR was selected on the basis of extremely promiscuous MR1 recognition. Structural and molecular dynamic analyses associated promiscuity to unique TCRß-chain features that were enriched within self-reactive MAIT cells of healthy individuals. Thus, self-reactive recognition of MR1 represents a functionally relevant indication of MAIT TCR crossreactivity, suggesting a potentially broader role of MAIT cells in immune homeostasis and diseases, beyond microbial immunosurveillance.

Interrogation of an Enzyme Library Reveals the Catalytic Plasticity of Naturally Evolved [4+2] Cyclases.

Stereoselective carbon-carbon bond forming reactions are quintessential transformations in organic synthesis. One example is the Diels-Alder reaction, a [4+2] cycloaddition between a conjugated diene and a dienophile to form cyclohexenes. The development of biocatalysts for this reaction is paramount for unlocking sustainable routes to a plethora of important molecules. To obtain a comprehensive understanding of naturally evolved [4+2] cyclases, and to identify hitherto uncharacterised biocatalysts for this reaction, we constructed a library comprising forty-five enzymes with reported or predicted [4+2] cycloaddition activity. Thirty-one library members were successfully produced in recombinant form. In vitro assays employing a synthetic substrate incorporating a diene and a dienophile revealed broad-ranging cycloaddition activity amongst these polypeptides. The hypothetical protein Cyc15 was found to catalyse an intramolecular cycloaddition to generate a novel spirotetronate. The crystal structure of this enzyme, along with docking studies, establishes the basis for stereoselectivity in Cyc15, as compared to other spirotetronate cyclases.

Rare and Common Variants in GALNT3 May Affect Bone Mass Independently of Phosphate Metabolism.

Anabolic treatment options for osteoporosis remain limited. One approach to discovering novel anabolic drug targets is to identify genetic causes of extreme high bone mass (HBM). We investigated a pedigree with unexplained HBM within the UK HBM study, a national cohort of probands with HBM and their relatives. Whole exome sequencing (WES) in a family with HBM identified a rare heterozygous missense variant (NM_004482.4:c.1657C > T, p.Arg553Trp) in GALNT3, segregating appropriately. Interrogation of data from the UK HBM study and the Anglo-Australasian Osteoporosis Genetics Consortium (AOGC) revealed an unrelated individual with HBM with another rare heterozygous variant (NM_004482.4:c.831 T > A, p.Asp277Glu) within the same gene. In silico protein modeling predicted that p.Arg553Trp would disrupt salt-bridge interactions, causing instability of GALNT3, and that p.Asp277Glu would disrupt manganese binding and consequently GALNT3 catalytic function. Bi-allelic loss-of-function GALNT3 mutations alter FGF23 metabolism, resulting in hyperphosphatemia and causing familial tumoral calcinosis (FTC). However, bone mineral density (BMD) in FTC cases, when reported, has been either normal or low. Common variants in the GALNT3 locus show genome-wide significant associations with lumbar, femoral neck, and total body BMD. However, no significant associations with BMD are observed at loci coding for FGF23, its receptor FGFR1, or coreceptor klotho. Mendelian randomization analysis, using expression quantitative trait loci (eQTL) data from primary human osteoblasts and genome-wide association studies data from UK Biobank, suggested increased expression of GALNT3 reduces total body, lumbar spine, and femoral neck BMD but has no effect on phosphate concentrations. In conclusion, rare heterozygous loss-of-function variants in GALNT3 may cause HBM without altering phosphate concentration. These findings suggest that GALNT3 may affect BMD through pathways other than FGF23 regulation, the identification of which may yield novel anabolic drug targets for osteoporosis. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Synthetic cannabinoid receptor agonists are monoamine oxidase-A selective inhibitors.

Synthetic cannabinoid receptor agonists (SCRAs) are one of the fastest growing classes of recreational drugs. Despite their growth in use, their vast chemical diversity and rapidly changing landscape of structures make understanding their effects challenging. In particular, the side effects for SCRA use are extremely diverse, but notably include severe outcomes such as cardiac arrest. These side effects appear at odds with the main putative mode of action, as full agonists of cannabinoid receptors. We have hypothesized that SCRAs may act as MAO inhibitors, owing to their structural similarity to known monoamine oxidase inhibitors (MAOI's) as well as matching clinical outcomes (hypertensive crisis) of 'monoaminergic toxicity' for users of MAOIs and some SCRA use. We have studied the potential for SCRA-mediated inhibition of MAO-A and MAO-B via a range of SCRAs used commonly in the UK, as well as structural analogues to prove the atomistic determinants of inhibition. By combining in silico and experimental kinetic studies we demonstrate that SCRAs are MAO-A-specific inhibitors and their affinity can vary significantly between SCRAs, most notably affected by the nature of the SCRA 'head' group. Our data allow us to posit a putative mechanism of inhibition. Crucially our data demonstrate that SCRA activity is not limited to just cannabinoid receptor agonism and that alternative interactions might account for some of the diversity of the observed side effects and that these effects can be SCRA-specific.

The Role of Cytochrome†P450 AbyV in the Final Stages of Abyssomicin†C Biosynthesis.

Abyssomicin C and its atropisomer are potent inhibitors of bacterial folate metabolism. They possess complex polycyclic structures, and their biosynthesis has been shown to involve several unusual enzymatic transformations. Using a combination of synthesis and in vitro assays we reveal that AbyV, a cytochrome P450 enzyme from the aby gene cluster, catalyses a key late-stage epoxidation required for the installation of the characteristic ether-bridged core of abyssomicin C. The X-ray crystal structure of AbyV has been determined, which in combination with molecular dynamics simulations provides a structural framework for our functional data. This work demonstrates the power of combining selective carbon-13 labelling with NMR spectroscopy as a sensitive tool to interrogate enzyme-catalysed reactions in vitro with no need for purification.

The Role of Cytochrome†P450 AbyV in the Final Stages of Abyssomicin†C Biosynthesis.

Abyssomicin C and its atropisomer are potent inhibitors of bacterial folate metabolism. They possess complex polycyclic structures, and their biosynthesis has been shown to involve several unusual enzymatic transformations. Using a combination of synthesis and in vitro assays we reveal that AbyV, a cytochrome P450 enzyme from the aby gene cluster, catalyses a key late-stage epoxidation required for the installation of the characteristic ether-bridged core of abyssomicin C. The X-ray crystal structure of AbyV has been determined, which in combination with molecular dynamics simulations provides a structural framework for our functional data. This work demonstrates the power of combining selective carbon-13 labelling with NMR spectroscopy as a sensitive tool to interrogate enzyme-catalysed reactions in vitro with no need for purification.

Structurally Informed Mutagenesis of a Stereochemically Promiscuous Aldolase Produces Mutants That Catalyze the Diastereoselective Syntheses of All Four Stereoisomers of 3-Deoxy-hexulosonic Acid.

A 2-keto-3-deoxygluconate aldolase from the hyperthermophile Sulfolobus solfataricus catalyzes the nonstereoselective aldol reaction of pyruvate and d-glyceraldehyde to produce 2-keto-3-deoxygluconate (d-KDGlc) and 2-keto-3-deoxy-d-galactonate (d-KDGal). Previous investigations into curing the stereochemical promiscuity of this hyperstable aldolase used high-resolution structures of the aldolase bound to d-KDGlc or d-KDGal to identify critical amino acids involved in substrate binding for mutation. This structure-guided approach enabled mutant variants to be created that could stereoselectively catalyze the aldol reaction of pyruvate and natural d-glyceraldehyde to selectively afford d-KDGlc or d-KDGal. Here we describe the creation of two further mutants of this Sulfolobus aldolase that can be used to catalyze aldol reactions between pyruvate and non-natural l-glyceraldehyde to enable the diastereoselective synthesis of l-KDGlc and l-KDGal. High-resolution crystal structures of all four variant aldolases have been determined (both unliganded and liganded), including Variant 1 with d-KDGlc, Variant 2 with pyruvate, Variant 3 with l-KDGlc, and Variant 4 with l-KDGal. These structures have enabled us to rationalize the observed changes in diastereoselectivities in these variant-catalyzed aldol reactions at a molecular level. Interestingly, the active site of Variant 4 was found to be sufficiently flexible to enable catalytically important amino acids to be replaced while still retaining sufficient enzymic activity to enable production of l-KDGal.

Therapeutic high affinity T cell receptor targeting a KRASG12D cancer neoantigen.

Neoantigens derived from somatic mutations are specific to cancer cells and are ideal targets for cancer immunotherapy. KRAS is the most frequently mutated oncogene and drives the pathogenesis of several cancers. Here we show the identification and development of an affinity-enhanced T cell receptor (TCR) that recognizes a peptide derived from the most common KRAS mutant, KRASG12D, presented in the context of HLA-A*11:01. The affinity of the engineered TCR is increased by over one million-fold yet fully able to distinguish KRASG12D over KRASWT. While crystal structures reveal few discernible differences in TCR interactions with KRASWT versus KRASG12D, thermodynamic analysis and molecular dynamics simulations reveal that TCR specificity is driven by differences in indirect electrostatic interactions. The affinity enhanced TCR, fused to a humanized anti-CD3 scFv, enables selective killing of cancer cells expressing KRASG12D. Our work thus reveals a molecular mechanism that drives TCR selectivity and describes a soluble bispecific molecule with therapeutic potential against cancers harboring a common shared neoantigen.

QM/MM Simulations Reveal the Determinants of Carbapenemase Activity in Class A ß-Lactamases.

ß-lactam antibiotic resistance in Gram-negative bacteria, primarily caused by ß-lactamase enzymes that hydrolyze the ß-lactam ring, has become a serious clinical problem. Carbapenems were formerly considered "last resort" antibiotics because they escaped breakdown by most ß-lactamases, due to slow deacylation of the acyl-enzyme intermediate. However, an increasing number of Gram-negative bacteria now produce ß-lactamases with carbapenemase activity: these efficiently hydrolyze the carbapenem ß-lactam ring, severely limiting the treatment of some bacterial infections. Here, we use quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reactions of acyl-enzyme complexes of eight ß-lactamases of class A (the most widely distributed ß-lactamase group) with the carbapenem meropenem to investigate differences between those inhibited by carbapenems (TEM-1, SHV-1, BlaC, and CTX-M-16) and those that hydrolyze them (SFC-1, KPC-2, NMC-A, and SME-1). QM/MM molecular dynamics simulations confirm the two enzyme groups to differ in the preferred acyl-enzyme orientation: carbapenem-inhibited enzymes favor hydrogen bonding of the carbapenem hydroxyethyl group to deacylating water (DW). QM/MM simulations of deacylation give activation free energies in good agreement with experimental hydrolysis rates, correctly distinguishing carbapenemases. For the carbapenem-inhibited enzymes, free energies for deacylation are significantly higher than for the carbapenemases, even when the hydroxyethyl group was restrained to prevent interaction with the DW. Analysis of these simulations, and additional simulations of mutant enzymes, shows how factors including the hydroxyethyl orientation, the active site volume, and architecture (conformations of Asn170 and Asn132; organization of the oxyanion hole; and the Cys69-Cys238 disulfide bond) collectively determine catalytic efficiency toward carbapenems.

Multiscale Simulations Identify Origins of Differential Carbapenem Hydrolysis by the OXA-48 ß-Lactamase.

OXA-48 ß-lactamases are frequently encountered in bacterial infections caused by carbapenem-resistant Gram-negative bacteria. Due to the importance of carbapenems in the treatment of healthcare-associated infections and the increasingly wide dissemination of OXA-48-like enzymes on plasmids, these ß-lactamases are of high clinical significance. Notably, OXA-48 hydrolyzes imipenem more efficiently than other commonly used carbapenems, such as meropenem. Here, we use extensive multiscale simulations of imipenem and meropenem hydrolysis by OXA-48 to dissect the dynamics and to explore differences in the reactivity of the possible conformational substates of the respective acylenzymes. Quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reaction for both substrates demonstrate that deacylation is favored when the 6α-hydroxyethyl group is able to hydrogen bond to the water molecule responsible for deacylation but disfavored by the increasing hydration of either oxygen of the carboxylated Lys73 general base. Differences in free energy barriers calculated from the QM/MM simulations correlate well with the experimentally observed differences in hydrolytic efficiency between meropenem and imipenem. We conclude that the impaired breakdown of meropenem, compared to imipenem, which arises from a subtle change in the hydrogen bonding pattern between the deacylating water molecule and the antibiotic, is most likely induced by the meropenem 1ß-methyl group. In addition to increased insights into carbapenem breakdown by OXA ß-lactamases, which may aid in future efforts to design antibiotics or inhibitors, our approach exemplifies the combined use of atomistic simulations in determining the possible different enzyme-substrate substates and their influence on enzyme reaction kinetics.

Path to Actinorhodin: Regio- and Stereoselective Ketone Reduction by a Type II Polyketide Ketoreductase Revealed in Atomistic Detail.

In type II polyketide synthases (PKSs), which typically biosynthesize several antibiotic and antitumor compounds, the substrate is a growing polyketide chain, shuttled between individual PKS enzymes, while covalently tethered to an acyl carrier protein (ACP): this requires the ACP interacting with a series of different enzymes in succession. During biosynthesis of the antibiotic actinorhodin, produced by Streptomyces coelicolor, one such key binding event is between an ACP carrying a 16-carbon octaketide chain (actACP) and a ketoreductase (actKR). Once the octaketide is bound inside actKR, it is likely cyclized between C7 and C12 and regioselective reduction of the ketone at C9 occurs: how these elegant chemical and conformational changes are controlled is not yet known. Here, we perform protein-protein docking, protein NMR, and extensive molecular dynamics simulations to reveal a probable mode of association between actACP and actKR; we obtain and analyze a detailed model of the C7-C12-cyclized octaketide within the actKR active site; and we confirm this model through multiscale (QM/MM) reaction simulations of the key ketoreduction step. Molecular dynamics simulations show that the most thermodynamically stable cyclized octaketide isomer (7R,12R) also gives rise to the most reaction competent conformations for ketoreduction. Subsequent reaction simulations show that ketoreduction is stereoselective as well as regioselective, resulting in an S-alcohol. Our simulations further indicate several conserved residues that may be involved in selectivity of C7-12 cyclization and C9 ketoreduction. Detailed insights obtained on ACP-based substrate presentation in type II PKSs can help design ACP-ketoreductase systems with altered regio- or stereoselectivity.

Reliable In Silico Ranking of Engineered Therapeutic TCR Binding Affinities with MMPB/GBSA.

Accurate and efficient in silico ranking of protein-protein binding affinities is useful for protein design with applications in biological therapeutics. One popular approach to rank binding affinities is to apply the molecular mechanics Poisson-Boltzmann/generalized Born surface area (MMPB/GBSA) method to molecular dynamics (MD) trajectories. Here, we identify protocols that enable the reliable evaluation of T-cell receptor (TCR) variants binding to their target, peptide-human leukocyte antigens (pHLAs). We suggest different protocols for variant sets with a few (≤4) or many mutations, with entropy corrections important for the latter. We demonstrate how potential outliers could be identified in advance and that just 5-10 replicas of short (4 ns) MD simulations may be sufficient for the reproducible and accurate ranking of TCR variants. The protocols developed here can be applied toward in silico screening during the optimization of therapeutic TCRs, potentially reducing both the cost and time taken for biologic development.