Charge transfer in FeO: a combined molecular-dynamics and ab initio study

Molecular-dynamics simulations and ab initio electronic structure calculations were carried out to determine the rate of charge transfer in stoichiometric wustite (FeO). The charge transfer of interest occurs by II/III valence interchange between nearest-neighbor Fe atoms, with the Fe(III) constituting a "hole" electronic defect. There are two possible nearest-neighbor charge transfers in the FeO lattice, which occur between edge-sharing or corner-sharing FeO(6) octahedra. Molecular-dynamics simulations predict charge-transfer rates of 3.7 x 10(11) and 1.9 x 10(9) s(-1) for the edge and corner transfers, respectively, in good agreement with those calculated using an ab initio cluster approach (1.6 x 10(11) and 8.0 x 10(8) s(-1), respectively). The calculated rates are also similar to those along the basal and c-axis directions in hematite (alpha-Fe(2)O(3)) determined previously. Therefore, as is the case for hematite, wustite is predicted to show anisotropic electrical conductivity. Our findings indicate that a rigid-ion model does not give acceptable results, thus showing the need to account for the change in polarizability of the system upon charge transfer. Our model achieves this by using a simple mechanical shell model. By calculating the electronic coupling matrix elements for many transition state configurations obtained from the molecular-dynamics simulations, we found evidence that the position of the bridging oxygen atoms can greatly affect the amount of electronic coupling between the donor and acceptor states. Finally, we address the effect of oxygen vacancies on the charge transfer. It was found that an oxygen vacancy not only creates a driving force for holes to transport away from the vacancy (or equivalently for electrons to diffuse toward the vacancy) but also lowers the free-energy barriers for charge transfer. In addition, the reorganization energy significantly differed from the nondefective case in a small radius around the defect.     

Adduct formation in electrospray ionization mass spectrometry II. Benzoic acid derivatives

This work serves as a follow-up to Part I of experiments designed to determine the underlying principles in the formation of pseudomolecular, or adduct, ions during electrospray ionization. Aromatic acids were studied by flow injection analysis in the negative ionization mode of electrospray ionization mass spectrometry. Part I dealt with common acidic anti-inflammatory pharmaceuticals. such as ibuprofen and related analogues. Part II deals with functionally less complex molecules, namely benzoic acid (BA) and substituted benzoic acids. Halide-substituted molecules are investigated to deduce the effect of electron-withdrawing substituents (bromo-, chloro-, and fluoro-) and ring position (ortho-, meta- and para-) on the response of a traditional deprotonated molecular ion ([M-H]-) and a sodium-bridged dimer ion ([2M-2 H+Na]-). Amino-substituted benzoic acids are also analyzed in order to study the effect of an additional ionizable group on the molecule, and para-tert.-butyl-BA is analyzed to study the effect of increased hydrophobicity, as they relate to the formation of pseudomolecular ions. This study shows that solution character [octanol-water partition coefficient (or log P) and pKa] of the model compounds controls the relative efficiency of formation of [M-H]- and [2M-2H+Na]- ions. However the relative gas phase character (gas phase basicity and proton affinity) also has a significant effect on the formation of the sodium-bridged dimer ion. For the halide-substituted species, placement of the electron-withdrawing atom at the meta-position gives the greatest enhancement in sensitivity. Observations also show that as the structural complexity of the model compound increases, predictions relating analyte acidity to sodium-bridged dimer ion formation give way to a stronger dependence between log P values and ionization efficiency. Supporting this hypothesis is the nearly ten-fold enhancement in signal for tert.-butyl BA relative to BA. due to the greater hydrophobicity, and consequently, increased surface activity in an electrosprayed droplet of the analyte molecule.     

Effect of amine ligand bulk on the interaction of methionine with platinum(II) diamine complexes

