Computational modeling of the dizinc-ferroxidase complex of human H ferritin: direct comparison of the density functional theory calculated and experimental structures.

Density functional theory optimizations of structures of dizinc(II) complexes with a six-residue model of the ferroxidase center of human H ferritin have been performed and the results compared with the crystallographically determined structure of the complex as presented in Protein Data Bank file 2CEI. The model employs the full structures of Glu27, Glu62, His65, Glu107, Gln141, and Ala144, and the structural effect of Tyr34 is also examined. The mean absolute deviation from experiment of atomic positions in the best calculated structures is less than 0.3 A. The experimental structure is reproduced well enough to determine the coordination environment of the metal ions. Each zinc(II) center is pentacoordinate with a single water ligand, and the two centers are bridged by a hydroxide ion. Ala144 interacts weakly and repulsively with the rest of the complex. Tyr34 displays a weak attraction through a hydrogen bond to Glu107 that affects the orientation of that group.

DFT comparison of Fe2+ hydration in the binding sites of the ferroxidase center of bullfrog M ferritin.

Density functional theory optimizations have been conducted on structures of complexes of Fe(2+) with (H(2)O)(n) (n = 0-3) in three-residue models of binding sites A and B of the ferroxidase center of bullfrog M ferritin. Each site is modeled by the full structures of its three active amino acids. The potential surface at each site in the presence of water molecules is complex; coordination numbers of iron from three to six are seen. Water contributes to the complexity through its ability to hydrogen bond, to coordinate to iron, and to displace the neutral ligands glutamine and histidine. Intrinsic differences are noted at each site; at site B, the most stable complexes are found to favor tetracoordinate iron, while pentacoordination is preferred at site A in the two- and three-water complexes. While each incremental addition of a water molecule results in increased stability, successive binding energies are found to decrease.

Effects of cluster formation on spectra of benzo[a]pyrene and benzo[e]pyrene.

Absorption and fluorescence emission spectra of the polycyclic aromatic hydrocarbons benzo[a]pyrene (BaP) and benzo[e]pyrene (BeP) in solution and adsorbed on silica have been obtained and compared to examine the spectroscopic effects of clustering. Molecular mechanics calculations with the UFF potential were done to optimize monomer, dimer and trimer geometries, and energy differences were determined by MP2/6-31G* calculations. Fluorescence emission spectra of adsorbed BeP and BaP display a red shift that progresses with increased loading, and the two differ in their photodegradation kinetics. The experimental and theoretical results are found to be consistent.

High-spin versus broken symmetry-Effect of DFT spin density representation on the geometries of three diiron (III) model compounds.

Unrestricted density functional theory calculations have been conducted on three diiron(III) synthetic model compounds containing antiferromagnetically coupled high-spin (HS) irons for which crystallographic structures and Raman spectral data are available. Three density functionals have been employed: BPW91, PWC, and BOP. The study compares the effects on optimized geometries and harmonic vibrational frequencies of spin-paired (SP) low-spin, HS, and broken symmetry antiferromagnetically coupled singlet representations of the spin density distribution. The geometries around the diiron centers in the HS and broken symmetry (BS) representations are found to be similar, both markedly different from those arising from the SP representation. Small differences between the HS and BS results are seen in bond lengths, angles, Raman frequencies, and spin densities associated with oxo and peroxo bridges between the irons.

Computational study of iron(II) and -(III) complexes with a simple model human H ferritin ferroxidase center.

Interaction of iron ions with a six-amino acid model of the ferroxidase center of human H chain ferritin has been examined in density functional theory calculations. The model, based on experimental studies of oxidation of Fe2+ at the center, consists of Glu27, Glu62, His65, Glu107, Gln141, and Ala144. Reasonable structures are obtained in a survey of types of iron complexes inferred to occur in the ferroxidase reaction. Structures of complexes of the model center with one and two Fe2+ ions, with diiron(III) bridged by peroxide and bridged by oxide-peroxide combinations, have been optimized. Calculations on diiron(III) complexes confirm that stable peroxide-bridged complexes can form and that the Fe-Fe distance in at least one is consistent with the experimental Fe-Fe distance observed in the blue peroxodiferric complex of ferritin.

Sulfide-binding hemoglobins: Effects of mutations on active-site flexibility.

The dynamics of Hemoglobin I (HbI) from the clam Lucina pectinata, from wild-type sperm whale (SW) myoglobin, and from the L29F/H64Q/V68F triple mutant of SW, both unligated and bound to hydrogen sulfide (H2S), have been studied in molecular dynamics simulations. Features that account for differences in H2S affinity among the three have been examined. Our results verify the existence of an unusual heme rocking motion in unligated HbI that can promote the entrance of large ligands such as H2S. The FQF-mutant partially reproduces the amplitude and relative orientation of the motion of HbI's heme group. Therefore, besides introducing favorable electrostatic interactions with H2S, the three mutations in the distal pocket change the dynamic properties of the heme group. The active-site residues Gln-64(E7), Phe-43(CD1), and His-93(F8) are also shown to be more flexible in unligated HbI than in FQF-mutant and SW. Further contributions to H2S affinity come from differences in hydrogen bonding between the heme propionate groups and nearby amino acid residues.

Structures and energetics of Be(n)Si(n) and Be(2n)Si(n) (n = 1-4) clusters.

The structures and energies of Be(n)Si(n) and Be(2n)Si(n) (n = 1-4) clusters have been examined in ab initio theoretical electronic structure calculations. Cluster geometries have been established in B3LYP/6-31G(2df) calculations and accurate relative energies determined by the G3XMP2 method. The two atoms readily bond to each other and to other atoms of their own kind. The result is a great variety of low-energy clusters in a variety of structural types.