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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

EF-Hand Battle Royale: Hetero-ion Complexation in Lanmodulin

The lanmodulin (LanM) protein has emerged as an effective means for rare earth element (REE) extraction and separation from complex feedstocks without the use of organic solvents. Whereas the binding of LanM to individual REEs has been well characterized, little is known about the thermodynamics of mixed metal binding complexes (i.e., heterogeneous ion complexes), which limits the ability to accurately predict separation performance for a given metal ion mixture. In this paper, we employ the law of mass action to establish a theory of perfect cooperativity for LanM-REE complexation at the two highest-affinity binding sites. The theory is then used to derive an equation that explains the nonintuitive REE binding behavior of LanM, where separation factors for binary pairs of ions vary widely based on the ratio of ions in the aqueous phase, a phenomenon that is distinct from single-ion-binding chemical chelators. We then experimentally validate this theory and perform the first quantitative characterization of LanM complexation with heterogeneous ion pairs using resin-immobilized LanM. Importantly, the resulting homogeneous and heterogeneous constants enable accurate prediction of the equilibrium state of LanM in the presence of mixtures of up to 10 REEs, confirming that the perfect cooperativity model is an accurate mechanistic description of REE complexation by LanM. We further employ the model to simulate separation performance over a range of homogeneous and heterogeneous binding constants, revealing important insights into how mixed binding differentially impacts REE separations based on the relative positioning of the ion pairs within the lanthanide series. In addition to informing REE separation process optimization, these results provide mathematical and experimental insight into competition dynamics in other ubiquitous and medically relevant, cooperative binding proteins, such as calmodulin.

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

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↗
DOE OSTI2022

Helium interaction with solutes and impurities in neutron-irradiated nanostructured ferritic alloys: A first principles study

Density functional theory calculations are performed to explore the binding between He and alloying solutes, impurities, and transmutation products expected in neutron irradiated nanostructured ferritic alloys (NFAs), here 14YWT is taken as an example. Elements that exhibit significant binding (attraction) with an interstitial He are Y (binding energy = 0.46 eV), Mg (0.32), O (0.33), Ti (0.16), and C (0.15). Those that provide significant binding to a substitutional He are O (1.44), Y (1.24), N (0.73), H (0.56), Mg (0.52), Ti (0.34), Si (0.34), C (0.33), Al (0.32), Ni (0.26), Ta (0.23), and Mn (0.16). The presence of these elements in Fe matrix could reduce the transport of He towards oxide particles, dislocations, and internal boundaries, and could promote He bubble nucleation in the matrix. For convenience, we compile existing binding energy data of He with He n and He n V (He-vacancy) clusters. Dissociation pathway analysis reveals that, in general, the most likely dissociation of a He n V cluster is by a sequential emission of individual He atoms. Furthermore, larger bubbles are more prone to dissociation than smaller ones. In addition, we estimate the binding energy (segregation energy) of He in bulk Y 2 Ti 2 O 7 (YTO) single crystal, YTO/Fe interface, and YTO particle embedded in Fe, with respect to interstitial He in Fe, from existing formation energies of He in these structures. We also compile available data of He binding with Fe self-interstitial atom (SIA), SIA clusters, and edge and screw dislocations. Note that given the absence of DFT data, the binding with SIA clusters and dislocations are gathered from simulations with empirical potentials. Finally, the data presented in this paper is important to inform multiscale simulations of He bubble accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bacterial hemophilin homologs and their specific type eleven secretor proteins have conserved roles in heme capture and are diversifying as a family

