Atomistic Characterization of Healthy and Damaged Hair Surfaces: A Molecular Dynamics Simulation Study of Fatty Acids on Protein Layer.

Amid the bourgeoning demand for in-silico designed, environmentally sustainable, and highly effective hair care formulations, a growing interest is evident in the exploration of realistic computational model for the hair surface. In this work, we present an atomistic model for the outermost layer of the hair surface derived through molecular dynamics simulations, which comprises 18-Methyleicosanoic acid (18-MEA) fatty acid chains covalently bound onto the keratin-associated protein 10-4 (KAP10-4) at a spacing distance of ~1 nm. Remarkably, this hair surface model facilitates the inclusion of free fatty acids (free 18-MEA) into the gaps between chemically bound 18-MEA chains, up to a maximum number that results in a packing density of 0.22 nm2 per fatty acid molecule, consistent with the optimal spacing identified through free energy analysis. Atomistic insights are provided for the organization of fatty acid chains, structural features, and interaction energies on protein-inclusive hair surface models with varying amounts of free 18-MEA (FMEA) depletion, as well as varying degrees of anionic cysteic acid from damaged bound 18-MEA (BMEA), under both dry and wet conditions. In the presence of FMEA and water, the fatty acid chains in a pristine hair surface prefers to adopt a thermodynamically favored extended chain conformation, forming a thicker protective layer (~3 nm) on the protein surface. Our simulation results reveal that, while the depletion of FMEA can induce a pronounced impact on the thickness, tilt angle, and order parameters of fatty acid chains, the removal of BMEA has a marked effect on water penetration. There is a "sweet spot" spacing between the 18-MEA whereby damaged hair surface properties can be reinstated by replenishing FMEA. Through the incorporation of the protein layer and free fatty acids, the hair surface models presented in this study enables a realistic representation of the intricate details within the hair epicuticle, facilitating a molecular scale assessment of surface properties during the formulation design process.

Mobility and association of ions in aqueous solutions: the case of imidazolium based ionic liquids.

The mobility and the mechanism of ion pairing of 1,1 electrolytes in aqueous solutions were investigated systematically on nine imidazolium based ionic liquids (ILs) from 1-methylimidazolium chloride, [MIM][Cl], to 1-dodecyl-3-methylimidazolium chloride, [1,3-DoMIM][Cl], with two isomers 1,2-dimethylimidazolium chloride, [1,2-MMIM][Cl], and 1,3-dimethylimidazolium chloride, [1,3-MMIM][Cl]. Molecular dynamics (MD) simulations, statistical mechanics calculations in the framework of the integral equation theory using one-dimensional (1D-) and three-dimensional (3D-) reference interaction site model (RISM) approaches as well as conductivity measurements were applied. From experiment and MD simulations it was found that the mobility/diffusion coefficients of cations in the limit of infinite dilution decrease with an increasing length of the cation alkyl chain, but not linearly. The aggregation tendency of cations with long alkyl chains at higher IL concentrations impedes their diffusivity. Binding free energies of imidazolium cations with the chloride anion estimated by RISM calculations, MD simulations and experiments reveal that the association of investigated ILs as model 1,1 electrolytes in water solutions is weak but evidently dependent on the molecular structure (alkyl chain length), which also strongly affects the mobility of cations.

How the spontaneous insertion of amphiphilic imidazolium-based cations changes biological membranes: a molecular simulation study.

The insertion of 1-octyl-3-methylimidazolium cations (OMIM(+)) from a diluted aqueous ionic liquid (IL) solution into a model of a bacterial cell membrane is investigated. Subsequently, the mutual interactions of cations inside the membrane and their combined effect on membrane properties are derived. The ionic liquid solution and the membrane model are simulated using molecular dynamics in combination with empirical force fields. A high propensity of OMIM(+) for membrane insertion is observed, with a cation concentration at equilibrium inside the membrane 47 times larger than in the solvent. Once inserted, cations exhibit a weak effective attraction inside the membrane at a distance of 1.3 nm. At this free energy minimum, negatively charged phosphates of the phospholipids are sandwiched between two OMIM(+) to form energetically favorable OMIM(+)-phosphate-OMIM(+) types of coordination. The cation-cation association free energy is 5.9 kJ mol(-1), whereas the activation barrier for dissociation is 10.1 kJ mol(-1). Subsequently, OMIM(+) are inserted into the leaflet of the membrane bilayer that represents the extracellular side. The cations are evenly distributed with mutual cation distances according to the found optimum distance of 1.3 nm. Because of the short length of the cation alkyl chains compared to lipid fatty acids, voids are generated in the hydrophobic core of the membrane. These voids disorder the fatty acids, because they enable fatty acids to curl into these empty spaces and also cause a thinning of the membrane by 0.6 nm. Additionally, the membrane density increases at its center. The presence of OMIM(+) in the membrane facilitates the permeation of small molecules such as ammonia through the membrane, which is chosen as a model case for small polar solutes. The permeability coefficient of the membrane with respect to ammonia increases substantially by a factor of seven. This increase is caused by a reduction of the involved free energy barriers, which is effected by the cations through the thinning of the membrane and favorable interactions of the delocalized OMIM(+) charge with ammonia inside the membrane. Overall, the results indicate the antimicrobial effect of amphiphilic imidazolium-based cations that are found in various common ILs. This effect is caused by an alteration of the permeability of the bacterial membrane and other property changes.

