Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Transition from Cassie to impaled state during drop impact on groove-textured solid surfaces

Liquid drops impacted on textured surfaces undergo a transition from the Cassie state characterized by the presence of air pockets inside the roughness valleys below the drop to an impaled state with at least one of the roughness valleys filled with drop liquid. This occurs when the drop impact velocity exceeds a particular value referred to as the critical impact velocity. The present study investigates such a transition process during water drop impact on surfaces textured with unidirectional parallel grooves referred to as groove-textured surfaces. The process of liquid impalement into a groove in the vicinity of drop impact through de-pinning of the three-phase contact line (TPCL) beneath the drop as well as the critical impact velocity were identified experimentally from high speed video recordings of water drop impact on six different groove-textured surfaces made from intrinsically hydrophilic (stainless steel) as well as intrinsically hydrophobic (PDMS and rough aluminum) materials. The surface energy of various 2-D configurations of liquid-vapor interface beneath the drop near the drop impact point was theoretically investigated to identify the locally stable configurations and establish a pathway for the liquid impalement process. A force balance analysis performed on the liquid-vapor interface configuration just prior to TPCL de-pinning provided an expression for the critical drop impact velocity, U-o,U-cr, beyond which the drop state transitions from the Cassie to an impaled state. The theoretical model predicts that Uo, cr increases with the increase in pillar side angle, a, and intrinsic hydrophobicity whereas it decreases with the increase in groove top width, w, of the groove-textured surface. The quantitative predictions of the theoretical model were found to show good agreement with the experimental measurements of U-o,U-cr plotted against the surface texture geometry factor in our model, {tan(alpha/2)/w}(0.5).

Transition in vortex breakdown modes in a coaxial isothermal unconfined swirling jet

This paper reports first observations of transition in recirculation pattern from an open-bubble type axisymmetric vortex breakdown to partially open bubble mode through an intermediate, critical regime of conical sheet formation in an unconfined, co-axial isothermal swirling flow. This time-mean transition is studied for two distinct flow modes which are characterized based on the modified Rossby number (Ro(m)), i.e., Ro(m) <= 1 and Ro(m) > 1. Flow modes with Ro(m) <= 1 are observed to first undergo cone-type breakdown and then to partially open bubble state as the geometric swirl number (S-G) is increased by similar to 20% and similar to 40%, respectively, from the baseline open-bubble state. However, the flow modes with Ro(m) > 1 fail to undergo such sequential transition. This distinct behavior is explained based on the physical significance associated with Ro(m) and the swirl momentum factor (xi). In essence, xi represents the ratio of angular momentum distributed across the flow structure to that distributed from central axis to the edge of the vortex core. It is observed that xi increases by similar to 100% in the critical swirl number band where conical breakdown occurs as compared to its magnitude in the S-G regime where open bubble state is seen. This results from the fact that flow modes with Ro(m) <= 1 are dominated by radial pressure gradient due to swirl/rotational effect when compared to radial pressure deficit arising from entrainment (due to the presence of co-stream). Consequently, the imparted swirl tends to penetrate easily towards the central axis causing it to spread laterally and finally undergo conical sheet breakdown. However, the flow modes with Ro(m) > 1 are dominated by pressure deficit due to entrainment effect. This blocks the radial inward penetration of imparted angular momentum thus preventing the lateral spread of these flow modes. As such these structures fail to undergo cone mode of vortex breakdown which is substantiated by a mere 30%-40% rise in xi in the critical swirl number range. (C) 2014 AIP Publishing LLC.

Studying electrons on curved surfaces by trapping and manipulating multielectron bubbles in liquid helium

Investigations of two-dimensional electron systems (2DES) have been achieved with two model experimental systems, covering two distinct, non-overlapping regimes of the 2DES phase diagram, namely the quantum liquid phase in semiconducting heterostructures and the classical phases observed in electrons confined above the surface of liquid helium. Multielectron bubbles in liquid helium offer an exciting possibility to bridge this gap in the phase diagram, as well as to study the properties of electrons on curved flexible surfaces. However, this approach has been limited because all experimental studies have so far been transient in nature. Here we demonstrate that it is possible to trap and manipulate multielectron bubbles in a conventional Paul trap for several hundreds of milliseconds, enabling reliable measurements of their physical properties and thereby gaining valuable insight to various aspects of curved 2DES that were previously unexplored.

Effect of wettability and surface roughness on ice-adhesion strength of hydrophilic, hydrophobic and superhydrophobic surfaces

