Probing the transition from hydrophilic to hydrophobic solvation with atomic scale resolution

Picosecond and femtosecond X-ray absorption spectroscopy is used to probe the changes of the solvent shell structure upon electron abstraction of aqueous iodide using an ultrashort laser pulse. The transient L(1,3) edge EXAFS at 50 ps time delay points to the formation of an expanded water cavity around the iodine atom, in good agreement with classical and quantum mechanical/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. These also show that while the hydrogen atoms pointed toward iodide, they predominantly point toward the bulk solvent in the case of iodine, suggesting a hydrophobic behavior. This is further confirmed by quantum chemical (QC) calculations of I(-)/I(0)(H(2)O)(n=1-4) clusters. The L(1) edge sub-picosecond spectra point to the existence of a transient species that is not present at 50 ps. The QC calculations and the QM/MM MD simulations identify this transient species as an I(0)(OH(2)) complex inside the cavity. The simulations show that upon electron abstraction most of the water molecules move away from iodine, while one comes closer to form the complex that lives for 3-4 ps. This time is governed by the reorganization of the main solvation shell, basically the time it takes for the water molecules to reform an H-bond network. Only then is the interaction with the solvation shell strong enough to pull the water molecule of the complex toward the bulk solvent. Overall, much of the behavior at early times is determined by the reorientational dynamics of water molecules and the formation of a complete network of hydrogen bonded molecules in the first solvation shell.     

Dynamics of proton transfer and vibrational relaxation in dilute hydrofluoric acid

The molecular mechanisms in both vibrational relaxation and proton transfer (PT) associated with infrared (IR)-induced PT in a dilute hydrofluoric acid solution at ambient temperature are studied by molecular dynamics (MD) simulations with the multistate empirical valence bond model. To investigate the solvation dynamics, a collective solvent coordinate and its perpendicular bath modes are defined from the diabatic energy gap and their motions are examined by the generalized Langevin equation (GLE) formalism. The GLE analysis using the equilibrium MD simulation shows that the major solvent reorganizations in the PT are represented by the libration and hindered translation. In particular, the libration gives the stronger coupling to the solvent reorganization and the faster relaxation. The nonequilibrium MD simulation demonstrated that both the HF stretching vibration and the solvent reorganization relax on a similar time scale and thus compete in the PT. It also supported the "presolvation mechanism" for the PT in this system.     

Molecular dynamics with the united-residue model of polypeptide chains. II. Langevin and Berendsen-bath dynamics and tests on model alpha-helical systems

The implementation of molecular dynamics (MD) with our physics-based protein united-residue (UNRES) force field, described in the accompanying paper, was extended to Langevin dynamics. The equations of motion are integrated by using a simplified stochastic velocity Verlet algorithm. To compare the results to those with all-atom simulations with implicit solvent in which no explicit stochastic and friction forces are present, we alternatively introduced the Berendsen thermostat. Test simulations on the Ala(10) polypeptide demonstrated that the average kinetic energy is stable with about a 5 fs time step. To determine the correspondence between the UNRES time step and the time step of all-atom molecular dynamics, all-atom simulations with the AMBER 99 force field and explicit solvent and also with implicit solvent taken into account within the framework of the generalized Born/surface area (GBSA) model were carried out on the unblocked Ala(10) polypeptide. We found that the UNRES time scale is 4 times longer than that of all-atom MD simulations because the degrees of freedom corresponding to the fastest motions in UNRES are averaged out. When the reduction of the computational cost for evaluation of the UNRES energy function is also taken into account, UNRES (with hydration included implicitly in the side chain-side chain interaction potential) offers about at least a 4000-fold speed up of computations relative to all-atom simulations with explicit solvent and at least a 65-fold speed up relative to all-atom simulations with implicit solvent. To carry out an initial full-blown test of the UNRES/MD approach, we ran Berendsen-bath and Langevin dynamics simulations of the 46-residue B-domain of staphylococcal protein A. We were able to determine the folding temperature at which all trajectories converged to nativelike structures with both approaches. For comparison, we carried out ab initio folding simulations of this protein at the AMBER 99/GBSA level. The average CPU time for folding protein A by UNRES molecular dynamics was 30 min with a single Alpha processor, compared to about 152 h for all-atom simulations with implicit solvent. It can be concluded that the UNRES/MD approach will enable us to carry out microsecond and, possibly, millisecond simulations of protein folding and, consequently, of the folding process of proteins in real time.     

Coarse-grained molecular dynamics simulations of the sphere to rod transition in surfactant micelles

