Predição dos fracionamentos de proteínas de Brachiaria brizantha cv Marandu utilizando RNA

Prediction of protein fractioning of Brachiaria brizantha cv Marandu through Artificial Neural Networks

Brennecke, K.¹, Tech, A.R.B.^2*, Arce, A.I.C.^2A, Luz, P.H.C.^2B, Herling, V.R.^2C e Costa, E.J.X.^2D

¹Mestrado Stricto-Senso em Produção Animal. UNICASTELO Campus Descalvado. São Paulo. Brasil. katheryb@usp.br
²ZAB-FZEA. University of São Paulo. Pirassununga, SP. Brazil. *adriano.tech@usp.br, ^Aacaid@usp.br, ^Bphcerluz@usp.br, ^Cvrherlin@usp.br, ^Dernane@usp.br

RESUMO

Este artigo tem como objetivo relacionar variáveis morfogenéticas da forragem e dados climáticos com as frações de proteínas (A, B1, B2, B3 e C), através da rede neural aritifical conhecida como Multi-Layer Perceptron, com três camadas e algoritmo de treinamento baseado na retropropagação do gradiente do erro, a fim de criar um modelo capaz de predizer as frações protéicas de forragens a partir de suas características morfogenéticas e climáticas. Os dados utilizados para o treinamento e teste da rede foram coletados em um experimento que foi realizado em uma área experimental de 25,2 ha, com coordenadas 47^o 26' W, 21^o 59' S, formada com capim braquiarão, em um delineamento em blocos completos e casualizados, com 4 repetições e uma oferta de forragem de 5% (5 kg de massa seca por 100 kg de peso animal.dia^-1). Cada bloco foi dividido em quatro unidades experimentais de 1,575 ha, com cinco piquetes de 0,315 ha cada. As amostras foram colhidas dois dias antes da entrada dos animais, sendo a análise das frações de proteínas realizadas em laboratório para posteriormente serem comparadas com os valores estimados pela rede. Assim, através da comparação entre os dados de saída da rede e os obtidos pela análise laboratorial foi possível calcular o erro médio para as frações A, B1, B2, B3 e C de proteínas e com isso, pode-se concluir que o modelo MLP é capaz de predizer com eficiência as frações de proteínas de Brachiaria brizantha.

Palavras chave: Forragicultura.

SUMMARY

This paper aims to connect morphogenetic variables of forage and climatic data, with protein fractions (A, B1, B2, B3 and C) through the artificial neural networks known as Multi-Layer Perceptron, with three layers and algorithm training based on back-propagation of error gradient, in order to create a model to predict the protein fractions of fodder from their morphogenetic characteristics and climate. The data used for training and the test were collected in an experiment that was conducted on a 25.2 ha experimental area, located at 47^o26'W, 21^o59'S, and with pasture composed of Brachiaria brizantha cv Marandu, in a completely randomized block, with four replicates and a forage allowance of 5% (5 kg of dry matter per 100 kg animal.day^-1). Each block was divided into four experimental units of 1.575 ha, with five paddocks of 0.315 ha each. Samples were taken two days before the entry of animals, protein was analysed at laboratory to subsequently be compared with the values estimated by the network. Thus, by comparing the output of the network and those obtained by laboratory analysis, it was possible to calculate the average error for fractions A, B1, B2, B3 and C proteins and, thus, can be concluded that the model MLP is able to efficiently predict protein fractions of Brachiaria brizantha.

Key words: Forage.

Introdução

O uso de técnicas avançadas de processamento digital de sinais, principalmente o uso de modelos, baseados em redes neurais artificiais, permite inferir relações existentes entre variáveis do meio ambiente e características qualitativas de plantas. Possibilitando assim, a avaliação antecipada das características morfofisiológicas e químicas das plantas em função das condições ambientais em menor de tempo e com menor custo de análise.

