Viewpoint 9--molecular structure of aqueous interfaces

Chemical Reaction Dynamics at Liquid Interfaces: A Computational Approach

Recent advances in experimental and theoretical studies of liquid interfaces provide remarkable evidence for the unique properties of these systems. In this review we examine how these properties affect the thermodynamics and kinetics of chemical reactions which take place at the liquid/vapor interface and at the liquid/liquid interface. We demonstrate how the rapidly varying density and viscosity, the marked changes in polarity and the surface roughness manifest themselves in isomerization, electron transfer and photodissociation reactions.

Miscibility at the immiscible liquid/liquid interface: A molecular dynamics study of thermodynamics and mechanism

Molecular dynamics simulations are used to study the dissolution of water into an adjacent, immiscible organic liquid phase. Equilibrium thermodynamic and structural properties are calculated during the transfer of water molecule(s) across the interface using umbrella sampling. The net free energy of transfer agrees reasonably well with experimental solubility values. We find that water molecules "prefer" to transfer into the adjacent phase one-at-a-time, without co-transfer of the hydration shell, as in the case of evaporation. To study the dynamics and mechanism of transfer of water to liquid nitrobenzene, we collected over 400 independent dissolution events. Analysis of these trajectories suggests that the transfer of water is facilitated by interfacial protrusions of the water phase into the organic phase, where one water molecule at the tip of the protrusion enters the organic phase by the breakup of a single hydrogen bond.

Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited

The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between -0.9 and -0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.

Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: an accurate correction scheme for electrostatic finite-size effects

The calculation of a protein-ligand binding free energy based on molecular dynamics (MD) simulations generally relies on a thermodynamic cycle in which the ligand is alchemically inserted into the system, both in the solvated protein and free in solution. The corresponding ligand-insertion free energies are typically calculated in nanoscale computational boxes simulated under periodic boundary conditions and considering electrostatic interactions defined by a periodic lattice-sum. This is distinct from the ideal bulk situation of a system of macroscopic size simulated under non-periodic boundary conditions with Coulombic electrostatic interactions. This discrepancy results in finite-size effects, which affect primarily the charging component of the insertion free energy, are dependent on the box size, and can be large when the ligand bears a net charge, especially if the protein is charged as well. This article investigates finite-size effects on calculated charging free energies using as a test case the binding of the ligand 2-amino-5-methylthiazole (net charge +1 e) to a mutant form of yeast cytochrome c peroxidase in water. Considering different charge isoforms of the protein (net charges -5, 0, +3, or +9 e), either in the absence or the presence of neutralizing counter-ions, and sizes of the cubic computational box (edges ranging from 7.42 to 11.02 nm), the potentially large magnitude of finite-size effects on the raw charging free energies (up to 17.1 kJ mol(-1)) is demonstrated. Two correction schemes are then proposed to eliminate these effects, a numerical and an analytical one. Both schemes are based on a continuum-electrostatics analysis and require performing Poisson-Boltzmann (PB) calculations on the protein-ligand system. While the numerical scheme requires PB calculations under both non-periodic and periodic boundary conditions, the latter at the box size considered in the MD simulations, the analytical scheme only requires three non-periodic PB calculations for a given system, its dependence on the box size being analytical. The latter scheme also provides insight into the physical origin of the finite-size effects. These two schemes also encompass a correction for discrete solvent effects that persists even in the limit of infinite box sizes. Application of either scheme essentially eliminates the size dependence of the corrected charging free energies (maximal deviation of 1.5 kJ mol(-1)). Because it is simple to apply, the analytical correction scheme offers a general solution to the problem of finite-size effects in free-energy calculations involving charged solutes, as encountered in calculations concerning, e.g., protein-ligand binding, biomolecular association, residue mutation, pKa and redox potential estimation, substrate transformation, solvation, and solvent-solvent partitioning.

Computation of methodology-independent single-ion solvation properties from molecular simulations. III. Correction terms for the solvation free energies, enthalpies, entropies, heat capacities, volumes, compressibilities, and expansivities of solvated ions

