Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.

Simulation of Influence of Bilayer Melting on Dynamics and Thermodynamics of Interfacial Water

We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303

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.

An Efficient Probe for Rapid Detection of Cyanide in Water at Parts per Billion Levels and Naked-Eye Detection of Endogenous Cyanide

A new molecular probe based on an oxidized bis-indolyl skeleton has been developed for rapid and sensitive visual detection of cyanide ions in water and also for the detection of endogenously bound cyanide. The probe allows the naked-eye detection of cyanide ions in water with a visual color change from red to yellow ((max)=80nm) with the immediate addition of the probe. It shows high selectivity towards the cyanide ion without any interference from other anions. The detection of cyanide by the probe is ratiometric, thus making the detection quantitative. A Michael-type addition reaction of the probe with the cyanide ion takes place during this chemodosimetric process. In water, the detection limit was found to be at the parts per million level, which improved drastically when a neutral micellar medium was employed, and it showed a parts-per-billion-level detection, which is even 25-fold lower than the permitted limits of cyanide in water. The probe could also efficiently detect the endogenously bound cyanide in cassava (a staple food) with a clear visual color change without requiring any sample pretreatment and/or any special reaction conditions such as pH or temperature. Thus the probe could serve as a practical naked-eye probe for in-field experiments without requiring any sophisticated instruments.

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.

Structure, dynamics and thermodynamics of single-file water under confinement: effects of polarizability of water molecules

Various structural, dynamic and thermodynamic properties of water molecules confined in single-wall carbon nanotubes (CNTs) are investigated using both polarizable and non-polarizable water models. The inclusion of polarizability quantitatively affects the nature of hydrogen bonding, which governs many properties of confined water molecules. Polarizable water leads to tighter hydrogen bonding and makes the distance between neighboring water molecules shorter than that for non-polarizable water. Stronger hydrogen bonding also decreases the rotational entropy and makes the diffusion constant smaller than in TIP3P and TIP3PM water models. The reorientational dynamics of the water molecules is governed by a jump mechanism, the barrier for the jump being highest for the polarizable water model. Our results highlight the role of polarizability in governing the dynamics of confined water and demonstrate that the inclusion of polarizability is necessary to obtain agreement with the results of ab initio simulations for the distributions of waiting and jump times. The SPC/E water model is found to predict various water properties in close agreement with the results of polarizable water models with much lower computational costs.

Molecular mechanism of water permeation in a helium impermeable graphene and graphene oxide membrane

Layers of graphene oxide (GO) are found to be good for the permeation of water but not for helium (Science, 2012, 335(6067), 442-444) suggesting that the GO layers are dynamic in the formation of a permeation route depending on the environment they are in (i.e., water or helium). To probe the microscopic origin of this observation we calculate the potential of mean force (PMF) of GO sheets (with oxidized and reduced parts), with the inter-planar distance as a reaction coordinate in helium and water. Our PMF calculation shows that the equilibrium interlayer distance between the oxidized part of the GO sheets in helium is at 4.8 angstrom leaving no space for helium permeation. In contrast, the PMF of the oxidized part of the GO in water shows two minima, one at 4.8 angstrom and another at 6.8 angstrom, corresponding to no water and a water filled region, thus giving rise to a permeation path. The increased electrostatic interaction between water with the oxidized part of the sheet helps the sheet open up and pushes water inside. Based on the entropy calculations for water trapped between graphene sheets and oxidized graphene sheets at different inter-sheet spacings, we also show the thermodynamics of filling.

The physicochemical characteristics and anaerobic degradability of desiccated coconut industry waste water

Desiccated coconut industries (DCI) create various intermediates from fresh coconut kernel for cosmetic, pharmaceutical and food industries. The mechanized and non-mechanized DCI process between 10,000 and 100,000 nuts/day to discharge 6-150 m(3) of malodorous waste water leading to a discharge of 2646642 kg chemical oxygen demand (COD) daily. In these units, three main types of waste water streams are coconut kernel water, kernel wash water and virgin oil waste water. The effluent streams contain lipids (1-55 g/l), suspended solids (6-80 g/l) and volatile fatty acids (VFA) at concentrations that are inhibitory to anaerobic bacteria. Coconut water contributes to 20-50 % of the total volume and 50-60 % of the total organic loads and causes higher inhibition of anaerobic bacteria with an initial lag phase of 30 days. The lagooning method of treatment widely adopted failed to appreciably treat the waste water and often led to the accumulation of volatile fatty acids (propionic acid) along with long-chain unsaturated free fatty acids. Biogas generation during biological methane potential (BMP) assay required a 15-day adaptation time, and gas production occurred at low concentrations of coconut water while the other two streams did not appear to be inhibitory. The anaerobic bacteria can mineralize coconut lipids at concentrations of 175 mg/l; however; they are severely inhibited at a lipid level of = 350 mg/g bacterial inoculum. The modified Gompertz model showed a good fit with the BMP data with a simple sigmoid pattern. However, it failed to fit experimental BMP data either possessing a longer lag phase and/or diauxic biogas production suggesting inhibition of anaerobic bacteria.

