The energy landscape for solvent dynamics in electron transfer reactions: A minimalist model

Energy fluctuations of a solute molecule embedded in a polar solvent are investigated to depict the energy landscape for solvation dynamics. The system is modeled by a charged molecule surrounded by two layers of solvent dipolar molecules with simple rotational dynamics. Individual solvent molecules are treated as simple dipoles that can point toward or away from the central charge (Ising spins). Single-spin-flip Monte Carlo kinetics simulations are carried out in a two-dimensional lattice for different central charges, radii of outer shell, and temperatures. By analyzing the density of states as a function of energy and temperatures, we have determined the existence of multiple freezing transitions. Each of them can be associated with the freezing of a different layer of the solvent. (C) 2002 American Institute of Physics.

Smooth landscape solvent dynamics in electron transfer reactions

Solvent effects play a major role in controlling electron-transfer reactions. The solvent dynamics happens on a very high-dimensional surface, and this complex landscape is populated by a large number of minima. A critical problem is to understand the conditions under which the solvent dynamics can be represented by a single collective reaction coordinate. When this unidimensional representation is valid, one recovers the successful Marcus theory. In this study the approach used in a previous work [V. B. P. Leite and J. N. Onuchic; J. Phys. Chem. 100, 7680 (1996)] is extended to treat a more realistic solvent model, which includes energy correlation. The dynamics takes place in a smooth and well behaved landscape. The single shell of solvent molecules around a cavity is described by a two-dimensional system with periodic boundary conditions with nearest neighbor interaction. It is shown how the polarization-dependent effects can be inferred. The existence of phase transitions depends on a factor y proportional to the contribution from the two parameters of the model. For the present model, γ suggests the existence of weak kinetic phase transitions, which are used in the analysis of solvent effects in charge-transfer reactions. © 1999 American Institute of Physics.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Solvent induced NMR chemical shifts that arise from molecular encounters

