2025 Metals in Biology Gordon Research Conference and Bioinorganic Chemistry Gordon Research Seminar
Final Conference Program & Meeting Summary
SEARCH · Engineering Papers
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.
Final Conference Program & Meeting Summary
The overarching goal of this project was to elucidate biochemical mechanisms that ubiquitous microorganisms use to enable primary metabolism based on nitrification, the multi-electron oxidation of ammonia (NH 3 ) to nitrite (NO 2 – ). This chemistry enables life to substitute carbon based fuel with nitrogen. Understanding the underlying chemistry was pursued to help inform the use of ammonia as an alternative fuel source compared to petroleum, capitalizing on the widespread infrastructure for production and distribution of this substance in the US. This project specifically focused on characterizing ammonia monooxygenase, an integral membrane copper protein that initiates nitrification via the aerobic oxidation of NH 3 to hydroxylamine (NH 2 OH).
A machine-learning (ML) model that predicts metal–ligand binding constants was developed using the open-source Chemprop software. The model was trained on over 30,000 experimental log K 1 values, which include both protonation and metal–ligand stability constants, comprising over 3500 ligands and 10 2 metal ions from 73 total elements, thus generalizing beyond existing limited approaches, which focus only on specific metals or ligand families. The best-performing model included a combination of SMILES-based molecular representations along with descriptors for the metal ion and experimental conditions. It had an external test R 2 value of 0.942, and MAE value of 0.834. A “SMILES-only” simpler version also produced accurate predictions and preserved the binding trends, serving as a quick and easily accessible alternative for users without computational expertise. The SMILES-only model performed comparably to density functional theory (DFT) calculations but utilized a fraction of the computational resources. The model was successfully applied across diverse domains, including bioinorganic chemistry, heavy metal remediation, and sensor development and demonstrated its effectiveness as a rapid and reliable screening tool for both academic and industrial uses.
Metal-binding proteins have the exceptional ability to facilitate long-range electron transport in nature. Despite recent progress, the sequence-structure–function relationships governing electron transport in heme-binding peptides and protein assemblies are not yet fully understood. In this work, the electronic properties of a series of heme-binding peptides inspired by cytochrome bc1 are studied using a combination of molecular electronics experiments, molecular modeling, and simulation. Self-assembled monolayers (SAMs) are prepared using sequence-defined heme-binding peptides capable of forming helical secondary structures. Following monolayer formation, the structural properties and chemical composition of assembled peptides are determined using atomic force microscopy and X-ray photoelectron spectroscopy, and the electronic properties (current density–voltage response) are characterized using a soft contact liquid metal electrode method based on eutectic gallium–indium alloys (EGaIn). Our results show a substantial 1000-fold increase in current density across SAM junctions upon addition of heme compared to identical peptide sequences in the absence of heme, while maintaining a constant junction thickness. These findings show that amino acid composition and sequence directly control enhancements in electron transport in heme-binding peptides. Overall, this study demonstrates the potential of using sequence-defined synthetic peptides inspired by nature as functional bioelectronic materials.
Summary Our work reveals the structure of the active state of Methyl-Coenzyme M Reductase (MCR), the key and rate-limiting enzyme in biological methane formation. We find large differences between the active Ni(I) and inactive Ni(II) proteins and provide insight into how nature makes and breaks the C-H bond of methane. The Ni(II)-F430 center in inactive MCR contains four planar nitrogen ligands, a lower axial glutamine oxo, and an upper axial thiolate. The Ni(I)-enzyme replaces the axial ligands with a single water. The one-electron redox change results in movement of the Ni ion and upward swing of the β-lactam ring in the tetrapyrrole coupled to a domino-like protein quake through second sphere residues, inter-subunit interactions, a substrate tunnel, affecting even the dimensions of the unit cell. These structural changes lead Ni(I)-MCR to release a charge clamp that, in the Ni(II) state, locks down substrate Coenzyme B. Determining the Ni(I)-MCR structure required development of rigorous anaerobic crystallographic techniques. Validation of the MCR redox state was accomplished by in-line and parallel spectroscopic and unit cell analyses. This structure has large implications for developing technologies to limit methane emissions and efficiently produce biofuels. Methodology described here will enhance structural biology for other oxygen-sensitive enzymes.
The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.
High-valent Fe IV =O species are common intermediates in biological and artificial catalysts. Heme and nonheme S=1 Fe IV =O sites have been synthesized and studied for decades but little quantitative experimental comparison of their electronic structures has been available, due to the lack of direct methods focused on the iron. This study allows a rigorous determination of the electronic structure of a nonheme Fe IV =O center and its comparison to an Fe IV =O heme site using 1s2p resonant inelastic X-ray scattering (RIXS) and Fe L-edge X-ray absorption spectroscopy (XAS). Further, variable temperature magnetic circular dichroism (VT-MCD) of the ligand field transitions, combined with nuclear resonance vibrational spectroscopy of the two S=1 Fe IV =O systems show that the equatorial ligand field decreases from a nonheme to a heme Fe IV =O site. Alternatively, RIXS and Fe L-edge XAS combined with MCD show that the Fe dπ orbitals are unperturbed in the Fe IV =O heme relative to the nonheme site because the strong axial Fe-O bond uncouples the Fe dπ orbitals from the porphyrin π-system. As a consequence, the thermodynamics and kinetics of the H-atom abstraction reactions are actually very similar for heme compound II and nonheme Fe IV =O active sites.