Transformations in plasma membranes of cancerous cells and resulting consequences for cation insertion studied with molecular dynamics.

Structural and energetic transformations in the plasma membrane of a cancerous cell are investigated together with related consequences for the insertion of small cationic compounds. Molecular dynamics simulations are performed with an empirical force field on two membrane models that represent the membrane of a cancerous cell (M-Cancer) and of a healthy cell (M-Eukar), respectively. An eight-fold increase of negatively charged phosphatidylserine in the external membrane layer as well as a reduction of cholesterol concentration by half is taken into account to describe the membrane transformation. Three additional reference membranes are prepared and consist of pure phosphatidylcholine (M-PC), where 20% is replaced with phosphatidylserine (M-PC0.8S0.2), and where 34% is replaced with cholesterol (M-PC0.66Ch0.34), respectively. Moreover, the free energy released by inserting octadecylmethylimidazolium (OMIM(+)), a cation found in a class of common ionic liquids, into M-Eukar, M-Cancer as well as into the three reference model membranes is derived by applying thermodynamic integration. We find that the presence of serine improves the solvation of the membrane through favorable electrostatic interactions with solvated sodium ions, where a significant number of sodium ions are capable of penetrating the upper polar layer of the membrane. However, the insertion free energy of OMIM(+) does not seem to be influenced by serine in the membrane. Furthermore, a significant serine induced structural reorganization of the membrane is not observed. In contrast, a reduction of cholesterol in the membrane models leads to smaller lipid surface densities, thinner membranes as well as less ordered and less stretched lipids as expected. We also observe that cholesterol reduction leads to a rougher membrane surface and an increased solvent accessibility of the hydrophobic membrane core. Membrane insertion of OMIM(+) becomes significantly more favorable in the absence of cholesterol, with an increased insertion free energy release of 7.1 kJ mol(-1) in M-Cancer compared to M-Eukar. Overall, the results suggest only a minor influence of serine on membrane organization but do not rule out an influence on cation insertion through a stronger cation adsorption to the membrane surface. In contrast, cholesterol seems to impede OMIM(+) insertion by increasing the density of polar lipids on the membrane surface and by flattening the membrane surface. These observations are shedding some light on the previously observed selective disruption of cancerous cells induced by cationic compounds such as found in ionic liquids.

Variation and decomposition of the partial molar volume of small gas molecules in different organic solvents derived from molecular dynamics simulations.

The partial molar volumes, V(i), of the gas solutes H2, CO, and CO2, solvated in acetone, methanol, heptane, and diethylether are determined computationally in the limit of infinite dilution and standard conditions. Solutions are described with molecular dynamics simulations in combination with the OPLS-aa force field for solvents and customized force field for solutes. V(i) is determined with the direct method, while the composition of V(i) is studied with Kirkwood-Buff integrals (KBIs). Subsequently, the amount of unoccupied space and size of pre-formed cavities in pure solvents is determined. Additionally, the shape of individual solvent cages is analyzed. Calculated V(i) deviate only 3.4 cm(3) mol(-1) (7.1%) from experimental literature values. Experimental V(i) variations across solutions are reproduced qualitatively and also quantitatively in most cases. The KBI analysis identifies differences in solute induced solvent reorganization in the immediate vicinity of H2 (<0.7 nm) and solvent reorganization up to the third solvation shell of CO and CO2 (<1.6 nm) as the origin of V(i) variations. In all solutions, larger V(i) are found in solvents that exhibit weak internal interactions, low cohesive energy density and large compressibility. Weak internal interactions facilitate solvent displacement by thermal solute movement, which enhances the size of solvent cages and thus V(i). Additionally, attractive electrostatic interactions of CO2 and the solvents, which do not depend on internal solvent interactions only, partially reversed the V(i) trends observed in H2 and CO solutions where electrostatic interactions with the solvents are absent. More empty space and larger pre-formed cavities are found in solvents with weak internal interactions, however, no evidence is found that solutes in any considered solvent are accommodated in pre-formed cavities. Individual solvent cages are found to be elongated in the negative direction of solute movement. This wake behind the moving solute is more pronounced in case of mobile H2 and in solvents with weaker internal interactions. However, deviations from a spherical solvent cage shape do not influence solute-solvent radial distribution functions after averaging over all solvent cage orientations and hence do not change V(i). Overall, the applied methodology reproduces V(i) and its variations reliably and the used V(i) decompositions identify the underlying reasons behind observed V(i) variations.

