Design of LNA Analogues Using a Combined Density Functional Theory and Molecular Dynamics Approach for RNA Therapeutics.

Antisense therapeutics treat a wide spectrum of diseases, many of which cannot be addressed with the current drug technologies. In the quest to design better antisense oligonucleotide drugs, we propose five novel LNA analogues (A1-A5) for modifying antisense oligonucleotides and establishing each with the five standard nucleic acids: adenine (A), guanine (G), cytosine (C), thymine (T), and uracil (U). Monomer nucleotides of these modifications were considered for a detailed Density Functional Theory (DFT)-based quantum chemical analysis to determine their molecular-level structural and electronic properties. A detailed MD simulation study was done on a 14-mer ASO (5'-CTTAGCACTGGCCT-3') containing these modifications targeting PTEN mRNA. Results from both molecular- and oligomer-level analysis clearly depicted LNA-level stability of the modifications, the ASO/RNA duplexes maintaining stable Watson-Crick base pairing preferring RNA-mimicking A-form duplexes. Notably, monomer MO isosurfaces for both purines and pyrimidines were majorly distributed on the nucleobase region in modifications A1 and A2 and in the bridging unit in modifications A3, A4, and A5, suggesting that A3/RNA, A4/RNA, and A5/RNA duplexes interact more with the RNase H and solvent environment. Accordingly, solvation of A3/RNA, A4/RNA, and A5/RNA duplexes was higher compared to that of LNA/RNA, A1/RNA, and A2/RNA duplexes. This study has resulted in a successful archetype for creating advantageous nucleic acid modifications tailored for particular needs, fulfilling a useful purpose of designing novel antisense modifications, which may overcome the drawbacks and improve the pharmacokinetics of existing LNA antisense modifications.

Structural insight into locked nucleic acid based novel antisense modifications: A DFT calculations at monomer and MD simulations at oligomer level.

In the present study, five novel LNA based antisense modifications have been proposed. A conformational search was carried out using TANGO, followed by geometry optimization using MOPAC. Based on their electronic energies the most stable conformation for each modification was identified. Further, DFT based full geometry optimization on the most stable conformations at the gas phase B3LYP/6-31G(d,p) using a Gaussian03 and single point energy calculations on the optimized structures at the solvent phase B3LYP/6-311G(d,p) level of theory were done to derive their quantum chemical descriptors using the Gaussian09. A comparison of global reactivity descriptors confirmed that the LNA based modifications were the most reactive. Base-pair stability was recorded by observing the binding energies and base-pairing conformations of modified GC base pairs at the B3LYP/6-311G(d,p) level of theory. Molecular dynamics simulations have been performed at the oligomer duplex level by incorporating individual modifications on 20-mer RNA-RNA duplexes using AMBER16. Free energy calculations of duplex structures suggested that incorporation of A2 modification into the RNA-RNA duplex increased the duplex binding affinity similar to LNA. Whereas, the A3 modification showed less binding compared to LNA but improved binding compared to MOE. This computational approach using quantum chemical methods may be very useful to propose better modifications than the existing ones before performing the experiments in the area of antisense technology.

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations.

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6-41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA. Communicated by Ramaswamy H. Sarma.

Acetylcholinesterase and Aß Aggregation Inhibition by Heterometallic Ruthenium(II)-Platinum(II) Polypyridyl Complexes.

Two heteronuclear ruthenium(II)-platinum(II) complexes [Ru(bpy)2(BPIMBp)PtCl2]2+ (3) and [Ru(phen)2(BPIMBp)PtCl2]2+ (4), where bpy = 2,2'-bipyridine, phen = 1,10-phenanthroline, and BPIMBp = 1,4'-bis[(2-pyridin-2-yl)-1H-imidazol-1-ylmethyl]-1,1'-biphenyl, have been designed and synthesized from their mononuclear precursors [Ru(bpy)2(BPIMBp)]2+ (1) and [Ru(phen)2(BPIMBp)]2+ (2) as multitarget molecules for Alzheimer's disease (AD). The inclusion of the cis-PtCl2 moiety facilitates the covalent interaction of Ru(II) polypyridyl complexes with amyloid ß (Aß) peptide. These multifunctional complexes act as inhibitors of acetylcholinesterase (AChE), Aß aggregation, and Cu-induced oxidative stress and protect neuronal cells against Aß-toxicity. The study highlights the design of metal based anti-Alzheimer's disease (AD) systems.