Surfactant molecules self-assemble in aqueous solutions to form various micellar structures such as spheres, rods, or lamellae. Although phase transitions in surfactant solutions have been studied experimentally, their molecular mechanisms are still not well understood. In this work, we show that molecular dynamics (MD) simulations using the coarse-grained (CG) MARTINI force field and explicit CG solvent, validated against atomistic MD studies, can accurately represent micellar assemblies of cetyltrimethylammonium chloride (CTAC). The effect of salt on micellar structures is studied for aromatic anionic salts, e.g., sodium salicylate, and simple inorganic salts, e.g., sodium chloride. Above a threshold concentration, sodium salicylate induces a sphere to rod transition in the micelle. CG MD simulations are shown to capture the dynamics of this shape transition and support a mechanism based on the reduction in the micelle-water interfacial tension induced by the adsorption of the amphiphilic salicylate ions. At the threshold salt concentration, the interface is nearly saturated with adsorbed salicylate ions. Predictions of the effect of salt on the micelle structure in different CG solvent models, namely, single-site standard water and three-site polarizable water, show qualitative agreement. This suggests that phase transitions in aqueous micelle solutions could be investigated by using standard CG water models which allow for 3 orders of magnitude reduction in the computational time as compared to that required for atomistic MD simulations.     

Redox entropy of plastocyanin: developing a microscopic view of mesoscopic polar solvation

We report applications of analytical formalisms and molecular dynamics (MD) simulations to the calculation of redox entropy of plastocyanin metalloprotein in aqueous solution. The goal of our analysis is to establish critical components of the theory required to describe polar solvation at the mesoscopic scale. The analytical techniques include a microscopic formalism based on structure factors of the solvent dipolar orientations and density and continuum dielectric theories. The microscopic theory employs the atomistic structure of the protein with force-field atomic charges and solvent structure factors obtained from separate MD simulations of the homogeneous solvent. The MD simulations provide linear response solvation free energies and reorganization energies of electron transfer in the temperature range of 280-310 K. We found that continuum models universally underestimate solvation entropies, and a more favorable agreement is reported between the microscopic calculations and MD simulations. The analysis of simulations also suggests that difficulties of extending standard formalisms to protein solvation are related to the inhomogeneous structure of the solvation shell at the protein-water interface combining islands of highly structured water around ionized residues along with partial dewetting of hydrophobic patches. Quantitative theories of electrostatic protein hydration need to incorporate realistic density profile of water at the protein-water interface.     

Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.     

Adsorption of acridine orange at a C8,18/water/acetonitrile interface

Fully atomistic simulations are used to characterize the molecular dynamics (MD) of acridine orange (3,6-dimethylaminoacridine) at a chromatographic interface. Multiple 1 ns MD simulations were performed for acridine orange at the interface between three different acetonitrile/water mixtures (0/100, 20/80, and 50/50) with C8 and C18 alkyl chains. The diffusion coefficient, D, of acridine orange in pure solvent was found to be 4 times smaller at the water/C18 interface (D = 0.022 x 10(-4) cm2/s) than in bulk water (D = 0.087 x 10(-4) cm2/s), in qualitative agreement with experiment. Rotational reorientation times were 20 and 700 ps, which also agree favorably with the measured time scales of 130 and 740 ps. Contrary to experiment, the simulations found that for increasing surface coverage, the diffusion coefficient for acridine decreased. Detailed analysis of the solvent structure showed that the transport properties of acridine were primarily governed by the solvent distribution above the functionalized surface. The solvent structure, in turn, was largely determined by the surface consisting of the silica layer, the alkyl chains, and their functionalization.     

Molecular dynamics of ��-CD in water/co-solvent mixtures

Molecular dynamics (MD) simulations on ??-cyclodextrin (??-CD) in water, ethanol (EtOH), methanol (MeOH) and mixtures of these solvents have been carried out at 300 K over a time period of 15 ns using the AMBER force field. The hydrated X-ray crystallographic structure has four water molecules inside the cavity, defined by a more precise boundary for the ??-CD cavity. From the simulations, 2?C4 encapsulated water molecules are most probably found. In an ethanol co-solvent system, the ??-CD cavity is occupied with one ethanol molecule located in two discrete sites: below and above the O4(n) plane, which is in agreement with experimental results. In all systems, the average values of tilt angles of the obtained structures are higher than the tilt angles of the X-ray structures. The investigations of the alcohol orientations in co-solvent mixtures reveal the hydrophobic environment of the cavity and the hydrophilic atmosphere at both rims of ??-CD.     

Incorporating the excluded solvent volume and surface charges for computing solvation free energy

Gauss's law or Poisson's equation is conventionally used to calculate solvation free energy. However, the near‐solute dielectric polarization from Gauss's law or Poisson's equation differs from that obtained from molecular dynamics (MD) simulations. To mimic the near‐solute dielectric polarization from MD simulations, the first‐shell water was treated as two layers of surface charges, the densities of which are proportional to the electric field at the solvent molecule that is modeled as a hard sphere. The intermediate water was treated as a bulk solvent. An equation describing the solvation free energy of ions using this solvent scheme was derived using the TIP3P water model. © 2013 Wiley Periodicals, Inc.     