As redes multi-layer percetron são um tipo de redes neurais que tem a habilidade de identificar padrões ou características após um treinamento. Essa estrutura de aprendizagem está ligada a um algoritmo de aprendizagem que é utilizado para treinar a rede e gerar informações sobre determinados eventos. Ela é formada basicamente por uma camada de neurônios de entrada, um ou mais camadas escondidas e uma camada de saída, responsável por gerar a saída da rede (Hertz et al., 1991).

Há estudos que mostram a viabilidade da RNA em programas agronômicos e zootécnicos, mostrando a capacidade das RNA's de mapear relações complexas entre variáveis ambientais e características intrínsecas das plantas e animais. Como exemplos da aplicação de RNAs na modelagem de ambientes de ambientes biológicos podem-se citar o estudo de Elizondo et al. (1994) que utilizaram uma RNA para predizer o florescimento e a maturidade fisiológica em soja, onde usaram trinta e três campos experimentais durante cinco anos; Edwards e Coob (1997) que estudaram o efeito da temperatura no acúmulo de clorofila e glicoalcalóide em batatas (Solanun tuberosum L. cv King Edward), armazenadas em 5, 10, 20 e 25^oC por 8 dias, em baixa densidade de fluxo de fóton, incluindo um modelo preliminar usando uma RNA; Parmar et al. (1997) estimaram a contaminação de aflatoxina em amendoim usando RNA, e concluíram que a RNA tem melhor performance que as técnicas tradicionais. Tambén, Yanget al. (1997a), ao utilizarem a RNA para avaliar as concentrações de pesticida no solo e Yanget al. (1997b) aplicaram RNA para simular temperatura de solos; Zee (2003), trabalhando com o programa NASA's Advance Life Support, utilizou redes neurais artificiais para caracterizar o processo de fotossíntese em soja.

Estes trabalhos mostram que a capacidade de predição e análise das RNAs em sistemas, onde as relações entre o clima x planta x animal associadas ao seu manejo, são complexas e todos os modelos comprovaram a viabilidade do emprego das RNAs.

Uma vez que a qualidade nutricional de uma forrageira está, entre outros fatores de ordem nutricionais, relacionada as frações da proteína, surge a seguinte questão: Pode se estimar essas frações de proteínas a partir de uma correlação das condições ambientais e morfofisiológicas da forrageira, a campo?

Nos ecossistemas de pastagens existe uma interação complexa entre os componentes solo-planta-animal. Portanto, o desenvolvimento, o crescimento e a senescência de folhas e perfilhos constituem processos fisiológicos que determinam mudanças significativas nos ecossistemas. Crescimento e desenvolvimento são processos distintos, porém correlacionados (Nascimento Júnior et al., 2002).

As análises desses parâmetros constituem importante ferramenta, tanto para a caracterização do potencial de produção da espécie como para a definição do potencial de uso de dado ecossistema na produção animal. Tão importante quanto saber qual será a produção é saber como esta produção está distribuída no tempo e espaço. Deste modo, a interação da planta com o meio ambiente constitui-se de fundamental importância para suportar, tanto o crescimento quanto a manutenção da capacidade produtiva da pastagem e com isso garantir um produto final de qualidade para o animal.

Em relação ao fracionamento de proteínas, o grande mérito desse sistema é caracterizar quimicamente os alimentos de forma compatível aos processos de fermentação ruminal e digestão pós-ruminal, de forma que torna possível a utilização de modelos dinâmicos a partir dos quais o desempenho animal pode ser predito, no entanto, não estabelecem diretamente as relações entre as variáveis do ambiente, as características morfológicas e fisiológicas e a qualidade nutricional, que são consideradas, em geral, complexas, necessitando assim, serem analisadas em laboratórios.