How ion properties determine the stability of a lipase enzyme in ionic liquids: a molecular dynamics study.

The influence of eight different ionic liquid (IL) solvents on the stability of the lipase Candida antarctica lipase B (CAL-B) is investigated with molecular dynamics (MD) simulations. Considered ILs contain cations that are based either on imidazolium or guanidinium as well as nitrate, tetrafluoroborate or hexafluorophosphate anions. Partial unfolding of CAL-B is observed during high-temperature MD simulations and related changes of CAL-B regarding its radius of gyration, surface area, secondary structure, amount of solvent close to the backbone and interaction strength with the ILs are evaluated. CAL-B stability is influenced primarily by anions in the order NO(3)(-)≪ BF(4)(-) < PF(6)(-) of increasing stability, which agrees with experiments. Cations influence protein stability less than anions but still substantially. Long decyl side chains, polar methoxy groups and guanidinium-based cations destabilize CAL-B more than short methyl groups, other non-polar groups and imidazolium-based cations, respectively. Two distinct causes for CAL-B destabilization are identified: a destabilization of the protein surface is facilitated mostly by strong Coulomb interactions of CAL-B with anions that exhibit a localized charge and strong polarization as well as with polar cation groups. Surface instability is characterized by an unraveling of α-helices and an increase of surface area, radius of gyration and protein-IL total interaction strength of CAL-B, all of which describe a destabilization of the folded protein state. On the other hand, a destabilization of the protein core is facilitated when direct core-IL interactions are feasible. This is the case when long alkyl chains are involved or when particularly hydrophobic ILs induce major conformational changes that enable ILs direct access to the protein core. This core instability is characterized by a disintegration of ß-sheets, diffusion of ions into CAL-B and increasing protein-IL van der Waals interactions. This process describes a stabilization of the unfolded protein state. Both of these processes reduce the folding free energy and thus destabilize CAL-B. The results of this work clarify the impact of ions on CAL-B stabilization. An extrapolation of the observed trends enables proposing novel ILs in which protein stability could be enhanced further.

On the different roles of anions and cations in the solvation of enzymes in ionic liquids.

The solvation of the enzyme Candida antarctica lipase B (CAL-B) was studied in eight different ionic liquids (ILs). The influence of enzyme-ion interactions on the solvation of CAL-B and the structure of the enzyme-IL interface are analyzed. CAL-B and ILs are described with molecular dynamics (MD) simulations in combination with an atomistic empirical force field. The considered cations are based on imidazolium or guanidinium that are paired with nitrate, tetrafluoroborate or hexafluorophosphate anions. The interactions of CAL-B with ILs are dominated by Coulomb interactions with anions, while the second largest contribution stems from van der Waals interactions with cations. The enzyme-ion interaction strength is determined by the ion size and the magnitude of the ion surface charge. The solvation of CAL-B in ILs is unfavorable compared to water because of large formation energies for the CAL-B solute cages in ILs. The internal energy in the IL and of CAL-B increases linearly with the enzyme-ion interaction strength. The average electrostatic potential on the surface of CAL-B is larger in ILs than in water, due to a weaker screening of charged enzyme residues. Ion densities increased moderately in the vicinity of charged residues and decreased close to non-polar residues. An aggregation of long alkyl chains close to non-polar regions and the active site entrance of CAL-B are observed in one IL that involved long non-polar decyl groups. In ILs that contain 1-butyl-3-methylimidazolium cations, the diffusion of one or two cations into the active site of CAL-B occurs during MD simulations. This suggests a possible obstruction of the active site in these ILs. Overall, the results indicate that small ions lead to a stronger electrostatic screening within the solvent and stronger interactions with the enzyme. Also a large ion surface charge, when more hydrophilic ions are used, increases enzyme-IL interactions. An increase of these interactions destabilizes the enzyme and impedes enzyme solvation due to an increase in solute cage formation energies.

