Microscopic wetting of self-assembled monolayers with different surfaces: a combined molecular dynamics and quantum mechanics study

Molecular dynamics simulations are used to study the micronature of the organization of water molecules on the flat surface of well-ordered self-assembled monolayers (SAMs) of 18-carbon alkanethiolate chains bound to a silicon (111) substrate. Six different headgroups (-CH(3), -C═C, -OCH(3), -CN, -NH(2), -COOH) are used to tune the character of the surface from hydrophobic to hydrophilic, while the level of hydration is consistent on all six SAM surfaces. Quantum mechanics calculations are employed to optimize each alkyl chain of the different SAMs with one water molecule and to investigate changes in the configuration of each headgroup under hydration. We report the changes of the structure of the six SAMs with different surfaces in the presence of water, and the area of the wetted surface of each SAM, depending on the terminal group. Our results suggest that a corrugated and hydrophobic surface will be formed if the headgroups of SAM surface are not able to form H-bonds either with water molecules or between adjacent groups. In contrast, the formation of hydrogen bonds not only among polar heads but also between polar heads and water may enhance the SAM surface hydrophilicity and corrugation. We explicitly discuss the micromechanisms for the hydration of three hydrophilic SAM (CN-, NH(2)- and COOH-terminated) surfaces, which is helpful to superhydrophilic surface design of SAM in biomimetic materials.     

The hydrophobic effect: molecular dynamics simulations of water confined between extended hydrophobic and hydrophilic surfaces

Structural and dynamic properties of water confined between two parallel, extended, either hydrophobic or hydrophilic crystalline surfaces of n-alkane C(36)H(74) or n-alcohol C(35)H(71)OH, are studied by molecular dynamics simulations. Electron density profiles, directly compared with corresponding experimental data from x-ray reflectivity measurements, reveal a uniform weak de-wetting characteristic for the extended hydrophobic surface, while the hydrophilic surface is weakly wetted. These microscopic data are consistent with macroscopic contact angle measurements. Specific water orientation is present at both surfaces. The ordering is characteristically different between the surfaces and of longer range at the hydrophilic surface. Furthermore, the dynamic properties of water are different at the two surfaces and different from the bulk behavior. In particular, at the hydrophobic surface, time-correlation functions reveal that water molecules have characteristic diffusive behavior and orientational ordering due to the lack of hydrogen bonding interactions with the surface. These observations suggest that the altered dynamical properties of water in contact with extended hydrophobic surfaces together with a partial drying of the surfaces are more indicative of the hydrophobic effect than structural ordering, which we suggest to be independent of surface topology.     

Polypeptide friction and adhesion on hydrophobic and hydrophilic surfaces: a molecular dynamics case study

Using all-atomistic MD simulations including explicit water, the mobility and adhesion of a mildly hydrophobic single polypeptide chain adsorbed on hydrophobic and hydrophilic diamond surfaces is investigated by application of lateral and vertical pulling forces. Forced motion on the hydrophilic surface exhibits stick-slip due to breaking and reformation of hydrogen bonds; in contrast, on the hydrophobic surface, the motion is smooth. By carefully tuning the driving force magnitude, the linear-response regime is reached on a hydrophobic surface and equilibrium values for mobility and adhesive strength are obtained. On the hydrophilic surface, on the other hand, slow hydrogen-bond kinetics prevents equilibration and only upper bounds for adhesion force and mobility can be estimated. Whereas the desorption force is rather comparable on the two surfaces and differs at most by a factor of 2, the mobility on the hydrophilic surface is at least 30-fold reduced compared to the hydrophobic one. A simple model based on a single particle diffusing in a corrugated potential landscape suggests that cooperativity is rather limited and that the small mobility on a hydrophilic surface can be rationalized in terms of incoherently moving monomers. The experimentally well-known peptide mobility in bulk water is quantitatively reproduced in our simulations, which serves as a sensitive test on our methodology employed.     

Chemical Force Microscopy Study of Adhesion and Friction between Surfaces Functionalized with Self-Assembled Monolayers and Immersed in Solvents