Chickpea-Fusarium oxysporum interaction transcriptome reveals differential modulation of plant defense strategies.

Fusarium wilt is one of the major biotic stresses reducing chickpea productivity. The use of wilt-resistant cultivars is the most appropriate means to combat the disease and secure productivity. As a step towards understanding the molecular basis of wilt resistance in chickpea, we investigated the transcriptomes of wilt-susceptible and wilt-resistant cultivars under both Fusarium oxysporum f.sp. ciceri (Foc) challenged and unchallenged conditions. Transcriptome profiling using LongSAGE provided a valuable insight into the molecular interactions between chickpea and Foc, which revealed several known as well as novel genes with differential or unique expression patterns in chickpea contributing to lignification, hormonal homeostasis, plant defense signaling, ROS homeostasis, R-gene mediated defense, etc. Similarly, several Foc genes characteristically required for survival and growth of the pathogen were expressed only in the susceptible cultivar with null expression of most of these genes in the resistant cultivar. This study provides a rich resource for functional characterization of the genes involved in resistance mechanism and their use in breeding for sustainable wilt-resistance. Additionally, it provides pathogen targets facilitating the development of novel control strategies.

Ruthenium(II) polypyridyl complexes with hydrophobic ancillary ligand as Aß aggregation inhibitors.

The synthesis, spectral and electrochemical characterization of the complexes of the type [Ru(NN)2(txbg)](2+) where NN is 2,2'-bipyridine (bpy) (1), 1,10-phenanthroline (phen) (2), dipyrido [3,2-d:2',3f] quinoxaline (dpq) (3), and dipyrido[3,2-a:2',3'-c]phenazine (dppz) (4) which incorporate the tetra-xylene bipyridine glycoluril (txbg) as the ancillary ligand are described in detail. Crystal structures of ligand txbg and complex 2 were solved by single crystal X-ray diffraction. Thioflavin T (ThT) fluorescence and Transmission Electron Microscopy (TEM) results indicated that at micromolar concentration all complexes exhibit significant potential of Aß aggregation inhibition, while the ligand txbg displayed weak activity towards Aß aggregation. Complex 1 showed relatively low inhibition (70%) while complexes 2-4 inhibited nearly 100% Aß aggregation after 240 h of incubation. The similar potential of complexes 2-4 and absence of any trend in their activity with the planarity of polypyridyl ligands suggests there is no marked effect of planarity of coligands on their inhibitory potential. Further studies on acetylcholinesterase (AChE) inhibition indicated very weak activity of these complexes against AChE. Detailed interactions of Aß with both ligand and complex 2 have been studied by molecular modeling. Complex 2 showed interactions involving all three polypyridyl ligands with hydrophobic region of Aß. Furthermore, the toxicity of these complexes towards human neuroblastoma cells was evaluated by MTT assay and except complex 4, the complexes displayed very low toxicity.

Auto-regulation of Slug mediates its activity during epithelial to mesenchymal transition.

Slug, a five C2H2 zinc finger (ZF) motif transcription factor mediates cell migration in development, adult tissue repair and regeneration, as well as during tumor metastases through epithelial to mesenchymal transition. At the molecular level, this involves interactions with E-box (CACC/GGTG) consensus elements within target gene promoters to achieve transcriptional repression. However, precise elucidation of events involved in this DNA recognition and binding of specific promoters to regulate target genes have not been achieved. In the present study, we show that besides transcriptional repression, Slug can also directly activate its own expression by preferential binding to specific E-box elements in the distal binding region of its promoter. Our findings suggest that while the first ZF does not contribute to the transcription-associated functions of Slug, all the remaining four ZFs are involved in regulating the expression of target genes with ZF3 and ZF4 being more crucial than ZF2 or ZF5. We also report that recognition and binding preferences of ZFs are defined through intrinsic differences in the E-box core base pairs and/or flanking sequences, with the S2 E-box element being most critical during autoregulation. However, specific target E-box recognition and binding are also defined by the cellular context, which implies that in silico and/or biochemical DNA binding preferences may not necessarily be able to accurately predict in situ events. Our studies thus constitute a novel understanding of transcriptional regulation.

Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study.

Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.