The influence of aqueous versus glassy solvents on protein dynamics: vibrational echo experiments and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo measurements are used to measure the vibrational dephasing of the CO stretching mode of carbonmonoxy-hemoglobin (HbCO), a myoglobin mutant (H64V), and a bacterial cytochrome c(552) mutant (Ht-M61A) in aqueous solution and trehalose glasses. The vibrational dephasing of the heme-bound CO is significantly slower for all three proteins embedded in trehalose glasses compared to that of aqueous protein solutions. All three proteins exhibit persistent but notably slower spectral diffusion when the protein surface is fixed by the glassy solvent. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to reveal that the structural dynamics, as sensed by the CO, of the three proteins in trehalose and aqueous solution are dominated by fast (tens of femtoseconds), motionally narrowed fluctuations. MD simulations of H64V in dynamic and "static" water are presented as models of the aqueous and glassy environments. FFCFs are calculated from the H64V simulations and qualitatively reproduce the important features of the experimentally extracted FFCFs. The suppression of long time scale (picoseconds to tens of picoseconds) frequency fluctuations (spectral diffusion) in the glassy solvent is the result of a damping of atomic displacements throughout the protein structure and is not limited to structural dynamics that occur only at the protein surface. The analysis provides evidence that some dynamics are coupled to the hydration shell of water, supporting the idea that the bioprotection offered by trehalose is due to its ability to immobilize the protein surface through a thin, constrained layer of water.     

Parametrization, molecular dynamics simulation, and calculation of electron spin resonance spectra of a nitroxide spin label on a polyalanine alpha-helix

The nitroxide spin label 1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl-methanethiosulfonate (MTSSL), commonly used in site-directed spin labeling of proteins, is studied with molecular dynamics (MD) simulations. After developing force field parameters for the nitroxide moiety and the spin label linker, we simulate MTSSL attached to a polyalanine alpha-helix in explicit solvent to elucidate the factors affecting its conformational dynamics. Electron spin resonance spectra at 9 and 250 GHz are simulated in the time domain using the MD trajectories and including global rotational diffusion appropriate for the tumbling of T4 Lysozyme in solution. Analysis of the MD simulations reveals the presence of significant hydrophobic interactions of the spin label with the alanine side chains.     

Double-stranded DNA dissociates into single strands when dragged into a poor solvent

DNA displays a richness of biologically relevant supramolecular structures, which depend on both sequence and ambient conditions. The effect of dragging double-stranded DNA (dsDNA) from water into poor solvent on the double-stranded structure is still unclear because of condensation. Here, we employed single molecule techniques based on atomic force microscopy and molecular dynamics (MD) simulations to investigate the change in structure and mechanics of DNA during the ambient change. We found that the two strands are split apart when the dsDNA is pulled at one strand from water into a poor solvent. The findings were corroborated by MD simulations where dsDNA was dragged from water into poor solvent, revealing details of the strand separation at the water/poor solvent interface. Because the structure of DNA is of high polarity, all poor solvents show a relatively low polarity. We speculate that the principle of spontaneous unwinding/splitting of dsDNA by providing a low-polarity (in other word, hydrophobic) micro-environment is exploited as one of the catalysis mechanisms of helicases.     

Improvement of dissolution properties of a new Helicobacter pylori eradicating agent (TG44) by inclusion complexation with beta-cyclodextrin

