Hydration of l-tyrosine in aqueous medium. An experimental and theoretical study by mixed quantum mechanical/molecular mechanics methods

Quantum mechanical/molecular mechanics (QM/MM) calculations were carried out in order to study the theoretical structures of l-tyrosine in both gas phase and in aqueous solution and observe the changes that occur on the structural and vibrational properties in two phases. Therefore, the molecule was characterized by infrared and Raman spectroscopy in solid phase and aqueous solution. Optimized geometries and relative stabilities for the zwitterion l-tyrosine derivatives have been calculated taking into account the solvent effects by using the self-consistent reaction field (SCRF) theory. For a complete assignment of the IR and Raman spectra of l-tyrosine in solid and aqueous solution phases, density functional theory (DFT) calculations were combined with Pulay's scaled quantum mechanical force field (SQMFF) methodology in order to fit the theoretical wavenumber values to the experimental ones. A good agreement between theoretical and available experimental results is found.     

Ab initio quantum mechanical/molecular mechanical simulation of electron transfer process: fractional electron approach

Electron transfer (ET) reactions are one of the most important processes in chemistry and biology. Because of the quantum nature of the processes and the complicated roles of the solvent, theoretical study of ET processes is challenging. To simulate ET processes at the electronic level, we have developed an efficient density functional theory (DFT) quantum mechanical (QM)/molecular mechanical (MM) approach that uses the fractional number of electrons as the order parameter to calculate the redox free energy of ET reactions in solution. We applied this method to study the ET reactions of the aqueous metal complexes Fe(H(2)O)(6)(2+/3+) and Ru(H(2)O)(6)(2+/3+). The calculated oxidation potentials, 5.82 eV for Fe(II/III) and 5.14 eV for Ru(II/III), agree well with the experimental data, 5.50 and 4.96 eV, for iron and ruthenium, respectively. Furthermore, we have constructed the diabatic free energy surfaces from histogram analysis based on the molecular dynamics trajectories. The resulting reorganization energy and the diabatic activation energy also show good agreement with experimental data. Our calculations show that using the fractional number of electrons (FNE) as the order parameter in the thermodynamic integration process leads to efficient sampling and validate the ab initio QM/MM approach in the calculation of redox free energies.     

Experimental evidence for a variable first coordination shell of the cadmium(II) ion in aqueous, dimethyl sulfoxide, and N,N'-dimethylpropyleneurea solution

A combined extended X-ray absorption fine structure (EXAFS) and large angle X-ray scattering (LAXS) investigation has been performed to evaluate the coordination structure of the cadmium(II) ion in aqueous, dimethyl sulfoxide, and N,N'-dimethylpropyleneurea (dmpu) solutions. This approach has singled out the existence of a flexible coordination shell around the cadmium(II) ion in aqueous and dimethyl sulfoxide solutions, whereas a regular octahedral complex is detected in dmpu. The EXAFS and LAXS techniques provide different values of the Cd-O first shell distance (2.27(1) A and 2.302(5) A, respectively) for the hydrated and dimethyl sulfoxide solvated complexes, and this discrepancy is originated by the simultaneous presence of hexa- and heptacoordinated complexes in solution, giving rise to a broad distribution of distances around the ion. These findings demonstrate that, in solution, the cadmium(II) ion forms quite flexible hydration and dimethyl sulfoxide solvate complexes undergoing a solvent exchange with unusually stable seven-coordinated intermediate complexes, and therefore the mean ion-solvent distance is longer in solution than in the solid state. In the dmpu solution, due to the bulkiness of the solvent molecules, the octahedral cadmium(II) solvate is extremely crowded and it is not possible for a seventh ligand to enter the inner-coordination shell. This investigation shows that the combined analysis of the EXAFS and LAXS data allows a reliable determination of the structural properties of electrolyte solutions, also in the presence of flexible coordination shell with a variable number of coordinating molecules.     

Structure of Zinc and Nickel Histidine Complexes in Solution Revealed by Molecular Dynamics and Raman Optical Activity

The histidine residue has an exceptional affinity for metals, but solution structure of its complexes are difficult to study. For zinc and nickel complexes, Raman and Raman optical activity (ROA) spectroscopy methods to investigate the link between spectral shapes and the geometry were used. The spectra were recorded and interpreted on the basis of ionic equilibria, molecular dynamics, ab initio molecular dynamics, and density functional theory. For zwitterionic histidine the dominant tautomer was determined by the decomposition of experimental spectra into calculated subspectra. An octahedral structure was found to prevail for the ZnHis₂ complex in solution, in contrast to a tetrahedral arrangement in the crystal phase. The solution geometry of NiHis₂ is more similar to the octahedral structure found by X-ray. The Raman and ROA structural determinations of metal complexes are dependent on extensive computations, but reveal unique information about the studied systems.     