Cellular life relies on enzymes that require metals, which must be acquired from extracellular sources. Bacteria utilize surface and secreted proteins to acquire such valuable nutrients from their environment. These include the cargo proteins of the type eleven secretion system (T11SS), which have been connected to host specificity, metal homeostasis, and nutritional immunity evasion. This Sec-dependent, Gram-negative secretion system is encoded by organisms throughout the phylum Proteobacteria, including human pathogens Neisseria meningitidis, Proteus mirabilis, Acinetobacter baumannii, and Haemophilus influenzae. Experimentally verified T11SS-dependent cargo include transferrin-binding protein B (TbpB), the hemophilin homologs heme receptor protein C (HrpC), hemophilin A (HphA), the immune evasion protein factor-H binding protein (fHbp), and the host symbiosis factor nematode intestinal localization protein C (NilC). Here, we examined the specificity of T11SS systems for their cognate cargo proteins using taxonomically distributed homolog pairs of T11SS and hemophilin cargo and explored the ligand binding ability of those hemophilin cargo homologs. In vivo expression in Escherichia coli of hemophilin homologs revealed that each is secreted in a specific manner by its cognate T11SS protein. Sequence analysis and structural modeling suggest that all hemophilin homologs share an N-terminal ligand-binding domain with the same topology as the ligand-binding domains of the Haemophilus haemolyticus heme binding protein (Hpl) and HphA. We term this signature feature of this group of proteins the hemophilin ligand-binding domain. Network analysis of hemophilin homologs revealed five subclusters and representatives from four of these showed variable heme-binding activities, which, combined with sequence-structure variation, suggests that hemophilins are diversifying in function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparison of PsbQ and Psb27 in photosystem II provides insight into their roles

Photosystem II (PSII) catalyzes the oxidation of water at its active site that harbors a high-valent inorganic Mn 4 CaO x cluster called the oxygen-evolving complex (OEC). Extrinsic subunits generally serve to protect the OEC from reductants and stabilize the structure, but diversity in the extrinsic subunits exists between phototrophs. Recent cryo-electron microscopy experiments have provided new molecular structures of PSII with varied extrinsic subunits. We focus on the extrinsic subunit PsbQ, that binds to the mature PSII complex, and on Psb27, an extrinsic subunit involved in PSII biogenesis. PsbQ and Psb27 share a similar binding site and have a four-helix bundle tertiary structure, suggesting they are related. Here, we use sequence alignments, structural analyses, and binding simulations to compare PsbQ and Psb27 from different organisms. We find no evidence that PsbQ and Psb27 are related despite their similar structures and binding sites. Evolutionary divergence within PsbQ homologs from different lineages is high, probably due to their interactions with other extrinsic subunits that themselves exhibit vast diversity between lineages. This may result in functional variation as exemplified by large differences in their calculated binding energies. Psb27 homologs generally exhibit less divergence, which may be due to stronger evolutionary selection for certain residues that maintain its function during PSII biogenesis which is consistent with their more similar calculated binding energies between organisms. Previous experimental inconsistencies, low confidence binding simulations, and recent structural data suggest that Psb27 is likely to exhibit flexibility that may be an important characteristic of its activity. Furthermore, the analysis provides insight into the functions and evolution of PsbQ and Psb27, and an unusual example of proteins with similar tertiary structures and binding sites that probably serve different roles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

All-Atom Molecular Dynamics Simulations of Cationic Polyelectrolyte Brushes in the Presence of Halide Counterions

Understanding the response of the charged polyelectrolyte (PE) brushes and brush-supported water and ions to the changes in the nature of screening counterions is significant in developing strategies for utilizing PE brushes in various applications. In this paper, we employ all-atom molecular dynamics (MD) simulations for studying the behavior of cationic [poly(2-(methacryloyloxy)ethyl) trimethylammonium] (PMETA) brushes and brush-supported water and ions in the presence of different halide (X: I – , Br – , Cl – , and F – ) screening counterions. We find that despite the F– ion having the largest charge density, the extent of binding of the counterions on the PMETAX brushes varies as I – > Br – > Cl – > F – , leading to PMETAX brush height being least with I– counterions and greatest with F – counterions. This trend in the binding of the halide ions matches the previous experimental result and can be explained by identifying the chaotropic nature of the I – and Br – ions that promote a disruption of water structure and a more favorable binding of the ions to the polymer chains. Such a binding trend also ensures that the order of the water molecules around the PMETAX chains or the counterions as well as the number of water–water hydrogen bonds inside the brush layer increases in the following order of the counterion-specific PMETAX brushes: F – > Cl – > Br – > I – . Furthermore, halide-ion-PMETAX-chain binding takes place via both interchain and intrachain bridging: intrachain bridging dominates for the case of counterions that show enhanced binding (I – and Br – ), while interchain bridging is more favored for the case of counterions that show weakened binding (F – and Cl – ). Also, a greater degree of intrachain bridging leads to greater compressibility and flexibility of the brush layer. Lastly, we show that the mobility of the halide ions follows a nonmonotonic trend with the charge density: the mobility decreases as I – < Br – < Cl – (as their binding to the PMETAX chains varies as I – > Br – > Cl – ), but the mobility of F – ions is in between that of I – and Br – ions. We argue that the strongly attached hydration layer and the ensuing friction inside the brush layer lead to such a reduced mobility of the F – ions.

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

