Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

Solvatochromism in Binary Mixtures: First Report on a Solvation Free Energy Relationship between Solvent Exchange Equilibrium Constants and the Properties of the Medium

We have employed UV-vis spectroscopy in order to investigate details of the solvation of six solvatochromic indicators, hereafter designated as ""probes"", namely, 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr,, respectively. These can be divided into three pairs, each includes two probes of similar pK(a) in water and different lipophilicity. Solvation has been studied in binary mixtures, BMs, of water, W, with 12 protic organic solvents, S, including mono- and bifunctional alcohols (2-alkoxyethanoles, unsaturated and chlorinated alcohols). Each medium was treated as a mixture of S, W, and a complex solvent, S-W, formed by hydrogen bonding. Values of lambda(max) (of the probe intramolecular charge transfer) were converted into empirical polarity scales, E(T)(probe) in kcal/mol, whose values were correlated with the effective mole fraction of water in the medium, chi w(effective). This correlation furnished three equilibrium constants for the exchange of solvents in the probe solvation shell; phi(W/S) (W substitutes S): phi(S-W/W) (S-W substitutes W), and phi(S-W/S) (S-W substitutes S), respectively. The values of these constants depend on the physicochemical properties of the probe and the medium. We tested, for the first time, the applicability of a new solvation free energy relationship: phi = constant + a alpha(BM) + b beta(BM) + s(pi*(BM) + d delta) + p log P(BM), where a, b, s, and p are regression coefficients alpha(BM), beta(BM), and pi*(BM) are solvatochromic parameters of the BM, delta is a correction term for pi*, and log P is an empirical scale of lipophilicity. Correlations were carried out with two-, three-, and four-medium descriptors. In all cases, three descriptors gave satisfactory correlations; use of four parameters gave only a marginal increase of the goodness of fit. For phi(W/S), the most important descriptor was found to be the lipophilicity of the medium; for phi(S-W/W) and phi(S-W/S), solvent basicity is either statistically relevant or is the most important descriptor. These responses are different from those of E(T)(probe) of many solvatochromic indicators in pure solvents, where the importance of solvent basicity is usually marginal, and can be neglected.

Solvation in Pure Liquids: What Can Be Learned from the Use of Pairs of Indicators?

The solvation of six solvatochromic probes in a large number of solvents (33-68) was examined at 25 degrees C. The probes employed were the following: 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenyl pyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), respectively. Of these, MePMBr is a novel compound. They can be grouped in three pairs, each with similar pK(a) in water but with different molecular properties, for example, lipophilicity and dipole moment. These pairs are formed by RB and MePM; QB and MePMBr; WB and MePMBr(2), respectively. Theoretical calculations were carried out in order to calculate their physicochemical properties including bond lengths, dihedral angles, dipole moments, and wavelength of absorption of the intramolecular charge-transfer band in four solvents, water, methanol, acetone, and DMSO, respectively. The data calculated were in excellent agreement with available experimental data, for example, bond length and dihedral angles. This gives credence to the use of the calculated properties in explaining the solvatochromic behaviors observed. The dependence of an empirical solvent polarity scale E(T)(probe) in kcal/mol on the physicochemical properties of the solvent (acidity, basicity, and dipolarity/polarizability) and those of the probes (pK(a), and dipole moment) was analyzed by using known multiparameter solvation equations. For each pair of probes, values of E(T)(probe) (for example, E(T)(MePM) versus E(T)(RB)) were found to be linearly correlated with correlation coefficients, r, between 0.9548 and 0.9860. For the mercyanine series, the values of E(T)(probe) also correlated linearly, with (r) of 0.9772 (MePMBr versus MePM) and 0.9919 (MePMBr(2) versus MePM). The response of each pair of probes (of similar pK(a)) to solvent acidity is the same, provided that solute-solvent hydrogen-bonding is not seriously affected by steric crowding (as in case of RB). We show, for the first time, that the response to solvent dipolarity/polarizability is linearly correlated to the dipole moment of the probes. The successive introduction of bromine atoms in MePM (to give MePMBr, then MePMBr(2)) leads to the following linear decrease: pK(a) in water, length of the phenolate oxygen-carbon bond, length of the central ethylenic bond, susceptibility to solvent acidity, and susceptibility to solvent dipolarity/polarizability. Thus studying the solvation of probes whose molecular structures are varied systematically produces a wealth of information on the effect of solute structure on its solvation. The results of solvation of the present probes were employed in order to test the goodness of fit of two independent sets of solvent solvatochromic parameters.

First study on the thermo-solvatochromism in aqueous 1-(1-butyl)-3-methylimidazolium tetrafluoroborate: A comparison between the solvation by an ionic liquid and by aqueous alcohols

