Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

GPU-based normalized cuts for road extraction using satellite imagery

This paper presents a GPU implementation of normalized cuts for road extraction problem using panchromatic satellite imagery. The roads have been extracted in three stages namely pre-processing, image segmentation and post-processing. Initially, the image is pre-processed to improve the tolerance by reducing the clutter (that mostly represents the buildings, vegetation,. and fallow regions). The road regions are then extracted using the normalized cuts algorithm. Normalized cuts algorithm is a graph-based partitioning `approach whose focus lies in extracting the global impression (perceptual grouping) of an image rather than local features. For the segmented image, post-processing is carried out using morphological operations - erosion and dilation. Finally, the road extracted image is overlaid on the original image. Here, a GPGPU (General Purpose Graphical Processing Unit) approach has been adopted to implement the same algorithm on the GPU for fast processing. A performance comparison of this proposed GPU implementation of normalized cuts algorithm with the earlier algorithm (CPU implementation) is presented. From the results, we conclude that the computational improvement in terms of time as the size of image increases for the proposed GPU implementation of normalized cuts. Also, a qualitative and quantitative assessment of the segmentation results has been projected.

Surface Freshwater Storage Variations in the Orinoco Floodplains Using Multi-Satellite Observations

Variations in surface water extent and storage are poorly characterized from regional to global scales. In this study, a multi-satellite approach is proposed to estimate the water stored in the floodplains of the Orinoco Basin at a monthly time-scale using remotely-sensed observations of surface water from the Global Inundation Extent Multi-Satellite (GIEMS) and stages from Envisat radar altimetry. Surface water storage variations over 2003-2007 exhibit large interannual variability and a strong seasonal signal, peaking during summer, and associated with the flood pulse. The volume of surface water storage in the Orinoco Basin was highly correlated with the river discharge at Ciudad Bolivar (R = 0.95), the closest station to the mouth where discharge was estimated, although discharge lagged one month behind storage. The correlation remained high (R = 0.73) after removing seasonal effects. Mean annual variations in surface water volume represented similar to 170 km(3), contributing to similar to 45% of the Gravity Recovery and Climate Experiment (GRACE)-derived total water storage variations and representing similar to 13% of the total volume of water that flowed out of the Orinoco Basin to the Atlantic Ocean.

Genome wide chromatin occupancy of mrhl RNA and its role in gene regulation in mouse spermatogonial cells

Mrhl RNA is a nuclear lncRNA encoded in the mouse genome and negatively regulates Wnt signaling in spermatogonial cells through p68/Ddx5 RNA helicase. Mrhl RNA is present in the chromatin fraction of mouse spermatogonial Gc1-Spg cells and genome wide chromatin occupancy of mrhl RNA by ChOP (Chromatin oligo affinity precipitation) technique identified 1370 statistically significant genomic loci. Among these, genes at 37 genomic loci also showed altered expression pattern upon mrhl RNA down regulation which are referred to as GRPAM (Genes Regulated by Physical Association of Mrhl RNA). p68 interacted with mrhl RNA in chromatin at these GRPAM loci. p68 silencing drastically reduced mrhl RNA occupancy at 27 GRPAM loci and also perturbed the expression of GRPAM suggesting a role for p68 mediated mrhl RNA occupancy in regulating GRPAM expression. Wnt3a ligand treatment of Gc1-Spg cells down regulated mrhl RNA expression and also perturbed expression of these 27 GRPAM genes that included genes regulating Wnt signaling pathway and spermatogenesis, one of them being Sox8, a developmentally important transcription factor. We also identified interacting proteins of mrhl RNA associated chromatin fraction which included Pc4, a chromatin organizer protein and hnRNP A/B and hnRNP A2/B1 which have been shown to be associated with lincRNA-Cox2 function in gene regulation. Our findings in the Gc1-Spg cell line also correlate with the results from analysis of mouse testicular tissue which further highlights the in vivo physiological significance of mrhl RNA in the context of gene regulation during mammalian spermatogenesis.

Demonstration of helicase activity in the nonstructural protein, NSs, of the negative-sense RNA virus, Groundnut bud necrosis virus

The nonstructural protein NSs, encoded by the S RNA of groundnut bud necrosis virus (GBNV) (genus Tospovirus, family Bunyaviridae) has earlier been shown to possess nucleic-acid-stimulated NTPase and 50 a phosphatase activity. ATP hydrolysis is an essential function of a true helicase. Therefore, NSs was tested for DNA helicase activity. The results demonstrated that GBNV NSs possesses bidirectional DNA helicase activity. An alanine mutation in the Walker A motif (K189A rNSs) decreased DNA helicase activity substantially, whereas a mutation in the Walker B motif resulted in a marginal decrease in this activity. The parallel loss of the helicase and ATPase activity in the K189A mutant confirms that NSs acts as a non-canonical DNA helicase. Furthermore, both the wild-type and K189A NSs could function as RNA silencing suppressors, demonstrating that the suppressor activity of NSs is independent of its helicase or ATPase activity. This is the first report of a true helicase from a negative-sense RNA virus.

RNA nucleosides as chiral sensing agents in NMR spectroscopy

The study reports chiral sensing properties of RNA nucleosides. Adenosine, guanosine, uridine and cytidine are used as chiral derivatizing agents to differentiate chiral 1 degrees-amines. A three component protocol has been adopted for complexation of nucleosides and amines. The chiral differentiating ability of nucleosides is examined for different amines based on the H-1 NMR chemical shift differences of diastereomers (Delta delta(R,S)). Enantiomeric differentiation has been observed at multiple chemically distinct proton sites. Adenosine and guanosine exhibit large chiral differentiation (Delta delta(R,S)) due to the presence of a purine ring. The diastereomeric excess (de) measured by using adenosine is in good agreement with the gravimetric values.

Investigation of orographic features influencing spatial distribution of rainfall over the Western Ghats of India using satellite data

Using remotely sensed Tropical Rainfall Measuring Mission (TRMM) 3B42 rainfall and topographic data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Digital Elevation Model (DEM), the impact of oroghraphical aspects such as topography, spatial variability of elevation and altitude of apexes are examined to investigate capacious summer monsoon rainfall over the Western Ghats (WG) of India. TRMM 3B42 v7 rainfall data is validated with Indian Meteorological Department (IMD) gridded rainfall data at 0.5 degrees resolution over the WG. The analysis of spatial pattern of monsoon rainfall with orography of the WG ascertains that the grade of orographic precipitation depends mainly on topography of the mountain barrier followed by steepness of windward side slope and altitude of the mountain. Longer and broader, i.e. cascaded topography, elevated summits and gradually increasing slopes impel the enhancement in precipitation. Comparing topography of various states of the WG, it has been observed that windward side of Karnataka receives intense rainfall in the WG during summer monsoon. It has been observed that the rainfall is enhanced before the peak of the mountain and confined up to the height about 800m over the WG. In addition to this, the spatial distribution of heavy and very heavy rainfall events in the last 14 years has also been explored. Heavy and very heavy rain events on this hilly terrain are categorized with a threshold of precipitation (R) in the range 150>R>120mmday(-1) and exceeding 150mmday(-1) using probability distribution of TRMM 3B42 v7 rainfall. The areas which are prone to heavy precipitation are identified. The study would help policy makers to manage the hazard scenario and, to improve weather predictions on mountainous terrain of the WG.

A carboxy terminal domain of the L protein of rinderpest virus possesses RNA triphosphatase activity - The first enzyme in the viral mRNA capping pathway

The large protein L of negative-sense RNA viruses is a multifunctional protein involved in transcription and replication of genomic RNA. It also possesses enzymatic activities involved in capping and methylation of viral mRNAs. The pathway for mRNA capping followed by the L protein of the viruses in the Morbillivirus genus has not been established, although it has been speculated that these viruses may follow the unconventional capping pathway as has been shown for some viruses of Rhabdoviridae family. We had earlier shown that the large protein L of Rinderpest virus expressed as recombinant L-P complex in insect cells as well as the ribonucleoprotein complex from purified virus possesses RNA triphosphatase (RTPase) and guanylyltransferase activities, in addition to RNA dependent RNA polymerase activity. In the present work, we demonstrate that RTPase as well as nucleoside triphosphatase (NTPase) activities are exhibited by a subdomain of the L protein in the C terminal region (a.a. 1640 1840). The RTPase activity depends absolutely on a divalent cation, either magnesium or manganese. Both the RTPase and NTPase activities of the protein show dual metal specificity. Two mutant proteins having alanine mutations in the glutamic acid residues in motif-A of the RTPase domain did not show RTPase activity, while exhibiting reduced NTPase activity suggesting overlapping active sites for the two enzymatic functions. The RTPase and NTPase activities of the L subdomain resemble those of the Vaccinia capping enzyme D1 and the baculovirus LEF4 proteins. (C) 2015 Elsevier Inc. All rights reserved.

Evidence for N-7 guanine methyl transferase activity encoded within the modular domain of RNA-dependent RNA polymerase L of a Morbillivirus

Post-transcriptional modification of viral mRNA is essential for the translation of viral proteins by cellular translation machinery. Due to the cytoplasmic replication of Paramyxoviruses, the viral-encoded RNA-dependent RNA polymerase (RdRP) is thought to possess all activities required for mRNA capping and methylation. In the present work, using partially purified recombinant RNA polymerase complex of rinderpest virus expressed in insect cells, we demonstrate the in vitro methylation of capped mRNA. Further, we show that a recombinant C-terminal fragment (1717-2183 aa) of L protein is capable of methylating capped mRNA, suggesting that the various post-transcriptional activities of the L protein are located in independently folding domains.

Differential binding of ppGpp and pppGpp to E-coli RNA polymerase: photo-labeling and mass spectral studies

(p) ppGpp, a secondary messenger, is induced under stress and shows pleiotropic response. It binds to RNA polymerase and regulates transcription in Escherichia coli. More than 25 years have passed since the first discovery was made on the direct interaction of ppGpp with E. coli RNA polymerase. Several lines of evidence suggest different modes of ppGpp binding to the enzyme. Earlier cross-linking experiments suggested that the beta-subunit of RNA polymerase is the preferred site for ppGpp, whereas recent crystallographic studies pinpoint the interface of beta'/omega-subunits as the site of action. With an aim to validate the binding domain and to follow whether tetra-and pentaphosphate guanosines have different location on RNA polymerase, this work was initiated. RNA polymerase was photo-labeled with 8-azido-ppGpp/8-azido-pppGpp, and the product was digested with trypsin and subjected to mass spectrometry analysis. We observed three new peptides in the trypsin digest of the RNA polymerase labeled with 8-azido-ppGpp, of which two peptides correspond to the same pocket on beta'-subunit as predicted by X-ray structural analysis, whereas the third peptide was mapped on the beta-subunit. In the case of 8-azido-pppGpp-labeled RNA polymerase, we have found only one cross-linked peptide from the beta'-subunit. However, we were unable to identify any binding site of pppGpp on the beta-subunit. Interestingly, we observed that pppGpp at high concentration competes out ppGpp bound to RNA polymerase more efficiently, whereas ppGpp cannot titrate out pppGpp. The competition between tetraphosphate guanosine and pentaphosphate guanosine for E. coli RNA polymerase was followed by gel-based assay as well as by a new method known as DRaCALA assay.

Micro-satellite angular rate estimation using second-order sliding mode observer

Development of computationally efficient and accurate attitude rate estimation algorithm using low-cost commercially available star sensor arrays and processing unit for micro-satellite mission is presented. Our design reduces the computational load of least square (LS)-based rate estimation method while maintaining the same accuracy compared to other rate estimation approaches. Furthermore, rate estimation accuracy is improved by using recently developed fast and accurate second-order sliding mode observer (SOSMO) scheme. It also gives robust estimation in the presence of modeling uncertainties, unknown disturbances, and measurement noise. Simulation study shows that rate estimation accuracy achieved by our LS-based method is comparable with other methods for a typical commercially available star sensor array. The robustness analysis of SOSMO with respect to measurement noise is also presented in this paper. Simulation test bench for a practical scenario of satellite rate estimation uses moment-of-inertia variation and environmental disturbances affecting a typical micro-satellite at 500km circular orbit. Comparison studies of SOSMO with 1-SMO and pseudo-linear Kalman filter show that satisfactory estimation accuracy is achieved by SOSMO.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Anthropogenic aerosol fraction over the Indian region: model simulations versus multi-satellite data analysis

Anthropogenic aerosols play a crucial role in our environment, climate, and health. Assessment of spatial and temporal variation in anthropogenic aerosols is essential to determine their impact. Aerosols are of natural and anthropogenic origin and together constitute a composite aerosol system. Information about either component needs elimination of the other from the composite aerosol system. In the present work we estimated the anthropogenic aerosol fraction (AF) over the Indian region following two different approaches and inter-compared the estimates. We espouse multi-satellite data analysis and model simulations (using the CHIMERE Chemical transport model) to derive natural aerosol distribution, which was subsequently used to estimate AF over the Indian subcontinent. These two approaches are significantly different from each other. Natural aerosol satellite-derived information was extracted in terms of optical depth while model simulations yielded mass concentration. Anthropogenic aerosol fraction distribution was studied over two periods in 2008: premonsoon (March-May) and winter (November-February) in regard to the known distinct seasonality in aerosol loading and type over the Indian region. Although both techniques have derived the same property, considerable differences were noted in temporal and spatial distribution. Satellite retrieval of AF showed maximum values during the pre-monsoon and summer months while lowest values were observed in winter. On the other hand, model simulations showed the highest concentration of AF in winter and the lowest during pre-monsoon and summer months. Both techniques provided an annual average AF of comparable magnitude (similar to 0.43 +/- 0.06 from the satellite and similar to 0.48 +/- 0.19 from the model). For winter months the model-estimated AF was similar to 0.62 +/- 0.09, significantly higher than that (0.39 +/- 0.05) estimated from the satellite, while during pre-monsoon months satellite-estimated AF was similar to 0.46 +/- 0.06 and the model simulation estimation similar to 0.53 +/- 0.14. Preliminary results from this work indicate that model-simulated results are nearer to the actual variation as compared to satellite estimation in view of general seasonal variation in aerosol concentrations.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

A Novel Approach for Multispectral Satellite Image Classification Based on the Bat Algorithm

Among the multiple advantages and applications of remote sensing, one of the most important uses is to solve the problem of crop classification, i.e., differentiating between various crop types. Satellite images are a reliable source for investigating the temporal changes in crop cultivated areas. In this letter, we propose a novel bat algorithm (BA)-based clustering approach for solving crop type classification problems using a multispectral satellite image. The proposed partitional clustering algorithm is used to extract information in the form of optimal cluster centers from training samples. The extracted cluster centers are then validated on test samples. A real-time multispectral satellite image and one benchmark data set from the University of California, Irvine (UCI) repository are used to demonstrate the robustness of the proposed algorithm. The performance of the BA is compared with two other nature-inspired metaheuristic techniques, namely, genetic algorithm and particle swarm optimization. The performance is also compared with the existing hybrid approach such as the BA with K-means. From the results obtained, it can be concluded that the BA can be successfully applied to solve crop type classification problems.