Manifolds of low energy structures for a magic number of hydrated sulfate: SO42-(H2O)24.

Low energy structures of SO42-(H2O)24 have been obtained using a combination of classical molecular dynamics simulations and refinement of structures and energies by quantum chemical calculations. Extensive exploration of the potential energy surface led to a number of low-energy structures, confirmed by accurate calibration calculations. An overall analysis of this large set was made after devising appropriate structural descriptors such as the numbers of cycles and their combinations. Low energy structures bear common motifs, the most prominent being fused cycles involving alternatively four and six water molecules. The latter adopt specific conformations which ensure the appropriate surface curvature to form a closed cage without dangling O-H bonds and at the same time provide 12-coordination of the sulfate ion. A prominent feature to take into account is isomerism via inversion of hydrogen bond orientations along cycles. This generates large families of ca. 100 isomers for this cluster size, spanning energy windows of 10-30 kJ mol-1. This relatively ignored isomerism must be taken into account to identify reliably the lowest energy minima. The overall picture is that the magic number cluster SO42-(H2O)24 does not correspond to formation of a single, remarkable structure, but rather to a manifold of structural families with similar stabilities. Extensive calculations on isomerization mechanisms within a family indicate that large barriers are associated to direct inversion of hydrogen bond networks. Possible implications of these results for magic number clusters of other anions are discussed.

Infrared Spectra of Deprotonated Dicarboxylic Acids: IRMPD Spectroscopy and Empirical Valence-Bond Modeling.

Experimental infrared multiple-photon dissociation (IRMPD) spectra recorded for a series of deprotonated dicarboxylic acids, HO2 (CH2 )n CO 2 - (n=2-4), are interpreted using a variety of computational methods. The broad bands centered near 1600 cm-1 can be reproduced neither by static vibrational calculations based on quantum chemistry nor by a dynamical description of individual structures using the many-body polarizable AMOEBA force field, strongly suggesting that these molecules experience dynamical proton sharing between the two carboxylic ends. To confirm this assumption, AMOEBA was combined with a two-state empirical valence-bond (EVB) model to allow for proton transfer in classical molecular dynamics simulations. Upon suitable parametrization based on ab initio reference data, the EVB-AMOEBA model satisfactorily reproduces the experimental infrared spectra, and the finite temperature dynamics reveals a significant amount of proton sharing in such systems.

Strategy for Modeling the Infrared Spectra of Ion-Containing Water Drops.

Hydrated ions are ubiquitous in environmental and biological media. Understanding the perturbation exerted by an ion on the water hydrogen bond network is possible in the nanodrop regime by recording vibrational spectra in the O-H bond stretching region. This has been achieved experimentally in recent years by forming gaseous ions containing tens to hundreds of water molecules and recording their infrared photodissociation spectra. In this paper, we demonstrate the capabilities of a modeling strategy based on an extension of the AMOEBA polarizable force field to implement water atomic charge fluctuations along with those of intramolecular structure along the dynamics. This supplementary flexibility of nonbonded interactions improves the description of the hydrogen bond network and, therefore, the spectroscopic response. Finite temperature IR spectra are obtained from molecular dynamics simulations by computing the Fourier transform of the dipole moment autocorrelation function. Simulations of 1-2 ns are required for extensive sampling in order to reproduce the experimental spectra. Furthermore, bands are assigned with the driven molecular dynamics approach. This method package is shown to compare successfully with experimental spectra for 11 ions in water drops containing 36-100 water molecules. In particular, band frequency shifts of the free O-H stretching modes at the cluster surface are well reproduced as a function of both ion charge and drop size.

Structure and thermodynamics of Mg:phosphate interactions in water: a simulation study.

The association of Mg(2+) and H2 PO4 (-) in water can give insights into Mg:phosphate interactions in general, which are very widespread, but for which experimental data is surprisingly sparse. It is studied through molecular dynamics simulations (>100 ns) by using the polarizable AMOEBA force field, and the association free energy is computed for the first time. Explicit consideration of outer-sphere and two types of inner-sphere association provides considerable insight into the dynamics and thermodynamics of ion pairing. After careful assessment of the computational approximations, the agreement with experimental values indicates that the methodology can be extended to other inorganic and biological Mg:phosphate interactions in solution.