Structural Dynamics and Catalytic Mechanism of ATP13A2 (PARK9) from Simulations.

ATP13A2 is a gene encoding a protein of the P5B subfamily of ATPases and is a PARK gene. Molecular defects of the gene are mainly associated with variations of Parkinson's disease (PD). Despite the established importance of the protein in regulating neuronal integrity, the three-dimensional structure of the protein currently remains unresolved crystallographically. We have modeled the structure and reactivity of the full-length protein in its E1-ATP state. Using molecular dynamics (MD), quantum cluster, and quantum mechanical/molecular mechanical (QM/MM) methods, we aimed at describing the main catalytic reaction, leading to the phosphorylation of Asp513. Our MD simulations suggest that two positively charged Mg2+ cations are present at the active site during the catalytic reaction, stabilizing a specific triphosphate binding mode. Using QM/MM calculations, we subsequently calculated the reaction profiles for the phosphoryl transfer step in the presence of one and two Mg2+ cations. The calculated barrier heights in both cases are found to be ∼12.5 and 7.5 kcal mol-1, respectively. We elucidated details of the catalytically competent ATP conformation and the binding mode of the second Mg2+ cofactor. We also examined the role of the conserved Arg686 and Lys859 catalytic residues. We observed that by significantly lowering the barrier height of the ATP cleavage reaction, Arg686 had major effect on the reaction. The removal of Arg686 increased the barrier height for the ATP cleavage by more than 5.0 kcal mol-1 while the removal of key electrostatic interactions created by Lys859 to the γ-phosphate and Asp513 destabilizes the reactant state. When missense mutations occur in close proximity to an active site residue, they can interfere with the barrier height of the reaction, which can halt the normal enzymatic rate of the protein. We also found large binding pockets in the full-length structure, including a transmembrane domain pocket, which is likely where the ATP13A2 cargo binds.

Magnesium-cationic dummy atom molecules enhance representation of DNA polymerase beta in molecular dynamics simulations: improved accuracy in studies of structural features and mutational effects.

Human DNA polymerase beta (pol beta) fills gaps in DNA as part of base excision DNA repair. Due to its small size it is a convenient model enzyme for other DNA polymerases. Its active site contains two Mg(2+) ions, of which one binds an incoming dNTP and one catalyzes its condensation with the DNA primer strand. Simulating such binuclear metalloenzymes accurately but computationally efficiently is a challenging task. Here, we present a magnesium-cationic dummy atom approach that can easily be implemented in molecular mechanical force fields such as the ENZYMIX or the AMBER force fields. All properties investigated here, namely, structure and energetics of both Michaelis complexes and transition state (TS) complexes were represented more accurately using the magnesium-cationic dummy atom model than using the traditional one-atom representation for Mg(2+) ions. The improved agreement between calculated free energies of binding of TS models to different pol beta variants and the experimentally determined activation free energies indicates that this model will be useful in studying mutational effects on catalytic efficiency and fidelity of DNA polymerases. The model should also have broad applicability to the modeling of other magnesium-containing proteins.

A model for self-diffusion of guanidinium-based ionic liquids: a molecular simulation study.

We propose a novel self-diffusion model for ionic liquids on an atomic level of detail. The model is derived from molecular dynamics simulations of guanidinium-based ionic liquids (GILs) as a model case. The simulations are based on an empirical molecular mechanical force field, which has been developed in our preceding work, and it relies on the charge distribution in the actual liquid. The simulated GILs consist of acyclic and cyclic cations that were paired with nitrate and perchlorate anions. Self-diffusion coefficients are calculated at different temperatures from which diffusive activation energies between 32-40 kJ/mol are derived. Vaporization enthalpies between 174-212 kJ/mol are calculated, and their strong connection with diffusive activation energies is demonstrated. An observed formation of cavities in GILs of up to 6.5% of the total volume does not facilitate self-diffusion. Instead, the diffusion of ions is found to be determined primarily by interactions with their immediate environment via electrostatic attraction between cation hydrogen and anion oxygen atoms. The calculated average time between single diffusive transitions varies between 58-107 ps and determines the speed of diffusion, in contrast to diffusive displacement distances, which were found to be similar in all simulated GILs. All simulations indicate that ions diffuse by using a brachiation type of movement: a diffusive transition is initiated by cleaving close contacts to a coordinated counterion, after which the ion diffuses only about 2 A until new close contacts are formed with another counterion in its vicinity. The proposed diffusion model links all calculated energetic and dynamic properties of GILs consistently and explains their molecular origin. The validity of the model is confirmed by providing an explanation for the variation of measured ratios of self-diffusion coefficients of cations and paired anions over a wide range of values, encompassing various ionic liquid classes as well as the simulated GILs. The proposed diffusion model facilitates the qualitative a priori prediction of the impact of ion modifications on the diffusive characteristics of new ionic liquids.