The raw single-ion solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions (finite or periodic system, system or box size) and treatment of electrostatic interactions (Coulombic, lattice-sum, or cutoff-based) used during these simulations. However, as shown by Kastenholz and Hünenberger [J. Chem. Phys. 124, 224501 (2006)], correction terms can be derived for the effects of: (A) an incorrect solvent polarization around the ion and an incomplete or/and inexact interaction of the ion with the polarized solvent due to the use of an approximate (not strictly Coulombic) electrostatic scheme; (B) the finite-size or artificial periodicity of the simulated system; (C) an improper summation scheme to evaluate the potential at the ion site, and the possible presence of a polarized air-liquid interface or of a constraint of vanishing average electrostatic potential in the simulated system; and (D) an inaccurate dielectric permittivity of the employed solvent model. Comparison with standard experimental data also requires the inclusion of appropriate cavity-formation and standard-state correction terms. In the present study, this correction scheme is extended by: (i) providing simple approximate analytical expressions (empirically-fitted) for the correction terms that were evaluated numerically in the above scheme (continuum-electrostatics calculations); (ii) providing correction terms for derivative thermodynamic single-ion solvation properties (and corresponding partial molar variables in solution), namely, the enthalpy, entropy, isobaric heat capacity, volume, isothermal compressibility, and isobaric expansivity (including appropriate standard-state correction terms). The ability of the correction scheme to produce methodology-independent single-ion solvation free energies based on atomistic simulations is tested in the case of Na(+) hydration, and the nature and magnitude of the correction terms for derivative thermodynamic properties is assessed numerically.

Quantum mechanical/molecular mechanical modeling finds Diels-Alder reactions are accelerated less on the surface of water than in water

Quantum and molecular mechanics calculations for the Diels-Alder reactions of cyclopentadiene with 1,4-naphthoquinone, methyl vinyl ketone, and acrylonitrile have been carried out at the vacuum-water interface and in the gas phase. In conjunction with previous studies of these cycloadditions in dilute solution, a more complete picture of aqueous environmental effects emerges with implications for the origin of observed rate accelerations using heterogeneous aqueous suspensions, "on water" conditions. The pure TIP4P water slab maintains the bulk density and hydrogen-bonding properties in central water layers. The bulk region merges to vacuum over a ca. 5 A band with progressive diminution of the density and hydrogen bonding. The relative free energies of activation and transition structures for the reactions at the interface are found to be intermediate between those calculated in the gas phase and in bulk water; i.e., for the reaction with 1,4-naphthoquinone, the DeltaDeltaG(++) values relative to the gas phase are -3.6 and -7.3 kcal/mol at the interface and in bulk water, respectively. Thus, the results do not support the notion that a water surface is more effective than bulk water for catalysis of such pericyclic reactions. The trend is in qualitative agreement with expectations based on density considerations and estimates of experimental rate constants for the gas phase, a heterogeneous aqueous suspension, and a dilute aqueous solution for the reaction of cyclopentadiene with methyl vinyl ketone. Computed energy pair distributions reveal a uniform loss of 0.5-1.0 hydrogen bond for the reactants and transition states in progressing from bulk water to the vacuum-water interface. Orientational effects are apparent at the surface; e.g., the carbonyl group in the methyl vinyl ketone transition structure is preferentially oriented into the surface. Also, the transition structure for the 1,4-naphthoquinone case is buried more in the surface, and the free energy of activation for this reaction is most similar to the result in bulk water.

In silico prediction of drug solubility: 4. Will simple potentials suffice?

In view of the extreme importance of reliable computational prediction of aqueous drug solubility, we have established a Monte Carlo simulation procedure which appears, in principle, to yield reliable solubilities even for complex drug molecules. A theory based on judicious application of linear response and mean field approximations has been found to reproduce the computationally demanding free energy determinations by simulation while at the same time offering mechanistic insight. The focus here is on the suitability of the model of both drug and solvent, i.e., the force fields. The optimized potentials for liquid simulations all atom (OPLS-AA) force field, either intact or combined with partial charges determined either by semiempirical AM1/CM1A calculations or taken from the condensed-phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field has been used. The results illustrate the crucial role of the force field in determining drug solubilities. The errors in interaction energies obtained by the simple force fields tested here are still found to be too large for our purpose but if a component of this error is systematic and readily removed by empirical adjustment the results are significantly improved. In fact, consistent use of the OPLS-AA Lennard-Jones force field parameters with partial charges from the COMPASS force field will in this way produce good predictions of amorphous drug solubility within 1 day on a standard desktop PC. This is shown here by the results of extensive new simulations for a total of 47 drug molecules which were also improved by increasing the water box in the hydration simulations from 500 to 2000 water molecules.