The interaction between poly(vinylpyrrolidone) and reversed micelles of water/AOT/n-heptane

The interactions between poly(vinylpyrrolidone) (PVP) and the reversed micelles composed of water, AOT, and n-heptane are investigated with the aid of phase diagram, measurements of conductivity and viscosity, Fourier transform infrared (FTIR) spectrum, and dynamic light scattering (DLS). The phase diagrams of water/AOT/heptane in the presence of and absence of PVP are given. The conductivity of the water/AOT/heptane reversed micelle without PVP initially increases and then decreases with the increase of water content, ω0 (the molar ratio of water to AOT), while the plots of conductivity (K) versus ω0 of the reversed micelle in the presence of PVP depend on the PVP concentrations. The plot of K versus ω0 with 2.0%wt PVP is similar to that without PVP. Only the ω0,max (the water content that the maximum conductivity corresponds to) is larger than that without PVP. Nevertheless, the conductivity of the reversed micelle containing more than 4%wt PVP always rises with the increase of the water content in the measured range. The DLS results indicate that the hydrodynamic radius (Rh) in the presence and absence of PVP rises with the increase of ω0. The plots with PVP and without PVP have almost the same value when ω0<17; and after that, it quickly increases with the increase of ω0. It is interesting to find that there is almost no effect of the PVP concentration on the viscosity and Rh of the reversed micelle at ω0 = 15. The FTIR results suggest that the contents of SO3--bound water and Na+-bound water both decrease with PVP added, while the content of the bulky-like water increases. However, the trapped water in the hydrophobic chain of the surfactant is nearly unaffected by PVP. It is also found from the FTIR that the carbonyl group stretching vibration of AOT is fitted into two sub-peaks, which center at 1740 and 1729 cm-1, corresponding to the trans and cis conformations of AOT, respectively.

Determination and modelling of water quality parameter in Abuja metropolis

This research work involves the determination and modelling of water parameter such as pH, temperature, turbidity, chloride, hardness. The result of the analysis was used as important operating variables to generate a model equation of pH, hardness, temperature, turbidity and chloride. The values obtained from the model equation were compared with those from experiment. On an average bases the values were close. These parameters can be used to monitor the extent of pollution of pond water and to monitor stress and diseases of fish. The experimental data of pH was in the range of 6.7 to 6.9 while the modelled result was also between 6.7 to 7.0. The turbidity experimental value was close to the modelled value also. The chloride value for the experimental data was in the range of 25.32 to 35.0. The total hardness value ranges between 4.5 to 65.1 mg/l while the modelled result ranges between 11.025 to 68.402 mg/l. The result was within the acceptable limit of world health organization standard on water quality parameter.

Oxidation of Volatile Organic Compounds in Aqueous Solution and at the Air-water Interface of Aqueous Microdroplets

Isoprene (ISO),the most abundant non-methane VOC, is the major contributor to secondary organic aerosols (SOA) formation. The mechanisms involved in such transformation, however, are not fully understood. Current mechanisms, which are based on the oxidation of ISO in the gas-phase, underestimate SOA yields. The heightened awareness that ISO is only partially processed in the gas-phase has turned attention to heterogeneous processes as alternative pathways toward SOA.

During my research project, I investigated the photochemical oxidation of isoprene in bulk water. Below, I will report on the λ > 305 nm photolysis of H₂O₂ in dilute ISO solutions. This process yields C₁₀H₁₅OH species as primary products, whose formation both requires and is inhibited by O₂. Several isomers of C₁₀H₁₅OH were resolved by reverse-phase high-performance liquid chromatography and detected as MH⁺ (m/z = 153) and MH⁺-18 (m/z = 135) signals by electrospray ionization mass spectrometry. This finding is consistent with the addition of ·OH to ISO, followed by HO-ISO· reactions with ISO (in competition with O₂) leading to second generation HO(ISO)₂· radicals that terminate as C₁₀H₁₅OH via β-H abstraction by O₂.

