Hydrophobic interaction of tetrapropylammonium ions in water: A neutron diffraction and reverse Monte Carlo study

Hydration of biologically relevant tetramethylammonium cation by neutron scattering and molecular dynamics

Neutron scattering and molecular dynamics studies were performed on a concentrated aqueous tetramethylammonium (TMA) chloride solution to gain insight into the hydration shell structure of TMA, which is relevant for understanding its behavior in biological contexts of, e.g., properties of phospholipid membrane headgroups or interactions between DNA and histones. Specifically, neutron diffraction with isotopic substitution experiments were performed on TMA and water hydrogens to extract the specific correlation between hydrogens in TMA (HTMA) and hydrogens in water (HW). Classical molecular dynamics simulations were performed to help interpret the experimental neutron scattering data. Comparison of the hydration structure and simulated neutron signals obtained with various force field flavors (e.g. overall charge, charge distribution, polarity of the CH bonds and geometry) allowed us to gain insight into how sensitive the TMA hydration structure is to such changes and how much the neutron signal can capture them. We show that certain aspects of the hydration, such as the correlation of the hydrogen on TMA to hydrogen on water, showed little dependence on the force field. In contrast, other correlations, such as the ion-ion interactions, showed more marked changes. Strikingly, the neutron scattering signal cannot discriminate between different hydration patterns. Finally, ab initio molecular dynamics was used to examine the three-dimensional hydration structure and thus to benchmark force field simulations. Overall, while neutron scattering has been previously successfully used to improve force fields, in the particular case of TMA we show that it has only limited value to fully determine the hydration structure, with other techniques such as ab initio MD being of a significant help.

Physico-chemical properties of aqueous drug solutions: From the basic thermodynamics to the advanced experimental and simulation results

The physical chemical properties of aqueous solutions of model compounds are illustrated in relation to hydration and solubility issues by using three perspectives: thermodynamic, spectroscopic and molecular dynamics simulations. The thermodynamic survey of the fundamental backgrounds of concentration dependence and experimental solubility results show some peculiar behavior of aqueous solutions with several types of similar solutes. Secondly, the use of a variety of experimental spectroscopic devices, operating under different experimental conditions of dimension and frequency, has produced a large amount of structural and dynamic data on aqueous solutions showing the richness of the information produced, depending on where and how the experiment is carried out. Finally, the use of molecular dynamics computational work is presented to highlight how the different types of solute functional groups and surface topologies organize adjacent water molecules differently. The highly valuable contribution of computer simulation studies in providing molecular explanations for experimental deductions, either of a thermodynamic or spectroscopic nature, is shown to have changed the current knowledge of many aqueous solution processes. While this paper is intended to provide a collective view on the latest literature results, still the presentation aims at a tutorial explanation of the potentials of the three methodologies in the field of aqueous solutions of pharmaceutical molecules.

Do guanidinium and tetrapropylammonium ions specifically interact with aromatic amino acid side chains?

Many ions are known to affect the activity, stability, and structural integrity of proteins. Although this effect can be generally attributed to ion-induced changes in forces that govern protein folding, delineating the underlying mechanism of action still remains challenging because it requires assessment of all relevant interactions, such as ion-protein, ion-water, and ion-ion interactions. Herein, we use two unnatural aromatic amino acids and several spectroscopic techniques to examine whether guanidinium chloride, one of the most commonly used protein denaturants, and tetrapropylammonium chloride can specifically interact with aromatic side chains. Our results show that tetrapropylammonium, but not guanidinium, can preferentially accumulate around aromatic residues and that tetrapropylammonium undergoes a transition at ∼1.3 M to form aggregates. We find that similar to ionic micelles, on one hand, such aggregates can disrupt native hydrophobic interactions, and on the other hand, they can promote α-helix formation in certain peptides.

A neutron scattering and modelling study of aqueous solutions of tetramethylammonium and tetrapropylammonium bromide

Tetrapropylammonium ions cluster and have less structured water arrangements whereas tetramethylammonium ions occur as single molecules with tetrahedral water arrangements.

