Simulation of interfaces between room temperature ionic liquids and other liquids

The structure and properties of the interfaces between the room temperature ionic liquid dimethylimidazolium chloride ([dmim]Cl) and different Lennard-Jones fluids and between ionic liquid and water have been studied by molecular dynamics simulations, and compared to the ionic liquid-vapour interface. Two contrasting types of interface were investigated, thermodynamically stable interfaces between ionic liquid and vapour and between ionic liquid and Lennard-Jones fluids, and diffusing interfaces between miscible phases of different compositions involving water. The density profiles of different species through the interface are presented. The cations and water molecules near the former type of interface are aligned relative to the surface, but no orientational preference was found near or in the broad diffusing interface. The ionic liquid has a negative electrostatic potential relative to vapour or Lennard-Jones fluid, but is more positive than pure water. This contrast is explained in terms of the relative importance of orientation and concentration differences in the two types of interface.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

Structural basis of substrate selectivity in the glycerol-3-phosphate: phosphate antiporter GlpT.

Major facilitators represent the largest superfamily of secondary active transporter proteins and catalyze the transport of an enormous variety of small solute molecules across biological membranes. However, individual superfamily members, although they may be architecturally similar, exhibit strict specificity toward the substrates they transport. The structural basis of this specificity is poorly understood. A member of the major facilitator superfamily is the glycerol-3-phosphate (G3P) transporter (GlpT) from the Escherichia coli inner membrane. GlpT is an antiporter that transports G3P into the cell in exchange for inorganic phosphate (Pi). By combining large-scale molecular-dynamics simulations, mutagenesis, substrate-binding affinity, and transport activity assays on GlpT, we were able to identify key amino acid residues that confer substrate specificity upon this protein. Our studies suggest that only a few amino acid residues that line the transporter lumen act as specificity determinants. Whereas R45, K80, H165, and, to a lesser extent Y38, Y42, and Y76 contribute to recognition of both free Pi and the phosphate moiety of G3P, the residues N162, Y266, and Y393 function in recognition of only the glycerol moiety of G3P. It is the latter interactions that give the transporter a higher affinity to G3P over Pi.

Interaction of Room Temperature Ionic Liquid Solutions with a Cholesterol Bilayer

Water solutions of representative (IC(4)mim][Cl] and [C(4)mim][Tf2N] room temperature ionic liquids (ILs) in contact with a neutral lipid bilayer made of cholesterol molecules has been investigated by molecular dynamics simulations based on an empirical force field model. The results show that both ILs display selective adsorption at the water-cholesterol interface, with partial inclusion of ions into the bilayer. In the case Of [C(4)mim][Cl], the adsorption of ions at the water-cholesterol interface is limited by a sizable bulk solubility of the IL, driven by the high water affinity of [Cl](-). The relatively low Solubility Of [C(4)mim][Tf2N], instead, gives rise to a nearly complete segregation of the IL component on the bilayer, altering its volume, compressibility, and electrostatic environment. The computational results display important similarities to the results of recent experimental measurements for ILs in contact with phospholipid model membranes (see Evans, K. O. Int. J. Mol. Sci. 2008, 9, 498-511 and references therein).

Structural motifs of cholesterol nanoparticles

The growth sequence of gas-phase cholesterol clusters (Ch(N)) with up to N=36 molecules has been investigated by atomistic simulation based on an empirical force field model. The results of long annealings from high temperature show that the geometric motifs characterizing the structure of pure cholesterol crystals already appear in nanometric aggregates. In all clusters molecules tend to align along a common direction. For cluster sizes above the smallest ones, dispersion interactions among the hydrocarbon body and tails of cholesterol cooperate with hydrogen bonding to give rise to a bilayer structure. Analysis of snapshots from the annealing shows that the condensation of hydrogen bonds into a connected network of rings and chains is an important step in the self-organization of cholesterol clusters. The effect of solvation on the equilibrium properties of medium-size aggregates is investigated by short molecular dynamics simulations for the N=30 and N=40 clusters in water at near ambient conditions and in supercritical carbon dioxide at T=400 K.

