Attractions, water structure, and thermodynamics of hydrophobic polymer collapse

We explore the prospects of a perturbation approach for predicting how weak attractive interactions affect collapse thermodynamics of hydrophobic polymers in water. Specifically, using molecular dynamics simulations of model polymers in explicit water, we show that the hydration structure is sensitive to the strength of the van der Waals attractions but that the hydration contribution to the potential of mean force for collapse is not. We discuss how perturbation theory ideas developed for small spherical apolar solutes need to be modified in order to account for the effect of attractions on the conformational equilibria of polymers.     

Hydration and dewetting near graphite-CH(3) and graphite-COOH plates

The dynamics of water near the nanoscale hydrophobic (graphite-CH(3)) and hydrophilic (graphite-COOH) plates has been studied in detail with molecular dynamics simulations in this paper. It is shown that these designed surfaces (by growing a layer of methyl or carboxyl groups on top of graphite) can have a significant impact on the neighboring water dynamics, with the hydrophilic carboxyl surface having even more profound effects. The water hydrogen bond lifetime is much longer near both types of surfaces than that in the bulk, while on the other hand the water diffusion constant is much smaller than that in the bulk. The difference in the diffusion constant can be as large as a factor of 8 and the difference in the hydrogen bond lifetime can be as large as a factor of 2, depending on the distance from the surface. Furthermore, the water molecules in the first solvation shell of surface atoms show a strong bias in hydroxyl group orientation near the surface, confirming some of the previous findings. Finally, the possible water dewetting transition between two graphite-CH(3) plates and the effect of the strength of the solute-solvent attractions on the water drying transition are investigated. The relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on the graphite-CH(3) surface is also analyzed on the basis of a simple macroscopic theory, which can be used to predict the dewetting transition critical distance.     

Role of water in mediating the assembly of Alzheimer amyloid-beta Abeta16-22 protofilaments

The role of water in promoting the formation of protofilaments (the basic building blocks of amyloid fibrils) is investigated using fully atomic molecular dynamics simulations. Our model protofilament consists of two parallel beta-sheets of Alzheimer Amyloid-beta 16-22 peptides (Ac-K(16)-L(17)-V(18)-F(19)-F(20)-A(21)-E(22)-NH2). Each sheet presents a distinct hydrophobic and hydrophilic face and together self-assemble to a stable protofilament with a core consisting of purely hydrophobic residues (L(17), F(19), A(21)), with the two charged residues (K(16), E(22)) pointing to the solvent. Our simulations reveal a subtle interplay between a water mediated assembly and one driven by favorable energetic interactions between specific residues forming the interior of the protofilament. A dewetting transition, in which water expulsion precedes hydrophobic collapse, is observed for some, but not all molecular dynamics trajectories. In the trajectories in which no dewetting is observed, water expulsion and hydrophobic collapse occur simultaneously, with protofilament assembly driven by direct interactions between the hydrophobic side chains of the peptides (particularly between F-F residues). For those same trajectories, a small increase in the temperature of the simulation (on the order of 20 K) or a modest reduction in the peptide-water van der Waals attraction (on the order of 10%) is sufficient to induce a dewetting transition, suggesting that the existence of a dewetting transition in simulation might be sensitive to the details of the force field parametrization.     

Hydration and dewetting near fluorinated superhydrophobic plates