Diante disso, uma resposta a essa questão é adequar a interação à produção quantitativa e qualitativa da planta forrageira e para isso pode se usar uma ferramenta que permita uma forma de predizer as condições nutricionais da planta forrageira, ainda no campo, com a utilização de dados históricos, da planta, previamente obtidos, sem a necessidade de avaliações laboratoriais, a fim de facilitar a aquisição de dados qualitativos de pastagens e aprimoramento de manejo, além de reduzir o tempo de análises e custos.

Portanto, o propósito deste trabalho foi o de investigar a eficiência de utilização de um modelo de redes neurais do tipo MLP (Multilayer Perceptron) com o objetivo de predizer os fracionamentos de proteínas de Brachiaria brizantha cv Marandu.

Material e métodos

ÁREA DOESTUDO

O experimento foi conduzido em área pertencente à Faculdade de Zootecnia e Engenharia de Alimentos FZEA/USP, Pirassununga/SP, no período de janeiro a julho de 2004. A posição geográfica é 21^o59' de latitude sul e 47^o26' de longitude oeste, com uma altitude de 634 metros, sendo o clima subtropical do tipo Cwa (Köppen e Geiger, 1928), com inverno seco e verão quente e chuvoso.

DELINEAMENTO EXPERIMENTAL

A área formada com capim braquiarão foi dividida em 80 piquetes de 0,315 ha e posteriormente foi realizada a análise química do solo, para efeito de correção do solo.

O experimento constou de 4 ciclos de pastejo com duração de 28 dias cada, isto é, os animais ficavam durante um período total de 28 dias para completar o tempo de pastejo na área; em cada área haviam quatro unidades onde os animais ficavam durante sete dias, e o manejo dos animais estabelecido foi o método rotacionado, sendo o período de ocupação de 7 dias. Os tratamentos experimentais foram representados pela oferta de forragem a 5% (kg MS/100 kg de peso animal/dia), O delineamento experimental utilizado foi em blocos completos e casualizados com quatro repetições. Cada bloco foi dividido em quatro unidades experimentais de 1,575 ha, com cinco piquetes de 0,315 ha cada. As amostras da forragem para as análises laboratoriais foram colhidas dois dias antes da entrada dos animais, à altura do resíduo do pastejo anterior.

No período do experimento, foram realizadas análises morfofisiológicas da pastagem. As amostras de forragem colhidas foram embaladas em sacos de polietileno, identificadas e enviadas ao laboratório para pesagem da matéria original. No laboratório, essas amostras foram levadas a estufa de circulação forçada a temperatura de 55^oC, até que se atingisse uma massa constante, armazenadas para posteriores análises.

As frações protéicas foram obtidas conforme método descrito por Licitra et al. (1996) e foi realizada uma estatística básica para verificação da média dos resultados e do desvio padrão para cada ciclo de pastejo. Esses resultados serviram como valores de alimentação do banco de conhecimento, no caso os ajustes dos pesos sinápticos da rede neural, na fase de treino.

REDE NEURAL ARTIFICIAL

Os experimentos computacionais foram feitos em três etapas:

1) Preparação e digitalização dos dados obtidos pelo fracionamento de proteínas e os demais dados de ambiente;

2) Treinamento da RNA e;

3) Teste e validação da RNA.

Os dados que compuseram a etapa de preparação e digitalização foram:

A -Dados fixos: latitude, longitude, ciclos de pastejo (fevereiro a julho -1 ao 4), oferta de forragens (5, 10, 15 e 20%) e unidades experimentais (1 ao 16).

