Melting of aqueous NaCl solutions in porous materials: shifted phase transition distribution (SIDI) approach for determining NMR cryoporometry pore size distributions.

Nuclear magnetic resonance cryoporometry (NMRC) and differential scanning calorimetry thermoporometry (DSC-TPM) are powerful methods for measuring mesopore size distributions. The methods are based on the fact that, according to the Gibbs-Thomson equation, the melting point depression of a liquid confined to a pore is inversely proportional to the pore size. However, aqueous salt solutions, which inherently exist in a broad range of biological porous materials as well as technological applications such as electrolytes, do not melt at a single temperature. This causes artefacts in the pore size distributions extracted by traditional Gibbs-Thomson analysis of NMRC and DSC-TPM data. Bulk aqueous NaCl solutions are known to have a broad distribution of melting points between the eutectic and pure water phase transition points (252-273 K). Here, we hypothesize that, when aqueous NaCl solution (saline) is confined to a small pore, the whole melting point distribution is shifted toward lower temperatures by the value predicted by the Gibbs-Thomson equation. We show that this so-called shifted phase transition distribution (SIDI) approach removes the artefacts arising from the traditional Gibbs-Thomson analysis and gives correct pore size distributions for saline saturated mesoporous silica gel and controlled pore materials analyzed by NMR cryoporometry. Furthermore, we demonstrate that the method can be used for determining pore sizes in collagen-chondroitin sulphate hydrogels resembling the composition of the extracellular matrix of articular cartilage. It is straightforward to apply the SIDI analysis for DSC-TMP data as well.

SPICY: a method for single scan rotating frame relaxometry.

T 1ρ is an NMR relaxation mode that is sensitive to low frequency molecular motions, making it an especially valuable tool in biomolecular research. Here, we introduce a new method, SPICY, for measuring T1ρ relaxation times. In contrast to conventional T1ρ experiments, in which the sequence is repeated many times to determine the T1ρ time, the SPICY sequence allows determination of T1ρ within a single scan, shortening the experiment time remarkably. We demonstrate the method using 1H T1ρ relaxation dispersion experiments. Additionally, we combine the sequence with spatial encoding to produce 1D images in a single scan. We show that T1ρ relaxation times obtained using the single scan approach are in good agreement with those obtained using the traditional experiments.

Molecular-scale evidence of aerosol particle formation via sequential addition of HIO3.

Homogeneous nucleation and subsequent cluster growth leads to the formation of new aerosol particles in the atmosphere. The nucleation of sulfuric acid and organic vapours is thought to be responsible for the formation of new particles over continents, whereas iodine oxide vapours have been implicated in particle formation over coastal regions. The molecular clustering pathways that are involved in atmospheric particle formation have been elucidated in controlled laboratory studies of chemically simple systems, but direct molecular-level observations of nucleation in atmospheric field conditions that involve sulfuric acid, organic or iodine oxide vapours have yet to be reported. Here we present field data from Mace Head, Ireland, and supporting data from northern Greenland and Queen Maud Land, Antarctica, that enable us to identify the molecular steps involved in new particle formation in an iodine-rich, coastal atmospheric environment. We find that the formation and initial growth process is almost exclusively driven by iodine oxoacids and iodine oxide vapours, with average oxygen-to-iodine ratios of 2.4 found in the clusters. On the basis of this high ratio, together with the high concentrations of iodic acid (HIO3) observed, we suggest that cluster formation primarily proceeds by sequential addition of HIO3, followed by intracluster restructuring to I2O5 and recycling of water either in the atmosphere or on dehydration. Our study provides ambient atmospheric molecular-level observations of nucleation, supporting the previously suggested role of iodine-containing species in the formation of new aerosol particles, and identifies the key nucleating compound.

A potential for molecular simulation of compounds with linear moieties.

