Dynamics of 4-benzylamino-7-nitrobenzofurazan in the 1-propanol/water binary solvent system. Evidence for composition-dependent solvent organization

We report on the fluorescence lifetime and rotational diffusion dynamics of 4-benzylamino-7-nitrobenzofurazan (BBD) in a series of 1-propanol/water binary solvent systems. The fluorescence lifetime of BBD increases monotonically with increasing 1-propanol concentration. The rotational diffusion dynamics of BBD also vary with solution 1-propanol content, but this variation is not monotonic. Comparison of the BBD rotational diffusion time constant to solution viscosity and 1-propanol composition reveals the presence of a solution composition dependence of solvent-solute interactions, with a relative decrease in solvent-solute interaction strength for solvent system compositions where the 1-propanol/water azeotrope is known to exist. These data point collectively to the existence of microscopic heterogeneity in these binary solvent systems.     

Isomerization reaction dynamics and equilibrium at the liquid-vapor interface of water. A molecular-dynamics study

The gauche-trans isomerization reaction of 1,2-dichloroethane at the liquid-vapor interface of water is studied using molecular-dynamics computer simulations. The solvent bulk and surface effects on the torsional potential of mean force and on barrier recrossing dynamics are computed. The isomerization reaction involves a large change in the electric dipole moment, and as a result the trans/gauche ratio is considerably affected by the transition from the bulk solvent to the surface. Reactive flux correlation function calculations of the reaction rate reveal that deviation from the transition-state theory due to barrier recrossing is greater at the surface than in the bulk water. This suggests that the system exhibits non-Rice-Ramsperger-Kassel-Marcus behavior due to the weak solvent-solute coupling at the water liquid-vapor interface.     

Photodissociation of ICN at the liquid/vapor interface of water

The photodissociation of ICN adsorbed at the liquid/vapor interface of water is studied using classical molecular dynamics with nonadiabatic surface hopping. The cage escape, geminate recombination to form ICN and INC and the subsequent vibrational relaxation of these two molecules (on their ground electronic states) is compared with the same process in bulk water and with previous photodissociation studies at liquid interfaces. We find that the reduced surface density and weaker solvent-solute interactions give rise to reduced rate of nonadiabatic transitions and that the probability for cage escape at the interface is significantly enhanced due to the possibility that one or both of the photodissociation fragments desorb into the gas phase. The overall desorption probability varies from 75% to 92% for ICN initially located just below the Gibbs surface (50% bulk density) to ICN located just above the Gibbs surface, respectively. The corresponding geminate recombination probabilities are 18% and 9%, respectively. The vibrational relaxation rate of the recombined ICN is slower than in the bulk by a factor of 2.3.     

Including anharmonicity in the calculation of rate constants. 1. The HCN/HNC isomerization reaction

A method for calculating anharmonic vibrational energy levels in asymmetric top and linear systems that is based on second-order perturbation theory in curvilinear coordinates is extended to the bound generalized normal modes at nonstationary points along a reaction path. Explicit formulas for the anharmonicity coefficients, x(ij), and the constant term, E0, are presented, and the necessary modifications for resonance cases are considered. The method is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the HCN/HNC isomerization reaction. Although the results for this system are not very sensitive to the choice of coordinates, we find that the inclusion of anharmonicity leads to a substantial improvement in the vibrational energy levels. We also present detailed comparisons of rate constants computed with and without anharmonicity, with various approximations for incorporating tunneling along the reaction path, and with a more practical approach to calculating the vibrational partition functions needed for larger systems.     

The dissociation rate of water molecules in systems with cation- and anion-exchange membranes

Within the framework of the mathematical model of Nernst-Planck-Poisson, an attempt is undertaken to theoretically describe the electrodiffusion of ions in the system diffusion layer/monopolar ionexchange membrane, which is accompanied by dissociation of water molecules. The formulas for estimating the current density transferred through a monopolar membrane by hydrogen or hydroxyl ions formed in dissociation of water in the space-charge region are derived. The rate constants and other parameters of dissociation of water molecules in the space-charge region of monopolar membranes under conditions of stabilization of the diffusion layer thickness are calculated. Their comparative analysis with the similar characteristics of bipolar membranes is carried out. For the phosphoric-acid heterogeneous membrane MK-41 in which the polarization conditions in the current density range under study are not so severe and the reaction layer is not being depleted as in the bipolar membrane MB-3 (contains the same phosphoric-acid groups), it is shown that only single-charged phosphoric-acid groups are involved in the water dissociation reaction. For MK-41, the calculated constants of the heterolytic reaction of water molecule dissociation are lower than for the heterogeneous membrane MA-40 containing ternary and quaternary amino groups. It is confirmed that the nature of ionogenic groups in membranes is a factor that determines the rate of water dissociation in systems with ion-exchange membranes.     