B -Dados variáveis (morfofisiológicos): média da altura pré pastejo, média da altura pós pastejo, diferenças entre as alturas pré e pós pastejo, perfilhos decapitados, perfilhos remanescentes, perfilhos novos, perfilhos aéreos remanescentes, perfilhos aéreos novos, perfilhos resprodutivos, perdas no solo, perdas na planta, perdas total, produção de massa (kg/ha) pré pastejo, produção de massa (kg/ha), pós pastejo, densidade da folha pré pastejo (<15, 15-30, 30-45, >45 cm), densidade da folha pós pastejo (<15, 15-30, 30-45, >45 cm), densidade total pré pastejo (<15, 15-30, 30-45, >45 cm), densidade total pós pastejo (<15, 1530, 30-45, >45 cm), estratificação pré pastejo folha verde (<15, 15-30, 30-45, >45 cm), estratificação pré pastejo colmo + bainha (<15, 15-30, 30-45, >45 cm), estratificação pré pastejo senescência (<15, 15-30, 30-45, >45 cm), estratificação pós pastejo folha verde (<15, 15-30, 30-45, >45 cm), estratificação pós pastejo colmo + bainha (<15, 15-30, 30-45, >45 cm), estratificação pós pastejo senescência (<15, 15-30, 30-45, >45 cm), idade no período do experimento, peso de matéria verde em função da simulação de pastejo, peso de matéria seca em função da simulação de pastejo.

C-Dados variáveis (químicos): frações de proteínas (fração A, B1, B2, B3 e C).

Para o treinamento da rede MLP foram utilizados 2/3 dos dados.

O modelo de rede utilizado por ser visualizado através da figura 1.

A figura 1 ilustra os elementos de entrada da rede, onde foram utilizados dados variáveis e dados fixos, sendo a saída da mesma representada pelas respectivas respostas, A, B1, B2, B3 e C.

Na etapa de teste e validação da RNA foi utilizado o 1/3 restante da base de dados. Foram processados pela rede já treinada todos os valores das variáveis analisadas, e foram obtidas as estimativas de fracionamento de proteínas (A, B1, B2, B3 e C). Os valores de todas as variáveis de entrada foram dispostos na rede através de médias.

A validação estatística da RNA foi efetuada segundo o método de matriz de confusão ou método Kapa (Vuolo, 1996), conforme verificado na tabela I, para cada uma das frações de proteínas. Essas matrizes de confusão foram desenvolvidas a partir da média dos resultados obtidos em laboratório (demonstrados pelos valores da coluna) com os valores absolutos obtidos na RNA (demonstrados pelos valores da linha), conforme exemplificado na figura 1. Os valores atribuídos para o teste da rede foram aleatórios e estão marcados em negrito na tabela I.

Resultados e discussões

O resultado da fração A, bem como a matriz de erro podem ser observados na tabela II. A fração A é composta de aminoácidos, peptídeos e nitrogênio não protéico, sendo de disponibilidade imediata no rúmen. Essa fração depende diretamente do teor de proteína bruta, e isso irá depender da idade da planta, bem como do manejo adotado. Como um elemento biológico, cada planta irá possuir um valor de proteína e conseqüentemente de frações de proteínas. Diante disso, uma maneira de aproximar os valores obtidos em laboratório é compondo a média das amostras e seu desvio padrão.

No caso do presente experimento, os valores da fração A de proteína obtidos em laboratório tiveram um desvio padrão médio de 0,05; 0,13; 0,13 e 0,06.

A RNA tem como resultado um único valor absoluto, porém, vindo de dados inseridos no seu treinamento. Esses dados foram coletados do sistema solo x planta ao longo do experimento.

Pode-se observar que, considerando os erros, estatisticamente, os valores médios obtidos em laboratório e o valor encontrado pela RNA, são iguais.

A tabela III mostra a média dos resultados dos valores da fração B1 de proteína no laboratório, os valores encontrados pela rede e a matriz de erro.

A subdivisão da fração B (B1, B2 e B3) é baseada nas taxas de degradação, no rumem, de cada fração.

A taxa de degradação desta fração é muito superior à taxa de pesagem da dieta no rumem, sendo totalmente degradada. Em silagens, a fração B1 corresponde à porção muito pequena da proteína solúvel (por volta de 5%), porém, em forragens verdes, a maior parte da proteína solúvel é proteína verdadeira (B1).