Intrinsic structure of hydrophobic surfaces: the oil-water interface

We investigate the water-oil interface using molecular dynamics simulations of realistic models of alkanes and water. The intrinsic density profiles are computed using a methodology that removes the smoothing effect of the capillary waves. We show that at 300 K the intrinsic width of the gap separating the oil and water phases spans little more than one water molecule diameter, and undergoes very weak short-ranged fluctuations, indicating that the water-oil interface is a rigid molecular structure at ambient temperature. Only near the drying transition (above 500 K for dodecane), the gap features uncoupled fluctuations of the oil and water surfaces, as expected in a typical drying structure. We find that the intrinsic structure of water next to the oil phase is remarkably similar to the bare water-vapor interface.

On the nature of ions at the liquid water surface

A qualitatively new understanding of the nature of ions at the liquid water surface is emerging. Traditionally, the characterization of liquid surfaces has been limited to macroscopic experimental techniques such as surface tension and electrostatic potential measurements, wherein the microscopic picture then has to be inferred by applying theoretical models. Because the surface tension of electrolyte solutions generally increases with ion concentration, all inorganic ions have been thought to be repelled from the air-water interface, leaving the outermost surface layer essentially devoid of ions. This oversimplified picture has recently been challenged: first by chemical kinetics measurements, then by theoretical molecular dynamics simulations using polarizable models, and most recently by new surface sensitive experimental observations. Here we present an overview of the nature of the interfacial structure of electrolyte solutions and give a detailed description of the new picture that is emerging.

In Silico Prediction of Drug Solubility: 1. Free Energy of Hydration

As a first step in the computational prediction of drug solubility the free energy of hydration, DeltaG*(vw) in TIP4P water has been computed for a data set of 48 drug molecules using the free energy of perturbation method and the optimized potential for liquid simulations all-atom force field. The simulations were performed in two steps, where first the Coulomb and then the Lennard-Jones interactions between the solute and the water molecules were scaled down from full to zero strength to provide physical understanding and simpler predictive models. The results have been interpreted using a theory assuming DeltaG*(vw) = A(MS)gamma + E(LJ) + E(C)/2 where A(MS) is the molecular surface area, gamma is the water-vapor surface tension, and E(LJ) and E(C) are the solute-water Lennard-Jones and Coulomb interaction energies, respectively. It was found that by a proper definition of the molecular surface area our results as well as several results from the literature were found to be in quantitative agreement using the macroscopic surface tension of TIP4P water. This is in contrast to the surface tension for water around a spherical cavity that previously has been shown to be dependent on the size of the cavity up to a radius of approximately 1 nm. The step of scaling down the electrostatic interaction can be represented by linear response theory.

Wetting Properties of Passivated Metal Nanocrystals at Liquid−Vapor Interfaces: A Computer Simulation Study

Molecular dynamics simulations have been employed to determine the contact angles of alkylthiol passivated gold nanocrystals adsorbed at the air-water interface. Simulations were performed using butane-, dodecane-, and octadecanethiol passivated nanoparticles. We demonstrate how the length of the surfactant chain can profoundly influence the wetting behavior of these nanoparticles. All particles were found to be stable at the air-water interface, possessing large, well-defined contact angles. We find that the shape of the dodecane- and octadecanethiol particles is strongly perturbed by the interface. We also present an analysis of the orientational ordering of water molecules at the dodecane-water interface and around butane- and dodecanethiol passivated nanoparticles. The orientational ordering translates into an electrostatic field around the nanoparticles, the magnitude of which corresponds with that of the water liquid-vapor interface.

Computation of methodology-independent ionic solvation free energies from molecular simulations. II. The hydration free energy of the sodium cation