Differential roles of kinetic on- and off-rates in T-cell receptor signal integration revealed with a modified Fab’-DNA ligand

Antibody-derived T-cell receptor (TCR) agonists are commonly used to activate T cells. While antibodies can trigger TCRs regardless of clonotype, they bypass native T cell signal integration mechanisms that rely on monovalent, membrane-associated, and relatively weakly binding ligand in the context of cellular adhesion. Commonly used antibodies and their derivatives bind much more strongly than native peptide major histocompatibility complex (pMHC) ligands bind their cognate TCRs. Because ligand dwell time is a critical parameter that tightly correlates with physiological function of the TCR signaling system, there is a general need, both in research and therapeutics, for universal TCR ligands with controlled kinetic binding parameters. To this end, we have introduced point mutations into recombinantly expressed α-TCRβ H57 Fab to modulate the dwell time of monovalent Fab binding to TCR. When tethered to a supported lipid bilayer via DNA complementation, these monovalent Fab’-DNA ligands activate T cells with potencies well-correlated with their TCR binding dwell time. Single-molecule tracking studies in live T cells reveal that individual binding events between Fab'-DNA ligands and TCRs elicit local signaling responses closely resembling native pMHC. The unique combination of high on- and off-rates of the H57 R97L mutant enables direct observations of cooperative interplay between ligand binding and TCR-proximal condensation of the linker for activation of T cells, which is not readily visualized with pMHC. This work provides insights into how T cells integrate kinetic information from TCR ligands and introduces a method to develop affinity panels for polyclonal T cells, such as cells from a human patient.

Science & Technology - Other Topics↗
DOE OSTI2020

Recognition of nuclear export signals by CRM1 carrying the oncogenic E571K mutation

The E571K mutation of CRM1 is highly prevalent in some cancers, but its mechanism of tumorigenesis is unclear. Glu571 of CRM1 is located in its nuclear export signal (NES)-binding groove, suggesting that binding of select NESs may be altered. We generated HEK 293 cells with either monoallelic CRM1WT/E571K or biallelic CRM1E571K/E571K using CRISPR/Cas9. We also combined analysis of binding affinities and structures of 27 diverse NESs for wild-type and E571K CRM1 with structure-based bioinformatics. While most NESs bind the two CRM1 similarly, NESs from Mek1, eIF4E-transporter, and RPS2 showed >10-fold affinity differences. These NESs have multiple charged side chains binding close to CRM1 position 571, but this feature alone was not sufficient to predict different binding to CRM1(E571K). Consistent with eIF4E-transporter NES binding weaker to CRM1(E571K), eIF4E-transporter was mislocalized in tumor cells carrying CRM1(E571K). This serves as proof of concept that understanding how CRM1(E571K) affects NES binding provides a platform for identifying cargoes that are mislocalized in cancer upon CRM1 mutation. Finally, we showed that large affinity changes seen with some NES peptides (of Mek1 and RPS2) do not always translate to the full-length cargoes, suggesting limitations with current NES prediction methods. Therefore, comprehensive studies like ours are imperative to identify CRM1 cargoes with real pathogenic potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The Structure of Maltooctaose-Bound Escherichia coli Branching Enzyme Suggests a Mechanism for Donor Chain Specificity

Glycogen is the primary storage polysaccharide in bacteria and animals. It is a glucose polymer linked by α-1,4 glucose linkages and branched via α-1,6-linkages, with the latter reaction catalyzed by branching enzymes. Both the length and dispensation of these branches are critical in defining the structure, density, and relative bioavailability of the storage polysaccharide. Key to this is the specificity of branching enzymes because they define branch length. Herein, we report the crystal structure of the maltooctaose-bound branching enzyme from the enterobacteria E. coli. The structure identifies three new malto-oligosaccharide binding sites and confirms oligosaccharide binding in seven others, bringing the total number of oligosaccharide binding sites to twelve. In addition, the structure shows distinctly different binding in previously identified site I, with a substantially longer glucan chain ordered in the binding site. Using the donor oligosaccharide chain-bound Cyanothece branching enzyme structure as a guide, binding site I was identified as the likely binding surface for the extended donor chains that the E. coli branching enzyme is known to transfer. Furthermore, the structure suggests that analogous loops in branching enzymes from a diversity of organisms are responsible for branch chain length specificity. Together, these results suggest a possible mechanism for transfer chain specificity involving some of these surface binding sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

