Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Systems Analysis of Tank Irrigation: II. Delayed Start and Water Deficit

Tank irrigation systems in the semiarid regions of India are discussed in this paper. To optimize the grain yield of rice, it is essential to start the agricultural operations in the second week of July so that favorable climatic conditions will prevail during flowering and yield formation stages. Because of low inflow during the initial few weeks of the crop season, often farmers are forced to delay planting until sufficient sowing rain and inflow have occurred or to adopt deficit irrigation during this period. The delayed start affects the grain yield, but will lead to an improved irrigation efficiency. A delayed start of agricultural operations with increased irrigation efficiency leads to the energy resources becoming critical during the peak requirement week, particularly those of female labor and animal power. This necessitates augmenting these resources during weeks of their peak use, either by reorganizing the traditional methods of cultivation or by importing from outside the system.

Active and break spells of the Indian summer monsoon

In this paper, we suggest criteria for the identification of active and break events of the Indian summer monsoon on the basis of recently derived high resolution daily gridded rainfall dataset over India (1951-2007). Active and break events are defined as periods during the peak monsoon months of July and August, in which the normalized anomaly of the rainfall over a critical area, called the monsoon core zone exceeds 1 or is less than -1.0 respectively, provided the criterion is satisfied for at least three consecutive days. We elucidate the major features of these events. We consider very briefly the relationship of the intraseasonal fluctuations between these events and the interannual variation of the summer monsoon rainfall. We find that breaks tend to have a longer life-span than active spells.While, almost 80% of the active spells lasted 3-4 days, only 40% of the break spells were of such short duration. A small fraction (9%) of active spells and 32% of break spells lasted for a week or longer. While active events occurred almost every year, not a single break occurred in 26% of the years considered. On an average, there are 7 days of active and break events from July through August. There are no significant trends in either the days of active or break events. We have shown that there is a major difference between weak spells and long intense breaks. While weak spells are characterized by weak moist convective regimes, long intense break events have a heat trough type circulation which is similar to the circulation over the Indian subcontinent before the onset of the monsoon. The space-time evolution of the rainfall composite patterns suggests that the revival from breaks occurs primarily from northward propagations of the convective cloud zone. There are important differences between the spatial patterns of the active/break spells and those characteristic of interannual variation, particularly those associated with the link to ENSO. Hence, the interannual variation of the Indian monsoon cannot be considered as primarily arising from the interannual variation of intraseasonal variation. However, the signature over the eastern equatorial Indian Ocean on intraseasonal time scales is similar to that on the interannual time scales.

Susceptibility of sperm chromatin to acid denaturation in situ: A study in endogenous FSH-deprived adult male bonnet monkeys (Macaca radiata)

Acid denaturation of calf thymus DNA in vitro followed by acridine orange (AO) binding induced a 112% increase in the emission of red, a 58% decrease in green, and a consequential decrease in the ratio of green:red fluorescences from 1.7 to 0.9. This metachromatic property of AO on binding to DNA following acid denaturation was utilized to study the susceptibility of normal and ovine follicle-stimulating hormone (oFSH) actively immunized bonnet monkey spermatozoa voided throughout the year. For analyses, the scattergram generated by the emission of red and green fluorescences by 10,000 AO-bound sperm from each semen sample was divided into 4 quadrant zones representing percentage cells fluorescing high green-low red (Q1), high green-high red (Q2), low green-low red (Q3) and low green-high red. (Q4). Normal monkey sperm obtained during the months of July-December exhibited 76, 13, and 11% cells in Q2, Q3, and Q4 quadrants, respectively. However, during January-June, when the females of the species are markedly subfertile, noncycling, and amenorrhoeic, the spermatozoa ejaculated by the male monkeys exhibited 38, 39, and 23% sperm in Q2, Q3, and Q4, respectively, the differences being highly significant (p < .01-.001). FSH deprivation induced significant shifts in fluorescence emissions, from respective controls, with 39, 33, and 28% cells in Q2, Q3, and Q4, respectively, during July-December, and 15, 48, and 37% sperm in Q2, Q3, and Q4 quadrants, respectively, during January-June. It is postulated that the altered kinetics of germ cell transformations and the deficient spermiogenesis observed earlier following FSH deprivation in these monkeys may have induced the enhanced susceptibility to acid denaturation in sperm.

Direct numerical simulation of turbulent spots

A group of high-order finite-difference schemes for incompressible flow was implemented to simulate the evolution of turbulent spots in channel flows. The long-time accuracy of these schemes was tested by comparing the evolution of small disturbances to a plane channel flow against the growth rate predicted by linear theory. When the perturbation is the unstable eigenfunction at a Reynolds number of 7500, the solution grows only if there are a comparatively large number of (equispaced) grid points across the channel. Fifth-order upwind biasing of convection terms is found to be worse than second-order central differencing. But, for a decaying mode at a Reynolds number of 1000, about a fourth of the points suffice to obtain the correct decay rate. We show that this is due to the comparatively high gradients in the unstable eigenfunction near the walls. So, high-wave-number dissipation of the high-order upwind biasing degrades the solution especially. But for a well-resolved calculation, the weak dissipation does not degrade solutions even over the very long times (O(100)) computed in these tests. Some new solutions of spot evolution in Couette flows with pressure gradients are presented. The approach to self-similarity at long times can be seen readily in contour plots.

Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

