Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Covalently Linked, Water-Dispersible, Cyciodextrin: Reduced-Graphene Oxide Sheets

Reduced-graphene oxide (rGO) sheets have been functionalized by covalently linking beta-cyclodextrin (beta CD) cavities to the sheets via an amide linkage. The functionalized beta-CD:rGO sheets, in contrast to rGO, are dispersible over a wide range of pH values (2-13). Zeta potential measurements indicate that there is more than one factor responsible for the dispersibility. We show here that planar aromatic molecules adsorbed on the rGO sheet as well as nonplanar molecules included in the tethered beta-CD cavities have their fluorescence effectively quenched by the beta-CD:rGO sheets. The beta-CD:rGO sheets combine the hydrophobicity associated with rGO along with the hydrophobicity of the cyclodextrin cavities in a single water-dispersible material.

Effect of pressure on the ionic conductivity of Li+ and Cl- ions in water

A molecular dynamics simulation study of aqueous solution of LiCl is reported as a function of pressure. Experimental measurements of conductivity of Li+ ion as a function of pressure shows an increase in conductivity with pressure. Our simulations are able to reproduce the observed trend in conductivity. A number of relevant properties have been computed in order to understand the reasons for the increase in conductivity with pressure. These include radial distribution function, void and neck distributions, hydration or coordination numbers, diffusivity, velocity autocorrelation functions, angles between ion-oxygen and dipole of water as well as OH vector, mean residence time for water in the hydration shell, etc. These show that the increase in pressure acts as a structure breaker. The decay of the self part of the intermediate scattering function at small wave number k shows a bi-exponential decay at 1 bar which changes to single exponential decay at higher pressures. The k dependence of the ratio of the self part of the full width at half maximum of the dynamic structure factor to 2Dk(2) exhibits trends which suggest that the void structure of water is playing a role. These support the view that the changes in void and neck distributions in water can account for changes in conductivity or diffusivity of Li+ with pressure. These results can be understood in terms of the levitation effect. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4756909]

Prediction of Ground Water Levels in the Uplands of a Tropical Coastal Riparian Wetland using Artificial Neural Networks

Artificial Neural Networks (ANNs) have been found to be a robust tool to model many non-linear hydrological processes. The present study aims at evaluating the performance of ANN in simulating and predicting ground water levels in the uplands of a tropical coastal riparian wetland. The study involves comparison of two network architectures, Feed Forward Neural Network (FFNN) and Recurrent Neural Network (RNN) trained under five algorithms namely Levenberg Marquardt algorithm, Resilient Back propagation algorithm, BFGS Quasi Newton algorithm, Scaled Conjugate Gradient algorithm, and Fletcher Reeves Conjugate Gradient algorithm by simulating the water levels in a well in the study area. The study is analyzed in two cases-one with four inputs to the networks and two with eight inputs to the networks. The two networks-five algorithms in both the cases are compared to determine the best performing combination that could simulate and predict the process satisfactorily. Ad Hoc (Trial and Error) method is followed in optimizing network structure in all cases. On the whole, it is noticed from the results that the Artificial Neural Networks have simulated and predicted the water levels in the well with fair accuracy. This is evident from low values of Normalized Root Mean Square Error and Relative Root Mean Square Error and high values of Nash-Sutcliffe Efficiency Index and Correlation Coefficient (which are taken as the performance measures to calibrate the networks) calculated after the analysis. On comparison of ground water levels predicted with those at the observation well, FFNN trained with Fletcher Reeves Conjugate Gradient algorithm taken four inputs has outperformed all other combinations.

Short- and long-term temporal variations in salinity and the oxygen, carbon and hydrogen isotopic compositions of the Hooghly Estuary water, India

The Hooghly River estuary provides a unique experimental site to understand the effect of monsoonal river discharge on freshwater and seawater mixing. Water samples collected bi-weekly for a duration of 17 months were analyzed for salinity, delta O-18,delta C-13(DIC), as well as delta D to investigate the differential mixing of freshwater and seawater. The differences in salinity and delta O-18 of samples collected during low and high tides on the same day are strongly correlated suggesting a well mixed water column at our sampling site. Low salinity and depleted delta O-18 during monsoon is consistent with increased river discharge as well as high rainfall. We identified different slopes in a delta O-18 versus salinity plot for the estuary water samples collected during monsoon and non-monsoon seasons. This is driven by composition of the freshwater source which is dominated by rainwater during monsoon and rivers during non-monsoon months. Selected delta D analyses of samples indicate that groundwater contributes significantly to the Hooghly Estuary during low rainfall times of the year. delta C-13(DIC) measured in the water recorded low values towards the end of monsoon indicating low productivity (i.e. increased organic respiration) while progressively increasing delta C-13(DIC) values from October till January as well as during some of the pre-monsoon months can be explained by increasing productivity. Very low delta C-13(DIC) (similar to-20%0) suggests involvement of carbon derived from anaerobic oxidation of organics and/or methane with potential contribution from increased anthropogenic water supply. An estimate of seawater incursion into the Hooghly Estuary at different times of the year is obtained by using salinity data in a two-component mixing model. Presence of seawater was found maximum (31-37%) during February till July and lowest (less than or equal to 6%) from September till November. We notice a temporal offset between Ganges River discharge farther upstream at Farakka and salinity variation at the Hooghly Estuary. We believe that this time lag is a result of the physical distance between Farakka and Kakdweep (our sampling location) and put constraints on the travel time of river water during early monsoon. (c) 2012 Published by Elsevier B.V.

