Use of genetic algorithms and solvation potential to study peptides structure

In this work, genetic algorithms concepts along with a rotamer library for proteins side chains and implicit solvation potential are used to optimize the tertiary structure of peptides. We starting from the known PDB structure of its backbone which is kept fixed while the side chains allowed adopting the conformations present in the rotamer library. It was used rotamer library independent of backbone and a implicit solvation potential. The structure of Mastoporan-X was predicted using several force fields with a growing complexity; we started it with a field where the only present interaction was Lennard-Jones. We added the Coulombian term and we considered the solvation effects through a term proportional to the solvent accessible area. This paper present good and interesting results obtained using the potential with solvation term and rotamer library. Hence, the algorithm (called YODA) presented here can be a good tool to the prediction problem. (c) 2007 Elsevier B.V. All rights reserved.

Predicting peptides structure with solvation potential and rotamer library dependent of the backbone

In this work, genetic algorithms concepts along with a rotamer library dependent of backbone and implicit solvation potential are used to study the tertiary structure of peptides. We starting from known primary sequence and optimize the structure of the backbone while the side chains allowed adopting the conformations present in a rotamer library. The GA, implemented with two force fields with a growing complexity, was used predict the structure of a polyalanine and a poly-isolueucine. This paper presents good and interesting results about the study of peptides structures and about the development of computational tools to study peptides structures. (C) 2007 Elsevier B.V. All rights reserved.

Ab initio protein folding simulations using atomic burials as informational intermediates between sequence and structure

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Bothrops moojeni myotoxin-II, a Lys49-phospholipase A(2) homologue: An example of function versatility of snake venom proteins

MjTX-II, a myotoxic phospholipase A(2) (PLA(2)) homologue from Bothrops moojeni venom, was functionally and structurally characterized. The MjTX-II characterization included: (i) functional characterization (antitumoral, antimicrobial and antiparasitic effects); (ii) effects of structural modifications by 4-bromophenacyl bromide (BPB), cyanogen bromide (CNBr), acetic anhydride and 2-nitrobenzenesulphonyl fluoride (NBSF); (iii) enzymatic characterization: inhibition by low molecular weight heparin and EDTA; and (iv) molecular characterization: cDNA sequence and molecular structure prediction. The results demonstrated that MjTX-II displayed antimicrobial activity by growth inhibition against Escherichia coli and Candida albicans, antitumoral activity against Erlich ascitic tumor (EAT), human breast adenocarcinoma (SK-BR-3) and human T leukemia cells (JURKAT) and antiparasitic effects against Schistosoma mansoni and Leishmania spp., which makes MjTX-II a promising molecular model for future therapeutic applications, as well as other multifunctional homologous Lys49-PLA(2)S or even derived peptides. This work provides useful insights into the structural determinants of the action of Lys49-PLA2 homologues and, together with additional strategies, supports the concept of the presence of others bioactive sites distinct from the catalytic site in snake venom myotoxic PLA(2)s. (c) 2005 Elsevier B.V. All rights reserved.

Using genetic algorithm to design protein sequence

In this work, genetic algorithms concepts along with a rotamer library for proteins side chains are used to optimize the tertiary structure of the hydrophobic core of Cytochrome b(562) starting from the known PDB structure of its backbone which is kept fixed while the side chains of the hydrophobic core are allowed to adopt the conformations present in the rotamer library. The atoms of the side chains forming the core interact via van der Waals energy. Besides the prediction of the native core structure, it is also suggested a set of different amino acid sequences for this core. Comparison between these new cores and the native are made in terms of their volumes, van der Waals energies values and the numbers of contacts made by the side chains forming the cores. This paper proves that genetic algorithms area efficient to design new sequence for the protein core. (C) 2007 Elsevier B.V. All rights reserved.

O uso de rede neural artificial MLP na predição de estruturas secundárias de proteínas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Alternativas para redução de efeitos de multicolinearidade em modelos de avaliação de efeitos genéticos em bovinos de corte

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Cupin: A candidate molecular structure for the Nep1-like protein family

Background: NEP1-like proteins (NLPs) are a novel family of microbial elicitors of plant necrosis. Some NLPs induce a hypersensitive-like response in dicot plants though the basis for this response remains unclear. In addition, the spatial structure and the role of these highly conserved proteins are not known.Results: We predict a 3d-structure for the beta-rich section of the NLPs based on alignments, prediction tools and molecular dynamics. We calculated a consensus sequence from 42 NLPs proteins, predicted its secondary structure and obtained a high quality alignment of this structure and conserved residues with the two Cupin superfamily motifs. The conserved sequence GHRHDWE and several common residues, especially some conserved histidines, in NLPs match closely the two cupin motifs. Besides other common residues shared by dicot Auxin-Binding Proteins (ABPs) and NLPs, an additional conserved histidine found in all dicot ABPs was also found in all NLPs at the same position.Conclusion: We propose that the necrosis inducing protein class belongs to the Cupin superfamily. Based on the 3d-structure, we are proposing some possible functions for the NLPs.

Competition and phylogeny determine community structure in Mullerian co-mimics

