Search for: All records

The transfer of lanthanide–ligand complexes across aqueous–organic interfaces was studied with rare event molecular dynamics simulations. Relative solubilities were quantified from potentials of mean force. 

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison toin vivostructure and mechanism. This work introduces an MD to QM-cluster model workflow. 

Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu3+ ion complexed with EDTA in all its protonation states in aqueous solution. Eu-EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X-ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu3+ aqua ion and Eu-EDTA complexes at pH 3 and 11. This work details how Eu-EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu3+ ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation. 

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. 

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.