Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reversible domain closure modulates GlnBP ligand binding affinity

Glutamine binding protein (GlnBP) is an Escherichia Coli periplasmic binding protein, which binds and carries glutamine to the inner membrane ATP-binding cassette (ABC) transporter. GlnBP binds the ligand with affinity around 0.1μM measured by isothermal titration calorimetry (ITC) and ligand binding stabilizes protein structure shown by its increase in thermodynamic stability. However, the molecular determinant of GlnBP ligand binding is not known. Electrostatic and hydrophobic interaction between GlnBP and glutamine are critical factors. We propose that the freedome of closure movement is also vital for ligand binding. In order to approve this hypothesis, we generate a series of mutants with different linker length that has different magnitude of domain closure. Mutants show different ligand binding affinity, which indicates that the propensity of domain closure determines the ligand binding affinity. Ligand binding triggers gradual ensemble conformational change. Structural changes upon ligand binding are monitored by combination of small angle x-ray scattering (SAXS) and NMR spectroscopy. Detailed structure characterization of GlnBP contributes to a better understanding of ligand binding and provides the structural basis for biosensor design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structure-guided approach to modulate small molecule binding to a promiscuous ligand-activated protein

Ligand-binding promiscuity in detoxification systems protects the body from toxicological harm but is a roadblock to drug development due to the difficulty in optimizing small molecules to both retain target potency and avoid metabolic events. Immense effort is invested in evaluating metabolism of molecules to develop safer, more effective treatments, but engineering specificity into or out of promiscuous proteins and their ligands is a challenging task. To better understand the promiscuous nature of detoxification networks, we have used X-ray crystallography to characterize a structural feature of pregnane X receptor (PXR), a nuclear receptor that is activated by diverse molecules (with different structures and sizes) to up-regulate transcription of drug metabolism genes. We found that large ligands expand PXR’s ligand-binding pocket, and the ligand-induced expansion occurs through a specific unfavorable compound-protein clash that likely contributes to reduced binding affinity. Removing the clash by compound modification resulted in more favorable binding modes with significantly enhanced binding affinity. We then engineered the unfavorable ligand-protein clash into a potent, small PXR ligand, resulting in marked reduction in PXR binding and activation. Structural analysis showed that PXR is remodeled, and the modified ligands reposition in the binding pocket to avoid clashes, but the conformational changes result in less favorable binding modes. Thus, ligand-induced binding pocket expansion increases ligand-binding potential of PXR but is an unfavorable event; therefore, drug candidates can be engineered to expand PXR’s ligand-binding pocket and reduce their safety liability due to PXR binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Stable binding to phosphatidylserine-containing membranes requires conserved arginine residues in tandem C domains of blood coagulation factor VIII

