Sub-genomic RNA of defective interfering (D.I.) dengue viral particles is replicated in the same manner as full length genomes

The predicted secondary structure of sub-genomic RNA in dengue virus defective interfering (D.I.) particles from patients, or generated in vitro, resembled that of the 3′ and 5′ regions of wild type dengue virus (DENV) genomes. While these structures in the sub-genomic RNA were found to be essential for its replication, their nucleotide sequences were not, so long as any new sequences maintained wild type RNA secondary structure. These observations suggested that these sub-genomic fragments of RNA from dengue viruses were replicated in the same manner as the full length genomes of their wild type, “helper”, viruses and that they probably represent the smallest fragments of DENV RNA that can be replicated during a natural infection. While D.I. particles containing sub-genomic RNA are completely parasitic, the relationship between wild type and D.I. DENV may be symbiotic, with the D.I. particles enhancing the transmission of infectious DENV.

RNASeqBrowser: A genome browser for simultaneous visualization of raw strand specific RNAseq reads and UCSC genome browser custom tracks

Background Strand specific RNAseq data is now more common in RNAseq projects. Visualizing RNAseq data has become an important matter in Analysis of sequencing data. The most widely used visualization tool is the UCSC genome browser that introduced the custom track concept that enabled researchers to simultaneously visualize gene expression at a particular locus from multiple experiments. Our objective of the software tool is to provide friendly interface for visualization of RNAseq datasets. Results This paper introduces a visualization tool (RNASeqBrowser) that incorporates and extends the functionality of the UCSC genome browser. For example, RNASeqBrowser simultaneously displays read coverage, SNPs, InDels and raw read tracks with other BED and wiggle tracks -- all being dynamically built from the BAM file. Paired reads are also connected in the browser to enable easier identification of novel exon/intron borders and chimaeric transcripts. Strand specific RNAseq data is also supported by RNASeqBrowser that displays reads above (positive strand transcript) or below (negative strand transcripts) a central line. Finally, RNASeqBrowser was designed for ease of use for users with few bioinformatic skills, and incorporates the features of many genome browsers into one platform. Conclusions The features of RNASeqBrowser: (1) RNASeqBrowser integrates UCSC genome browser and NGS visualization tools such as IGV. It extends the functionality of the UCSC genome browser by adding several new types of tracks to show NGS data such as individual raw reads, SNPs and InDels. (2) RNASeqBrowser can dynamically generate RNA secondary structure. It is useful for identifying non-coding RNA such as miRNA. (3) Overlaying NGS wiggle data is helpful in displaying differential expression and is simple to implement in RNASeqBrowser. (4) NGS data accumulates a lot of raw reads. Thus, RNASeqBrowser collapses exact duplicate reads to reduce visualization space. Normal PC’s can show many windows of NGS individual raw reads without much delay. (5) Multiple popup windows of individual raw reads provide users with more viewing space. This avoids existing approaches (such as IGV) which squeeze all raw reads into one window. This will be helpful for visualizing multiple datasets simultaneously. RNASeqBrowser and its manual are freely available at http://www.australianprostatecentre.org/research/software/rnaseqbrowser webcite or http://sourceforge.net/projects/rnaseqbrowser/ webcite

Crystal structures of the endoplasmic reticulum aminopeptidase-1 (ERAP1) reveal the molecular basis for N-terminal peptide trimming

Endoplasmatic reticulum aminopeptidase 1 (ERAP1) is a multifunctional enzyme involved in trimming of peptides to an optimal length for presentation by major histocompatibility complex (MHC) class I molecules. Polymorphisms in ERAP1 have been associated with chronic inflammatory diseases, including ankylosing spondylitis (AS) and psoriasis, and subsequent in vitro enzyme studies suggest distinct catalytic properties of ERAP1 variants. To understand structure-activity relationships of this enzyme we determined crystal structures in open and closed states of human ERAP1, which provide the first snapshots along a catalytic path. ERAP1 is a zinc-metallopeptidase with typical H-E-X-X-H-(X)18-E zinc binding and G-A-M-E-N motifs characteristic for members of the gluzincin protease family. The structures reveal extensive domain movements, including an active site closure as well as three different open conformations, thus providing insights into the catalytic cycle. A K 528R mutant strongly associated with AS in GWAS studies shows significantly altered peptide processing characteristics, which are possibly related to impaired interdomain interactions.

Watson-Crick pairing, the Heisenberg group and Milnor invariants

We study the secondary structure of RNA determined by Watson-Crick pairing without pseudo-knots using Milnor invariants of links. We focus on the first non-trivial invariant, which we call the Heisenber invariant. The Heisenberg invariant, which is an integer, can be interpreted in terms of the Heisenberg group as well as in terms of lattice paths. We show that the Heisenberg invariant gives a lower bound on the number of unpaired bases in an RNA secondary structure. We also show that the Heisenberg invariant can predict allosteric structures for RNA. Namely, if the Heisenberg invariant is large, then there are widely separated local maxima (i.e., allosteric structures) for the number of Watson-Crick pairs found.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Characterization of the restriction enzyme-like endonuclease encoded by the Entamoeba histolytica non-long terminal repeat retrotransposon EhLINE1

