Abstract

Hoogsteen-sugar base pattern of RNA base pairing is studied using density functional theory. Hydrogen bonding patterns of these base apirs are characterized using NBO analysis and AIM analysis. Correlation between strenghth of base pairing and nature of donor-acceptor combinations is also carried out.

Abstract

Predicting the 3D native conformation of a protein given the amino acid sequence is known as protein structure prediction (PSP)problem and is one of the greatest challenge of computational biology. The current work uses a real valued Genetic Algorithm (GA), a powerful variant of conventional GA to simulate the PSP problem. It consists of three evolutionary operators to manipulate the genes and a fitness function based on a simple force field. The individuals of the population correspond to different conformations of the same polypeptide chain, represented by a string of (phi,psi)torsion angles. The conformations are generated under the constraints of Ramachandran plot along with secondary structure information, which are then screened through a set of knowledge based biophysical filters, viz. persistence length and radius of gyration. Only the structures satisfying the filtering criteria are considered for energy minimization using GA. The algorithm has been validated on a set of known globular protein containing 2-4 secondary structure elements. For all the test proteins the algorithm converges rapidly ad the converged structure shows a backbone RMSD (root mean square deviation) or 3-6A as compared to the native structure.

Abstract

Motivation: Various types of functional RNAs-tRNA, Ribosomes, Introns, Aptamers, Ribozymes Diverse and Dynamic molecules Recent availability of RNA crystal structures has facilitated study and analysis

Abstract

Schistosoma mansoni, a trematode parasite causes schistosomiasis, a tropical disease affecting more than 200 million people. It lives in an aerobic environment and therefore has effective redox mechanisms for its survival as well as to evade reactive oxygen species from its host. Although, the host performs these activities by two different systems: thioredoxin and glutaredoxin separately, the parasite has only one unique multifunctional enzyme, thioredoxin glutathione reductase (TGR) for both the activities which is its main bottleneck. Therefore this emerges as a promising target for schistosomiasis. TGR is a homodimer of 65kDa sub units. each consisting 598 amino acids with a selenocysteine residue in the C-terminus. It shows an unprecedented fusion of two proteins: glutaredoxin (Grx) and thioredoxin reductase (TR). Although many biochemical studies have been reported for TGR, the three dimensional structural information is still missing. Herein, we present some of our efforts towards development of a homology model for S. mansoni TGR (SmTGR) using Modeller9.1. A BLAST search with SmTGR sequence revealed about 60% homology with different mammalian TRs and the initial 100 residues showed close similarity with Grx. Therefore, homology modeling was carried out using both TR and Grx templates. The developed models were dimerized and validated by Procheck analysis. The validated structure was then minimized by MD simulations using heated simulated annealing method.

Abstract

Purine riboswitches constitute a class of highly structured regions within the 5'-untranaslated sequence stretches of mRNAs that specifically bind purines to regulate the gene expression. A unique feature of the ligand binding in these riboswitches is the larger role of hydrogen bonding compared to stacking [1]. In order to study the inherent strength and specificity of molecular recognition involved in these interactions, we have carried out gas phase quantum chemical studies on the metabolite binding pockets of adenine [2] and guanine [2] binding riboswitches, with adenine, guanine and a modified pyrimidine respectively, at the B3LYP/6-31G level of theory. An analysis of optimized structures reveals that for three conserved residues - C74, U47 and U51, the ligand-base hydrogen bonding pattern remains consereved. On the other hand, the U22 residue deviates from its crystal geometary, suggesting that its is held in the crystal geometry through strong backbone and local environmental constraints. We have carried out Morokuma decomposition and interaction energy evaluations as well as NBO and AIM analysis to understand the strength and nature of binding of the individual aptamer residues with their respective purine metabolites. The dominant features of these non-covalent interactions have also been analyazed. Acknowledgements: AM thanks DBT for research grant. PS and SS thank CSIR, New Delhi for JRFs. [1] Gilbert, S. D.; Stoddard, C. D.; Wise, S. J.; Batey, R. T. J. Mol. Biol. 2006, 359, 754-768 [2] Serganov, A.; Yuan, Y.; Pikovskaya, O.; Polonskaia, A.; Malinina, L.; Phan, A. T.; Hobartner, C.; Micura, R.; Breakar, R. R.; Patel, D. J. Chem. Biol. 2004, 11, 1729-1741.

Abstract