Molecular mechanics and dynamics calculations have been used in conjunction with experimental data to study the effects of amine ligand bulk on the formation of both guanine and methionine complexes with platinum diamine compounds. The AMBER force field has been supplemented with previous modifications [Yao; et al. Inorg. Chem. 1994, 33, 6061-6077. Cerasino; et al. Inorg. Chem. 1997, 36, 6070-6079] and has been further modified to include parameters for platinum bound to the sulfur atom of methionine. Molecular mechanics calculations with this modified AMBER force field have suggested that a platinum complex with two sulfur-bound methionine ligands and a bulky diamine ligand (N,N,N',N'-tetramethylethylenediamine, Me(4)en) would have severe interligand clashes; such interligand clashes are less pronounced in bis(9-ethylguanine) complexes. Consistent with these observations, NMR studies with [Pt(Me(4)en)(D(2)O)(2)](2+) have indicated that guanine 5'-monophosphate reacts in a 2:1 guanine:platinum ratio while both methionine and N-acetylmethionine react with only a 1:1 stoichiometry. Methionine forms a chelate via the sulfur and nitrogen atoms whereas N-acetylmethionine forms a chelate via the sulfur and oxygen atoms. The oxygen of the latter chelate can be displaced by the addition of guanosine 5'-monophosphate, although complete displacement of the N-acetylmethionine was not observed.     

Hydroxyl radical probe of the calmodulin-melittin complex interface by electrospray ionization mass spectrometry

The calcium-dependent interaction of calmodulin and melittin is studied through the application of a radical probe approach in which solutions of the protein and peptide and protein alone are subjected to high fluxes of hydroxyl and other oxygen radicals on millisecond timescales. These radicals are generated by an electrical discharge within an electrospray ion source of a mass spectrometer. Condensation of the electrosprayed droplets followed by proteolytic digestion of both calmodulin and melittin has identified residues in both which participate in the interaction and/or are shielded from solvent within the protein complex. Consistent with other theoretical models and available experimental data, the tryptophan residue of melittin at position 19 is shown to be critical to the formation of the complex with the C-terminal domain of peptide enveloped by and protected from oxidation upon binding to the protein. Furthermore, the N-terminal domain (to residue 36) and tyrosine at position 99 in calmodulin are significantly protected from limited oxidation upon the binding of melittin while exposing the phenylalanine residue at position 92 of the flexible loop domain. The N-terminus (through residue 36) of calmodulin is shown to lie in closer proximity to the melittin helix than its C-terminal counterpart (residues 127-148) based upon the protection levels measured at reactive residues within these segments of the protein.     

Rydberg-London potential for diatomic molecules and unbonded atom pairs

We propose and test a pair potential that is accurate at all relevant distances and simple enough for use in large-scale computer simulations. A combination of the Rydberg potential from spectroscopy and the London inverse-sixth-power energy, the proposed form fits spectroscopically determined potentials better than the Morse, Varnshi, and Hulburt-Hirschfelder potentials and much better than the Lennard-Jones and harmonic potentials. At long distances, it goes smoothly to the London force appropriate for gases and preserves van der Waals's "continuity of the gas and liquid states," which is routinely violated by coefficients assigned to the Lennard-Jones 6-12 form.     

Using molecular dynamics simulations to identify the key factors responsible for chiral recognition by an amino acid-based molecular micelle

Molecular dynamics (MD) simulations were used to investigate the binding of six chiral compounds to the amino acid-based molecular micelle (MM) poly-(sodium undecyl-(L)-leucine-leucine) or poly(SULL). The MM investigated is used as a chiral selector in capillary electrophoresis. The project goal was to characterize the chiral recognition mechanism in these separations and to move toward predictive models to identify the best amino acid-based MM for a given separation. Poly(SULL) was found to contain six binding sites into which chiral compounds could insert. Four sites had similar sizes, shapes, and electrostatic properties. Enantiomers of alprenolol, propranolol, 1,1′-bi-2-naphthyl-2,2′-diyl hydrogen phosphate, 1,1′-bi-2-naphthol, chlorthalidone, or lorazepam were separately docked into each binding pocket and MD simulations with the resulting intermolecular complexes were performed. Solvent-accessible surface area calculations showed the compounds preferentially associated with binding sites where they penetrated into the MM core and shielded their non-polar atoms from solvent. Furthermore, with five of the six compounds the enantiomer with the most favorable free energy of MM association also experienced the most favorable intermolecular hydrogen bonding interactions with the MM. This result suggests that stereoselective intermolecular hydrogen bonds play an important role in chiral discrimination in separations using amino acid-based MMs.GRAPHICAL ABSTRACT     

