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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Engineered selective biotoxin‐binding hydrogels for toxin sequestration

The development of synthetic selective membranes that separate materials of similar sizes, charges, and/or polarities remains a difficult challenge, and looking towards biology provides inspiration for new designs. In this work, a series of cholera toxin binding peptides (CTBPs) are identified, spanning a range of binding inhibitions, and integrated into chemically cross-linked cholera toxin binding gels (CTBGs) via thiol-Michael polycondensation reactions. All gels demonstrate rheological profiles consistent with elastic solids. The CTBGs are probed via small-angle neutron scattering and exhibit a correlation length, ξ, smaller than most proteins (1.3–2.5 nm). Thus, an effective entropic mesh is formed to block non-targeted proteins. However, the CTBGs have a dynamic mesh size, Ξ, that is larger than cholera toxin (CT) to allow the transport of target proteins. The CTBGs with the highest binding inhibitions both show high selectivity and permeation of CT, rejecting all other tested proteins. In total, two new highly selective CTBGs are synthesized and validated for use in cholera toxin remediation. Together, this platform demonstrates the wide applicability of selectively-diffusive materials for difficult separations.

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

Determination of affinities of lanthanide-binding proteins using chelator-buffered titrations

The recent discoveries of the first proteins that bind lanthanides as part of their biological function not only are relevant to the emerging field of lanthanide-dependent biology, but also hold promise to revolutionize the technologically critical rare earths industry. Although protocols to assess the thermodynamics of metal–protein interactions are well established for “traditional” metal ions in biology, the characterization of lanthanide-binding proteins presents a challenge to biochemists due to the lanthanides' Lewis acidity, propensity for hydrolysis, and high-affinity complexes with biological ligands. These properties necessitate the preparation of metal stock solutions with very low buffered “free” metal concentrations (e.g., femtomolar to nanomolar) for such determinations. Here in this paper we describe several protocols to overcome these challenges. First, we present standardization methods for the preparation of chelator-buffered solutions of lanthanide ions with easily calculated free metal concentrations. We also describe how these solutions can be used in concert with analytical methods including UV–visible spectrophotometry, circular dichroism spectroscopy, Förster resonance energy transfer (FRET), and sensitized terbium luminescence, in order to accurately determine dissociation constants (K d s) of lanthanide–protein complexes. Finally, we highlight how application of these methods to lanthanide-binding proteins, such as lanmodulin, has yielded insights into selective recognition of lanthanides in biology. We anticipate that these protocols will facilitate discovery and characterization of additional native lanthanide-binding proteins, will motivate the understanding of their biological context, and will prompt their applications in biotechnology.

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

A Trichomonas vaginalis C2-XYPPX-repeat protein with a structured C2 domain displaying dampened flexibility upon binding calcium

C2 domains are ubiquitous membrane-binding modules of ∼130 residues in eukaryotes that are often associated with proteins involved in membrane trafficking and lipid modification. The genome of Trichomonas vaginalis, the most common, non-viral, sexually transmitted human pathogen, encodes eight genes that contain a N-terminal C2 module linked to a XYPPX-repeat domain of more than four XYPPX repeats (C2-XYPPX). While the function of the XYPPX-repeat domain remains unknown, its multiple association with C2 domains in T. vaginalis suggests it is important. Here, the C2 domain from one of these C2-XYPPX-repeat proteins, Tv-C2-1, was structurally and physically characterized using X-ray crystallography and NMR spectroscopy. The crystal structure for Tv-C2-1 shows that this domain shares a fold common to all C2 domains, a compact Greek-key motif composed of eight anti-parallel β-strands in the type-2 topology. An NMR chemical shift perturbation study with Ca 2+ showed that Tv-C2-1 bound two Ca 2+ atoms primarily via two loops (loop-1 and loop-3) on the predicted calcium binding face of the protein with K d s of 58.0 ± 0.1 μM and 232 ± 6 μM. Estimations of the overall rotational correlation time, τ c , in the apo (11.1 ns) and Ca 2+ -bound (9.2 ns) state suggests the protein becomes more compact upon Ca 2+ binding, consistent with a decrease in dynamics in loop-3 and marginally in loop-1 suggested by amide 15 N heteronuclear steady-state { 1 H}- 15 N NOEs. Showing Tv-C2-1 binds calcium and adopts a compact Greek-key motif structure, two primary features of C2 domains, suggests understanding the function of the XYPPX-repeat domain may be warranted.

NMR spectroscopy↗
DOE OSTI2021

