Classification of Protein Kinases on the Basis of Both Kinase and Non-Kinase Regions

Background: Protein phosphorylation is a generic way to regulate signal transduction pathways in all kingdoms of life. In many organisms, it is achieved by the large family of Ser/Thr/Tyr protein kinases which are traditionally classified into groups and subfamilies on the basis of the amino acid sequence of their catalytic domains. Many protein kinases are multidomain in nature but the diversity of the accessory domains and their organization are usually not taken into account while classifying kinases into groups or subfamilies. Methodology: Here, we present an approach which considers amino acid sequences of complete gene products, in order to suggest refinements in sets of pre-classified sequences. The strategy is based on alignment-free similarity scores and iterative Area Under the Curve (AUC) computation. Similarity scores are computed by detecting common patterns between two sequences and scoring them using a substitution matrix, with a consistent normalization scheme. This allows us to handle full-length sequences, and implicitly takes into account domain diversity and domain shuffling. We quantitatively validate our approach on a subset of 212 human protein kinases. We then employ it on the complete repertoire of human protein kinases and suggest few qualitative refinements in the subfamily assignment stored in the KinG database, which is based on catalytic domains only. Based on our new measure, we delineate 37 cases of potential hybrid kinases: sequences for which classical classification based entirely on catalytic domains is inconsistent with the full-length similarity scores computed here, which implicitly consider multi-domain nature and regions outside the catalytic kinase domain. We also provide some examples of hybrid kinases of the protozoan parasite Entamoeba histolytica. Conclusions: The implicit consideration of multi-domain architectures is a valuable inclusion to complement other classification schemes. The proposed algorithm may also be employed to classify other families of enzymes with multidomain architecture.

High-quality annotation of promoter regions for 913 bacterial genomes

Motivation: The number of bacterial genomes being sequenced is increasing very rapidly and hence, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The present work addresses this requirement and presents a generic method applicable across organisms. Results: Relative stability of the DNA double helical sequences has been used to discriminate promoter regions from non-promoter regions. Based on the difference in stability between neighboring regions, an algorithm has been implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content. Threshold values to identify promoter regions have been derived using sequences flanking a subset of translation start sites from all microbial genomes and then used to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset.

Intraspecific variation in brain size and architecture : population divergence and phenotypic plasticity

Brain size and architecture exhibit great evolutionary and ontogenetic variation. Yet, studies on population variation (within a single species) in brain size and architecture, or in brain plasticity induced by ecologically relevant biotic factors have been largely overlooked. Here, I address the following questions: (i) do locally adapted populations differ in brain size and architecture, (ii) can the biotic environment induce brain plasticity, and (iii) do locally adapted populations differ in levels of brain plasticity? In the first two chapters I report large variation in both absolute and relative brain size, as well as in the relative sizes of brain parts, among divergent nine-spined stickleback (Pungitius pungitius) populations. Some traits show habitat-dependent divergence, implying natural selection being responsible for the observed patterns. Namely, marine sticklebacks have relatively larger bulbi olfactorii (chemosensory centre) and telencephala (involved in learning) than pond sticklebacks. Further, I demonstrate the importance of common garden studies in drawing firm evolutionary conclusions. In the following three chapters I show how the social environment and perceived predation risk shapes brain development. In common frog (Rana temporaria) tadpoles, I demonstrate that under the highest per capita predation risk, tadpoles develop smaller brains than in less risky situations, while high tadpole density results in enlarged tectum opticum (visual brain centre). Visual contact with conspecifics induces enlarged tecta optica in nine-spined sticklebacks, whereas when only olfactory cues from conspecifics are available, bulbus olfactorius become enlarged.Perceived predation risk results in smaller hypothalami (complex function) in sticklebacks. Further, group-living has a negative effect on relative brain size in the competition-adapted pond sticklebacks, but not in the predation-adapted marine sticklebacks. Perceived predation risk induces enlargement of bulbus olfactorius in pond sticklebacks, but not in marine sticklebacks who have larger bulbi olfactorii than pond fish regardless of predation. In sum, my studies demonstrate how applying a microevolutionary approach can help us to understand the enormous variation observed in the brains of wild animals a point-of-view which I high-light in the closing review chapter of my thesis.