Hydration of the sulfate dianion in cold nanodroplets: SO4(2-)(H2O)12 and SO4(2-)(H2O)13.

The structures, energetics and infrared spectra of SO4(2-)(H2O)12 and SO4(2-)(H2O)13 have been investigated by a combination of classical polarizable molecular dynamics and static quantum chemical calculations. Snapshots extracted from MD trajectories were used as inputs for local DFT optimization. Energies of the most stable structures were further refined at the ab initio level. A number of new low energy structures have thus been identified. The most stable structures of SO4(2-)(H2O)12 have the sulfate on the surface of the water cluster, while it may be slightly more burried in SO4(2-)(H2O)13, however still with an incomplete first hydration shell. Differences in the infrared spectra arise in part from mixing of sulfate stretching and water librational modes in the 900-1100 cm(-1) region, leading to some sensitivity of the IR spectrum to the structure. Second shell water molecules however do not generate signatures that are specific enough to relate spectra to structures straightforwardly, at least in this frequency range. Thus the emergence of a new band at 970 cm(-1) in the SO4(2-)(H2O)13 spectrum cannot be taken as a clue as to the number of water molecules which is necessary for a cluster to close the first hydration shell of sulfate. This number is at least 14 and possibly larger. However the density of low energy isomers is large enough that individual structures may loose meaning at all but the lowest temperatures.

Vibrational mode assignment of finite temperature infrared spectra using the AMOEBA polarizable force field.

The calculation of infrared spectra by molecular dynamics simulations based on the AMOEBA polarizable force field has recently been demonstrated [Semrouni et al., J. Chem. Theory Comput., 2014, 10, 3190]. While this approach allows access to temperature and anharmonicity effects, band assignment requires additional tools, which we describe in this paper. The Driven Molecular Dynamics approach, originally developed by Bowman, Kaledin et al. [Bowman et al. J. Chem. Phys., 2003, 119, 646, Kaledin et al. J. Chem. Phys., 2004, 121, 5646] has been adapted and associated with AMOEBA. Its advantages and limitations are described. The IR spectrum of the Ac-Phe-Ala-NH2 model peptide is analyzed in detail. In addition to differentiation of conformations by reproducing frequency shifts due to non-covalent interactions, DMD allows visualizing the temperature-dependent vibrational modes.

Interaction modes and absolute affinities of α-amino acids for Mn2+: a comprehensive picture.

Manganese is involved as a cofactor in the activation of numerous enzymes as well as the oxygen-evolving complex of photosystem II. Full understanding of the role played by the Mn(2+) ion requires detailed knowledge of the interaction modes and energies of manganese with its various environments, a knowledge that is far from complete. To bring detailed insight into the local interactions of Mn in metallopeptides and proteins, theoretical studies employing first-principles quantum mechanical calculations are carried out on [Mn-amino acid](2+) complexes involving all 20 natural α-amino acids (AAs). Detailed investigation of [Mn-serine](2+), [Mn-cysteine](2+), [Mn-phenylalanine](2+), [Mn-tyrosine](2+), and [Mn-tryptophan](2+) indicates that with an electron-rich side chain, the most stable species involves interaction of Mn(2+) with carbonyl oxygen, amino nitrogen, and an electron-rich section of the side chain of the AA in its canonical form. This is in sharp contrast with aliphatic side chains for which a salt bridge is formed. For aromatic AAs, complexation to manganese leads to partial oxidation as well as aromaticity reduction. Despite multisite binding, AAs do not generate strong enough ligand fields to switch the metal to a low- or even intermediate-spin ground state. The affinities of Mn(2+) for all AAs are reported at the B3LYP and CCSD(T) levels of theory, thereby providing the first complete series of affinities for a divalent metal ion. The trends are compared with those of other cations for which affinities of all AAs have been previously obtained.

Structure of sodiated polyglycines.

