Progressive structure-based alignment of homologous proteins: Adopting sequence comparison strategies

Comparison of multiple protein structures has a broad range of applications in the analysis of protein structure, function and evolution. Multiple structure alignment tools (MSTAs) are necessary to obtain a simultaneous comparison of a family of related folds. In this study, we have developed a method for multiple structure comparison largely based on sequence alignment techniques. A widely used Structural Alphabet named Protein Blocks (PBs) was used to transform the information on 3D protein backbone conformation as a ID sequence string. A progressive alignment strategy similar to CLUSTALW was adopted for multiple PB sequence alignment (mulPBA). Highly similar stretches identified by the pairwise alignments are given higher weights during the alignment. The residue equivalences from PB based alignments are used to obtain a three dimensional fit of the structures followed by an iterative refinement of the structural superposition. Systematic comparisons using benchmark datasets of MSTAs underlines that the alignment quality is better than MULTIPROT, MUSTANG and the alignments in HOMSTRAD, in more than 85% of the cases. Comparison with other rigid-body and flexible MSTAs also indicate that mulPBA alignments are superior to most of the rigid-body MSTAs and highly comparable to the flexible alignment methods. (C) 2012 Elsevier Masson SAS. All rights reserved.

Targeting ribosome assembly on the HCV RNA using a small RNA molecule

Translation initiation of hepatitis C virus (HCV) RNA is the initial obligatory step of the viral life cycle, mediated through the internal ribosome entry site (IRES) present in the 5'-untranslated region (UTR). Initiation on the HCV IRES is mediated by multiple structure-specific interactions between IRES RNA and host 40S ribosomal subunit. In the present study we demonstrate that the SLIIIef domain, in isolation from other structural elements of HCV IRES, retain the ability to interact with 40S ribosome subunit. A small RNA SLRef, mimicking the SLIIIef domain was found to interact specifically with human La protein and the ribosomal protein S5 and selectively inhibit HCV RNA translation. More importantly, SLRef RNA showed significant suppression of replication in HCV monocistronic replicon and decrease of negative strand synthesis in HCV cell culture system. Finally, using Sendai virus based virosome, the targeted delivery of SLRef RNA into mice liver succeeded in selectively inhibiting HCV IRES mediated translation in vivo.

Sequence and structural basis for chromosomal fragility during translocations in cancer

Chromosomal aberration is considered to be one of the major characteristic features in many cancers. Chromosomal translocation, one type of genomic abnormality, can lead to deregulation of critical genes involved in regulating important physiological functions such as cell proliferation and DNA repair. Although chromosomal translocations were thought to be random events, recent findings suggest that certain regions in the human genome are more susceptible to breakage than others. The possibility of deviation from the usual B-DNA conformation in such fragile regions has been an active area of investigation. This review summarizes the factors that contribute towards the fragility of these regions in the chromosomes, such as DNA sequences and the role of different forms of DNA structures. Proteins responsible for chromosomal fragility, and their mechanism of action are also discussed. The effect of positioning of chromosomes within the nucleus favoring chromosomal translocations and the role of repair mechanisms are also addressed.

Draft Genome Sequence of Staphylococcus aureus ST672, an Emerging Disease Clone from India

We report the draft genome sequence of methicillin-resistant Staphylococcus aureus (MRSA) strain ST672, an emerging disease clone in India, from a septicemia patient. The genome size is about 2.82 Mb with 2,485 open reading frames (ORFs). The staphylococcal cassette chromosome mec (SCCmec) element (type V) and immune evasion cluster appear to be different from those of strain ST772 on preliminary examination.

Generation and analysis of drought stressed subtracted expressed sequence tags from safflower (Carthamus tinctorius L.)