Aqueous solutions of tetraalkylammonium halides: ion hydration, dynamics and ion-ion interactions in light of steric effects

Molecular simulations have allowed us to probe the atomic details of aqueous solutions of tetramethylammonium (TMA) and tetrabutylammonium (TBA) bromide, across a wide range of concentrations (0.5 to 3-4 molal). We highlight the space-filling (TMA(+)) versus penetrable (TBA(+)) nature of these polyatomic cations and its consequence for ion hydration, ion dynamics and ion-ion interactions. A well-established hydration is seen for both TMA(+) and TBA(+) throughout the concentration range studied. A clear penetration of water molecules, as well as counterions, between the hydrocarbon arms of TBA(+), which remain in an extended configuration, is seen. Global rotation of individual TBA(+) points towards isolated rather than aggregated ions (from dilute up to 1 m concentration). Only for highly concentrated solutions, in which inter-penetration of adjacent TBA(+)s cannot be avoided, does the rotational time increase dramatically. From both structural and dynamic data we conclude that there is absence of hydrophobicity-driven cation-cation aggregation in both TMABr and TBABr solutions studied. The link between these real systems and the theoretical predictions for spherical hydrophobic solutes of varying size does not seem straightforward.

Simulation and Neutron Diffraction Studies of Small Biomolecules in Water

Modern biophysics has benefited greatly from the use of X-ray and neutron diffraction from ordered single crystals of proteins and other macromolecules to give highly detailed pictures of these molecules in the solid state. However, the most biologically relevant environments for these molecules are liquid solutions, and their liquid state properties are sensitive to details of the liquid structuring. The best experimental method for studying such structuring is also neutron diffraction, but of course, the inherent disorder of the liquid state means that these experiments cannot hope to achieve the level of informational detail available from single crystal diffraction. Nonetheless, recent advances in neutron beam intensity, beam stability, and detector sensitivity mean that it should be possible, at least in principle, to use such measurements to extract information about structuring in much more complex systems than have previously been studied. We describe a series of neutron diffraction studies of isotopically labeled molecules in aqueous solution which, when combined with results from computer simulations, can be used to extract conformational information of the hydration of the molecules themselves, essentially opening up new avenues of investigation in structural biology.

Observation of pyridine aggregation in aqueous solution using neutron scattering experiments and MD simulations

Neutron diffraction with isotopic substitution (NDIS) experiments have been used to examine the structuring of aqueous solutions of pyridine. A new method is described for extracting the structure factors relating to intermolecular correlations from neutron scattering experiments on liquid solutions of complex molecular species. This approach performs experiments at different concentrations and exploits the intramolecular coordination number concentration invariance (ICNCI) to separate the internal and intermolecular contributions to the total intensities. The ability of this method to deconvolute molecular and intermolecular correlations is tested and demonstrated using simulated neutron scattering results predicted from molecular dynamics simulations of aqueous solutions of the polyatomic solute pyridine in which the inter- and intramolecular terms are known. The method is then implemented using neutron scattering measurements on solutions of pyridine. The results confirm that pyridine shows a significant propensity to aggregate in solution and demonstrate the prospects for the future application of the ICNCI approach to the study of large polyatomic solutes using next-generation neutron sources and detectors.

Modeling Spontaneous Formation of Precursor Nanoparticles in Clear-Solution Zeolite Synthesis

We present a lattice model describing the formation of silica nanoparticles in the early stages of the clear-solution templated synthesis of silicalite-1 zeolite. Silica condensation/hydrolysis is modeled by a nearest-neighbor attraction, while the electrostatics are represented by an orientation-dependent, short-range interaction. Using this simplified model, we show excellent qualitative agreement with published experimental observations. The nanoparticles are identified as a metastable state, stabilized by electrostatic interactions between the negatively charged silica surface and a layer of organic cations. Nanoparticle size is controlled mainly by the solution pH, through nanoparticle surface charge. The size and concentration of the charge-balancing cation are found to have a negligible effect on nanoparticle size. Increasing the temperature allows for further particle growth by Ostwald ripening. We suggest that this mechanism may play a role in the growth of zeolite crystals.