At sites of vascular damage, factor VIII (fVIII) is proteolytically activated by thrombin and binds to activated platelet surfaces with activated factor IX (fIXa) to form the intrinsic “tenase” complex. Previous structural and mutational studies of fVIII have identified the C1 and C2 domains in binding to negatively charged membrane surfaces through β-hairpin loops with solvent-exposed hydrophobic residues and a ring of positively charged basic residues. Several hemophilia A-associated mutations within the C domains are suggested to disrupt lipid binding, preventing formation of the intrinsic tenase complex. In this study, we devised a novel platform for generating recombinant C1, C2, and C1C2 domain constructs and performed mutagenesis of several charged residues proximal to the putative membrane binding region of each C domain. Binding measurements between phosphatidylserine (PS)-containing lipid membrane surfaces and fVIII C domains demonstrated an ionic strength dependence on membrane binding affinity. Mutations to basic residues adjacent to the surface-exposed hydrophobic regions of C1 and C2 differentially disrupted membrane binding, with abrogation of binding occurring for mutations to conserved arginine residues in the C1 (R2163) and C2 (R2320) domains. Lastly, we determined the X-ray crystal structure of the porcine fVIII C2 domain bound to o-phospho-L-serine, the polar headgroup of PS, which binds to a basic cleft and makes charge-charge contact with R2320. We conclude that basic clefts in the fVIII C domains bind to PS-containing membranes through conserved arginine residues via a C domain modularity, where each C domain possesses modest electrostatic-dependent affinity and tandem C domains are required for high affinity binding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Enhanced Binding of SARS-CoV-2 Spike Protein to Receptor by Distal Polybasic Cleavage Sites

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein plays a crucial role in binding the human cell receptor ACE2 that is required for viral entry. Many studies have been conducted to target the structures of RBD–ACE2 binding and to design RBD-targeting vaccines and drugs. Nevertheless, mutations distal from the SARS-CoV-2 RBD also impact its transmissibility and antibody can target non-RBD regions, suggesting the incomplete role of the RBD region in the spike protein–ACE2 binding. Here, in order to elucidate distant binding mechanisms, we analyze complexes of ACE2 with the wild-type spike protein and with key mutants via large-scale all-atom explicit solvent molecular dynamics simulations. We find that though distributed approximately 10 nm away from the RBD, the SARS-CoV-2 polybasic cleavage sites enhance, via electrostatic interactions and hydration, the RBD–ACE2 binding affinity. A negatively charged tetrapeptide (GluGluLeuGlu) is then designed to neutralize the positively charged arginine on the polybasic cleavage sites. We find that the tetrapeptide GluGluLeuGlu binds to one of the three polybasic cleavage sites of the SARS-CoV-2 spike protein lessening by 34% the RBD–ACE2 binding strength. This significant binding energy reduction demonstrates the feasibility to neutralize RBD–ACE2 binding by targeting this specific polybasic cleavage site. Our work enhances understanding of the binding mechanism of SARS-CoV-2 to ACE2, which may aid the design of therapeutics for COVID-19 infection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Tandem sialoglycan-binding modules in a Streptococcus sanguinis serine-rich repeat adhesin create target dependent avidity effects

Sialic acid–binding immunoglobulin-like lectins (Siglec)–like domains of streptococcal serine-rich repeat (SRR) adhesins recognize sialylated glycans on human salivary, platelet, and plasma glycoproteins via a YTRY sequence motif. The SRR adhesin from Streptococcus sanguinis strain SK1 has tandem sialoglycan-binding domains and has previously been shown to bind sialoglycans with high affinity. However, both domains contain substitutions within the canonical YTRY motif, making it unclear how they interact with host receptors. To identify how the S. sanguinis strain SK1 SRR adhesin affects interactions with sialylated glycans and glycoproteins, we determined high-resolution crystal structures of the binding domains alone and with purified trisaccharides. These structural studies determined that the ligands still bind at the noncanonical binding motif, but with fewer hydrogen-bonding interactions to the protein than is observed in structures of other Siglec-like adhesins. Complementary biochemical studies identified that each of the two binding domains has a different selectivity profile. Interestingly, the binding of SK1 to platelets and plasma glycoproteins identified that the interaction to some host targets is dominated by the contribution of one binding domain, whereas the binding to other host receptors is mediated by both binding domains. These results provide insight into outstanding questions concerning the roles of tandem domains in targeting host receptors and suggest mechanisms for how pathogens can adapt to the availability of a range of related but nonidentical host receptors. They further suggest that the definition of the YTRY motif should be changed to ΦTRX, a more rigorous description of this sialic acid–recognition motif given recent findings.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Quantifying metal-binding specificity of Cc NikZ-II from Clostridium carboxidivorans in the presence of competing metal ions

Many proteins bind transition metal ions as cofactors to carry out their biological functions. Despite binding affinities for divalent transition metal ions being predominantly dictated by the Irving-Williams series for wild-type proteins, in vivo metal ion binding specificity is ensured by intracellular mechanisms that regulate free metal ion concentrations. However, a growing area of biotechnology research considers the use of metal-binding proteins in vitro to purify specific metal ions from wastewater, where specificity is dictated by the protein's metal binding affinities. A goal of metalloprotein engineering is to modulate these affinities to improve a protein's specificity towards a particular metal; however, the quantitative relationship between the affinities and the equilibrium metal-bound protein fractions depends on the underlying binding mechanisms. Here we demonstrate a high-throughput intrinsic tryptophan fluorescence quenching method to validate binding models in multi-metal solutions for CcNikZ-II, a nickel-binding protein from Clostridium carboxidivorans. Using our validated models, we quantify the relationship between binding affinity and specificity in different classes of metal-binding models for CcNikZ-II. In conclusion, we further illustrate the potential relevance of data-informed models to predicting engineering targets for improved specificity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structures of the DarR transcription regulator reveal unique modes of second messenger and DNA binding