A force field for guanidinium-based ionic liquids that utilizes the electron charge distribution of the actual liquid: a molecular simulation study.

We propose a new all-atom force field for guanidinium-based ionic liquids (GILs) which is based on the charge distribution in the actual liquid. It comprises all cations that can be built by attaching alkyl chains of variable length to an acyclic or cyclic guanidinium compound and that are paired with nitrate or perchlorate anions. We based the parametrization of the force field on liquid-phase charge distributions to improve the prediction of energetic and dynamic properties of GILs. The impact of electron charge transfer and polarization on various properties of GILs is systematically assessed. A significant average electron charge transfer between -0.12 and -0.06 e from anions to the central guanidinium group of the cations and a strong polarization of acyclic cations are observed by applying a combined quantum mechanical/molecular mechanical (QM/MM) approach. Molecular dynamics simulations of GILs are performed, utilizing the proposed force field. Derived structures approach the accuracy of QM/MM structures, and a previously reported crystal structure remains stable throughout the simulations. Mass densities are reproduced with a deviation of only 1.4% from experimental data. The calculated melting point of a GIL crystal deviates only 8% from the measured value. Self-diffusion coefficients of various GILs are reported, and a comparison with a diffusion coefficient derived from experimental data indicates that the values are within a reasonable range. We observe that the melting point of a GIL crystal was lowered up to 60 K and that diffusion coefficients are substantially increased by a factor of up to 3.5 upon consideration of charge transfer and polarization. The results demonstrate that liquid-phase partial charges are capable of improving the quality of ionic liquid force fields substantially and that their utilization led to a model that can be applied to predict structural, energetic, and dynamic properties of GILs.

On the Stability of Proteins Solvated in Imidazolium-Based Ionic Liquids Studied with Replica Exchange Molecular Dynamics.

The stability of two small proteins, one composed of three α-helices (α-peptide) and another composed of a ß-sheet (ß-peptide) solvated in five different ionic liquids (ILs), is analyzed using replica exchange molecular dynamics (REMD) simulations. ILs are composed of 1-butyl-3-methylimidazolium (BMIM) cations, paired with five different anions of varying hydrophilicity and size, namely, Cl-, NO3-, BF4-, PF6-, and NTf2-. REMD simulations greatly improve structure sampling and mitigate bias toward the initial folded peptide structure, thereby providing more adequate simulations to study protein stability. Cluster analysis, DSSP analysis and derivation of radius of gyration, interaction energies, and hydrogen bonding are used to quantify structural peptide changes in a large temperature range from 250 to 650 K. α-Peptides are least stable in ILs that contain small anions with localized negative charge, such as in BMIM-Cl and BMIM-NO3. Destabilization is caused by direct electrostatic interactions of anions with α-helices that are exposed to the solvent. This destabilization is characterized not by unfolded but instead by compact misfolded structures. Also, ß-peptides retain compact structures up to at least 400 K, below which unfolding hardly occurs. However, intrapeptide hydrogen bonds that constitute the ß-sheet are not exposed to the solvent. Therefore, ß-peptides are generally more stable than α-peptides in all considered ILs. Moreover, on contrary to α-peptides, ß-peptides are least stable in less polar ILs, such as BMIM-PF6 and BMIM-NTf2, because dissolving ß-sheets requires large structural changes of the peptide. Such transitions are energetically less opposed in ILs with weaker mutual ion coordination. A large interaction density within ILs, for example, in BMIM-Cl, is thus kinetically trapping ß-peptides in the original folded state. Additionally, in BMIM-BF4, interactions with ß-peptides are so weak, compared to an aqueous solvent, resulting in stronger interactions within the peptide, which extend ß-sheets, hence causing misfolding of a different kind. The results reveal how direct ion-peptide interactions and solvent reorganization energy in ILs are both crucial in determining protein stability. These insights could translate into guidelines for the design of new IL solvents with improved protein stability.