The harmonic angle bending potential is used in many force fields for (bio)molecular simulation. The force associated with this potential is discontinuous at angles close to 180°, which can lead to numeric instabilities. Angle bending of linear groups, such as alkynes or nitriles, or linear molecules, such as carbon dioxide, can be treated by a simple harmonic potential if we describe the fluctuations as a deviation from a reference position of the central atom, the position of which is determined by the flanking atoms. The force constant for the linear angle potential can be derived analytically from the corresponding force constant in the traditional potential. The new potential is tested on the properties of alkynes, nitriles, and carbon dioxide. We find that the angles of the linear groups remain about 2° closer to 180° using the new potential. The bond and angle force constants for carbon dioxide were tuned to reproduce the experimentally determined frequencies. An interesting finding was that simulations of liquid carbon dioxide under pressure with the new flexible model were stable only when explicitly modeling the long-range Lennard-Jones (LJ) interactions due to the very long-range nature of the LJ interactions (>1.7 nm). In the other tested liquids, we find that a Lennard-Jones cutoff of 1.1 nm yields similar results as the particle mesh Ewald algorithm for LJ interactions. Algorithmic factors influencing the stability of liquid simulations are discussed as well. Finally, we demonstrate that the linear angle potential can be used in free energy perturbation calculations.

Mechanism of the Oxidation of Heptafulvenes to Tropones Studied by Online Mass Spectrometry and Density Functional Theory Calculations.

We have identified the most likely reaction mechanism for oxidizing heptafulvenes to the corresponding tropones by experimental and theoretical investigations. The experimental studies were done by coupling a three-dimensional printed miniaturized reactor with an integrated electrospray ionization needle to a mass spectrometer. Using the experimentally observed ions as a basis, nine alternative reaction pathways were investigated with density functional theory calculations. The lowest energy reaction pathway starts with the formation of an epoxide that is opened upon the addition of a second equivalent of the oxidizing species meta-chloroperoxybenzoic acid. The adduct formed then undergoes a Criegee-like rearrangement to yield a positively charged hemiketal, which on deprotonation dissociates into acetone and tropone. Overall, the reaction mechanism resembles a Hock-like rearrangement.

Molecular understanding of sulphuric acid-amine particle nucleation in the atmosphere.

Nucleation of aerosol particles from trace atmospheric vapours is thought to provide up to half of global cloud condensation nuclei. Aerosols can cause a net cooling of climate by scattering sunlight and by leading to smaller but more numerous cloud droplets, which makes clouds brighter and extends their lifetimes. Atmospheric aerosols derived from human activities are thought to have compensated for a large fraction of the warming caused by greenhouse gases. However, despite its importance for climate, atmospheric nucleation is poorly understood. Recently, it has been shown that sulphuric acid and ammonia cannot explain particle formation rates observed in the lower atmosphere. It is thought that amines may enhance nucleation, but until now there has been no direct evidence for amine ternary nucleation under atmospheric conditions. Here we use the CLOUD (Cosmics Leaving OUtdoor Droplets) chamber at CERN and find that dimethylamine above three parts per trillion by volume can enhance particle formation rates more than 1,000-fold compared with ammonia, sufficient to account for the particle formation rates observed in the atmosphere. Molecular analysis of the clusters reveals that the faster nucleation is explained by a base-stabilization mechanism involving acid-amine pairs, which strongly decrease evaporation. The ion-induced contribution is generally small, reflecting the high stability of sulphuric acid-dimethylamine clusters and indicating that galactic cosmic rays exert only a small influence on their formation, except at low overall formation rates. Our experimental measurements are well reproduced by a dynamical model based on quantum chemical calculations of binding energies of molecular clusters, without any fitted parameters. These results show that, in regions of the atmosphere near amine sources, both amines and sulphur dioxide should be considered when assessing the impact of anthropogenic activities on particle formation.

Computational Study on the Effect of Hydration on New Particle Formation in the Sulfuric Acid/Ammonia and Sulfuric Acid/Dimethylamine Systems.