The intrinsic folding of peptides about a sodium ion has been investigated in detail by using infrared multiple photon dissociation (IRMPD) spectroscopy and a combination of theoretical methods. IRMPD spectroscopy was carried out on sodiated polyglycines G(n)-Na(+) (n=2-8), in both the fingerprint and N-H/O-H stretching regions. Interplay between experimental and computational approaches (classical and quantum) enables us to decipher most structural details. The most stable structures of the small peptides up to G(6)-Na(+) maximize metal-peptide interactions with all peptidic C=O groups bound to sodium. In addition, direct interactions between peptide termini are possible for G(6)-Na(+) and larger polyglycines. The increased flexibility of larger peptides leads to more complex folding and internal peptide structuration through γ or ß turns. A structural transition is found to occur between G(6)-Na(+) and G(7)-Na(+), leading to a structure with sodium coordination that becomes tri-dimensional for the latter. This transition was confirmed by H/D exchange experiments on G(n)-Na(+) (n=3-8). The most favorable hydrogen-bonding pattern in G(8)-Na(+) involves direct interactions between the peptide termini and opens the way to salt-bridge formation; however, there is only good agreement between experimental and computational data over the entire spectral range for the charge solvation isomer.

Structural, energetic and dynamical properties of sodiated oligoglycines: relevance of a polarizable force field.

Oligoglycine peptides (from two to ten residues) complexed to the sodium ion were studied by quantum chemical and molecular mechanics calculations to understand their structural and energetic properties. Modeling such systems required the use of a polarizable force field and AMOEBA, as developed by Ren and Ponder [J. Comput. Chem., 2002, 23, 1497], was chosen. Some electrostatic and torsional parameters were re-optimized using a rigorous procedure and validated against both geometric and energetic ab initio data in the gas phase. Molecular dynamics simulations were performed on seven sodiated octa-glycine (G(8)) structures. Structural transitions were generally observed (with the notable exception of the a-helix), leading to new structures that were further proved by ab initio calculations to be of low energies. The main result is that for G(8)-Na(+), there is a compromise between sodium peptide interactions and multiple hydrogen bonding. The accuracy achieved with AMOEBA demonstrates the potential of this force field for the realistic modeling of gaseous peptides.

Cu(II)-catalyzed reactions in ternary [Cu(AA)(AA - H)]+ complexes (AA = Gly, Ala, Val, Leu, Ile, t-Leu, Phe).