The anti-icing properties of hydrophilic, hydrophobic and superhydrophobic surfaces/coatings were evaluated using a custom-built apparatus based on zero-degree cone test method. The ice-adhesion reduction factor (ARF) of these coatings has been evaluated using bare aluminium alloy as a reference. The wettability of the surfaces was evaluated by measuring water contact angle (WCA) and sliding angle. It was found that the ice-adhesion strength (tau) on silicone based hydrophobic surfaces was similar to 43 times lower than compared to bare polished aluminium alloy indicating excellent anti-icing property of these coatings. Superhydrophobic coatings displayed poor anti-icing property in spite of their high water repellence. Field Emission Scanning Electron Microscope reveal that Silicone based hydrophobic coatings exhibited smooth surface whereas the superhydrophobic coatings had a rough surface consisting of microscale bumps and protrusions superimposed with nanospheres. Both surface roughness and surface energy play a major role on the ice-adhesion strength of the coatings. The 3D surface roughness profiles of the coatings also indicated the same trend of roughness. An attempt is made to correlate the observed ice-adhesion strength of different surfaces with their wettability and surface roughness. (C) 2014 Elsevier B.V. All rights reserved.

Impact dynamics of high Weber number drops on chemically modified metallic surfaces

Hydrophobic/superhydrophobic metallic surfaces prepared via chemical treatment are encountered in many industrial scenarios involving the impingement of spray droplets. The effectiveness of such surfaces is understood through the analysis of droplet impact experiments. In the present study, three target surfaces with aluminum (Al-6061) as base material-acid-etched, Octadecyl Trichloro Silane (OTS) coated, and acid-etched plus OTS-coated-were prepared. Experiments on the impact of inertia dominated water drops on these chemically modified aluminum surfaces were carried out with the objective to highlight the effect of chemical treatment on the target surfaces on key sub-processes occurring in drop impact phenomenon. High speed videos of the entire drop impact dynamics were captured at three Weber number (We) conditions representative of high We (We > 200) regime. During the early stages of drop spreading, the drop impact resulted in ejection of secondary droplets from spreading drop front on the etched surfaces resembling prompt splash on rough surfaces whereas no such splashing was observable on untreated aluminum surface. Prominent development of undulations (fingers) were observed at the rim of drop spreading on the etched surfaces; between the etched surfaces the OTS-coated surface showed a subdued development of fingers than the uncoated surface. The impacted drops showed intense receding on OTS-coated surfaces whereas on the etched surface a highly irregular receding, with drop liquid sticking to the surface, was observed. Quantitative analyses were performed to reveal the effect of target surface characteristics on drop impact parameters such as temporal variation of spread factor of drop lamella, temporal variation of average finger length during spreading phase, maximum drop spreading, time taken to attain maximum spreading, sensitivity of maximum spreading to We, number of fingers at maximum spreading, and average receding velocity of drop lamella. Existing models for maximum drop spreading showed reasonably good agreement with the experimental measurements on the target surfaces except the acid-etched surface. (C) 2014 Elsevier B.V. All rights reserved.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Spinning strings and minimal surfaces in AdS(3) with mixed 3-form fluxes

Motivated by the recent proposal for the S-matrix in AdS(3) x S-3 with mixed three form fluxes, we study classical folded string spinning in AdS(3) with both Ramond and Neveu-Schwarz three form fluxes. We solve the equations of motion of these strings and obtain their dispersion relation to the leading order in the Neveu-Schwarz flux b. We show that dispersion relation for the spinning strings with large spin S acquires a term given by -root lambda/2 pi b(2) log(2) S in addition to the usual root lambda/pi log S term where root lambda is proportional to the square of the radius of AdS(3). Using SO(2, 2) transformations and re-parmetrizations we show that these spinning strings can be related to light like Wilson loops in AdS(3) with Neveu-Schwarz flux b. We observe that the logarithmic divergence in the area of the light like Wilson loop is also deformed by precisely the same coefficient of the b(2) log(2) S term in the dispersion relation of the spinning string. This result indicates that the coefficient of b(2) log(2) S has a property similar to the coefficient of the log S term, known as cusp-anomalous dimension, and can possibly be determined to all orders in the coupling lambda using the recent proposal for the S-matrix.

Hydrophobic hydration driven self-assembly of curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. (C) 2014 AIP Publishing LLC.

Transitory second-order reciprocal connection for two surfaces in point contact

This paper investigates the instantaneous spatial higher pair to lower pair substitute-connection which is kinematically equivalent up to acceleration analysis for two smooth surfaces in point contact. The existing first-order equivalent substitute-connection consisting of a Hooke's joint (U-joint) and a spherical joint (S-joint) connected by an additional link is extended up to second-order. A two step procedure is chalked out for achieving this equivalence. First, the existing method is employed for velocity equivalence. In the second step, the two centers of substitution are obtained as a conjugate relationship involving the principal normal curvatures of the surfaces at the contact point and the screw coordinates of the instantaneous screw axis (ISA) of the first-order relative motion. Unlike the classical planar replacement, this particular substitution cannot be done by merely examining the profiles of the contacting surfaces. An illustrative example of a three-link direct-contact mechanism is presented. (C) 2014 Elsevier Ltd. All rights reserved.