The raw ionic solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions (finite or periodic system, system shape, and size) and treatment of electrostatic interactions (Coulombic, lattice sum, or cutoff based) used during these simulations. In the present article, it is shown that correction terms can be derived for the effect of (A) an incorrect solvent polarization around the ion due to the use of an approximate (not strictly Coulombic) electrostatic scheme; (B) the finite size or artificial periodicity of the simulated system; (C) an improper summation scheme to evaluate the potential at the ion site and the possible presence of a liquid-vacuum interface in the simulated system. Taking the hydration free energy of the sodium cation as a test case, it is shown that the raw solvation free energies obtained using seven different types of boundary conditions and electrostatic schemes commonly used in explicit-solvent simulations (for a total of 72 simulations differing in the corresponding simulation parameters) can be corrected so as to obtain a consistent value for this quantity.

Computation of methodology-independent ionic solvation free energies from molecular simulations. I. The electrostatic potential in molecular liquids

The computation of ionic solvation free energies from atomistic simulations is a surprisingly difficult problem that has found no satisfactory solution for more than 15 years. The reason is that the charging free energies evaluated from such simulations are affected by very large errors. One of these is related to the choice of a specific convention for summing up the contributions of solvent charges to the electrostatic potential in the ionic cavity, namely, on the basis of point charges within entire solvent molecules (M scheme) or on the basis of individual point charges (P scheme). The use of an inappropriate convention may lead to a charge-independent offset in the calculated potential, which depends on the details of the summation scheme, on the quadrupole-moment trace of the solvent molecule, and on the approximate form used to represent electrostatic interactions in the system. However, whether the M or P scheme (if any) represents the appropriate convention is still a matter of on-going debate. The goal of the present article is to settle this long-standing controversy by carefully analyzing (both analytically and numerically) the properties of the electrostatic potential in molecular liquids (and inside cavities within them). Restricting the discussion to real liquids of "spherical" solvent molecules (represented by a classical solvent model with a single van der Waals interaction site), it is concluded that (i) for Coulombic (or straight-cutoff truncated) electrostatic interactions, the M scheme is the appropriate way of calculating the electrostatic potential; (ii) for non-Coulombic interactions deriving from a continuously differentiable function, both M and P schemes generally deliver an incorrect result (for which an analytical correction must be applied); and (iii) finite-temperature effects, including intermolecular orientation correlations and a preferential orientational structure in the neighborhood of a liquid-vacuum interface, must be taken into account. Applications of these results to the computation methodology-independent ionic solvation free energies from molecular simulations will be the scope of a forthcoming article.

Electrochemical Studies of the Lateral Diffusion of TEMPO in the Aqueous Liquid/Vapor Interfacial Region

Surface partitioning and lateral mobility of TEMPO (2,2,6,6-tetramethyl-1-piperidynyloxy free radical) in the aqueous liquid/gas interfacial region were investigated electrochemically with 100 nm wide, 1.0 cm long microband electrodes positioned at the air/water interface. For redox active amphiphiles such as TEMPO, the electrochemical current is the sum of the surface and solution components representing the diffusive transport of TEMPO in both domains as well as the dynamics of equilibration at the air/water interface. Interpretation of the recorded current-voltage curves was aided by a FEMLAB simulation code developed to analyze transport processes in this class of systems. TEMPO and TEMPO(+) partition constants (K(T), K(T)+) and solution diffusivities (D(sol), equal for the two species) were obtained experimentally yielding K(T) = 5.0 +/- 0.7 x 10(2) M(-1), K(T+) = 41 +/- 3 M(-1), and D(sol) = 7.7 +/- 0.35 x 10(-6) cm(2)/s. In view of the low value of K(T)+, TEMPO(+) was assumed not to partition to the air/water interface. We further assumed that the desorption rate constants (k(des)) of both TEMPO and TEMPO(+) were the same. Good fits between the recorded and simulated cyclic voltammograms were obtained using two correlated, adjustable parameters, k(des) and the TEMPO lateral, surface diffusion constant (D(surf)). Detailed analysis of the behavior of this class of systems was obtained for a broad range of D(surf) and k(des) values. In addition, a calibration curve of k(des) versus D(surf) was obtained. Assuming that TEMPO k(des) is in a likely range of 10-100 s(-1), its lateral diffusion constant is in the range of 7.9-3.6 x 10(-5) cm(2)/s. In view of our earlier work (Wu, D. G.; Malec, A. D.; Head-Gordon, M.; Majda, M. J. Am. Chem. Soc. 2005, 127, 4490-4496) showing that at the air/water interface TEMPO is unimmersed, and that its interactions with water are limited to hydrogen bonding with one or two water molecules, D(surf) can be related to the viscosity of the aqueous interfacial region.