Water treatment potential in relation to phytoplankton dynamics in tropical sewage-fed urban water bodies

Urban water bodies frequently receive untreated sewage and water levels in such water bodies are maintained by daily inputs of sewage. They function as “variable-zone” anaerobic-aerobic lagoons suffering several macrophyte, biotic and abiotic stresses. We have studied two such lakes in Bangalore (Bellandur-360 ha and Varthur-220 ha) to understand whether such an occurrence could be made beneficial (maintaining water levels as well as treatment). Such hypertrophic water body receives sewage at 180-250mg/L and is discharged at 25-80mg/L COD/BOD in different seasons. In an earlier study we reported macrophyte altering the purification function of the water body. In this paper we studied the impact of phytoplankton dynamics and macrophyte cover on the functions such as organic load removal. Algal community analysis, algal biomass, macrophyte cover, water quality, nutrient status was studied seasonally during 2009-2010. Oxygen deficiency and sometimes anoxia, recorded from surface samples resulted in high quantities of NH4+-N (30-40mg/L) and phosphate (0.5-4mg/L)-characteristics of anoxic hypertrophic urban lakes. The productiveness favoured high phytoplanktonic community characterized by small cells (<10μm; Chlorella sp. - highest numbers). The lake could be clearly demarcated into an initial anaerobic zone (40% area), a facultative zone (20%) and an aerobic zone (40%) based on redox values and GIS/bathymetry. During summer the lake is covered by floating macrophytes converting the lake into an anoxic/anaerobic water pool subduing the water purification function as well as aesthetics. When macrophytes are controlled such sewage fed water bodies can be used for treating urban wastewater while also maintaining water sustainability in these semi-arid ecosystems. This paper reports the community dynamics of phytoplankton, their function and competition with macrophytes.

Microscopic dynamics of nanoparticle monolayers at air-water interface