DFT-based investigation of solvation of Nd(III)/Yb(III) cations in room-temperature ionic-liquid

We report a DFT study on the solvation of bare rare-earth (RE) cations [Nd(III) and Yb(III)] in the room-temperature ionic-liquid (RTIL) based on 1,3-dimethyl-imidazolium + , PF 6 - , [MMI][PF 6 ]. In the present article, we aim to investigate whether the chosen RTIL can preferentially solvate any of the RE cations under study. First, the RTIL binding to the RE cation is studied in the gas phase, and then explored with continuum solvation model employing RTIL parameters. First, we considered the formation of the first solvation shell, where the PF 6 - anions directly bond with the Ln(III) cation. We then further surrounded the first solvation shell with [MMI]+ cations to form the second solvation shell. The stability of PF 6 - coordinated complexes is studied by calculating the binding energy values. The number of PF 6 - anions for both Ln(III) cations starts with three ligands and goes until the maximum binding energy in the RTIL medium is reached. The coordination number and binding modes of PF 6 - depend on the size of the RE cation. Comparing the binding energy values between the RE cations for the first solvation shell complexes, it is noticed that Yb-complexes have relatively greater binding energy than Nd-complexes. We have further performed a Boltzmann population analysis to find out the most dominant conformer among each category of PF 6 - coordinated complexes. A study using Natural Population Analysis shows that after complexation with PF 6 - the charge on the Yb-centers lower to a greater extent than Nd-centers, implying stronger bonding between Yb and PF 6 - , which is in harmony with the binding energy analysis. Overall, the computational investigation provides an in-depth understanding of the complexation phenomena of RE cations in RTIL, which can be useful for the practical separation of RE cations using RTIL.

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

Classification of Adsorbed Hydrocarbons Based on Bonding Configurations of the Adsorbates and Surface Site Stabilities

The design of heterogeneous catalysts can be accelerated by identifying relevant descriptors that accurately and effectively link the binding and activation energies to reactivity. Herein, we investigated scaling relations between binding energies of various hydrocarbon-based adsorbates on three different Pt surfaces and metal binding energies estimated via the recently developed α-scheme model. In this study, we find that the scaling slopes are similar for certain groups of adsorbates, which then can be classified based on their spatial and electronic structure enabling fast description of binding strengths for each member of the class. Hence, our findings show that the binding energies of simple hydrocarbons CH x , x = {0,1,2,3,4}, and CHCH 2 can be used to identify the binding energies of more complex hydrocarbon-based adsorbates. We introduce this classification to establish a generalizable scheme in which complex hydrogenation/dehydrogenation processes of higher hydrocarbons can be predicted via the binding energies of simpler hydrocarbon-based species and ultimately through surface site stabilities.

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

Differentiating Ligand Tailoring and Cation Incorporation as Strategies for Tuning Heterobimetallic Cerium Complexes

Tuning of redox-active complexes featuring metals with high coordination numbers by incorporation of secondary redox-inactive cations has received far less attention than it deserves. Here, appending moderate steric bulk to a tripodal ligand framework has been tested for its influence on secondary-cation-driven structural and electrochemical tuning of cerium, a lanthanide that tends to adopt high coordination numbers. A quasi -C 3 -symmetric cerium(III) complex denoted [Ce] has been prepared that features pendant benzyloxy groups, and this work demonstrates that this species offers a site capable of binding single Na + or Ca 2+ ions. Electrochemical and UV-visible spectroscopic studies reveal equilibrium binding affinity of [Ce] for Na + in acetonitrile solvent, contrasting with tight binding of all cations in all other previously studied systems of this type. The modulated cation binding can be attributed to the bulky benzyloxy groups, which impact the thermodynamics of cation binding but do not impede the formation of cerium centers with coordination number 8 upon binding of either Na + or Ca2 + . Finally, the Ce(IV/III) reduction potential was found to be tunable under the equilibrium binding conditions, highlighting the potentially significant role that controlled structural changes can play in modulating the solution properties of heterobimetallic complexes.

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

