Modeling protein evolution using secondary structures

L’évolution des protéines est un domaine important de la recherche en bioinformatique et catalyse l'intérêt de trouver des outils d'alignement qui peuvent être utilisés de manière fiable et modéliser avec précision l'évolution d'une famille de protéines. TM-Align (Zhang and Skolnick, 2005) est considéré comme l'outil idéal pour une telle tâche, en termes de rapidité et de précision. Par conséquent, dans cette étude, TM-Align a été utilisé comme point de référence pour faciliter la détection des autres outils d'alignement qui sont en mesure de préciser l'évolution des protéines. En parallèle, nous avons élargi l'actuel outil d'exploration de structures secondaires de protéines, Helix Explorer (Marrakchi, 2006), afin qu'il puisse également être utilisé comme un outil pour la modélisation de l'évolution des protéines.

Origin and evolution of 3'UTR of flaviviruses: long direct repeats as a basis for the formation of secondary structures and their significance for virus transmission

The 3' untranslated regions (3'UTRs) of flaviviruses are reviewed and analyzed in relation to short sequences conserved as direct repeats (DRs). Previously, alignments of the 3'UTRs have been constructed for three of the four recognized flavivirus groups, namely mosquito-borne, tick-borne, and nonclassified flaviviruses (MBFV, TBFV, and NCFV, respectively). This revealed (1) six long repeat sequences (LRSs) in the 3'UTR and open-reading frame (ORF) of the TBFV, (2) duplication of the 3'UTR of the NCFV by intramolecular recombination, and (3) the possibility of a common origin for all DRs within the MBFV. We have now extended this analysis and review it in the context of all previous published analyses. This has been achieved by constructing a robust alignment between all flaviviruses using the published DRs and secondary RNA structures as "anchors" to reveal additional homologies along the 3'UTR. This approach identified nucleotide regions within the MBFV, NKV (no-known vector viruses), and NCFV 3'UTRs that are homologous to different LRSs in the TBFV 3'UTR and ORF. The analysis revealed that some of the DRs and secondary RNA structures described individually within each flavivirus group share common evolutionary origins. The 3'UTR of flaviviruses, and possibly the ORF, therefore probably evolved through multiple duplication of an RNA domain, homologous to the LRS previously identified only in the TBFV. The short DRs in all virus groups appear to represent the evolutionary remnants of these domains rather than resulting from new duplications. The relevance of these flavivirus DRs to evolution, diversity, 3'UTR enhancer function, and virus transmission is reviewed.

Investigating the effect of different positioning of lysine residues along the peptide chain of mastoparans for their secondary structures and biological activities

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

De Novo Design of secondary structures: Synthesis and conformational studies

