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.

Formation of a G-quadruplex at the BCL2 major breakpoint region of the t(14;18) translocation in follicular lymphoma

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Most breaks on chromosome 18 are located at the 3'-UTR of the BCL2 gene and are mainly clustered in the major breakpoint region (MBR). Recently, we found that the BCL2 MBR has a non-B DNA character in genomic DNA. Here, we show that single-stranded DNA modeled from the template strand of the BCL2 MBR, forms secondary structures that migrate faster on native PAGE in the presence of potassium, due to the formation of intramolecular G-quadruplexes. Circular dichroism shows evidence for a parallel orientation for G-quadruplex structures in the template strand of the BCL2 MBR. Mutagenesis and the DMS modification assay confirm the presence of three guanine tetrads in the structure. 1H nuclear magnetic resonance studies further confirm the formation of an intramolecular G-quadruplex and a representative model has been built based on all of the experimental evidence. We also provide data consistent with the possible formation of a G-quadruplex structure at the BCL2 MBR within mammalian cells. In summary, these important features could contribute to the single-stranded character at the BCL2 MBR, thereby contributing to chromosomal fragility.

Formation of a G-quadruplex at the BCL2 major breakpoint region of the t(14;18) translocation in follicular lymphoma

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Most breaks on chromosome 18 are located at the 3'-UTR of the BCL2 gene and are mainly clustered in the major breakpoint region (MBR). Recently, we found that the BCL2 MBR has a non-B DNA character in genomic DNA. Here, we show that single-stranded DNA modeled from the template strand of the BCL2 MBR, forms secondary structures that migrate faster on native PAGE in the presence of potassium, due to the formation of intramolecular G-quadruplexes. Circular dichroism shows evidence for a parallel orientation for G-quadruplex structures in the template strand of the BCL2 MBR. Mutagenesis and the DMS modification assay confirm the presence of three guanine tetrads in the structure. 1H nuclear magnetic resonance studies further confirm the formation of an intramolecular G-quadruplex and a representative model has been built based on all of the experimental evidence. We also provide data consistent with the possible formation of a G-quadruplex structure at the BCL2 MBR within mammalian cells. In summary, these important features could contribute to the single-stranded character at the BCL2 MBR, thereby contributing to chromosomal fragility.

Analysis of the effects of amino-acid-sequence on the structure of proteins

The conformation of amino acid side chains as observed in well-determined structures of globular proteins has earlier been extensively investigated. In contrast, the structural features of the polypeptide backbone that result from the occurrence of specific amino acids along the polypeptide have not been analysed. In this article, we present the statistically significant features in the backbone geometry that appear to be a consequence of the occurrence of rotamers of different amino acid side chains by analysing 102 well-refined structures that form a random collection of proteins. It is found that the persistence of helical segments around each residue is influenced by the residue type. Several residues exert asymmetrical influence between the carboxyl and amino terminal polypeptide segments. The degree to which secondary structures depart from an average geometry also appears to depend on residue type. These departures are correlated to the corresponding Chou and Fasman parameters of amino acid residues. The frequency distribution of the side chain rotamers is influenced by polypeptide secondary structure. In turn, the rotamer conformation of side chain affects the extension of the secondary structure of the backbone. The strongest correlation is found between the occurrence of g+ conformation and helix propagation on the carboxyl side of many residues.

Investigating Diproline Segments in Proteins: Occurrences, Conformation and Classification

The covalent linkage between the side-chain and the backbone nitrogen atom of proline leads to the formation of the five-membered pyrrolidine ring and hence restriction of the backbone torsional angle phi to values of -60 degrees +/- 30 degrees for the L-proline. Diproline segments constitute a chain fragment with considerably reduced conformational choices. In the current study, the conformational states for the diproline segment ((L)Pro-(L)Pro) found in proteins has been investigated with an emphasis on the cis and trans states for the Pro-Pro peptide bond. The occurrence of diproline segments in turns and other secondary structures has been studied and compared to that of Xaa-Pro-Yaa segments in proteins which gives us a better understanding on the restriction imposed on other residues by the diproline segment and the single proline residue. The study indicates that P(II)-P(II) and P(II)-alpha are the most favorable conformational states for the diproline segment. The analysis on Xaa-Pro-Yaa sequences reveals that the XaaPro peptide bond exists preferably as the trans conformer rather than the cis conformer. The present study may lead to a better understanding of the behavior of proline occurring in diproline segments which can facilitate various designed diproline-based synthetic templates for biological and structural studies. (C) 2011 Wiley Periodicals, Inc. Biopolymers 97: 54-64, 2012.