Ruthenium(II) polypyridyl complex as inhibitor of acetylcholinesterase and Aß aggregation.

Two ruthenium(II) polypyridyl complexes [Ru(phen)3](2+) (1) and [Ru(phen)2(bxbg)](2+) (2) (where phen = 1,10 phenanthroline, bxbg = bis(o-xylene)bipyridine glycoluril) have been evaluated for acetylcholinesterase (AChE) and Amyloid-ß peptide (Aß) aggregation inhibition. Complex 2 exhibits higher potency of AChE inhibition and kinetics and molecular modeling studies indicate that ancillary ligand plays significant role in inhibitory potency exhibited by complex 2. The inhibitory effect of these complexes on Aß (1-40) aggregation is investigated using Thioflavin T fluorescence and Transmission Electron Microscopy. Both complexes efficiently inhibit Aß (1-40) aggregation and are negligibly toxic to human neuroblastoma cells. This is the first demonstration that ruthenium(II) polypyridyl complexes simultaneously inhibit AChE and Aß aggregation.

MD simulations of HIV-1 RT primer-template complex: effect of modified nucleosides and antisense PNA oligomer.

Human immunodeficiency virus type 1 (HIV-1) requires the human tRNA(3)(Lys) as a reverse transcriptase (RT) primer. The annealing of 3' terminal 18 nucleotides of tRNA(3)(Lys) with the primer binding site (PBS) of viral RNA (vRNA) is crucial for reverse transcription. Additional contacts between the A rich (A-loop) region of vRNA and the anticodon domain of tRNA(3)(Lys) are necessary, which show the specific requirement of tRNA(3)(Lys). The importance of modified nucleosides, present in tRNA(3)(Lys), in giving stability to the primer-template complex has been determined in earlier experiments. It has been observed that the PNA oligomer targeted to PBS of vRNA destabilized the crucial interactions between primer and template due to which the reverse transcription is inhibited. Molecular dynamics simulations have been carried out to study the effect of modified nucleosides on the vRNA-tRNA(3)(Lys) complex stability and the destabilization effect of PNA oligomer on the vRNA-tRNA(3)(Lys)-PNA complex. The root-mean-square deviation, hydrogen bonding, tertiary interactions, and free energy calculations of the simulation data support the experimental results. The analyses have revealed the structural changes in PBS region of vRNA which might be another strong reason for the inability of RT binding to 7F helix for its normal functioning of reverse transcription.

Genome sequences of Salmonella enterica serovar typhimurium, Choleraesuis, Dublin, and Gallinarum strains of well- defined virulence in food-producing animals.

Salmonella enterica is an animal and zoonotic pathogen of worldwide importance and may be classified into serovars differing in virulence and host range. We sequenced and annotated the genomes of serovar Typhimurium, Choleraesuis, Dublin, and Gallinarum strains of defined virulence in each of three food-producing animal hosts. This provides valuable measures of intraserovar diversity and opportunities to formally link genotypes to phenotypes in target animals.

Mixed ligand cobalt(II) picolinate complexes: synthesis, characterization, DNA binding and photocleavage.