Towards accurate ab initio QM/MM calculations of free-energy profiles of enzymatic reactions.

Reliable studies of enzymatic reactions by combined quantum mechanical/molecular mechanics (QM/MM) approaches, with an ab initio description of the quantum region, presents a major challenge to computational chemists. The main problem is the need for a large amount of computer time to evaluate the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. This work presents major progress toward the evaluation of ab initio QM/MM free-energy surfaces and activation free energies of reactions in enzymes and in solutions. This is done by exploiting our previous idea of using the empirical valence bond (EVB) method as a reference potential and then using the linear response approximation (LRA) approach to evaluate the free energies of transfer from the EVB to the QM/MM surfaces in the reactant and product state. However, the new crucial step involves the use of a constraint at the transition state that fixes the system at a given value of the reaction coordinate and allows us to use the LRA at the transition state. The advance offered by the present approach is particularly significant because it evaluates the free energy associated with both the substrate and the solvent motions. This evaluation appeared to be a relatively simple task once one uses a classical reference potential. The main problem has been using the reference potential for the evaluation of the free-energy contributions associated with the solute motions where the difference between the reference EVB potential and the QM/MM potential can be large. The present refinement finally allows us to overcome the problems with the solute fluctuations and therefore to obtain, for the first time, a free-energy barrier that reflects the solute entropy properly. Thus, we present a way to evaluate the complete QM/MM activation free energy with an equal footing treatment of the solute and the solvent. This provides a general consistent and effective strategy for evaluating the QM/MM activation free energies in proteins and in solution. Our advance allows one to explore consistently various mechanistic and catalytic proposals while using ab initio (ai) QM/MM approaches.

On possible pitfalls in ab initio quantum mechanics/molecular mechanics minimization approaches for studies of enzymatic reactions.

Reliable studies of enzymatic reactions by combined quantum mechanics/molecular mechanics (QM/MM) approaches, with an ab initio description of the quantum region, presents a major challenge to computational chemists. The main problem is the need for a very large computer time for the evaluation of the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. A seemingly reasonable alternative is to perform energy minimization studies of the type used in gas-phase ab initio studies. However, it is hard to see why such an approach should give reliable results in protein active sites. To examine the problems with energy minimization QM/MM approaches, we chose the hypothetical reaction of a metaphosphate ion with water in the Ras.GAP complex. This hypothetical reaction served as a simple benchmark reaction. The possible problems with the QM/MM minimization were explored by generating several protein configurations from long MD simulations and using energy minimization and scanning of the reaction coordinates to evaluate the corresponding potential energy surfaces of the reaction for each of these different protein configurations. Comparing these potential energy surfaces, we found major variations of the corresponding minima. Furthermore, the reaction energies and activation energies also varied significantly even for similar protein configurations. The specific coordination of a magnesium ion, present in the active center of the protein complex, turned out to influence the energetics of the reaction in a major way, where a direct coordination to the reactant leads to an increase of the activation energy by 17 kcal/mol. Apparently, using energy minimization to generate potential surfaces for an enzymatic reaction, while starting from a single protein structure, could lead to major errors in calculations of activation free energies and binding free energies. Thus we believe that extensive samplings of the configurational space of the protein are essential for meaningful determination of the energetics of enzymatic reactions. The possible relevance of our conclusion with regard to a recent study of the RasGAP reaction is discussed.

What Determines CO2 Solubility in Ionic Liquids? A Molecular Simulation Study.

Molecular dynamics (MD) simulations of 10 different pure and CO2-saturated ionic liquids are performed to identify the factors that determine CO2 solubility. Imidazolium-based cations with varying alkyl chain length and functionalization are paired with anions of different hydrophobicity and size. Simulations are carried out with an empirical force field based on liquid-phase charges. The partial molar volume of CO2 in ionic liquids (ILs) varies from 30 to 40 cm(3)/mol. This indicates that slight ion displacements are necessary to enable CO2 insertions. However, the absorption of CO2 does not affect the overall organization of ions in the ILs as demonstrated by almost equal cation-anion radial distribution functions of pure ILs and ILs saturated with CO2. The solubility of CO2 in ILs is not influenced by direct CO2-ion interactions. Instead, a strong correlation between the ratio of unoccupied space in pure ILs and their ability to absorb CO2 is found. This preformed unoccupied space is regularly dispersed throughout the ILs and needs to be expanded by slight ion displacements to accommodate CO2. The amount of preformed unoccupied space is a good indicator for ion cohesion in ILs. Weak electrostatic cation-anion interaction densities in ILs, i.e., weak ion cohesion, leads to larger average distances between ions and hence to more unoccupied space. Weak ion cohesion facilitates ion displacement to enable an expansion of empty space to accommodate CO2. Moreover, it is demonstrated that the strength of ion cohesion is primarily determined by the ion density, which in turn is given by the ion sizes. Ion cohesion is influenced additionally to a smaller extent by local electrostatic interactions among ion moieties between which CO2 is inserted and which do not depend on the ion density. Overall, the factors that determine the solubility of CO2 in ILs are identified consistently across a large variety of constituting ions through MD simulations.

