Search for: All records

Abstract Human carbonic anhydrase (CA) metalloenzymes utilize a Zn 2+ -containing active site to catalyze the interconversion of carbon dioxide to bicarbonate. The Zn 2+ ion may be replaced with other divalent transition metals, though the catalytic efficiency of the enzyme will be reduced. In this work, quantum mechanical cluster models of the active site are used to map the reaction profile for the hydration mechanism of carbon dioxide. The Lipscomb proton transfer and Lindskog rotation mechanisms were examined for the native Zn 2+ -enzyme along with variants where the metal was substituted with Cd 2+ , Ni 2+ , Fe 2+ , and Fe 3+ . The findings highlight the impact the metal coordination geometry has on the reaction profile. The results also suggest Fe 2+ , which is the functional metal for a prototypical CA of an anaerobic bacterium, might also be functional for human CA if cultured within an anaerobic environment. 

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

 

Diethylenetriaminepentaacetic acid (DTPA) is a chelating agent whose complex with the Gd³⁺ion is used in medical imaging. DTPA is also used in lanthanide‐actinide separation processes. As protonation of the DTPA ligand can facilitate dissociation of the Gd³⁺ion from the Gd‐DTPA complex, this work investigates the coordination structures of the aqueous Gd³⁺ion and its environment when chelated by DTPA in eight different DTPA protonation states. Both classical and ab initio molecular dynamics (MD) simulations are conducted to model the solvated complexes. Extended X‐ray absorption fine structure (EXAFS) measurements of the Gd³⁺aqua ion, and the Gd‐DTPA complex at pH 1 and 11, are compared to EXAFS spectra predicted from the MD simulations to verify the accuracy of the MD structures. The findings of this work provide atomic‐level details into the fluctuating Gd‐DTPA complex environment as the DTPA ligand gradually detaches from the Gd³⁺ion with increased protonation.

 

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ∼7 and ∼11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ∼11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ∼7 and ∼11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. Increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature. 

The restricted rotation of chemical bonds may lead to the formation of stable, conformationally chiral molecules. While the asymmetry in chiral molecules is generally observed in the presence of one or more stereocenters, asymmetry exhibited by conformational chirality in compounds lacking stereocenters, called atropisomerism, depends on structural and temperature factors that are still not fully understood. This atropisomerism is observed in natural diarylether heptanoids where the length of the intramolecular tether constrains the compounds to isolable enantiomers at room temperature. In this work, we examine the impact tether length has on the activation free energies to isomerization of a diarylether cyclophane substructure with a tether ranging from 6 to 14 carbons. Racemization activation energies are observed to decay from 48 kcal/mol for a 7-carbon tether to 9.2 kcal/mol for a 14-carbon tether. Synthetic efforts to experimentally test these constraints are also presented. This work will likely guide the design and synthesis of novel asymmetric cyclophanes that will be of interest in the catalysis community given the importance of atropisomeric ligands in the field of asymmetric catalysis.