A combined Raman, DFT and MD study of the solvation dynamics and the adsorption process of pyridine in silver hydrosols

The adsorption of pyridine onto silver colloids has been investigated by Raman spectroscopy experiments and by ab initio DFT and MP2 calculations. The solvation dynamics of the pyridine in water has been studied by a molecular dynamics simulation. The results are compared with the latest available experimental and theoretical data. It is found that the pyridine is essentially hydrogen bonded to one solvent molecule. Calculations based on pyridine-water and pyridine-Ag(+) complexes allow the reproduction of the experimentally observed Raman features and explain the adsorption process of the ligand in silver hydrosols.     

Computational study of the Zr4+ tetranuclear polymer, [Zr4(OH)8(H2O)16]8+

The Zr(4+) tetramer, [Zr(4)(OH)(8)(H(2)O)(16)](8+), is thought to be the major component of the Zr(4+) polymer system in aqueous solution, present as a dominant ionic cluster species compared to other Zr(4+) clusters under various experimental conditions. Despite widespread applications of zirconium, the structure and dynamics of the tetramer in aqueous solution are not well understood. We conducted a combination of ab initio molecular dynamics and quantum mechanical studies in the gas phase and aqueous solution and related our results to the available experimental data to provide atom-level information on the behavior of this species in aqueous solution. Our simulations indicate that the tetramer structure is stable on the picosecond time scale in an aqueous environment and that it is of a planar form, comprising eight-coordinated Zr(4+) ions with an antiprism/irregular dodecahedron ligand arrangement. In combination with our studies of Zr(4+) dimer and trimer clusters, our results provide detailed geometrical information on structural motifs for building zirconium polymers and suggest a possible polymerization path.     

DFT study on Raman frequencies of molten lead (II) chloride

The study concerns ab inito calculations of the essential Raman frequencies of molten lead (II) chloride near the melting point. Modelling of topologically disordered lead (II) chloride was carried out within the framework of the density functional theory using the Perdew–Burke-Ernzerhof (PBE) functional and optimised basis sets. Calculations were performed for a cluster containing 24 formula units. The optimum geometry of the cluster was determined and the local structure parameters were found. Nano-size effect leads to the picture of damped oscillations on radial atomic density distributions typical of molten salts. Distorted octahedra of PbCl₆ were detected inside the cluster. Ab initio calculation of the Raman spectrum of the octahedral complexes inside the cluster structure was implemented. It was shown that the spectrum has a peak at 192 cm^?1, which agrees well with the experimental Raman spectra of lead (II) chloride melt near the melting point.     

Characterization of hydrated Na+(phenol) and K+(phenol) complexes using infrared spectroscopy

Hydrated alkali metal ion-phenol complexes were studied to model these species in aqueous solution for M=Na and K. IR predissociation spectroscopy in the O-H stretch region was used to analyze the structures of M+(Phenol)(H2O)n cluster ions, for n = 1-4. The onset of hydrogen bonding was observed to occur at n=4. Ab initio calculations were used to qualitatively explore the types of hydrogen-bonded structures of the M+(Phenol)(H2O)4 isomers. By combining the ab initio calculations and IR spectra, several different structures were identified for each metal ion. In contrast to benzene, detailed in a previous study of Na+(Benzene)n(H2O)m [J. Chem. Phys. 110, 8429 (1999)], phenol is able to bind directly to Na+ even in the presence of four waters. This is likely the result of the sigma-type interaction between the phenol oxygen and the ion. With K+, the dominant isomers are those in which the phenol O-H group is involved in a hydrogen bond with the water molecules, while with Na+, the dominant isomers are those in which the phenol O-H group is free and the water molecules are hydrogen-bonded to each other. Spectra and ab initio calculations for the M+(Phenol)Ar cluster ions for M=Na and K are reported to characterize the free phenol O-H stretch in the M+(Phenol) complex. While pi-type configurations were observed for binary M+(Phenol) complexes, sigma-type configurations appear to dominate the hydrated cluster ions.     

First principles and experimental 1H NMR signatures of solvated ions: The case of HCl(aq).