Reactions of self-healing agents and the chemical binding of aggressive ions in sea water: Thermodynamics and kinetics

New self-healing agents that can chemically bind seawater ions invading cracked cementitious materials were proposed. The potential of self-healing and binding of seawater ions were investigated by thermodynamic modeling. It was found that CaO-NaAlO{sub 2} and CaO-metakaolin agents can have Cl{sup −}, SO{sub 4}{sup 2−} and Mg{sup 2+} chemically bound by reacting with sea water to form Friedel's salt, Kuzel's salt, ettringite and hydrotalcite. The removal of Cl{sup −} from seawater firstly increased and then decreased with the increase of Ca/Al molar ratio in both agents, while the removal of Mg{sup 2+} and SO{sub 4}{sup 2−} were hardly influenced and approximated 100%. Because NaAlO{sub 2} dissolves and releases Al(OH){sub 4}{sup −} rapidly, precipitates binding Cl{sup −}, SO{sub 4}{sup 2−} and Mg{sup 2+} were formed fast. In comparison, the reaction of metakaolin binding aggressive ions occurred after 3 days. Because of the faster reaction and the capacity to make [Cl{sup −}]/[OH{sup −}] lower in the solution, CaO-NaAlO{sub 2} would be more efficient for self-healing and mitigating reinforcement corrossion than CaO-metakaolin.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neutralizing monoclonal antibodies elicited by mosaic RBD nanoparticles bind conserved sarbecovirus epitopes

Increased immune evasion by SARS-CoV-2 variants of concern highlights the need for new therapeutic neutralizing antibodies. Immunization with nanoparticles co-displaying spike receptor-binding domains (RBDs) from eight sarbecoviruses (mosaic-8 RBD-nanoparticles) efficiently elicits cross-reactive polyclonal antibodies against conserved sarbecovirus RBD epitopes. Here, we identified monoclonal antibodies (mAbs) capable of cross-reactive binding and neutralization of animal sarbecoviruses and SARS-CoV-2 variants by screening single mouse B cells secreting IgGs that bind two or more sarbecovirus RBDs. Single-particle cryo-EM structures of antibody-spike complexes, including a Fab-Omicron complex, mapped neutralizing mAbs to conserved class 1/4 RBD epitopes. Structural analyses revealed neutralization mechanisms, potentials for intra-spike trimer cross-linking by IgGs, and induced changes in trimer upon Fab binding. In addition, we identified a mAb-resembling Bebtelovimab, an EUA-approved human class 3 anti-RBD mAb. These results support using mosaic RBD-nanoparticle vaccination to generate and identify therapeutic pan-sarbecovirus and pan-variant mAbs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The crystal structure of Vibrio cholerae (6-4) photolyase reveals interactions with cofactors and a DNA-binding region

Photolyases (PLs) reverse UV-induced DNA damage using blue light as an energy source. Of these PLs, (6-4) PLs repair (6-4)-lesioned photoproducts. We recently identified a gene from Vibrio cholerae (Vc) encoding a (6-4) PL, but structural characterization is needed to elucidate specific interactions with the chromophore cofactors. Here, we determined the crystal structure of Vc (6-4) PL at 2.5 Å resolution. Our high-resolution structure revealed that the two well-known cofactors, flavin adenine dinucleotide and the photoantenna 6,7-dimethyl 8-ribityl-lumazin (DMRL), stably interact with an α-helical and an α/β domain, respectively. Additionally, the structure has a third cofactor with distinct electron clouds corresponding to a [4Fe-4S] cluster. Moreover, we identified that Asp106 makes a hydrogen bond with water and DMRL, which indicates further stabilization of the photoantenna DMRL within Vc (6-4) PL. Further analysis of the Vc (6-4) PL structure revealed a possible region responsible for DNA binding. The region located between residues 478 to 484 may bind the lesioned DNA, with Arg483 potentially forming a salt bridge with DNA to stabilize further the interaction of Vc (6-4) PL with its substrate. Our comparative analysis revealed that the DNA lesion could not bind to the Vc (6-4) PL in a similar fashion to the Drosophila melanogaster (Dm, (6-4)) PL without a significant conformational change of the protein. The 23rd helix of the bacterial (6-4) PLs seems to have remarkable plasticity, and conformational changes facilitate DNA binding. In conclusion, our structure provides further insight into DNA repair by a (6-4) PL containing three cofactors.

(6-4) lesioned photoproduct↗
DOE OSTI2023

