Bioinformatics tools and databases for the study of human growth hormone

Advances in novel molecular biological diagnostic methods are changing the way of diagnosis and study of metabolic disorders like growth hormone deficiency. Faster sequencing and genotyping methods require strong bioinformatics tools to make sense of the vast amount of data generated by modern laboratories. Advances in genome sequencing and computational power to analyze the whole genome sequences will guide the diagnostics of future. In this chapter, an overview of some basic bioinformatics resources that are needed to study metabolic disorders are reviewed and some examples of bioinformatics analysis of human growth hormone gene, protein and structure are provided.

Stereochemistry of Amino Acid Spacers Determines the Pharmacokinetics of (111)In-DOTA-Minigastrin Analogues for Targeting the CCK2/Gastrin Receptor.

The metabolic instability and high kidney retention of minigastrin (MG) analogues hamper their suitability for use in peptide-receptor radionuclide therapy of CCK2/gastrin receptor-expressing tumors. High kidney retention has been related to N-terminal glutamic acids and can be substantially reduced by coinjection of polyglutamic acids or gelofusine. The aim of the present study was to investigate the influence of the stereochemistry of the N-terminal amino acid spacer on the enzymatic stability and pharmacokinetics of (111)In-DOTA-(d-Glu)6-Ala-Tyr-Gly-Trp-Met-Asp-Phe-NH2 ((111)In-PP11-D) and (111)In-DOTA-(l-Glu)6-Ala-Tyr-Gly-Trp-Met-Asp-Phe-NH2 ((111)In-PP11-L). Using circular dichroism measurements, we demonstrate the important role of secondary structure on the pharmacokinetics of the two MG analogues. The higher in vitro serum stability together with the improved tumor-to-kidney ratio of the (d-Glu)6 congener indicates that this MG analogue might be a good candidate for further clinical study.

Double-target antisense U7 snRNAs promote efficient skipping of an aberrant exon in three human beta-thalassemic mutations.

We have used three beta-thalassemic mutations, IVS2-654, -705 and -745, that create aberrant 5' splice sites (5' ss) and activate a common cryptic 3' ss further upstream in intron 2 of the human beta-globin gene to optimize a generally applicable exon-skipping strategy using antisense derivatives of U7 small nuclear RNA (snRNA). Introducing a modified U7 snRNA gene carrying an antisense sequence against the cryptic 3' ss into cultured cells expressing the mutant beta-globin genes, restored correct beta-globin mRNA splicing for all three mutations, but the efficiency was much weaker for IVS2-654 than for the other mutations. The length of antisense sequence influenced the efficiency with an optimum of approximately 24 nucleotides. Combining two antisense sequences directed against different target sites in intron 2, either on separate antisense RNAs or, even better, on a single U7 snRNA, significantly enhanced the efficiency of splicing correction. One double-target U7 RNA was expressed on stable transformation resulting in permanent and efficient suppression of the IVS2-654 mutation and production of beta-globin. These results suggest that forcing the aberrant exon into a looped secondary structure may strongly promote its exclusion from the mRNA and that this approach may be used generally to induce exon skipping.

NMR determination of the major solution conformation of a peptoid pentamer with chiral side chains

Polymers of N-substituted glycines (“peptoids”) containing chiral centers at the α position of their side chains can form stable structures in solution. We studied a prototypical peptoid, consisting of five para-substituted (S)-N-(1-phenylethyl)glycine residues, by NMR spectroscopy. Multiple configurational isomers were observed, but because of extensive signal overlap, only the major isomer containing all cis-amide bonds was examined in detail. The NMR data for this molecule, in conjunction with previous CD spectroscopic results, indicate that the major species in methanol is a right-handed helix with cis-amide bonds. The periodicity of the helix is three residues per turn, with a pitch of ≈6 Å. This conformation is similar to that anticipated by computational studies of a chiral peptoid octamer. The helical repeat orients the amide bond chromophores in a manner consistent with the intensity of the CD signal exhibited by this molecule. Many other chiral polypeptoids have similar CD spectra, suggesting that a whole family of peptoids containing chiral side chains is capable of adopting this secondary structure motif. Taken together, our experimental and theoretical studies of the structural properties of chiral peptoids lay the groundwork for the rational design of more complex polypeptoid molecules, with a variety of applications, ranging from nanostructures to nonviral gene delivery systems.

Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15–25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

Chrysanthemum chlorotic mottle viroid: Unusual structural properties of a subgroup of self-cleaving viroids with hammerhead ribozymes