The mycobacterial repressor, DarR, a TetR family regulator (TFR), was the first transcription regulator shown to bind c-di-AMP. However, the molecular basis for this interaction and the mechanism involved in DNA binding by DarR remain unknown. Here we describe DarR-c-di-AMP and DarR-DNA structures and complementary biochemical assays. The DarR-c-di-AMP structure reveals a unique effector binding site for a TFR, located between DarR dimer subunits. Strikingly, we show this motif also binds cAMP. The location of the adenine nucleotide binding site between subunits suggests this interaction may facilitate dimerization and hence DNA binding. Indeed, biochemical assays show cAMP enhances DarR DNA binding. Finally, DarR-DNA structures reveal a distinct TFR DNA-binding mechanism involving two interacting dimers on the DNA. Thus, the combined data unveil a newly described second messenger binding motif and DNA binding mode for this important family of regulators.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Intra- and inter-molecular regulation by intrinsically-disordered regions governs PUF protein RNA binding

PUF proteins are characterized by globular RNA-binding domains. They also interact with partner proteins that modulate their RNA-binding activities. Caenorhabditis elegans PUF protein fem-3 binding factor-2 (FBF-2) partners with intrinsically disordered Lateral Signaling Target-1 (LST-1) to regulate target mRNAs in germline stem cells. Here, we report that an intrinsically disordered region (IDR) at the C-terminus of FBF-2 autoinhibits its RNA-binding affinity by increasing the off rate for RNA binding. Moreover, the FBF-2 C-terminal region interacts with its globular RNA-binding domain at the same site where LST-1 binds. This intramolecular interaction restrains an electronegative cluster of amino acid residues near the 5' end of the bound RNA to inhibit RNA binding. LST-1 binding in place of the FBF-2 C-terminus therefore releases autoinhibition and increases RNA-binding affinity. This regulatory mechanism, driven by IDRs, provides a biochemical and biophysical explanation for the interdependence of FBF-2 and LST-1 in germline stem cell self-renewal.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Affinity of Nicotinoids to a Model Nicotinic Acetylcholine Receptor (nAChR) Binding Pocket in the Human Brain

The binding affinity of nicotinoids to the binding residues of the a4ß2 variant of the nicotinic acetylcholine receptor (nAChR) was identified as a strong predictor of the nicotinoid’s addictive character. Using ab-iniito calculations for model binding pockets of increasing size comprising of 3, 6, and 14 amino acids (3AA, 6AA, and 14AA) that are derived from the crystal structure, the differences in binding affinity of 6 nicotinoids, namely nicotine (NIC), nornicotine (NOR), anabasine (ANB), anatabine (ANT), myosmine (MYO), and cotinine (COT) were correlated to their previously reported doses required for increases in intracranial self-stimulation (ICSS) thresholds, a metric for their addictive function. By employing the many body decomposition, the differences in the binding affinities of the various nicotinoids could be attributed mainly to the proton exchange energy between the Pyridine and non-Pyridine rings of the nicotinoids and the interactions between them and a handful of proximal amino acids, namely Trp156, Trpß57, Tyr100, and Tyr204. Interactions between the guest nicotinoid and the amino acids of the binding pocket were found to be mainly classical in nature, except for those between the nicotinoid and Trp156. The larger pockets were found to model binding structures more accurately and predicted the addictive character of all nicotinoids while smaller models, which are more computationally feasible, would only predict the addictive character of nicotinoids that are similar to nicotine. Here, the present study identifies the binding affinity of the guest nicotinoid to the host binding pocket as a strong descriptor of the nicotinoid’s addiction potential and as such it can be employed as a fast screening technique for the potential addiction of nicotine analogs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

