Biological water: Femtosecond dynamics of macromolecular hydration

The unique features of a macromolecule and water as a solvent make the issue of solvation unconventional, with questions about the static versus dynamic nature of hydration and the, physics of orientational and translational diffusion at the boundary. For proteins, the hydration shell that covers the surface is critical to the stability of its structure and function. Dynamically speaking, the residence time of water at the surface is a signature of its mobility and binding. With femtosecond time resolution it is possible to unravel the shortest residence times which are key for the description of the hydration layer, static or dynamic. In this article we review these issues guided by experimental studies, from this laboratory, of polar hydration dynamics at the surfaces of two proteins (Subtilisin Carlsberg (SC) and Monellin). The natural probe tryptophan amino acid was used for the interrogation of the dynamics, and for direct comparison we also studied the behavior in bulk water - a complete hydration in 1 ps. We develop a theoretical description of solvation and relate the theory to the experimental observations. In this - theoretical approach, we consider the dynamical equilibrium in the hydration shell, defining the rate processes for breaking and making the transient hydrogen bonds, and the effective friction in the layer which is defined by the translational and orientational motions of water molecules. The relationship between the residence time of water molecules and the observed slow component in solvation dynamics is a direct one. For the two proteins studied, we observed a "bimodal decay" for the hydration correlation function, with two primary relaxation times: ultrafast, typically 1 ps or less, and longer, typically 15-40 ps, and both are related to the residence time at the protein surface, depending on the binding energies. We end by making extensions to studies of the denatured state of the protein, random coils, and the biomimetic micelles, and conclude with our thoughts on the relevance of the dynamics of native structures to their functions.

Performance analysis of the random early detection algorithm

In this article we consider a finite queue with its arrivals controlled by the random early detection algorithm. This is one of the most prominent congestion avoidance schemes in the Internet routers. The aggregate arrival stream from the population of transmission control protocol sources is locally considered stationary renewal or Markov modulated Poisson process with general packet length distribution. We study the exact dynamics of this queue and provide the stability and the rates of convergence to the stationary distribution and obtain the packet loss probability and the waiting time distribution. Then we extend these results to a two traffic class case with each arrival stream renewal. However, computing the performance indices for this system becomes computationally prohibitive. Thus, in the latter half of the article, we approximate the dynamics of the average queue length process asymptotically via an ordinary differential equation. We estimate the error term via a diffusion approximation. We use these results to obtain approximate transient and stationary performance of the system. Finally, we provide some computational examples to show the accuracy of these approximations.

Assessment of treatment capabilities of Varthur Lake, Bangalore, India

Manmade waterbodies have traditionally been used for domestic and irrigation purposes. Unplanned urbanization and ad-hoc approaches have led to these waterbodies receiving untreated sewage. This enriches and eutrophies the waterbody. A physicochemical and biological analysis of sewage-fed Varthur Lake in Bangalore was carried out and its treatment capabilities in terms of BOD removal, nutrient assimilation and self-remediation were assessed. Anaerobic conditions (0 mg/L) prevail at the inlet which improves towards the outlets due to algal aeration. This removed > 50% BOD in the monsoon season but was inhibited by floating macrophytes in all other seasons. Alkalinity, TDS, conductivity and hardness values were higher when compared to earlier studies. This study shows the lake behaves as an anaerobic~aerobic lagoon with a residence time of 4.8 d treating the wastewater to a considerable extent. Further research is required to optimise the system performance.

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.

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]

Catalysis of tRNA Aminoacylation: Single Turnover to Steady-State Kinetics of tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRS) catalyze the bimolecular association reaction between amino acid and tRNA by specifically and unerringly choosing the cognate amino acid and tRNA. There are two classes of such synthetases that perform tRNA-aminoacylation reaction. Interestingly, these two classes of aminoacyl-tRNA synthetases differ not only in their structures but they also exhibit remarkably distinct kinetics under pre-steady-state condition. The class I synthetases show initial burst of product formation followed by a slower steady-state rate. This has been argued to represent the influence of slow product release. In contrast, there is no burst in the case of class H enzymes. The tight binding of product with enzyme for class I enzymes is correlated with the enhancement of rate in presence of elongation factor. EF-TU. In spite of extensive experimental studies, there is no detailed theoretical analysis that can provide a quantitative understanding of this important problem. In this article, we present a theoretical investigation of enzyme kinetics for both classes of aminoacyl-tRNA synthetases. We present an augmented kinetic scheme and then employ the methods of time-dependent probability statistics to obtain expressions for the first passage time distribution that gives both the time-dependent and the steady-state rates. The present study quantitatively explains all the above experimental observations. We propose an alternative path way in the case of class II enzymes showing the tRNA-dependent amino acid activation and the discrepancy between the single-turnover and steady-state rate.

Treatment efficacy of algae-based sewage treatment plants

