Integration of experiment, theory and modeling to understand the interaction type and kinetics of limonene on silica surfaces.

The high membrane permeability of H2S was studied using polarizable molecular dynamics simulations of a DPPC lipid bilayer. The solubility-diffusion model predicts permeability coefficients of H2S and H2O that are in good agreement with experiment. The computed diffusion coefficient profile shows H2S to diffuse at a lower rate than H2O, but the barrier for H2S permeation on the Gibbs energy profile is negligible. The hydrophobicity of H2S allows it to partition into the paraffinic interior of the membrane readily.

A polarizable force field for liquid hydrogen sulfide (H2S) has been developed based on the Drude oscillator model. This force field has been designed to be analogous to the SWM4-NDP water model; the model is rigid with point charges assigned to the H and S atoms and a lone pair on the bisector of ∠HSH in the molecular plane. Positions of the lone pair and the charges have been defined such that the model has a static dipole moment of 0.98 D, equal to the experimental value. Polarizability is incorporated by a charged (Drude) particle attached to the S atom through a harmonic potential. Intermolecular nonbonded forces are included by use of a Lennard-Jones potential between S atoms. The model was parametrized to reproduce the density, enthalpy of vaporization, and dielectric constant of pure H2S at 212 K and 1 atm. The calculated density, enthalpy of vaporization, shear viscosity coefficient, and self-diffusion coefficient are in good agreement with experiment over the temperature range 212-298 K along the liquid-vapor coexistence curve of liquid H2S. The radial distribution function calculated from this model is in good agreement with experimental diffraction data and ab initio molecular dynamics simulations.

The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM-TURBOMOLE is well-suited for extended QM/MM MD simulations using first principles methods with large (triple-ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg(2+) (aq), where the MM outer sphere water molecules are represented using the SWM4-NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg(2+) and Zn(2+) were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time-averaged solution absorption spectrum of indole, the emission spectrum of the indole 1La excited state, and the electronic circular dichroism spectrum of an oxacepham.

Molecular dynamics (MD) simulations using the Drude polarizable force field are used to study the solution and interfacial properties of hydrogen sulfide (H2S) in water. Pairwise H2O-H2S Lennard-Jones interactions were optimized to the experimental H2S gas solubility at 298 K. These parameters yield hydration free energies and diffusion coefficients for H2S that are in good agreement with the experiment over 273-323 K and 298-368 K, respectively. H2S is sparingly soluble in water, with a ΔG(hydr)° of -0.5 kcal mol(-1). The free energy perturbation (FEP) calculations and analysis of the radial distribution functions show that H2S has limited hydrogen bonding and electrostatic interactions with the water solvent and generally behaves like a hydrophobic solute. These features were confirmed by ab initio MD simulations. Umbrella sampling simulations were used to calculate the free energy profile of the transition of H2S across the water-vapor interface, which showed that H2S has a sizable surface excess, with a ΔG(surf) of 1.3 kcal mol(-1). This high surface excess is consistent with our calculations of the surface tension, which decreases to 20 dyn cm(-1) under high densities of H2S (g). The dipole moment of H2S increases from its gas phase value of 0.98 to 1.25 D in bulk water as it moves across the interface. Adsorbed H2S tends to be oriented perpendicular to the interface, with the sulfur atom pointing toward the vapor phase.

The indoor environment is a dynamic one with many variables impacting indoor air quality and indoor air chemistry. These include relative humidity (RH) and the presence of different surfaces. Although it has been suggested that the indoor concentrations of gas-phase compounds increase at higher relative humidity, because of displacement of these compounds from indoor surfaces, little is known from a molecular perspective about how RH and adsorbed water impact the adsorption of indoor relevant organic compounds such as limonene with indoor relevant surfaces. Herein, we investigate the effects of RH on the adsorption of limonene, a hydrophobic molecule, on hydroxylated SiO 2 surfaces, a model for glass surfaces. Experimental data using infrared spectroscopy to directly measure limonene adsorption are combined with both force field-based molecular dynamics (MD) and ab initio molecular dynamics (AIMD) simulations to understand the competitive interactions between limonene, water, and the SiO 2 surface. The spectroscopic data provide evidence that adsorbed limonene is not completely displaced by adsorbed water, even at high RH (∼80%) when the water layer coverage is close to three monolayers (MLs). These experimental data are supported by AIMD and MD simulations, which indicate that limonene is present at the adsorbed water interface but displaced from direct interactions with SiO 2 . This study shows that although some limonene can desorb from the surface, even at the highest RH, more than half the limonene remains adsorbed on the surface that can undergo continued surface reactivity. A complex network of π-hydrogen bonds, water−water hydrogen bonds, and SiO 2 −water hydrogen bonds explains these interactions at the air/adsorbed water/SiO 2 interface that hold the hydrophobic limonene molecule at the interface. Importantly, these interactions are most likely present for a range of other systems involving organic compounds and solid surfaces at ambient relative humidity and may be important in a range of scientific areas, from sensor development to cultural heritage science.

The adsorption of limonene, a common organic compound found in indoor air, on hydrophilic surfaces such as glass (SiO 2 ), a prevalent surface in the indoor environment, is poorly understood. In this study, we have investigated the interaction of limonene and three other cyclic hydrocarbons (cyclohexane, cyclohexene, and benzene) on hydroxylated SiO 2 using infrared spectroscopy and ab initio molecular dynamics (AIMD) simulations. Experimental results show that there is an interaction between these cyclic hydrocarbons and surface hydroxyl groups. AIMD simulations demonstrate that all of the cyclic molecules, except for cyclohexane, π-hydrogen bond with surface hydroxyl groups while cyclohexane interacts with the surface OH groups through dispersion forces. According to experiments and simulations, the intermolecular interaction between limonene and SiO 2 is significantly stronger than those of other compounds explored. This study provides an understanding of some of the driving forces behind the formation of organic coatings on glass surfaces important in indoor environments.

The human voltage-gated proton channel Hv1 is a drug target for cancer, ischemic stroke, and neuroinflammation. It resides on the plasma membrane and endocytic compartments of a variety of cell types, where it mediates outward proton movement and regulates the activity of NOX enzymes. Its voltage-sensing domain (VSD) contains a gated and proton-selective conduction pathway, which can be blocked by aromatic guanidine derivatives such as 2-guanidinobenzimidazole (2GBI). Mutation of Hv1 residue F150 to alanine (F150A) was previously found to increase 2GBI apparent binding affinity more than two orders of magnitude. Here, we explore the contribution of aromatic interactions between the inhibitor and the channel in the presence and absence of the F150A mutation, using a combination of electrophysiological recordings, classic mutagenesis, and site-specific incorporation of fluorinated phenylalanines via nonsense suppression methodology. Our data suggest that the increase in apparent binding affinity is due to a rearrangement of the binding site allowed by the smaller residue at position 150. We used this information to design new arginine mimics with improved affinity for the nonrearranged binding site of the wild-type channel. The new compounds, named “Hv1 Inhibitor Flexibles” (HIFs), consist of two “prongs,” an aminoimidazole ring, and an aromatic group connected by extended flexible linkers. Some HIF compounds display inhibitory properties that are superior to those of 2GBI, thus providing a promising scaffold for further development of high-affinity Hv1 inhibitors.