Can equipartition fields produce the tilts of bipolar magnetic regions?

The effect of turbulence on the nonaxisymmetric flux rings of equipartition field strength in bipolar magnetic regions is studied on the basis of the small-scale momentum exchange mechanism and the giant cell drag combined with the Kelvin-Helmholtz drag mechanism. It is shown that the giant cell drag and small-scale momentum exchange mechanism can make equipartition flux loops emerge at low latitudes, in addition to making them exhibit the observed tilts. However, the sizes of the flux tubes have to be restricted to a couple of hundred kilometers. An ad hoc constraint on the footpoints of the flux loops is introduced by not letting them move in the phi direction, and it is found that equipartition fields of any size can be made to emerge at sunspot latitudes with the observed tilts by suitably adjusting the footpoint separations.

Nucleotide sequence of the first cosmid from the Mycobacterium leprae genome project: structure and function of the Rif-Str regions

The nucleotide sequence of cosmid B1790, carrying the Rif-Str regions of the Mycobacterium leprae chromosome, has been determined. Twelve open reading frames were identified in the 36716bp sequence, representing 40% of the coding capacity. Five ribosomal proteins, two elongation factors and the β and β'subunits of RNA polymerase have been characterized and two novel genes were found. One of these encodes a member of the so-called ABC family of ATP-binding proteins while the other appears to encode an enzyme involved in repairing genomic lesions caused by free radicals. This finding may well be significant as M. leprae, an intracellular pathogen, lives within macrophages.

An ensemble of B-DNA dinucleotide geometries lead to characteristic nucleosomal DNA structure and provide plasticity required for gene expression

