Smoothness within ruggedness: the role of neutrality in adaptation.

RNA secondary structure folding algorithms predict the existence of connected networks of RNA sequences with identical structure. On such networks, evolving populations split into subpopulations, which diffuse independently in sequence space. This demands a distinction between two mutation thresholds: one at which genotypic information is lost and one at which phenotypic information is lost. In between, diffusion enables the search of vast areas in genotype space while still preserving the dominant phenotype. By this dynamic the success of phenotypic adaptation becomes much less sensitive to the initial conditions in genotype space.

How are close residues of protein structures distributed in primary sequence?

Structurally neighboring residues are categorized according to their separation in the primary sequence as proximal (1-4 positions apart) and otherwise distal, which in turn is divided into near (5-20 positions), far (21-50 positions), very far ( > 50 positions), and interchain (from different chains of the same structure). These categories describe the linear distance histogram (LDH) for three-dimensional neighboring residue types. Among the main results are the following: (i) nearest-neighbor hydrophobic residues tend to be increasingly distally separated in the linear sequence, thus most often connecting distinct secondary structure units. (ii) The LDHs of oppositely charged nearest-neighbors emphasize proximal positions with a subsidiary maximum for very far positions. (iii) Cysteine-cysteine structural interactions rarely involve proximal positions. (iv) The greatest numbers of interchain specific nearest-neighbors in protein structures are composed of oppositely charged residues. (v) The largest fraction of side-chain neighboring residues from beta-strands involves near positions, emphasizing associations between consecutive strands. (vi) Exposed residue pairs are predominantly located in proximal linear positions, while buried residue pairs principally correspond to far or very far distal positions. The results are principally invariant to protein sizes, amino acid usages, linear distance normalizations, and over- and underrepresentations among nearest-neighbor types. Interpretations and hypotheses concerning the LDHs, particularly those of hydrophobic and charged pairings, are discussed with respect to protein stability and functionality. The pronounced occurrence of oppositely charged interchain contacts is consistent with many observations on protein complexes where multichain stabilization is facilitated by electrostatic interactions.

Mg-SINE: a short interspersed nuclear element from the rice blast fungus, Magnaporthe grisea.

A short interspersed nuclear element, Mg-SINE, was isolated and characterized from the genome of the rice blast fungus, Magnaporthe grisea. Mg-SINE was isolated as an insertion element within Pot2, an inverted-repeat transposon from M. grisea and shows typical features of a mammalian SINE. Mg-SINE is present as a 0.47-kb interspersed sequence at approximately 100 copies per haploid genome in both rice and non-rice isolates of M. grisea, indicating a common evolutionary origin. Secondary structure analysis of Mg-SINE revealed a tRNA-related region at the 5' end which folds into a cloverleaf structure. Genomic fusions resulting in chimeric Mg-SINEs (Ch-SINEs) composed of a sequence homologous to Mg-SINE at the 3' end and an unrelated sequence at its 5' end were also isolated, indicating that this and other DNA rearrangements mediated by these elements may have a major effect on the genomic architecture of this fungus.

Amyloid beta peptide potentiates cytokine secretion by interleukin-1 beta-activated human astrocytoma cells.

Neurodegenerative processes in Alzheimer disease (AD) are thought to be driven in part by the deposition of amyloid beta (A beta), a 39- to 43-amino acid peptide product resulting from an alternative cleavage of amyloid precursor protein. Recent descriptions of in vitro neurotoxic effects of A beta support this hypothesis and suggest toxicity might be mediated by A beta-induced neuronal calcium disregulation. In addition, it has been reported that "aging" A beta results in increased toxic potency due to peptide aggregation and formation of a beta-sheet secondary structure. In addition, A beta might also promote neuropathology indirectly by activating immune/inflammatory pathways in affected areas of the brain (e.g., cortex and hippocampus). Here we report that A beta can modulate cytokine secretion [interleukins 6 and 8 (IL-6 and IL-8)] from human astrocytoma cells (U-373 MG). Freshly prepared and aged A beta modestly stimulated IL-6 and IL-8 secretion from U-373 MG cells. However, in the presence of interleukin-1 beta (IL-1 beta), aged, but not fresh, A beta markedly potentiated (3- to 8-fold) cytokine release. In contrast, aged A beta did not potentiate substance P (NK-1)- or histamine (H1)-stimulated cytokine production. Further studies showed that IL-1 beta-induced cytokine release was potentiated by A beta-(25-35), while A beta-(1-16) was inactive. Calcium disregulation may be responsible for the effects of A beta on cytokine production, since the calcium ionophore A23187 similarly potentiated IL-1 beta-induced cytokine secretion and EGTA treatment blocked either A beta or A23187 activity. Thus, chronic neurodegeneration in AD-affected brain regions may be mediated in part by the ability of A beta to exacerbate inflammatory pathways in a conformation-dependent manner.