The unimolecular chemistry of [Cu(II)AA(AA - H)](+) complexes, composed of an intact and a deprotonated amino acid (AA) ligand, have been probed in the gas phase by tandem and multistage mass spectrometry in an electrospray ionization quadrupole ion trap mass spectrometer. The amino acids examined include Gly, Ala, Val, Leu, Ile, t-Leu and Phe. Upon collisionally-activated dissociation (CAD), the [Cu(II)AA(AA - H)](+) complexes undergo decarboxylation with simultaneous reduction of Cu(II) to Cu(I); during this process, a radical site is created at the alpha-carbon of the decarboxylated ligand (H(2)N(1) - (*)C(alpha)H - C(beta)H(2) - R; R = side chain substituent). The radical site is able to move along the backbone of the decarboxylated amino acid to form two new radicals (HN(1)(*) - C(alpha)H(2) - C(beta)H(2) - R and H(2)N(1) - C(alpha)H(2) - (*)C(beta)H - R). From the complexes of Gly and t-Leu, only C(alpha) and N(1) radicals can be formed. The whole radical ligand can be lost to form [Cu(I)AA](+) from these three isomeric radicals. Alternatively, further radical induced dissociations can take place along the backbone of the decarboxylated amino acid ligand to yield [Cu(II)AA(AA - 2H - CO(2))](+), [Cu(I)AA((*)NH(2))](+), [Cu(I)AA(HN = C(alpha)H(2))](+), or [Cu(I)AA(H(2)N - C(alpha)H = C(beta)H - R'](+) (R' = partial side chain substituent). The sodiated copper complexes, [Cu(II)(AA - H + Na)(AA - H)](+), show the same fragmentation patterns as their non-sodiated counterparts; sodium ion is retained on the intact amino acid ligand and is not involved in the CAD pathways. The amino groups of both AA units, the carbonyl group of the intact amino acid, and the deprotonated hydroxyl oxygen coordinate Cu(II) in square-planar fashion. Ab initio calculations indicate that the metal ion facilitates hydrogen atom shuttling between the N(1), C(alpha) and C(beta) atoms of the decarboxylated amino acid ligand. The dissociations of the decarboxylated radical ions unveil important insight about the so far largely unknown intrinsic chemistry of alpha-amino acid and peptide radicals, which are implicated as intermediates in numerous pathogenic biological processes.

Vibrational signatures of protonated, phosphorylated amino acids in the gas phase.

Structural characterization of protonated phosphorylated serine, threonine, and tyrosine was performed using mid-infrared multiple photon dissociation (IRMPD) spectroscopy and density functional theory (DFT) calculations. The ions were generated and analyzed by an external electrospray source coupled to a Paul ion-trap type mass spectrometer. Their fragmentation was induced by the resonant absorption of multiple photons from a tunable free electron laser (FEL) beam. IRMPD spectra were recorded in the 900-1850 cm(-1) energy range and compared to the corresponding computed IR spectra. On the basis of the frequency and intensity of two independent bands in the 900-1400 cm(-1) energy range, it is possible to identify the phosphorylated residue. IRMPD spectra for a 12-residue fragment of stathmin in its phosphorylated and nonphosphorylated forms were also recorded in the 800-1400 cm(-1) energy range. The lack of spectral congestion in the 900-1300 cm(-1) region makes their distinction facile. Our results show that IRMPD spectroscopy may became a valuable tool for structural characterization of small phosphorylated peptides.

A comparative study of semiempirical, ab initio, and DFT methods in evaluating metal-ligand bond strength, proton affinity, and interactions between first and second shell ligands in Zn-biomimetic complexes.

Although theoretical methods are now available which give very accurate results, often comparable to the experimental ones, modeling chemical or biological interesting systems often requires less demanding and less accurate theoretical methods, mainly due to computer limitations. Therefore, it is crucial to know the precision of such less reliable methods for relevant models and data. This has been done in this work for small zinc-active site models including O- (H(2)O and OH(-)) and N-donor (NH(3) and imidazole) ligands. Calculations using a number of quantum mechanical methods were carried out to determine their precision for geometries, coordination number relative stability, metal-ligand bond strengths, proton affinities, and interaction energies between first and second shell ligands. We have found that obtaining chemical accuracy can be as straightforward as HF geometry optimization with a double-zeta plus polarization basis followed by a B3LYP energy calculation with a triple-zeta quality basis set including diffuse and polarization functions. The use of levels as low as PM3 geometry optimization followed by a B3LYP single-point energy calculation with a double-zeta quality basis including polarization functions already yields useful trends in bond length, proton affinities or bond dissociation energies, provided that appropriate caution is taken with the optimized structures. The reliability of these levels of calculation has been successfully demonstrated for real biomimetic cases.

IRMPD spectroscopy of a protonated, phosphorylated dipeptide.

The protonated, phosphorylated dipeptide [GpY+H](+) is characterized by mid-infrared multiple-photon dissociation (IRMPD) spectroscopy and quantum-chemical calculations. The ions are generated in an external electrospray source and analyzed in a Fourier transform ion cyclotron resonance mass spectrometer, and their fragmentation is induced by resonant absorption of multiple photons emitted by a tunable free-electron laser. The IRMPD spectra are recorded in the 900-1730 cm(-1) range and compared to the absorption spectra computed for the lowest energy structures. A detailed calibration of computational levels, including B3LYP-D and coupled cluster, is carried out to obtain reliable relative energies of the low-energy conformers. It turns out that a single structure can be invoked to assign the IRMPD spectrum. Protonation at the N terminus leads to the formation of a strong ionic hydrogen bond with the phosphate P=O group in all low-energy structures. This leads to a P=O stretching frequency for [GpY+H](+) that is closer to that of [pS+H](+) than to that of [pY+H](+) and thus demonstrates the sensitivity of this mode to the phosphate environment. The COP phosphate ester stretching mode is confirmed to be an intrinsic diagnostic for identification of which type of amino acid is phosphorylated.

The sodium ion affinities of cytosine and its methylated derivatives.

The sodium ion affinities of cytosine (Cyt), 5-methylcytosine (5MeCyt) and 1-methylcytosine (1MeCyt) have been determined by experimental and quantum chemical methods. Na(+)-bound heterodimers were produced carrying one cytosine or methylated cytosine ligand (designated as C) and one peptide or amino acid reference base (designated as Pep); the Pep molecules included the peptides GlyLeu, GlyPhe, SerGly, and PheGly, and the amino acid His. The dissociation kinetics of these C--Na(+)--Pep ions were determined by collisionally activated dissociation (CAD) and converted to relative and absolute Na(+) affinities via kinetic method approaches. Relative Na(+) affinities increase in the order (kJ/mol): GlyLeu (0) < Cyt (3) < GlyPhe (4) < SerGly (6) < 5MeCyt (8) < PheGly (11) < 1MeCyt (13) < His (17). Anchoring the relative values of the nucleobases to the absolute affinities of the reference bases leads to absolute Na(+) affinities of 214 +/- 8, 219 +/- 8, and 224 +/- 8 kJ/mol for Cyt, 5MeCyt, and 1MeCyt, respectively. Ab initio calculations were used to confirm these results. The computed affinities of Cyt (213 kJ/mol) and 1MeCyt (217 kJ/mol) are in very good agreement with the experiments. These values unambiguously correspond to Na(+) complexes with the keto form of cytosine and its methyl derivatives. Ab initio calculations on tautomerization mechanisms in the gas versus condensed phase are used to discuss why the sodiated keto isomers were formed in the present electrospray ionization (ESI) experiments, but the enol isomers in previous fast atom bombardment (FAB) experiments.

[Nepsilon-(gamma-glutamyl) lysine] as a potential biomarker in neurological diseases: new detection method and fragmentation pathways.

Protein aggregates are characteristic of a number of diseases of the central nervous system such as diseases of polyQ expansion. Covalent bonds formed by the action of transglutaminase are thought to participate in the stabilization of these aggregates. Transglutaminase catalyzes the formation of cross-links between the side chains of glutaminyl and lysyl residues of polypeptides. Identification of the isodipeptide N(epsilon)-(gamma-glutamyl) lysine (iEK) in terminal proteolytic digests of neuronal aggregates would demonstrate participation of transglutaminase in neurological diseases. In order to identify and quantify the iEK present in the brain of patients with neurological disease, a method combining liquid chromatography and multistep mass spectrometry was developed. Because isobaric peptides of iEK could be present in the digest of aggregated proteins, the choice of fragment diagnostic ions was crucial. These ions were identified by mass spectrometry on sodiated iEK, which was derivatized on the carboxylic functions and terminal amines in order to improve sensitivity. Deuterated molecules as well as (13)C(6)- and (15)N(2)-isotopomers were used to derive filiations in the multistep fragmentations. The main fragmentation patterns have been identified, so that two ions (m/z 396 [MH - 56-42 u](+) and 350 [MH - 56-88 u](+)) are shown to be adequate markers for quantitation experiments. In order to gain a better understanding of the fragmentation processes, detailed quantum chemical calculations have been performed at levels which are expected to provide good accuracy. A thorough study has been carried out with a reduced model in which only the 'active' part of the molecule is retained. This allowed obtaining full mechanistic details on the pathways leading to a number of observed fragments. In particular, it has been shown that losses of 87 and 88 u from A(+) = [MH - 56 u](+) are competitive. Computations on the entire derivatized isodipeptide have been used to validate the use of the smaller model in order to obtain reliable energetics and mechanisms.

The alkylation mechanism of zinc-bound thiolates depends upon the zinc ligands.

Alkylation of zinc-bound thiolates occurs in both catalytic and structural zinc sites of enzymes. Recent biomimetic studies have led to a controversy as to which mechanism is operative in thiolate alkylation. Building on one of these biomimetic complexes, we have devised a series of models that allow for an appraisal of the roles of charge, ligand nature, and hydrogen bonding to sulfur on reactivity. The reactions of these complexes with methyl iodide, leading to thioethers and zinc iodide complexes, have been examined by density functional theory calculations, in the gas phase as well as in an aqueous solution. In all cases, a S(N)2 reaction is favored over sigma-bond metathesis. Both the net electronic charge and the hydrogen bond play a significant role in the nucleophilicity of the thiolate. We find that the mechanistic diversity observed experimentally can be explained by the difference in the net charge of the complexes. A dianionic complex follows a dissociative pathway, whereas an associative one is preferred for a neutral system.

The sodium ion affinities of simple di-, tri-, and tetrapeptides.

The sodium ion affinities (binding energies) of nineteen peptides containing 2-4 residues have been determined by experimental and computational approaches. Na(+)-bound heterodimers with amino acid and peptide ligands (Pep(1), Pep(2)) were produced by electrospray ionization. The dissociations of these Pep(1)-Na(+)-Pep(2) ions to Pep(1)-Na(+) and Pep(2)-Na(+) were examined by collisionally activated dissociation to construct a ladder of relative affinities via the kinetic method. The accuracy of this ladder was subsequently ascertained by experiments using several excitation energies for four peptide pairs. The relative scale was converted to absolute affinities by anchoring the relative values to the known Na(+) affinity of GlyGly. The Na(+) affinities of AlaAla, HisGly, GlyHis, GlyGlyGly, AlaAlaAla, GlyGlyGlyGly, and AlaAlaAlaAla were also calculated at the MP2(full)/6-311 + G(2d,2p) level of ab initio theory using geometries that were optimized at the MP2(full)/6-31G(d) level for AlaAla or HF/6-31G(d) level for the other peptides; the resulting values agree well with experimental Na(+) affinities. Increasing the peptide size is found to dramatically augment the Na(+) binding energy. The calculations show that in nearly all cases, all available carbonyl oxygens are sodium binding sites in the most stable structures. Whenever side chains are available, as in HisGly and GlyHis, specific additional binding sites are provided to the cation. Oligoglycines and oligoalanines have similar binding modes for the di- and tripeptides, but differ significantly for the tetrapeptides: while the lowest energy structure of GlyGlyGlyGly-Na(+) has the peptide folded around the ion with all four carbonyl oxygens in close contact with Na(+), that of AlaAlaAlaAla-Na(+) involves a pseudo-cyclic peptide in which the C and N termini interact via hydrogen bonding, while Na(+) sits on top of the oxygens of three nearly parallel C=O bonds.

Electrolytic electrospray ionization mass spectrometry of quaterthiophene-bridged bisporphyrins: beyond the identification tool.

Several quaterthiophene-bridged bisporphyrins were analyzed by electrospray ionization mass spectrometry (ESI-MS). The active centers of these molecular assemblies are two porphyrins moieties complexed (Z) or not (H) with a metal ion, typically Zn(2+), and the spacer is a quaterthiophene. The two end-groups were chemically linked to the quaterthiophene spacers by (i) a C--C single bond, (ii) a trans double bond or (iii) a triple bond. The formation of charged species either by protonation ([M + H](+) and [M + 2H](2+)) or electron(s) loss (M(+) and M(2+)), account for the occurrence of electrochemical processes in the basic operation of an electrospray source acting in a non-aqueous solvent. The nature of the observed charged species is correlated with the electro-oxidation properties and proton production by electro-oxidation of residual water. The occurence of these electrochemical reaction is proposed when the electroactivity of the electrosprayed substrates is not sufficient to support the current demand of the ESI source. In this way, the results obtained from the analysed series suggest the occurrence of such a process when the interfacial potential of the metal capillary reaches a value of 0.75 V vs Ag/AgCl. The results of theoretical calculations confirm the importance of the ionization energy with regard to the protonation energy in the course of the ionization reaction. The structural differences at the porphyrin-linker junctions lead to significantly smaller ionization energy in the case of the trans double bond. The MS observation of discharged dimers from molecular assemblies, including two complexed porphyrins ZZ or two free bases HH as end-group and a triple bond as the quaterthiophene-bisporphyrin junction, indicates together with molecular modelling (carried out at the semi-empirical PM3 level), that the planar and symmetric structures favour stacking.