Effect of Ionic Strength and Anion Identity on Rare‐Earth Element Chelation by Poly(acrylic acid)

Abstract Rare‐earth elements (REEs) are critical components of many modern technologies, but their efficient separation and purification continue to pose significant challenges. Metal‐binding polymers may enable the development of promising alternative approaches to conventional separation techniques (e.g., solvent extraction), and insights into the interactions between the metal ions and polymers in solution are crucial to such development. In this work, the effect of ionic strength on the chelation of REEs by poly(acrylic acid) is investigated using isothermal titration calorimetry (ITC) to measure the metal‐binding thermodynamics. ITC experiments revealed that when ionic strength increases from increasing the sodium acetate buffer strength, polymer‐REE binding becomes increasingly disfavored. However, when ionic strength is increased by the addition of sodium chloride, binding again becomes increasingly disfavored, but to a lesser extent than when ionic strength increases from increasing the buffer strength. Additionally, the reduction in binding affinity is enthalpically driven when ionic strength is raised with sodium chloride, while it is entropically driven when raised with sodium acetate buffer, suggesting competive binding by the acetate ions. These results highlight the importance of ionic strength in polymer‐metal binding thermodynamics and will guide the development of new processes for the purification of critical metals.

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

Structural and biochemical analysis of the metallo-β-lactamase L1 from emerging pathogen Stenotrophomonas maltophilia revealed the subtle but distinct di-metal scaffold for catalytic activity

Emergence of Enterobacteriaceae harboring metallo-beta-lactamases (MBL) has raised global threats due to their broad antibiotic resistance profiles and the lack of effective inhibitors against them. We have been studied one of the emerging environmental MBL, the L1 from Stenotrophomonas maltophilia K279a . We determined several crystal structures of L1 complexes with three different classes of beta-lactam antibiotics (penicillin G, moxalactam, meropenem, and imipenem), with the inhibitor captopril and different metal ions (Zn +2 , Cd +2 , and Cu +2 ). All hydrolyzed antibiotics and the inhibitor were found binding to two Zn +2 ions mainly through the opened lactam ring and some hydrophobic interactions with the binding pocket atoms. Without a metal ion, the active site is very similarly maintained as that of the native form with two Zn +2 ions, however, the protein does not bind the substrate moxalactam. When two Zn +2 ions were replaced with other metal ions, the same di-metal scaffold was maintained and the added moxalactam was found hydrolyzed in the active site. Differential scanning fluorimetry and isothermal titration calorimetry were used to study thermodynamic properties of L1 MBL compared with New Deli Metallo-beta-lactamase-1 (NDM-1). Both enzymes are significantly stabilized by Zn +2 and other divalent metals but showed different dependency. These studies also suggest that moxalactam and its hydrolyzed form may bind and dissociate with different kinetic modes with or without Zn +2 for each of L1 and NDM-1. Our analysis implicates metal ions, in forming a distinct di-metal scaffold, which is central to the enzyme stability, promiscuous substrate binding and versatile catalytic activity. Statement: The L1 metallo-beta-lactamase from an environmental multidrug-resistant opportunistic pathogen Stenotrophomonas maltophilia K279a has been studied by determining 3D structures of L1 enzyme in the complexes with several beta-lactam antibiotics and different divalent metals and characterizing its biochemical and ligand binding properties. We found that the two-metal center in the active site is critical in the enzymatic process including antibiotics recognition and binding, which explains the enzyme's activity toward diverse antibiotic substrates. This study provides the critical information for understanding the ligand recognition and for advanced drug development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Integrated structural model of the palladin–actin complex using XL ‐ MS , docking, NMR , and SAXS

