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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Localization and functional exploration of leiomodin-2’s C-terminal binding sites

Striated muscle contraction occurs through interactions between overlapping myosin-based thick and actin-based thin filaments within the sarcomere. For effective contraction to occur, the length of the thin filament must be maintained to allow for sufficient overlap with the thick filament. The proteins leiomodin and tropomodulin compete for binding at the pointed end of thin filaments to regulate their length, utilizing their homologous N-terminal actin and tropomyosin binding sites. Leiomodin also has a region called the C-terminal extension, absent in tropomodulin. In this region, the cardiac isoform (leiomodin-2) contains additional actin-binding sites that enable it to bind along the sides of thin filaments in a Ca2+-dependent manner. Here, using nuclear magnetic resonance spectroscopy, we localize the regions of the C-terminal extension that contain residues involved in thin filament side-binding. Using co-sedimentation assays, we reveal that these regions can independently bind thin filaments and discover that the poly-proline region plays a role as a linker, maintaining an adequate distance between two of the regions required for effective interaction to occur. In addition to its role in side-binding, we provide direct evidence that the poly-proline region interacts with profilin and propose a new mechanism by which leiomodin-2 may assist in the polymerization of profilin-bound actin at thin filament pointed ends.

ACTIN-BINDING PROTEINS↗
DOE OSTI2023

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro WT ) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ~400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (T m ), respectively, relative to MPro WT . The H41A mutation results in a 20-fold increase in the dissociation constant (K d ), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T m , respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro H41A , whereas no significant change is observed in binding to MPro WT . Structures of the four inhibitor complexes with MPro 1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro C145A also results in a significant, but smaller, increase in K d and decrease in ΔG and T m , relative to NMV.

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

Specific iron binding to natural sphingomyelin membrane induced by non-specific co-solutes

Sphingomyelin (SPM), a crucial phospholipid in the myelin sheath, plays a vital role in insulating nerve fibers. We hypothesize that iron ions selectively bind to the phosphatidylcholine (PC) template within the SPM membrane under near-physiological conditions, resulting in disruptions to membrane organization. These interactions could potentially contribute to the degradation of the myelin sheath, thereby playing a role in the development of neurodegenerative diseases. We utilized synchrotron-based X-ray spectroscopy and diffraction techniques to study the interaction of iron ions with a bovine spinal-cord SPM monolayer (ML) at the liquid-vapor interface under physiological conditions. The SPM ML serves as a model system, representing localized patches of lipids within a more complex membrane structure. The experiments assessed iron binding to the SPM membrane both in the presence of salts and with additional evaluation of the effects of various ion species on membrane behavior. Grazing incidence X-ray diffraction was employed to analyze the impact of iron binding on the structural integrity of the SPM membrane. Furthermore, our results demonstrate that iron ions in dilute solution selectively bind to the PC template of the SPM membrane exclusively at near-physiological salt concentrations (e.g., NaCl, KCl, KI, or CaCl 2 ) and are pH-dependent. In-significant binding was detected in the absence of these salts or at near-neutral pH with salts. The surface adsorption of iron ions is correlated with salt concentration, reaching saturation at physiological levels. In contrast, multivalent ions such as La 3+ and Ca 2+ do not bind to SPM under similar conditions. Notably, iron binding to the SPM membrane disrupts its in-plane organization, suggesting that these interactions may compromise membrane integrity and contribute to myelin sheath damage associated with neurological disorders.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single–amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single– and multiple–amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

DeepGRN: prediction of transcription factor binding site across cell-types using attention-based deep neural networks

