Abstract

The Independent Component Analysis with Reference (ICA - R) also called as constrained ICA (cICA) extracts only the desired source signals from the mixture of source signals by incorporating some prior information into the separation process. To overcome the problem of designing the reference signal when there is no prior information about the desired signal in the cICA, an improved method is proposed to use a different speech signal generated by the same physical source. The cICA is extended to use Bessel coefficients of the observed signals and the reference signal for processing as they converge faster than the other transformations. Since the Bessel functions provide the desired properties, efficient in representing speech signals, less memory storage they have been exploited in speech processing [1]. The results demonstrate the efficiency of the proposed method

Abstract

In many circumstances in image processing, image resizing is done, either to magnify or to reduce the size of a digital image. The spatial domain based resizing methods such as bilinear interpolation and bi-cubic interpolation are simple and work better for image size magnification. The main drawback in using them is that they are not suitable for image size reduction. In this paper, we propose a new method for image resizing based on Bessel transform (BT). The performance of image resizing based on BT is compared to that of spatial domain based resizing techniques, the results are viewed in terms of Peak Signal to Noise ratio (PSNR) and Mean Square Error (MSE). Experimental results confirm that the proposed method maintains better image quality when image size is enlarged and also when image size is reduced.

Abstract

This paper describes the challenges faced in language learning and use, with a focus on domain-specific vocabulary learning. We focus on the architecture and the presentation model of domain-specific knowledge components. The paper first discusses the significance of assistive mechanisms for vocabulary learning in an academic domain and then outlines two architectures for the presentation model. We suggest how the design of this environment for vocabulary use and learning in the academic domain must adopt different presentations and search strategies to confront the specific challenges in the domain. As this problem can only be solved on a long term basis, we also briefly describe a set of games as a first step in the implementation of the assistive mechanisms.

Abstract

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at theMP2/aug-ccpVDZ level of theory. The geometries for each of the base pairs were characterized in terms of Hbonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from-24 to-49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98%of the total interaction energy values.Onthe basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

Abstract

Riboswitches, typically an integral part of several mRNAs, consist of two domains – a ‘cognate metabolite’ sensing aptamer domain and an expression platform – regulating gene expression at the transcription or in the translation level. Because of their small size and high discriminative ability, p urine riboswitches have been studied extensively both experimentally as well as computationally. Crystal structures of ligand bound aptamer domain of purine riboswitches (Fig. 1) reveal a compact three-way junction architecture consisting of two stem-loops (P2-L2 and P3-L3), held together in parallel orientation by tertiary base pairing interactions between the loops (L2 and L3), converging to a base stem (P1) which communicates with the expression platform. Nucleotides belonging to the three junction loops (J1-2, J2-3 and J3-1) constitute the binding pocket which is lined above and below with double layers of stacked base triples. Though crystal structures of corresponding ligand-free aptamer domains are not available, biophysical and computational studies, including our earlier MD simulations of the add adenine riboswitch1, have provided insights into the functional dynamics connecting the two mechanistic events: discriminative ligand binding and communication with the expression Fig. 1 Ligand bound add platform. Adenine riboswitch In order to gain better insights into the ligand binding induced conformational changes, at the molecular level, and an understanding of the general principles underlying the switching event, we have carried out all-atom explicit solvent molecular dynamics simulations of add adenine riboswitch and xpt-pbuX guanine riboswitch, in ligand-free (APO) state and ligand-bound (HOLO) state, using ff99parmbsc0. Detailed analysis, including principal component analysis and conformational clustering, of the trajectories, reveal interesting insights into changes in conformational dynamics accompanying ligand binding. While confirming the existence of a preformed and ‘binding ready’ APO form and ligand binding induced stabilization of the expression platform linking P1 helices in both aptamers, our studies indicate possible role of stronger L2-L3 interactions in giving rise to fine differences in the functional dynamics of the guanine riboswitch, vis-a-vis the adenine riboswitch.

Abstract