Adhesive and frictional forces between surfaces modified with self-assembled monolayers (SAMs) and immersed in solvents were measured with chemical force microscopy as functions of surface functionality and solvent. Si/SiO2 substrates were modified with SAMs of alkylsiloxanes (SiCl3(CH2)n-X), and gold-coated AFM tips were modified with SAMs of alkylthiolates (HS-(CH2)n-X). SAMs of alkylsiloxanes terminated in a methyl or oxidized vinyl group; SAMs of alkanethiolates terminated in a methyl or carboxyl group. Adhesive and frictional forces were measured in hexadecane, ethanol, 1,2-propanediol, 1,3-propanediol, and water. The work of adhesion (W) was calculated with the Johnson-Kendall-Roberts theory of adhesive contact. The JKR values agreed well with values derived from the Fowkes-van Oss-Chaudhury-Good surface tension model and from contact angle results. Calculated values of W for all combinations of contacting surfaces and solvents spanned two orders of magnitude. W correlated with the surface tension of the solvent for hydrophobic/hydrophobic interactions; hydrophilic/hydrophilic and hydrophobic/hydrophilic interactions were more complex. Friction forces were fit to a modified form of Amonton's law. For any solvent, friction coefficients were largest for the hydrophilic/hydrophilic contacting surfaces. The friction coefficient for any contacting pair was largest in hexadecane. In polar solvents, friction coefficients scaled with solvent polarity only for hydrophobic/hydrophobic contacting pairs. Copyright 1999 Academic Press.     

State and structure of water in vicinity of hydrophobic surfaces

Criterial values of the specific heat of water wetting, surface pressure, and contact angle classifying surfaces into hydrophilic and hydrophobic are proposed based on the analysis of own and published data. The most characteristic properties of hydrophobic surfaces, i.e., large surface area per water molecule in the conventional adsorption monolayer and the absence of continuous two-layer water film on the adsorbent surface at vapor pressure close to saturation, are discussed using nonporous carbon-based materials as example. The presence of residual hydrophilic groups that act as sites of the clusterization of polar molecules on the surface of graphitized carbon black is confirmed by gas chromatography and the concentration of these sites is calculated. The amount of water molecules in the surface cluster is determined at different stages of adsorption. Procedures for preparing organically modified layered silicates and silica as basic objects of the study of the interaction between water molecules and hydrophobic surfaces are considered. It is proven that the boundary water layer in the vicinity of hydrophobic surface consists of a thin (∼0.5 nm) depletion layer with a density of 0.4 g/cm³ and a considerable amount (25–30%) of water molecules with free OH groups and thicker (∼35 nm) layer, which is characterized by a more ordered network of hydrogen bonds compared to liquid water. Data obtained by X-ray scattering and neutron and reflection methods, and sum-frequency vibrational spectroscopy are compared with the results of calorimetric study of the interaction between water and hydrophobic surface, as well as with the data of molecular-statistical calculations of the state of water molecules in the surface layer.     

Water hydrogen bond analysis on hydrophilic and hydrophobic biomolecule sites

Elastic and quasielastic neutron scattering experiments have been used to investigate the hydrogen bonding network dynamics of hydration water on hydrophilic and hydrophobic sites. To this end the evolution of hydration water dynamics of a prototypical hydrophobic amino acid with polar backbone, N-acetyl-leucine-methylamide (NALMA), and hydrophilic amino acid, N-acetyl-glycine-methylamide (NAGMA), has been investigated as a function of the molecular ratio water : peptide. The results suggest that the dynamical contribution of the intrinsic and low hydration molecules of water is characteristic of pure librational/rotational movement. The water molecule remains attached to the hydrophilic site with only the possibility of hindered rotations that eventually break the bond with the peptide and reform it immediately after. A gradual evolution from librational motions to hindered rotations is observed as a function of temperature. When the hydration increases, we observe (together with the hindered rotations of hydrogen bonds) a slow diffusion of water molecules on the surface of the peptides.     

Unusual properties of water at hydrophilic/hydrophobic interfaces

The behaviour of water at mosaic hydrophilic/hydrophobic surfaces of different silicas and in biosystems (biomacromolecules, yeast cells, wheat seeds, bone and muscular tissues) was studied in different dispersion media over wide temperature range using 1H NMR spectroscopy with layer-by-layer freezing-out of bulk water (close to 273 K) and interfacial water (180 < T < 273 K), thermally stimulated depolarization current (TSDC) (90 < T < 270 K), infrared (IR) spectroscopy, and quantum chemical methods. Bulk water and water bound to hydrophilic/hydrophobic interfaces can be assigned to different structural types. There are (i) weakly associated interfacial water (1H NMR chemical shift delta(H) = 1.1-1.7 ppm) that can be assigned to high-density water (HDW) with collapsed structure (CS), representing individual molecules in hydrophobic pockets, small clusters and interstitial water with strongly distorted hydrogen bonds or without them, and (ii) strongly associated interfacial water (delta(H) = 4-5 ppm) with larger clusters, nano- and microdomains, and continuous interfacial layer with both HDW and low-density water (LDW). The molecular mobility of weakly associated bound water is higher (because hydrogen bonds are distorted and weakened and their number is smaller than that for strongly associated water) than that of strongly associated bound water (with strong hydrogen bonds but nevertheless weaker than that in ice Ih) that results in the difference in the temperature dependences of the 1H NMR spectra at T < 273 K. These different waters are also appear in changes in the IR and TSDC spectra.     