Recently Homer and Percival have postulated that intermolecular van der Waals dispersion forces can be characterized by three mechanisms. The first arises via the mean square reaction field < R1; 2> due to the transient dipole of a particular solute molecule that is considered situated in a cavity surrounded by solvent molecules; this was characterized by an extended Onsager approach. The second stems from the extra cavity mean square reaction field < R2; 2> of the near neighbour solvent molecules. The third originates from square field electric fields E2BI due to a newly characterized effect in which solute atoms are `buffeted' by the peripheral atoms of adjacent solvent molecules. The present work concerns more detailed studies of the buffeting screening, which is governed by sterically controlled parameter (2T - T)2, where and are geometric structural parameters. The original approach is used to characterise the buffeting shifts induced by large solvent molecules and the approach is found to be inadequate. Consequently, improved methods of calculating and are reported. Using the improved approach it is shown that buffeting is dependent on the nature of the solvent as well as the nature of the solute molecule. Detailed investigation of the buffeting component of the van der Waals chemical shifts of selected solutes in a range of solvents containing either H or Cl as peripheral atoms have enabled the determination of a theoretical acceptable value for the classical screening coefficient B for protons. 1H and 13C resonance studies of tetraethylmethane and 1H, 13C and 29Si resonance studies of TMS have been used to support the original contention that three (< R1; 2> , < R2; 2> and E2BI) components of intermolecular van der Waals dispersion fields are required to characterise vdW chemical shifts.

The molecular design of a new solvent for the absorption of carbon dioxide

Gas absorption, the removal of one or more constitutents from a gas mixture, is widely used in chemical processes. In many gas absorption processes, the gas mixture is already at high pressure and in recent years organic solvents have been developed for the process of physical absorption at high pressure followed by low pressure regeneration of the solvent and recovery of the absorbed gases. Until now the discovery of new solvents has usually been by expensive and time consuming trial and error laboratory tests. This work describes a new approach, whereby a solvent is selected from considerations of its molecular structure by applying recently published methods of predicting gas solubility from the molecular groups which make up the solvent molecule. The removal of the acid gases of carbon dioxide and hydrogen sulfide from methane or hydrogen was used as a commercially important example. After a preliminary assessment to identify promising moecular groups, more than eighty new solvent molecules were designed and evaluated by predicting gas solubility. The other important physical properties were also predicted by appropriate theoretical procedures, and a commercially promising new solvent was chosen to have a high solubility for acid gases, a low solubility for methane and hydrogen, a low vapour pressure, and a low viscosity. The solvent chosen, of molecular structure Ch3-COCH2-CH2-CO-CH3, was tested in the laboratory and shown to have physical properties, except for vapour pressures, close to those predicted. That is gas solubilities were within 10% but lower than predicted. Viscosity within 10% but higher than predicted and a vapour pressure significantly lower than predicted. A computer program was written to predict gas solubility in the new solvent at the high pressures (25 bar) used in practice. This is based on the group contribution method of Skold Jorgensen (1984). Before using this with the new solvent, Acetonyl acetone, the method was show to be sufficiently accurate by comparing predicted values of gas solubility with experimental solubilities from the literature for 14 systems up to 50 bar. A test of the commercial potential of the new solvent was made by means of two design studies which compared the size of plant and approximate relative costs of absorbing acid gases by means of the new solvent with other commonly used solvents. These were refrigerated methanol(Rectisol process) and Dimethyl Ether or Polyethylene Glycol(Selexol process). Both studies showed in terms of capital and operating cost some significant advantage for plant designed for the new solvent process.

2-[(3,3-Dimethylindolin-2-ylidene)methyl]-4-[(3,3-dimethyl-3H-indol-1-ium-2-yl)methylidene]-3-oxocyclobut-1-en-1-olate chloroform disolvate

In the title squaraine dye solvate, C26H24N2O2·2CHCl3, the dye molecule is essentially planar, except for the methyl groups, having a maximum deviation over the 26-membered delocalized bond system of 0.060 (2) Å. It possesses crystallographic twofold rotational symmetry with the indole ring systems adopting a syn conformation. The molecular structure features intramolecular N-HO hydrogen bonds enclosing conjoint S7 ring motifs about one of the dioxocyclobutene O atoms, while the two chloroform solvent molecules are linked to the second O atom through C-HO hydrogen bonds.

Ultrafast Solvation Dynamics of an Ion in the gamma-Cyclodextrin Cavity: The Role of Restricted Environment

A detailed study of the solvation dynamics of a charged coumarin dye molecule in gamma-cyclodextrin/water has been carried out by using two different theoretical approaches. The first approach is based on a multishell continuum model (MSCM). This model predicts the time scales of the dynamics rather well, provided an accurate description of the frequency-dependent dielectric function is supplied. The reason for this rather surprising agreement is 2-fold. First, there is a cancellation of errors, second, the two-zone model mimics the heterogeneous microenvironment surrounding the ion rather well. The second approach is based on the molecular hydrodynamics theory (MI-IT). In this molecular approach, the solvation dynamics has been studied by restricting the translational motion of the solvent molecules enclosed within the cavity. The results from the molecular theory are also in good agreement with the experimental results. Our study indicates that, in the present case, the restricted environment affects only the long time decay of the solvation time correlation function. The short time dynamics is still governed by the librational (and/or vibrational) modes present in bulk water.

The Influence of Bilayer Composition on the Gel to Liquid Crystalline Transition

We report molecular dynamics simulations of bilayers using a united atom model with explicit solvent molecules. The bilayer consists of the single tail cationic surfactant behenyl trimethyl ammonium chloride (BTMAC) with stearyl alcohol (SA) as the cosurfactant. We study the gel to liquid crystalline transitions in the bilayer by varying the amount of water at fixed BTMAC to SA ratio as well as by varying the BTMAC to SA ratio at fixed water content. The bilayer is found to exist in the tilted, Lβ′ phase at low temperatures, and for the compositions investigated in this study, the Lβ′ to Lα melting transition occurred in the temperature range 330−338 K. For the highest BTMAC to SA composition (2:3 molar ratio), a diffuse headgroup−water interface is observed at lower temperatures, and an increase in the d-spacing occurs prior to the melting transition. This pretransition swelling is accompanied by a sharpening in the water density variation across the headgroup region of the bilayer. Signatures of this swelling effect which can be observed in the alkane density distributions, area per headgroup, and membrane thickness are attributed to the hydrophobic effect. At a fixed bilayer composition, the transition temperature (>338 K) from the Lβ′ to Lα transition obtained for the high water content bilayer (80 wt %) is similar to that obtained with low water content (54.3 wt %), confirming that the melting transition at these water contents is dominated by chain melting.

Self-assembly of three new coordination complexes: Formation of 2-D square grid, 1-D chain and tape structures

Three distinct coordination complexes, viz., [Co(imi)(2)(tmb)(2)] (1) [where imi = imidazole], {[Ni(tmb)(2)(H2O)(3)]center dot 2H(2)O}(n) (2) and [Cu-2(mu-tmb)(4)(CH3OH)(2)] (3), have been synthesized hydrothermally by the reactions of metal acetates,2,4,6-trimethylbenzoic acid (Htmb) and with or without appropriate amine. The Ni analogue of 1 and the Co analogue of 2 have also been synthesized. X-ray single-crystal diffraction suggests that complex 1 represents discrete mononuclear species and complex 2 represents a 1D chain coordination polymer in which the Ni(H) ions are connected by the bridging water molecules. Complex 3 represents a neutral dinuclear complex. In 1, the central metal ions are associated by the carboxylate moiety and imidazole ligands, whereas the central metal atom is coordinated to the carboxylate moiety and the respective solvent molecules in 2 and 3. In 3, the four 2,4,6-trimethylbenzoate moieties act as a bridge connecting two copper (11) ions and the 0 atoms of methanol coord geometry, with the methanol molecule at the apical position. In all the three structures the central metal atom sits on a crystallographic inversion centre. In all the cases, the coordination entities are further organized via hydrogen bonding interactions to generate multifarious supramolecular networks. Complexes 1, 2 and 3 have also been characterized by spectroscopic (UV/Vis and IR) and thermal analysis (TGA). In addition, the complexes were found to exhibit antimicrobial activity. The magnetic susceptibility measurements, measured from 8 to 300 K, revealed antiferromagnetic interactions between the Co(II) ions in compound 1 and the Ni(II) ions in la, respectively.

1-D hydrogen-bonded organization of hexanuclear {3d-4f-5d} complexes: evidence for slow relaxation of the magnetization for [{LMe2Ni(H2O)Ln(H2O)4.5}2{W(CN)8}2] with Ln = Tb and Dy

Heterometallic {3d-4f-5d} aggregates with formula [{LMe2Ni(H2O)Ln(H2O)4.5}2{W(CN)8}2]·15H2O, (LMe2 stands for N,N-2,2-dimethylpropylenedi(3-methoxysalicylideneiminato) Schiff-base ligand) with Ln = Gd, Tb, Dy, have been obtained by reacting bimetallic [LMe2Ni(H2O)2Ln(NO3)3] and Cs3{W(CN)8} in H2O. The hexanuclear complexes are organized in 1-D arrays by means of hydrogen bonds established between the solvent molecules coordinated to Ln and the CN ligands of an octacyanometallate moiety. The X-ray structure was solved for the Tb derivative. Magnetic behavior indicates ferromagnetic {W–Ni} and {Ni–Ln} interactions (JNiW = 18.5 cm-1, JNiGd = 1.85 cm-1) as well as ferromagnetic intermolecular interactions mediated by the H-bonds. Dynamic magnetic susceptibility studies reveal slow magnetic relaxation processes for the Tb and Dy derivatives, suggesting SMM type behavior for these compounds.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Vibrational phase relaxation of O–H stretch in bulk water: Role of large amplitude angular jumps and negative cross-correlations among the forces on the O–H bond

The ultrafast vibrational phase relaxation of O–H stretch in bulk water is investigated in molecular dynamics simulations. The dephasing time (T2) of the O–H stretch in bulk water calculated from the frequency fluctuation time correlation function (Cω(t)) is in the range of 70–80 femtosecond (fs), which is comparable to the characteristic timescale obtained from the vibrational echo peak shift measurements using infrared photon echo [W.P. de Boeij, M.S. Pshenichnikov, D.A. Wiersma, Ann. Rev. Phys. Chem. 49 (1998) 99]. The ultrafast decay of Cω(t) is found to be responsible for the ultrashort T2 in bulk water. Careful analysis reveals the following two interesting reasons for the ultrafast decay of Cω(t). (A) The large amplitude angular jumps of water molecules (within 30–40 fs time duration) provide a large scale contribution to the mean square vibrational frequency fluctuation and gives rise to the rapid spectral diffusion on 100 fs time scale. (B) The projected force, due to all the atoms of the solvent molecules on the oxygen (FO(t)) and hydrogen (FH(t)) atom of the O–H bond exhibit a large negative cross-correlation (NCC). We further find that this NCC is partly responsible for a weak, non-Arrhenius temperature dependence of the dephasing rate.

Aggregation studies in crystals of apolar helical peptides: Boc-Aib-Val-Ala-Leu-Aib-Val-Ala-Leu-Aib-OMe

In the crystal, the backbone of Boc-(Aib-Val-Ala-Leu)2-Aib-OMe adopts a helical form with four alpha-type hydrogen bonds in the middle, flanked by 3(10)-type hydrogen bonds at either end. The helical molecules stack in columns with head-to-tail hydrogen bonds, either directly between NH and CO, or bridged by solvent molecules. The packing of the helices is parallel, even in space group P2(1). Cell parameters are a = 9.837(2) A, b = 15.565(3) A, c = 20.087(5) A, beta = 96.42(2) degrees, dcalc = 1.091 g/cm3 for C46H83N9O12.1.5H2O.0.67CH3OH. There appears to be some hydration of the backbone in this apolar helix.

Microcalorimetric indications for ligand binding as a function of the protein for galactoside-specific plant and avian lectins

The process cascade leading to the final accommodation of the carbohydrate ligand in the lectin’s binding site comprises enthalpic and entropic contributions of the binding partners and solvent molecules. With emphasis on lactose, N-acetyllactosamine, and thiodigalactoside as potent inhibitors of binding of galactoside-specific lectins, the question was addressed to what extent these parameters are affected as a function of the protein. The microcalorimetric study of carbohydrate association to the galectin from chicken liver (CG-16) and the agglutinin from Viscum album (VAA) revealed enthalpy–entropy compensation with evident protein type-dependent changes for N-acetyllactosamine. Reduction of the entropic penalty by differential flexibility of loops or side chains and/or solvation properties of the protein will have to be reckoned with to assign a molecular cause to protein type-dependent changes in thermodynamic parameters for lectins sharing the same monosaccharide specificity.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.