The characterization of trace elements in minerals using extended X-ray absorption fine structure (EXAFS) spectroscopy constitutes a first step toward understanding how impurities and contaminants interact with the host phase and the environment. However, limitations to EXAFS interpretation complicate the analysis of trace concentrations of impurities that are distributed across multiple phases in a heterogeneous system. Ab initio molecular dynamics (AIMD)-informed EXAFS analysis was employed to investigate the immobilization of trace uranium associated with nanophase iron (oxyhydr)oxides, a model system for the geochemical sequestration of radiotoxic actinides. The reductive transformation of ferrihydrite [Fe(OH)] to nanoparticulate iron oxyhydroxide minerals in the presence of uranyl (UO) resulted in the preferential incorporation of U into goethite (α-FeOOH) over lepidocrocite (γ-FeOOH), even though reaction conditions favored the formation of excess lepidocrocite. This unexpected result is supported by atomically resolved transmission electron microscopy. We demonstrate how AIMD-informed EXAFS analysis lifts the strict statistical limitations and uncertainty of traditional shell-by-shell EXAFS fitting, enabling the detailed characterization of the local bonding environment, charge compensation mechanisms, and oxidation states of polyvalent impurities in complex multiphase mineral systems.