Lagoons have been traditionally used in India for decentralized treatment of domestic sewage. These are cost effective as they depend mainly on natural processes without any external energy inputs. This study focuses on the treatment efficiency of algae-based sewage treatment plant (STP) of 67.65 million liters per day (MLD) capacity considering the characteristics of domestic wastewater (sewage) and functioning of the treatment plant, while attempting to understand the role of algae in the treatment. STP performance was assessed by diurnal as well as periodic investigations of key water quality parameters and algal biota. STP with a residence time of 14.3 days perform moderately, which is evident from the removal of total chemical oxygen demand (COD) (60 %), filterable COD (50 %), total biochemical oxygen demand (BOD) (82 %), and filterable BOD (70 %) as sewage travels from the inlet to the outlet. Furthermore, nitrogen content showed sharp variations with total Kjeldahl nitrogen (TKN) removal of 36 %; ammonium N (NH4-N) removal efficiency of 18 %, nitrate (NO3-N) removal efficiency of 22 %, and nitrite (NO2-N) removal efficiency of 57.8 %. The predominant algae are euglenoides (in facultative lagoons) and chlorophycean members (maturation ponds). The drastic decrease of particulates and suspended matter highlights heterotrophy of euglenoides in removing particulates.

Analysis of an queue operated under the FCFS policy with exact admission control

In this article, we present an exact theoretical analysis of an system, with arbitrary distribution of relative deadline for the end of service, operated under the first come first served scheduling policy with exact admission control. We provide an explicit solution to the functional equation that must be satisfied by the workload distribution, when the system reaches steady state. We use this solution to derive explicit expressions for the loss ratio and the sojourn time distribution. Finally, we compare this loss ratio with that of a similar system operating without admission control, in the cases of some common distributions of the relative deadline.

Influence of surface area to volume ratio of fuel particles on gasification process in a fixed bed

The paper addresses the effect of particle size on tar generation in a fixed bed gasification system. Pyrolysis, a diffusion limited process, depends on the heating rate and the surface area of the particle influencing the release of the volatile fraction leaving behind residual char. The flaming time has been estimated for different biomass samples. It is found that the flaming time for wood flakes is almost one fourth than that of coconut shells for same equivalent diameter fuel samples. The particle density of the coconut shell is more than twice that of wood spheres, and almost four times compared with wood flakes; having a significant influence on the flaming time. The ratio of the particle surface area to that of an equivalent diameter is nearly two times higher for flakes compared with wood pieces. Accounting for the density effect, on normalizing with density of the particle, the flaming rate is double in the case of wood flakes or coconut shells compared with the wood sphere for an equivalent diameter. This is due to increased surface area per unit volume of the particle. Experiments are conducted on estimation of tar content in the raw gas for wood flakes and standard wood pieces. It is observed that the tar level in the raw gas is about 80% higher in the case of wood flakes compared with wood pieces. The analysis suggests that the time for pyrolysis is lower with a higher surface area particle and is subjected to fast pyrolysis process resulting in higher tar fraction with low char yield. Increased residence time with staged air flow has a better control on residence time and lower tar in the raw gas. (C) 2014 International Energy Initiative. Published by Elsevier Inc. All rights reserved.

The interaction of N-trifluoroacetylgalactosamine and its derivatives with winged bean (Psophocarpus tetragonolobus) basic agglutinin reveals differential mechanism of their recognition: a fluorine-19 nuclear magnetic resonance study

Here, we show the binding results of a leguminosae lectin, winged bean basic agglutinin (WBA I) to N-trifluoroacetylgalactosamine (NTFAGalN), methyl-alpha-N-trifluoroacetylgalactosamine (Me alpha NTFAGalN) and methyl-beta-tifluoroacetylgalactosamine (Me beta NTFAGalN) using (19) F NMR spectroscopy. No chemical shift difference between the free and bound states for NTFAGalN and Me beta NTFAGalN, and 0.01-ppm chemical shift change for Me alpha NTFAGalN, demonstrate that the Me alpha NTFAGalN has a sufficiently long residence time on the protein binding site as compared to Me beta NTFAGalN and the free anomers of NTFAGalN. The sugar anomers were found in slow exchange with the binding site of agglutinin. Consequently, we obtained their binding parameters to the protein using line shape analyses. Aforementioned analyses of the activation parameters for the interactions of these saccharides indicate that the binding of alpha and beta anomers of NTFAGalN and Me alpha NTFAGalN is controlled enthalpically, while that of Me beta NTFAGalN is controlled entropically. This asserts the sterically constrained nature of the interaction of the Me beta NTFAGalN with WBA I. These studies thus highlight a significant role of the conformation of the monosaccharide ligands for their recognition by WBA I.

Electron transfer statistics and thermal fluctuations in molecular junctions

