Search

The interaction of chloroquine (CQ) and the μ-oxo dimer of iron(III) protoporphyrin IX (ferriheme) in aqueous solution was modeled using molecular dynamics (MD) simulations. Two models of the CQ-(μ-oxo ferriheme) complex were investigated, one involving CQ π-stacked with an unligated porphyrin face of μ-oxo ferriheme and the other in which CQ was docked between the two porphyrin rings. The feasibility of both models was tested by fitting computed structures to the experimental extended X-ray absorption fine structure (EXAFS) spectrum of the CQ-(μ-oxo ferriheme) complex in frozen aqueous solution. The docked model produced better agreement with experimental data, suggesting that this is the more likely structure in aqueous solution. The EXAFS fit indicated a longer than expected Fe-O bond of 1.87Å, accounting for the higher than expected magnetic moment of the complex. As a consequence, the asymmetric Fe-O-Fe stretch shifts much lower in frequency and was identified in the precipitated solid at 744cm(-1) with the aid of the O(18) isomer shift. Three important CQ-ferriheme interactions were identified in the docked structure. These were a hydrogen bond between the oxide bridge of μ-oxo ferriheme and the protonated quinolinium nitrogen atom of CQ; π-stacking between the quinoline ring of CQ and the porphyrin rings; and a close contact between the 7-chloro substituent of CQ and the porphyrin methyl hydrogen atoms. These interactions can be used to rationalize previously observed structure-activity relationships for quinoline-ferriheme association., (Copyright © 2015 Elsevier Inc. All rights reserved.)

CHARMM force field parameters have been developed to model nonprotein bound five-coordinate ferriheme (ferriprotoporphyrin IX) species in aqueous solution. Structures and solvation were determined from molecular dynamics (MD) simulations at 298 K of monomeric [HO-ferriheme](2-), [H2O-ferriheme](-), and [H2O-ferriheme](0); π-π dimeric [(HO-ferriheme)2](4-), [(H2O-ferriheme)(HO-ferriheme)](3-), [(H2O-ferriheme)2](2-), and [(H2O-ferriheme)2](0); and μ-oxo dimeric [μ-(ferriheme)2O](4-). Solvation of monomeric species predominated around the axial ligand, meso-hydrogen atoms of the porphyrin ring (Hmeso), and the unligated face. Existence of π-π ferriheme dimers in aqueous solution was supported by MD calculations where such dimers remained associated over the course of the simulation. Porphyrin rings were essentially coplanar. In these dimers major and minor solvation was observed around the axial ligand and Hmeso positions, respectively. In μ-oxo ferriheme, strong solvation of the unligated face and bridging oxide ligand was observed. The solution structure of the μ-oxo dimer was investigated using extended X-ray absorption fine structure (EXAFS) spectroscopy. The EXAFS spectrum obtained from frozen solution was markedly different from that recorded on dried μ-oxo ferriheme solid. Inclusion of five solvent molecules obtained from spatial distribution functions in the structure generated from MD simulation was required to produce acceptable fits to the EXAFS spectra of the dimer in solution, while the solid was suitably fitted using the crystal structure of μ-oxo ferriheme dimethyl ester which included no solvent molecules.

Speciation of ferriprotoporphyrin IX, Fe(III)PPIX, in aqueous solution is complex. Despite the use of its characteristic spectroscopic features for identification, the theoretical basis of the unique UV-visible absorbance spectrum of μ-[Fe(III)PPIX](2)O has not been explored. To investigate this and to establish a structural and spectroscopic model for Fe(III)PPIX species, density functional theory (DFT) calculations were undertaken for H(2)O-Fe(III)PPIX and μ-[Fe(III)PPIX](2)O. The models agreed with related Fe(III)porphyrin crystal structures and reproduced vibrational spectra well. The UV-visible absorbance spectra of H(2)O-Fe(III)PPIX and μ-[Fe(III)PPIX](2)O were calculated using time-dependent DFT and reproduced major features of the experimental spectra of both. Transitions contributing to calculated excitations have been identified. The features of the electronic spectrum calculated for μ-[Fe(III)PPIX](2)O were attributed to delocalization of electron density between the two porphyrin rings of the dimer, the weaker ligand field of the axial ligand, and antiferromagnetic coupling of the Fe(III) centers. Room temperature magnetic circular dichroism (MCD) spectra have been recorded and are shown to be useful in distinguishing between these two Fe(III)PPIX species. Bands underlying major spectroscopic features were identified through simultaneous deconvolution of UV-visible and MCD spectra. Computed UV-visible spectra were compared to deconvoluted spectra. Interpretation of the prominent bands of H(2)O-Fe(III)PPIX largely conforms to previous literature. Owing to the weak paramagnetism of μ-[Fe(III)PPIX](2)O at room temperature and the larger number of underlying excitations, interpretation of its experimental UV-visible spectrum was necessarily tentative. Nonetheless, comparison with the calculated spectra of antiferromagnetically coupled and paramagnetic forms of the μ-oxo dimer of Fe(III)porphine suggested that the composition of the Soret band involves a mixture of π→π* and π→d(π) charge transfer transitions. The Q-band and charge transfer bands appear to amalgamate into a mixed low energy envelope consisting of excitations with heavily admixed π→π* and charge transfer transitions.