DFT Analysis of the Binding of Rare Earth Nitrates at Internal and External Surfaces of MCM-22

MCM-22-type zeolites constitute a well-characterized tunable class of aluminosilicates suitable for elucidating the fundamental aspects of the binding of rare earth elements (REEs) in layered materials. Starting from the pure silica version ITQ-1 with a unit cell of Si 72 O 144 , a model for periodic bulk crystalline MCM-22 with a finite Al concentration is provided by replacing a Si atom with an Al atom in the unit cell at suitable tetrahedral sites near an internal pore surface. Then, a H atom is added to an O atom bridging Si and Al atoms to create a Brønsted acid site (BAS). There are no internal silanol groups in this bulk model. To generate a model for an external surface, we adopt the fully hydroxylated surface structure of a layer within the ITQ-1 precursor with two silanols per lateral unit cell. A BAS on the external surface can be generated by replacing a near-surface Si atom with an Al atom and adding a H atom, as above. The strength of binding at a BAS of REE, X, taken to be present in the solution phase as nitrates, is determined from the energy change in the reaction X(NO 3 ) 3 + ≡Si–{OH}–Al≡ → ≡Si–{OX(NO 3 ) 2 }–Al≡+ HNO 3 . The strength of binding at the silanols is determined similarly. Binding energies are determined from two approaches. The first performs periodic plane-wave density functional theory (DFT) total energy analysis for an entire unit cell of MCM-22. The second utilizes cluster models capturing the local environment of REE binding sites and performs DFT analysis with localized basis sets. The two approaches yield consistent results for Nd, revealing similarly strong binding at either an external or internal BAS, but much weaker binding at a silanol site. This is consistent with the picture deduced from recent experiments. Here, we also comment on binding at Al-bridged siloxane sites, which have been suggested as alternative binding sites to BAS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A novel assay to measure low-density lipoproteins binding to proteoglycans

Background: The binding of low-density lipoprotein (LDL) to proteoglycans (PGs) in the extracellular matrix (ECM) of the arterial intima is a key initial step in the development of atherosclerosis. Although many techniques have been developed to assess this binding, most of the methods are labor-intensive and technically challenging to standardize across research laboratories. Thus, sensitive, and reproducible assay to detect LDL binding to PGs is needed to screen clinical populations for atherosclerosis risk. Objectives: The aim of this study was to develop a quantitative, and reproducible assay to evaluate the affinity of LDL towards PGs and to replicate previously published results on LDL-PG binding. Methods: Immunofluorescence microscopy was performed to visualize the binding of LDL to PGs using mouse vascular smooth muscle (MOVAS) cells. An in-cell ELISA (ICE) was also developed and optimized to quantitatively measure LDL-PG binding using fixed MOVAS cells cultured in a 96-well format. Results: We used the ICE assay to show that, despite equal APOB concentrations, LDL isolated from adults with cardiovascular disease bound to PG to a greater extent than LDL isolated from adults without cardiovascular disease (p<0.05). Conclusion: We have developed an LDL-PG binding assay that is capable of detecting differences in PG binding affinities despite equal APOB concentrations. Future work will focus on candidate apolipoproteins that enhance or diminish this interaction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Binding of Sulfates and Water to Monovalent Cations

The binding of the sulfate ligand group to monovalent cations in the presence of water is important for many systems. To understand the structure and energetics of sulfate complexes, we use density functional theory to study ethyl sulfate binding to the monovalent cations Li + , Na + , and K + , and to water. The free energies of binding and optimal structures are calculated for a range of the number of ethyl sulfates and waters. Without water, the most optimal structure for all the cations is bidentate binding by two ethyl sulfates, yielding a 4-fold coordination. With water, the lowest free energy structures also have two ethyl sulfates, but the coordination varies with cations. For complexes with water, the four oxygen atoms in the sulfate group enable multiple binding geometries for the cations and for hydrogen bonding with water. Many of these geometries differ in free energy by only a small amount (1–2 kcal/mol), meaning there will be multiple binding configurations in bulk solution. In comparison to the optimal structures for binding to the carboxylate group, there is more variation for binding to the sulfate group as a function of cation type and the number of waters. Further, the polarization of the atoms is significant and varies among the sulfate oxygen atoms. The water oxygen charge is often larger than that of sulfate oxygen, which plays a role in the preference for monodentate ligand binding to cations in the presence of water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cyclic Peptides for Lanthanide Binding