The water dynamics near nanoscale fluorinated (CF(3)(CF(2))(7)(CH(2))(2)SiH(3)) monolayers (plates) as well as possible dewetting transitions in-between two such plates have been studied with molecular dynamics simulations in this paper. A "weak water depletion" is found near the single fluorinated surface, with an average water density in the first solvation shells 6-8% lower than its hydrogenated counterpart. The fluorinated molecules are also found to be water impermeable, consistent with experimental findings. More surprisingly, a dewetting transition is found in the interplate region with a critical distance D(c) of 10 A (3-4 water diameters) for double plates with 8 x 8 molecules each (plate size approximately 4 nm x 4 nm). This transition, although occurring on a microscopic length scale, is reminiscent of a first-order phase transition from liquid to vapor. The unusual superhydrophobicity of fluorocarbons is found to be related to their larger size (or surface area) as compared to hydrocarbons, which "dilutes" their physical interactions with water. The water-plate interaction profile shows that the fluorinated carbons have a 10-12% weaker water-plate interaction than their hydrogenated counterparts in the nearest solvation shell, even though the fluorocarbons do have a stronger electrostatic interaction with water due to their larger partial charges. However, the van der Waals interactions dominate the water-plate interaction within the nearest shell, with up to 90% contributions to the total interaction energy, and fluorocarbons have a noticeably weaker (by 10-15%) van der Waals interaction with water in the nearest shell than do hydrocarbons. Both the slightly weaker water-plate interaction and larger surface area contribute to the stronger dewetting transition in the current fluorinated carbon plates.     

Studies on the interaction mechanism between CXS and papain by spectroscopic and molecular modeling methods

We have studied the interaction between cefuroxime sodium (CXS) and papain at different temperatures by a fluorescence method, and confirmed that the mechanism of fluorescence quenching of CXS to papain is mainly static quenching. We also determined the binding constant K. Based on the thermodynamic functions at different temperatures, the results show that the major forces between CXS and papain are van der Waals’ forces and H bond. According to the Forster non-radiation energy transfer mechanism, we determined the binding distance between CXS and papain, and studied the confirmation effect of CXS to papain by synchronous fluorescence and UV–Vis spectroscopy. Molecular simulations show that the binding types of CXS and papain are van der Waals’ forces, hydrophobic interaction, and H-bond.     

Hydrophobic drying and hysteresis at different length scales by molecular dynamics simulations

We performed molecular dynamics simulations to investigate hydrophobic interactions between two parallel hydrophobic plates immersed in water. The two plates are separated by a distance D ranging from contact to a few nanometers. To mimic the attractive hydrophobic force measurement in a surface force experiment, a driving spring is used to measure the hydrophobic force between two hydrophobic plates. The force-distance curves, in particular the force variations from spontaneous drying to hydrophobic collapse are obtained. These details are usually not accessible in the surface force measurement due to the unstable jump into contact. The length-scale effect on the hydrophobic drying during normal approach and the hydrophobic hysteresis during retraction has been studied. We find that the critical distance at which a spontaneous drying occurs is determined by the shorter characteristic dimension of the plate, whereas the hydrophobic hysteresis is determined by the longer characteristic dimension of the plate. The variations of the potential of mean force versus separation during approach and retraction are also calculated. The results show that water confined between two parallel hydrophobic plates is in a thermodynamic metastable state. This comparably high energy state leads to the spontaneous drying at some critical distance.     

Corresponding-states laws for protein solutions

The solvent around protein molecules in solutions is structured and this structuring introduces a repulsion in the intermolecular interaction potential at intermediate separations. We use Monte Carlo simulations with isotropic, pair-additive systems interacting with such potentials. We test if the liquid-liquid and liquid-solid phase lines in model protein solutions can be predicted from universal curves and a pair of experimentally determined parameters, as done for atomic and colloid materials using several laws of corresponding states. As predictors, we test three properties at the critical point for liquid-liquid separation: temperature, as in the original van der Waals law, the second virial coefficient, and a modified second virial coefficient, all paired with the critical volume fraction. We find that the van der Waals law is best obeyed and appears more general than its original formulation: A single universal curve describes all tested nonconformal isotropic pair-additive systems. Published experimental data for the liquid-liquid equilibrium for several proteins at various conditions follow a single van der Waals curve. For the solid-liquid equilibrium, we find that no single system property serves as its predictor. We go beyond corresponding-states correlations and put forth semiempirical laws, which allow prediction of the critical temperature and volume fraction solely based on the range of attraction of the intermolecular interaction potential.     

Communication: Length scale dependent oil-water energy fluctuations