Abstract Background Due to the complexity of the biological systems, the prediction of the potential DNA binding sites for transcription factors remains a difficult problem in computational biology. Genomic DNA sequences and experimental results from parallel sequencing provide available information about the affinity and accessibility of genome and are commonly used features in binding sites prediction. The attention mechanism in deep learning has shown its capability to learn long-range dependencies from sequential data, such as sentences and voices. Until now, no study has applied this approach in binding site inference from massively parallel sequencing data. The successful applications of attention mechanism in similar input contexts motivate us to build and test new methods that can accurately determine the binding sites of transcription factors. Results In this study, we propose a novel tool (named DeepGRN) for transcription factors binding site prediction based on the combination of two components: single attention module and pairwise attention module. The performance of our methods is evaluated on the ENCODE-DREAM in vivo Transcription Factor Binding Site Prediction Challenge datasets. The results show that DeepGRN achieves higher unified scores in 6 of 13 targets than any of the top four methods in the DREAM challenge. We also demonstrate that the attention weights learned by the model are correlated with potential informative inputs, such as DNase-Seq coverage and motifs, which provide possible explanations for the predictive improvements in DeepGRN. Conclusions DeepGRN can automatically and effectively predict transcription factor binding sites from DNA sequences and DNase-Seq coverage. Furthermore, the visualization techniques we developed for the attention modules help to interpret how critical patterns from different types of input features are recognized by our model.

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

Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Isolation and characterization of a novel calmodulin-binding protein from potato

Tuberization in potato is controlled by hormonal and environmental signals. Ca(2+), an important intracellular messenger, and calmodulin (CaM), one of the primary Ca(2+) sensors, have been implicated in controlling diverse cellular processes in plants including tuberization. The regulation of cellular processes by CaM involves its interaction with other proteins. To understand the role of Ca(2+)/CaM in tuberization, we have screened an expression library prepared from developing tubers with biotinylated CaM. This screening resulted in isolation of a cDNA encoding a novel CaM-binding protein (potato calmodulin-binding protein (PCBP)). Ca(2+)-dependent binding of the cDNA-encoded protein to CaM is confirmed by (35)S-labeled CaM. The full-length cDNA is 5 kb long and encodes a protein of 1309 amino acids. The deduced amino acid sequence showed significant similarity with a hypothetical protein from another plant, Arabidopsis. However, no homologs of PCBP are found in nonplant systems, suggesting that it is likely to be specific to plants. Using truncated versions of the protein and a synthetic peptide in CaM binding assays we mapped the CaM-binding region to a 20-amino acid stretch (residues 1216-1237). The bacterially expressed protein containing the CaM-binding domain interacted with three CaM isoforms (CaM2, CaM4, and CaM6). PCBP is encoded by a single gene and is expressed differentially in the tissues tested. The expression of CaM, PCBP, and another CaM-binding protein is similar in different tissues and organs. The predicted protein contained seven putative nuclear localization signals and several strong PEST motifs. Fusion of the N-terminal region of the protein containing six of the seven nuclear localization signals to the reporter gene beta-glucuronidase targeted the reporter gene to the nucleus, suggesting a nuclear role for PCBP.

NASA Discipline Plant Biology↗
DOE OSTI2023

Characterization of promoter elements of isoprene‐responsive genes and the ability of isoprene to bind START domain transcription factors

Abstract Isoprene has recently been proposed to be a signaling molecule that can enhance tolerance of both biotic and abiotic stress. Not all plants make isoprene, but all plants tested to date respond to isoprene. We hypothesized that isoprene interacts with existing signaling pathways rather than requiring novel mechanisms for its effect on plants. We analyzed the cis‐regulatory elements (CREs) in promoters of isoprene‐responsive genes and the corresponding transcription factors binding these promoter elements to obtain clues about the transcription factors and other proteins involved in isoprene signaling. Promoter regions of isoprene‐responsive genes were characterized using the Arabidopsis cis‐regulatory element database. CREs bind ARR1, Dof, DPBF, bHLH112, GATA factors, GT‐1, MYB, and WRKY transcription factors, and light‐responsive elements were overrepresented in promoters of isoprene‐responsive genes; CBF‐, HSF‐, WUS‐binding motifs were underrepresented. Transcription factors corresponding to CREs overrepresented in promoters of isoprene‐responsive genes were mainly those important for stress responses: drought‐, salt/osmotic‐, oxidative‐, herbivory/wounding and pathogen‐stress. More than half of the isoprene‐responsive genes contained at least one binding site for TFs of the class IV (homeodomain leucine zipper) HD‐ZIP family, such as GL2, ATML1, PDF2, HDG11, ATHB17. While the HD‐zipper‐loop‐zipper (ZLZ) domain binds to the L1 box of the promoter region, a special domain called the steroidogenic acute regulatory protein‐related lipid transfer, or START domain, can bind ligands such as fatty acids (e.g., linolenic and linoleic acid). We tested whether isoprene might bind in such a START domain. Molecular simulations and modeling to test interactions between isoprene and a class IV HD‐ZIP family START‐domain‐containing protein were carried out. Without membrane penetration by the HDG11 START domain, isoprene within the lipid bilayer was inaccessible to this domain, preventing protein interactions with membrane bound isoprene. The cross‐talk between isoprene‐mediated signaling and other growth regulator and stress signaling pathways, in terms of common CREs and transcription factors could enhance the stability of the isoprene emission trait when it evolves in a plant but so far it has not been possible to say what how isoprene is sensed to initiate signaling responses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2 ) n , n = 1,2