We present results of surface mechanical and particle tracking measurements of nanoparticles trapped at the air-water interface as a function of their areal density. We monitor both the surface pressure (II) and isothermal compression modulus (epsilon) as well as the dynamics of nanoparticle clusters, using fluorescence confocal microscopy while they are compressed to very high density near the two dimensional close packing density Phi similar to 0.82. We observe non-monotonic variation in both epsilon and the dynamic heterogeneity, characterized by the dynamical susceptibility chi(4) with Phi, in such high density monolayers. We provide insight into the underlying nature of such transitions in close packed high density nanoparticle monolayers in terms of the morphology and flexibility of these soft colloidal particles.. We discuss the significance our results in the context of related studies on two dimensional granular or colloidal systems. (C) 2013 Elsevier Inc. All rights reserved.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Kinetic proofreading at single molecular level: aminoacylation of tRNA(Ile) and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA(Ile) which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case, respectively), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8x10(2)).(4.33x10(2)) = 7.8x10(4), (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA(Ile). The present study shows a non-michaelis type dependence of rate of reaction on tRNA(Ile) concentration in case of valine. The overall editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA(Ile) whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low. We find that the presence of additional editing domain in class I editing enzyme makes the kinetic proofreading more efficient through enhanced hydrolysis of wrong product at the editing CP1 domain.

Albumin-mediated incorporation of water-insoluble therapeutics in layer-by-layer assembled thin films and microcapsules

The present study demonstrates a method to deliver hydrophobic drugs by incorporation into thin films and microcapsules fabricated via a layer-by-layer assembly approach. The hydrophobic molecule binding properties of albumin have been exploited for solubilization of a water-insoluble molecule, pyrene (model drug), by preparation of non-covalent conjugates with bovine serum albumin (BSA). Conjugation with BSA renders a highly negative zeta potential to the previously uncharged pyrene which favors the assembly formation by electrostatic interaction with a positively charged polyelectrolyte, chitosan (at acidic pH). The growth of the assembly was followed by monitoring pyrene absorbance with successive layer deposition. The thin film assembly was demonstrated to be capable of releasing its hydrophobic cargo under physiological conditions. We demonstrated the applicability of this approach by encapsulating a water-insoluble drug, curcumin. These assemblies were further loaded with the anti-cancer drug Doxorubicin. Biocompatible calcium carbonate microparticles were used for capsule preparation. The porous nature of the microparticles allows for the pre-encapsulation of therapeutic macromolecules like protein. The fabrication of protein encapsulated stable microcapsules with hydrophobic molecules incorporated into the shell of the microcapsules has been demonstrated. The microcapsules were further capable of loading hydrophilic molecules like Rhodamine B. Thus, using the approach described, a multi-agent carrier for hydrophobic and hydrophilic drugs as well as therapeutic macromolecules can be envisioned.

Purifying water containing both anionic and cationic species using a (Zn, Cu)O, ZnO, and Cobalt Ferrite based multiphase adsorbent system

For the purpose of water purification, novel and low-cost adsorbents which are promising replacements for activated carbon are being actively pursued. However, a single-phase material that adsorbs both cationic and anionic species remains elusive. Hence, a low-cost, multiphase adsorbent bed that purifies water containing both anionic and cationic pollutants is a desirable alternative. We choose anionic (Congo red, Orange G) and cationic (methylene blue, malachite green) dyes as model pollutants. These dyes are chosen since they are widely found in effluents from textile, leather, fishery, and pharmaceutical industries, and their carcinogenic, mutagenic, genotoxic, and cytotoxic impact on mammalian cells is well-established. We show that ZnO, (Zn0.24Cu0.76)O and cobalt ferrite based multiphase fixed adsorbent bed efficiently adsorbs model anionic (Congo red, Orange G) and cationic (methylene blue and malachite green) pollutants, and their complex mixtures. All adsorbent phases are synthesized using room-temperature, high-yield (similar to 96-100%), green chemical processes. The nanoadsorbents are characterized by using X-ray powder diffraction (XRD), scanning electron microscopy (SEM), Brunauer-Emmett-Teller (BET) surface area analysis, and zeta potential measurements. The constituent nanophases are deliberately chosen to be beyond 50 nm, in order to avoid the nanotoxic size regime of oxides. Adsorption characteristics of each of the phases are examined. Isotherm based analysis shows that adsorption is both spontaneous and highly favorable. zeta potential measurements indicate that electrostatic interactions are the primary driving force for the observed adsorption behavior. The isotherms obtained are best described using a composite Langmuir-Freundlich model. Pseudo-first-order, rapid kinetics is observed (with adsorption rate constants as high as 0.1-0.2 min(-1) in some cases). Film diffusion is shown to be the primary mechanism of adsorption.

Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

Remarkable role of positional isomers in the design of sensors for the ratiometric detection of copper and mercury ions in water

Cation sensing properties of the three positional isomers of rhodamine based sensors (1-3) are studied in water. The sensors differ only in the position of pyridine's nitrogen. The chemosensor 1, with pyridine nitrogen at ortho-position, showed a selective colorimetric detection of Cu(II) ions in water, at physiological pH 7.4 and also in medium containing BSA (bovine serum albumin) and blood serum. Notably the compound 2 and 3, with pyridine end located at meta-and para-positions did not show any color change with Cu(II) ions, although both the compounds showed turn-on change both in color and fluorescence with Hg(II) ions specifically. All the probes showed ratiometric changes with the specific metal ions. The changing position of nitrogen also changed the complexation pattern of the sensors with the metal ions. Probe 1 showed 2 : 1 complexation with Cu(II), whereas 2 and 3 showed 1 : 1 complexation with Hg(II) ions. The mechanism investigation showed that the change in color upon addition of metal ions is due to the ring-opening of the spirolactam ring of the probes. Cu(II) interacted with ligand 1 through a three-point interaction mode comprising carbonyl oxygen, amido nitrogen and pyridine nitrogen end. But in case of 2 and 3, Hg2+ only interacted through pyridine nitrogen ends. Quantitative estimation of Cu2+ and Hg2+ in complex biological media such as bovine albumin protein (BSA) and human blood serum were performed using these sensors. Rapid on-site detection as well as discrimination of these toxic ions was demonstrated using easily prepared portable test-strips.

Selective and Efficient Detection of Nitro-Aromatic Explosives in Multiple Media including Water, Micelles, Organogel, and Solid Support

Selective detection of nitro-aromatic compounds (NACs) at nanomolar concentration is achieved for the first time in multiple media including water, micelles or in organogels as well as using test strips. Mechanism of interaction of NACs with highly fluorescent p-phenylenevinylene-based molecules has been described as the electron transfer phenomenon from the electron-rich chromophoric probe to the electron deficient NACs. The selectivity in sensing is guided by the pK(a) of the probes as well as the NACs under consideration. TNP-induced selective gel-to-sol transition in THF medium is also observed through the reorganization of molecular self-assembly and the portable test trips are made successfully for rapid on-site detection purpose.

Modeling Equitable Distribution of Water: Dynamic Inversion-Based Controller Approach

It is a well-known fact that most of the developing countries have intermittent water supply and the quantity of water supplied from the source is also not distributed equitably among the consumers. Aged pipelines, pump failures, and improper management of water resources are some of the main reasons for it. This study presents the application of a nonlinear control technique to overcome this problem in different zones in the city of Bangalore. The water is pumped to the city from a large distance of approximately 100km over a very high elevation of approximately 400m. The city has large undulating terrain among different zones, which leads to unequal distribution of water. The Bangalore, inflow water-distribution system (WDS) has been modeled. A dynamic inversion (DI) nonlinear controller with proportional integral derivative (PID) features (DI-PID) is used for valve throttling to achieve the target flows to different zones of the city. This novel approach of equitable water distribution using DI-PID controllers that can be used as a decision support system is discussed in this paper.