Drought is the most crucial environmental factor that limits productivity of many crop plants. Exploring novel genes and gene combinations is of primary importance in plant drought tolerance research. Stress tolerant genotypes/species are known to express novel stress responsive genes with unique functional significance. Hence, identification and characterization of stress responsive genes from these tolerant species might be a reliable option to engineer the drought tolerance. Safflower has been found to be a relatively drought tolerant crop and thus, it has been the choice of study to characterize the genes expressed under drought stress. In the present study, we have evaluated differential drought tolerance of two cultivars of safflower namely, A1 and Nira using selective physiological marker traits and we have identified cultivar A1 as relatively drought tolerant. To identify the drought responsive genes, we have constructed a stress subtracted cDNA library from cultivar A1 following subtractive hybridization. Analysis of similar to 1,300 cDNA clones resulted in the identification of 667 unique drought responsive ESTs. Protein homology search revealed that 521 (78 %) out of 667 ESTs showed significant similarity to known sequences in the database and majority of them previously identified as drought stress-related genes and were found to be involved in a variety of cellular functions ranging from stress perception to cellular protection. Remaining 146 (22 %) ESTs were not homologous to known sequences in the database and therefore, they were considered to be unique and novel drought responsive genes of safflower. Since safflower is a stress-adapted oil-seed crop this observation has great relevance. In addition, to validate the differential expression of the identified genes, expression profiles of selected clones were analyzed using dot blot (reverse northern), and northern blot analysis. We showed that these clones were differentially expressed under different abiotic stress conditions. The implications of the analyzed genes in abiotic stress tolerance are discussed in our study.

The defect sequence for contractive tuples

We introduce the defect sequence for a contractive tuple of Hilbert space operators and investigate its properties. The defect sequence is a sequence of numbers, called defect dimensions associated with a contractive tuple. We show that there are upper bounds for the defect dimensions. The tuples for which these upper bounds are obtained, are called maximal contractive tuples. The upper bounds are different in the non-commutative and in the commutative case. We show that the creation operators on the full Fock space and the coordinate multipliers on the Drury-Arveson space are maximal. We also study pure tuples and see how the defect dimensions play a role in their irreducibility. (C) 2012 Elsevier Inc. All rights reserved.

Plasmodium falciparum Prp16 homologue and its role in splicing

Large numbers of Plasmodium genes have been predicted to have introns. However, little information exists on the splicing mechanisms in this organism. Here, we describe the DExD/DExH-box containing Pre-mRNA processing proteins (Prps), PfPrp2p, PfPrp5p, PfPrp16p, PfPrp22p, PfPrp28p, PfPrp43p and PfBrr2p, present in the Plasmodium falciparum genome and characterized the role of one of these factors, PfPrp16p. It is a member of DEAH-box protein family with nine collinear sequence motifs, a characteristic of helicase proteins. Experiments with the recombinantly expressed and purified PfPrp16 helicase domain revealed binding to RNA, hydrolysis of ATP as well as catalytic helicase activities. Expression of helicase domain with the C-terminal helicase-associated domain (HA2) reduced these activities considerably, indicating that the helicase-associated domain may regulate the PfPrp16 function. Localization studies with the PfPrp16 GFP transgenic lines suggested a role of its N-terminal domain (1-80 amino acids) in nuclear targeting. Immunodepletion of PfPrp16p, from nuclear extracts of parasite cultures, blocked the second catalytic step of an in vitro constituted splicing reaction suggesting a role for PfPrp16p in splicing catalysis. Further we show by complementation assay in yeast that a chimeric yeast-Plasmodium Prp16 protein, not the full length PfPrp16, can rescue the yeast prp16 temperature-sensitive mutant. These results suggest that although the role of Prp16p in catalytic step II is highly conserved among Plasmodium, human and yeast, subtle differences exist with regards to its associated factors or its assembly with spliceosomes. (C) 2012 Elsevier B.V. All rights reserved.

Translocation and encapsulation of siRNA inside carbon nanotubes

We report spontaneous translocation of small interfering RNA (siRNA) inside carbon nanotubes (CNTs) of various diameters and chirality using all atom molecular dynamics simulations with explicit solvent. We use umbrella sampling method to calculate the free energy landscape of the siRNA entry and translocation event. Free energy profiles show that siRNA gains free energy while translocating inside CNT, and barrier for siRNA exit from CNT ranges from 40 to 110 kcal/mol depending on CNT chirality and salt concentration. The translocation time tau decreases with the increase of CNT diameter with a critical diameter of 24 angstrom for the translocation. In contrast, double strand DNA of the same sequence does not translocate inside CNT due to large free energy barrier for the translocation. This study helps in understanding the nucleic acid transport through nanopores at microscopic level and may help designing carbon nanotube based sensor for siRNA. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4773302]