Pyrimidine has two in-plane CH(δ+)/N̈(δ–)/CH(δ+) binding sites that are complementary to the (δ–/2δ+/δ–) quadrupole moment of CO 2 . For this study, we recorded broadband microwave spectra over the 7.5–17.5 GHz range for pyrimidine-(CO 2 ) n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14 N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO 2 complex is planar, with CO 2 occupying one of the quadrupolar binding sites, forming a structure in which the CO 2 is stabilized in the plane by interactions with the C–H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO 2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191–195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C 2v symmetry in which the second CO 2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO 2 is –13.7 kJ mol –1 , including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO 2 ) 2 is almost exactly double that size, indicating little interaction between the two CO 2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

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

Revealing the Influence of Binding Motifs on Electron Transfer and Recombination Kinetics for CdSe Quantum Dots Functionalized with a Modified Viologen

Anchoring of molecules to the surfaces of semiconductor nanocrystals (NCs) presents an opportunity to leverage the precise synthetic tunability of molecular function and the remarkable light harvesting properties of NCs to drive photochemical reactions. However, charge transfer between the two species depends not only on the energy level alignments but also on the details of their binding interactions, which are difficult to probe. Here, in this work, we characterize the binding between CdSe quantum dots (QDs) and a new phosphonated derivative of the electron acceptor methyl viologen, designed to attach to the QD surface via the phosphonate group. We use isothermal titration calorimetry to probe the thermodynamics of the QD–molecule interaction and use the parameters determined therein to analyze transient absorption spectroscopy measurements of forward and back electron transfer from QDs to the viologen. We find that the ligand-like phosphonate binding leads to an electron-transfer rate constant that is 3 orders of magnitude smaller than that for the face-on binding of the bipyridine ring of methyl viologen. Back electron transfer is also significantly slower in the derivative. Interestingly, a minor fraction of the phosphonated derivative also binds in the face-on configuration, with similar forward and back electron transfer kinetics as methyl viologen. Numerical simulations show that the ligand-like binding will lead to significantly improved quantum yields of photocatalysis over a wide range of reaction rates. By independently characterizing binding thermodynamics and charge transfer kinetics, this work reveals how the complexities underlying electron transfer at the NC–molecule interface determine photocatalytic outcomes. This work also represents a step toward controlling forward and back electron transfer kinetics via rational molecular design.

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

Molecular Dynamics Simulation of Atomic Interactions in the Vancomycin Binding Site

