Computational Antigen Optimization through Symbolic Optimization and Affinity Maturation Simulation
Explore the source record for details and available documents.
Biological and medical sciences, Mathematics and C↗
SEARCH · Engineering Papers
Results for “affinity”
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.
Strategies to Scale Up Antibody Affinity Assays
Explore the source record for details and available documents.
Biological and medical sciences↗
Strategies to Scale Up Antibody Affinity Assays
Explore the source record for details and available documents.
Biological and medical sciences↗
Base-Level Optimization of Library Design Facilitates Cost-Effective Scale Up of Antibody Affinity Assays
Explore the source record for details and available documents.
Biological and medical sciences↗
Rational Design of Lanmodulin Variants for Size-Based Selectivity of Individual Rare Earth Elements
Rare earth elements (REEs) are essential to modern technologies, yet their high physical and chemical similarity makes separation of individual REEs difficult and environmentally taxing. Metalloproteins offer a promising alternative for selective REE binding, as they tend to have high metal ion affinity and specificity. Lanmodulin (LanM), in particular, has arisen as a potential candidate for REE separation as it exhibits picomolar affinity for elements in the REE family. Prior work has shown that the single point mutation D9N can shift LanM’s preference away from lanthanides toward actinides, motivating efforts to tune selectivity of LanM through targeted mutagenesis. Here, we tested the hypothesis that introducing selective aspartic acid to glutamic acid substitutions in the metal coordinating EF hands of LanM would impose steric constraints that would drive LanM affinity away from larger ions, such as La3+, to smaller ions, such as Y3+. To test this hypothesis, a combination of computational and experimental approaches were employed to evaluate the signal mutations LanM D5E and LanM D3E and the double mutants LanM D1ED5E and LanM D3ED9E. Surprisingly, increasing the number of mutations within the metal center did not enhance affinity for smaller REEs, or decrease affinity for larger ions. Only the single point mutation LanM D5E weakened La3+ binding by one order of magnitude relative to LanM wild type (WT), and pairing it with a second mutation to produce LanM D1ED5E drove La3+ affinity to be stronger than that seen for LanM WT. The D3E mutation alone prevented proper expression and folding, but paring it with D9E to produce LanM D3ED9E rescued expression and yielded La3+ affinities comparable to LanM WT. All variants that expressed (LanM D5E, LanM D1ED5E, LanM D3ED9E) displayed Y3+ affinities comparable to LanM WT. Overall, these results highlight the tunability of LanM’s metal-binding environment but also expose current limitations in predicting structural responses to point mutations within a protein sequence. This work establishes a foundation that can be used for refining computational and experimental strategies to engineer metalloproteins with tailored REE selectivity.
Close, Emily [Pacific Northwest National Laborator↗
Author Correction: Genome-guided isolation of the hyperthermophilic aerobe Fervidibacter sacchari reveals conserved polysaccharide metabolism in the Armatimonadota
Correction to: Nature Communicationshttps://doi.org/10.1038/s41467-024-53784-3, published online 4 November 2024 In the version of this article initially published, Table 1 did not include the properties of the taxa being proposed or refer directly to another location in the main manuscript describing the properties. As such, the original manuscript did not comply with Rule 27 (2)(c) of the ICNP. Also, Table 1 listed the order Fervidibacterales as the nomenclatural type for the class Fervidibacteria, which violates latest emended version of Rule 15 stating that the nomenclatural type for a class must be a genus. Below we provide a modification of Table 1 containing protologues with these errors corrected. We have also changed the order of the taxa in the table to meet the most common ordering. (Table presented.) Taxon names proposed under the ICNP Proposed taxon Etymology Description Genus Fervidibacter Fer.vi.di.bac’ter. L. masc. adj. fervidus, hot, steaming; N.L. masc. n. bacter, a rod; N.L. masc. n. Fervidibacter, a hot rod Thermophilic or hyperthermophilic inhabitants of freshwater thermal environments. All members are likely polysaccharide-degrading chemoheterotrophs with numerous carbohydrate-active enzymes encoded in their genomes. Aerobic, with high-affinity and/or low-affinity terminal oxidases present in the genomes. The oxidative pentose phosphate pathway and the tricarboxylic acid cycle are complete in genomes belonging to the genus. Gram-stain-negative and diderm cell envelope structure. Ovoid- to rod-shaped morphology. Spores are not formed. The genus is a distinct phylogenetic lineage in the family Fervidibacteraceae, the order Fervidibacterales, and the class Fervidibacteria in the phylum Armatimonadota. The type species is Fervidibacter sacchariT. Species Fervidibacter sacchari sac’cha.ri. N.L. gen. n. sacchari, of sugar Hyperthermophilic, microaerophilic, facultatively anaerobic, and grows chemoheterotrophically on monosaccharides and polysaccharides. Cells are ovoid- to rod-shaped, Gram-stain negative, and are 0.9–1.3 µm in width and 1.6–3.6 µm in length. Grows between 65 and 87.5 °C and an optimum temperature of 80 °C, and a pH range of 6.5–8.6 with an optimum pH of 7.5. Grows at an optimum O2 concentration of 5–10%. Grows on D-arabinose, D-galactose, D-glucose, D-rhamnose, D-ribose, D-xylose, chondroitin sulfate, colloidal chitin, galactan, gellan gum, guar gum, karaya gum, locust bean gum, xantham gum, xyloglucan, β-glucan, glycogen, starch, AFEX-pretreated corn stover, miscanthus, sugarcane bagasse, acetate and casamino acids. Grows weakly on xyloglucan under fermentation conditions. The major fatty acids (>10%) are C16:0, C18:0 and/or cyclo-C17:0, and iso-C16:0. The major respiratory quinones (>10%) are MK-8 and MK-9. The isolate and genomes of the species have been recovered from geothermal springs in the Great Basin, Nevada, USA. GC content of genomes range between 51–52%. Subunits for both the high-affinity and low-affinity terminal oxidases are encoded in the genomes. Genomes also encode a Group 3d [NiFe] hydrogenase, which produces hydrogen as an electron sink for NAD+ regeneration. The type strain PD1T (= JCM 39283T = DSM 113467T) was isolated from Great Boiling Spring in Nevada, USA. Family Fervidibacteraceae Fer.vi.di.bac.te.ra’ce.ae. N.L. masc. n. Fervidibacter type genus of the family; L. suff. -aceae ending to denote a family; N.L. fem. pl. n. Fervidibacteraceae the family of the genus Fervidibacter Thermophilic or hyperthermophilic inhabitants of freshwater thermal environments. All members are likely polysaccharide-degrading chemoheterotrophs with numerous carbohydrate-active enzymes encoded in their genomes. Aerobic, with high-affinity and/or low-affinity terminal oxidases present in the genomes. The oxidative pentose phosphate pathway and the tricarboxylic acid cycle are complete in genomes belonging to the family. The family is a distinct phylogenetic lineage in the order Fervidibacterales and the class Fervidibacteria in the phylum Armatimonadota. The type genus is Fervidibacter. Order Fervidibacterales Fer.vi.di.bac.te.ra’les. N.L. masc. n. Fervidibacter type genus of the order; L. suff. -ales ending to denote an order; N.L. fem. pl. n. Fervidibacterales the order of the genus Fervidibacter Thermophilic or hyperthermophilic inhabitants of freshwater thermal environments. All members are likely polysaccharide-degrading chemoheterotrophs with numerous carbohydrate-active enzymes encoded in their genomes. Aerobic or strictly anaerobic. Phylogenomic placement of this lineage within the Fervidibacteria and relative evolutionary divergence supports delineation of this lineage as an order within the class Fervidibacteria and phylum Armatimonadota. The type genus is Fervidibacter. Class Fervidibacteria Fer.vi.di.bac.te’ri.a. N.L. masc. n. Fervidibacter type genus of the type order of the class; L. suff. -ia ending to denote a class; N.L. neut. pl. n. Fervidibacteria the class of the order Fervidibacterales Thermophilic or hyperthermophilic inhabitants of freshwater thermal environments. All members are likely polysaccharide-degrading chemoheterotrophs with numerous carbohydrate-active enzymes encoded in their genomes. Aerobic or strictly anaerobic. Phylogenomic placement of this lineage within the Armatimonadota and relative evolutionary divergence supports delineation of this lineage as a class within the Armatimonadota. The type genus is Fervidibacter. The error has not been corrected in the PDF or HTML versions of the Article.
Nou, Nancy O↗
Machine-learning-informed scattering correlation analysis of sheared colloids
We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.
Gaussian process regression↗
Cross-species analysis of FcγRIIa/b (CD32a/b) polymorphisms at position 131: structural and functional insights into the mechanism of IgG- mediated phagocytosis in human and macaque
Introduction Antibodies play a critical role in immunity in part by mediating clearance of pathogens and infected cells by antibody-dependent cellular phagocytosis (ADCP) through engagement of Fc gamma receptors (FcγRs) on innate immune cells. Among these, FcγRIIa (CD32a) is a key activating receptor expressed on macrophages, dendritic cells, and other antigen-presenting cells. Its affinity for IgG and ability to mediate ADCP is influenced by allelic polymorphisms. In humans, a single amino acid polymorphism at position 131, where histidine (H) is substituted with arginine (R), leads to decreased IgG1 and IgG2 subclass binding affinity and, consequently, lower efficiency of phagocytic responses. Rhesus macaques ( Macaca mulatta ), which are widely used as nonhuman primate models, exhibit a similar polymorphism at position 131 of FcγRIIa, but with arginine replaced by proline (P). Here, we investigated structure-function relationships associated with the FcγRIIa polymorphism at position 131 in both species, specifically with respect to IgG1 and IgG2. Methods We determined the structures of complexes formed by each variant with IgG1 Fc and those formed by the higher affinity variant with IgG2 Fc for both species by x-ray crystallography and linked these structures to affinity and activity using SPR and an ADCP assay. We also determined the structure of human inhibitory FcγRIIb (CD32b) in complex with IgG1 Fc by x-ray crystallography. Results Through analysis of these structures, our studies reveal that FcγRIIa engagement is minimally influenced by Fc glycan composition, distinguishing it from FcγRIIIa whose affinity is strongly influenced by glycan-composition. Comparative structures of human and macaque FcγRIIa variants demonstrate species- and allele-specific differences in Fc binding, but our functional assays showed only minimal allele-specific effects in humans. In contrast, allele-specific effects in macaques were highly significant; the macaque P 131 variant showing uniformly reduced IgG affinity. Conclusion These insights highlight fundamental interspecies and allelic distinctions that are critical for interpreting FcγRIIa-mediated effector functions in macaque models and for optimizing translational antibody and vaccine design.
Tolbert, William D.↗
Mechanistic Similarities and Differences for Hydrogenation of Aromatic Heterocycles and Aliphatic Carbonyls on Sulfided Ru Nanoparticles
This study establishes the contrasting reactivity trends for the hydrogenation of aromatic hydrocarbons (AHCs = C n H m X, X = N, S, O, and C) and aliphatic carbonyls [R a C(O)R b ; R a = alkyl group, R b = CH 3 or H] on sulfided Ru clusters arising from the difference in the elementary proton and hydride attack catalytic sequence. Both reactions require sequential additions of a proton from either Ru n+ –(SH 2 ) or S 2 –(H δ+ ) species and a hydride from Ru n+ –(H δ– ) species to the unsaturated C=X bonds. For the five-membered-ring aromatic heterocycles (AHCs = pyrrole, thiophene, and furan), an initial proton addition limits the catalytic turnovers; thus, their hydrogenation reactivity increases with increasing gas-phase proton affinities of the AHCs. Pyridine as the more basic six-membered N-AHC is more susceptible to protonation; therefore, it is more reactive, and its initial proton addition is quasi-equilibrated, followed by the kinetically relevant hydride addition. Conversely, aliphatic carbonyls prefer to undergo hydrogen additions in a reverse sequence, where a Ru n+ –(H δ– ) hydride initially attacks the electron-deficient carbonyl C atom as the kinetically relevant step before a subsequent rapid S 2 –(H δ+ ) proton addition on the electron-rich O atom, as confirmed by isotopic exchange studies with butanal-D 2 and 1-butanol-D 2 mixtures and density functional theory calculations on S-deficient RuS 2 (100) surfaces. For these reasons, their hydrogenation reactivity increases with increasing gas-phase hydride affinities of the carbonyls. On metallic Ru surfaces without sulfur, hydrogen adatoms (H*) are the only reactive hydrogen species; the reactivity of their attack on aromatic heterocycles increases with increasing reactant proton affinities much more sensitively than that on sulfided Ru surfaces. This work illustrates the distinct catalytic roles of the diverse hydrogen species in hydrogenation—the interplay between proton and hydride additions has marked catalytic consequences in shaping the free energy landscape of the reactions, which in turn leads to the observed kinetic dependences, kinetic parameters, reactivity trends, and scaling relations between the measured barriers and the appropriate kinetic descriptors, that is, proton affinity of AHCs and hydride affinity of carbonyls in hydrogenation catalysis. Here, these mechanistic similarities and differences provide explanations of the observed reactivity trends and thus have profound implications for industrial hydrogenation and hydrotreating catalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
Probing surface-adsorbate interactions through active particle dynamics
Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structural basis for breadth development in the HIV-1 V3-glycan targeting DH270 antibody clonal lineage
Antibody affinity maturation enables adaptive immune responses to a wide range of pathogens. In some individuals broadly neutralizing antibodies develop to recognize rapidly mutating pathogens with extensive sequence diversity. Vaccine design for pathogens such as HIV-1 and influenza has therefore focused on recapitulating the natural affinity maturation process. Here, we determine structures of antibodies in complex with HIV-1 Envelope for all observed members and ancestral states of the broadly neutralizing HIV-1 V3-glycan targeting DH270 antibody clonal B cell lineage. These structures track the development of neutralization breadth from the unmutated common ancestor and define affinity maturation at high spatial resolution. By elucidating contacts mediated by key mutations at different stages of antibody development we identified sites on the epitope-paratope interface that are the focus of affinity optimization. Thus, our results identify bottlenecks on the path to natural affinity maturation and reveal solutions for these that will inform immunogen design aimed at eliciting a broadly neutralizing immune response by vaccination.
60 APPLIED LIFE SCIENCES↗
Structural and kinetic characterization of an acetoacetyl-Coenzyme A: acetate Coenzyme A transferase from the extreme thermophile Thermosipho melanesiensis
Family 1 Coenzyme A transferases (CtfAB) from the extremely thermophilic bacterium, Thermosipho melanesiensis, has been used for in vivo acetone production up to 70°C. This enzyme has tentatively been identified as the rate-limiting step, due to its relatively low-binding affinity for acetate. However, existing kinetic and mechanistic studies on this enzyme are insufficient to evaluate this hypothesis. Here, kinetic analysis of purified recombinant T. melanesiensis CtfAB showed that it has a ping-pong bi-bi mechanism typical of Coenzyme A (CoA) transferases with Km values for acetate and acetoacetyl-CoA of 85 mM and 135 μM, respectively. Product inhibition by acetyl-CoA was competitive with respect to acetoacetyl-CoA and non-competitive with respect to acetate. Crystal structures of wild-type and mutant T. melanesiensis CtfAB were solved in the presence of acetate and in the presence or absence of acetyl-CoA. These structures led to a proposed structural basis for the competitive and non-competitive inhibition of acetyl-CoA: acetate binds independently of acetyl-CoA in an apparent low-affinity binding pocket in CtfA that is directly adjacent to a catalytic glutamate in CtfB. Similar to other CoA transferases, acetyl-CoA is bound in an apparent high-affinity binding site in CtfB with most interactions occurring between the phospho-ADP of CoA and CtfB residues far from the acetate binding pocket. This structural-based mechanism also explains the organic acid promiscuity of CtfAB. High-affinity interactions are predominantly between the conserved phospho-ADP of CoA, and the variable organic acid binding site is a low-affinity binding site with few specific interactions.
Biochemistry & Molecular Biology↗
The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy
An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Modeling receptor flexibility in the structure-based design of KRAS G12C inhibitors
KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS G12C —sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS G12C . In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS G12C , we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.
59 BASIC BIOLOGICAL SCIENCES↗
A machine learning approach for identifying variables associated with risk of developing neutralizing antidrug antibodies to factor VIII
A key unmet need in the management of hemophilia A (HA) is the lack of clinically validated markers that are associated with the development of neutralizing antibodies to Factor VIII (FVIII) (commonly referred to as inhibitors). This study aimed to identify relevant biomarkers for FVIII inhibition using Machine Learning (ML) and Explainable AI (XAI) using the My Life Our Future (MLOF) research repository. The dataset includes biologically relevant variables such as age, race, sex, ethnicity, and the variants in the F8 gene. In addition, we previously carried out Human Leukocyte Antigen Class II (HLA-II) typing on samples obtained from the MLOF repository. Using this information, we derived other patient-specific biologically and genetically important variables. These included identifying the number of foreign FVIII derived peptides, based on the alignment of the endogenous FVIII and infused drug sequences, and the foreign-peptide HLA-II molecule binding affinity calculated using NetMHCIIpan. The data were processed and trained with multiple ML classification models to identify the top performing models. The top performing model was then chosen to apply XAI via SHAP, (SHapley Additive exPlanations) to identify the variables critical for the prediction of FVIII inhibitor development in a hemophilia A patient. Using XAI we provide a robust and ranked identification of variables that could be predictive for developing inhibitors to FVIII drugs in hemophilia A patients. These variables could be validated as biomarkers and used in making clinical decisions and during drug development. The top five variables for predicting inhibitor development based on SHAP values are: (i) the baseline activity of the FVIII protein, (ii) mean affinity of all foreign peptides for HLA DRB 3, 4, & 5 alleles, (iii) mean affinity of all foreign peptides for HLA DRB1 alleles), (iv) the minimum affinity among all foreign peptides for HLA DRB1 alleles, and (v) F8 mutation type.
60 APPLIED LIFE SCIENCES↗
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↗
Effects of SARS-CoV-2 Main Protease Mutations at Positions L50, E166, and L167 Rendering Resistance to Covalent and Noncovalent Inhibitors
SARS-CoV-2 propagation under nirmatrelvir and ensitrelvir pressure selects for main protease (MPro) drug-resistant mutations E166V (DRM2), L50F/E166V (DRM3), E166A/L167F (DRM4), and L50F/E166A/L167F (DRM5). DRM2-DRM5 undergoes N-terminal autoprocessing to produce mature MPro with dimer dissociation constants (K dimer ) 2–3 times larger than that of the wildtype. Co-selection of L50F restores catalytic activity of DRM2 and DRM4 from ~10 to 30%, relative to that of the wild-type enzyme, without altering K dimer . Binding affinities and thermodynamic profiles that parallel the drug selection pressure, exhibiting significant decreases in affinity through entropy/enthalpy compensation, were compared with GC373. Reorganization of the active sites due to mutations observed in the inhibitor-free DRM3 and DRM4 structures as compared to MPro WT may account for the reduced binding affinities, although DRM2 and DRM3 complexes with ensitrelvir are almost identical to MPro WT -ensitrelvir. In conclusion, chemical reactivity changes of the mutant active sites due to differences in electrostatic and protein dynamics effects likely contribute to losses in binding affinities.
60 APPLIED LIFE SCIENCES↗