The thermo-solvatochrornic behaviors of 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate, RB; 2,6-dichloro-4-(2,4,6-triphenyloyridinium-1-yl) phenolate, WB; 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2); 2,6-dibromo-4-[(E)-2-(1-n-octylpyridinium-4-yl)ethenyl] phenolate, OcPMBr(2), have been investigated in binary mixtures of the ionic liquid, IL, 1-(1-butyl)-3-methylimidazolium tetrafluorborate, [BuMeIm][BF(4)], and water (W), in the temperature range from 10 to 60 degrees C. Plots of the empirical solvent polarities, ET (probe) in kcal mol(-1), versus the mole fraction of water in the binary mixture, chi(w) showed nonlinear, i.e., nonideal behavior. Solvation by these IL-W mixtures shows the following similarities to that by aqueous aliphatic alcohols: The same solvation model can be conveniently employed to treat the data obtained; it is based on the presence in the system-bulk medium and probe solvation shell of IL, W, and the ""complex"" solvent 1:1 IL-W. The origin of the nonideal solvation behavior appears to be the same, preferential solvation of the probe, in particular by the complex solvent. The strength of association of the IL-W complex, and the polarity of the IL are situated between the corresponding values of aqueous methanol and aqueous ethanol. Temperature increase causes a gradual desolvation of all probes employed. A difference between solvation by IL-W and aqueous alcohols is that probe-solvent hydrophobic interactions appear to play a minor role in case of the former mixture, probably because solvation is dominated by hydrogen-bonding and Coulombic interactions between the ions of the IL and the zwitterionic probes.

Selectivity in the mechanism of action of antimicrobial mastoparan peptide Polybia-MP1

Many potent antimicrobial peptides also present hemolytic activity, an undesired collateral effect for the therapeutic application. Unlike other mastoparan peptides, Polybia-MP1 (IDWKKLLDAAKQIL), obtained from the venom of the social wasp Polybia paulista, is highly selective of bacterial cells. The study of its mechanism of action demonstrated that it permeates vesicles at a greater rate of leakage on the anionic over the zwitterionic, impaired by the presence of cholesterol or cardiolipin; its lytic activity is characterized by a threshold peptide to lipid molar ratio that depends on the phospholipid composition of the vesicles. At these particular threshold concentrations, the apparent average pore number is distinctive between anionic and zwitterionic vesicles, suggesting that pores are similarly formed depending on the ionic character of the bilayer. To prospect the molecular reasons for the strengthened selectivity in Polybia-MP1 and its absence in Mastoparan-X, MD simulations were carried out. Both peptides presented amphipathic alpha-helical structures, as previously observed in Circular Dichroism spectra, with important differences in the extension and stability of the helix; their backbone solvation analysis also indicate a different profile, suggesting that the selectivity of Polybia-MP1 is a consequence of the distribution of the charged and polar residues along the peptide helix, and on how the solvent molecules orient themselves according to these electrostatic interactions. We suggest that the lack of hemolytic activity of Polybia-MP1 is due to the presence and position of Asp residues that enable the equilibrium of electrostatic interactions and favor the preference for the more hydrophilic environment.

Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Solvent effects on the electronic absorption spectrum of camphor using continuum, discrete or explicit approaches

We address the effect of solvation on the lowest electronic excitation energy of camphor. The solvents considered represent a large variation in-solvent polarity. We consider three conceptually different ways of accounting for the solvent using either an implicit, a discrete or an explicit solvation model. The solvatochromic shifts in polar solvents are found to be in good agreement with the experimental data for all three solvent models. However, both the implicit and discrete solvation models are less successful in predicting solvatochromic shifts for solvents of low polarity. The results presented suggest the importance of using explicit solvent molecules in the case of nonpolar solvents. (C) 2009 Elsevier B.V. All rights reserved.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Dipole polarizability and Rayleigh light scattering by the hydrated electron

Assuming the existence of a confined state of the electron in bulk water the polarizability of the hydrated electron is analyzed. Statistically uncorrelated supermolecular structures composed of seven water molecules (first solvation shell) with an extra electron were extracted from classical Monte Carlo simulation and used in quantum mechanical second-order Moller-Plesset calculations. It is found that the bound excess electron contributes with 274 a.u. to the total dipole polarizability of 345 a.u. for (H(2)O)(7)(-). From the calculated polarizabilities the Rayleigh elastic light scattering properties are inferred and found to considerably enhance activity and light depolarization. (C) 2009 Elsevier B.V. All rights reserved.

Experimental and theoretical investigation of first hyperpolarizability in aminophenols

This work reports theoretical and experimental studies on the first hyperpolarizability (beta) of aminophenols, evaluating the influence of the NH(2) group position relative to the OH group on the hyperpolarizability. A new extension of hyper-Rayleigh scattering technique using picosecond pulse trains was employed to obtain the experimental absolute values of (beta). The theoretical static beta(0) values were calculated using AMI method implemented in the AMPAC program. The theoretical and experimental data show a clear dependence between beta and the relative position of the electron donor (D) and acceptor (A) groups, presenting the 2-aminophenol the higher values. Moreover, calculations show excellent qualitative agreement between theoretical and experimental data, which are improved when the simulations considering the solvated molecule in a combination of discrete solvent molecules interacting with the solute and the application of continuous dielectric model. Besides, the study indicates that the experimental hyperpolarizabilities seem to be a property of the solute-solvation shell system. These facts have affirmed that the theoretical approach employed can be successfully used to foresee the variation in beta due to modifications in the D/A position. Moreover, a theoretical study of the ground state absorption is performed and compared with experimental data. (C) 2008 Elsevier B.V. All rights reserved.