Interfacial fluctuations in the cohesive (van der Waals) interaction energy of spherical oil-drops with water provide evidence of a length scale dependent transition from linear to non-linear response behavior. For sub-nanometer oil-drop sizes, energy fluctuations are found to be independent of the van der Waals coupling strength, while nanometer (and larger) size oil drops experience highly non-linear energy fluctuations. The latter behavior is linked to enhanced hydrophobic density fluctuations and the emergence of entropic contributions to oil-water cohesive interaction free energies.     

Modeling van der Waals interactions between proteins and inorganic surfaces from time-dependent density functional theory calculations

In this work we show how the ab initio determination of van der Waals coefficients within time-dependent density functional theory can be used to build efficient and accurate atomistic models that describe the long-range interactions of proteins with other proteins and of proteins with semi-conducting surfaces. The model parameters are fitted so that they reproduce the ab initio van der Waals coefficients of amino acids and dipeptides. We then assess the quality of our results by comparing ab initio van der Waals coefficients for larger peptides with the coefficients yielded by the models. The different sets of parameters can be easily incorporated in current empirical force field methods, thus providing an essential ingredient for molecular dynamics simulations of proteins close to surfaces.     

Specific ion-dependent attraction and phase behavior of polymer-coated colloids

This paper presents results demonstrating the role of temperature and specific ions in mediating attraction between polymer-coated colloids and determining their equilibrium phase behavior. In particular, theoretical predictions of continuum van der Waals attraction between poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO)-coated polystyrene colloids are used to explain measured temperature and specific ion-dependent fluid-gel transitions in dispersions of these particles. Building on previous studies of PEO-PPO-PEO-coated polystyrene colloids dispersed in aqueous NaCl media, this work reports rheologically measured fluid-gel transitions as a function of temperature and NaCl/MgSO4 composition. Adhesive-sphere predictions of percolation thresholds are fit to measured fluid-gel data by allowing the adsorbed copolymer layer thickness as a single adjustable parameter. This allows the attraction between the PEO-PPO-PEO layers to be interpreted as a function of temperature and NaCl/MgSO4 composition. Quantitative predictions of a polymeric van der Waals attraction associated with the layer collapse in diminishing solvent conditions provides a simple mechanism for explaining the measured fluid-gel data as a dynamic percolation transition. Ultimately, this work identifies the importance of continuum polymeric van der Waals attraction for explaining specific ion-dependent phenomena.     

Cracking in drying latex films

Thin films of latex dispersions containing particles of high glass transition temperature generally crack while drying under ambient conditions. Experiments with particles of varying radii focused on conditions for which capillary stresses normal to the film deform the particles elastically and generate tensile stresses in the plane of the film. Irrespective of the particle size, the drying film contained, simultaneously, domains consisting of a fluid dispersion, a fully dried packing of deformed spheres, and a close packed array saturated with water. Interestingly, films cast from dispersions containing 95-nm sized particles developed tensile stresses and ultimately became transparent even in the absence of water, indicating that van der Waals forces can deform the particles. Employing the stress-strain relation for a drying latex film along with the well-known Griffith's energy balance concept, we calculate the critical stress at cracking and the accompanying crack spacing, in general agreement with the observed values.     

长程范德华力导向作用下胶体凝聚的计算机模拟

采用计算机模拟方法研究了长程范德华力在胶体凝聚过程中的作用, 发现由于胶粒间的范德华力是长程力, 它对胶粒或团簇运动将产生导向作用. 与不考虑导向作用的扩散控制团簇凝聚(DLCA)模型比较, 这种导向作用不仅加速了胶体的凝聚过程, 而且形成了更致密、分形维数更大的结构体. 研究还发现, 长程范德华力导向作用对胶粒的初始浓度非常敏感, 不论是在凝聚物的结构还是凝聚速率方面, 只有在胶粒初始浓度较低时, 该导向作用效应才明显. 其可能的原因是,在胶粒初始浓度较高时, 由于胶粒布朗运动的平均自由程很短而且位阻效应大, 从而使导向作用效应未能反映出来.     