Theoretical Studies of Solute Vibrational Energy Relaxation at Liquid Interfaces

Recent advances in the theoretical understanding of solute vibrational energy relaxation at liquid interfaces and surfaces are described. Non-equilibrium molecular dynamics simulations of the relaxation of an initially excited solute molecule are combined with equilibrium force autocorrelation calculations to gain insight into the factors that influence the vibrational relaxation rate. Diatomic and triatomic nonpolar, polar, and ionic solute molecules adsorbed at the liquid/vapor interface of several liquids as well as at the water/CCl(4) liquid/liquid interface are considered. In general, the vibrational relaxation rate is significantly slower (a factor of 3 to 4) at the liquid/vapor and liquid/liquid interface than in the bulk due to the reduced density, which gives rise to a reduced contribution of the repulsive solvent-solute forces on the vibrational mode. The surface effects on the ionic solutes are much smaller (50% or less slower relaxation relative to the bulk). This is due to the fact that ionic solutes at the interface are able to keep part of their solvation shell to a degree that depends on their size. Thus, a significant portion of the repulsive forces is maintained. A high degree of correlation is found between the peak height of the solvent-solute radial distribution function and the vibrational relaxation rate. The relaxation rate at the liquid/liquid interface strongly depends on the location of the solute across the interface and correlates with the change in the density and polarity profile of the interface.

Hydration Shell Exchange Dynamics during Ion Transfer Across the Liquid/Liquid Interface

We examine using molecular dynamics simulations the rate and mechanism of water molecules exchange around the Li(+) and Na(+) ions during ion transfer across the interface between water and nitrobenzene. As the ions are transferred from the water to the organic phase, they keep their first hydration shell and an incomplete second shell. The rate of water exchange between the first shell and the rest of the interfacial water molecule decreases during the transfer, which is consistent with an increase in the barrier along the ion-water potential of mean force. While in bulk water the exchange of water molecules around the Li(+) follows an associative (A) or associative interchange (I(a)) type mechanism, the fraction of exchange events of type A increases at the interface. In contrast, while in bulk water the exchange of water molecules around the six coordinated Na(+) hydrated species mainly follows a dissociative mechanism, the situation at the interface involves an equilibrium interchange between the four- and five-coordinated hydrated ion. Simulation of the reversed process, in which the hydrated Li(+) ion is transferred to the aqueous phase, shows the same general behavior as a function of location from the interface.

Viscosity of the Aqueous Liquid/Vapor Interfacial Region: 2D Electrochemical Measurements with a Piperidine Nitroxy Radical Probe

Surface partitioning of 2,2,6,6-tetramethyl-1-piperidynyloxy radical (Tempo) to the air/water interface follows a Langmuir isotherm. The partition constant was obtained by the surface tension measurements in the concentration range of 1.0 x 10(-4)-2.4 x 10(-3) M yielding K = 640 +/- 99 M(-1). The lateral mobility of Tempo at the air/water interface was measured electrochemically in the surface concentration range of 2.0 x 10(-11)-1.4 x 10(-10) mol/cm2, corresponding to ca. 7.3-50% full monolayer coverage. The measurements employed cyclic voltammetry with line microelectrodes touching the air/water interface. The Tempo lateral diffusion constant of (1.5 +/- 0.7) x 10(-4) cm2/s is independent of surface concentration below 4.0 x 10(-11) mol/cm2. The extent of Tempo water interactions was assessed by the electronic structure calculations. These calculations showed that, at most, two water molecules can hydrogen bond with the oxygen atom of the nitroxyl group of Tempo, and that a single water molecule forms a hydrogen bond that is ca. 30% stronger than the H2O-H2O hydrogen bond. These calculations led to a postulate that Tempo diffuses along the interface largely unimmersed, and that it is coupled to the interfacial water via hydrogen bonding with H2O. In view of this postulate, the viscosity of the aqueous liquid/vapor interfacial region obtained by interpreting the Tempo diffusion constant in the low concentration region is as much as 4 times smaller than that of bulk liquid water.