Pode-se observar diante dos resultados obtidos que a RNA conseguiu predizer os valores da fração B1, onde o erro foi menor que o desvio padrão encontrado no laboratório, que foi de 0,04; 0,07; 0,06 e 0,04 para os ciclos de pastejos 1, 2, 3 e 4.

Os resultados e a matriz de erro da fração B2 podem ser observados na tabela IV.

A fração B2 é calculada subtraindo-se da proteína bruta a proteína solúvel (frações A e B1), a fração B3 e a fração C são calculadas através da fórmula: B2= [PB - (A + B1 + B3 + C)].

A taxa de degradação da fração B2 é semelhante à taxa de pesagem, de forma que a quantidade efetivamente degradada no rúmen é função direta da taxa de pesagem, isto é, dependerá da relação entre taxa de degradação/taxa de passagem. Em alimentos cuja fração B2 é significativa, a correta avaliação destes parâmetros é fundamental à avaliação dos alimentos e das dietas.

O desvio padrão das frações B2 obtidas em laboratórios foram 0,04; 0,05; 0,13 e 0,06, valores aquém ao encontrado na matriz de erros, quando comparadas às médias dos valores obtidas em laboratório com o resultado da RNA, para a fração B2.

Como se pode observar a fração B2 é considerada como um resíduo e obtida por cálculo matemático, portanto, essa fração é completamente dependente dos resultados das demais frações da proteína, e a proteína, por sua vez, dependente dos processos biológicos do trinômio planta x solo x clima.

A fração B3 é determinada pela subtração entre o nitrogênio insolúvel em detergente neutro e o nitrogênio insolúvel em detergente ácido (B3= NIDN - NIDA), sendo caracterizada pela lenta degradação no rúmen, uma vez que está associada à parede celular.

Os resultados obtidos da fração B3 e a matriz de erro podem ser observados na tabela V.

Os resultados dos desvios padrão dos valores da fração B3 obtidas em laboratórios foram respectivamente para os cilos de pastejos 1, 2, 3 e 4 de 0,02; 0,01; 0,05 e 0,01.

A fração C corresponde à proteína indisponível, sendo insolúvel em detergente ácido, abreviado por NIDA (nitrogênio insolúvel em detergente ácido). Essa fração contém proteínas associadas à lignina, taninos e a compostos de Maillard altamente resistentes à degradação microbiana e enzimática, sendo considerada inaproveitável, tanto no rúmen como pós-ruminal (Sniffen et al., 1992). No entanto, sabe-se que o NIDA não é totalmente indigestível. Em forragens, os autores apontaram digestibilidade da ordem de 30%, ao passo que em concentrados esta digestibilidade atingiu até 60%.

A porcentagem resultante da fração C constata a tendência de aumento da proporção indigerível da parede celular, constituída basicamente de lignina, com o avanço da idade da planta forrageira.

Os valores da fração C da proteína e a matriz de erro podem ser observados através da tabela VI. Os valores do desvio padrão obtidos nas análises de laboratório em função dos ciclos de pastejos foram 0,01; 0,01; 0,01 e 0,01, esses valores são próximos aos erros obtidos quando comparados com o erro mostrado na matriz de erro.

Os resultados mostram que as predições realizadas pela rede MLP são semelhantes aos resultados obtidos pela análise laboratorial, ou seja, a rede MLP, após o treinamento efetuado, conseguiu determinar valores próximos aos encontrados em laboratório e com a margem de erro dentro do esperado.

Conclusões

O estudo permitiu verificar através dos testes que é possível realizar as predições referentes às frações A, B1, B2, B3 e C sem a necessidade de análises laboratoriais, diminuindo assim, o tempo de resposta e os custos operacionais das mesmas.

Observou-se que todos os teores das frações protéicas obtidas pela RNA, mesmo não estando dentro do erro esperado, estão dentro da variação encontrada em pesquisas com plantas forrageiras de clima tropical.