Stopped-flow NMR spectroscopy: real-time unfolding studies of 6-19F-tryptophan-labeled Escherichia coli dihydrofolate reductase.

Escherichia coli dihydrofolate reductase (DHFR; EC 1.5.1.3) contains five tryptophan residues that have been replaced with 6-19F-tryptophan. The 19F NMR assignments are known in the native, unliganded form and the unfolded form. We have used these assignments with stopped-flow 19F NMR spectroscopy to investigate the behavior of specific regions of the protein in real time during urea-induced unfolding. The NMR data show that within 1.5 sec most of the intensities of the native 19F resonances of the protein are lost but only a fraction (approximately 20%) of the intensities of the unfolded resonances appears. We postulate that the early disappearance of the native resonances indicates that most of the protein rapidly forms an intermediate in which the side chains have considerable mobility. Stopped-flow far-UV circular dichroism measurements indicate that this intermediate retains native-like secondary structure. Eighty percent of the intensities of the NMR resonances assigned to the individual tryptophans in the unfolded state appear with similar rate constants (k approximately 0.14 sec-1), consistent with the major phase of unfolding observed by stopped-flow circular dichroism (representing 80% of total amplitude). These data imply that after formation of the intermediate, which appears to represent an expanded structural form, all regions of the protein unfold at the same rate. Stopped-flow measurements of the fluorescence and circular dichroism changes associated with the urea-induced unfolding show a fast phase (half-time of about 1 sec) representing 20% of the total amplitude in addition to the slow phase mentioned above. The NMR data show that approximately 20% of the total intensity for each of the unfolded tryptophan resonances is present at 1.5 sec, indicating that these two phases may represent the complete unfolding of the two different populations of the native protein.

Prediction of protein folding class using global description of amino acid sequence.

We present a method for predicting protein folding class based on global protein chain description and a voting process. Selection of the best descriptors was achieved by a computer-simulated neural network trained on a data base consisting of 83 folding classes. Protein-chain descriptors include overall composition, transition, and distribution of amino acid attributes, such as relative hydrophobicity, predicted secondary structure, and predicted solvent exposure. Cross-validation testing was performed on 15 of the largest classes. The test shows that proteins were assigned to the correct class (correct positive prediction) with an average accuracy of 71.7%, whereas the inverse prediction of proteins as not belonging to a particular class (correct negative prediction) was 90-95% accurate. When tested on 254 structures used in this study, the top two predictions contained the correct class in 91% of the cases.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

Catalytically critical nucleotide in domain 5 of a group II intron.

Domain 5 (D5) is a small hairpin structure within group II introns. A bimolecular assay system depends on binding by D5 to an intron substrate for self-splicing activity. In this study, mutations in D5 identify two among six nearly invariant nucleotides as being critical for 5' splice junction hydrolysis but unimportant for binding. A mutation at another site in D5 blocks binding. Thus, mutations can distinguish two D5 functions: substrate binding and catalysis. The secondary structure of D5 may resemble helix I formed by the U2 and U6 small nuclear RNAs in the eukaryotic spliceosome. Our results support a revision of the previously proposed correspondence between D5 and helix I on the basis of the critical trinucleotide 5'-AGC-3' present in both. We suggest that this trinucleotide plays a similar role in promoting the chemical reactions for both splicing systems.

On the distribution of amino acid residues in transmembrane alpha-helix bundles.

The periodic distribution of residues in the sequence of 469 putative transmembrane alpha-helices from eukaryotic plasma membrane polytopic proteins has been analyzed with correlation matrices. The method does not involve any a priori assumption about the secondary structure of the segments or about the physicochemical properties of individual amino acid residues. Maximal correlation is observed at 3.6 residues per period, characteristic of alpha-helices. A scale extracted from the data describes the propensity of the various residues to lie on the same or on opposite helix faces. The most polar face of transmembrane helices, presumably that buried in the protein core, shows a strong enrichment in aromatic residues, while residues likely to face the fatty acyl chains of lipids are largely aliphatic.