Structure of warm dense matter via angularly resolved x-ray scatter

In this paper we describe experimental results on angularly resolved x-ray scatter from a sample of warm dense aluminium that has been created by double sided laser-driven shock compression. The experiment was carried out on the Central Laser Facility of the Rutherford Appleton Laboratory, using the VULCAN laser. The form of the angularly resolved scatter cross-section was compared with predictions based on a series of molecular dynamics simulations with an embedded atom potential, a Yukakwa potential and a bare Coulomb potential. The importance of screening is evident from the comparison and the embedded atom model seems to match experiment better than the Yukawa potential.

Two Arginine-Glutamate Ionic Locks Near the Extracellular Surface of FFAR1 Gate Receptor Activation

Activation of a number of class A G protein-coupled receptors (GPCRs) is thought to involve two molecular switches, a rotamer toggle switch within the transmembrane domain and an ionic lock at the cytoplasmic surface of the receptor; however, the mechanism by which agonist binding changes these molecular interactions is not understood. Importantly, 80% of GPCRs including free fatty acid receptor 1 (FFAR1) lack the complement of amino acid residues implicated in either or both of these two switches; the mechanism of activation of these GPCRs is therefore less clear. By homology modeling, we identified two Glu residues (Glu-145 and Glu-172) in the second extracellular loop of FFAR1 that form putative interactions individually with two transmembrane Arg residues (Arg-183(5.39) and Arg-258(7.35)) to create two ionic locks. Molecular dynamics simulations showed that binding of agonists to FFAR1 leads to breakage of these Glu-Arg interactions. In mutagenesis experiments, breakage of these two putative interactions by substituting Ala for Glu-145 and Glu-172 caused constitutive receptor activation. Our results therefore reveal a molecular switch for receptor activation present on the extracellular surface of FFAR1 that is broken by agonist binding. Similar ionic locks between the transmembrane domains and the extracellular loops may constitute a mechanism common to other class A GPCRs also.

Atomistic Insights into Rhodopsin Activation from a Dynamic Model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Mechanism of Activation of a G Protein-coupled Receptor, the Human Cholecystokinin-2 Receptor

G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.

Solid and liquid charge-transfer complex formation between 1-methylnaphthalene and 1-alkyl-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide ionic liquids

Liquid charge-transfer (CT) complexes were observed to form on contacting electron-rich aromatics with electron withdrawing group appended 1-alkyl-4-cyanopyridinium ionic liquids (ILs). Cooling below the melting point of the ionic liquid resulted in crystallisation of ionic liquid from the complex for 2-cyano and 3-cyano pyridinium isomers and in the formation of a 1 : 1 IL : aromatic crystalline CT-complex with the 4-cyanopyridinium isomer. The liquid structure of a 1 : 1 mixture of 1-methyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl} imide with 1-methylnaphthalene has been probed by neutron diffraction experiments and molecular dynamics simulations. A high degree of correlation between the experimental data and the simulations was found with a significant displacement of the anions from around the cation by the aromatic species and the resulting structure having pi-pi stacks between the cations and the aromatic.

A crystallographic and modelling study of a human telomeric RNA (TERRA) quadruplex

DNA telomeric repeats in mammalian cells are transcribed to guanine-rich RNA sequences, which adopt parallel-stranded G-quadruplexes with a propeller-like fold. The successful crystallization and structure analysis of a bimolecular human telomeric RNA G-quadruplex, folded into the same crystalline environment as an equivalent DNA oligonucleotide sequence, is reported here. The structural basis of the increased stability of RNA telomeric quadruplexes over DNA ones and their preference for parallel topologies is described here. Our findings suggest that the 2'-OH hydroxyl groups in the RNA quadruplex play a significant role in redefining hydration structure in the grooves and the hydrogen bonding networks. The preference for specific nucleotides to populate the C3'-endo sugar pucker domain is accommodated by alterations in the phosphate backbone, which leads to greater stability through enhanced hydrogen bonding networks. Molecular dynamics simulations on the DNA and RNA quadruplexes are consistent with these findings. The computations, based on the native crystal structure, provide an explanation for RNA G-quadruplex ligand binding selectivity for a group of naphthalene diimide ligands as compared to the DNA G-quadruplex.