Conclui-se assim, que as redes neurais baseado no modelo MLP e com retropropagação de erros, consegue predizer com eficiência todas as frações de proteínas em relação às estimativas obtidas por análises laboratoriais, desde que alteradas os padrões que as relacionam nas frações que são determinadas basicamente por cálculos matemáticos excluindo elementos nutricionais.

Bibliografia

Edwards, E.J. and Cobb, A.H. 1997. Effect of temperature on glicoalkaloid and chlorophyll accumulation in potatoes (Solanum tuberosum L. cv King Edward) stored at low photon flux density, including preliminary modeling using an artificial neural network. Am. Chem. Soc., 45: 1032-1038.         [ Links ]

Elizondo, D.A., McClendon R.W. and Hoogenboom G. 1994. Neural network models for predicting flowering and phisiological maturity of soybean. Amer. Soc. Agric. Eng., 37: 981-988.         [ Links ]

Hertz, J., Krogh, A. and Palmer, R.G. 1991. Introduction to the theory of neural computation. Santa Fe Institute studies in the sciences of complexity. Computation and neural system series 1. 280 pp.         [ Links ]

Köppen, W. and Geiger, R. 1928. Klimate der erde. Wall-map 150cmx200cm. Justus Perthes, Verlag. Gotha. 1928.         [ Links ]

Krishnamoorthy, U., Sniffen, C.J. and Stern, M.D. 1983. Evaluation of a mathematical model of rumen digestion and an in vitro simulation of rumen proteolysis to estimate the rumenundegraded nitrogen content of feedstuffs. Brit. J. Nutr., 50: 555-568.         [ Links ]

Licitra, G., Hernández, T.M. and Van Soest, P.J. 1996. Nitrogen fraction in selected feedtuffs. J. Dairy Sci., 65: 217-225.         [ Links ]

Nascimento Junior, D., Neto Garces, A. F., Barbosa, R.A. e Andrade, C.M. 2002. Fundamentos sobre manejo de pastagem: Evolução e atualidade. Simpósio sobre Manejo de Pastagem. Universidade Federal de Viçosa. Viçosa. pp. 149-187.         [ Links ]

Parmar, R.S., McClendon R.W., Hoogenboom, G., Blankenship, P.D., Cole, R.J. and Dorner, J.W. 1997. Estimation of aflatoxinin contamination in preharvest peanuts using neural networks. Amer. Soc. Agric. Eng., 40: 809-813.         [ Links ]

Rumelhart, D.E., Hilton, G.E. and Willians, R.J. 1986. Learning representation by backpropagation errors. Nature, 323: 533-536.         [ Links ]

Sniffen, C.J., O'connor, J.D. and Van Soest, P.J. 1992. A net carboydrate and protein system for evaluating cattle diets: II. Carbohydrate and protein availability. J. Anim. Sci., 70: 3562-3577.         [ Links ]

Vuolo, J.H. 1996. Fundamentos da teoria dos erros. 2^a ed. Ed. Edigard Blücher Ltda. São Paulo, SP. 240 pp.         [ Links ]

Yang, C.C., Prasher, S.O., Sreekanth, S., Patni, N.K. and Masse, L. 1997a. An artificial neural network model for simulating pesticide concentrations in soil. T. ASAE, 40: 1285-1294.         [ Links ]

Yang, C.C., Prasher, S.O., Mehuys, G.R. and Patni, N.K. 1997b. Application of artificial neural networks for simulation of soil temperature. T. ASAE, 40: 649-656.         [ Links ]

Zee, F. 2003. Modeling of assimilation and transpiration in soybean plants using artificial neural networks. Sponsored by NASA. Office of life and microgravity Sciences and Applications. Jet Propulsion Laboratory. California Institute of Technology. <http://trsnew.jpl.nasa.gov/dspace/bitstream/2014/16369/1/00-2368.pdf> (30/09/2009).         [ Links ]

Recibido: 2-3-10
Aceptado: 13-4-11.