The formation of new particles through condensation from the gas phase is an important source of atmospheric aerosols. The properties of the electrically neutral clusters formed in the very first steps of the condensation process are, however, not directly observable by experimental means. We present here electronic structure calculations on the hydrates of clusters of three molecules of sulfuric acid and three molecules of ammonia or dimethylamine. On the basis of the results of these new calculations together with previously published material we simulate the influence of hydration on the dynamic processes involved in particle formation. Most strongly affected by hydration and most important as a mediator for the effect on particle formation rates are the evaporation rates of clusters. The results give an estimate of the sensitivity of the atmospheric particle formation rate for humidity. The particle formation rate can change approximately two orders of magnitude in either direction due to hydration; the net effect, however, is highly dependent on the exact conditions.

Effect of Conformers on Free Energies of Atmospheric Complexes.

In this article we show how to calculate free energies for atmospherically relevant complexes when multiple conformers and/or isomers are present. We explain why the thermal averaging methods used in several published works are incorrect. On the basis of our two sample cases, the sulfuric acid-pinic acid complex and the (H2SO4)3(NH3)3(H2O)4 cluster, we provide numerical evidence that the use of these incorrect formulas can result in errors larger than 1 kcal/mol. We recommend that if vibrational frequencies and thus Gibbs free energies of the individual conformers are unavailable, one should not attempt to correct for the presence of multiple conformers and instead use only the global minimum conformers for both reactants and products. On the contrary, if the free energies for the conformers are calculated for both reactants and products, their effect can be accounted for by the statistical mechanical methods presented in this article.

Structures, Hydration, and Electrical Mobilities of Bisulfate Ion-Sulfuric Acid-Ammonia/Dimethylamine Clusters: A Computational Study.

Despite the well-established role of small molecular clusters in the very first steps of atmospheric particle formation, their thermochemical data are still not completely available due to limitation of the experimental techniques to treat such small clusters. We have investigated the structures and the thermochemistry of stepwise hydration of clusters containing one bisulfate ion, sulfuric acid, base (ammonia or dimethylamine), and water molecules using quantum chemical methods. We found that water facilitates proton transfer from sulfuric acid or the bisulfate ion to the base or water molecules, and depending on the hydration level, the sulfate ion was formed in most of the base-containing clusters. The calculated hydration energies indicate that water binds more strongly to ammonia-containing clusters than to dimethylamine-containing and base-free clusters, which results in a wider hydrate distribution for ammonia-containing clusters. The electrical mobilities of all clusters were calculated using a particle dynamics model. The results indicate that the effect of humidity is negligible on the electrical mobilities of molecular clusters formed in the very first steps of atmospheric particle formation. The combination of the results of this study with those previously published on the hydration of neutral clusters by our group provides a comprehensive set of thermochemical data on neutral and negatively charged clusters containing sulfuric acid, ammonia, or dimethylamine.

Hydration of atmospherically relevant molecular clusters: computational chemistry and classical thermodynamics.

Formation of new particles through clustering of molecules from condensable vapors is a significant source for atmospheric aerosols. The smallest clusters formed in the very first steps of the condensation process are, however, not directly observable by experimental means. We present here a comprehensive series of electronic structure calculations on the hydrates of clusters formed by up to four molecules of sulfuric acid, and up to two molecules of ammonia or dimethylamine. Though clusters containing ammonia, and certainly dimethylamine, generally exhibit lower average hydration than the pure acid clusters, populations of individual hydrates vary widely. Furthermore, we explore the predictions obtained using a thermodynamic model for the description of these hydrates. The similar magnitude and trends of hydrate formation predicted by both methods illustrate the potential of combining them to obtain more comprehensive models. The stabilization of some clusters relative to others due to their hydration is highly likely to have significant effects on the overall processes that lead to formation of new particles in the atmosphere.

Dipolar Relaxation of Water Protons in the Vicinity of a Collagen-like Peptide.