We derive analytical expressions for probability distribution function (PDF) for electron transport in a simple model of quantum junction in presence of thermal fluctuations. Our approach is based on the large deviation theory combined with the generating function method. For large number of electrons transferred, the PDF is found to decay exponentially in the tails with different rates due to applied bias. This asymmetry in the PDF is related to the fluctuation theorem. Statistics of fluctuations are analyzed in terms of the Fano factor. Thermal fluctuations play a quantitative role in determining the statistics of electron transfer; they tend to suppress the average current while enhancing the fluctuations in particle transfer. This gives rise to both bunching and antibunching phenomena as determined by the Fano factor. The thermal fluctuations and shot noise compete with each other and determine the net (effective) statistics of particle transfer. Exact analytical expression is obtained for delay time distribution. The optimal values of the delay time between successive electron transfers can be lowered below the corresponding shot noise values by tuning the thermal effects. (C) 2015 AIP Publishing LLC.

An ultrastable conjugate of silver nanoparticles and protein formed through weak interactions

In recent years, silver nanoparticles (AgNPs) have attracted significant attention owing to their unique physicochemical, optical, conductive and antimicrobial properties. One of the properties of AgNPs which is crucial for all applications is their stability. In the present study we unravel a mechanism through which silver nanoparticles are rendered ultrastable in an aqueous solution in complex with the protein ubiquitin (Ubq). This involves a dynamic and reversible association and dissociation of ubiquitin from the surface of AgNP. The exchange occurs at a rate much greater than 25 s(-1) implying a residence time of <40 ms for the protein. The AgNP-Ubq complex remains stable for months due to steric stabilization over a wide pH range compared to unconjugated AgNPs. NMR studies reveal that the protein molecules bind reversibly to AgNP with an approximate dissociation constant of 55 mu M and undergo fast exchange. At pH > 4 the positively charged surface of the protein comes in contact with the citrate capped AgNP surface. Further, NMR relaxation-based experiments suggest that in addition to the dynamic exchange, a conformational rearrangement of the protein takes place upon binding to AgNP. The ultrastability of the AgNP-Ubq complex was found to be useful for its anti-microbial activity, which allowed the recycling of this complex multiple times without the loss of stability. Altogether, the study provides new insights into the mechanism of protein-silver nanoparticle interactions and opens up new avenues for its application in a wide range of systems.

A Simulation Model for Unified Interrill Erosion and Rill Erosion on Hillslopes

A mathematical model was developed for simulating runoff generation and soil erosion on hillslopes. The model is comprised of three modules: one for overland flow, one for soil infiltration, and one for soil erosion including rill erosion and interrill er

Impact of the commercial fishery on the population of bait shrimp (Penaeus spp.) in Biscayne Bay, 1986

Monthly population size of bait shrimp in the Bay was estimated from December 1984 to July 1985. Growth rates for male and female P. duorarum showed that pink shrimp exhibit a mean residence time in the nursery area (Biscayne Bay) of approximately 21 weeks. Monthly mortality rates were determined for each sex of pink shrimp. It was estimated that 23% and 26% of the male and female monthly population size, respectively, was absorbed by both the fishery and ecosystem monthly. Monthly proportion of the standing stock expected to die exclusively through fishing was 6.5% and 6.0% for males and females respectively. Estimates of emigration rates showed that approximately 4.0% of the population was lost from the Bay system each month. This surplus production was about 50% of the average monthly catch by the fleet. Fishing mortality represents only 8 - 9% of the losses to the shrimp population. The biggest source of loss is emigration, suggesting that most shrimp beyond the size at recruitment (to the fishery) are not utilized for food while in the Bay. Thus, it appears that the direct impact of the fishery on the bait shrimp population is relatively small. (PDF contains 46 pages)

Some aspects of the temperature, oxygen and nutrient distributions in Monterey Bay, California. Annual Report, Part 1, 1973

Seasonal variations in temperature, dissolved oxygen, and nutrients in the nearshore areas and in the canyon area of Monterey Bay, California during 1971-1972 were similar~ During upwelling periods, however, water in the nearshore areas was higher in temperature and oxygen and lower in nutrients than water in the canyon area~ This was caused by upwelled water moving north and south of the canyon into counterclockwise and clockwise flow in the northern and southern ends of the bay respectively. The water was heated by insolation and depleted of its nutrients by photosynthesis during this movement. The residence time of water in the nearshore northern and southern bay during upwelling is estimated to be 3 to 8 days, and this fits well into the above circulation pattern and average measured current velocities of 10 to 15 cm/sec~ There is sorne evidence that this circulation pattern and the estimated residence time may be also valid for on-upwelling periods. Upwelling apparently occurred in Monterey Submarine Canyon at rates of 0.4 to 2.9 m/day and was stronger in 1971 than 1972. (PDF contains 107 pages)