Polysaccharide charge density regulating protein adsorption to air/water interfaces by protein/polysaccharide complex formation

Because the formation of protein/polysaccharide complexes is dominated by electrostatic interaction, polysaccharide charge density is expected to play a major role in the adsorption behavior of the complexes. In this study, pullulan (a non-charged polysaccharide) carboxylated to four different charge densities (fraction of carboxylated subunits: 0.1, 0.26, 0.51, and 0.56) was used to investigate the effect of charge density on the properties of mixed protein/polysaccharide adsorbed layers at air/water interfaces. With all pullulan samples, soluble complexes with beta-lactoglobulin could be formed at low ionic strength, pH 4.5. It was shown that the higher was the pullulan charge density, the more the increase of surface pressure in time was retarded as compared to that for pure beta-lactoglobulin. The retardation was even more pronounced for the development of the dilatational modulus. The lower dilatational modulus can be explained by the ability of the polysaccharides to prevent the formation of a compact protein layer at the air/water interface due to electrostatic repulsion. This ability of the polysaccharides to prevent "layer compactness" increases with the net negative charge of the complexes. If charge density is sufficient (> or = 0.26), polysaccharides may enhance the cohesion between complexes within the adsorbed layer. The charge density of polysaccharides is shown to be a dominant regulator of both the adsorption kinetics as well as the resulting surface rheological behavior of the mixed layers formed. These findings have significant value for the application of complex protein-polysaccharide systems.     

Quantitative determination of alpha-tocopherol distribution in a tributyrin/Brij 30/water model food emulsion

Until recently, determining the distribution of antioxidants, AOs, between the oil, interfacial and aqueous regions of opaque emulsions has not worked well because the concentrations of AOs in interfacial regions cannot be determined separately from their concentrations in the oil and water phases. However, our novel kinetic method based on the reaction between an arenediazonium ion and vitamin E, or alpha-tocopherol, provides the first good estimates for the two partition constants that describe alpha-tocopherol distribution between the oil/interfacial and water/interfacial regions of tributyrin/Brij 30/water emulsions without physical isolation of any phase. The reaction is monitored by a new derivatization method based on trapping unreacted arenediazonium ion as an azo dye and confirmed by linear sweep voltammetry, LSV. The results by both derivatization and LSV methods are in good agreement and show that alpha-tocopherol distributes strongly in favor of the interfacial region when the oil is tributyrin, e.g., ca. 90% when the surfactant volume fraction is Phi I=0.01. The second-order rate constant for reaction in the interfacial region is also obtained from the results. Our kinetic method provides a robust approach for determining antioxidant distributions in emulsions and should help develop a quantitative interpretation of antioxidant efficiency in emulsions.     

Interactions between phosphate and water in solution: a natural bond orbital based analysis in a QM/MM framework

The natural bond orbital (NBO) and natural energy decomposition analysis (NEDA) calculations are used to analyze the interaction between mono-methyl phosphate-ester (MMP) and its solvation environment in a combined quantum mechanical/molecular mechanical (QM/MM) framework. The solute-solvent configurations are generated using a specific parametrization of the self-consistent-charge density functional tight-binding (SCC-DFTB) model for the MMP and TIP3P for water. The NBO and NEDA calculations are done with several QM/MM partitioning schemes with HF/6-31+G** as the QM level. Regardless of the size of the QM region, a notable amount of charge transfer is observed between MMP and the neighboring water molecules and the charge-transfer interactions are, in the NEDA framework, as important as the electric (electrostatic and polarization) components. This work illustrates that NBO based analyses are effective tools for probing intermolecular interactions in condensed phase systems.     

Magnitude of CH/O interactions between carbohydrate and water