A combined experimental and ab initio study is presented of the 1H NMR chemical shift distribution of aqueous hydrogen chloride solution as a function of acid concentration, based on Car-Parrinello molecular dynamics simulations and fully periodic NMR chemical-shift calculations. The agreement of computed and experimental spectra is very good. From first-principles calculations, we can show that the individual contributions of Eigen and Zundel ions, regular water molecules, and the chlorine solvation shell to the NMR line are very distinct and almost independent of the acid concentration. From the computed instantaneous NMR distributions, it is further possible to characterize the average variation in hydrogen-bond strength of the different complexes.     

Toward a DFT-based molecular dynamics description of Co(II) binding in sulfur-rich peptides

In this paper, we investigated the reliability of a Car-Parrinello molecular dynamics (CPMD) approach to characterize the binding of Co(II) metal cation to peptide molecules containing cysteine. To this end, we compared pseudo-potentials and DFT plane wave expansion, which are used as key ingredients in the CPMD method, with standard all-electron Gaussian basis set DFT calculations. The simulations presented here are the first attempts to characterize interactions and dynamics of Co(II) metal with the building blocks of phytochelatin peptide molecules. Benchmark calculations are performed on [Co(Cys-H)]+ and [Co(Glutathione-H)]+ complexes, since they are the main fragments of the Co(II)-Cys and Co(II)-glutathione systems found in gas phase electrospray ionisation mass spectrometry (ESI-MS) experiments done in our laboratory. We also present benchmark calculations on the [Co(H2O)6)]2+ cluster with direct comparisons to highly correlated ab initio calculations and experiments. In particular, we investigated the dissociation path of one water molecule from the first hydration shell of Co(II) with CPMD. Overall, our molecular dynamics simulations shed some light on the nature of the Co(II) interaction and reactivity in Co(II)-phytochelatin building block systems related to the biological and environmental activity of the metal, either in the gas or liquid phase.     

Zinc(II) and cadmium(II) <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>'-Bis(2-<Emphasis Type="Italic">N</Emphasis>-Tosylaminobenzylidene) diaminodipropyliminates: Syntheses,structures, and photoluminescence properties

Chemical and electrochemical syntheses of the zinc(II) and cadmium(II) complexes with the tetradentate Schiff base (H₂L), the condensation product of 2-N-tosylaminobenzaldehyde with diaminodipropylamine, are carried out. The structures, compositions, and properties of the synthesized metal complexes are studied by elemental analyses, IR spectroscopy, ¹H NMR, UV spectroscopy, X-ray absorption spectroscopy, and quantum-chemical calculations. The structure of the cadmium(II) complex is determined by X-ray diffraction analysis (CIF file CCDC no. 1446393). The cadmium(II) and zinc(II) complexes exhibit luminescence in a CH₂Cl₂ solution in the blue spectral range (λ_PL = 425–428 nm) with the photoluminescence quantum yields ? = 0.20 and 0.75, respectively.     

Solvent effects on peroxynitrite structure and properties from QM/MM simulations

We have investigated the structure and the vibrational spectrum of peroxynitrite anion in aqueous solution by means of combined quantum-classical (QM/MM) molecular dynamics simulations. In our QM/MM scheme, the reactant was modeled using density functional theory with a Gaussian basis set and the solvent was described using the mean-field TIP4P and the polarizable TIP4P-FQ force fields. The choice of basis sets, functionals and force field parameters has been validated by performing calculations on isolated peroxynitrite and on small peroxynitrite-water complexes. Poor values for isolated peroxynitrite structural properties and vibrational frequencies are found for most ab initio methods, particularly regarding the ON-OO(-) bond distance and the harmonic stretching frequency, probably due to the singlet-triplet instability found in the HF wave function. On the other hand, DFT methods yield results in better agreement with high level CCSD(T) ab initio calculations. We have computed the vibrational spectrum for aqueous peroxynitrite by calculating the Fourier transform of the velocity autocorrelation function extracted from the QM-MM molecular dynamics simulations. Our computational scheme, which allows for the inclusion of both anharmonicity and solvent effects, is able to clarify previous discrepancies regarding the experimental spectra assignments and to shed light on the subtle interplay between solvation and peroxynitrite structure and properties.     

Solvent effects and dynamic averaging of 195Pt NMR shielding in cisplatin derivatives