Complexes of the type [Co(pic)(2)(NN)], where pic = picolinate, NN = dipyrido[3,2-d:2',3'-f]quinoxaline (dpq) (4) and 4b,5,7,7a-tetrahydro-4b,7a-epiminomethanoimino-6H-imidazo[4,5-f][1,10]-phenanthroline-6,13-dione (bipyridyl-glycoluril) (bpg) (6) have been synthesized and characterized by elemental analysis, IR, UV-vis, NMR and ESI-MS spectroscopy and thermogravimetic analysis (TGA). Their physicochemical properties are compared with previously synthesized complexes, where NN = (H(2)O)(2) (1), 2,2'-bipyridine (bpy) (2), 1,10-phenanthroline (phen) (3) and dipyrido[3,2-a:2',3'-c]phenazine (dppz) (5). The crystal structures of the complexes 4-6 were solved by single-crystal X-ray diffraction. The complexes 4 and 5 crystallize from a mixture of chloroform and methanol in monoclinic and orthorhombic crystal systems, respectively, whereas complex 6 crystallizes from dimethyl sulfoxide (DMSO) in a tetragonal crystal system. The coordination sphere consists of two oxygen atoms and two nitrogen atoms from the two picolinates and two nitrogen atoms from the dpq, dppz or bpg ligand, respectively. Co(ii)/Co(iii) oxidation potentials have been determined by cyclic voltammetry. The DNA binding of complexes 1-5 has been investigated using thermal melting, fluorescence quenching and viscosity measurements, which indicate the partial intercalation of complex 5 with an apparent binding constant (k(app)) of 8.3 × 10(5) M(-1). DNA cleavage studies of complexes 1-5 have been investigated using gel electrophoresis in the presence of H(2)O(2) as an oxidizing agent and also by photoirradiation at 365 nm. The mechanistic investigations suggest that singlet oxygen ((1)O(2)) is the major species involved in the DNA cleavage by these complexes. The structures of complexes 2-6 were optimized with density functional theory (DFT) method (B3LYP/6-31G(d,p)). The low vertical ionization potential values indicate photoredox pathways for the DNA cleavage activity by complexes 4 and 5, which is corroborated by DNA cleavage experiments.

Molecular dynamics simulations of cyclohexyl modified peptide nucleic acids (PNA).

Peptide Nucleic Acids (PNA) that bind sequence specifically to DNA/RNA are of major interest in the field of molecular biology and could form the basis for gene-targeted drugs. Molecular dynamics simulations are aimed to characterize the structural and dynamical features to understand the effect of backbone modification on the structure and dynamics along with the stability of the resulting 10mer complexes of PNA with DNA/RNA. Twelve Molecular Dynamics (MD) simulations of duplexes and triplexes with and without cyclohexyl modification were carried out for 10ns each. The simulations indicate that the cyclohexyl modification with different stereoisomers has influenced all the PNA-DNA/RNA complexes. Modification has added rigidity to backbone by restricting beta to +60 in case of (1R,2S) cyclohexyl PNA and to -60 in case of (1S,2R) cyclohexyl PNA. The results of MD simulations were able to show the backbone rigidification and preference for RNA complexes over DNA due to presence of cyclohexyl ring in the PNA backbone.

Generation, annotation, and analysis of ESTs from midgut tissue of adult female Anopheles stephensi mosquitoes.

BACKGROUND: Malaria is a tropical disease caused by protozoan parasite, Plasmodium, which is transmitted to humans by various species of female anopheline mosquitoes. Anopheles stephensi is one such major malaria vector in urban parts of the Indian subcontinent. Unlike Anopheles gambiae, an African malaria vector, transcriptome of A. stephensi midgut tissue is less explored. We have therefore carried out generation, annotation, and analysis of expressed sequence tags from sugar-fed and Plasmodium yoelii infected blood-fed (post 24 h) adult female A. stephensi midgut tissue. RESULTS: We obtained 7061 and 8306 ESTs from the sugar-fed and P. yoelii infected mosquito midgut tissue libraries, respectively. ESTs from the combined dataset formed 1319 contigs and 2627 singlets, totaling to 3946 unique transcripts. Putative functions were assigned to 1615 (40.9%) transcripts using BLASTX against UniProtKB database. Amongst unannotated transcripts, we identified 1513 putative novel transcripts and 818 potential untranslated regions (UTRs). Statistical comparison of annotated and unannotated ESTs from the two libraries identified 119 differentially regulated genes. Out of 3946 unique transcripts, only 1387 transcripts were mapped on the A. gambiae genome. These also included 189 novel transcripts, which were mapped to the unannotated regions of the genome. The EST data is available as ESTDB at http://mycompdb.bioinfo-portal.cdac.in/cgi-bin/est/index.cgi. CONCLUSION: 3946 unique transcripts were successfully identified from the adult female A. stephensi midgut tissue. These data can be used for microarray development for better understanding of vector-parasite relationship and to study differences or similarities with other malaria vectors. Mapping of putative novel transcripts from A. stephensi on the A. gambiae genome proved fruitful in identification and annotation of several genes. Failure of some novel transcripts to map on the A. gambiae genome indicates existence of substantial genomic dissimilarities between these two potent malaria vectors.