The interaction energy potentials for six orientations of the fucose–water complex were calculated for evaluating the magnitude of the CH/O interactions in the complex. The calculations show that the C–H bonds of the nonpolar surface of fucose prefer to have contact with the oxygen atom of the water. The substantial attraction exists between the C–H bonds of fucose and water. The interaction energy calculated for the fucose–water complex at the MP2/aug-cc-pVTZ level is −2.55 kcal/mol. The CH/O interactions in the fucose–water complex are significantly larger than those in the cyclohexane–water complex (−1.13 kcal/mol), which shows that the oxygen atoms of fucose enhance the CH/O interactions. The electrostatic and dispersion interactions are responsible for the attraction in the CH/O interactions in the fucose–water complex, while the electrostatic contributions to the attraction in the CH/O interactions in the cyclohexane–water complex is small. The DFT-SAPT calculations also show that the electrostatic interactions are responsible for the larger attraction in the fucose–water complex. These results suggest that the nature of the CH/O interactions between carbohydrate and water is significantly different from that of the CH/O interactions between saturated hydrocarbon and water.     

A theoretical thermochemical study of solute-solvent dielectric effects in the displacement of codon-anticodon base pairs

Quantum-chemical solvent effect theories describe the electronic structure of a molecular subsystem embedded in a solvent or other molecular environment. The solvation of biomolecules is important in molecular biology, since numerous processes involve proteins interacting in changing solvent-solute systems. In this theoretical study, we focus on mRNA-tRNA base pairs as a fundamental step in protein synthesis influenced by hydrogen bonding between two antiparallel trinucleotides, namely, the mRNA codon and tRNA anticodon. We use the mean reaction field theories, which describe electrostatic and polarization interactions between solute and solvent in the AAA, UUU, AAG, and UUC triplex sequences optimized in various solvent media such as water, dimethylsulfoxide, methanol, ethanol, and cyclopean using the self-consistent reaction field model. This process depends on either the reaction potential function of the solvent or charge transfer operators that appear in solute-solvent interaction. Because of codon and anticodon biological criteria, we performed nonempirical quantum-mechanical calculations at the BLYP and B3LYP/3-21G, 6-31G, and 6-31G* levels of theory in the gas phase and five solvents at three temperatures. Finally, to obtain more information, we calculated thermochemical parameters to find that the dielectric constant of solvents plays an important role in the displacement of amino acid sequences on codon-anticodon residues in proteins, which can cause some mutations in humans.     

Preferential solvation of some antiepileptic drugs in {cosolvent (1) + water (2)} mixtures at 298.15 K

The solubility of lamotrigine (LTG), clonazepam (CZP) and diazepam (DZP) in some {cosolvent (1) + water (2)} mixtures expressed in mole fraction at 298.15 K was calculated from reported solubility values expressed in molarity by using the densities of the saturated solutions. Aqueous binary mixtures of ethanol, propylene glycol and N-methyl-2-pyrrolidone were considered. From mole fraction solubilities and some thermodynamic properties of the solvent mixtures, the preferential solvation of these drugs by both solvents in the mixtures was analysed by using the inverse Kirkwood–Buff integrals. It is observed that LTG, CZP and DZP are preferentially solvated by water in water-rich mixtures in all the three binary systems analysed. In {ethanol (1) + water (2)} mixtures, preferential solvation by water is also observed in ethanol-rich mixtures. Nevertheless, in {propylene glycol (1) + water (2)} and {N-methyl-2-pyrrolidone (1) + water (2)} mixtures preferential solvation by the cosolvent was observed in cosolvent-rich mixtures.     

Dynamics of an enzymatic substitution reaction in haloalkane dehalogenase

Reactive flux molecular dynamics simulations have been carried out using a combined QM/MM potential to study the dynamics of the nucleophilic substitution reaction of dichloroethane by a carboxylate group in haloalkane dehalogenase and in water. We found that protein dynamics accelerates the reaction rate by a factor of 2 over the uncatalyzed reaction. Compared to the thermodynamic effect in barrier reduction, protein dynamic contribution is relatively small. However, analyses of the friction kernel reveal that the origins of the reaction dynamics in water and in the enzyme are different. In aqueous solution, there is significant electrostatic solvation effect, which is reflected by the slow reorganization relaxation of the solvent. On the other hand, there is no strong electrostatic coupling in the enzyme and the major effect on reaction coordinate motion is intramolecular energy relaxation.     