Structural basis of DNA binding by the WhiB-like transcription factor WhiB3 in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) WhiB3 is an iron–sulfur cluster-containing transcription factor belonging to a subclass of the WhiB-Like (Wbl) family that is widely distributed in the phylum Actinobacteria. WhiB3 plays a crucial role in the survival and pathogenesis of Mtb. It binds to the conserved region 4 of the principal sigma factor (σ A 4 ) in the RNA polymerase holoenzyme to regulate gene expression like other known Wbl proteins in Mtb. However, the structural basis of how WhiB3 coordinates with σ A 4 to bind DNA and regulate transcription is unclear. Here we determined crystal structures of the WhiB3:σ A 4 complex without and with DNA at 1.5 Å and 2.45 Å, respectively, to elucidate how WhiB3 interacts with DNA to regulate gene expression. These structures reveal that the WhiB3:σ A 4 complex shares a molecular interface similar to other structurally characterized Wbl proteins and also possesses a subclass-specific Arg-rich DNA-binding motif. We demonstrate that this newly defined Arg-rich motif is required for WhiB3 binding to DNA in vitro and transcriptional regulation in Mycobacterium smegmatis. Together, our study provides empirical evidence of how WhiB3 regulates gene expression in Mtb by partnering with σ A 4 and engaging with DNA via the subclass-specific structural motif, distinct from the modes of DNA interaction by WhiB1 and WhiB7.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

An alcove at the acetyl-CoA synthase nickel active site is required for productive substrate CO binding and anaerobic carbon fixation

One of the seven natural CO 2 fixation pathways, the anaerobic Wood-Ljungdahl pathway (WLP) is unique in generating CO as a metabolic intermediate, operating through organometallic intermediates, and in conserving (versus utilizing) net ATP. The key enzyme in the WLP is acetyl-CoA synthase (ACS), which uses an active site [2Ni-4Fe-4S] cluster (A-cluster), a CO tunnel, and an organometallic (Ni-CO, Ni-methyl, and Ni-acetyl) reaction sequence to generate acetyl-CoA. Here, we reveal that an alcove, which interfaces the tunnel and the A-cluster, is essential for CO 2 fixation and autotrophic growth by the WLP. In vitro spectroscopy, kinetics, binding, and in vivo growth experiments reveal that a Phe229A substitution at one wall of the alcove decreases CO affinity thirty-fold and abolishes autotrophic growth; however, a F229W substitution enhances CO binding 80-fold. Our results indicate that the structure of the alcove is exquisitely tuned to concentrate CO near the A-cluster; protect ACS from CO loss during catalysis, provide a haven for inhibitory CO, and stabilize the tetrahedral coordination at the Nip site where CO binds. The directing, concentrating, and protective effects of the alcove explain the inability of F209A to grow autotrophically. The alcove also could help explain current controversies over whether ACS binds CO and methyl through a random or ordered mechanism. Our work redefines what we historically refer to as the metallocenter "active site". The alcove is so crucial for enzymatic function that we propose it is part of the active site. The community should now look for such alcoves in all "gas handling" metalloenzymes.

bacterial metabolism↗
DOE OSTI2024

Functional diversification within the heme-binding split-barrel family

Due to neofunctionalization, a single fold can be identified in multiple proteins that have distinct molecular functions. Depending on the time that has passed since gene duplication and the number of mutations, the sequence similarity between functionally divergent proteins can be relatively high, eroding the value of sequence similarity as the sole tool for accurately annotating the function of uncharacterized homologs. Here, we combine bioinformatic approaches with targeted experimentation to reveal a large multifunctional family of putative enzymatic and nonenzymatic proteins involved in heme metabolism. This family (homolog of HugZ (HOZ)) is embedded in the “FMN-binding split barrel” superfamily and contains separate groups of proteins from prokaryotes, plants, and algae, which bind heme and either catalyze its degradation or function as nonenzymatic heme sensors. In prokaryotes these proteins are often involved in iron assimilation, whereas several plant and algal homologs are predicted to degrade heme in the plastid or regulate heme biosynthesis. In the plant Arabidopsis thaliana, which contains two HOZ subfamilies that can degrade heme in vitro (HOZ1 and HOZ2), disruption of AtHOZ1 (AT3G03890) or AtHOZ2A (AT1G51560) causes developmental delays, pointing to important biological roles in the plastid. In the tree Populus trichocarpa, a recent duplication event of a HOZ1 ancestor has resulted in localization of a paralog to the cytosol. Structural characterization of this cytosolic paralog and comparison to published homologous structures suggests conservation of heme-binding sites. This study unifies our understanding of the sequence-structure-function relationships within this multilineage family of heme-binding proteins and presents new molecular players in plant and bacterial heme metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