Vancomycin is a glycopeptide antibiotic produced by Amycolaptopsis orientalis used to treat serious infections by Grampositive pathogens including methicillin-resistant Staphylococcus aureus. Vancomycin inhibits cell wall biosynthesis by targeting lipid II, which is the membrane-bound peptidoglycan precursor. The heptapeptide aglycon structure of vancomycin binds to the D-Ala-D-Ala of the pentapeptide stem structure in lipid II. The third residue of vancomycin aglycon is asparagine, which is not directly involved in the dipeptide binding. Nonetheless, asparagine plays a crucial role in substrate recognition, as the vancomycin analogue with asparagine substituted by aspartic acid (V D ) shows a reduction in antibacterial activities. To characterize the function of asparagine, binding of vancomycin and its aspartic-acid-substituted analogue V D to L-Lys-D-Ala-D-Ala and L-Lys-D-Ala-D-Lac was investigated using molecular dynamic simulations. Binding interactions were analyzed using root-mean-square deviation (RMSD), two-dimensional (2D) contour plots, hydrogen bond analysis, and free energy calculations of the complexes. The analysis shows that the aspartate substitution introduced a negative charge to the binding cleft of V D , which altered the aglycon conformation that minimized the repulsive lone pair interaction in the binding of a depsipeptide. Our findings provide new insight for the development of novel glycopeptide antibiotics against the emerging vancomycin-resistant pathogens by chemical modification at the third residue in vancomycin to improve its binding affinity to the D-Ala-D-Lac-terminated peptidoglycan in lipid II found in vancomycin-resistant enterococci and vancomycinresistant S. aureus.

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

Structure of an open K ATP channel reveals tandem PIP 2 binding sites mediating the Kir6.2 and SUR1 regulatory interface

ATP-sensitive potassium (K ATP ) channels, composed of four pore-lining Kir6.2 subunits and four regulatory sulfonylurea receptor 1 (SUR1) subunits, control insulin secretion in pancreatic β-cells. K ATP channel opening is stimulated by PIP 2 and inhibited by ATP. Mutations that increase channel opening by PIP 2 reduce ATP inhibition and cause neonatal diabetes. Although considerable evidence has implicated a role for PIP 2 in K ATP channel function, previously solved open-channel structures have lacked bound PIP 2 , and mechanisms by which PIP 2 regulates K ATP channels remain unresolved. Here, we report the cryoEM structure of a K ATP channel harboring the neonatal diabetes mutation Kir6.2-Q52R, in the open conformation, bound to amphipathic molecules consistent with natural C18:0/C20:4 long-chain PI(4,5)P 2 at two adjacent binding sites between SUR1 and Kir6.2. The canonical PIP 2 binding site is conserved among PIP 2 -gated Kir channels. The non-canonical PIP 2 binding site forms at the interface of Kir6.2 and SUR1. Functional studies demonstrate both binding sites determine channel activity. Kir6.2 pore opening is associated with a twist of the Kir6.2 cytoplasmic domain and a rotation of the N-terminal transmembrane domain of SUR1, which widens the inhibitory ATP binding pocket to disfavor ATP binding. The open conformation is particularly stabilized by the Kir6.2-Q52R residue through cation-π bonding with SUR1-W51. Together, these results uncover the cooperation between SUR1 and Kir6.2 in PIP 2 binding and gating, explain the antagonistic regulation of K ATP channels by PIP 2 and ATP, and provide a putative mechanism by which Kir6.2-Q52R stabilizes an open channel to cause neonatal diabetes.

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

Structure of HIV-1 RRE stem-loop II identifies two conformational states of the high-affinity Rev binding site

During HIV infection, specific RNA-protein interaction between the Rev response element (RRE) and viral Rev protein is required for nuclear export of intron-containing viral mRNA transcripts. Rev initially binds the high-affinity site in stem-loop II, which promotes oligomerization of additional Rev proteins on RRE. Here, we present the crystal structure of RRE stem-loop II in distinct closed and open conformations. The high-affinity Rev-binding site is located within the three-way junction rather than the predicted stem IIB. The closed and open conformers differ in their non-canonical interactions within the three-way junction, and only the open conformation has the widened major groove conducive to initial Rev interaction. Rev binding assays show that RRE stem-loop II has high- and low-affinity binding sites, each of which binds a Rev dimer. We propose a binding model, wherein Rev-binding sites on RRE are sequentially created through structural rearrangements induced by Rev-RRE interactions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