The causal agent of chrysanthemum chlorotic mottle (CChM) disease has been identified, cloned, and sequenced. It is a viroid RNA (CChMVd) of 398–399 nucleotides. In vitro transcripts with the complete CChMVd sequence were infectious and induced the typical symptoms of the CChM disease. CChMVd can form hammerhead structures in both polarity strands. Plus and minus monomeric CChMVd RNAs self-cleaved during in vitro transcription and after purification as predicted by these structures, which are stable and most probably act as single hammerhead structures as in peach latent mosaic viroid (PLMVd), but not in avocado sunblotch viroid (ASBVd). Moreover, the plus CChMVd hammerhead structure also appears to be active in vivo, because the 5′ terminus of the linear plus CChMVd RNA isolated from infected tissue is that predicted by the corresponding hammerhead ribozyme. Both hammerhead structures of CChMVd display some peculiarities: the plus self-cleaving domain has an unpaired A after the conserved A9 residue, and the minus one has an unusually long helix II. The most stable secondary structure predicted for CChMVd is a branched conformation that does not fulfill the rod-like or quasi-rod-like model proposed for the in vitro structure of most viroids with the exception of PLMVd, whose proposed secondary structure of lowest free energy also is branched. The unusual conformation of CChMVd and PLMVd is supported by their insolubility in 2 M LiCl, in contrast to ASBVd and a series of representative non-self-cleaving viroids that are soluble under the same high salt conditions. These results support the classification of self-cleaving viroids into two subgroups, one formed by ASBVd and the other one by PLMVd and CChMVd.

Three novel families of miniature inverted-repeat transposable elements are associated with genes of the yellow fever mosquito, Aedes aegypti

Three novel families of transposable elements, Wukong, Wujin, and Wuneng, are described in the yellow fever mosquito, Aedes aegypti. Their copy numbers range from 2,100 to 3,000 per haploid genome. There are high degrees of sequence similarity within each family, and many structural but not sequence similarities between families. The common structural characteristics include small size, no coding potential, terminal inverted repeats, potential to form a stable secondary structure, A+T richness, and putative 2- to 4-bp A+T-biased specific target sites. Evidence of previous mobility is presented for the Wukong elements. Elements of these three families are associated with 7 of 16 fully or partially sequenced Ae. aegypti genes. Characteristics of these mosquito elements indicate strong similarities to the miniature inverted-repeat transposable elements (MITEs) recently found to be associated with plant genes. MITE-like elements have also been reported in two species of Xenopus and in Homo sapiens. This characterization of multiple families of highly repetitive MITE-like elements in an invertebrate extends the range of these elements in eukaryotic genomes. A hypothesis is presented relating genome size and organization to the presence of highly reiterated MITE families. The association of MITE-like elements with Ae. aegypti genes shows the same bias toward noncoding regions as in plants. This association has potentially important implications for the evolution of gene regulation.

The increased level of β1,4-galactosyltransferase required for lactose biosynthesis is achieved in part by translational control

β1,4-Galactosyltransferase (β4GalT-I) participates in both glycoconjugate biosynthesis (ubiquitous activity) and lactose biosynthesis (mammary gland-specific activity). In somatic tissues, transcription of the mammalian β4GalT-I gene results in a 4.1-kb mRNA and a 3.9-kb mRNA as a consequence of initiation at two start sites separated by ≈200 bp. In the mammary gland, coincident with the increased β4GalT-I enzyme level (≈50-fold) required for lactose biosynthesis, there is a switch from the 4.1-kb start site to the preferential use of the 3.9-kb start site, which is governed by a stronger tissue-restricted promoter. The use of the 3.9-kb start site results in a β4GalT-I transcript in which the 5′- untranslated region (UTR) has been truncated from ≈175 nt to ≈28 nt. The 5′-UTR of the 4.1-kb transcript [UTR(4.1)] is predicted to contain extensive secondary structure, a feature previously shown to reduce translational efficiency of an mRNA. In contrast, the 5′-UTR of the 3.9-kb mRNA [UTR(3.9)] lacks extensive secondary structure; thus, this transcript is predicted to be more efficiently translated relative to the 4.1-kb mRNA. To test this prediction, constructs were assembled in which the respective 5′-UTRs were fused to the luciferase-coding sequence and enzyme levels were determined after translation in vitro and in vivo. The luciferase mRNA containing the truncated UTR(3.9) was translated more efficiently both in vitro (≈14-fold) and in vivo (3- to 5-fold) relative to the luciferase mRNA containing the UTR(4.1). Consequently, in addition to control at the transcriptional level, β4GalT-I enzyme levels are further augmented in the lactating mammary gland as a result of translational control.

Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Skipping of internal exons during removal of introns from pre-mRNA must be avoided for proper expression of most eukaryotic genes. Despite significant understanding of the mechanics of intron removal, mechanisms that ensure inclusion of internal exons in multi-intron pre-mRNAs remain mysterious. Using a natural two-intron yeast gene, we have identified distinct RNA–RNA complementarities within each intron that prevent exon skipping and ensure inclusion of internal exons. We show that these complementarities are positioned to act as intron identity elements, bringing together only the appropriate 5′ splice sites and branchpoints. Destroying either intron self-complementarity allows exon skipping to occur, and restoring the complementarity using compensatory mutations rescues exon inclusion, indicating that the elements act through formation of RNA secondary structure. Introducing new pairing potential between regions near the 5′ splice site of intron 1 and the branchpoint of intron 2 dramatically enhances exon skipping. Similar elements identified in single intron yeast genes contribute to splicing efficiency. Our results illustrate how intron secondary structure serves to coordinate splice site pairing and enforce exon inclusion. We suggest that similar elements in vertebrate genes could assist in the splicing of very large introns and in the evolution of alternative splicing.