Comparison of solid-phase extraction methods for the determination of azaspiracids in shellfish by liquid chromatography-electrospray mass spectrometry

Azaspiracids have been identified as the cause of a new toxic syndrome called azaspiracid poisoning (AZP) that has led to incidents of human intoxications throughout Europe following the consumption of mussels. Although five AZP toxins have been structurally elucidated to-date, azaspiracid (AZA1), 8-methylazaspiracid (AZA2) and 22-demethylazaspiracid (AZA3) are the predominant toxins. Separation of the three main AZP toxins was achieved using reversed-phase liquid chromatography (LC) and coupled to an electrospray ionisation source of an ion-trap mass spectrometer. Five reversed-phase (C18) and three diol solid-phase extraction (SPE) cartridges were compared for their efficacy in the cleanup of shellfish matrix. The comparison was based on the optimum recoveries of AZA1, AZA2 and AZA3 from extracts of mussel tissues. LC-electrospray MS3 analysis was used to quantify the AZP toxins in wash and eluate fractions in the SPE studies. Good recovery and reproducibility data were obtained for one diol SPE cartridge and two C18 SPE cartridge types.     

Molecular dynamics investigation of ferrous-ferric electron transfer in a hydrolyzing aqueous solution: calculation of the pH dependence of the diabatic transfer barrier and the potential of mean force

We present a molecular model for ferrous-ferric electron transfer in an aqueous solution that accounts for electronic polarizability and exhibits spontaneous cation hydrolysis. An extended Lagrangian technique is introduced for carrying out calculations of electron-transfer barriers in polarizable systems. The model predicts that the diabatic barrier to electron transfer increases with increasing pH, due to stabilization of the Fe3+ by fluctuations in the number of hydroxide ions in its first coordination sphere, in much the same way as the barrier would increase with increasing dielectric constant in the Marcus theory. We have also calculated the effect of pH on the potential of mean force between two hydrolyzing ions in aqueous solution. As expected, increasing pH reduces the potential of mean force between the ferrous and ferric ions in the model system. The magnitudes of the predicted increase in diabatic transfer barrier and the predicted decrease in the potential of mean force nearly cancel each other at the canonical transfer distance of 0.55 nm. Even though hydrolysis is allowed in our calculations, the distribution of reorganization energies has only one maximum and is Gaussian to an excellent approximation, giving a harmonic free energy surface in the reorganization energy F(DeltaE) with a single minimum. There is thus a surprising amount of overlap in electron-transfer reorganization energies for Fe(2+)-Fe(H2O)6(3+), Fe(2+)-Fe(OH)(H2O)5(2+), and Fe(2+)-Fe(OH)2(H2O)+ couples, indicating that fluctuations in hydrolysis state can be viewed on a continuum with other solvent contributions to the reorganization energy. There appears to be little justification for thinking of the transfer rate as arising from the contributions of different hydrolysis states. Electronic structure calculations indicate that Fe(H2O)6(2+)-Fe(OH)n(H2O)(6-n)(3-n)+ complexes interacting through H3O2- bridges do not have large electronic couplings.     

An addition to the oxoacid family: H2B12(OH)12

The acid H(2)B(12)(OH)(12) can be isolated as a crystalline solid by protonation of the hydroxylated borane anion, B(12)(OH)(12)(2)(-). This acidic compound has low solubility in water, conducts protons in the solid state, and has thermal stability to a temperature of 400 degrees C. The conductivity mechanism is a Grotthuss mechanism with a low activation enthalpy (9-13 kcal/mol). This new acid represents an addition to the class of oxoacids, of which sulfuric and phosphoric acid are the most prominent examples.