Effect of nanostructure on the properties of water at the water-hydrophobic interface: a molecular dynamics simulation

The local structure of water near hydrophobic surfaces of different surface topographies has been analyzed by molecular dynamics simulation. An alkane crystal has been taken as the parent model for a hydrophobic surface. Surface structures were created by placing pits into it, which were half a nanometer deep and several nanometers wide. Around all structures, the water has a lower density, less orientational ordering, fewer water-water hydrogen bonds, and fewer surface contacts than for a flat unstructured surface. This indicates that the structured surfaces are more hydrophobic than the flat surface. Of the structures investigated, pits with a diameter of approximately 2.5 nm were effective in increasing the hydrophobic character of the surface.     

Influence of nanobubbles on the adsorption of nanoparticles

A quartz crystal microbalance was used to study the influence of nanobubbles on the adsorption of polystyrene nanoparticles onto surfaces coated with gold, or coated with dodecanethiol or mercaptoundecanoic acid self-assembled monolayers (SAMs). Adsorption of the nanoparticles onto the surface causes the resonant frequency of the quartz crystal to decrease. We found that particles were adsorbed onto the gold-coated quartz crystal in air-rich water, but not in degassed water. This finding supports the long-standing hypothesis that nanobubbles play a key role in the long-range attractive force between hydrophobic surfaces in aqueous solutions. When the experiments were conducted using quartz crystals coated with a hydrophobic dodecanethiol SAM, the nanoparticles were adsorbed onto the surface even in degassed water due to the short-range hydrophobic interactions between the nanoparticles and the dodecanethiol molecules. In contrast, the nanoparticles were adsorbed to a lesser degree onto the hydrophilic mercaptoundecanoic acid-coated crystals due to electrostatic repulsive forces.     

Aqueous sodium oleate–sodium dehydrocholate mixtures at low concentration

The aqueous mixed system sodium dehydrocholate (NaDHC)–sodium oleate (NaOL) was studied by several methods to determine the influence of the hydrophobic structure of both surfactants in the mixed micellization and the formation of the mixed monolayer adsorbed at the air–water interface. The molecular area at the critical micelle concentration in pure surfactant solutions suggests that the adsorbed oleate chain was folded to allow the double bond in the middle of the molecule to remain in contact with water, and that the NaDHC molecule was situated with its plane laying parallel to the water surface, allowing the three carbonyl groups in the hydrocarbon backbone to form hydrogen bonds with water. The interaction was repulsive at the surface, and in the mixed monolayer some molecules must move away the less hydrophilic groups from water (double bond of NaOL, carbonyl groups of NaDHC). The interaction in mixed micelles was strongly attractive, showing a preferential composition roughly equimolar. The hydrolysis in mixed micelles was augmented in comparison with pure surfactants systems, which could be explained by assuming the existence of a more hydrophobic mixed micelle core. The mixed micelle degree of ionization was below that of the pure micelles, thus indicating a high surface charge density.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Dynamics of water confined in the interdomain region of a multidomain protein

Molecular dynamics simulations are performed to study the dynamics of interfacial water confined in the interdomain region of a two-domain protein, BphC enzyme. The results show that near the protein surface the water diffusion constant is much smaller and the water-water hydrogen bond lifetime is much longer than that in bulk. The diffusion constant and hydrogen bond lifetime can vary by a factor of as much as 2 in going from the region near the hydrophobic domain surface to the bulk. Water molecules in the first solvation shell persist for a much longer time near local concave sites than near convex sites. Also, the water layer survival correlation time shows that on average water molecules near the extended hydrophilic surfaces have longer residence times than those near hydrophobic surfaces. These results indicate that local surface curvature and hydrophobicity have a significant influence on water dynamics.     

Roughness and ordering at the interface of oxidized polystyrene and water

