Contact angles and their hysteresis as a measure of liquid-solid adhesion

The wetting behavior of a series of aliphatic polyamides was examined. Polyamides and polyethylene were molded against glass to produce smooth surfaces. After cleaning, chemical composition of the surfaces was verified with X-ray photoelectron spectroscopy. Advancing and receding contact angles were measured from small sessile water drops. Contact angles decreased with amide content while contact angle hysteresis increased. Wetting free energies calculated from contact angles were equal to those from dewetting, suggesting that contact angle hysteresis did not arise from surface anomalies, but from hydrogen bonding between water and the amide groups in the polyamide surfaces.     

Adsorption hysteresis for a slit-like pore model

The Frenkel-Halsey-Hill equation is used to describe the adsorption branch of a hysteresis loop upon polylayer adsorption with an H3 loop according to IUPAC nomenclature. The equation for the desorption branch of a hysteresis loop is derived from a combined solution to the equation for the Gibbs potential change, given the adsorbent swelling and pore connectivity function, and the Laplace equation taken for the conditions of infinitely elongated meniscus. This equation is shown to connect the adsorbate relative pressure in a bulk phase for the desorption branch of a hysteresis loop with the key parameters of the adsorption system. The equation obtained was verified by a water adsorption isotherm on natural mineral schungite.     

Contact angle hysteresis at the nanoscale: a molecular dynamics simulation study

In this paper, the contact angle hysteresis (CAH) of nanodroplets on both rigid and flexible substrates with different wettabilities was studied using molecular dynamics (MD) simulations. The critical shear stress (CSS) that determines the motion of the contact line (CL) was investigated. A theoretical correlation between CAH and CSS was proposed. Both CAH and CSS reflect the energy dissipation at the CL of the droplet in response to the exerted force. MD results of CAH are qualitatively consistent with the theoretical model. Simulation results also show that, for the same liquid–solid interactions, CAH on the flexible substrate is larger than that on the rigid substrate. These findings aim to enhance our understanding of the mechanism of the CAH at the nanoscale.     

Capillary condensation in a geometrically and a chemically heterogeneous pore: a molecular simulation study

A computer simulation study has been carried out, using an extended Gibbs ensemble Monte Carlo technique, to examine the influence of so-called geometric and chemical disorder on the thermodynamic behavior of simple fluids confined in porous media. The technique allows the equilibrium coexistence of gas and liquid phases to be calculated in a single run. The phase diagram of Lennard-Jones fluid has been calculated in a perfectly cylindrical pore as a reference. Some disorder is then introduced in the porous material, first by spatially modifying the external potential of the initially cylindrical pore, to imitate the geometric disorder of a more realistic pore (undulation, constrictions, etc.) and second by modulating the amplitude of the same initially cylindrical potential to reproduce the energetic disorder of realistic pores due to chemical variations along it. It is shown that the chemical disorder has a much stronger effect on the phase diagram of the confined fluid. The complete adsorption/desorption isotherms are also calculated to help in understanding the large effects of chemical disorder.     

Contact angle hysteresis, adhesion, and marine biofouling

Adhesive and marine biofouling release properties of coatings containing surface-oriented perfluoroalkyl groups were investigated. These coatings were prepared by cross-linking a copolymer of 1H,1H,2H,2H-heptadecafluorodecyl acrylate and acrylic acid with a copolymer of poly(2-isopropenyl-2-oxazoline) and methyl methacrylate at different molar ratios. The relationships between contact angle, contact angle hysteresis, adhesion, and marine biofouling were studied. Adhesion was determined by peel tests using pressure-sensitive adhesives. The chemical nature of the surfaces was studied by using X-ray photoelectron spectroscopy. Resistance to marine biofouling of an optimized coating was studied by immersion in seawater and compared to previous, less optimized coatings. The adhesive release properties of the coatings did not correlate well with the surface energies of the coatings estimated from the static and advancing contact angles nor with the amount of fluorine present on the surface. The adhesive properties of the surfaces, however, show a correlation with water receding contact angles and contact angle hysteresis (or wetting hysteresis) resulting from surface penetration and surface reconstruction. Coatings having the best release properties had both the highest cross-link density and the lowest contact angle hysteresis. An optimized coating exhibited unprecedented resistance to marine biofouling. Water contact angle hysteresis appears to correlate with marine biofouling resistance.     

Generalized Born model with a simple, robust molecular volume correction

