Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

Lyophilization-induced reversible changes in the secondary structure of proteins.

Changes in the secondary structure of some dozen different proteins upon lyophilization of their aqueous solutions have been investigated by means of Fourier-transform infrared spectroscopy in the amide III band region. Dehydration markedly (but reversibly) alters the secondary structure of all the proteins studied, as revealed by both the quantitative analysis of the second derivative spectra and the Gaussian curve fitting of the original infrared spectra. Lyophilization substantially increases the beta-sheet content and lowers the alpha-helix content of all proteins. In all but one case, proteins become more ordered upon lyophilization.

Genetic and comparative analyses reveal an alternative secondary structure in the region of nt 912 of Escherichia coli 16S rRNA.

Mutations at position 912 of Escherichia coli 16S rRNA result in two notable phenotypes. The C-->U transition confers resistance to streptomycin, a translational-error-inducing antibiotic, while a C-->G transversion causes marked retardation of cell growth rate. Starting with the slow-growing G912 mutant, random mutagenesis was used to isolate a second site mutation that restored growth nearly to the wild-type rate. The second site mutation was identified as a G-->C transversion at position 885 in 16S rRNA. Cells containing the G912 mutation had an increased doubling time, abnormal sucrose gradient ribosome/subunit profile, increased sensitivity to spectinomycin, dependence upon streptomycin for growth in the presence of spectinomycin, and slower translation rate, whereas cells with the G912/C885 double mutation were similar to wild type in these assays. Comparative analysis showed there was significant covariation between positions 912 and 885. Thus the second-site suppressor analysis, the functional assays, and the comparative data suggest that the interaction between nt 912 and nt 885 is conserved and necessary for normal ribosome function. Furthermore, the comparative data suggest that the interaction extends to include G885-G886-G887 pairing with C912-U911-C910. An alternative secondary structure element for the central domain of 16S rRNA is proposed.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Periodicity of polar and nonpolar amino acids is the major determinant of secondary structure in self-assembling oligomeric peptides.

The tendency of a polypeptide chain to form alpha-helical or beta-strand secondary structure depends upon local and nonlocal effects. Local effects reflect the intrinsic propensities of the amino acid residues for particular secondary structures, while nonlocal effects reflect the positioning of the individual residues in the context of the entire amino acid sequence. In particular, the periodicity of polar and nonpolar residues specifies whether a given sequence is consistent with amphiphilic alpha-helices or beta-strands. The importance of intrinsic propensities was compared to that of polar/nonpolar periodicity by a direct competition. Synthetic peptides were designed using residues with intrinsic propensities that favored one or the other type of secondary structure. The polar/nonpolar periodicities of the peptides were designed either to be consistent with the secondary structure favored by the intrinsic propensities of the component residues or in other cases to oppose these intrinsic propensities. Characterization of the synthetic peptides demonstrated that in all cases the observed secondary structure correlates with the periodicity of the peptide sequence--even when this secondary structure differs from that predicted from the intrinsic propensities of the component amino acids. The observed secondary structures are concentration dependent, indicating that oligomerization of the amphiphilic peptides is responsible for the observed secondary structures. Thus, for self-assembling oligomeric peptides, the polar/nonpolar periodicity can overwhelm the intrinsic propensities of the amino acid residues and serves as the major determinant of peptide secondary structure.

Identifying sequence regions undergoing conformational change via predicted continuum secondary structure

Motivation: Conformational flexibility is essential to the function of many proteins, e.g. catalytic activity. To assist efforts in determining and exploring the functional properties of a protein, it is desirable to automatically identify regions that are prone to undergo conformational changes. It was recently shown that a probabilistic predictor of continuum secondary structure is more accurate than categorical predictors for structurally ambivalent sequence regions, suggesting that such models are suited to characterize protein flexibility. Results: We develop a computational method for identifying regions that are prone to conformational change directly from the amino acid sequence. The method uses the entropy of the probabilistic output of an 8-class continuum secondary structure predictor. Results for 171 unique amino acid sequences with well-characterized variable structure (identified in the 'Macromolecular movements database') indicate that the method is highly sensitive at identifying flexible protein regions, but false positives remain a problem. The method can be used to explore conformational flexibility of proteins (including hypothetical or synthetic ones) whose structure is yet to be determined experimentally.