Chemists have long sought to extrapolate the power of biological catalysis and recognition to synthetic systems. These efforts have focused largely on low molecular weight catalysts and receptors; however, biological systems themselves rely almost exclusively on polymers, proteins and RNA, to perform complex chemical functions. Proteins and RNA are unique in their ability to adopt compact, well-ordered conformations, and specific folding provides precise spatial orientation of the functional groups that comprise the “active site”. These features suggest that identification of new polymer backbones with discrete and predictable folding propensities (“foldamers”) will provide a basis for design of molecular machines with unique capabilities. The foldamer approach complements current efforts to design unnatural properties into polypeptides and polynucleotides. The aim of this thesis is the synthesis and conformational studies of new classes of foldamers, using a peptidomimetic approach. Moreover their attitude to be utilized as ionophores, catalysts, and nanobiomaterials were analyzed in solution and in the solid state. This thesis is divided in thematically chapters that are reported below. It begins with a very general introduction (page 4) which is useful, but not strictly necessary, to the expert reader. It is worth mentioning that paragraph I.3 (page 22) is the starting point of this work and paragraph I.5 (page 32) isrequired to better understand the results of chapters 4 and 5. In chapter 1 (page 39) is reported the synthesis and conformational analysis of a novel class of foldamers containing (S)-β3-homophenylglycine [(S)-β3-hPhg] and D- 4-carboxy-oxazolidin-2-one (D-Oxd) residues in alternate order is reported. The experimental conformational analysis performed in solution by IR, 1HNMR, and CD spectroscopy unambiguously proved that these oligomers fold into ordered structures with increasing sequence length. Theoretical calculations employing ab initio MO theory suggest a helix with 11-membered hydrogenbonded rings as the preferred secondary structure type. The novel structures enrich the field of peptidic foldamers and might be useful in the mimicry of native peptides. In chapter 2 cyclo-(L-Ala-D-Oxd)3 and cyclo-(L-Ala-DOxd) 4 were prepared in the liquid phase with good overall yields and were utilized for bivalent ions chelation (Ca2+, Mg2+, Cu2+, Zn2+ and Hg2+); their chelation skill was analyzed with ESI-MS, CD and 1HNMR techniques and the best results were obtained with cyclo-(L-Ala-D-Oxd)3 and Mg2+ or Ca2+. Chapter 3 describes an application of oligopeptides as catalysts for aldol reactions. Paragraph 3.1 concerns the use of prolinamides as catalysts of the cross aldol addition of hydroxyacetone to aromatic aldeydes, whereas paragraphs 3.2 and 3.3 are about the catalyzed aldol addition of acetone to isatins. By means of DFT and AIM calculations, the steric and stereoelectronic effects that control the enantioselectivity in the cross-aldol addition of acetone to isatin catalysed by L-proline have been studied, also in the presence of small quantities of water. In chapter 4 is reported the synthesis and the analysis of a new fiber-like material, obtained from the selfaggregation of the dipeptide Boc-L-Phe-D-Oxd-OBn, which spontaneously forms uniform fibers consisting of parallel infinite linear chains arising from singleintermolecular N-H···O=C hydrogen bonds. This is the absolute borderline case of a parallel β-sheet structure. Longer oligomers of the same series with general formula Boc-(L-Phe-D-Oxd)n-OBn (where n = 2-5), are described in chapter 5. Their properties in solution and in the solid state were analyzed, in correlation with their attitude to form intramolecular hydrogen bond. In chapter 6 is reported the synthesis of imidazolidin-2- one-4-carboxylate and (tetrahydro)-pyrimidin-2-one-5- carboxylate, via an efficient modification of the Hofmann rearrangement. The reaction affords the desired compounds from protected asparagine or glutamine in good to high yield, using PhI(OAc)2 as source of iodine(III).

Comparative analyses of the secondary structures of synthetic and intracellular yeast MFA2 mRNAs

The overall folded (global) structure of mRNA may be critical to translation and turnover control mechanisms, but it has received little experimental attention. Presented here is a comparative analysis of the basic features of the global secondary structure of a synthetic mRNA and the same intracellular eukaryotic mRNA by dimethyl sulfate (DMS) structure probing. Synthetic MFA2 mRNA of Saccharomyces cerevisiae first was examined by using both enzymes and chemical reagents to determine single-stranded and hybridized regions; RNAs with and without a poly(A) tail were compared. A folding pattern was obtained with the aid of the mfold program package that identified the model that best satisfied the probing data. A long-range structural interaction involving the 5′ and 3′ untranslated regions and causing a juxtaposition of the ends of the RNA, was examined further by a useful technique involving oligo(dT)-cellulose chromatography and antisense oligonucleotides. DMS chemical probing of A and C nucleotides of intracellular MFA2 mRNA was then done. The modification data support a very similar intracellular structure. When low reactivity of A and C residues is found in the synthetic RNA, ≈70% of the same sites are relatively more resistant to DMS modification in vivo. A slightly higher sensitivity to DMS is found in vivo for some of the A and C nucleotides predicted to be hybridized from the synthetic structural model. With this small mRNA, the translation process and mRNA-binding proteins do not block DMS modifications, and all A and C nucleotides are modified the same or more strongly than with the synthetic RNA.

Function of RNA secondary structures in transcriptional attenuation of the Bacillus subtilis pyr operon