Improved Detection of Remote Homologues Using Cascade PSI-BLAST: Influence of Neighbouring Protein Families on Sequence Coverage

Background: Development of sensitive sequence search procedures for the detection of distant relationships between proteins at superfamily/fold level is still a big challenge. The intermediate sequence search approach is the most frequently employed manner of identifying remote homologues effectively. In this study, examination of serine proteases of prolyl oligopeptidase, rhomboid and subtilisin protein families were carried out using plant serine proteases as queries from two genomes including A. thaliana and O. sativa and 13 other families of unrelated folds to identify the distant homologues which could not be obtained using PSI-BLAST. Methodology/Principal Findings: We have proposed to start with multiple queries of classical serine protease members to identify remote homologues in families, using a rigorous approach like Cascade PSI-BLAST. We found that classical sequence based approaches, like PSI-BLAST, showed very low sequence coverage in identifying plant serine proteases. The algorithm was applied on enriched sequence database of homologous domains and we obtained overall average coverage of 88% at family, 77% at superfamily or fold level along with specificity of similar to 100% and Mathew's correlation coefficient of 0.91. Similar approach was also implemented on 13 other protein families representing every structural class in SCOP database. Further investigation with statistical tests, like jackknifing, helped us to better understand the influence of neighbouring protein families. Conclusions/Significance: Our study suggests that employment of multiple queries of a family for the Cascade PSI-BLAST searches is useful for predicting distant relationships effectively even at superfamily level. We have proposed a generalized strategy to cover all the distant members of a particular family using multiple query sequences. Our findings reveal that prior selection of sequences as query and the presence of neighbouring families can be important for covering the search space effectively in minimal computational time. This study also provides an understanding of the `bridging' role of related families.

A cyclic peptide mimic of an RNA recognition motif of human La protein is a potent inhibitor of hepatitis C virus

Due to limited available therapeutic options, developing new lead compounds against hepatitis C virus is an urgent need. Human La protein stimulates hepatitis C virus translation through interaction with the hepatitis C viral RNA. A cyclic peptide mimicking the beta-turn of the human La protein that interacts with the viral RNA was synthesized. It inhibits hepatitis C viral RNA translation significantly better than the corresponding linear peptide at longer post-treatment times. The cyclic peptide also inhibited replication as measured by replicon RNA levels using real time RT-PCR. The cyclic peptide emerges as a promising lead compound against hepatitis C.

Human la protein interaction with GCAC near the initiator AUG enhances hepatitis C virus RNA replication by promoting linkage between 5 ` and 3 ` untranslated regions

Human La protein has been implicated in facilitating the internal initiation of translation as well as replication of hepatitis C virus (HCV) RNA. Previously, we demonstrated that La interacts with the HCV internal ribosome entry site (IRES) around the GCAC motif near the initiator AUG within stem-loop IV by its RNA recognition motif (RRM) (residues 112 to 184) and influences HCV translation. In this study, we have deciphered the role of this interaction in HCV replication in a hepatocellular carcinoma cell culture system. We incorporated mutation of the GCAC motif in an HCV monocistronic subgenomic replicon and a pJFH1 construct which altered the binding of La and checked HCV RNA replication by reverse transcriptase PCR (RT-PCR). The mutation drastically affected HCV replication. Furthermore, to address whether the decrease in replication is a consequence of translation inhibition or not, we incorporated the same mutation into a bicistronic replicon and observed a substantial decrease in HCV RNA levels. Interestingly, La overexpression rescued this inhibition of replication. More importantly, we observed that the mutation reduced the association between La and NS5B. The effect of the GCAC mutation on the translation-to-replication switch, which is regulated by the interplay between NS3 and La, was further investigated. Additionally, our analyses of point mutations in the GCAC motif revealed distinct roles of each nucleotide in HCV replication and translation. Finally, we showed that a specific interaction of the GCAC motif with human La protein is crucial for linking 5' and 3' ends of the HCV genome. Taken together, our results demonstrate the mechanism of regulation of HCV replication by interaction of the cis-acting element GCAC within the HCV IRES with human La protein.