For the first time, atomistically detailed molecular dynamics calculations revealed molecular ordering of the water-oxidized atactic polystyrene (aPS) interface. Both ordering of the water molecules and the phenyl rings occur. In addition, the natural roughness of the surface has been simulated and compared to experimental values. The composition of the simulated aPS films is based on spin-coated aPS films that have been oxidized and characterized experimentally. The aPS surfaces are oxidized with ultraviolet-ozone radiation and have been characterized by XPS, AFM, and water contact angle measurements. XPS measurements show that the oxygen content in the sample increases rapidly with exposure and reaches saturation near 24 at. % of oxygen. The molecular dynamics simulations show smoothening of an hydrophobic aPS surface upon transition from vacuum to water. The smoothening decreases with increasing hydrophilicity. The calculations reveal ordering of oxidized phenyl rings for aPS surfaces in water. The order increases with increasing hydrophilicity. Additionally, we investigated the water structure near the aPS-water interface as a function of the surface hydrophilicity. With increasing hydrophilicity, the density of water at the aPS-water interface increases. The water density profile is steeper in the presence of hydrophobic aPS. The water shows an ordered layer near both the hydrophobic and hydrophilic surfaces; the position of this layer shifts toward the interface with increasing hydrophilicity.     

Computer simulation of adhesion between hydrophilic and hydrophobic self-assembled monolayers in water

The grand canonical Monte Carlo technique and atomistic force fields are used to calculate the force-distance relations and free energies of adhesion between carboxyl and methyl terminated alkanethiolate self-assembled monolayers (SAMs) in water. Both symmetric and asymmetric confinements are considered, as formed by like and unlike SAMs, respectively. As the confinement is increased, water confined by the hydrophobic methyl terminated SAMs experiences capillary evaporation. As a consequence, the adhesion energy is determined by the direct interaction between bare SAMs. In the asymmetric system, an incomplete capillary evaporation is observed, with the number of water molecules dropped by more than an order of magnitude. The remaining water molecules are all adsorbed on the hydrophilic SAM, while the hydrophobic SAM is separated from the rest of the system by a thin vapor layer. The calculated free energies of adhesion are in acceptable agreement with experiment.     

Molecular packing of lysozyme,fibrinogen, and bovine serum albumin on hydrophilic and hydrophobic surfaces studied by infrared-visible sum frequency generation and fluorescence microscopy

Infrared-visible sum frequency generation (SFG) vibrational spectroscopy, in combination with fluorescence microscopy, was employed to investigate the surface structure of lysozyme, fibrinogen, and bovine serum albumin (BSA) adsorbed on hydrophilic silica and hydrophobic polystyrene as a function of protein concentration. Fluorescence microscopy shows that the relative amounts of protein adsorbed on hydrophilic and hydrophobic surfaces increase in proportion with the concentration of protein solutions. For a given bulk protein concentration, a larger amount of protein is adsorbed on hydrophobic polystyrene surfaces compared to hydrophilic silica surfaces. While lysozyme molecules adsorbed on silica surfaces yield relatively similar SFG spectra, regardless of the surface concentration, SFG spectra of fibrinogen and BSA adsorbed on silica surfaces exhibit concentration-dependent signal intensities and peak shapes. Quantitative SFG data analysis reveals that methyl groups in lysozyme adsorbed on hydrophilic surfaces show a concentration-independent orientation. However, methyl groups in BSA and fibrinogen become less tilted with respect to the surface normal with increasing protein concentration at the surface. On hydrophobic polystyrene surfaces, all proteins yield similar SFG spectra, which are different from those on hydrophilic surfaces. Although more protein molecules are present on hydrophobic surfaces, lower SFG signal intensity is observed, indicating that methyl groups in adsorbed proteins are more randomly oriented as compared to those on hydrophilic surfaces. SFG data also shows that the orientation and ordering of phenyl rings in the polystyrene surface is affected by protein adsorption, depending on the amount and type of proteins.     

Recognition of Hydrophilic Cyclic Compounds by a Water-Soluble Cavitand

A water-soluble deep cavitand bearing amides on the upper rim and trimethyl ammonium groups on the feet was synthesized. The open-ended cavity is stabilized by the intramolecular hydrogen bonds formed between the adjacent amides, and the introduction of trimethylammonium imparts to the cavitand good solubility in water. The cavitand exhibits high binding affinity and selectivity to hydrophilic molecules in water. With certain guests, such as cyclohexyl alcohols, amines and acids, the recognition involves the synergistic action of hydrogen bonding with hydrophobic effects. The binding phenomena are interpreted in terms of a fixed solvent cage presented by the host to the guest.     