Ruthenium(II) complexes of bipyridine-glycoluril and their interactions with DNA.

Nine complexes of the type [Ru(N-N)(2)(BPG)]Cl(2) 1-4, [Ru(N-N)(BPG)(2)]Cl(2) 5-8, and [Ru(BPG)(3)]Cl(2) 9 where N-N is 2,2'-bipyridine (bpy), 1,10-phenanthroline (phen), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq), dipyrido[3,2-a:2',3'-c]phenazine (dppz), which incorporates bipyridine-glycoluril (BPG-4b,5,7,7a-tetrahydro-4b,7a-epiminomethanoimino-6H-imidazo[4,5-f][1,10]phenanthroline-6,13-dione) as the ancillary ligand, have been synthesized and characterized. These complexes with the peripheral polypyridyl ligands have the ability to form conjugates with DNA. The DNA binding (absorption spectroscopy, steady-state and time-resolved emission measurements, steady-state emission quenching measurements) and cleavage (under dark and irradiated conditions) by these complexes has been studied to investigate the influence of the ancillary ligand. The binding ability of these complexes to DNA is dependent on the planarity of the intercalative polypyridyl ligand, which is further affected by the ancillary bipyridine-glycoluril ligand. The complexes 3, 4, 7, and 8 bind to CT-DNA with binding constants on the order of 10(4) M(-1). Time-resolved emission measurements on the DNA-bound complexes 1, 3, 5-7, and 9 show monoexponential decay of the excited states, whereas complexes 2, 4, and 8 show biexponential decay with short- and long-lived components. Interaction of complexes 2-9 with plasmid pBR322 DNA studied by gel electrophoresis experiments reveals that all complexes cleave DNA efficiently at micromolar concentrations under dark and anaerobic conditions probably by a hydrolytic mechanism. Complexes 3, 4, 7, 8, and [Ru(bpy)(2)(dppz)](2+) show extensive DNA cleavage in the presence of light with a shift in mobility of form I of DNA probably due to the high molecular weight of DNA-complex conjugates. However, the extent of the cleavage is augmented on irradiation in the case of complexes 3, 4, 7, and 8, which include the planar dpq and dppz ligands, suggesting a combination of hydrolytic and oxidative mechanism for the DNA scission. Molecular mechanics calculations of these systems corroborate the DNA binding and cleavage mechanisms.

Study of early events in the protein folding of villin headpiece using molecular dynamics simulation.

Protein folding is scientifically and computationally challenging problem. The early phases of protein folding are interesting due to various events like nascent secondary structure formation, hydrophobic collapse leading to formation of non-native or meta-stable conformations. These events occur within a very short time span of 100 ns as compared to total folding time of few microseconds. It is highly difficult to observe these events experimentally due to very short lifetime. Molecular dynamics simulation technique can efficiently probe the detailed atomic level understanding about these events. In the present paper, all atom molecular dynamics simulation trajectory of nearly 200 ns was carried out for fully solvated villin headpiece with PME treatment using AMBER 7 package. Initial hydrophobic collapse along with secondary structure formation resulted into formation of partially stable non-native conformations. The formation of secondary structural elements and hydrophobic collapse takes place simultaneously in the folding process.