The interaction of a newly developed Helicobacter pylori eradicating agent (TG44, 4-methylbenzyl-4'-[trans-4-(guanidinomethyl)cyclohexylcarbonyloxy]biphenyl-4-carboxlylate monohydrochloride) with beta-cyclodextrin (beta-CyD) in aqueous solution and in solid state was studied to gain insight into the high in-vivo H. pylori eradicating activity of TG44/beta-CyD complex. The interaction was studied by the solubility method, spectroscopic methods, powder X-ray diffractometry and differential scanning colorimetry (DSC). TG44 gave A(L)-type phase solubility diagram with beta-CyD in water, showing a linear increase in solubility of the drug up to 8 mM beta-CyD concentration. The solubility of TG44 (0.04 mM in water at 25 degrees C) increased about 70-folds at 8 mM beta-CyD. Ultraviolet, circular dichroism, fluorescence and (1)H-nuclear magnetic resonance spectroscopic studies indicated that TG44 forms the inclusion complex with beta-CyD in a 1:1 stoichiometry and the biphenyl moiety of TG44 is preferably included in the beta-CyD cavity in water. The Giordano plot made by monitoring changes in the fusion enthalpy of TG44 (about 184 degrees C) suggested that TG44 forms the 1:1 complex with beta-CyD in the solid state. The TG44/beta-CyD solid complex in a 1:1 stoichiometry was prepared by the grinding and spray-drying methods and confirmed by powder X-ray diffractometry and DSC that the complex is in an amorphous state. The initial dissolution rate of TG44/beta-CyD complex was significantly faster than those of the drug alone and the physical mixture of both components, maintaining higher supersaturated concentrations of the drug for a long time. The results suggested that the higher eradicating activity of TG44/beta-CyD complex to Helicobacter pylori, compared with that of the drug alone, is attributable at least partly to the faster dissolving property of the complex and its ability to maintain the supersaturated state of the drug in the gastric fluid.     

Self-assembling cyclic peptides: molecular dynamics studies of dimers in polar and nonpolar solvents

The self-assembly of cyclic D,L-alpha-peptides into hollow nanotubes is a crucial mechanistic step in their application as antibacterial and drug-delivery agents. To understand this process, molecular dynamics (MD) simulations were performed on dimers of cyclic peptides formed from cyclo [(-L-Trp-D-N-MeLeu-)4-]2 and cyclo [(-L-Trp-D-Leu-)4-]2 subunits in nonpolar (nonane) and polar (water) solvent. The dimers were observed to be stable only in nonpolar solvent over the full 10 ns length of the MD trajectory. The behavior of the dimers in different solvents is rationalized in terms of the intersubunit hydrogen bonding, hydrogen bonding with the solvent, and planarity of the rings. It is shown that the phi and psi dihedral angles of a single uncapped ring in nonane lie in the beta-sheet region of the Ramachandran plot, and the ring stays in a flat conformation. Steered MD (SMD) simulations based on Jarzynski's equality were performed to obtain the potential of mean force as a function of the distance between the two rings of the capped dimer in nonane. It is also shown that a single peptide subunit prefers to reside close to the nonane/water interface rather than in bulk solvent because of the amphiphilic character of the peptide ring. The present MD results build the foundation for using MD simulations to study the mechanism of the formation of cyclic peptide nanotubes in lipid bilayers.     

Refinement of the NMR structures of alpha-conotoxin MI using molecular dynamics simulation with explicit solvent water and a full molecular force field

Three NMR structures of alpha-conotoxin MI, a potent antagonist of the nicotinic acetylcholine receptor, have been refined using molecular dynamics (MD) simulation with explicit water. Although the convergence of the NMR structures of alpha-conotoxin MI was not sufficient to provide detailed structural features, the average structures obtained from MD simulations converged to one conformation, providing structural characteristics. The resulting structure was also found to be consistent with the results of amide proton-exchange experiments. These results demonstrate that MD simulation with explicit solvent water is very useful in refining NMR structures.     

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes in polar and apolar solutions

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes nTs (n = 1 to 10) are investigated in solvents of n-hexane, 1,4-dioxane, carbon tetrachloride, chloroform, and water by using density functional theory (DFT) and molecular dynamics (MD) simulations. Both implicit and explicit solvent models are employed. The polarized continuum model (PCM) calculations and MD simulations demonstrate the weak solvent effects on the electronic structures of alpha-oligothiophenes. The lowest dipole-allowed vertical excitation energies of nTs, obtained from time-dependent DFT/PCM calculations at the B3LYP/6-31G(d) level, exhibit a red shift as the solvent polarity increases, in agreement with experiments. The studied solvents have little impact on the state order of the low-lying excited states provided that the nTs are kept in C2h or C2v symmetry. The MD simulations demonstrate that the chain conformations are distorted to some extent in polar and nonpolar solvents. A qualitative picture of the distribution of solvent molecules around the solvated nTs is drawn by means of radial and spatial distribution functions. The S...H-O and pi...H-O solute-solvent interactions are insignificant in aqueous solution.     

Impact of the solvent on the conformational isomerism of calix[4]arenes: a study based on continuum solvation models

The influence of solvation on the conformational isomerism of calix[4]arene and p-tert-butylcalix[4]arene has been investigated by using the continuum model reported by Miertus, Scrocco, and Tomasi (MST). The quantum mechanical (QM) and semiclassical (SC) formalisms of the MST model have been considered for two different solvents (chloroform and water). The suitability of the QM-MST and SC-MST methods has been examined by comparison with previous results derived from classical molecular dynamics (MD) simulations with explicit solvent molecules. The application of the continuum model to the solute configurations generated by using in vacuo classical MD simulations provides a fast strategy to evaluate the effects of the solvent on the conformational preferences of calixarenes. These encouraging results allow us to propose the use of continuum models to solutes with complex molecular structures, which are traditionally studied by MD simulations.