NMRof peptides

Developments and applications of NMR spectroscopy especially with biomolecules has taken big strides over the decades. This review gives a brief overview of peptide analysis by NMR as carried out in the author’s laboratory. A brief introduction to peptide biomolecules and NMR useful parameters are discussed in the beginning. This is followed by diagnostics features observed in NMR for identification of secondary structures. It further goes on to show how a three dimensional structure could be obtained by all-important NOE and hydrogen bond information. Use of heteronuclear experiments, which could be done at natural abundance is also highlighted in getting more details of peptide structures.Applications using Solid state NMR at natural abundance in connecting peptide solution and x-ray structures is demonstrated with couple of examples.

Local Structural Differences in Homologous Proteins: Specificities in Different SCOP Classes

The constant increase in the number of solved protein structures is of great help in understanding the basic principles behind protein folding and evolution. 3-D structural knowledge is valuable in designing and developing methods for comparison, modelling and prediction of protein structures. These approaches for structure analysis can be directly implicated in studying protein function and for drug design. The backbone of a protein structure favours certain local conformations which include alpha-helices, beta-strands and turns. Libraries of limited number of local conformations (Structural Alphabets) were developed in the past to obtain a useful categorization of backbone conformation. Protein Block (PB) is one such Structural Alphabet that gave a reasonable structure approximation of 0.42 angstrom. In this study, we use PB description of local structures to analyse conformations that are preferred sites for structural variations and insertions, among group of related folds. This knowledge can be utilized in improving tools for structure comparison that work by analysing local structure similarities. Conformational differences between homologous proteins are known to occur often in the regions comprising turns and loops. Interestingly, these differences are found to have specific preferences depending upon the structural classes of proteins. Such class-specific preferences are mainly seen in the all-beta class with changes involving short helical conformations and hairpin turns. A test carried out on a benchmark dataset also indicates that the use of knowledge on the class specific variations can improve the performance of a PB based structure comparison approach. The preference for the indel sites also seem to be confined to a few backbone conformations involving beta-turns and helix C-caps. These are mainly associated with short loops joining the regular secondary structures that mediate a reversal in the chain direction. Rare beta-turns of type I' and II' are also identified as preferred sites for insertions.

Multidimensional free energy surface of unfolding of HP-36: Microscopic origin of ruggedness

The protein folding funnel paradigm suggests that folding and unfolding proceed as directed diffusion in a multidimensional free energy surface where a multitude of pathways can be traversed during the protein's sojourn from initial to final state. However, finding even a single pathway, with the detail chronicling of intermediates, is an arduous task. In this work we explore the free energy surface of unfolding pathway through umbrella sampling, for a small globular a-helical protein chicken-villin headpiece (HP-36) when the melting of secondary structures is induced by adding DMSO in aqueous solution. We find that the unfolding proceeds through the initial separation or melting of aggregated hydrophobic core that comprises of three phenylalanine residues (Phe7, Phe11, and Phe18). This separation is accompanied by simultaneous melting of the second helix. Unfolding is found to be a multistage process involving crossing of three consecutive minima and two barriers at the initial stage. At a molecular level, Phe18 is observed to reorient itself towards other hydrophobic grooves to stabilize the intermediate states. We identify the configuration of the intermediates and correlate the intermediates with those obtained in our previous works. We also give an estimate of the barriers for different transition states and observe the softening of the barriers with increasing DMSO concentration. We show that higher concentration of DMSO tunes the unfolding pathway by destabilizing the third minimum and stabilizing the second one, indicating the development of a solvent modified, less rugged pathway. The prime outcome of this work is the demonstration that mixed solvents can profoundly transform the nature of the energy landscape and induce unfolding via a modified route. A successful application of Kramer's rate equation correlating the free energy simulation results shows faster rate of unfolding with increasing DMSO concentration. This work perhaps presents the first systematic theoretical study of the effect of a chemical denaturant on the microscopic free energy surface and rates of unfolding of HP-36. (C) 2014 AIP Publishing LLC.

Dissecting pi-helices: sequence, structure and function