Lanthanide ions are difficult to separate from one another due to their similar chemical properties. The discovery of lanthanide-binding peptides and proteins in nature has led to an increased interest in the possibility of utilizing the strong binding of peptides to lanthanide ions for their separations; as such, there has been an effort to identify or design peptides with improved lanthanide binding and selectivity toward particular lanthanide ions. Here, in this study, we designed and characterized lanthanide-binding cyclic peptides (LBCPs) with molecular dynamics simulations, electronic structure calculations, and emission spectroscopy. Luminescent decay measurements were done to determine the number of water molecules coordinated to the Eu 3+ ion in Eu-LBCP complexes and compare to the predicted number of water molecules by computation to assess the lanthanide-binding affinity of LBCPs. Measured stability constants show binding of the LBCPs to the Eu 3+ ion with stronger than micromolar affinity. We were able to identify multiple peptides that selectively bind to middle lanthanides. We describe the structural basis of the lanthanide-binding selectivity trend with strongest binding to the middle lanthanides, followed by the heavier lanthanides, and finally to the lighter ions.

ions↗
DOE OSTI2023

Inhibitor-3 inhibits Protein Phosphatase 1 via a metal binding dynamic protein–protein interaction

To achieve substrate specificity, protein phosphate 1 (PP1) forms holoenzymes with hundreds of regulatory and inhibitory proteins. Inhibitor-3 (I3) is an ancient inhibitor of PP1 with putative roles in PP1 maturation and the regulation of PP1 activity. Here, we show that I3 residues 27–68 are necessary and sufficient for PP1 binding and inhibition. In addition to a canonical RVxF motif, which is shared by nearly all PP1 regulators and inhibitors, and a non-canonical SILK motif, I3 also binds PP1 via multiple basic residues that bind directly in the PP1 acidic substrate binding groove, an interaction that provides a blueprint for how substrates bind this groove for dephosphorylation. Unexpectedly, this interaction positions a CCC (cys-cys-cys) motif to bind directly across the PP1 active site. Using biophysical and inhibition assays, we show that the I3 CCC motif binds and inhibits PP1 in an unexpected dynamic, fuzzy manner, via transient engagement of the PP1 active site metals. Together, these data not only provide fundamental insights into the mechanisms by which IDP protein regulators of PP1 achieve inhibition, but also shows that fuzzy interactions between IDPs and their folded binding partners, in addition to enhancing binding affinity, can also directly regulate enzyme activity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

CovS inactivation reduces CovR promoter binding at diverse virulence factor encoding genes in group A Streptococcus

The control of virulence gene regulator (CovR), also called caspsule synthesis regulator (CsrR), is critical to how the major human pathogen group A Streptococcus fine-tunes virulence factor production. CovR phosphorylation (CovR~P) levels are determined by its cognate sensor kinase CovS, and functional abrogating mutations in CovS can occur in invasive GAS isolates leading to hypervirulence. Presently, the mechanism of CovR-DNA binding specificity is unclear, and the impact of CovS inactivation on global CovR binding has not been assessed. Thus, we performed CovR chromatin immunoprecipitation sequencing (ChIP-seq) analysis in the emm1 strain MGAS2221 and its CovS kinase deficient derivative strain 2221-CovS-E281A. We identified that CovR bound in the promoter regions of nearly all virulence factor encoding genes in the CovR regulon. Additionally, direct CovR binding was observed for numerous genes encoding proteins involved in amino acid metabolism, but we found limited direct CovR binding to genes encoding other transcriptional regulators. The consensus sequence AATRANAAAARVABTAAA was present in the promoters of genes directly regulated by CovR, and mutations of highly conserved positions within this motif relieved CovR repression of the hasA and MGAS2221_0187 promoters. Analysis of strain 2221-CovS-E281A revealed that binding of CovR at repressed, but not activated, promoters is highly dependent on CovR~P state. CovR repressed virulence factor encoding genes could be grouped dependent on how CovR~P dependent variation in DNA binding correlated with gene transcript levels. Taken together, the data show that CovR repression of virulence factor encoding genes is primarily direct in nature, involves binding to a newly-identified DNA binding motif, and is relieved by CovS inactivation. These data provide new mechanistic insights into one of the most important bacterial virulence regulators and allow for subsequent focused investigations into how CovR-DNA interaction at directly controlled promoters impacts GAS pathogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Identification of preferred multimodal ligand‐binding regions on IgG1 F C using nuclear magnetic resonance and molecular dynamics simulations