Unlike structures of double stranded DNA, stabilized by canonical base pairings, single stranded RNA molecules are known to fold into complex functional 3D structures. Several of these structures provide sites for binding with cognate molecules that induce conformational changes related to their respective functionalities. Their structures and functional dynamics are scripted by a multitude of noncovalent interactions, involving base edges, 2’-OH of ribose, backbone phosphates and metal cations. In the context of the importance of the ribosome as a potential drug target [1], we have studied the geometries and interaction energies associated with hydrogen bonded base pairs and pseudo pairs observed in crystal structures of drug molecules complexed with their target RNAs. Previous studies, as reviewed recently [2], provide insights into the structural aspects of these RNA-ligand interactions. However, the molecular level understanding of their role in RNA functions require the analysis of their stabilities and interaction energies. We have shown earlier how appropriate quantum chemical calculations can be used to assess gas phase interaction energies and intrinsic stabilities of nucleic acid base pairs [3 – 9]. Here we report results of our calculations, at B3LYP/6-31G(d,p) and RIMP2/aug-cc-pVDZ level, for geometry optimization and interaction energy calculation respectively, for base pairs and pseudobase pairs observed in experimental structures of distinct classes of drug molecules such as aminoglycosides, oxazolidinones, tetracyclines, lincosamides etc., bound to RNA. Our studies show that on full geometry optimization, these base pairs tend to achieve greater stability compared to that in their respective experimental geometries. The increase in interaction energies is not due to any increase in the number of hydrogen bonds but due to increase in their strengths because of increased planarity and by approaching in-range donor-acceptor distances of hydrogen bonds. Their gas phase interaction energies range from -11.32 to -29.70 kcal/mol and compare well with those of standard Watson Crick base pairs. Except for the Puromycin - Guanine pair involving sugar-sugar interaction, the Hartree Fock component dominates the interaction energies, in all fully optimised geometries, with values ranging from 53.09% to 87.13%. The EDFT-D energies ranging from -9.68 to -66.94 kcal/mol, suggest that the interaction between the RNA-Drug pairs is further stabilised due to dispersion. We also report detailed studies on the discriminative ability of the theophylline aptamer with respect to its cognate ligand. Unlike in the context of purine riboswitches observed earlier [6], cooperativity is found to be negligible in multiple base-ligand interactions studied in case of Theophylline interacting with its aptamer. However, the binding affinities of theophylline, caffeine and theobromine with the theophylline aptamer correlates well with the experimental data [10] with high binding affinity for theophylline and least for caffeine

Abstract

Protonation and deprotonation of RNA bases, although energetically unfavorable, are known to play a key role in stabilization of the global structure of the RNA molecules and in regulating their functionalities. In this study an attempt has been made to look at the gas phase energies of the individual protonated and deprotonated bases, and the accompanying changes in their respective geometries and charge distribution. . Geometry optimization and energy calculations were carried out using QM methods at different levels of theory and approximations (B3LYP/6-31G++(2d,2p) and MP2/aug-ccpvdz). As expected, free energy calculations show that protonation/deprotonation at any possible site of a RNA base as a disfavored process characterized by a positive free energy change. It is however important to note that protonation at any exocyclic site or at ring imino nitrogen sites, and deprotonation from any site lead to planar and stable geometry of the corresponding conjugate acids and bases respectively. Protonation of secondary amino nitrogen atoms in the ring, such as N1of G or N3 of U, on the other hand, give significantly deformed structures where covalent bonds are also broken. Comparison of the geometrical parameters and Mulliken charge distribution in the neutral bases, with those of the corresponding charged species, showed interesting variations depending on the identity of the bases as well as on their respective protonation/deprotonation sites. It has been found that protonation at imino nitrogen atoms are not only energetically less costly, geometrical changes and charge redistribution that take place in the protonated species also help the protonated base to take part in non-covalent interactions like hydrogen bonding. Protonation at N7 position of Guanine is found to be most favorable among all possible protonation sites in any RNA base. Similar free energy calculation for nucleosides also showed the same trend and N7 protonation of Guanosine comes out to be the energetically least costly process. When we examined experimentally detected protonated base pairs involving N7 of guanine, and carried out geometry optimization for the two protonation possibilities, it was remarkable to note that it was the partner base-protonated configuration which turned out to be more stable. In fact, geometry optimization of the configuration, containing N7-protonated Guanine, resulted in the transfer of proton from N7 of Guanine to the nearest electronegative atom of the partner base. This may be because of enhanced hydrogen bonding in the alternate configuration. The next best possible protonation site among the nucleic acid bases is N3 of Cytosine, which also has many other biochemical roles. N7 protonation of Guanine, however, leads to variation of several theoretically calculated IR frequencies, some of which may destabilize the glycosidic bond. As the protonated bases are quite stable, they may play greater role in base pairing interactions in RNA.