A new procedure for the identification of regular secondary structures using a C-alpha trace has identified 659 pi-helices in 3582 protein chains, solved at high resolution. Taking advantage of this significantly expanded database of pi-helices, we have analysed the functional and structural roles of helices and determined the position-wise amino acid propensity within and around them. These helices range from 5 to 18 residues in length with the average twist and rise being 85.2 +/- 7.2 and 1.28 +/- 0.31 angstrom, respectively. A total of 546 (similar to 83%) out of 659 pi-helices occur in conjunction with alpha-helices, with 101 pi-helices being interspersed between two alpha-helices. The majority of interspersed pi-helices were found to be conserved across a large number of structures within a protein family and produce a significant bend in the overall helical segment as well as local distortions in the neighbouring a-helices. The presence of a pi-helical fragment leads to appropriate orientation of the constituent residues, so as to facilitate favourable interactions and also help in proper folding of the protein chain. In addition to intra helical 6 -> 1 N H center dot center dot center dot O hydrogen bonds, pi-helices are also stabilized by several other non-bonded interactions. pi-Helices show distinct positional residue preferences, which are different from those of a-helices.

Muscle wasting in myotonic dystrophies: a model of premature aging

Myotonic dystrophy type 1 (DM1 or Steinert's disease) and type 2 (DM2) are multisystem disorders of genetic origin. Progressive muscular weakness, atrophy and myotonia are the most prominent neuromuscular features of these diseases, while other clinical manifestations such as cardiomyopathy, insulin resistance and cataracts are also common. From a clinical perspective, most DM symptoms are interpreted as a result of an accelerated aging (cataracts, muscular weakness and atrophy, cognitive decline, metabolic dysfunction, etc.), including an increased risk of developing tumors. From this point of view, DM1 could be described as a progeroid syndrome since a notable age dependent dysfunction of all systems occurs. The underlying molecular disorder in DM1 consists of the existence of a pathological (CTG) triplet expansion in the 3' untranslated region (UTR) of the Dystrophia ll/Iyotonica Protein Kinase (DMPK) gene, whereas (CCTG)n repeats in the first intron of the Cellular Nucleic acid Binding Protein/Zinc Finger Protein 9 (CNBP/ZNF9) gene cause DM2. The expansions are transcribed into (CUG)n and (CCUG)n-containing RNA, respectively, which form secondary structures and sequester RNA binding proteins, such as the splicing factor muscleblind-like protein (MBNL), forming nuclear aggregates known as foci. Other splicing factors, such as CUGBP, are also disrupted, leading to a spliceopathy of a large number of downstream genes linked to the clinical features of these diseases. Skeletal muscle regeneration relies on muscle progenitor cells, known as satellite cells, which are activated after muscle damage, and which proliferate and differentiate to muscle cells, thus regenerating the damaged tissue. Satellite cell dysfunction seems to be a common feature of both age-dependent muscle degeneration (sarcopenia) and muscle wasting in DM and other muscle degenerative diseases. This review aims to describe the cellular, molecular and macrostructural processes involved in the muscular degeneration seen in DM patients, highlighting the similarities found with muscle aging.

Maximins S, a novel group of antimicrobial peptides from toad Bombina maxima

Amphibian skin secretions are rich in antimicrobial peptides acting as important components of innate defense system against invading microorganisms. A novel type of peptide, designated as maximin S, was deduced by random sequencing of 793 clones from a constructed Bombina maxima skin cDNA library. The putative primary structures of maximin S peptides can be grouped into five species, in which maximin S I has 14 amino acid residues and the rest of maximin S peptides (S2-S5) all have 18 amino acid residues. Unlike most of the amphibian antimicrobial peptides so far identified, the newly characterized four maximin S precursors are composed of maximin S I and different combinations of tandem repeated maximin S2-S5 linked by internal peptides. Except maximin S I, the predicted secondary structures of maximin S2-S5 show a similar amphipathic alpha-helical structure. MALDI-TOF mass spectrometry analysis of partially isolated skin secretions of the toad indicates that most of the deduced maximin S peptides are expressed. Two deduced maximin S peptides (S1, S4) were synthesized and their antimicrobial activities were tested. Maximin S4 only had an antibiotic activity against mycoplasma and had no antibacterial or antifungal activity toward tested strains. Maximin S1 had no activity under the same conditions. (C) 2004 Elsevier Inc. All rights reserved.

Is there a genetic relationship between mRNA and the encoded protein at three-dimensional level?

A comparative study on the structures of some mRNAs and their encoded proteins shows an intriguing correlation between the two foldings. Non-random distribution of codons in the secondary structures of mRNAs is also shown, which appears to be in accordance with the conformational properties of amino acids in protein structures to some extent. These results seem to suggest that there may be a kind of genetic relationship between mRNA and protein at three-dimensional level.

Cardioviral RNA structure logo analysis: entropy, correlations, and prediction

In recent years, there has been an increased number of sequenced RNAs leading to the development of new RNA databases. Thus, predicting RNA structure from multiple alignments is an important issue to understand its function. Since RNA secondary structures are often conserved in evolution, developing methods to identify covariate sites in an alignment can be essential for discovering structural elements. Structure Logo is a technique established on the basis of entropy and mutual information measured to analyze RNA sequences from an alignment. We proposed an efficient Structure Logo approach to analyze conservations and correlations in a set of Cardioviral RNA sequences. The entropy and mutual information content were measured to examine the conservations and correlations, respectively. The conserved secondary structure motifs were predicted on the basis of the conservation and correlation analyses. Our predictive motifs were similar to the ones observed in the viral RNA structure database, and the correlations between bases also corresponded to the secondary structure in the database.

Bead-probe complex capture a couple of SINE and LINE family from genomes of two closely related species of East Asian cyprinid directly using magnetic separation

Background: Short and long interspersed elements (SINEs and LINEs, respectively), two types of retroposons, are active in shaping the architecture of genomes and powerful tools for studies of phylogeny and population biology. Here we developed special protocol to apply biotin-streptavidin bead system into isolation of interspersed repeated sequences rapidly and efficiently, in which SINEs and LINEs were captured directly from digested genomic DNA by hybridization to bead-probe complex in solution instead of traditional strategy including genomic library construction and screening. Results: A new couple of SINEs and LINEs that shared an almost identical 3'tail was isolated and characterized in silver carp and bighead carp of two closely related species. These SINEs (34 members), designated HAmo SINE family, were little divergent in sequence and flanked by obvious TSD indicated that HAmo SINE was very young family. The copy numbers of this family was estimated to 2 x 10(5) and 1.7 x 10(5) per haploid genome by Real-Time qPCR, respectively. The LINEs, identified as the homologs of LINE2 in other fishes, had a conserved primary sequence and secondary structures of the 3'tail region that was almost identical to that of HAmo SINE. These evidences suggest that HAmo SINEs are active and amplified recently utilizing the enzymatic machinery for retroposition of HAmoL2 through the recognition of higher-order structures of the conserved 42-tail region. We analyzed the possible structures of HAmo SINE that lead to successful amplification in genome and then deduced that HAmo SINE, SmaI SINE and FokI SINE that were similar in sequence each other, were probably generated independently and created by LINE family within the same lineage of a LINE phylogeny in the genomes of different hosts. Conclusion: The presented results show the advantage of the novel method for retroposons isolation and a pair of young SINE family and its partner LINE family in two carp fishes, which strengthened the hypotheses containing the slippage model for initiation of reverse transcription, retropositional parasitism of SINEs on LINEs, the formation of the stem loop structure in 3'tail region of some SINEs and LINEs and the mechanism of template switching in generating new SINE family.

Mutations stabilize small subunit ribosomal RNA in desiccation-tolerant cyanobacteria Nostoc

The ribosomal RNA molecule is an ideal model for evaluating the stability of a gene product under desiccation stress. We isolated 8 Nostoc strains that had the capacity to withstand desiccation in habitats and sequenced their 16S rRNA genes. The stabilities of 16S rRNAs secondary structures, indicated by free energy change of folding, were compared among Nostoc and other related species. The results suggested that 163 rRNA secondary structures of the desiccation-tolerant Nostoc strains were more stable than that of planktonic Nostocaceae species. The stabilizing mutations were divided into two categories: (1) those causing GC to replace other types of base pairs in stems and (2) those causing extension of stems. By mapping stabilizing mutations onto the Nostoc phylogenetic tree based on 16S rRNA gene, it was shown that most of stabilizing mutations had evolved during adaptive radiation among Nostoc spp. The evolution of 16S rRNA along the Nostoc lineage is suggested to be selectively advantageous under desiccation stress.