Hydroxamate类抑制剂与MMP-2的分子对接研究

通过分子动力学模拟研究了MMP-2和hydroxamate抑制剂之间的作用模式。在分子动力学模拟中,对于催化区的锌离子和其共价结合的配体(包括抑制剂和组氨酸)采用了键合的模型。从模拟的结果可以看到,R^1取代基和MMP-2的S1疏水口袋中的部分残基能形成很好的几何匹配,从而可以产生很强的范德华和疏水相互作用。模拟结果也表明,两个抑制剂和MMP-2之间分别能形成5个和8个氢键,抑制剂B比A活性更高的原因就是能够形成更加有利氢键作用模式。在整个模拟过程中,催化锌都能保持好的五配位形式,配位键的长度也处于稳定的状态,预测得到的MMP-2和其抑制剂的相互作用模式对于全新抑制剂的设计提供了非常重要的结构信息。     

The bridging force between two plates by polyelectrolyte chains

The interaction between particles in a colloidal system can be significantly affected by their bridging by polyelectrolyte chains. In this paper, the bridging is investigated by using a self-consistent field approach which takes into account the van der Waals interactions between the segments of the polyelectrolyte molecules and the plates, as well as the electrostatic and volume exclusion interactions. A positive contribution to the force between two plates is generated by the van der Waals interactions between the segments and the plates. This positive (repulsive) contribution plays an important role in the force when the distances between the plates are small. With increasing van der Waals interaction strength between segments and plates, the force between the plates becomes more repulsive at small distances and more attractive at large distances. When the surfaces of the plates have a constant surface electrical potential and a charge sign opposite to that of the polyelectrolyte chains, the force between the two plates becomes less attractive as the bulk polyelectrolyte concentration increases. This behavior is due to a higher bulk counterion concentration dissociated from the polyelectrolyte molecules. At short distances, the force between plates is more repulsive for stiffer chains. A comparison between theoretical and experimental results regarding the contraction of the interlayer separation between the platelets of vermiculite clays against the concentration of poly(vinyl methyl ether) was made.     

Anisotropy of Atomic Van der Waals Radii in the Gas-Phase and Condensed Molecules\sp{1}

The anisotropic components of van der Waals radii for 5b, 6b and 7b Group elements were determined from the structures of gaseous van der Waals complexes and of crystalline molecular solids. The transition from the gaseous to the solid state reduces the anisotropy of a van der Waals atomic shape. This anisotropy is largely responsible for the changes of intermolecular distances, which are often misinterpreted as an effect of hydrogen bonds.     

表皮生长因子受体和抑制剂之间分子对接的研究

研究了表皮生长因子受体(EGFR)和4-苯胺喹唑啉类抑制剂之间的相互作用模式，表皮生长因子受体的三维结构通过同源蛋白模建的方法得到，而抑制剂和靶酶结合复合物结构则通过分子力学和分子动力学结合的方法计算得到。从模拟结果得到的抑制剂和靶酶之间的相互作用模式表明范德华相互作用、疏水相互作用以及氢键相互作用对抑制剂的活性都有重要的影响，抑制剂的苯胺部分位于活性口袋的底部，能够与受体残基的非极性侧链产生很强的范德华和疏水相互作用，抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏水性相互作用，而抑制剂喹唑啉环上的氮原子能和周围的残基形成较强的氢键相互作用，对抑制剂的活性有较大的影响，计算得到抑制剂和靶酶之间的非键相互作用能以及抑制剂和靶酶之间的相互作用信息能够很好地解释抑制剂活性和结构的关系，为全新抑制剂的设计提供了重要的结构信息。     

Van der Waals surface graphs and molecular shape