Abstract Palladin is an actin‐binding protein that accelerates actin polymerization and is linked to the metastasis of several types of cancer. Previously, three lysine residues in an immunoglobulin‐like domain of palladin have been identified as essential for actin binding. However, it is still unknown where palladin binds to F‐actin. Evidence that palladin binds to the sides of actin filaments to facilitate branching is supported by our previous study showing that palladin was able to compensate for Arp2/3 in the formation of Listeria actin comet tails. Here, we used chemical crosslinking to covalently link palladin and F‐actin residues based on spatial proximity. Samples were then enzymatically digested, separated by liquid chromatography, and analyzed by tandem mass spectrometry. Peptides containing the crosslinks and specific residues involved were then identified for input to the HADDOCK docking server to model the most likely binding conformation. Small‐angle x‐ray scattering was used to provide further insight into palladin flexibility and the binding interface, and NMR spectra identified potential interactions between palladin's Ig domains. Our final structural model of the F‐actin:palladin complex revealed how palladin interacts with and stabilizes F‐actin at the interface between two actin monomers. Three actin residues that were identified in this study also appear commonly in the actin‐binding interface with other proteins such as myotilin, myosin, and tropomodulin. An accurate structural representation of the complex between palladin and actin extends our understanding of palladin's role in promoting cancer metastasis through the regulation of actin dynamics.

Sargent, Rachel [Department of Chemistry and Bioch↗
DOE OSTI2025

Uncovering Sequence and Structural Characteristics of Fungal Expansin‐Related Proteins With Potential to Drive Substrate Targeting

Expansins loosen plant cell wall networks through disrupting non-covalent bonds between cellulose microfibrils and matrix polysaccharides. Whereas expansins were first discovered in plants, expansin-related proteins have since been identified in bacteria and fungi. The biological function of microbial expansins remains unclear; however, several studies have shown distinct binding preferences toward different structural polysaccharides. Earlier studies of bacterial expansin-related proteins uncovered sequence and structural features that correlate to substrate binding. Herein, 20 fungal expansin-related sequences were recombinantly produced in Komagataella phaffii, and the purified proteins were compared in terms of substrate binding to cellulosic and chitinous substrates. The impact of pH on the zeta potential of prioritized substrates was also measured, and Principal Component Analysis was performed to uncover correlations between protein characteristics (e.g., pI, hydrophobicity, surface charge distribution) and measured substrate binding preferences. Whereas acidic proteins with a predicted pI less than 5.0 preferentially bound to chitin, basic proteins with pI greater than 8.0 preferentially bound to xylan and xylan-containing fiber. Similar to many cellulases, binding to cellulose was correlated to relatively high aromatic amino acid content in the protein sequence and presence of a carbohydrate binding module (CBM), which in the case of expansins is a C-terminal CBM63. Whereas overall sequence characteristics could be correlated to substrate binding preference, the identity of amino acids occupying conserved positions that impact protein activity was better correlated with loosenin versus expansin classifications.

chitin↗
DOE OSTI2025

Ion-selective conformational stabilization of a disordered repeats-in-toxin protein domain

Ion-binding intrinsically disordered proteins (IDPs) recruit and bind to specific metal ions to perform critical biological functions. In proteins where ion binding and structural transitions are coupled, interactions with off-target toxic metals can dramatically disrupt protein structure and function, exemplified by lead and mercury poisoning. Understanding the complex mechanisms underlying how IDPs exclude or allow binding to different ionic species is crucial for addressing the origins of metal toxicity in biological systems. Here, we elucidate mechanisms of ion selectivity in an IDP that adopts a structure upon Ca 2+ binding. We probed ion-induced conformational changes of a repeats-in-toxin (RTX) protein domain in the presence of different ion ligands—Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ —with chemical similarities but drastically different ionic radii. RTX adopts ion-selective conformations measured by x-ray crystallography, small-angle x-ray scattering, and circular dichroism. High-resolution x-ray structures reveal that Sr 2+ induces a nearly identical RTX structure as natively binding Ca 2+ , enabled by the intrinsic flexibility and disorder of the protein. Small-angle x-ray scattering and circular dichroism indicate that smaller Mg 2+ does not induce a significant conformational change in RTX, whereas larger Ba 2+ induces a partially folded structure. These results highlight the importance of geometric constraints imposed by protein structure in determining metal ion selectivity, yielding insights into how off-target ion binding may result in protein misfolding and malfunction.

Gudinas, Alana P. [Stanford Univ., CA (United Stat↗