Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda)

ITS2 sequences are used extensively in molecular taxonomy and population genetics of arthropods and other animals yet little is known about the molecular evolution of ITS2. We studied the secondary structure of ITS2 in species from each of the six main lineages of hard ticks (family Ixodidae). The ITS2 of these ticks varied in length from 679 bp in Ixodes scapularis to 1547 bp in Aponomma concolor. Nucleotide content varied also: the ITS2 of ticks from the Prostriata lineage (Ixodes spp.) had 46-49% GC whereas ITS2 sequences of ticks from the Metastriata lineage (all other hard ticks) had 61-62% GC. Despite variation in nucleotide sequence, the secondary structure of the ITS2 of all of these ticks apparently has five domains. Stems 1, 3, 4 and 5 of this secondary structure were obvious in all of the species studied. However, stem 2 was not always obvious despite the fact that it is flanked by highly conserved sequence motifs in the adjacent stems, stems 1 and 3. The ITS2 of hard ticks has apparently evolved mostly by increases and decreases in length of the nucleotide sequences, which caused increases, and decreases in the length of stems of the secondary structure. This is most obvious when stems of the secondary structures of the Prostriata (Ixodes spp.) are compared to those of the Metastriata (all other hard ticks). Increases in the size of the ITS2 may have been caused by replication slippage which generated large repeats, like those seen in Haemaphysalis humerosa and species from the Rhipicepalinae lineage, and the small repeats found in species from the other lineages of ticks.

Secondary structure of the third extracellular loop responsible for ligand selectivity of a mammalian gonadotropin-releasing hormone receptor

The extracellular loop 3 (ECL3) of the mammalian gonadotropin-releasing hormone receptor (GnRH-R) contains an acidic amino acid (Glu(301) in the mouse GnRH-R,) that confers agonist selectivity for Are in mammalian GnRH. It is proposed that a specific conformation of ECL3 is necessary to orientate the carboxyl side chain of the acidic residue for interaction with Arg(8) of GnRH, which is supported by decreased affinity for Arg(8) GnRH but not Gln(8) GnRH when an adjacent Pro is mutated to Ala. To probe the structural contribution of the loop domain to the proposed presentation of the carboxyl side chain, we synthesized a model peptide (CGPEMLNRVSEPGC) representing residues 293-302 of mouse ECL3, where Cys and Gly residues are added symmetrically at the N and C termini, respectively, allowing the introduction of a disulfide bridge to simulate the distances at which the ECL3 is tethered to the transmembrane domains 6 and 7 of the receptor. The ability of the ECL3 peptide to bind GnRH with low affinity was demonstrated by its inhibition of GnRH stimulation of inositol phosphate production in cells expressing the GnRH-R. The CD bands of the ECL3 peptides exhibited a superposition of predominantly unordered structure and partial contributions from beta-sheet structure. Likewise, the analysis of the amide I and amide III bands from micro-Raman and FT Raman experiments revealed mainly unordered conformations of the cyclic and of the linear peptide. NMR data demonstrated the presence of a beta-hairpin among an ensemble of largely disordered structures in the cyclic peptide. The location of the turn linking the two strands of the hairpin was assigned to the three central residues L-296, N-297, and R-298. A small population of structured species among an ensemble of predominantly random coil conformation suggests that the unliganded receptor represents a variety of structural conformers, some of which have the potential to make contacts with the ligand. We propose a mechanism of receptor activation whereby binding of the agonist to the inactive receptor state induces and stabilizes a particular structural state of the loop domain, leading to further conformational rearrangements across the transmembrane domain and signal propagating interaction with G proteins. Interaction of the Glu(301) of the receptor with Arg(8) of GnRH induces a folded configuration of the ligand. Our proposal thus suggests that conformational changes of both ligand and receptor result from this interaction.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

Secondary structure prediction with support vector machines

Motivation: A new method that uses support vector machines (SVMs) to predict protein secondary structure is described and evaluated. The study is designed to develop a reliable prediction method using an alternative technique and to investigate the applicability of SVMs to this type of bioinformatics problem. Methods: Binary SVMs are trained to discriminate between two structural classes. The binary classifiers are combined in several ways to predict multi-class secondary structure. Results: The average three-state prediction accuracy per protein (Q3) is estimated by cross-validation to be 77.07 ± 0.26% with a segment overlap (Sov) score of 73.32 ± 0.39%. The SVM performs similarly to the 'state-of-the-art' PSIPRED prediction method on a non-homologous test set of 121 proteins despite being trained on substantially fewer examples. A simple consensus of the SVM, PSIPRED and PROFsec achieves significantly higher prediction accuracy than the individual methods. Availability: The SVM classifier is available from the authors. Work is in progress to make the method available on-line and to integrate the SVM predictions into the PSIPRED server.

Benchmarking secondary structure prediction for fold recognition