Prediction of protein continuum secondary structure with probabilistic models based on NMR solved structures

Background: The structure of proteins may change as a result of the inherent flexibility of some protein regions. We develop and explore probabilistic machine learning methods for predicting a continuum secondary structure, i.e. assigning probabilities to the conformational states of a residue. We train our methods using data derived from high-quality NMR models. Results: Several probabilistic models not only successfully estimate the continuum secondary structure, but also provide a categorical output on par with models directly trained on categorical data. Importantly, models trained on the continuum secondary structure are also better than their categorical counterparts at identifying the conformational state for structurally ambivalent residues. Conclusion: Cascaded probabilistic neural networks trained on the continuum secondary structure exhibit better accuracy in structurally ambivalent regions of proteins, while sustaining an overall classification accuracy on par with standard, categorical prediction methods.

Preliminary crystallographic characterization of an RNA helicase from Kunjin virus

Kunjin virus is a member of the Flavivirus genus and is an Australian variant of West Nile virus. The C-terminal domain of the Kunjin virus NS3 protein displays helicase activity. The protein is thought to separate daughter and template RNA strands, assisting the initiation of replication by unwinding RNA secondary structure in the 3' nontranslated region. Expression, purification and preliminary crystallographic characterization of the NS3 helicase domain are reported. It is shown that Kunjin virus helicase may adopt a dimeric assembly in absence of nucleic acids, oligomerization being a means to provide the helicases with multiple nucleic acid-binding capability, facilitating translocation along the RNA strands. Kunjin virus NS3 helicase domain is an attractive model for studying the molecular mechanisms of flavivirus replication, while simultaneously providing a new basis for the rational development of anti-flaviviral compounds.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

Background: Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei. Results: The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infraordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes. Conclusions: Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Représentation et recherche de motifs cycliques et structuraux d’ARN connus dans les structures secondaires

L'acide désoxyribonucléique (ADN) et l'acide ribonucléique (ARN) sont des polymères de nucléotides essentiels à la cellule. À l'inverse de l'ADN qui sert principalement à stocker l'information génétique, les ARN sont impliqués dans plusieurs processus métaboliques. Par exemple, ils transmettent l’information génétique codée dans l’ADN. Ils sont essentiels pour la maturation des autres ARN, la régulation de l’expression génétique, la prévention de la dégradation des chromosomes et le ciblage des protéines dans la cellule. La polyvalence fonctionnelle de l'ARN résulte de sa plus grande diversité structurale. Notre laboratoire a développé MC-Fold, un algorithme pour prédire la structure des ARN qu'on représente avec des graphes d'interactions inter-nucléotidiques. Les sommets de ces graphes représentent les nucléotides et les arêtes leurs interactions. Notre laboratoire a aussi observé qu'un petit ensemble de cycles d'interactions à lui seul définit la structure de n'importe quel motif d'ARN. La formation de ces cycles dépend de la séquence de nucléotides et MC-Fold détermine les cycles les plus probables étant donnée cette séquence. Mon projet de maîtrise a été, dans un premier temps, de définir une base de données des motifs structuraux et fonctionnels d'ARN, bdMotifs, en terme de ces cycles. Par la suite, j’ai implanté un algorithme, MC-Motifs, qui recherche ces motifs dans des graphes d'interactions et, entre autres, ceux générés par MC-Fold. Finalement, j’ai validé mon algorithme sur des ARN dont la structure est connue, tels que les ARN ribosomaux (ARNr) 5S, 16S et 23S, et l'ARN utilisé pour prédire la structure des riborégulateurs. Le mémoire est divisé en cinq chapitres. Le premier chapitre présente la structure chimique, les fonctions cellulaires de l'ARN et le repliement structural du polymère. Dans le deuxième chapitre, je décris la base de données bdMotifs. Dans le troisième chapitre, l’algorithme de recherche MC-Motifs est introduit. Le quatrième chapitre présente les résultats de la validation et des prédictions. Finalement, le dernier chapitre porte sur la discussion des résultats suivis d’une conclusion sur le travail.

Riboswitches : le cas des atténuateurs de la transcription du type terminateur/antiterminateur chez les bactéries