Until recently, the study of negative and antagonistic interactions (for example, competition and predation) has dominated our understanding of community structure, maintenance and assembly(1). Nevertheless, a recent theoretical model suggests that positive interactions (for example, mutualisms) may counterbalance competition, facilitating long-term coexistence even among ecologically undifferentiated species(2). Mullerian mimics are mutualists that share the costs of predator education(3) and are therefore ideally suited for the investigation of positive and negative interactions in community dynamics. The sole empirical test of this model in a Mullerian mimetic community supports the prediction that positive interactions outweigh the negative effects of spatial overlap(4) (without quantifying resource acquisition). Understanding the role of trophic niche partitioning in facilitating the evolution and stability of Mullerian mimetic communities is now of critical importance, but has yet to be formally investigated. Here we show that resource partitioning and phylogeny determine community structure and outweigh the positive effects of Mullerian mimicry in a species-rich group of neotropical catfishes. From multiple, independent reproductively isolated allopatric communities displaying convergently evolved colour patterns, 92% consist of species that do not compete for resources. Significant differences in phylogenetically conserved traits (snout morphology and body size) were consistently linked to trait-specific resource acquisition. Thus, we report the first evidence, to our knowledge, that competition for trophic resources and phylogeny are pivotal factors in the stable evolution of Mullerian mimicry rings. More generally, our work demonstrates that competition for resources is likely to have a dominant role in the structuring of communities that are simultaneously subject to the effects of both positive and negative interactions.

Relationship between global structural parameters and Enzyme Commission hierarchy: Implications for function prediction

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.

Uncertainties in the prediction of spatial variability of soil CO2 emissions and related properties

A emissão de CO2 do solo apresenta alta variabilidade espacial, devido à grande dependência espacial observada nas propriedades do solo que a influenciam. Neste estudo, objetivou-se: caracterizar e relacionar a variabilidade espacial da respiração do solo e propriedades relacionadas; avaliar a acurácia dos resultados fornecidos pelo método da krigagem ordinária e simulação sequencial gaussiana; e avaliar a incerteza na predição da variabilidade espacial da emissão de CO2 do solo e demais propriedades utilizando a simulação sequencial gaussiana. O estudo foi conduzido em uma malha amostral irregular com 141 pontos, instalada sobre a cultura de cana-de-açúcar. Nesses pontos foram avaliados a emissão de CO2 do solo, a temperatura do solo, a porosidade livre de água, o teor de matéria orgânica e a densidade do solo. Todas as variáveis apresentaram estrutura de dependência espacial. A emissão de CO2 do solo mostrou correlações positivas com a matéria orgânica (r = 0,25, p < 0,05) e a porosidade livre de água (r = 0,27, p <0,01) e negativa com a densidade do solo (r = -0,41, p < 0,01). No entanto, quando os valores estimados espacialmente (N=8833) são considerados, a porosidade livre de água passa a ser a principal variável responsável pelas características espaciais da respiração do solo, apresentando correlação de 0,26 (p < 0,01). As simulações individuais propiciaram, para todas as variáveis analisadas, melhor reprodução das funções de distribuição acumuladas e dos variogramas, em comparação à krigagem e estimativa E-type. As maiores incertezas na predição da emissão de CO2 estiveram associadas às regiões da área estudada com maiores valores observados e estimados, produzindo estimativas, ao longo do período estudado, de 0,18 a 1,85 t CO2 ha-1, dependendo dos diferentes cenários simulados. O conhecimento das incertezas gerado por meio dos diferentes cenários de estimativa pode ser incluído em inventários de gases do efeito estufa, resultando em estimativas mais conservadoras do potencial de emissão desses gases.

Prediction of cassava starch edible film properties by chemometric analysis of infrared spectra

Cassava starch has been shown to make transparent and colorless flexible films without any previous chemical treatment. The functional properties of edible films are influenced by starch properties, including chain conformation, molecular bonding, crystallinity, and water content. Fourier-transform infrared (FTIR) spectroscopy in combination with attenuated total reflectance (ATR) has been applied for the elucidation of the structure and conformation of carbohydrates. This technique associated with chemometric data processing could indicate the relationship between the structural parameters and the functional properties of cassava starch-based edible films. Successful prediction of the functional properties values of the starch-based films was achieved by partial least squares regression data. The results showed that presence of the hydroxyl group on carbon 6 of the cyclic part of glucose is directly correlated with the functional properties of cassava starch films.

A structure-activity relationship study of quinone compounds with trypanocidal activity

A set of 25 quinone compounds with anti-trypanocidal activity was studied by using the density functional theory (DFT) method in order to calculate atomic and molecular properties to be correlated with the biological activity. The chemometric methods principal component analysis (PCA), hierarchical cluster analysis (HCA), stepwise discriminant analysis (SDA), Kth nearest neighbor (KNN) and soft independent modeling of class analogy (SIMCA) were used to obtain possible relationships between the calculated descriptors and the biological activity studied and to predict the anti-trypanocidal activity of new quinone compounds from a prediction set. Four descriptors were responsible for the separation between the active and inactive compounds: T-5 (torsion angle), QTS1 (sum of absolute values of the atomic charges), VOLS2 (volume of the substituent at region B) and HOMO-1 (energy of the molecular orbital below HOMO). These descriptors give information on the kind of interaction that occurs between the compounds and the biological receptor. The prediction study was done with a set of three new compounds by using the PCA, HCA, SDA, KNN and SIMCA methods and two of them were predicted as active against the Trypanosoma cruzi. (c) 2005 Elsevier SAS. All rights reserved.

Prediction and spatial variability of soil dynamic properties in sugar cane fields of São Paulo State Brazil

The objective of this work was to model and diagnose the spatial variability of soil load support capacity (SLSC) in sugar cane crop fields, as well as to evaluate the management impact on São Paulo State soil structure. The investigated variables were: pressure preconsolidation (sigma(p)), apparent cohesion () and internal friction angle (). The conclusions from the results were that the models and spatial dependence maps constitute important tools in the prediction and location of the mechanical internal strength of soils cultivated with sugar cane. They will help future soil management decisions so that soil structure sustainability will not be compromised.