The genome of the human pathogen Entamoeba histolytica, a primitive protist, contains non-long terminal repeat retrotransposable elements called EhLINEs. These encode reverse transcriptase and endonuclease required for retrotransposition. The endonuclease shows sequence similarity with bacterial restriction endonucleases. Here we report the salient enzymatic features of one such endonuclease. The kinetics of an EhLINE1-encoded endonuclease catalyzed reaction, determined under steady-state and single-turnover conditions, revealed a significant burst phase followed by a slower steady-state phase, indicating that release of product could be the slower step in this reaction. For circular supercoiled DNA the K-m was 2.6 x 10-8 m and the k(cat) was 1.6 x 10-2 sec-1. For linear E. histolytica DNA substrate the K-m and k(cat) values were 1.3 x 10-8 m and 2.2 x 10-4 sec-1 respectively. Single-turnover reaction kinetics suggested a noncooperative mode of hydrolysis. The enzyme behaved as a monomer. While Mg2+ was required for activity, 60% activity was seen with Mn2+ and none with other divalent metal ions. Substitution of PDX12-14D (a metal-binding motif) with PAX(12-14)D caused local conformational change in the protein tertiary structure, which could contribute to reduced enzyme activity in the mutated protein. The protein underwent conformational change upon the addition of DNA, which is consistent with the known behavior of restriction endonucleases. The similarities with bacterial restriction endonucleases suggest that the EhLINE1-encoded endonuclease was possibly acquired from bacteria through horizontal gene transfer. The loss of strict sequence specificity for nicking may have been subsequently selected to facilitate spread of the retrotransposon to intergenic regions of the E. histolytica genome.

X-ray crystallographic and NMR studies of pantothenate synthetase provide insights into the mechanism of homotropic inhibition by pantoate

The structural basis for the homotropic inhibition of pantothenate synthetase by the substrate pantoate was investigated by X-ray crystallography and high-resolution NMR spectroscopic methods. The tertiary structure of the dimeric N-terminal domain of Escherichia coli pantothenate synthetase, determined by X-ray crystallography to a resolution of 1.7 Å, showed a second molecule of pantoate bound in the ATP-binding pocket. Pantoate binding to the ATP-binding site induced large changes in structure, mainly for backbone and side chain atoms of residues in the ATP binding HXGH(34–37) motif. Sequence-specific NMR resonance assignments and solution secondary structure of the dimeric N-terminal domain, obtained using samples enriched in 2H, 13C, and 15N, indicated that the secondary structural elements were conserved in solution. Nitrogen-15 edited two-dimensional solution NMR chemical shift mapping experiments revealed that pantoate, at 10 mm, bound at these two independent sites. The solution NMR studies unambiguously demonstrated that ATP stoichiometrically displaced pantoate from the ATP-binding site. All NMR and X-ray studies were conducted at substrate concentrations used for enzymatic characterization of pantothenate synthetase from different sources [Jonczyk R & Genschel U (2006) J Biol Chem 281, 37435–37446]. As pantoate binding to its canonical site is structurally conserved, these results demonstrate that the observed homotropic effects of pantoate on pantothenate biosynthesis are caused by competitive binding of this substrate to the ATP-binding site. The results presented here have implications for the design and development of potential antibacterial and herbicidal agents.

Influence of phosphate ligands in abolishing the conformational difference between ribonuclease A and its acid-denatured derivative

The initial structural alteration of RNAase A due to acid denaturation (0.5 N HCl, 30 degrees C) that accompanies deamidation (without altering enzymic activity) has been dectected by spectrophotometric titration, fluorescence and ORD/CD measurements. It is shown that acid treated RNAase A has an altered conformation at neutral pH, 25 degrees C. This is characterized by the increased accessibility of buried tyrosine residue(s) towards the solvent. The most altered conformation of RNAase A is found in the 10 h acid-treated derivative. This has about 1.5 additional exposed tyrosine residues and a lesser amount of secondary structure than RNAase A. All three methods (titration, fluorescence and CD) established that the structural transition of RNAase A is biphasic. The first phase occurs within 1 h and the resulting subtle conformational change is constant up to 7 h. Following this, after the release of 0.55 mol of ammonia, the major conformational change begins. The altered conformation of the acid-denatured RNAase A could be reversed completely to the native state through a conformational change induced by substrate analogs like 2'- or 3'-CMP. Thus the monodeamidated derivative isolated from the acid-denatured RNAase A by phosphate is very similar to RNAase A in over-all conformation. The results suggest the possibility of flexibility in the RNAase A molecule that does not affect its catalytic activity, as probed through the tyrosine residues.

A novel approach to design of cis-acting DNA structural element for regulation of gene expression in vivo

