Hybrid quantum and classical mechanical Monte Carlo simulations of the interaction of hydrogen chloride with solid water clusters

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

Estimating the rate constant of cyclic GMP hydrolysis by activated phosphodiesterase in photoreceptors

The early steps of light response occur in the outer segment of rod and cone photoreceptor. They involve the hydrolysis of cGMP, a soluble cyclic nucleotide, that gates ionic channels located in the outer segment membrane. We shall study here the rate by which cGMP is hydrolyzed by activated phosphodiesterase (PDE). This process has been characterized experimentally by two different rate constants beta(d) and beta(sub): beta(d) accounts for the effect of all spontaneously active PDE in the outer segment, and beta(sub) characterizes cGMP hydrolysis induced by a single light-activated PDE. So far, no attempt has been made to derive the experimental values of beta(d) and beta(sub) from a theoretical model, which is the goal of this work. Using a model of diffusion in the confined rod geometry, we derive analytical expressions for beta(d) and beta(sub) by calculating the flux of cGMP molecules to an activated PDE site. We obtain the dependency of these rate constants as a function of the outer segment geometry, the PDE activation and deactivation rates and the aqueous cGMP diffusion constant. Our formulas show good agreement with experimental measurements. Finally, we use our derivation to model the time course of the cGMP concentration in a transversally well-stirred outer segment.

Loss of chlorine in mass spectra of DCl picked up by water clusters in a beam

Water clusters (H(2)O)(n), n<or=16, are produced by supersonic expansion of water vapor into vacuum, and then pick up a DCl molecule. The resulting mixed clusters are analyzed by electron bombardment ionization mass spectrometry. In all cases observed, the chlorine atom is lost in the ionization process, producing a deuterated water cluster cation [(H(2)O)(n)D(+)]. This suggests that the chlorine atom stays on the surface and has a weaker bond to the host cluster.

Uptake of Phenol on Aerosol Particles

We present a study of the interaction between a phenol molecule and an aerosol particle. The aerosol particle is represented by a cluster of 128 water molecules. Using a classical approach, we present interaction energy surfaces for different relative distances and for three orientations of phenol relative to the particle. From the energy surfaces we find the reaction pathways with the largest interaction between the molecule and the particle. We use a quantum mechanics/molecular mechanics (QM/MM) method to calculate a potential energy curve for each reaction path. Coupled cluster methods are used for the part of the system described by quantum mechanics, while the part described by molecular mechanics is represented by a polarizable force field. We compare results obtained from the classical approach with the QM/MM results. Furthermore, we use the QM/MM results to calculate mass accommodation coefficients using a quantum-statistical (QM-ST) model and show how the mass accommodation coefficient depends on the relative orientation of phenol with respect to the aerosol particle.

Reactions on cell membranes: comparison of continuum theory and Brownian dynamics simulations

Biochemical transduction of signals received by living cells typically involves molecular interactions and enzyme-mediated reactions at the cell membrane, a problem that is analogous to reacting species on a catalyst surface or interface. We have developed an efficient Brownian dynamics algorithm that is especially suited for such systems and have compared the simulation results with various continuum theories through prediction of effective enzymatic rate constant values. We specifically consider reaction versus diffusion limitation, the effect of increasing enzyme density, and the spontaneous membrane association/dissociation of enzyme molecules. In all cases, we find the theory and simulations to be in quantitative agreement. This algorithm may be readily adapted for the stochastic simulation of more complex cell signaling systems.

Adsorption of HCl on the Water Ice Surface Studied by X-ray Absorption Spectroscopy

The adsorption state of HCl at 20 and 90 K on crystalline water ice films deposited under ultrahigh vacuum at 150 K has been studied by X-ray absorption spectroscopy at the O1s K-edge and Cl2p L-edge. We show that HCl dissociates at temperatures as low as 20 K, in agreement with the prediction of a spontaneous ionization of HCl on ice. Comparison between the rate of saturation of the "dangling" hydrogen bonds and the chlorine uptake indicates that hydrogen bonding of HCl with the surface native water "dangling" groups only accounts for a small part of the ionization events (20% at 90 K). A further mechanism drives the rest of the dissociation/solvation process. We suggest that the weakening of the ice surface hydrogen-bond network after the initial HCl adsorption phase facilitates the generation of new dissociation/solvation sites, which increases the uptake capacity of ice. These results also emphasize the necessity to take into account not only a single dissociation event but its catalyzing effect on the subsequent events when modeling the uptake of hydrogen-bonding molecules on the ice surface.