The HypA and HypB metallochaperones from Methanococcus maripaludis have unique metal-binding properties and a distinct nickel transfer mechanism

[NiFe] hydrogenases are widespread microbial metalloenzymes that catalyze the reversible conversion of hydrogen (H2) to protons and electrons, playing key roles in energy metabolism. The biosynthesis of the NiFe(CN) 2 CO cofactor involves a suite of maturation proteins, including the HypA and HypB nickel metallochaperones. Here, we define the metal-binding properties, nucleotide-dependent behavior, and functional interplay of HypA and HypB from the hydrogenotrophic methanogenic archaeon, Methanococcus maripaludis . Methanogens have multiple essential nickel-dependent enzymes, so they require efficient systems for nickel delivery that remain largely unexplored. Purified M. maripaludis HypA binds zinc or mononuclear iron at the C-terminal metal binding site, the latter of which has not been reported in other HypA proteins and may serve a unique regulatory role in methanogens. The G-protein metallochaperone HypB binds nickel at the G-domain, which stimulates GTPase activity. Size exclusion chromatography experiments reveal that HypA and HypB form complexes in the presence of nickel, and zinc-bound HypA is optimized for nickel transfer from HypB. The identity of the nucleotide bound to HypB (GDP or GTP) alters the oligomeric state of HypA-HypB complexes, supporting a GTPase-mediated nickel delivery pathway. The HypA-HypB 2 complex configuration is enriched and stable in the presence of GDP and nickel, indicating that this complex delivers nickel to the hydrogenase as opposed to HypA alone. Interestingly, affinity purification-mass spectrometry revealed that HypB interacts with several nickel-dependent proteins, suggesting that HypB may play a broader role in nickel homeostasis in M. maripaludis . Together, this work establishes a biochemical framework for HypAB-mediated nickel trafficking in methanogens.

[NiFe] hydrogenase↗
DOE OSTI2023

Multi-level Force-dependent Allosteric Enhancement of αE-catenin Binding to F-actin by Vinculin

Classical cadherins are transmembrane proteins whose extracellular domains link neighboring cells, and whose intracellular domains connect to the actin cytoskeleton via β-catenin and α-catenin. The cadherin-catenin complex transmits forces that drive tissue morphogenesis and wound healing. In addition, tension-dependent changes in αE-catenin conformation enables it to recruit the actin-binding protein vinculin to cell-cell junctions, which contributes to junctional strengthening. How and whether multiple cadherin-complexes cooperate to reinforce cell-cell junctions in response to load remains poorly understood. Here, we used single-molecule optical trap measurements to examine how multiple cadherin-catenin complexes interact with F-actin under load, and how this interaction is influenced by the presence of vinculin. We show that force oriented toward the (-) end of the actin filament results in mean lifetimes 3-fold longer than when force was applied towards the barbed (+) end. We also measured force-dependent actin binding by a quaternary complex comprising the cadherin-catenin complex and the vinculin head region, which cannot itself bind actin. Binding lifetimes of this quaternary complex increased as additional complexes bound F-actin, but only when load was oriented toward the (-) end. In contrast, the cadherin-catenin complex alone did not show this form of cooperativity. Furthermore, these findings reveal multi-level, force-dependent regulation that enhances the strength of the association of multiple cadherin/catenin complexes with F-actin, conferring positive feedback that may strengthen the junction and polarize F-actin to facilitate the emergence of higher-order cytoskeletal organization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Design and Function of α-Helix-Rich, Heme-Binding Peptide Materials

Peptide materials often employ short peptides that self-assemble into unique nanoscale architectures and have been employed across many fields relevant to medicine and energy. A majority of peptide materials are high in beta-sheet, secondary structure content, including heme-binding peptide materials. To broaden the structural diversity of heme-binding peptide materials, a small series of peptides were synthesized to explore the design criteria required for (1) folding into an alpha-helix structure, (2) assembling into a nanoscale material, (3) binding heme, and (4) demonstrating functions similar to that of heme proteins. One peptide was identified to meet all four criteria, including the heme protein function of CO binding and its microsecond-to-millisecond recombination rates, as measured by transient absorption spectroscopy. In conclusion, implications of new design criteria and peptide material function through heme incorporation are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Deeper Understanding of Mononuclear Manganese(IV)–Oxo Binding Brønsted and Lewis Acids and the Manganese(IV)–Hydroxide Complex