Quantitative magnetic resonance imaging is one of the few available methods for noninvasive diagnosis of degenerative changes in articular cartilage. The clinical use of the imaging data is limited by the lack of a clear association between structural changes at the molecular level and the measured magnetic relaxation times. In anisotropic, collagen-containing tissues, such as articular cartilage, the orientation dependency of nuclear magnetic relaxation can obscure the content of the images. Conversely, if the molecular origin of the phenomenon would be better understood, it would provide opportunities for diagnostics as well as treatment planning of degenerative changes in these tissues. We study the magnitude and orientation dependence of the nuclear magnetic relaxation due to dipole-dipole coupling of water protons in anisotropic, collagenous structures. The water-collagen interactions are modeled with molecular dynamics simulations of a small collagen-like peptide dissolved in water. We find that in the vicinity of the collagen-like peptide, the dipolar relaxation of water hydrogen nuclei is anisotropic, which can result in orientation-dependent relaxation times if the water remains close to the peptide. However, the orientation-dependency of the relaxation is different from the commonly observed magic-angle phenomenon in articular cartilage MRI.

Rational design of biomimetic molecularly imprinted materials: theoretical and computational strategies for guiding nanoscale structured polymer development.

In principle, molecularly imprinted polymer science and technology provides a means for ready access to nano-structured polymeric materials of predetermined selectivity. The versatility of the technique has brought it to the attention of many working with the development of nanomaterials with biological or biomimetic properties for use as therapeutics or in medical devices. Nonetheless, the further evolution of the field necessitates the development of robust predictive tools capable of handling the complexity of molecular imprinting systems. The rapid growth in computer power and software over the past decade has opened new possibilities for simulating aspects of the complex molecular imprinting process. We present here a survey of the current status of the use of in silico-based approaches to aspects of molecular imprinting. Finally, we highlight areas where ongoing and future efforts should yield information critical to our understanding of the underlying mechanisms sufficient to permit the rational design of molecularly imprinted polymers.

Warfarin: an environment-dependent switchable molecular probe.

The complex nature of the structure of the anticoagulant warfarin is reflected in the diversity of binding modes observed in warfarin-protein recognition systems. A series of theoretical, (1)H-NMR and steady state and time resolved fluorescence spectroscopic studies, have been used to establish correlations between the molecular environment provided by various solvent systems and the relative concentrations of the various members of warfarin's ensemble of isomers. A consequence of these observations is that the judicious choice of solvent system or molecular environment of warfarin allows for manipulation of the position of the equilibrium between isomeric structures such as the hemiacetal and open phenol-keto forms, the latter even possible in a deprotonated form, where in each case unique spectroscopic properties are exhibited by the respective structures. Collectively, warfarin's capacity to adapt its structure as a function of environment in conjunction with the fluorescence behaviours of the various isomers together provide an environment-dependent molecular switch with reporter properties, which allows for the simultaneous detection of warfarin in different states with lifetimes spanning the range < 0.10-5.5 ns. These characteristics are here used to examine warfarin binding domains in a series of materials (solvents, protein, inorganic matrix and synthetic polymer). Moreover, these studies demonstrate the potential for using warfarin, or other switchable analogues thereof, as a tool for studying molecular level characteristics, for example local dielectricity.

An Intuitively Understandable Quality Measure for Theoretical Vibrational Spectra.

Computational chemistry has become a central tool in spectroscopic studies in most of chemical science. The quality of a calculated vibrational spectrum is commonly expressed as the deviation of the peak position from the experimental reference. With the increasing application of vibrational spectroscopy to complex (biological) systems, this is likely not sustainable. Here we present a quality measure for theoretical vibrational spectra based on matching the spectra to a reference database with the help of correlation coefficients. This approach can easily be applied to large sets of data and complex spectra without easily identifiable peak positions. We demonstrate this on a database of infrared spectra of 670 compounds using six different theoretical (DFT and force field) methods. Most importantly, it is intuitively understandable by both theoreticians and experimentalists.