The sequence and structure of snake gourd (Trichosanthes anguina) seed lectin, a three-chain nontoxic homologue of type II RIPs

The sequence and structure of snake gourd seed lectin (SGSL), a nontoxic homologue of type II ribosome-inactivating proteins (RIPs), have been determined by mass spectrometry and X-ray crystallography, respectively. As in type II RIPs, the molecule consists of a lectin chain made up of two beta-trefoil domains. The catalytic chain, which is connected through a disulfide bridge to the lectin chain in type II RIPs, is cleaved into two in SGSL. However, the integrity of the three-dimensional structure of the catalytic component of the molecule is preserved. This is the first time that a three-chain RIP or RIP homologue has been observed. A thorough examination of the sequence and structure of the protein and of its interactions with the bound methyl-alpha-galactose indicate that the nontoxicity of SGSL results from a combination of changes in the catalytic and the carbohydrate-binding sites. Detailed analyses of the sequences of type II RIPs of known structure and their homologues with unknown structure provide valuable insights into the evolution of this class of proteins. They also indicate some variability in carbohydrate-binding sites, which appears to contribute to the different levels of toxicity exhibited by lectins from various sources.

Energy Hyperspace for Stacking Interaction in AU/AU Dinucleotide Step: Dispersion-Corrected Density Functional Theory Study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

Co-evolution of specific amino acid in sigma 1.2 region and nucleotide base in the discriminator to act as sensors of small molecule effectors of transcription initiation in mycobacteria

The transcription from rrn and a number of other promoters is regulated by initiating ribonucleotides (iNTPs) and guanosine tetra/penta phosphate (p)ppGpp], either by strengthening or by weakening of the RNA polymerase (RNAP)-promoter interactions during initiation. Studies in Escherichia coli revealed the importance of a sequence termed discriminator, located between -10 and the transcription start site of the responsive promoters in this mode of regulation. Instability of the open complex at these promoters is attributed to the lack of stabilizing interactions between the suboptimal discriminator and the 1.2 region of sigma 70 (Sig70) in RNAP holoenzyme. We demonstrate a different pattern of interaction between the promoters and sigma A (SigA) of Mycobacterium tuberculosis to execute similar regulation. Instead of cytosine and methionine, thymine at three nucleotides downstream to -10 element and leucine 232 in SigA are found to be essential for iNTPs and pppGpp mediated response at the rrn and gyr promoters of the organism. The specificity of the interaction is substantiated by mutational replacements, either in the discriminator or in SigA, which abolish the nucleotide mediated regulation in vitro or in vivo. Specific yet distinct bases and the amino acids appear to have co-evolved' to retain the discriminator-sigma 1.2 region regulatory switch operated by iNTPs/pppGpp during the transcription initiation in different bacteria.

Inactivation of the Bacterial RNA Polymerase Due to Acquisition of Secondary Structure by the omega Subunit

The widely conserved omega subunit encoded by rpoZ is the smallest subunit of Escherichia coli RNA polymerase (RNAP) but is dispensable for bacterial growth. Function of omega is known to be substituted by GroEL in omega-null strain, which thus does not exhibit a discernable phenotype. In this work, we report isolation of omega variants whose expression in vivo leads to a dominant lethal phenotype. Studies show that in contrast to omega, which is largely unstructured, omega mutants display substantial acquisition of secondary structure. By detailed study with one of the mutants, omega(6) bearing N60D substitution, the mechanism of lethality has been deciphered. Biochemical analysis reveals that omega(6) binds to beta ` subunit in vitro with greater affinity than that of omega. The reconstituted RNAP holoenzyme in the presence of omega(6) in vitro is defective in transcription initiation. Formation of a faulty RNAP in the presence of mutant omega results in death of the cell. Furthermore, lethality of omega(6) is relieved in cells expressing the rpoC2112 allele encoding beta ` (2112), a variant beta ` bearing Y457S substitution, immediately adjacent to the beta ` catalytic center. Our results suggest that the enhanced omega(6)-beta ` interaction may perturb the plasticity of the RNAP active center, implicating a role for omega and its flexible state.