Binding of Lewis acidic metal ions and Brønsted acid at the metal–oxo group of high-valent metal–oxo complexes enhances their reactivities significantly in oxidation reactions. However, such a binding of Lewis acids and proton at the metal–oxo group has been questioned in several cases and remains to be clarified. In this work, we report the synthesis, characterization, and reactivity studies of a mononuclear manganese(IV)–oxo complex binding triflic acid, {[(dpaq)Mn IV (O)]–HOTf} + (1–HOTf). First, 1–HOTf was synthesized and characterized using various spectroscopic techniques, including resonance Raman (rRaman) and X-ray absorption spectroscopy/extended X-ray absorption fine structure. In particular, in rRaman experiments, we observed a linear correlation between the Mn–O stretching frequencies of 1–HOTf (e.g., ν Mn–O at ~793 cm –1 ) and 1–M n+ (Mn + = Ca 2+ , Zn 2+ , Lu 3+ , Al 3+ , or Sc 3+ ) and the Lewis acidities of H + and M n+ ions, suggesting that H+ and Mn+ bind at the metal–oxo moiety of [(dpaq)Mn IV (O)] + . Interestingly, a single-crystal structure of 1–HOTf was obtained by X-ray diffraction analysis, but the structure was not an expected Mn(IV)–oxo complex but a Mn(IV)–hydroxide complex, [(dpaq)Mn IV (OH)](OTf) 2 (4), with a Mn–O bond distance of 1.8043(19) Å and a Mn–O stretch at 660 cm –1 . More interestingly, 4 reverted to 1–HOTf upon dissolution, demonstrating that 1–HOTf and 4 are interconvertible depending on the physical states, such as 1–HOTf in solution and 4 in isolated solid. The reactivity of 1–HOTf was investigated in hydrogen atom transfer (HAT) and oxygen atom transfer (OAT) reactions and then compared with those of 1–Mn+ complexes; an interesting correlation between the Mn–O stretching frequencies of 1–HOTf and 1–M n+ and their reactivities in the OAT and HAT reactions is reported for the first time in this study.

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

CAML: Commutative Algebra Machine Learning─A Case Study on Protein–Ligand Binding Affinity Prediction

Recently, Suwayyid and Wei introduced commutative algebra as an emerging paradigm for machine learning and data science. In this work, we propose commutative algebra machine learning (CAML) for the prediction of protein−ligand binding affinities. Specifically, we apply persistent Stanley−Reisner theory, a key concept in combinatorial commutative algebra, to the affinity predictions of protein−ligand binding and metalloprotein−ligand binding. We present three new algorithms, i.e., element-specific commutative algebra, category-specific commutative algebra, and commutative algebra on bipartite complexes, to tackle the complexity of data involved in (metallo) protein−ligand complexes. We show that the proposed CAML outperforms other state-of-theart methods in (metallo) protein−ligand binding affinity predictions, indicating the great potential of commutative algebra learning.

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

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

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

Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions

Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems.

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

Elucidating the Influence of Metal Surface Composition on Organic Adsorbate Binding Using Active Particle Dynamics

The adsorption strengths of organic compounds on metal surfaces are sensitive to the metal composition, and they play a central role in many catalytic reactions, helping to control the coverage of the reactant and altering the overall reaction rate. While adsorption energies are straightforward to measure and calculate in vacuum and gas-phase environments, adsorption energetics can be dramatically altered by the presence of solvent in liquid-phase reactions. However, the effects of metal composition on binding strengths in a liquid environment are less well understood, primarily due to the difficulty of accurate in situ measurements of organic binding on metal surfaces in the liquid phase. Here, we utilize the motion of active particles in water to probe the adsorption energies of an organic adsorbate (furfural) on a range of metal surfaces (pure Pd, pure Pt, and four PdAu alloy compositions) to elucidate the effect of metal composition. Janus particles with catalytic caps of particular metal compositions all exhibited active motion resulting from consumption of H 2 O 2 ; adsorbate binding was inferred through the decrease in the velocity of active motion and was modeled by a Langmuir adsorption isotherm. The measured adsorption affinities were used to extract the adsorption enthalpy of furfural on the different metals. The Pd surface was found to bind furfural more strongly than the Pt surface by some 10 kJ/mol. Furthermore, the adsorption of furfural on the alloys was found to increase monotonically in magnitude with Pd content. As a result, the data reported herein aid the development of accurate understanding of organic adsorption in the presence of solvent and the role of the metal surface in tuning adsorption strengths to optimize catalytic processes in the liquid phase.

36 MATERIALS SCIENCE↗