Estudos de proteômica, estruturais e funcionais de proteínas envolvidas na degradação da biomassa lignocelulósica: expansinas microbianas e hidrolases de glicosídeos termofílicas

O Brasil possui uma posição privilegiada quando se refere à produção de etanol. Por questões históricas e geográficas o país é responsável por mais de 30 % da produção mundial de etanol, com uma produção nacional de mais de 28 bilhões de litros em 2014. Para maximizar o rendimento desse processo, está em desenvolvimento a tecnologia associada ao etanol de segunda geração ou etanol lignocelulósico. Os principais desafios desta tecnologia são: melhorar a eficiência de conversão do substrato em produto e a produção em grande escala utilizando substratos de baixo custo. Com o objetivo de melhorar a eficiência do processo de conversão foram estudadas proteínas auxiliares (expansinas) que, em conjunto com celulases, melhoram a despolimerização de biomassa lignocelulósica em açúcares fermentescíveis. Além disso, realizou-se também a caracterização de enzimas ativas de carboidratos (CAZymes) de origem termofílica do organismo Thermogemmatispora sp. T81, devido a capacidade que estas proteínas apresentam de manter a atividade e conformação estrutural em altas temperaturas por um prolongado período de tempo. A partir de análises utilizando bioinformática, os genes que codificam para expansinas de Xanthomonas campestris, Bacillus licheniformis e Trichoderma reesei foram clonados e expressos em E. coli, e seus produtos gênicos (as expansinas) tiveram seus índices de sinergismo (devido atuação conjunta com coquetéis comerciais) e atividade catalítica determinados. Adicionalmente, dispondo de alinhamentos estruturais, foi proposto um mecanismo hidrolítico para elas. Em relação à bactéria Thermogemmatispora sp. T81, foram realizadas análises genômicas e proteômicas, a fim de selecionar enzimas superexpressas em meio celulósico. Seus genes foram clonados heterologamente em E. coli e o produto de expressão caracterizado bioquimicamente (cromatografia, ensaios de atividade e perfil de hidrólise) e estruturalmente (SAXS e dicroísmo circular). Os índices de sinergismo determinados foram de 2,47; 1,96 e 2,44 para as expansinas de Xanthomonas campestris, Bacillus licheniformis e Trichoderma reesei, respectivamente. A partir dos alinhamentos estruturais foi proposto a díade Asp/Glu como sitio catalítico em expansinas. As análises de proteômica possibilitaram a seleção de quatro alvos de clonagem, por apresentarem alto índice de expressão quando a bactéria foi cultivada em meio celulósico. Estas proteínas foram caracterizadas quanto a atividade e apresentaram um perfil comum: temperatura ótima de ação (de 70 a 75 °C), pH ótimo de 5, e hidrolisam preferencialmente substratos hemicelulósicos (xilano). A porcentagem de estruturais secundárias das proteínas em estudo foram confirmadas com predições teóricas ao se utilizar a técnica de dicroísmo circular. Desta maneira, os objetivos iniciais propostos neste projeto foram concluídos com a determinação do grau de sinergismo das proteínas expansinas em estudo e a proposição de um mecanismo de hidrólise para as mesmas, considerando que tais proteínas por mais de 20 anos tiveram sua atividade definida exclusivamente como acessória. Além disso, este estudo contribui com a identificação e seleção de genes para CAZymes termofilícas com aplicação biotecnológica devido às propriedades termoestáveis apresentadas.

Identification of the essential histidine residue for high-affinity binding of AlbA protein to albicidin antibiotics