Taking advantage of the degeneracy of the genetic code we have developed a novel approach to introduce, within a gene, DNA sequences capable of adopting unusual structures and to investigate the role of such sequences in regulation of gene expression in vivo. We used a computer program that generates alternative codon sequences for the same amino-acid sequence to convert a stretch of nucleotides into an inverted-repeat sequence with the potential to adopt cruciform structure. This approach was used to replace a 51-base-pair EcoRI-HindIII segment in the N-terminal region of the beta-galactosidase gene in plasmid pUC19 with a 51-bp synthetic oligonucleotide sequence with the potential to adopt a cruciform structure with 18 bp in the stem region. In selecting the 51-bp sequence, care was taken to include those codons that are preferred in E. coli. E. coli DH5-alpha cells harbouring the plasmid containing the redesigned sequence showed drastic reduction in expression of the beta-galactosidase gene compared to cells harbouring the plasmid with the native sequence. This approach demonstrates the possibility of introducing DNA secondary-structure elements to alter regulation of gene expression in vivo.

Biochemical and structural characterization of recombinant hyoscyamine 6 beta-hydroxylase from Datura metel L.

Hyoscyamine 6 beta-hydroxylase (H6H; EC 1.14.11.11), an important enzyme in the biosynthesis of tropane alkaloids, catalyzes the hydroxylation of hyoscyamine to give 6 beta-hydroxyhyoscyamine and its epoxidation in the biosynthetic pathway leading to scopolamine. Datura metel produces scopolamine as the predominant tropane alkaloid. The cDNA encoding H6H from D. mete! (DmH6H) was cloned, heterologously expressed and biochemically characterized. The purified recombinant His-tagged H6H from D. mete! (DmrH6H) was capable of converting hyoscyamine to scopolamine. The functionally expressed DmrH6H was confirmed by HPLC and ESI-MS verification of the products, 6 beta-hydroxyhyoscyamine and its derivative, scopolamine; the DmrH6H epoxidase activity was low compared to the hydroxylase activity. The K-m values for both the substrates, hyoscyamine and 2-oxoglutarate, were 50 mu M each. The CD (circular dichroism) spectrum of the DmrH6H indicated a preponderance of alpha-helicity in the secondary structure. From the fluorescence studies, Stern-Volmer constants for hyoscyamine and 2-oxoglutarate were found to be 0.14 M-1 and 0.56 M-1, respectively. These data suggested that the binding of the substrates, hyoscyamine and 2-oxoglutarate, to the enzyme induced significant conformational changes. (C) 2010 Elsevier Masson SAS. All rights reserved.

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

Beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.

Effect of constraints by threonine on proline containing αa-helix—A molecular dynamics approach

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12−Thr-Pro-(Ala)10−NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr- (Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

5-Fluorouracil induces structural changes in rat liver tRNA

5-fluorouracil (FUra) has been shown to modulate the aminoacylation function of rat liver tRNA. The present study was aimed at studying the structure-function relationship of FUra-substituted tRNA. Male Wistar rats (2-3 month old) were given a single i.p. injection of FUra at 50, 250, or 500 mg/kg body wt. and FUra-substituted total liver tRNA, i.e. tRNA(FUra50, 250, and 500, respectively, were isolated 3 h later. Normal tRNA (tRNA(N)) was isolated from saline-treated control rats. Thermal denaturation studies showed higher melting temperatures for tRNA(FUra) compared to tRNA(N). Heat denaturation followed by renaturation of total tRNA did not affect the activity of tRNA(N) and tRNA(FUra50), where as tRNA(FUra250 and 500) lost 35% and 72% of activity, respectively, compared to the corresponding group of non-denatured tRNA. Antibodies specific to rat liver tRNA recognized normal and FUra-substituted tRNA in the order of tRNA(N) > tRNA(FUra50) > or = tRNA(FUra250) > tRNA(FUra500) in an avidin-biotin micro-enzyme linked immunosorbant assay. tRNA(N) or tRNA(FUra50) preincubated with tRNA antiserum showed 74% and 59% of aminoacylation activity, respectively, compared to that of corresponding tRNA preincubated with normal rabbit IgG. However, activities of similarly treated tRNA(FUra250 and 500) were not affected. The observations of possible changes in the secondary structure of rat liver tRNA upon incorporation of FUra are discussed.

Studies on molecular evolution and structural features of double-headed inhibitors of ?-amylase and trypsin in plants

Plant seeds usually have high concentrations of proteinase and amylase inhibitors. These inhibitors exhibit a wide range of specificity, stability and oligomeric structure. In this communication, we report analysis of sequences that show statistically significant similarity to the double-headed alpha-amylase/trypsin inhibitor of ragi (Eleusine coracana). Our aim is to understand their evolutionary and structural features. The 14 sequences of this family that are available in the SWISSPROT database form three evolutionarily distinct branches. The branches relate to enzyme specificities and also probably to the oligomeric state of the proteins and not to the botanical class of the plant from which the enzymes are derived. This suggests that the enzyme specificities of the inhibitors evolved before the divergence of commercially cultivated cereals. The inhibitor sequences have three regions that display periodicity in hydrophobicity. It is likely that this feature reflects extended secondary structure in these segments. One of the most variable regions of the polypeptide corresponds to a loop, which is most probably exposed in the native structure of the inhibitors and is responsible for the inhibitory property.