If secondary structure predictions are to be incorporated into fold recognition methods, an assessment of the effect of specific types of errors in predicted secondary structures on the sensitivity of fold recognition should be carried out. Here, we present a systematic comparison of different secondary structure prediction methods by measuring frequencies of specific types of error. We carry out an evaluation of the effect of specific types of error on secondary structure element alignment (SSEA), a baseline fold recognition method. The results of this evaluation indicate that missing out whole helix or strand elements, or predicting the wrong type of element, is more detrimental than predicting the wrong lengths of elements or overpredicting helix or strand. We also suggest that SSEA scoring is an effective method for assessing accuracy of secondary structure prediction and perhaps may also provide a more appropriate assessment of the “usefulness” and quality of predicted secondary structure, if secondary structure alignments are to be used in fold recognition.

Rapid protein domain assignment from amino acid sequence using predicted secondary structure

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (±20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed

Secondary structure analysis of the dissimilatory sulphite reductase in Desulfovibrio desulfuricans

The complete sequences of the dsrA and dsrB genes coding for the α− and β−subunits, respectively, of the sulphite reductase enzyme in Desulfovibrio desulfuricans were determined. Analyses of the amino acid sequences indicated a number of serohaem/Fe4S4 binding consensus sequences whilst predictive secondary structure analysis revealed a similar pattern of α−helix and β−strand structures between the two subunits which was indicative of gene duplication.

Towards an assessment of character interdependence in avian RNA phylogenetics: A general secondary structure model for the avian mitochondrial 16S rRNA

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

gamma-zein secondary structure in solution by circular dichroism

The proline-rich N-terminal domain of gamma-zein has been reported in relevant process, which include its ability to cross the cell membranes. Evidences indicate that synthetic hexapeptide (PPPVHL), naturally found in N-terminal portion of gamma-zein, can adopt the polyproline II (PPII) conformation in aqueous solution. The secondary structure of gamma-zein in maize protein bodies had been analyzed by solid state Fourier transform infrared and nuclear magnetic resonance spectroscopies. However, it was not possible to measure PPII content in physiological environment since the beta-sheet and PPII signals overlap in both solid state techniques. Here, the secondary structure of gamma-zein has been analyzed by circular dichroism in SDS aqueous solution with and without ditiothreitol (DTT), and in 60% of 2-propanol and water with DTT The results show that gamma-zein has high helical content in all solutions. The PPII conformation was present at about 7% only in water/DTT solution. (c) 2007 Wiley Periodicals, Inc.

Influence of N-Terminus Modifications on the Biological Activity, Membrane Interaction, and Secondary Structure of the Antimicrobial Peptide Hylin-a1

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Electrochemical genosensors for the detection of Bonamia parasite. Selection of single strand-DNA (ssDNA) probes by simulation of the secondary structure folding

A post-PCR nucleic acid work by comparing experimental data, from electrochemical genosensors, and bioinformatics data, derived from the simulation of the secondary structure folding and prediction of hybridisation reaction, was carried out in order to rationalize the selection of ssDNA probes for the detection of two Bonamia species, B. exitiosa and B. ostreae, parasites of Ostrea edulis.Six ssDNA probes (from 11 to 25 bases in length, 2 thiolated and 4 biotinylated) were selected within different regions of B. ostreae and B. exitiosa PCR amplicons (300 and 304 bases, respectively) with the aim to discriminate between these parasite species. ssDNA amplicons and probes were analyzed separately using the "Mfold Web Server" simulating the secondary structure folding behaviour. The hybridisation of amplicon-probe was predicted by means of "Dinamelt Web Server". The results were evaluated considering the number of hydrogen bonds broken and formed in the simulated folding and hybridisation process, variance in gaps for each sequence and number of available bases. In the experimental part, thermally denatured PCR products were captured at the sensor interface via sandwich hybridisation with surface-tethered probes (thiolated probes) and biotinylated signalling probes. A convergence between analytical signals and simulated results was observed, indicating the possibility to use bioinformatic data for ssDNA probes selection to be incorporated in genosensors. (C) 2011 Elsevier B.V. All rights reserved.

Structure-activity relationship of mastoparan analogs: effects of the number and positioning of Lys residues on secondary structure, interaction with membrane-mimetic systems and biological activity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Water-in-Oil Micro-Emulsion Enhances the Secondary Structure of a Protein by Confinement

A scheme is presented in which an organic solvent environment in combination with surfactants is used to confine a natively unfolded protein inside an inverse microemulsion droplet. This type of confinement allows a study that provides unique insight into the dynamic structure of an unfolded, flexible protein which is still solvated and thus under near-physiological conditions. In a model system, the protein osteopontin (OPN) is used. It is a highly phosphorylated glycoprotein that is expressed in a wide range of cells and tissues for which limited structural analysis exists due to the high degree of flexibility and large number of post-translational modifications. OPN is implicated in tissue functions, such as inflammation and mineralisation. It also has a key function in tumour metastasis and progression. Circular dichroism measurements show that confinement enhances the secondary structural features of the protein. Small-angle X-ray scattering and dynamic light scattering show that OPN changes from being a flexible protein in aqueous solution to adopting a less flexible and more compact structure inside the microemulsion droplets. This novel approach for confining proteins while they are still hydrated may aid in studying the structure of a wide range of natively unfolded proteins.