Theoretical Infrared Spectra: Quantitative Similarity Measures and Force Fields.

Infrared spectroscopy can provide significant insight into the structures and dynamics of molecules of all sizes. The information that is contained in the spectrum is, however, often not easily extracted without the aid of theoretical calculations or simulations. We present here the calculation of the infrared spectra of a database of 703 gas phase compounds with four different force fields (CGenFF, GAFF-BCC, GAFF-ESP, and OPLS) using normal-mode analysis. Modern force fields increasingly use virtual sites to describe, e.g., lone-pair electrons or the σ-holes on halogen atoms. This requires some adaptation of code to perform normal-mode analysis of such compounds, the implementation of which into the GROMACS software is briefly described as well. For the quantitative comparison of the obtained spectra with experimental reference data, we discuss the application of two different statistical correlation coefficients, Pearson and Spearman. The advantages and drawbacks of the different methods of comparison are discussed, and we find that both methods of comparison give the same overall picture, showing that present force field methods cannot match the performance of quantum chemical methods for the calculation of infrared spectra.

How warfarin's structural diversity influences its phospholipid bilayer membrane permeation.

The role of the structural diversity of the widely used anticoagulant drug warfarin on its distribution in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayer membranes was investigated using a series of both restrained (umbrella sampling) and unrestrained molecular dynamics simulations. Data collected from unrestrained simulations revealed favorable positions for neutral isomers of warfarin, the open side chain form (OCO), and the cyclic hemiketal (CCO), along the bilayer normal close to the polar headgroup region and even in the relatively distant nonpolar lipid tails. The deprotonated open side chain form (DCO) was found to have lower affinity for the DOPC bilayer membrane relative to the neutral forms, with only a small fraction interacting with the membrane, typically within the polar headgroup region. The conformation of OCO inside the lipid bilayer was found to be stabilized by intramolecular hydrogen bonding thereby mimicking the structure of CCO. Differences in free energies, for positions of OCO and CCO inside the bilayer membrane, as compared to positions in the aqueous phase, were -97 and -146 kJ·mol(-1). Kinetic analysis based on the computed free energy barriers reveal that warfarin will diffuse through the membranes within hours, in agreement with experimental results on warfarin's accumulation in the plasma, thus suggesting a passive diffusion mechanism. We propose that this membrane transport may be an isomerization-driven process where warfarin adapts to the various local molecular environments encountered under its journey through the membrane. Collectively, these results improve our understanding of the influence of warfarin's structural diversity on the drug's distribution and bioavailability, which in turn may provide insights for developing new formulations of this important pharmaceutical to better address its narrow therapeutic window.

Theoretical and computational strategies for rational molecularly imprinted polymer design.

The further evolution of molecularly imprinted polymer science and technology necessitates the development of robust predictive tools capable of handling the complexity of molecular imprinting systems. A combination of the rapid growth in computer power over the past decade and significant software developments have opened new possibilities for simulating aspects of the complex molecular imprinting process. We present here a survey of the current status of the use of in silico-based approaches to aspects of molecular imprinting. Finally, we highlight areas where ongoing and future efforts should yield information critical to our understanding of the underlying mechanisms sufficient to permit the rational design of molecularly imprinted polymers.

(H2O)10 and (H2O)12 on a virtual metal surface: the growth of ice.

(H2O)10 and (H2O)12 are used to investigate the growth of ice on metal surfaces with hexagonal symmetry. The model of the virtual metal surface was used to separate the electronic structure of the metal from that of the water cluster while maintaining the geometric constraints imposed by the metal surface on the water cluster. To complement the ab initio calculations on the water cluster, an additional multicenter analysis was done to analyze the hydrogen bonds within the clusters. These calculations suggested that the water bilayer structure adjacent to the virtual metal surface effectively shields the growing ice crystal from the metal surface.