Ionic Recognition by 7-Nitro-1,3,5-triaza Adamantane: First Thermodynamic Study

A thermodynamic study involving 7-nitro-1,3,5-triaza adamantane, 1, and its interaction with metal cations in nonaqueous media is first reported. Solubility data of 1 in various solvents were used to derive the standard Gibbs energies of solution, Delta G(s)degrees in these solvents. The effect of solvation in the different media was assessed from the Gibbs energy of transfer taking acetonitrile as a reference solvent. (1)H NMR studies of the interaction of 1 and metal cations were carried out in CD(3)CN and CD(3)OD and the data are reported. Conductance measurements revealed that this ligand forms lead(II) or zinc complexes of 1: 1 stoichiometry in acetonitrile. It also revealed a stoichiometry of two molecules of 1 per mercury(II) and two cadmiu (II) ions per molecule of 1. The addition of silver salt to 1 led to the precipitation of the silver-1 complex which was isolated and characterized by X-ray crystallography. At variance with conductance measurements in solution, in the solid state the X-ray structure show`s a 1:1 stoichiometry in the Hg(II) complex. The themiodynamics of complexation of 1 and these cations provide a quantitative assessment of the selective behavior of this ligand for ions of environmental relevance.

Combined Crystallographic and Solution Molecular Dynamics Study of Allosteric Effects in Ester and Ketone p-tert-Butylcalix[4]arene Derivatives and Their Complexes with Acetonitrile, Cd(II), and Pb(II)

We describe here a procedure to bridge the gap in the field of calixarene physicochemistry between solid-state atomic-resolution structural information and the liquid-state low-resolution thermodynamics and spectroscopic data. We use MD simulations to study the kinetics and energetics involved in the complexation of lower rim calix[4]arene derivatives (L), containing bidentate ester (1) and ketone (2) pendant groups, with acetonitrile molecule (MeCN) and Cd2+ and Pb2+ ions (M2+) in acetonitrile solution. On one hand, we found that the prior inclusion of MeCN into the calix to form a L(MeCN) adduct has only a weak effect in preorganizing the hydrophilic cavity toward metal ion binding. On the other hand, the strong ion-hydrophilic cavity interaction produces a wide open calix which enhances the binding of one MeCN molecule (allosteric effect) to stabilize the whole (M2+)1(MeCN) bifunctional complex. We reach two major conclusions: (i) the MD results for the (M2+)1(MeCN) binding are in close agreement with the ""endo"", fully encapsulated, metal complex found by X-ray diffraction and in vacuo MD calculations, and (ii) the MD structure for the more flexible 2 ligand, however, differs from the also endo solid-state molecule. In fact, it shows strong solvation effects at the calixarene lower bore by competing MeCN molecules that share the metal coordination sphere with the four C=O oxygens of an ""exo"" (M2+)2(MeCN) complex.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Solvatochromic Shift of the pi -> pi* Transition in All-Trans, Cis-13, Cis-11, Cis-9, and Cis-7 Retinal Isomers Induced by Water and Methanol

The solvatochromic shift of the lowest singlet it pi -> pi* electronic transition in the all-trans, cis-13, cis-11, cis-9, and cis-7 retinal isomers were computed under the influence of water, methanol, and benzene solvents. Excitation energies were calculated in gas phase and in solution. The calculations in solution were performed considering the sequential Monte Carlo (MC) /Quantum Mechanical approach. The MC simulations were performed considering the full retinal isomer molecules and 900 water molecules, 900 methanol, or 400 benzene ones. The OPLS/AA parametrization was chosen for retinal, methanol, and benzene molecules and the SPC model was used for water one. From the MC calculations 100 independent configurations were selected, with 100 solvent molecules in thermodynamical equilibrium at T = 298.15 K. Average point-charges were obtained from those independent configurations for water, methanol, and benzene solvent. TDDFT and CASSCF//CASPT2 methodologies were used to compute the vertical excitation energy of the retinal isomers in different environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2076-2087, 2010

Raman spectra of a pseudo-oxocarbon anion in ionic liquids

Raman and electronic spectra of the [3,5-bis(dicyanomethylene)cyclopentane-1,2,4-trionate] dianion, the croconate violet (CV), are reported in solutions of ionic liquids based on imidazolium cations. Different normal modes of the CV anion, nu (C=O), nu (CO) + nu (CC) + nu (CCN), and nu(C N), were used as probes of solvation characteristics of ionic liquids, and were compared with spectra of CV in common solvents. The spectra of CV in ionic liquids are similar to those in dichloromethane solution, but distinct from those in protic solvents such as ethanol or water. The UV-vis spectra of CV in ionic liquids strongly suggest pi-pi interactions between the CV anion and the imidazolium cation. Copyright (C) 2009 John Wiley & Sons, Ltd.