Il est essentiel pour chaque organisme d’avoir la possibilité de réguler ses fonctions afin de permettre sa survie et d’améliorer sa capacité de se reproduire en divers habitats. Avec l’information disponible, il semble que les organismes consacrent une partie assez importante de leur matériel génétique à des fonctions de régulation. On peut envisager que certains mécanismes de régulation ont persisté dans le temps parce qu’ils remplissent bien leurs rôles. Les premières études sur les procaryotes ont indiqué qu’il y avait peu de mécanismes de régulation exerçant le contrôle des gènes, mais il a été démontré par la suite qu’une variété de ces mécanismes est utilisée pour la régulation de gènes et d’opérons. En particulier, les opérons bactériens impliqués dans la biosynthèse des acides aminés, l’ARNt synthétase, la dégradation des acides aminés, les protéines ribosomales et l’ARN ribosomal font l’objet d’un contrôle par l’atténuation de la transcription. Ce mécanisme d’atténuation de la transcription diffère d’autres mécanismes pour la génération de deux structures différentes de l’ARNm, où l’une de ces structures réprime le gène en aval, et l’autre permet de continuer la transcription/traduction. Dans le cadre de cette recherche, nous nous sommes intéressé au mécanisme d’atténuation de la transcription chez les procaryotes où aucune molécule ne semble intervenir comme facteur de régulation, en me concentrant sur la régulation des opérons bactériens. Le but principal de ce travail est de présenter une nouvelle méthode de recherche des riborégulateurs qui combine la recherche traditionnelle des riborégulateurs avec la recherche structurale. En incorporant l’étude du repliement de l’ARNm, nous pouvons mieux identifier les atténuateurs répondant à ce type de mécanisme d’atténuation. Ce mémoire est divisé en quatre chapitres. Le premier chapitre présente une revue de la littérature sur l’ARN et un survol sur les mécanismes de régulation de l’expression génétique chez les procaryotes. Les chapitres 2 et 3 sont consacrés à la méthodologie utilisée dans cette recherche et à l’implémentation du logiciel TA-Search. Enfin, le chapitre 4 expose les conclusions et les applications potentielles de la méthode.

Studies on Anopheles (Kerteszia) homunculus Komp (Diptera: Culicidae)

The present findings suggest that Anopheles (Kerteszia) homunculus may comprise more than one species. The rDNA ITS2 sequence data corroborate the presence of An. homunculus l.s. in Mata Atlantica, southern Brazil, and suggest that specimens from Trinidad may belong to an unnamed morphologically similar species. There is a need for additional studies to establish the geographical distribution of An. homunculus l.s. in continental South America and in Trinidad, especially in southern Mata Atlantica, Brazil.

Patterns and processes of diversification in a widespread and ecologically diverse avian group, the buteonine hawks (Aves, Accipitridae)

Buteonine hawks represent one of the most diverse groups in the Accipitridae, with 58 species distributed in a variety of habitats on almost all continents. Variations in migratory behavior, remarkable dispersal capability, and unusual diversity in Central and South America make buteonine hawks an excellent model for studies in avian evolution. To evaluate the history of their global radiation, we used an integrative approach that coupled estimation of the phylogeny using a large sequence database (based on 6411 bp of mitochondrial markers and one nuclear intron from 54 species), divergence time estimates, and ancestral state reconstructions. Our findings suggest that Neotropical buteonines resulted from a long evolutionary process that began in the Miocene and extended to the Pleistocene. Colonization of the Nearctic, and eventually the Old World, occurred from South America, promoted by the evolution of seasonal movements and development of land bridges. Migratory behavior evolved several times and may have contributed not only to colonization of the Holarctic, but also derivation of insular species. In the Neotropics, diversification of the buteonines included four disjunction events across the Andes. Adaptation of monophyletic taxa to wet environments occurred more than once, and some relationships indicate an evolutionary connection among mangroves, coastal and varzea environments. On the other hand, groups occupying the same biome, forest, or open vegetation habitats are not monophyletic. Refuges or sea-level changes or a combination of both was responsible for recent speciation in Amazonian taxa. In view of the lack of concordance between phylogeny and classification, we propose numerous taxonomic changes. (C) 2009 Elsevier Inc. All rights reserved.