Generalized Born (GB) models provide a computationally efficient means of representing the electrostatic effects of solvent and are widely used, especially in molecular dynamics (MD). A class of particularly fast GB models is based on integration over an interior volume approximated as a pairwise union of atom spheres-effectively, the interior is defined by a van der Waals rather than Lee-Richards molecular surface. The approximation is computationally efficient, but if uncorrected, allows for high dielectric (water) regions smaller than a water molecule between atoms, leading to decreased accuracy. Here, an earlier pairwise GB model is extended by a simple analytic correction term that largely alleviates the problem by correctly describing the solvent-excluded volume of each pair of atoms. The correction term introduces a free energy barrier to the separation of non-bonded atoms. This free energy barrier is seen in explicit solvent and Lee-Richards molecular surface implicit solvent calculations, but has been absent from earlier pairwise GB models. When used in MD, the correction term yields protein hydrogen bond length distributions and polypeptide conformational ensembles that are in better agreement with explicit solvent results than earlier pairwise models. The robustness and simplicity of the correction preserves the efficiency of the pairwise GB models while making them a better approximation to reality.     

Capillary condensation hysteresis in overlapping spherical pores: a monte carlo simulation study

The mechanisms of hysteretic phase transformations in fluids confined to porous bodies depend on the size and shape of pores, as well as their connectivity. We present a Monte Carlo simulation study of capillary condensation and evaporation cycles in the course of Lennard-Jones fluid adsorption in the system of overlapping spherical pores. This model system mimics pore shape and connectivity in some mesoporous materials obtained by templating cubic surfactant mesophases or colloidal crystals. We show different mechanisms of capillary hysteresis depending on the size of the window between the pores. For the system with a small window, the hysteresis cycle is similar to that in a single spherical pore: capillary condensation takes place upon achieving the limit of stability of adsorption film and evaporation is triggered by cavitation. When the window is large enough, the capillary condensation shifts to a pressure higher than that of the isolated pore, and the possibility for the equilibrium mechanism of desorption is revealed. These finding may have important implications for practical problems of assessment of the pore size distributions in mesoporous materials with cagelike pore networks.     

In silico model for P-glycoprotein substrate prediction: insights from molecular dynamics and in vitro studies

P-glycoprotein (P-gp) is a plasma membrane efflux transporter belonging to ATP-binding cassette superfamily, responsible for multidrug resistance in tumor cells. Over-expression of P-gp in cancer cells limits the efficacy of many anticancer drugs. A clear understanding of P-gp substrate binding will be advantageous in early drug discovery process. However, substrate poly-specificity of P-gp is a limiting factor in rational drug design. In this investigation, we report a dynamic trans-membrane model of P-gp that accurately identified the substrate binding residues of known anticancer agents. The study included homology modeling of human P-gp based on the crystal structure of C. elegans P-gp, molecular docking, molecular dynamics analyses and binding free energy calculations. The model was further utilized to speculate substrate propensity of in-house anticancer compounds. The model demonstrated promising results with one anticancer compound (NSC745689). As per our observations, the molecule could be a potential lead for anticancer agents devoid of P-gp mediated multiple drug resistance. The in silico results were further validated experimentally using Caco-2 cell lines studies, where NSC745689 exhibited poor permeability (P _app 1.03 ± 0.16 × 10^?6 cm/s) and low efflux ratio of 0.26.     

Optimal drug cocktail design: methods for targeting molecular ensembles and insights from theoretical model systems

Drug resistance is a significant obstacle in the effective treatment of diseases with rapidly mutating targets, such as AIDS, malaria, and certain forms of cancer. Such targets are remarkably efficient at exploring the space of functional mutants and at evolving to evade drug binding while still maintaining their biological role. To overcome this challenge, drug regimens must be active against potential target variants. Such a goal may be accomplished by one drug molecule that recognizes multiple variants or by a drug "cocktail"--a small collection of drug molecules that collectively binds all desired variants. Ideally, one wants the smallest cocktail possible due to the potential for increased toxicity with each additional drug. Therefore, the task of designing a regimen for multiple target variants can be framed as an optimization problem--find the smallest collection of molecules that together "covers" the relevant target variants. In this work, we formulate and apply this optimization framework to theoretical model target ensembles. These results are analyzed to develop an understanding of how the physical properties of a target ensemble relate to the properties of the optimal cocktail. We focus on electrostatic variation within target ensembles, as it is one important mechanism by which drug resistance is achieved. Using integer programming, we systematically designed optimal cocktails to cover model target ensembles. We found that certain drug molecules covered much larger regions of target space than others, a phenomenon explained by theory grounded in continuum electrostatics. Molecules within optimal cocktails were often dissimilar, such that each drug was responsible for binding variants with a certain electrostatic property in common. On average, the number of molecules in the optimal cocktails correlated with the number of variants, the differences in the variants' electrostatic properties at the binding interface, and the level of binding affinity required. We also treated cases in which a subset of target variants was to be avoided, modeling the common challenge of closely related host molecules that may be implicated in drug toxicity. Such decoys generally increased the size of the required cocktail and more often resulted in infeasible optimizations. Taken together, this work provides practical optimization methods for the design of drug cocktails and a theoretical, physics-based framework through which useful insights can be achieved.     