A van der Waals surface graph is the graph defined on a van der Waals surface by the intersections of the atomic van der Waals spheres. A van der Waals shape graph has a vertex for each atom with a visible face on the van der Waals surface, and edges between vertices representing atoms with adjacent faces on the van der Waals surface. These are discrete invariants of three‐dimensional molecular shape. Some basic properties of van der Waals surface graphs are studied, including their relationship with the Voronoi diagram of the atom centres, and a class of molecular embeddings is identified for which the dual of the van der Waals surface graph coincides with the van der Waals shape graph. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modeling competitive guest binding to beta-cyclodextrin molecular printboards

Anchoring of functionalized guest molecules to self-assembled monolayers (SAMs) is key to the development of molecular printboards for nanopatterning. One very promising system involves guest binding to immobilized beta-cyclodextrin (beta-CD) hosts, with guest:host recognition facilitated by a hydrophobic interaction between uncharged anchor groups on the guest molecule and beta-CD hosts self-assembled at gold surfaces. We use molecular dynamics free energy (MDFE) simulations to describe the specificity of guest:beta-CD association. We find good agreement with experimental thermodynamic measurements for binding enthalpy differences between three commonly used phenyl guests: benzene, toluene, and t-butylbenzene. van der Waals interaction with the inside of the host cavity accounts for almost all of the net stabilization of the larger phenyl guests in beta-CD. Partial and full methylation of the secondary rim of beta-CD decreases host rigidity and significantly impairs binding of both phenyl and larger adamantane guest molecules. The beta-CD cavity is also very intolerant of guest charging, penalizing the oxidized state of ferrocene by at least 7 kcal/mol. beta-CD hence expresses moderate specificity toward uncharged organic guest molecules by van der Waals recognition, with a much higher specificity calculated for electrostatic recognition of organometallic guests.     

Ionic Liquids: Dissecting the Enthalpies of Vaporization

We calculate the heats of vaporisation for imidazolium‐based ionic liquids [C_nmim][NTf₂] with n=1, 2, 4, 6, 8 by means of molecular dynamics (MD) simulations and discuss their behavior with respect to temperature and the alkyl chain length. We use a force field developed recently. The different cohesive energies contributing to the overall heats of vaporisations are discussed in detail. With increasing alkyl chain length, the Coulomb contribution to the heat of vaporisation remains constant at around 80 kJ mol^?1, whereas the van der Waals interaction increases continuously. The calculated increase of about 4.7 kJ mol^?1 per CH₂‐group of the van der Waals contribution in the ionic liquid exactly coincides with the increase in the heats of vaporisation for n‐alcohols and n‐alkanes, respectively. The results support the importance of van der Waals interactions even in systems completely composed of ions.     

EGFR和4-苯胺喹唑啉类抑制剂之间相互作用模式的研究

采用分子动力学和MM/PBSA相结合的方法预测了表皮生长因子受体和4-苯胺喹 啉类抑制剂的相互作用模式。在分子动力学采样的基础上,采用MM/PBSA的方法分 别预测了四种可能结合模式下表皮生长因子受体和4-苯胺喹唑啉类抑制剂间的结合 自由能。在MM/PBSA计算中,受体和抑制剂之间的非键相互作用能采用分子力学 (MM)的方法得到;溶剂效应中极性部分对自由能的贡献通过解Possion- Boltzmanne (PB)方程的方法得到;溶液效应中非极性部分对自由能的贡献则通过 分子表面积计算(SA)的方法得到。计算表明,在四种结合模式下,表皮生长因子受 体和4-苯胺喹唑啉类抑制剂之间的结合自由能有较大的差别。在最佳的相互作用模 式中,抑制剂的苯胺部分位于活性口袋的底部,能够与受体残基的非极性侧链产生 很强的范德华和疏水相互作用。抑制剂喹唑啉环上的N(1)原子能够和Met-769上的 NH形成稳定的氢键,而抑制剂上的N(3)原子则和周围的一个水分子形成氢键。同时 ,抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏 水相互作用。最佳结合模式能够很好地解释已有抑制剂结构和活性间的关系。