Synthesis of thiiranes from oxiranes in water using polymeric cosolvents

Poly(vinylamine) (PVA) and poly(allylamine) (PAA) were prepared and used as polymeric cosolvents. Oxiranes were converted efficiently to the corresponding thiiranes under mild reaction conditions in water with ammonium thiocyanate (NH₄SCN) using these polymeric cosolvents. The polymeric cosolvents are reusable.     

Enthalpy-entropy compensation in micellization of sodium dodecyl sulphate in water/methanol,water/ethylene glycol and water/glycerol binary mixtures

The enthalpy-entropy compensation in micellization of sodium dodecyl sulphate (SDS) in binary mixtures of water/methanol (MeOH), water/ethylene glycol (EG) and water/glycerol (GL) over a temperature range of 10–60°C was examined. When the cosolvent concentration was low, the critical micelle concentration (CMC) depended only on the total amount of the hydroxyl group added. When the cosolvent concentration was high, the increase in CMC followed the sequence: MeOH>EG>GL. Enthalpy and entropy changes were evaluated from which the compensation temperature was determined. Both enthalpy and entropy changes decreased on the addition of the cosolvents, indicating a lowering of solution hydrophobicity. The compensation temperature was found as a constant over the cosolvent concentration range, as a result, was not a good index for characterizing the solute/solvent interactions. The two reference temperatures at which the enthalpy-entropy change respectively became zero were strongly influenced by the cosolvent addition, therefore could serve as a proper index for solution hydrophobicity.     

Ion-selective Marangoni instability coupled with the nonlinear adsorption/desorption rate

An oil/water interface containing bis(2-ethylhexyl)phosphate and Ca(2+) or Fe(3+) exhibits spontaneous Marangoni instability associated with the fluctuation in interfacial tension. This instability rarely appears for oil/water systems with Mg(2+), Sr(2+), Ba(2+), Cu(2+), or Co(2+). The same ion selectivity is observed for n-heptane and nitrobenzene despite their significant differences in density, viscosity, and the dielectric constant of oil. We studied this instability under acidic pH conditions to avoid the neutralization reaction effects. The result of the equilibrium interfacial tension and the extraction ratio of cations indicates that a large number of oil-soluble complexes form at the interfaces of Ca(2+)-containing systems and probably for Fe(3+)-containing systems. The results obtained by oscillating drop tensiometry and Brewster angle microscopy indicate that desorption, rather than adsorption, is more significant to the onset of instability and that the resulting complex tends to form aggregates in the interface. This aggregation gives the nonlinear desorption rate of the oil-soluble complex. Then, exfoliation of the aggregating matter occurs, which triggers the Marangoni instability. The induced convection removes the oil-soluble complex accumulated at the interface, creating a renewed interface, which is necessary for the successive occurrence of the Marangoni instability. For the other cations, the oil-soluble compounds are insignificant, and they rarely form aggregates. In such cases, adsorption/desorption proceeds without instability.     

Protein-lipid interactions at the air/water interface

Surface pressure measurements and external reflection FTIR spectroscopy have been used to probe protein-lipid interactions at the air/water interface. Spread monomolecular layers of stearic acid and phosphocholine were prepared and held at different compressed phase states prior to the introduction of protein to the buffered subphase. Contrasting interfacial behaviour of the proteins, albumin and lysozyme, was observed and revealed the role of both electrostatic and hydrophobic interactions in protein adsorption. The rate of adsorption of lysozyme to the air/water interface increased dramatically in the presence of stearic acid, due to strong electrostatic interactions between the negatively charged stearic acid head group and lysozyme, whose net charge at pH 7 is positive. Introduction of albumin to the subphase resulted in solubilisation of the stearic acid via the formation of an albumin-stearic acid complex and subsequent adsorption of albumin. This observation held for both human and bovine serum albumin. Protein adsorption to a PC layer held at low surface pressure revealed adsorption rates similar to adsorption to the bare air/water interface and suggested very little interaction between the protein and the lipid. For PC layers in their compressed phase state some adsorption of protein occurred after long adsorption times. Structural changes of both lysozyme and albumin were observed during adsorption, but these were dramatically reduced in the presence of a lipid layer compared to that of adsorption to the pure air/water interface.     