Protein-like dynamics of polycarbonate polymers in water.

The dynamics of amphiphilic peptide-mimicking polycarbonate polymers are investigated, considering variations in polymer length, monomer sequence, and monomer modification. The polymers are simulated in aqueous solution with atomistic molecular dynamics simulations and an empirical force field. Various structural polymer properties, interaction strengths, and solvation free energies are derived. It is found that water is a less favorable solvent for these polymers than for peptides. Moreover, polymers readily adopt irreversibly a compact state that consists of a variety of distinct compact conformations that are adopted through frequent transitions. Furthermore, the polymers exhibit a strong propensity to form large aggregates. The driving forces for these processes appear to be a hydrophobic effect and more favorable polymer-solvent interactions of aggregates that overcome the otherwise strong mutual repulsion between the positively charged polymers. Replacing hydrophobic residues with polar side chains destabilizes the compact conformations of the polymers. Our results also indicate that the monomer sequence has little effect on the overall solvation properties of the polymer molecule. However, the sequence influences flexibility and compactness of the monomer in solution. Overall, the results of this work confirm the protein-like characteristics of these polymers and elucidate the role of single residues in influencing the structure and aggregation in aqueous solution.

Impact of ionic liquids in aqueous solution on bacterial plasma membranes studied with molecular dynamics simulations.

The impact of five different imidazolium-based ionic liquids (ILs) diluted in water on the properties of a bacterial plasma membrane is investigated using molecular dynamics (MD) simulations. Cations considered are 1-octyl-3-methylimidazolium (OMIM), 1-octyloxymethyl-3-methylimidazolium (OXMIM), and 1-tetradecyl-3-methylimidazolium (TDMIM), as well as the anions chloride and lactate. The atomistic model of the membrane bilayer is designed to reproduce the lipid composition of the plasma membrane of Gram-negative Escherichia coli. Spontaneous insertion of cations into the membrane is observed in all ILs. Substantially more insertions of OMIM than of OXMIM occur and the presence of chloride reduces cation insertions compared to lactate. In contrast, anions do not adsorb onto the membrane surface nor diffuse into the bilayer. Once inserted, cations are oriented in parallel to membrane lipids with cation alkyl tails embedded into the hydrophobic membrane core, while the imidazolium-ring remains mostly exposed to the solvent. Such inserted cations are strongly associated with one to two phospholipids in the membrane. The overall order of lipids decreased after OMIM and OXMIM insertions, while on the contrary the order of lipids in the vicinity of TDMIM increased. The short alkyl tails of OMIM and OXMIM generate voids in the bilayer that are filled by curling lipids. This cation induced lipid disorder also reduces the average membrane thickness. This effect is not observed after TDMIM insertions due to the similar length of cation alkyl chain and the fatty acids of the lipids. This lipid-mimicking behavior of inserted TDMIM indicates a high membrane affinity of this cation that could lead to an enhanced accumulation of cations in the membrane over time. Overall, the simulations reveal how cations are inserted into the bacterial membrane and how such insertions change its properties. Moreover, the different roles of cations and anions are highlighted and the fundamental importance of cation alkyl chain length and its functionalization is demonstrated.

Extraction of tryptophan with ionic liquids studied with molecular dynamics simulations.