Abstract In this study, the binding of multimodal chromatographic ligands to the IgG1 F C domain were studied using nuclear magnetic resonance and molecular dynamics simulations. Nuclear magnetic resonance experiments carried out with chromatographic ligands and a perdeuterated 15 N‐labeled F C domain indicated that while single‐mode ion exchange ligands interacted very weakly throughout the F C surface, multimodal ligands containing negatively charged and aromatic moieties interacted with specific clusters of residues with relatively high affinity, forming distinct binding regions on the F C . The multimodal ligand‐binding sites on the F C were concentrated in the hinge region and near the interface of the C H 2 and C H 3 domains. Furthermore, the multimodal binding sites were primarily composed of positively charged, polar, and aliphatic residues in these regions, with histidine residues exhibiting some of the strongest binding affinities with the multimodal ligand. Interestingly, comparison of protein surface property data with ligand interaction sites indicated that the patch analysis on F C corroborated molecular‐level binding information obtained from the nuclear magnetic resonance experiments. Finally, molecular dynamics simulation results were shown to be qualitatively consistent with the nuclear magnetic resonance results and to provide further insights into the binding mechanisms. An important contribution to multimodal ligand‐F C binding in these preferred regions was shown to be electrostatic interactions and π–π stacking of surface‐exposed histidines with the ligands. This combined biophysical and simulation approach has provided a deeper molecular‐level understanding of multimodal ligand–F C interactions and sets the stage for future analyses of even more complex biotherapeutics.

Gudhka, Ronak B.↗
DOE OSTI2021

Toward Guided Mutagenesis: Gaussian Process Regression Predicts MHC Class II Antigen Mutant Binding

Antigen-specific immunotherapies (ASI) require successful loading and presentation of antigen peptides into the major histocompatibility complex (MHC) binding cleft. One route of ASI design is to mutate native antigens for either stronger or weaker binding interaction to MHC. Exploring all possible mutations is costly both experimentally and computationally. To reduce experimental and computational expense, here we investigate the minimal amount of prior data required to accurately predict the relative binding affinity of point mutations for peptide-MHC class II (pMHCII) binding. Using data from different residue subsets, we interpolate pMHCII mutant binding affinities by Gaussian process (GP) regression of residue volume and hydrophobicity. We apply GP regression to an experimental data set from the Immune Epitope Database, and theoretical data sets from NetMHCIIpan and Free Energy Perturbation calculations. We find that GP regression can predict binding affinities of nine neutral residues from a six-residue subset with an average R 2 coefficient of determination value of 0.62 ± 0.04 (±95% CI), average error of 0.09 ± 0.01 kcal/mol (±95% CI), and with an receiver operating characteristic (ROC) AUC value of 0.92 for binary classification of enhanced or diminished binding affinity. Similarly, metrics increase to an R2 value of 0.69 ± 0.04, average error of 0.07 ± 0.01 kcal/mol, and an ROC AUC value of 0.94 for predicting seven neutral residues from an eight-residue subset. Our work finds that prediction is most accurate for neutral residues at anchor residue sites without register shift. This work holds relevance to predicting pMHCII binding and accelerating ASI design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