The influences of solvent effects and dynamic averaging on the (195)Pt NMR shielding and chemical shifts of cisplatin and three cisplatin derivatives in aqueous solution were computed using explicit and implicit solvation models. Within the density functional theory framework, these simulations were carried out by combining ab initio molecular dynamics (aiMD) simulations for the phase space sampling with all-electron relativistic NMR shielding tensor calculations using the zeroth-order regular approximation. Structural analyses support the presence of a solvent-assisted "inverse" or "anionic" hydration previously observed in similar square-planar transition-metal complexes. Comparisons with computationally less demanding implicit solvent models show that error cancellation is ubiquitous when dealing with liquid-state NMR simulations. After aiMD averaging, the calculated chemical shifts for the four complexes are in good agreement with experiment, with relative deviations between theory and experiment of about 5% on average (1% of the Pt(II) chemical shift range).     

A resonance Raman,surface-enhanced resonance Raman,IR, and ab initio vibrational spectroscopic study of nickel(II) tetraazaannulene complexes

The IR and resonance Raman spectra of the nickel(II) complexes of dibenzo[b,i][1,4,8,11]tetraaza[14]annulene (TAA) and 5,7,12,14-tetramethyldibenzo[b,i][1,4,8,11]tetraaza[14]annulene (TMTAA) have been measured and compared with ab initio calculations of the vibrational wavenumbers at the B3-LYP level using the LanL2DZ basis set. An excellent fit is found between the experimental and calculated data, enabling precise vibrational assignments to be made. Surface-enhanced resonance Raman spectra were obtained following adsorption on Ag electrodes, with potentials in the range -0.1 to -1.1 V vs Ag/AgCl. There is evidence for contributions from both the electromagnetic and charge transfer (CT) surface enhancement mechanisms. The data indicate that variations in band intensities with electrode potential can be interpreted in terms of the CT mechanism.     

Potential energy surfaces and dynamics of Ni2+ ion aqueous solution: molecular dynamics simulation of the electronic absorption spectrum

We develop a model effective Hamiltonian for describing the electronic structures of first-row transition metals in aqueous solutions using a quasidegenerate perturbation theory. All the states consisting of 3d(n) electronic configurations are determined by diagonalizing a small effective Hamiltonian matrix, where various intermolecular interaction terms such as the electrostatic, polarization, exchange, charge transfer, and three-body interactions are effectively incorporated. This model Hamiltonian is applied to constructing the ground and triplet excited states potential energy functions of Ni(2+) in aqueous solution, based on the ab initio multiconfiguration quasidegenerate perturbation theory calculations. We perform molecular dynamics simulation calculations for the ground state of Ni(2+) aqueous solution to calculate the electronic absorption spectral shape as well as the ground state properties. Agreement between the simulation and experimental spectra is satisfactory, indicating that the present model can well describe the Ni(2+) excited state potential surfaces in aqueous solution.     

Hyperfine coupling tensors of the benzosemiquinone radical anion from Car-Parrinello molecular dynamics.

Based on Car-Parrinello ab initio molecular dynamics simulations of the benzosemiquinone radical anion in both aqueous solution and the gas phase, density functional calculations provide the currently most refined EPR hyperfine coupling (HFC) tensors of semiquinone nuclei and solvent protons. For snapshots taken at regular intervals from the molecular dynamics trajectories, cluster models with different criteria for inclusion of water molecules and an additional continuum solvent model are used to analyse the HFCs. These models provide a detailed picture of the effects of dynamics and of different intermolecular interactions on the spin-density distribution and HFC tensors. Comparison with static calculations allows an assessment of the importance of dynamical effects, and of error compensation in static DFT calculations. Solvent proton HFCs depend characteristically on the position relative to the semiquinone radical anion. A point-dipolar model works well for in-plane hydrogen-bonded protons but deviates from the quantum chemical values for out-of-plane hydrogen bonding.     

Perturbation of local solvent structure by a small dication: a theoretical study on structural, vibrational, and reactive properties of beryllium ion in water

A molecular based understanding of beryllium chemistry in the context of biomolecules is necessary for gaining progress in prevention and treatment of chronic beryllium disease. One aspect that has hindered the theoretical progress has been the lack of a simple classical two-body potential for the aqueous beryllium ion (Be2+) to be used with biomolecular simulations. We provide new parameters for Be2+ that capture the structural and reactive properties of this small dication. Using classical molecular dynamics simulations, we show that these parameters reproduce the correct radial distribution function and coordination numbers for this cation in explicit aqueous solution when compared to published diffraction and NMR measurements. The geometrical parameters obtained using classical simulations are also in agreement with ab initio calculations. We successfully predict the vibrational modes of the tetra aqua Be2+ dication from ab initio calculations on solvated structures obtained from the simulations. The calculated vibrational modes show better agreement with experiments compared to any published work. This new potential also produces a well-established hydrogen bonding between the first and second solvation shells. More importantly, when the molecular dynamics (MD) and ab initio results are interpreted in concert, the dynamics and nature of interactions between the first and second shells capture the pivotal role they play on the reactivity of aqua-Be complexes.     