The need for computational characterization of DNA structure in the context of DNA bending/binding, B ->A transitions in protein-DNA complexes and also in the context of strand orientation especially in DNA triplexes and quadruplexes have given rise to several fine grained parameters, such as backbone conformations, sugar pucker and chi angle, local strand orientation, base pair and base pair step parameters backbone base inclination etc. With their unpaired bases, extensive non helical regions, non-canonical base pairs and tertiary interactions, RNA molecules host a large variety of recurrent structural and functional motifs and parameters such as Zp, defined in the context of double helical DNA base pairs and base pair steps and effectively used for A and B DNA differentiations, can not be used for computationally driven automated mining of complex RNA motifs. Though several algorithms have been proposed based on backbone conformational parameters and though algorithms for computational detection of base pairs and analysis of base pair parameters, are available, there are no reports of motif mining studies involving parameters which, as with Zp, simultaneously involve backbone conformation as well as base pairing geometry. In this work, we propose a novel set of parameters, Omega Pseudo-Torsions and Omega Distances, which capture backbone geometries at base(i) - base(j) interaction points and which can be used to computationally characterize structural features both in RNA and DNA. Pseudo torsion angles: Omega eta: P(i)-C4'(j)-P(j); Omega theta: P(i+1)-C4'(j)-P(j+1) Omega 1: P(i)-C4'(i)-C4'(j)-P(j+1); Omega 2: P(i+1)-C4'(i)-C4'(j)-P(j) Pseudo bond distance: Omega distance: C4'(i)-C4'(j) Where 'i' and 'j' are bases interacting through hydrogen bonds. We report the results of our benchmarking studies of Omega parameters with standard DNA structure classification parameters and of our preliminary studies on their effectiveness in classifying local geometries of RNA strands in different base pairing (canonical and non-canonical) contexts.

Abstract

DNA structures and their variations are important in several contexts such as chromatin remodeling and transcriptional initiation and regulation etc. Theneed for computational characterization of DNA structures associated with DNA bending/binding, B ? A transitions in protein-DNA complexes, molecular dynamics trajectory analysis and also in the context of strand orientation, especially in DNA triplexes and quadruplexes, have given rise to several fine grained parameters, such as backbone conformations, sugar pucker and ? angle, local strand orientation, base pair and base pair step parameters, backbone base inclination etc. With their unpaired bases, extensive non helical regions, non-canonical base pairs and tertiary interactions, RNA molecules host a large variety of recurrent structural and functional motifs and parameters such as zp, defined in the context of double helical DNA base pairs and base pair steps and effectively used for A and B DNA differentiations, can not be used for computationally driven automated mining of complex RNA motifs. Thus for RNA structures, though several motif mining algorithms based on backbone conformational parameters and involving computational detection and analysis of base pairs and their parameters are available, there are no reports of studies involving parameters which simultaneously involve backbone conformation as well as base pairing geometry. In this work, we propose a novel set of parameters, Omega Torsions and Omega Distances, which captures backbone geometries at base(i) ? base(j) interaction points and which can be used to computationally characterize structural features both in RNA and DNA. Pseudo torsion angles: Omega ?: P(i)-C4'(i)-C4'(j)-P(j); Omega ?: P(i+1)-C4'(i)-C4'(j)-P(j+1) Omega 1: P(i)-C4'(i)-C4'(j)-P(j+1); Omega 2: P(i+1)-C4'(i)-C4'(j)-P(j) Pseudo bond distance: Omega distance: C4'(i)-C4'(j) Where 'i' and 'j' are bases interacting through hydrogen bonds. We report the results of our benchmarking studies of Omega parameters, as applied to protein DNA complexes including e.g. TATA box-TBP complexes and nucleosomal DNA, with standard DNA structure classification parameters, such as zp, and of our preliminary studies on their effectiveness in classifying local geometries of RNA strands in different base pairing contexts. We conclude that omega parameters can potentially lead to the development of a single set of parameters which can be used to characterize both DNA as well as RNA structures.

Abstract

Projection of a Markov process with constant rates of transition to a small number of observable aggregated states can result in complex kinetics with memory. Here, we define the entropy production along a single sequence of aggregated states and show that it obeys detailed and integral fluctuation theorems. More importantly, we prove that projection shifts the distribution of entropy production over the ensemble of paths for a nonequilibrium process toward one characteristic of a system at equilibrium. This statement represents an analog of the second law of thermodynamics for path-dependent entropies and thus a new form of constraint of irreversible systems. Numeric examples are presented to illustrate these ideas.

Abstract

Capillary Non-Perfusion (CNP) is a condition in diabetic retinopathy where blood ceases to flow to certain parts of the retina, potentially leading to blindness. This paper presents a solution for automatically detecting and segmenting CNP regions from fundus fluorescein angiograms (FFAs). CNPs are modelled as valleys, and a novel multiresolution technique for trough-based valley detection is presented. The proposed algorithm has been tested on 40 images and validated against expert-marked ground truth. Obtained results are presented as a receiver operating characteristic (ROC) curve. The area under this curve is 0.842 and the distance of ROC from the ideal point (0,1) is 0.31.

Abstract

In this paper we present a method for fovea localization which does not use organization information of other retinal structures like optic disk and arcades. The main advantage of this method is that it does not require segmentation/localization of other retinal structures which are required as prior knowledge in existing fovea localization methods. The key idea is to enhance the relative contrast between the fovea and its surrounding such that it is well-separated in a retinal image. The approach was tested on 520 retinal images which include normal and pathological color retinal images. An overall accuracy of 90.57% is reported independent of left eye, right eye, normal or pathological retinal images.