Deblurred Observation of the Molecular Structure of an Oil−Water Interface

We simulated the interface between liquid water and a stationary phase of tethered n-C18 alkyl chains at a thermodynamic state of low pressure and water vapor-liquid coexistence. The interfacial water (oxygen atom) density profile so obtained is compared with a precisely defined proximal density of water molecules (oxygen atoms) conditional on the alkyl chain configurations. Though the conventional interfacial density profile takes a traditional monotonic form, the proximal radial distribution of oxygen atoms around a specific methyl (methylene) group closely resembles that for a solitary methane solute in liquid water. Moreover, this proximal radial distribution function is sufficient to accurately reconstruct the water oxygen density profile of the oil-water interface. These observations provide an alternative interpretation to collective drying or vaporization interpretations of commonly observed oil-water interfacial profiles for which water penetration into the interfacial region plays a role.

Rational determination of charge distributions for free energy calculations

Point charges derived from RHF/6-31G* electrostatic potentials are attractive because they tend to exaggerate the polarity of solvated molecules, thereby compensating in an average fashion missing induction effects. In the context of free energy calculations, wherein the molecule is transferred from a polar environment to a nonpolar one, we propose a more rational approach based on a self-consistent reaction field computation at a higher level of theory, supplemented by an estimation of the corresponding distortion energy to account for the change of polarity of the surroundings. Application of this method to the test cases acetamide, acetic acid, methyl acetate and phenol, using multinanosecond molecular dynamics/"umbrella sampling" simulations, yields consistent hydration free energies in reasonably good agreement with experiment.

Correlation between hydrophobic properties and efficiency of carrier-mediated membrane transduction and apoptosis of a p53 C-terminal peptide

Two membrane transporters, the 17 amino acid (aa) oligopeptide penetratin derived from the homeodomain of Antennapedia (Ant) and an analogue of the basic domain of TAT (aa 47-57) (TAT-a) from HIV-1, were tested as carriers for a p53 C-terminal peptide (aa 361-382) into human breast cancer cells. The studies were performed to determine whether the membrane-transduction efficiency of membrane carriers: Ant, TAT or TAT analogue (TAT-a) correlated with peptide hydrophobic features. Peptide-sequence analysis clearly demonstrated that the Ant sequence and p53 peptide sequence (p53p) together created a peptide with enhanced hydrophobic characteristics; while the TAT or TAT analogue (TAT-a) and p53p sequence together created a peptide with significantly less hydrophobic qualities. The degree of hydrophobic moment and helical wheel plots for these peptides correlated directly with their ability to transduce the p53 peptide. Western blot analysis revealed that Ant was able to transduce p53 C-terminal peptide into human breast cancer cells as a highly efficient membrane transporter. Compared to Ant, TAT-a fused to the C-terminus of p53 peptide (p53p-TAT-a) was a less efficient carrier into these cells under the conditions of our study. Additionally, N-terminal linked TAT-a to p53p (TAT-a-p53p) showed even lower efficiency as a transporter than p53-TAT-a. Apoptosis assays showed that the p53 peptide, fused at its C-terminus to Ant (p53p-Ant), induced a higher percentage of apoptotic cells in human breast cancer cell lines expressing mutant or wild-type p53 as compared to p53 peptide fused at its C-terminus to the TAT-a sequence (p53p-TAT-a) or when fused at the N-terminus to TAT-a (TAT-a-p53p). These data suggested a direct correlation between hydrophobic characteristics and efficiency as a transporter. Sequence study, using hydrophobic moment and helical wheel analyses, may be useful predictive tools for choosing the best carrier for a peptide.

Molecular theory of hydrophobic effects: "She is too mean to have her name repeated."

This paper reviews the molecular theory of hydrophobic effects relevant to biomolecular structure and assembly in aqueous solution. Recent progress has resulted in simple, validated molecular statistical thermodynamic theories and clarification of confusing theories of decades ago. Current work is resolving effects of wider variations of thermodynamic state, e.g., pressure denaturation of soluble proteins, and more exotic questions such as effects of surface chemistry in treating stability of macromolecular structures in aqueous solution.