Computational investigation of order, structure, and dynamics in modified water models

Model liquids have been constructed to study the role of local structure in the anomalous properties of liquid water. The intermolecular potentials were modified by increasing the weight of the Lennard-Jones term relative to the electrostatic term in the SPC/E model for water. The resulting family of liquids varies from SPC/E water to a Lennard-Jones-like liquid. Properties were measured as a function of density and temperature. The local structure was described by two order parameters, one measuring the tetrahedral order and the other measuring the translational order. The translational order parameter was found to be large for both tetrahedral and Lennard-Jones liquids, but to go through a minimum as the potentials were modified, demonstrating that the two types of structure are incompatible. Just as in water several properties (e.g., the translational diffusion coefficient, entropy) exhibit anomalous density dependence as a result of the breakdown of local tetrahedrality, we observed nonmonotonic behavior of the translational diffusion constant and reorientational relaxation rate as the fluids were transformed from tetrahedral to Lennard-Jones-like. This is also an indication of the incompatibility between Lennard-Jones and water-like structure.     

Self-assembly and influence of the organic counterion in the ternary systems dodecylamine/acrylic acid/water and dodecylamine/methacrylic acid/water

Due to complete proton transfer from the acid to the amine, a reaction between an equimolar mixture of dodecylamine and (meth)acrylic acid leads to the formation of dodecylammonium (meth)acrylate. The latter can be considered as a surfactant with a polymerizable organic counterion. The ternary phase diagrams of the two systems dodecylamine/acrylic acid/water and dodecylamine/methacrylic acid/water are described. Both systems can form isotropic solutions and lyotropic liquid crystalline lamellar phases. Moreover, the system with the methacrylate counterion can also form a cubic phase in the water-rich part of the phase diagram. The difference in the self-organization observed for the two systems is explained by the greater bulkiness and hydrophobicity of the methacrylate. Whereas the acrylate counterion behaves rather like a classic inorganic counterion, the methacrylate counterion resides in the outermost part of the aggregates, giving rise to a change in the surface curvature.     

Prediction of the mechanisms of enzyme-catalysed reactions using hybrid quantum mechanical/molecular mechanical methods

The use of hybrid methods, involving both quantum mechanics and molecular mechanics, to model the mechanism of enzyme-catalysed reactions, is discussed. Two alternative approaches to treating the electrostatic interactions between the quantum mechanical and molecular mechanical regions are studied, involving either the inclusion of this term in the electronic Hamiltonian (QM/MM), or evaluating it purely classically (MO + MM). In the latter scheme, possible problems of using force fields that are standard for macromolecular modelling are identified. The use of QM/MM schemes to investigate the mechanism of the enzymes thymidine phosphorylase (ThdPase) and protein tyrosine phosphatase (PTP) is described. For both systems, transition states have been identified using a PM3 Hamiltonian. For ThdPase, concerted motion of the enzyme during the course of the reaction is suggested and, for PTP, a two-step dephosphorylation reaction is indicated, both with quite low barriers.     

Polymer/surfactant interactions at the air/water interface

The development of neutron reflectometry has transformed the study and understanding of polymer/surfactant mixtures at the air/water interface. A critical assessment of the results from this technique is made by comparing them with the information available from other techniques used to investigate adsorption at this interface. In the last few years, detailed information about the structure and composition of adsorbed layers has been obtained for a wide range of polymer/surfactant mixtures, including neutral polymers and synthetic and naturally occurring polyelectrolytes, with single surfactants or mixtures of surfactants. The use of neutron reflectometry together with surface tensiometry, has allowed the surface behaviour of these mixtures to be related directly to the bulk phase behaviour. We review the broad range of systems that have been studied, from neutral polymers with ionic surfactants to oppositely charged polyelectrolyte/ionic surfactant mixtures. A particular emphasis is placed upon the rich pattern of adsorption behaviour that is seen in oppositely charged polyelectrolyte/surfactant mixtures, much of which had not been reported previously. The strong surface interactions resulting from the electrostatic attractions in these systems have a very pronounced effect on both the surface tension behaviour and on adsorbed layers consisting of polymer/surfactant complexes, often giving rise to significant surface ordering.