In this work, the collective structure of aqueous solutions of ionic liquids was studied by means of molecular dynamics simulations. Various concentrations of 1-butyl-3-methyl-imidazolium tetrafluoroborate and TIP3P water were simulated at the very same size of the simulation box. For the analysis, the ternary system cation/anion/water was subdivided into binary networks. The local structure of each of these six networks is investigated by atom-atom radial distribution functions as well as by the so-called g coefficients, which reveal the mutual orientation of the network constituting partners. Furthermore, the collective structure of the whole samples was characterized by the contribution of each species to the static dielectric constant epsilon(omega=0) and to the Kirkwood G(K) factor. The combination of the analysis tools mentioned above provides knowledge about the cross-linking of the ionic species with the dipolar water. Thereby, the interplay between charge-charge and hydrogen bond networks is analyzed in detail.

The ionic liquids BMIM+I-, BMIM+BF4-, and BMIM+PF6- were simulated by means of the molecular dynamics method over a time period of more than 100 ns. Besides the common structural analysis, e.g., radial distribution functions and three dimensional occupancy plots, a more sophisticated orientational analysis was performed. The angular correlation functions g(00)110(r) and g(00)101(r) are the first distance dependent coefficients of the pairwise orientational distribution function g(rij,Omega1,Omega2,Omega12). These functions help to interpret the three dimensional plot and reveal interesting insights into the local structure of the analyzed ionic liquids. Furthermore, the collective network of ionic liquids can be characterized by the Kirkwood factor Gkappa(r) [J. Chem. Phys. 7, 911 (1939)]. The short-range behavior (r<10 A) of this factor may be suitable to predict the water miscibility of the ionic liquid. The long-range limit of Gkinfinity is below 1 which demonstrates the strongly coupled nature of the ionic liquid networks. In addition, this factor relates the orientational structure and the dielectric properties of the ionic liquids. The static dielectric constant epsilon(omega=0) for the simulated system is 8.9-9.5. Since in ionic liquids the very same molecule contributes to the total dipole moment as well as carries a net charge, a small, but significant contribution of the cross term between the total dipole moment and the electric current to epsilon(omega=0) is observed.

We report molecular dynamics simulations of three globular proteins: ubiquitin, apo-calbindin D(9K), and the C-terminal SH2 domain of phospholipase C-gamma1 in explicit water. The proteins differ in their overall charge and fold type and were chosen to represent to some degree the structural variability found in medium-sized proteins. The length of each simulation was at least 15 ns, and larger than usual solvent boxes were used. We computed radial distribution functions, as well as orientational correlation functions about the surface residues. Two solvent shells could be clearly discerned about charged and polar amino acids. Near apolar amino acids the water density near such residues was almost devoid of structure. The mean residence time of water molecules was determined for water shells about the full protein, as well as for water layers about individual amino acids. In the dynamic properties, two solvent shells could be characterized as well. However, by comparison to simulations of pure water it could be shown that the influence of the protein reaches beyond 6 A, i.e., beyond the first two shells. In the first shell (r < or =3.5 A), the structural and dynamical properties of solvent waters varied considerably and depended primarily on the physicochemical properties of the closest amino acid side chain, with which the waters interact. By contrast, the solvent properties seem not to depend on the specifics of the protein studied (such as the net charge) or on the secondary structure element in which an amino acid is located. While differing considerably from the neat liquid, the properties of waters in the second solvation shell (3.5< r < or =6 A) are rather uniform; a direct influence from surface amino acids are already mostly shielded.

We report molecular dynamics (MD) simulations of three protein-water systems (ubiquitin, apo-calbindin D(9K), and the C-terminal SH2 domain of phospholipase C-gamma1), from which we compute the dielectric properties of the solutions. Since two of the proteins studied have a net charge, we develop the necessary theory to account for the presence of charged species in a form suitable for computer simulations. In order to ensure convergence of the time correlation functions needed for the analysis, the minimum length of the MD simulations was 20 ns. The system sizes (box length, number of waters) were chosen so that the resulting protein concentrations are comparable to experimental conditions. A dielectric component analysis was carried out to analyze the contributions from protein and water to the frequency-dependent dielectric susceptibility chi(omega) of the solutions. Additionally, an even finer decomposition into protein, two solvation shells, and the remaining water (bulk water) was carried out. The results of these dielectric decompositions were used to study protein solvation at mesoscopic resolution, i.e., in terms of protein, first and second solvation layers, and bulk water. This study, therefore, complements the structural and dynamical analyses at molecular resolution that are presented in the companion paper. The dielectric component contributions from the second shell and bulk water are very similar in all three systems. We find that the proteins influence the dielectric properties of water even beyond the second solvation shell, in agreement with what was observed for the mean residence times of water molecules in protein solutions. By contrast, the protein contributions, as well as the contributions of the first solvation shell, are system specific. Most importantly, the protein and the first water shell around ubiquitin and apo-calbindin are anticorrelated, whereas the first water shell around the SH2 domain is positively correlated.