YinYang atom: a simple combined ab initio quantum mechanical molecular mechanical model

A simple interface is proposed for combined quantum mechanical (QM) molecular mechanical (MM) calculations for the systems where the QM and MM regions are connected through covalent bonds. Within this model, the atom that connects the two regions, called YinYang atom here, serves as an ordinary MM atom to other MM atoms and as a hydrogen-like atom to other QM atoms. Only one new empirical parameter is introduced to adjust the length of the connecting bond and is calibrated with the molecule propanol. This model is tested with the computation of equilibrium geometries and protonation energies for dozens of molecules. Special attention is paid on the influence of MM point charges on optimized geometry and protonation energy, and it is found that it is important to maintain local charge-neutrality in the MM region in order for the accurate calculation of the protonation and deprotonation energies. Overall the simple YinYang atom model yields comparable results to some other QM/MM models.     

Hydrogen thermal desorption spectra: insights from molecular simulation

Thermal desorption spectra of a number of metal-organic frameworks were studied using grand canonical Monte Carlo simulation. Our simulation results are in qualitative agreement with experimental results but also show that great care must be taken when choosing the force field to describe the hydrogen/framework interaction. As the simulations additionally yield the positions and potential energies of the adsorbed molecules it is straightforward to assign the peaks and features in the thermal desorption spectra to specific adsorption sites. We show that the location of the peaks is directly related to the hydrogen-framework interaction which is a complex function of the chemical and topological environment of the pore space, the pore size and the presence of specific interaction sites such as open metal sites. Finally, we demonstrate that an IRMOF-8 sample used to obtain an experimental thermal desorption spectrum must have indeed been catenated as previously suspected. Overall, molecular simulation is a useful tool to complement the interpretation of experimental thermal desorption spectra.     

Study of the advancing and receding contact angles: liquid sorption as a cause of contact angle hysteresis

Two types of experiments were used to study the behavior of both advancing and receding contact angles, namely the dynamic one-cycle contact angle (DOCA) and the dynamic cycling contact angle (DCCA) experiments. For the preliminary study, DOCA measurements of different liquids on different solids were performed using an automated axisymmetric drop shape analysis-profile (ADSA-P). From these experimental results, four patterns of receding contact angle were observed: (1) time-dependent receding contact angle; (2) constant receding contact angle; (3) 'stick/slip'; (4) no receding contact angle. For the purpose of illustration, results from four different solid surfaces are shown. These solids are: FC-732-coated surface; poly(methyl methacrylate/n-butyl methacrylate) [P(MMA/nBMA)]; poly(lactic acid) (DL-PLA); and poly(lactic/glycolic acid) 50/50 (DL-PLGA 50/50). Since most of the surfaces in our studies exhibit time dependence in the receding contact angle, a more extended study was conducted using only FC-732-coated surfaces to better understand the possible causes of decreasing receding contact angle and contact angle hysteresis. Contact angle measurements of 21 liquids from two homologous series (i.e. n-alkanes and 1-alcohols) and octamethylcyclotetrasiloxane (OCMTS) on FC-732-coated surfaces were performed. It is apparent that the contact angle hysteresis decreases with the chain length of the liquid. It was found that the receding contact angle equals the advancing angle when the alkane molecules are infinitely large. These results strongly suggest that the chain length and size of the liquid molecule could contribute to contact angle hysteresis phenomena. Furthermore, DCCA measurements of six liquids from the two homologous series on FC-732-coated surfaces were performed. With these experimental results, one can construe that the time dependence of contact angle hysteresis on relatively smooth and homogeneous surfaces is mainly caused by liquid retention/sorption. The results also suggested that the contact angle hysteresis will eventually approach a steady state, where the rate of liquid retention-evaporation or sorption process would balance out each other. If the existence of contact angle hysteresis can be attributed to liquid sorption/retention, one should only use the advancing contact angles (measured on a dry surface) in conjunction with Young's equation for surface energetic calculations.     

Biomaterial immobilization in nanoporous carbon molecular sieves: influence of solution pH, pore volume, and pore diameter

The adsorption of lysozyme (Lz) onto nanoporous carbon molecular sieves with various pore diameters has been studied at different solution pH values. All the adsorption isotherms have successfully been correlated by the Langmuir equation. The amount of adsorbed Lz depends on the solution pH as well as on the specific pore volume and pore diameter of the adsorbents. The maximum adsorption was observed near the isoelectric point of the Lz (pI approximately 11), suggesting that suppression of electric repulsion between the enzymes plays an important role in the adsorption process. Moreover, the amount adsorbed depends on the pore size and pore volume of the nanoporous carbon adsorbents, indicating that the Lz molecules are adsorbed inside the mesopores. CMK-3-150 shows a larger amount of Lz adsorption as compared to CMK-3. The increased Lz adsorption capacity of CMK-3-150 may be due to the larger pore volume and pore diameter as compared to that of CMK-3. The unaltered structural order of the nanoporous adsorbents after the adsorption has been confirmed by the physicochemical characterization techniques such as XRD and N(2) adsorption. In addition, FT-IR spectroscopic studies confirm that the Lz used in this study is stable even after the adsorption on the nanoporous carbon. These results indicate that nanoporous carbon has superior water stability and thus is a more appropriate adsorbent for biomaterials than nanoporous silica.     