Extraction of amino acids from aqueous solutions with ionic liquids (ILs) in biphasic systems is analyzed with molecular dynamics (MD) simulations. Extraction of tryptophan (TRP) with the imidazolium-based ILs [C(4)mim][PF(6)], [C(8)mim][PF(6)], and [C(8)mim][BF(4)] are considered as model cases. Solvation free energies of TRP are calculated with MD simulations and thermodynamic integration in combination with an empirical force field, whose parametrization is based on the liquid-phase charge distribution of the ILs. Calculated solvation free energies reproduce successfully all observed experimental trends according to the previously reported partition of TRP between water and IL phases. Water is present in ILs as a cosolvent, due to direct contact with the aqueous phase during extraction, and is found to play a major role in the extraction of TRP. Water improves solvation of cationic TRP by 7.8 and 5.1 kcal/mol in [C(4)mim][PF(6)] and [C(8)mim][PF(6)], respectively, which is in the case of [C(4)mim][PF(6)] sufficient to extract TRP. Extraction in [C(8)mim][PF(6)] is not feasible, since the hydrophobic octyl groups of the cations limit the water concentration in the IL. The solvation of cationic TRP is 2.4 kcal/mol less favorable in [C(8)mim][PF(6)] than in [C(4)mim][PF(6)]. Water improves the solvation of TRP in ILs mostly through dipole-dipole interactions with the polar backbone of TRP. Extraction is most efficient with [C(8)mim][BF(4)], where hydrophilic BF(4)(-) anions substantially increase the water concentration in the IL. Additionally, stronger direct electrostatic interactions of TRP with BF(4)(-) anions improve its solvation in the IL further. The solvation of cationic TRP in [C(8)mim][BF(4)] is 3.4 kcal/mol more favorable than in [C(8)mim][PF(6)]. Overall, the extractive power of the ILs correlates with the water saturation concentration of the IL phase, which in turn is determined by the hydrophilicity of the constituting ions. The results of this work identify relations between the extraction performance of ILs and the basic chemical properties of the ions, which provide guidelines that could contribute to the design of improved novel ILs for amino acid extraction.

Proton transfer between tryptophan and ionic liquid solvents studied with molecular dynamics simulations.

The reaction free energies and associated pK(a) values for proton transfer from positively charged tryptophan (HTrp(+)) to the two pure ionic liquids (ILs) BMIM-PF6 and BMIM-BF4 are derived from molecular simulations. IL solvation effects are examined with molecular dynamics simulations together with an empirical force field in which the average charge distribution in the actual IL is taken into account. A combination of molecular mechanical and quantum mechanical description (QM/MM) is used to examine the protonation of the anion constituents of the ILs. A dissociation of the protonated anions is observed into hydrogen fluoride and BF3 or PF5. Finally, pK(a) values of 16.5 and 21.5 in BMIM-BF4 and BMIM-PF6, respectively, are found for proton transfer from HTrp(+) to PF6(-) and BF4(-) anions, which indicates that a deprotonation of HTrp(+) is highly unfavorable compared to aqueous solutions. An examination of the contributions to the reaction free energies demonstrates that a deprotonation of tryptophan is impeded because two ions need to be annihilated for the reaction to occur: HTrp(+) and an anion. While the solvation effects induced by the two ILs are similar, the low proton acceptance of PF6(-) anions leads to the larger pK(a) value in BMIM-PF6. Also, estimates suggest that IL-induced pK(a) shifts are comparably small in proton transfer reactions where the total number of ions remains unchanged. For the first time, pK(a) values of acids were determined computationally in ILs. The obtained results elucidate the role of solvation effects on proton transfer between amino acids and ILs and improve our understanding of the observed pH memory of proteins that are solvated in ILs.

What determines the miscibility of ionic liquids with water? Identification of the underlying factors to enable a straightforward prediction.

Whether an ionic liquid (IL) is water-miscible or immiscible depends on the particular ions that constitute it. We propose an explanation, based on molecular simulations, how ions determine the miscibility of ILs and suggest a straightforward and computationally inexpensive method to predict the miscibility of arbitrary new ILs. The influence of ions on the solvation of water is analyzed by comparing molecular dynamics simulations of water in 9 different ILs with varying cation and anion constituents. The solvation of water in ILs is found to depend primarily on the electrostatic water-ion interaction strength, which, in turn, is determined mainly by two factors: primarily, by the size of the ions and secondarily by the amount of charge on the ion surface that is coordinated with water. It is demonstrated that large ions lead to weaker interactions with water, due to the involved delocalization of the ion charge. A large charge on the ion surface, which is determined by the chemical structure of the ion, strengthens water-ion interactions. We observe that whenever the interaction strength of water with ions exceeds a certain threshold, an IL becomes water-miscible. On the basis of these findings, a simple equation is derived that estimates the water-ion interaction strength. With this equation it is possible to predict most of the observed water-miscibilities of a sample of 83 ILs correctly. A linear increase of the water saturation concentration with the estimated water-ion interaction strength is observed in water-immiscible ILs, which can be utilized to predict the water concentration in new ILs.