Synthesis, characterization, X-ray structure and DNA photocleavage by cis-dichloro bis(diimine) Co(III) complexes.

Complexes of the type [Co(LL)2Cl2]Cl, where LL = N,N'-ethylenediamine (en), 2,2'-bipyridine (bpy), 1,10-phenanthroline (phen), 1,10-phenanthroline-5,6-dione (phendione) and dipyrido[3,2-a:2',3'-c]phenazine (dppz) have been synthesized and characterized by elemental analyses, IR, UV-visible and NMR spectroscopy. Crystal structure of [Co(phendione)2Cl2]Cl x 0.5 HCl x 3.5 H2O has been solved and refined to R = 0.0552. The crystal is monoclinic with space group C2/c; a = 25.730(2) A, b = 12.375(1) A, c = 18.979(2) A, beta = 119.925(1) degrees and Z = 8. The DNA binding characteristics of the complexes, investigated by covalent binding assay, viscosity measurements and competitive binding fluorescence measurements show that the complexes interact with DNA covalently except the complex containing the planar dppz ligand which intercalates within the base pairs of DNA. The complexes containing en, phen and phendione cleave plasmid pBR 322 DNA upon irradiation under aerobic conditions while the complex containing the dppz ligand cleaves DNA upon irradiation under inert atmosphere. Molecular modeling studies show that the minimized structure of [Co(phendione)2Cl2]+, maintained the octahedral structure while binding to the N7 of guanines and the ligand fits into the major groove without disrupting the helical structure of the B-DNA.

Motif detection in Arabidopsis: correlation with gene expression data.

The genes having similar expression profiles are considered to have common regulatory mechanisms and are controlled by the binding of transcription factors to the regulatory elements present in their upstream regions. The detection of cis-regulatory elements can help in further understanding of co-expression of genes. This paper deals with the detection of motifs in the upstream regions of genes involved in diurnal rhythms of Arabidopsis and also deals with the correlation of expression data with sequence information. We detected motifs in the upstream regions of genes involved in diurnal cycles and checked for their presence in circadian regulated, dark induced and in light induced genes of Arabidopsis. Ten motifs were reported in this study, out of which five were already reported in available transcription factor databases as the elements involved in light responsiveness. Significance study of ten motifs was done by taking random sets of same data size. One of the ten motifs namely GGCCCA, which was found without any base variations in 62 genes, was further studied by analyzing the expression profiles of its respective genes within the set of diurnal regulated genes using SOM clustering method. It was found that the genes were clustered together into two major groups, out of which one group had glycine rich proteins and the second group had genes belonging to dehydrogenase and oxidoreductase family.

Accelerating comparative genomics using parallel computing.

In the past decade there has been an increase in the number of completely sequenced genomes due to the race of multibillion-dollar genome-sequencing projects. The enormous biological sequence data thus flooding into the sequence databases necessitates the development of efficient tools for comparative genome sequence analysis. The information deduced by such analysis has various applications viz. structural and functional annotation of novel genes and proteins, finding gene order in the genome, gene fusion studies, constructing metabolic pathways etc. Such study also proves invaluable for pharmaceutical industries, such as in silico drug target identification and new drug discovery. There are various sequence analysis tools available for mining such useful information of which FASTA and Smith-Waterman algorithms are widely used. However, analyzing large datasets of genome sequences using the above codes seems to be impractical on uniprocessor machines. Hence there is a need for improving the performance of the above popular sequence analysis tools on parallel cluster computers. Performance of the Smith-Waterman (SSEARCH) and FASTA programs were studied on PARAM 10000, a parallel cluster of workstations designed and developed in-house. FASTA and SSEARCH programs, which are available from the University of Virginia, were ported on PARAM and were optimized. In this era of high performance computing, where the paradigm is shifting from conventional supercomputers to the cost-effective general-purpose cluster of workstations and PCs, this study finds extreme relevance. Good performance of sequence analysis tools on a cluster of workstations was demonstrated, which is important for accelerating identification of novel genes and drug targets by screening large databases.