Vibrational spectra and structures of the anions of urazole and 4-methylurazole: DFT calculations of the normal modes and the influence of hydrogen bonding

Solid state IR and Raman as well as aqueous solution state Raman spectra are reported for the anions of urazole and 4-methylurazole, and their N-deuterated derivatives. DFT calculations, at the B3-LYP/cc-pVTZ level, established that the structures and vibrational spectra of both anions can be interpreted using a model that incorporates hydrogen-bonded water molecules, in conjunction with the polarizable continuum solvation method. In the case of the urazole anion it is shown that deprotonation occurs primarily at N1 rather than N4, but there is also evidence for the second tautomer both in the solid state and in aqueous solution. The vibrational spectra were computed at the optimised molecular geometry in each case, enabling normal coordinate analysis, which yielded satisfactory agreement with the experimental IR and Raman data. Computed potential energy distributions of the normal modes provided detailed vibrational assignments.     

Infrared and ab initio studies of hydrogen bonding and proton transfer in the complexes formed by pyrazoles

The results of experimental studies and quantum mechanical calculations of vibrational spectra and structure of hydrogen bonded complexes formed by pyrazole (P) and 3,5-dimethylpyrazole (DMP) are presented. IR spectra of pyrazoles in solutions, gas phase, and solid state have been investigated in wide range of concentrations and temperatures. It has been found that in the gas phase both P and DMP reveal the equilibrium between monomers, dimers, and trimers. In solutions the equilibrium between monomers and trimers dominates, no bands, which can be attributed to dimers were detected. DMP retains the trimer structure in solid state, while in the case of pyrazole P, formation of the crystal provides another type of association. Geometrical and spectral characteristics of dimers and trimers, obtained by ab initio calculations, are presented and compared with experimental data.

IR spectra of solutions containing P and DMP with a number of acids (acetic and trifluoroacetic acids, pentachlorophenol, HBr) have been studied in parallel with ab initio calculations. It has been found that pentachlorophenol forms with pyrazoles complexes with one strong hydrogen bond O–HN, while NH pyrazole group remains unbonded. With carboxylic acids DMP forms 1:1 cyclic complexes with two hydrogen bonds. In the case of acetic acid, the complex in CH₂Cl₂ solution reveals molecular structure with OHN and C=OHN bonds, in accordance with results of the calculations. For trifluoroacetic acid, the calculations predict the molecular structure to be energetically more stable in the case of the isolated binary complex (in gas phase), while the experimental spectrum of CH₂Cl₂ solution gives an evidence of the proton transfer with formation of the cyclic ionic pair with two NH⁺O⁻ bonds. The agreement with experimental results can be improved by taking into account the influence of environment in the framework of Onsager or Tomasi models. The shape of proton potential function of the complexes and medium effect on its parameters, resulted from experimental data and calculations, are discussed. It has been found that the number of potential minima and their relative depth depend strongly on the method of calculations and the basic set. Under excess of trifluoroacetic acid, the formation of 2:1 acid–DMP complex has been detected. Spectral characteristics and results of calculations point to the cyclic structure of this complex, which includes homoconjugated bis-trifluoroacetate anion and DMPH⁺ cation. With HBr both studied pyrazoles were found to form ionic complexes including one or two pyrazole molecules per one acid molecule and correspondingly monocation or homoconjugated cation BHB⁺.     

Fine speciation of active sites in zeolites by a CO probe: Dynamics and IR frequencies

The influence of system dynamics onto infrared spectra calculated and measured for probe molecules weakly interacting with active sites is discussed. OH stretching frequencies for Brønsted sites in Y zeolite and their shift on CO adsorption (measured by IR spectroscopy at ca. 170 K) are simulated from AIMD trajectory, based on BO DFT calculations. We present here very good agreement between experimental and computed OH frequency red‐shift on CO adsorption in dealuminated zeolite Y. It is shown that extracting nonharmonic OH frequencies from ab initio molecular dynamics trajectory provides better estimate of experimental IR results than standard vibrational analyses due to realistic computational model and better correspondence with experimental conditions.