Computational Verification of Two Universal Relations for Simple Ionic Liquids. Kinetic Properties of a Model 2:1 Molten Salt

Two semianalytical relations [Nature, 1996, 381, 137 and Phys. Rev. Lett. 2001, 87, 245901] predicting dynamical coefficients of simple liquids on the basis of structural properties have been tested by extensive molecular dynamics simulations for an idealized 2:1 model molten salt. In agreement with previous simulation studies, our results support the validity of the relation expressing the self-diffusion coefficient as a Function of the radial distribution functions for all thermodynamic conditions such that the system is in the ionic (ie., fully dissociated) liquid state. Deviations are apparent for high-density samples in the amorphous state and in the low-density, low-temperature range, when ions condense into AB(2) molecules. A similar relation predicting the ionic conductivity is only partially validated by our data. The simulation results, covering 210 distinct thermodynamic states, represent an extended database to tune and validate semianalytical theories of dynamical properties and provide a baseline for the interpretation of properties of more complex systems such as the room-temperature ionic liquids.

Effect of hydrophobic nanopatches within an ionic surface on the structure of liquids

The structures of liquid water and isopropanol have been studied as a function of the size of a hydrophobic patch present in a model hydrophilic surface via molecular dynamics simulations. A significant anisotropy extending into the first few solvent layers is found over the patch which suggests implications for many real-world systems in which nanoscale heterogeneity is found.

Nonadiabatic forces in ion-solid interactions:The initial stages of radiation damage

The Born-Oppenheimer approximation is the keystone for molecular dynamics simulations of radiation damage processes; however, actual materials response involves nonadiabatic energy exchange between nuclei and electrons. In this work, time dependent density functional theory is used to calculate the electronic excitations produced by energetic protons in Al. We study the influence of these electronic excitations on the interatomic forces and find that they differ substantially from the adiabatic case, revealing a nontrivial connection between electronic and nuclear stopping that is absent in the adiabatic case. These results unveil new effects in the early stages of radiation damage cascades.

On thin ice: surface order and disorder during pre-melting

The effect of temperature on the structure of the ice Ih (0001) surface is considered through a series of molecular dynamics simulations on an ice slab. At relatively low temperatures (200K) a small fraction of surface self-interstitials (i.e. admolecules) appear that are formed exclusively from molecules leaving the outermost bilayer. At higher temperatures (ca. 250 K), vacancies start to appear in the inner part of the outermost bilayer exposing the underlying bilayer and providing sites with a high concentration of the dangling hydrogen bonds. Around 250-260 K aggregates of molecules formed on top of the outermost bilayer from self-interstitials become more mobile and have diffusivities approaching that of liquid water. At similar to 270-280 K the inner bilayer of one surface noticeably destructures and it appears that at above 285 K both surfaces are melting. The observed disparity in the onset of melting between the two sides of the slab is rationalised by considering the relationship between surface energy and the spatial distribution of protons at the surface; thermodynamic stability is conferred on the surface by maximising separations between dangling protons at the crystal exterior. Local hotspots associated with a high dangling proton density are suggested to be susceptible to pre-melting and may be more efficient at trapping species at the external surface than regions with low concentrations of protons thus potentially helping ice particles to catalyse reactions. A preliminary conclusion of this work is that only about 10-20 K below the melting temperature of the particular water potential employed is major disruption of the crystalline lattice noted which could be interpreted as being "liquid", the thickness of this film being about a nanometre.