It is not generally realized that chemistry on the surface of water cannot be deduced, extrapolated or translated to those in bulk gas and liquid phases. The water density drops a thousand-fold within a few Angstroms through the gas-liquid interfacial region and therefore hydrophobic VOCs such as ISO will likely remain in these relatively 'dry' interfacial water layers rather than proceed into bulk water. In previous experiments from our laboratory, it was found that gas-phase olefins can be protonated on the surface of pH < 4 water. This phenomenon increases the residence time of gases at the interface, an event that makes them increasingly susceptible to interaction with gaseous atmospheric oxidants such as ozone and hydroxyl radicals.

In order to test this hypothesis, I carried out experiments in which ISO(g) collides with the surface of aqueous microdroplets of various compositions. Herein I report that ISO(g) is oxidized into soluble species via Fenton chemistry on the surface of aqueous Fe(II)Cl₂ solutions simultaneously exposed to H₂O₂(g). Monomer and oligomeric species (ISO)1-8H⁺ were detected via online electrospray ionization mass spectrometry (ESI-MS) on the surface of pH ~ 2 water, and were then oxidized into a suite of products whose combined yields exceed ~ 5% of (ISO)1-8H⁺. MS/MS analysis revealed that products mainly consisted of alcohols, ketones, epoxides and acids. Our experiments demonstrated that olefins in ambient air may be oxidized upon impact on the surface of Fe-containing aqueous acidic media, such as those of typical to tropospheric aerosols.

Related experiments involving the reaction of ISO(g) with ·OH radicals from the photolysis of dissolved H₂O₂ were also carried out to test the surface oxidation of ISO(g) by photolyzing H₂O₂(aq) at 266 nm at various pH. The products were analyzed via online electrospray ionization mass spectrometry. Similar to our Fenton experiments, we detected (ISO)1-7H⁺ at pH < 4, and new m/z⁺ = 271 and m/z^- = 76 products at pH > 5.

Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

Salmon and freshwater fisheries in the Solway Firth. Report of the Working Party of North West Water Authority and the Annan District Salmon Fishery Board.

The Workinq Party of North West Water Authority and the Annan District Salmon Fishery Board's aims were to identify the current problems relating to the salmon and freshwater fisheries in the Solway Firth and its tributaries (embracing the catchments of the Annan, Kirtle Water, Sark, Esk, Eden, Wampool and Waver - and possibly also the Lochar Water and the Nith) to produce agreed outline solutions, and to consider the financial aspects of the future administration of these fisheries.

Effects of land use and physicochemical water quality on grass shrimp, Palaemonetes pugio, and its parasitic isopod, Probopyrus pandalicola, in South Carolina, USA tidal creeks

Grass shrimp, Palaemonetes pugio, are a common inhabitant of US East and Gulf coast salt marshes and are a food source for recreationally and economically important fish and crustacean species. Due to the relationship of grass shrimp with their ecosystem, any significant changes in grass shrimp population may have the potential to affect the estuarine system. Land use is a crucial concern in coastal areas where increasing development impacts the surrounding estuaries and salt marshes and has made grass shrimp population studies a logical choice to investigate urbanization effects. Any impact on tidal creeks will be an impact on grass shrimp populations and their associated micro-environment whether predator, prey or parasitic symbiont. Anthropogenic stressors introduced into the grass shrimp ecosystem may even change the intensity of infections from parasitic symbionts. An ectoparasite found on P. pugio is the bopyrid isopod Probopyrus pandalicola. Little is known about factors that may affect the occurrence of this isopod in grass shrimp populations. The goal was to analyze the prevalence of P. pandalicola in grass shrimp in relation to land use classifications, water quality parameters, and grass shrimp population metrics. Eight tidal creeks in coastal South Carolina were sampled monthly over a three year period. The occurrence of P. pandalicola ranged from 1.2% to 5.7%. Analysis indicated that greater percent water and marsh coverage resulted in a higher incidence of bopyrid occurrence. Analysis also indicated that higher bopyrid incidence occurred in creeks with higher salinity, temperature, and pH but lower dissolved oxygen. The land use characteristics found to limit bopyrid incidence were limiting to grass shrimp (definitive host) populations and probably copepod (intermediate host) populations as well.