Abstract

The detailed molecular level understanding of physiochemical interactions responsible for the overall architecture and dynamics of functional RNA molecules embodies a major challenge to the scientific community. RNA 3D structure is currently analyzed in terms of strongly hydrogenbonded secondary structural elements interacting through primarily noncanonical tertiary base pairing interactions. Identification of recurring combinations of such elements as motifs, important in terms of architecture or as anchors for association, enhances our understanding of the RNA tertiary structures. Current approaches in motif mining are constrained by their dependence on identification of strongly hydrogen bonded base pairs. In our work, we have attempted a bottom-up approach towards identifying and characterizing structural features responsible for the complex architecture and dynamics of functional RNA molecules. At the lower end we have studied, noncovalent interactions such as inter-nucleotide hydrogen bonding interactions, water and ion mediated interactions; and at the higher end we have tried to capture the recurrent structural motifs made of several base-pairs and higher order. To study all the important weak and 'other than base base' interactions, within a non-redundant dataset of functional RNA molecules, we have developed a python based automated tool Inter-Nucleotide Contact Annotator in RNA (INCAR). INCAR efficiently manages space and time complexity to capture, categorize and suitably annotate all significant non-covalent interactions, such as hydrogen bonded and ion mediated interactions. Results are with high specificity, without compromising on sensitivity, with respect to currently available benchmarks like BPFIND and FR3D. Key features of our analysis include a) classification of bifurcated geometries, b) Identification of various structural elements such as platforms, bulges, loops, stems etc. c) significance of weak hydrogen bonds d) classification of higher order interactions such as triplets and quartets and e) use of a new automated centroid based approach which doesn't require a prior knowledge on standard base-pair geometries and atoms involving in hydrogen bonding is designed and implemented for detection of glycosidic bond orientation of various standard and non-standard base-pairing geometries (bifurcated, single hydrogen bonded, deformed). We have identified a new structurally important A-A Di-loop and an interesting tertiary interaction in a ribozyme. In addition, INCAR has been applied in several areas such as comparative analysis to investigate the structural conservancy of a class of sequence variants, investigation of motifs in terms of 'other than base-base' contacts and analysis of molecular dynamics simulation trajectories of RNA molecules.

Abstract

Protein domains are minimal structural units that can independently fold and carry out discrete biological functions. Evolutionary divergence amongst proteins not only cause considerable sequence changes of protein domains of similar folds and functions, but can also give rise to remarkable length variations. Rapid and heuristic sequence search algorithms are generally sensitive and effective in recognizing protein domains that are distantly related within large sequence databases, but are not well-suited to identify remote homologues of varying lengths. It is also challenging to distinguish reliable hits from a vast number of putative false positives that could have suboptimal sequence similarities. Here, we present a data-mining approach that provides stage-specific filters in sequence searches to reliably accumulate remote homologues which encourages sampling of length variations albeit no compensation on the validity of hitherto identified distant relationships. Realization of remote homologues with vivid length variations could contribute to better understanding of functional variety within protein domain superfamilies.

Abstract

Advances in sensing and data gathering technologies have resulted in an explosion in the volume of data that is being generated, processed, and archived. In particular, this information overload calls for new methods for querying large spatial datasets, since users are often not interested in merely retrieving a list of all data items satisfying a query, but would, instead, like a more informative \summary" of the retrieved items. An example is the so-called top-k problem, where the goal is to retrieve from a set of n weighted points in IRd the k most signicant points, ranked by their weights, that lie in an orthogonal query box in IRd (rather than get a list of all points lying in the query box). In this paper, ecient and output-sensitive solutions are presented for this problem in two settings. In the rst setting, the k points are reported in arbitrary order and the underlying set can be updated dynamically through insertions and deletions of points. In the second setting, the k points are reported in sorted order of their weights.