The identity of metal ions surrounding DNA is key to its biological function and materials applications. In this work, we compare atomistic molecular dynamics simulations of double strand DNA (dsDNA) with four alkaline earth metal ions (Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ ) to elucidate the physical interactions that govern DNA–ion binding. Simulations accurately model the ion–phosphate distance of Mg 2+ and reproduce ion counting experiments for Ca 2+ , Sr 2+ , and Ba 2+ . Our analysis shows that alkaline earth metal ions prefer to bind at the phosphate backbone compared to the major groove and negligible binding occurs in the minor groove. Larger alkaline earth metal ions with variable first solvation shells (Ca 2+ , Sr 2+ , and Ba 2+ ) show both direct and indirect binding, where indirect binding increases with ion size. Mg 2+ does not fit this trend because the strength of its first solvation shell predicts indirect binding only. Ions bound to the phosphate backbone form fewer contacts per ion compared to the major groove. Within the major groove, metal ions preferentially bind to guanine–cystosine base pairs and form simultaneous contacts with the N7 and O6 atoms of guanine. Overall, we find that the interplay among ion size, DNA–ion interaction, and the size and flexibility of the first solvation shell are key to predicting how alkaline earth metal ions interact with DNA.

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

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

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

Allosteric autoregulation of DNA binding via a DNA-mimicking protein domain: a biophysical study of ZNF410–DNA interaction using small angle X-ray scattering

ZNF410 is a highly-conserved transcription factor, remarkable in that it recognizes a 15-base pair DNA element but has just a single responsive target gene in mammalian erythroid cells. ZNF410 includes a tandem array of five zinc-fingers (ZFs), surrounded by uncharacterized N- and C-terminal regions. Unexpectedly, full-length ZNF410 has reduced DNA binding affinity, compared to that of the isolated DNA binding ZF array, both in vitro and in cells. AlphaFold predicts a partially-folded N-terminal subdomain that includes a 30-residue long helix, preceded by a hairpin loop rich in acidic (aspartate/glutamate) and serine/threonine residues. This hairpin loop is predicted by AlphaFold to lie against the DNA binding interface of the ZF array. In solution, ZNF410 is a monomer and binds to DNA with 1:1 stoichiometry. Surprisingly, the single best-fit model for the experimental small angle X-ray scattering profile, in the absence of DNA, is the original AlphaFold model with the N-terminal long-helix and the hairpin loop occupying the ZF DNA binding surface. For DNA binding, the hairpin loop presumably must be displaced. After combining biophysical, biochemical, bioinformatic and artificial intelligence-based AlphaFold analyses, we suggest that the hairpin loop mimics the structure and electrostatics of DNA, and provides an additional mechanism, supplementary to sequence specificity, of regulating ZNF410 DNA binding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

N-Glycolylneuraminic Acid Binding of Avian and Equine H7 Influenza A Viruses

Influenza A viruses (IAV) initiate infection by binding to glycans with terminal sialic acids on the cell surface. Hosts of IAV variably express two major forms of sialic acid, N-acetylneuraminic acid (NeuAc) and N-glycolylneuraminic acid (NeuGc). NeuGc is produced in most mammals, including horses and pigs, but is absent in humans, ferrets, and birds. The only known naturally occurring IAV that exclusively bind NeuGc are extinct highly pathogenic equine H7N7 viruses. We determined the crystal structure of a representative equine H7 hemagglutinin (HA) in complex with NeuGc and observed high similarity in the receptor-binding domain with an avian H7 HA. To determine the molecular basis for NeuAc and NeuGc specificity, we performed systematic mutational analyses, based on the structural insights, on two distant avian H7 HAs and an H15 HA. We found that the A135E mutation is key for binding α2,3-linked NeuGc but does not abolish NeuAc binding. The additional mutations S128T, I130V, T189A, and K193R converted the specificity from NeuAc to NeuGc. We investigated the residues at positions 128, 130, 135, 189, and 193 in a phylogenetic analysis of avian and equine H7 HAs. This analysis revealed a clear distinction between equine and avian residues. The highest variability was observed at key position 135, of which only the equine glutamic acid led to NeuGc binding. These results demonstrate that genetically distinct H7 and H15 HAs can be switched from NeuAc to NeuGc binding and vice versa after the introduction of several mutations, providing insights into the adaptation of H7 viruses to NeuGc receptors.

59 BASIC BIOLOGICAL SCIENCES↗