Background: A nucleosome is the fundamental repeating unit of the eukaryotic chromosome. It has been shown that the positioning of a majority of nucleosomes is primarily controlled by factors other than the intrinsic preference of the DNA sequence. One of the key questions in this context is the role, if any, that can be played by the variability of nucleosomal DNA structure. Results: In this study, we have addressed this question by analysing the variability at the dinucleotide and trinucleotide as well as longer length scales in a dataset of nucleosome X-ray crystal structures. We observe that the nucleosome structure displays remarkable local level structural versatility within the B-DNA family. The nucleosomal DNA also incorporates a large number of kinks. Conclusions: Based on our results, we propose that the local and global level versatility of B-DNA structure may be a significant factor modulating the formation of nucleosomes in the vicinity of high-plasticity genes, and in varying the probability of binding by regulatory proteins. Hence, these factors should be incorporated in the prediction algorithms and there may not be a unique `template' for predicting putative nucleosome sequences. In addition, the multimodal distribution of dinucleotide parameters for some steps and the presence of a large number of kinks in the nucleosomal DNA structure indicate that the linear elastic model, used by several algorithms to predict the energetic cost of nucleosome formation, may lead to incorrect results.

Inefficient excision of uracil from loop regions of DNA oligomers by E.coli uracil DNA glycosylase

Kinetic parameters for uracil DNA glycosylase (E.coli)-catalysed excision of uracil from DNA oligomers containing dUMP in different structural contexts were determined. Our results show that single-stranded oligonucleotides (unstructured) are used as somewhat better substrates than the double-stranded oligonucleotides. This is mainly because of the favourable V-max value of the enzyme for single-stranded substrates. More interestingly, however, we found that uracil release from loop regions of DNA hairpins is extremely inefficient. The poor efficiency with which uracil is excised from loop regions is a result of both increased K-m and lowered V-max values. This observation may have significant implications in uracil DNA glycosylase-directed repair of DNA segments that can be extruded as hairpins. In addition, these studies are useful in designing oligonucleotides for various applications in DNA research where the use of uracil DNA glycosylase is sought.

Star formation in giant extragalactic H II regions

Star formation properties in Giant Extragalactic H II Regions (GEHRs) are investigated using optical photometry and evolutionary population synthesis models. Photometric data in $BVR$ bands and in the emission line of H-alpha are obtained by CCD imaging at Vainu Bappu Observatory, Kavalur. Aperture photometry is performed for 180 GEHRs in galaxies NGC 1365, 1566, 2366, 2903, 2997, 3351, 4303, 4449, 4656 and 5253. Thirty six of these GEHRs having published spectroscopic data are studied for star formation properties. The population synthesis model is constructed based on Maeder's stellar evolutionary and Kurucz stellar atmosphere models, to synthesize observational quantities of embedded clusters in GEHRs. The observed H-alpha luminosity is a measure of the number of massive stars while the contribution to BVR bands is from intermediate mass (5-15 solar mass) stars when the cluster is young and from evolving supergiants when the cluster is old (age >/= 6~Myr). Differential reddening between gas and embedded stars is essential to constrain the dereddened cluster colors within the range of youngest clusters. Obscuring dust closely associated with gas, which is distributed in filaments and clumps, as in the case of 30 Doradus, is the most likely configuration giving rise to net reduction of extinction towards stars. The fraction of the stellar photons escaping the nebula unattenuated is estimated to be 50%. GEHRs are rarely found to be simple systems containing stars from single generation. In the present sample such regions in addition to being older than 3~Myr, have their Lyman continuum luminosity reduced by as much as 60%, compared to the observed $B$ band luminosity for a normal IMF. The missing ionizing photons may be escaping the nebula, leading to the ionization of extra-H II region ionized medium. Co-existence of young (age

Interactions of 3' terminal and 5' terminal regions of physalis mottle virus genomic RNA with its replication complex

Physalis mottle virus (PhMV) belongs to the tymogroup of positive-strand RNA viruses with a genome size of 6 kb. Crude membrane preparations from PhMV-infected Nicotiana glutinosa plants catalyzed the synthesis of PhMV genomic RNA from endogenously bound template. Addition of exogenous genomic RNA enhanced the synthesis which was specifically inhibited by the addition of sense and antisense transcripts corresponding to 3' terminal 242 nucleotides as well as the 5' terminal 458 nucleotides of PhMV genomic RNA while yeast tRNA or ribosomal RNA failed to inhibit the synthesis. This specific inhibition suggested that the 5' and 3' non-coding regions of PhMV RNA might play an important role in viral replication.

Interfacial Regions Governing Internal Cavities of Dendrimers. Studies of Poly(alkyl aryl ether) Dendrimers Constituted with Linkers of Varying Alkyl Chain Length

This report deals with a study of the properties of internal cavities of dendritic macromolecules that are capable Of encapsulating and mediating photoreactions of guest molecules. The internal cavity structures of dendrimers are determined by the interfacial regions between the aqueous exterior and hydrocarbon like interior constituted by the linkers that connect symmetrically sited branch points constituting the dendrimer and head groups that cap the dendrimers. Phloroglucinol-based poly(alkyl aryl ether) dendrimers constituted with a homologous series of alkyl linkers were undertaken for the current study. Twelve dendrimers within first, second, and third generations, having ethyl, n-propyl, n-butyl, and n-pentyl groups as the linkers and hydroxyl groups at peripheries in each generation, were synthesized. Encapsulation of pyrene and coumarins by aqueous basic solutions of dendrimers were monitored-by UV-vis and fluorescence spectroscopies, which showed that a lower generation dendrimer with an optimal alkyl linker presented better encapsulation abilities than a higher generation dendrimer. Norrish type I photoreaction of dibenzyl ketone was carried out within the above: series of dendrimers to probe their abilities to hold guests and reactive inthermediate radical pairs within themselves. The extent of cage effect from the series of third generation dendrimers was observed to be higher with dendrimers having an n-pentyl group as the linker.

Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

On the large-scale diffuse magnetic field of the sun - II. The contribution of active regions

Dikpati and Choudhuri (1993, 1995) developed a model for the poleward migration of the weak diffuse magnetic field on the Sun's surface. This field was identified with the poloidal component produced by the solar dynamo operating at the base of the convection zone, and its evolution was studied by considering the effects of meridional circulation and turbulent diffusion. The earlier model is extended in this paper by incorporating the flux from, the decay of tilted active regions near the solar surface as an additional source of the poloidal field. This extended model can now explain various low-latitude features in the time-latitude diagram of the weak diffuse fields. These low-latitude features could not be accounted for in the earlier model, which was very successful in modeling the behavior at high latitudes. The time-latitude diagrams show that regions of a particular polarity often have 'tongues' of opposite polarity. Such tongues can be produced in the theoretical model by incorporating fluctuations in the source term arising out of the decaying active regions.

The orientational relaxation of bipolar active regions

In the mean, bipolar active regions are oriented nearly toroidally, according to Hale's polarity law, with a latitude-dependent tilt known as Joy's Law. The tilt angles of individual active regions deviate from this mean behavior and change over time. It has been found that on average the change is toward the mean angle at a rate characteristic of 4.37 days (Howard, 1996). We show that this orientational relaxation is consistent with the standard model of flux tube emergence from a deep dynamo layer. Under this scenario Joy's law results from the Coriolis effect on the rising flux tube (D'Silva and Choudhuri, 1993), and departures from it result from turbulent buffeting of the tubes (Longcope and Fisher, 1996). We show that relaxation toward Joy's angle occurs because the turbulent perturbations relax on shorter time scales than the perturbations from the Coriolis force. The turbulent perturbations relax more rapidly because they are localized to the topmost portion of the convection zone while the Coriolis perturbations are more widely distributed. If a fully-developed active region remains connected to the strong toroidal magnetic field at the base of the convection zone, its tilt will eventually disappear, leaving it aligned perfectly toroidally. On the other hand, if the flux becomes disconnected from the toroidal field the bipole will assume a tilt indicative of the location of disconnection. We compare models which are connected and disconnected from the toroidal field. Only those disconnected at points very deep in the convection zone a-re consistent with observed time scale of orientational relaxation.

Anthropogenic aerosol fraction over Afro-Asian regions inferred using Kalpana-I and MISR data

Mineral dust constitutes the single largest contributor of natural aerosols over continents. The first step towards separating natural aerosol radiative impact from its anthropogenic counterparts over continents is to gather information on dust aerosols. The infrared (IR) radiance (10.5-12.5 mu m) acquired from the Kalpana-I satellite (similar to 8-km resolution) was used to retrieve regional characteristics of dust aerosols over the Afro-Asian region during the winter of 2004, coinciding with a national aerosol campaign. Here, we used aerosol-induced IR radiance depression as an index of dust load. The regional distribution of dust over various arid and semi-arid regions of India and adjacent continents has been estimated, and these data in conjunction with regional maps of column aerosol optical depth (AOD) are used to infer anthropogenic aerosol fraction. Surprisingly, even over desert locations in India and Saudi Arabia, the anthropogenic fraction was relatively high (similar to 0.3 to 0.4) and the regionally averaged anthropogenic fraction over India was 0.62 +/- 0.06.

Analysis of shear bands in slow granular flows using a frictional Cosserat model

The slow flow of granular materials is often marked by the existence of narrow shear layers, adjacent to large regions that suffer little or no deformation. This behaviour, in the regime where shear stress is generated primarily by the frictional interactions between grains, has so far eluded theoretical description. In this paper, we present a rigid-plastic frictional Cosserat model that captures thin shear layers by incorporating a microscopic length scale. We treat the granular medium as a Cosserat continuum, which allows the existence of localised couple stresses and, therefore, the possibility of an asymmetric stress tensor. In addition, the local rotation is an independent field variable and is not necessarily equal to the vorticity. The angular momentum balance, which is implicitly satisfied for a classical continuum, must now be solved in conjunction with the linear momentum balances. We extend the critical state model, used in soil plasticity, for a Cosserat continuum and obtain predictions for flow in plane and cylindrical Couette devices. The velocity profile predicted by our model is in qualitative agreement with available experimental data. In addition, our model can predict scaling laws for the shear layer thickness as a function of the Couette gap, which must be verified in future experiments. Most significantly, our model can determine the velocity field in viscometric flows, which classical plasticity-based model cannot.