The importance of the shape of the protein-water interface of a kinesin motor domain for dynamics of the surface atoms of the protein

A single kinesin motor domain immersed in water has been investigated using molecular dynamics. It has been found that local properties of water in the solvation shell change along with the nature of the neighboring protein surface. However, a detailed analysis leads to the conclusion that the geometrical features of hydrogen bonds and overall structure of kinesin hydration water are not very different from bulk water. The local values of diffusion coefficients (translational and rotational) of water adjacent to specific patches on the protein surface seem not to be correlated to the orientational ordering of hydration water, but instead they depend on spatial roughness and degree of exposure of the patch to the solvent. Finally, a relationship between the mobility of various surface atoms of the protein and the mean values of the diffusion coefficient of the adjacent water molecules has been observed. The latter finding suggests a close relationship between the dynamics of the inner kinesin movements and the behavior of solvation water which is in turn determined by the topography of the contact surface between the protein and the surrounding water molecules.     

Poly[triaqua(μ4‐4‐carboxybenzenesulfonato‐κ4O:O′:O′′:O′′′)(4‐carboxybenzenesulfonato‐κO)strontium(II)]

This study presents the coordination modes and two‐dimensional network of a novel strontium(II) coordination polymer, [Sr(C₇H₅O₅S)₂(H₂O)₃]_n. The eight‐coordinate Sr²⁺ ion is in a distorted bis‐disphenoidal coordination environment, surrounded by four sulfonate and one carboxyl O atom from five benzenesulfonate ligands, two of which are symmetry unique, and by three O atoms from three independent aqua ligands. The compound exhibits a monolayer structure with coordination bonds within and hydrogen bonds between the layers. The μ₄ acid ligand bridges the metal ions in two dimensions to form a thick undulating monolayer with a hydrophobic interior and hydrophilic surfaces. A second independent monoanion is arranged outward from both sides of the monolayer and serves to link adjacent monolayers via carboxyl–water and water–carboxyl hydrogen bonds.     

Structure and dynamics of hydrogen bonds in the interface of a C12E6 spherical micelle in water solution: a MD study at various temperatures

The temperature dehydration of a C(12)E(6) spherical micelle is characterized through the study of the structure and dynamics of the hydrogen bonds formed by water within the micellar interface. Water molecules in proximity of the hydrophilic fragment of the C(12)E(6) surfactants form strong H-bonds with the oxyethilene units E and with the polar alcoholic heads. The activation energies of such H-bonds fall in the range 2-3 Kcal mol(-1). On the exposed oil core, the number of water-water H-bonds decreases as an effect of dehydration. The dynamics of such bonds exhibits a slow relaxation with respect to the bulk, and two time scales can be discerned: the first one, tau approximately 3-6 ps, is typical of water-water H-bonds around small hydrophobic molecules, whereas the second one, tau approximately 40-80 ps, is probably due to the confining effect of the long hydrophilic fragments which reduces the probability of a water molecule to leave the hydration layer of the exposed oil core. Water molecules around the core form H-bond clusters whose size and distribution change with temperature. From a cluster analysis, the system appears to be below the percolation threshold, suggesting that the exposed oily surface is formed by disconnected patches of size around 1 nm(2), close to the estimate of the solvated hydrophobic patches on protein surfaces. The network connectivity is also considered for concentric hydration shells along the interface: it turns out that near the oil core, the cluster size is larger than elsewhere in the interface demonstrating a strong structural effect induced by the exposed hydrocarbon tails. Temperature affects the cluster size only in the innermost shell.     

Water-induced ethanol dewetting transition

The dewetting transitions of two hydrophobic plates immersed in pure water, aqueous ethanol solutions with concentrations from 25% to 90%, and pure ethanol were investigated by molecular dynamics simulations, where the dewetting transition was analogous to a first-order phase transition from liquid to vapor. It was found that the dewetting transitions occurred except that in the pure ethanol system. Although the ethanol molecules prefer to locate in the vicinity of the two plates, the inter-plate region is unfavorable for water molecules, due to losing more than one hydrogen bond. Moreover, each inter-plate water molecule forms hydrogen bonds on average with about two ethanol molecules. These intermolecular hydrogen bonds cause water and ethanol to cooperatively fill or exit the inter-plate region. Thus, water molecules play a more important role in the inter-plate filling/empty process, and induce the ethanol dewetting transition. Our results provide insight into the effect of water on the ethanol dewetting phenomena.