Formation of Prestellar Cores via Non-Isothermal Gas Fragmentation

Sheet-like clouds are common in turbulent gas and perhaps form via collisions between turbulent gas flows. Having examined the evolution of an isothermal shocked slab in an earlier contribution, in this work we follow the evolution of a sheet-like cloud confined by (thermal) pressure and gas in it is allowed to cool. The extant purpose of this endeavour is to study the early phases of core-formation. The observed evolution of this cloud supports the conjecture that molecular clouds themselves are three-phase media (comprising viz. a stable cold and warm medium, and a third thermally unstable medium), though it appears, clouds may evolve in this manner irrespective of whether they are gravitationally bound. We report, this sheet fragments initially due to the growth of the thermal instability (TI) and some fragments are elongated, filament-like. Subsequently, relatively large fragments become gravitationally unstable and sub-fragment into smaller cores. The formation of cores appears to be a three stage process: first, growth of the TI leads to rapid fragmentation of the slab; second, relatively small fragments acquire mass via gas-accretion and/or merger and third, sufficiently massive fragments become susceptible to the gravitational instability and sub-fragment to form smaller cores. We investigate typical properties of clumps (and smaller cores) resulting from this fragmentation process. Findings of this work support the suggestion that the weak velocity field usually observed in dense clumps and smaller cores is likely seeded by the growth of dynamic instabilities. Simulations were performed using the smooth particle hydrodynamics algorithm.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

Calcium binding properties of calcium dependent protein kinase 1 (CaCDPK1) from Cicer arietinum

Calcium plays a crucial role as a secondary messenger in all aspects of plant growth, development and survival. Calcium dependent protein kinases (CDPKs) are the major calcium decoders, which couple the changes in calcium level to an appropriate physiological response. The mechanism by which calcium regulates CDPK protein is not well understood. In this study, we investigated the interactions of Ca2+ ions with the CDPK1 isoform of Cicer arietinum (CaCDPK1) using a combination of biophysical tools. CaCDPK1 has four different EF hands as predicted by protein sequence analysis. The fluorescence emission spectrum of CaCDPK1 showed quenching with a 5 nm red shift upon addition of calcium, indicating conformational changes in the tertiary structure. The plot of changes in intensity against calcium concentrations showed a biphasic curve with binding constants of 1.29 mu M and 120 mu M indicating two kinds of binding sites. Isothermal calorimetric (ITC) titration with CaCl2 also showed a biphasic curve with two binding constants of 0.027 mu M and 1.7 mu M. Circular dichroism (CD) spectra showed two prominent peaks at 208 and 222 nm indicating that CaCDPK1 is a alpha-helical rich protein. Calcium binding further increased the alpha-helical content of CaCDPK1 from 75 to 81%. Addition of calcium to CaCDPK1 also increased fluorescence of 8-anilinonaphthalene-1-sulfonic acid (ANS) indicating exposure of hydrophobic surfaces. Thus, on the whole this study provides evidence for calcium induced conformational changes, exposure of hydrophobic surfaces and heterogeneity of EF hands in CaCDPK1. (C) 2015 Elsevier GmbH. All rights reserved.

A new approach to depict bone surfaces in finger imaging using photoacoustic tomography

Imaging the vasculature close around the finger joints is of interest in the field of rheumatology. Locally increased vasculature in the synovial membrane of these joints can be a marker for rheumatoid arthritis. In previous work we showed that part of the photoacoustically induced ultrasound from the epidermis reflects on the bone surface within the finger. These reflected signals could be wrongly interpreted as new photoacoustic sources. In this work we show that a conventional ultrasound reconstruction algorithm, that considers the skin as a collection of ultrasound transmitters and the PA tomography probe as the detector array, can be used to delineate bone surfaces of a finger. This can in the future assist in the localization of the joint gaps. This can provide us with a landmark to localize the region of the inflamed synovial membrane. We test the approach on finger mimicking phantoms.

Recent advances in engineering topography mediated antibacterial surfaces

The tendency of bacterial cells to adhere and colonize a material surface leading to biofilm formation is a fundamental challenge underlying many different applications including microbial infections associated with biomedical devices and products. Although, bacterial attachment to surfaces has been extensively studied in the past, the effect of surface topography on bacteria-material interactions has received little attention until more recently. We review the recent progress in surface topography based approaches for engineering antibacterial surfaces. Biomimicry of antibacterial surfaces in nature is a popular strategy. Whereas earlier endeavors in the field aimed at minimizing cell attachment, more recent efforts have focused on developing bactericidal surfaces. However, not all such topography mediated bactericidal surfaces are necessarily cytocompatible thus underscoring the need for continued efforts for research in this area for developing antibacterial and yet cytocompatible surfaces for use in implantable biomedical applications. This mini-review provides a brief overview of the current strategies and challenges in the emerging field of topography mediated antibacterial surfaces.