The albA gene from Klebsiella oxytoca encodes a protein that binds albicidin phytotoxins and antibiotics with high affinity. Previously, it has been shown that shifting pH from 6 to 4 reduces binding activity of AlbA by about 30%, indicating that histidine residues might be involved in substrate binding. In this study, molecular analysis of the albA coding region revealed sequence discrepancies with the albA sequence reported previously, which were probably due to sequencing errors. The albA gene was subsequently cloned from K oxytoca ATCC 13182(T) to establish the revised sequence. Biochemical and molecular approaches were used to determine the functional role of four histidine residues (His(78), HiS(125), HiS(141) and His(189)) in the corrected sequence for AlbA. Treatment of AlbA with diethyl pyrocarbonate (DEPC), a histidine-specific alkylating reagent, reduced binding activity by about 95%. DEPC treatment increased absorbance at 240-244 nm by an amount indicating conversion to N-carbethoxyhistidine of a single histidine residue per AlbA molecule. Pretreatment with albicidin protected AlbA against modification by DEPC, with a 1 : 1 molar ratio of albicidin to the protected histidine residues. Based on protein secondary structure and amino acid surface probability indices, it is predicted that HiS125 might be the residue required for albicidin binding. Mutation of HiS125 to either alanine or leucine resulted in about 32% loss of binding activity, and deletion of HiS125 totally abolished binding activity. Mutation of HiS125 to arginine and tyrosine had no effect. These results indicate that HiS125 plays a key role either in an electrostatic interaction between AlbA and albicidin or in the conformational dynamics of the albicidin-binding site.

Alternative polyadenylation and splicing of mRNAs transcribed from the human Sin1 gene

Limited but significant sequence similarity has been observed between an uncharacterized human protein, SIN1, and the S. pombe SIN1, Dictyostelium RIP3 and S. cerevisiae AVO1 proteins. The human Sin1 gene has been automatically predicted (MAPKAP1; GenBank accession number NM_024117); however, this sequence appears to be incomplete. In this study, we have cloned and characterized the full-length human Sin1 mRNA and identified a highly conserved domain that defines the family of SIN1 orthologues, members of which are widely distributed in the fungal and metazoan kingdoms. We demonstrate that Sin1 transcripts can use alternative polyadenylation signals and describe a number of Sin1 splice variants that potentially encode functionally different isoforms. (C) 2004 Elsevier B.V. All rights reserved.

Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA Microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design microarray probes with high specificities and sensitivities. Highly specific probes correspond to probes having unique DNA sequences; whereas highly sensitive probes correspond to those with melting temperature within a desired range and having no secondary structure. The selection of these probes from a set of functional DNA sequences (exons) constitutes a computationally expensive discrete non-linear search problem. We delegate the search task to a simple yet effective Evolution Strategy algorithm. The computational efficiency is also greatly improved by making use of an available bioinformatics tool.

Metal Clips Induce Folding of a Short Unstructured Peptide into an a-Helix via Turn Conformations in Water. Kinetic versus Thermodynamic Products

Short peptides corresponding to two to four a-helical turns of proteins are not thermodynamically stable helices in water. Unstructured octapeptide Ac-His1*-Ala2-Ala3-His4*-His5*-Glu6-Leu7-His8*-NH2 (1) reacts with two [Pd ((NH2)-N-15(CH2)(2) (NH2)-N-15)(NO3)(2)] in water to form a kinetically stable intermediate, [{Pden}(2)-{(1,4)(5,8)-peptide}](2), in which two 19-membered metallocyclic rings stabilize two peptide turns. Slow subsequent folding to a thermodynamically more stable two-turn a-helix drives the equilibrium to [{Pden}(2)-{(1,5)(4,8)-peptide}] (3), featuring two 22-membered rings. This transformation from unstructured peptide via turns to an a-helix suggests that metal clips might be useful probes for investigating peptide folding.

Protein contact prediction using patterns of correlation

We describe a new method for using neural networks to predict residue contact pairs in a protein. The main inputs to the neural network are a set of 25 measures of correlated mutation between all pairs of residues in two windows of size 5 centered on the residues of interest. While the individual pair-wise correlations are a relatively weak predictor of contact, by training the network on windows of correlation the accuracy of prediction is significantly improved. The neural network is trained on a set of 100 proteins and then tested on a disjoint set of 1033 proteins of known structure. An average predictive accuracy of 21.7% is obtained taking the best L/2 predictions for each protein, where L is the sequence length. Taking the best L/10 predictions gives an average accuracy of 30.7%. The predictor is also tested on a set of 59 proteins from the CASP5 experiment. The accuracy is found to be relatively consistent across different sequence lengths, but to vary widely according to the secondary structure. Predictive accuracy is also found to improve by using multiple sequence alignments containing many sequences to calculate the correlations. (C) 2004 Wiley-Liss, Inc.