The Bacillus subtilis pyr operon is regulated by exogenous pyrimidines by a transcriptional attenuation mechanism. Transcription in vitro from pyr DNA templates specifying attenuation regions yielded terminated and read-through transcripts of the expected lengths. Addition of the PyrR regulatory protein plus UMP led to greatly increased termination. Synthetic antisense deoxyoligonucleotides were used to probe possible secondary structures in the pyr mRNA that were proposed to play roles in controlling attenuation. Oligonucleotides predicted to disrupt terminator structures suppressed termination, whereas oligonucleotides predicted to disrupt the stem of antiterminator stem-loops strongly promoted termination at the usual termination site. Oligonucleotides that disrupt a previously unrecognized stem-loop structure, called the anti-antiterminator, the formation of which interferes with formation of the downstream antiterminator, suppressed termination. We propose that transcriptional attenuation of the pyr operon is governed by switching between alternative antiterminator versus anti-antiterminator plus terminator structures, and that PyrR acts by UMP-dependent binding to and stabilization of the anti-antiterminator.

Targeting nucleic acid secondary structures by antisense oligonucleotides designed through in vitro selection.

Using an in vitro selection approach, we have isolated oligonucleotides that can bind to a DNA hairpin structure. Complex formation of these oligonucleotides with the target hairpin involves some type of triple-stranded structure with noncanonical interaction, as indicated by bandshift assays and footprinting studies. The selected oligomers can block restriction endonuclease cleavage of the target hairpin in a sequence-specific manner. We demonstrate that in vitro selection can extend the antisense approach to functional targeting of secondary structure motifs. This could provide a basis for interfering with regulatory processes mediated by a variety of nucleic acid structures.

Analysis of the relation between the sequence and secondary and three-dimensional structures of immunoglobulin molecules.

Methods of structural and statistical analysis of the relation between the sequence and secondary and three-dimensional structures are developed. About 5000 secondary structures of immunoglobulin molecules from the Kabat data base were predicted. Two statistical analyses of amino acids reveal 47 universal positions in strands and loops. Eight universally conservative positions out of the 47 are singled out because they contain the same amino acid in > 90% of all chains. The remaining 39 positions, which we term universally alternative positions, were divided into five groups: hydrophobic, charged and polar, aromatic, hydrophilic, and Gly-Ala, corresponding to the residues that occupied them in almost all chains. The analysis of residue-residue contacts shows that the 47 universal positions can be distinguished by the number and types of contacts. The calculations of contact maps in the 29 antibody structures revealed that residues in 24 of these 47 positions have contacts only with residues of antiparallel beta-strands in the same beta-sheet and residues in the remaining 23 positions always have far-away contacts with residues from other beta-sheets as well. In addition, residues in 6 of the 47 universal positions are also involved in interactions with residues of the other variable or constant domains.

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.

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.

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

Molecular archaeology of Flaviviridae untranslated regions: duplicated RNA structures in the replication enhancer of flaviviruses and pestiviruses emerged via convergent evolution.

RNA secondary structures in the 3'untranslated regions (3'UTR) of the viruses of the family Flaviviridae, previously identified as essential (promoters) or beneficial (enhancers) for replication, have been analysed. Duplicated enhancer elements are revealed as a global feature in the evolution of the 3'UTR of distantly related viruses within the genera Flavivirus and Pestivirus. For the flaviviruses, duplicated structures occur in the 3'UTR of all four distantly related ecological virus subgroups (tick-borne, mosquito-borne, no known vector and insect-specific flaviviruses (ISFV). RNA structural differences distinguish tick-borne flaviviruses with discrete pathogenetic characteristics. For Aedes- and Culex-associated ISFV, secondary RNA structures with different conformations display numerous short ssRNA direct repeats, exposed as loops and bulges. Long quadruplicate regions comprise almost the entire 3'UTR of Culex-associated ISFV. Extended duplicated sequence and associated RNA structures were also discovered in the 3'UTR of pestiviruses. In both the Flavivirus and Pestivirus genera, duplicated RNA structures were localized to the enhancer regions of the 3'UTR suggesting an adaptive role predominantly in wild-type viruses. We propose sequence reiteration might act as a scaffold for dimerization of proteins involved in assembly of viral replicase complexes. Numerous nucleotide repeats exposed as loops/bulges might also interfere with host immune responses acting as a molecular sponge to sequester key host proteins or microRNAs.