Super-Resolution Reconstruction in Frequency-Domain Optical-Coherence Tomography Using the Finite-Rate-of-Innovation Principle

The standard approach to signal reconstruction in frequency-domain optical-coherence tomography (FDOCT) is to apply the inverse Fourier transform to the measurements. This technique offers limited resolution (due to Heisenberg's uncertainty principle). We propose a new super-resolution reconstruction method based on a parametric representation. We consider multilayer specimens, wherein each layer has a constant refractive index and show that the backscattered signal from such a specimen fits accurately in to the framework of finite-rate-of-innovation (FRI) signal model and is represented by a finite number of free parameters. We deploy the high-resolution Prony method and show that high-quality, super-resolved reconstruction is possible with fewer measurements (about one-fourth of the number required for the standard Fourier technique). To further improve robustness to noise in practical scenarios, we take advantage of an iterated singular-value decomposition algorithm (Cadzow denoiser). We present results of Monte Carlo analyses, and assess statistical efficiency of the reconstruction techniques by comparing their performance against the Cramer-Rao bound. Reconstruction results on experimental data obtained from technical as well as biological specimens show a distinct improvement in resolution and signal-to-reconstruction noise offered by the proposed method in comparison with the standard approach.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

On the Maximum Cloud Zone and the ITCZ over Indian, Longitudes during the Southwest Monsoon

An investigation is presented of the daily variation of the maximum cloud zone (MCZ) and the 7W mb trough in the Northern Hemisphere over the Indian longitudes 70–90°E during April–October for 1973–77. It is found that during June–September there are two favorable locations for a MCZ over these longitudes–on a majority of days the MCZ is present in the monsoon zone north of 15°N, and often a secondary MCZ occurs in the equatorial region (0–10°N). The monsoon MCZ gets established by northward movement of the MCZ occurring over the equatorial Indian ocean in April and May. The secondary MCZ appears intermittently, and is characterized by long spells of persistence only when the monsoon MCZ is absent. In each of the seasons studied, the MCZ temporarily disappeared from the mean summer monsoon location (15–28°N) about four weeks after it was established near the beginning of July. It is reestablished by the northward movement of the secondary MCZ, which becomes active during the absence of the monsoon MCZ, in a manner strikingly similar to that observed in the spring to summer transition. A break in monsoon conditions prevails just prior to the temporary disappearance of the monsoon MCZ. Thus we conclude that the monsoon MCZ cannot survive for longer than a month without reestablishment by the secondary MCZ. Possible underlying mechanisms are also discussed.

Spectroscopic Observation of Oscillations in the Corona During the Total Solar Eclipse of 22 July 2009

We performed high resolution spectroscopy of the solar corona during the total solar eclipse of 22 July 2009 in two emission lines: the green line at 5303 due to Fe xiv and the red line at 6374 due to Fe x, simultaneously from Anji (latitude 30A degrees 28.1' N; longitude 119A degrees 35.4' E; elevation 890 m), China. A two-mirror coelostat with 100 cm focal length lens produced a 9.2 mm image of the Sun. The spectrograph using 140 cm focal length lens in Littrow mode and a grating with 600 lines per millimeter blazed at 2 mu m provided a dispersion of 30 m and 43 m per pixel in the fourth order around the green line and third order around the red line, respectively. Two Peltier cooled 1k x 1k CCD cameras, with a pixel size of 13 mu m square and 14-bit readout at 10 MHz operated in frame transfer mode, were used to obtain the time sequence spectra in two emission lines simultaneously. The duration of totality was 341 s, but we could get spectra for 270 s after a trial exposure at an interval of 5 s. We report here on the detection of intensity, velocity, and line width oscillations with periodicity in the range of 25 -50 s. These oscillations can be interpreted in terms of the presence of fast magnetoacoustic waves or torsional Alfv,n waves. The intensity ratios of green to red emission lines indicate the temperature of the corona to be 1.65 MK in the equatorial region and 1.40 MK in the polar region, relatively higher than the expected temperature during the low activity period. The width variation of the emission lines in different coronal structures suggests different physical conditions in different structures.

Invited paper Application of the least mean fourth (LMF) adaptive algorithm to signals associated with gaussian noise

This paper considers the applicability of the least mean fourth (LM F) power gradient adaptation criteria with 'advantage' for signals associated with gaussian noise, the associated noise power estimate not being known. The proposed method, as an adaptive spectral estimator, is found to provide superior performance than the least mean square (LMS) adaptation for the same (or even lower) speed of convergence for signals having sufficiently high signal-to-gaussian noise ratio. The results include comparison of the performance of the LMS-tapped delay line, LMF-tapped delay line, LMS-lattice and LMF-lattice algorithms, with the Burg's block data method as reference. The signals, like sinusoids with noise and stochastic signals like EEG, are considered in this study.

A QUADRATIC C degrees INTERIOR PENALTY METHOD FOR LINEAR FOURTH ORDER BOUNDARY VALUE PROBLEMS WITH BOUNDARY CONDITIONS OF THE CAHN-HILLIARD TYPE

We develop a quadratic C degrees interior penalty method for linear fourth order boundary value problems with essential and natural boundary conditions of the Cahn-Hilliard type. Both a priori and a posteriori error estimates are derived. The performance of the method is illustrated by numerical experiments.