Pentameric assembly of a neuronal glutamate transporter

Freeze-fracture electron microscopy was used to study the structure of a human neuronal glutamate transporter (EAAT3). EAAT3 was expressed in Xenopus laevis oocytes, and its function was correlated with the total number of transporters in the plasma membrane of the same cells. Function was assayed as the maximum charge moved in response to a series of transmembrane voltage pulses. The number of transporters in the plasma membrane was determined from the density of a distinct 10-nm freeze-fracture particle, which appeared in the protoplasmic face only after EAAT3 expression. The linear correlation between EAAT3 maximum carrier-mediated charge and the total number of the 10-nm particles suggested that this particle represented functional EAAT3 in the plasma membrane. The cross-sectional area of EAAT3 in the plasma membrane (48 ± 5 nm2) predicted 35 ± 3 transmembrane α-helices in the transporter complex. This information along with secondary structure models (6–10 transmembrane α-helices) suggested an oligomeric state for EAAT3. EAAT3 particles were pentagonal in shape in which five domains could be identified. They exhibited fivefold symmetry because they appeared as equilateral pentagons and the angle at the vertices was 110°. Each domain appeared to contribute to an extracellular mass that projects ≈3 nm into the extracellular space. Projections from all five domains taper toward an axis passing through the center of the pentagon, giving the transporter complex the appearance of a penton-based pyramid. The pentameric structure of EAAT3 offers new insights into its function as both a glutamate transporter and a glutamate-gated chloride channel.

Escherichia coli RNA polymerase terminates transcription efficiently at rho-independent terminators on single-stranded DNA templates

Several models have been proposed for the mechanism of transcript termination by Escherichia coli RNA polymerase at rho-independent terminators. Yager and von Hippel (Yager, T. D. & von Hippel, P. H. (1991) Biochemistry 30, 1097–118) postulated that the transcription complex is stabilized by enzyme–nucleic acid interactions and the favorable free energy of a 12-bp RNA–DNA hybrid but is destabilized by the free energy required to maintain an extended transcription bubble. Termination, by their model, is viewed simply as displacement of the RNA transcript from the hybrid helix by reformation of the DNA helix. We have proposed an alternative model where the RNA transcript is stably bound to RNA polymerase primarily through interactions with two single-strand specific RNA-binding sites; termination is triggered by formation of an RNA hairpin that reduces binding of the RNA to one RNA-binding site and, ultimately, leads to its ejection from the complex. To distinguish between these models, we have tested whether E. coli RNA polymerase can terminate transcription at rho-independent terminators on single-stranded DNA. RNA polymerase cannot form a transcription bubble on these templates; thus, the Yager–von Hippel model predicts that intrinsic termination will not occur. We find that transcript elongation on single-stranded DNA templates is hindered somewhat by DNA secondary structure. However, E. coli RNA polymerase efficiently terminates and releases transcripts at several rho-independent terminators on such templates at the same positions as termination occurs on duplex DNAs. Therefore, neither the nontranscribed DNA strand nor the transcription bubble is essential for rho-independent termination by E. coli RNA polymerase.

The European Large Subunit Ribosomal RNA Database

The European Large Subunit Ribosomal RNA Database compiles all complete or nearly complete large subunit ribosomal RNA sequences available from public sequence databases. These are provided in aligned format and the secondary structure, as derived by comparative sequence analysis, is included. Additional information about the sequences such as literature references and taxonomic information is also included. The database is available from our WWW server at http://rrna.uia.ac.be/lsu/.

RNA editing in Trypanosoma brucei: characterization of gRNA U-tail interactions with partially edited mRNA substrates