Drops down the hill: theoretical study of limiting contact angles and the hysteresis range on a tilted plate

The limiting inclination angle (slip angle), for which a two-dimensional water drop may be at equilibrium on a chemically heterogeneous surface, is exactly calculated for a variety of cases. The main conclusion is that, in the cases studied, the contact angles at the upper and lower contact line do not always simultaneously equal the receding and advancing contact angles, respectively. On a hydrophobic surface, the lowest contact angle (at the upper contact line) tends to be approximately equal to the receding contact angle, while the highest contact angle (at the lower contact line) may be much lower than the advancing contact angle. For hydrophilic surfaces, the opposite is true. These conclusions imply that the hysteresis range cannot in general be measured by analyzing the shape of a drop on an inclined plane. Also, the limiting inclination angle cannot in general be calculated from the classical equation based only on the advancing and receding contact angles.     

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.     

Free energy of a trans-membrane pore calculated from atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations of a lipid bilayer were performed to calculate the free energy of a trans-membrane pore as a function of its radius. The free energy was calculated as a function of a reaction coordinate using a potential of mean constraint force. The pore radius was then calculated from the reaction coordinate using Monte Carlo particle insertions. The main characteristics of the free energy that comes out of the simulations are a quadratic shape for a radius less than about 0.3 nm, a linear shape for larger radii than this, and a rather abrupt change without local minima or maxima between the two regions. In the outer region, a line tension can be calculated, which is consistent with the experimentally measured values. Further, this line tension can be rationalized and understood in terms of the energetic cost for deforming a part of the lipid bilayer into a hydrophilic pore. The region with small radii can be described and understood in terms of statistical mechanics of density fluctuations. In the region of crossover between a quadratic and linear free energy there was some hysteresis associated with filling and evacuation of the pore with water. The metastable prepore state hypothesized to interpret the experiments was not observed in this region.     

Aqua ions-graphene interfacial and confinement behavior: insights from isobaric-isothermal molecular dynamics

We carry out a systematic microstructural characterization of the solid-fluid interface (SFI) of water and simple metal chloride aqueous solutions in contact with a free-standing plate or with two such plates separated by an interplate distance 0 ≤ h (?) ≤ 30 at ambient conditions via isothermal-isobaric molecular dynamics. With this characterization, we target the interrogation of the system in search for answers to fundamental questions regarding the structure of the "external" and "internal" (confined) SFIs, the effect of the differential hydration behavior among species, and its link to species expulsion from confinement. For water at ambient conditions, we found that the structure of the "external" SFIs is independent of the interplate distance h in the range 0 ≤ h (?) ≤ 30, that is, the absence of wall-mediated correlation effects between "external" and "internal" SFIs, and that for h < 9 ? the slit-pores dewet. Moreover, we observed a selective expulsion of ions caused by the differential hydration between the anion and the cations with a consequent charging of the slit-pore. All these observations were interpreted in terms of the axial profiles for precisely defined order parameters, including tetrahedral configuration, hydrogen bonding, and species coordination numbers.     

A simple, efficient polarizable coarse-grained water model for molecular dynamics simulations

The development of coarse-grained (CG) models that correctly represent the important features of compounds is essential to overcome the limitations in time scale and system size currently encountered in atomistic molecular dynamics simulations. Most approaches reported in the literature model one or several molecules into a single uncharged CG bead. For water, this implicit treatment of the electrostatic interactions, however, fails to mimic important properties, e.g., the dielectric screening. Therefore, a coarse-grained model for water is proposed which treats the electrostatic interactions between clusters of water molecules explicitly. Five water molecules are embedded in a spherical CG bead consisting of two oppositely charged particles which represent a dipole. The bond connecting the two particles in a bead is unconstrained, which makes the model polarizable. Experimental and all-atom simulated data of liquid water at room temperature are used for parametrization of the model. The experimental density and the relative static dielectric permittivity were chosen as primary target properties. The model properties are compared with those obtained from experiment, from clusters of simple-point-charge water molecules of appropriate size in the liquid phase, and for other CG water models if available. The comparison shows that not all atomistic properties can be reproduced by a CG model, so properties of key importance have to be selected when coarse graining is applied. Yet, the CG model reproduces the key characteristics of liquid water while being computationally 1-2 orders of magnitude more efficient than standard fine-grained atomistic water models.