Guide RNAs (gRNAs), key components of the RNA editing reaction in Trypanosoma brucei, direct the insertion and deletion of uridylate (U) residues. Analyses of gRNAs reveal three functional elements. The 5′-end of the gRNA contains the anchor, which is responsible for selection and binding to the pre-edited mRNA. The second element (the guiding region) provides the information required for editing. At the 3′-end of the gRNA is a non-encoded U-tail, whose function remains unclear. However, the cleavage–ligation model for editing proposes that the U-tail binds to purine-rich regions upstream of editing sites, thereby strengthening the interaction and holding onto the 5′ cleavage product. Our previous studies demonstrated that the U-tail interacts with upstream sequences and may play roles in both stabilization and tethering. These studies also indicated that the U-tail interactions involved mRNA regions that were to be subsequently edited. This raised the question of what happens to the mRNA–U-tail interaction as editing proceeds in the 3′→5′ direction. We examined gCYb-558 and its U-tail interaction with 5′CYbUT and two partially edited 5′CYb substrates. Our results indicate that the 3′-end of the U-tail interacts with the same sequence in all three mRNAs. Predicted secondary structures using crosslinking data suggest that a similar structure is maintained as editing proceeds. These results indicate that the role of the U-tail may also involve maintenance of important secondary structure motifs.

Formation of a GNRA tetraloop in P5abc can disrupt an interdomain interaction in the Tetrahymena group I ribozyme

The secondary structure of a truncated P5abc subdomain (tP5abc, a 56-nucleotide RNA) of the Tetrahymena thermophila group I intron ribozyme changes when its tertiary structure forms. We have now used heteronuclear NMR spectroscopy to determine its conformation in solution. The tP5abc RNA that contains only secondary structure is extended compared with the tertiary folded form; both forms coexist in slow chemical exchange (the interconversion rate constant is slower than 1 s−1) in the presence of magnesium. Kinetic experiments have shown that tertiary folding of the P5abc subdomain is one of the earliest folding transitions in the group I intron ribozyme, and that it leads to a metastable misfolded intermediate. Previous mutagenesis studies suggest that formation of the extended P5abc structure described here destabilize a misfolded intermediate. This study shows that the P5abc RNA subdomain containing a GNRA tetraloop in P5c (in contrast to the five-nucleotide loop P5c in the tertiary folded ribozyme) can disrupt the base-paired interdomain (P14) interaction between P5c and P2.

Polynitrosylated proteins: characterization, bioactivity, and functional consequences.

Chemical modification of proteins is a common theme in their regulation. Nitrosylation of protein sulfhydryl groups has been shown to confer nitric oxide (NO)-like biological activities and to regulate protein functions. Several other nucleophilic side chains -- including those with hydroxyls, amines, and aromatic carbons -- are also potentially susceptible to nitrosative attack. Therefore, we examined the reactivity and functional consequences of nitros(yl)ation at a variety of nucleophilic centers in biological molecules. Chemical analysis and spectroscopic studies show that nitrosation reactions are sustained at sulfur, oxygen, nitrogen, and aromatic carbon centers, with thiols being the most reactive functionality. The exemplary protein, BSA, in the presence of a 1-, 20-, 100-, or 200-fold excess of nitrosating equivalents removes 0.6 +/- 0.2, 3.2 +/- 0.4, 18 +/- 4, and 38 +/- 10, respectively, moles of NO equivalents per mole of BSA from the reaction medium; spectroscopic evidence shows the proportionate formation of a polynitrosylated protein. Analogous reaction of tissue-type plasminogen activator yields comparable NO protein stoichiometries. Disruption of protein tertiary structure by reduction results in the preferential nitrosylation of up to 20 thus-exposed thiol groups. The polynitrosylated proteins exhibit antiplatelet and vasodilator activity that increases with the degree of nitrosation, but S-nitroso derivatives show the greatest NO-related bioactivity. Studies on enzymatic activity of tissue-type plasminogen activator show that polynitrosylation may lead to attenuated function. Moreover, the reactivity of tyrosine residues in proteins raises the possibility that NO could disrupt processes regulated by phosphorylation. Polynitrosylated proteins were found in reaction mixtures containing interferon-gamma/lipopolysaccharide-stimulated macrophages and in tracheal secretions of subjects treated with NO gas, thus suggesting their physiological relevance. In conclusion, multiple sites on proteins are susceptible to attack by nitrogen oxides. Thiol groups are preferentially modified, supporting the notion that S-nitrosylation can serve to regulate protein function. Nitrosation reactions sustained at additional nucleophilic centers may have (patho)physiological significance and suggest a facile route by which abundant NO bioactivity can be delivered to a biological system, with specificity dictated by protein substrate.