Análise de resíduos conservados e correlacionados das cisteíno proteases

As proteínas evoluem ao longo do tempo, devido ao acúmulo de mutações, inserções e deleções nos genes que as codificam. Posições fortemente conservadas geralmente referem-se a resíduos que são estritamente necessários para a estrutura da proteína ou para sua função. A metodologia DRCN (Decomposition of Residue Coevolution Networks) permite estudar coevolução de resíduos em famílias de proteínas. Para tanto, consiste em cinco etapas básicas: filtragem do alinhamento, cálculo de correlações, decomposição da rede, geração de arquivos de visualização auxiliares e anotação de posições por busca automática no UniProt. As cisteíno proteases são um grande e diverso grupo de enzimas que pertencem ao clã CA e à família C1 (família da papaína) e podem ser encontradas em todos os animais e reinos de plantas bem como em muitos vírus e organismos procariotas. O presente trabalho tem como objetivo analisar e discutir a importância dos resíduos correlacionados e conservados da família das cisteíno proteases a partir do alinhamento múltiplo de sequências. Para tanto utilizamos, o software PFSTATS, que emprega a metodologia DRCN. Foram encontrados 34 resíduos correlacionados dispostos em 5 comunidades de resíduos que coevoluem. A comunidade 1 é a maior comunidade, sendo formada por 22 resíduos que coevoluem tanto por questões estruturais, como por exemplo Cys22 e Cys63 e Cys56 e Cys101, que realizam pontes dissulfeto; quanto por questões funcionais, como Pro2, Phe28, Gly35 e Tyr89, que se localizam em sítios alostéricos preditos da cruzaína. A comunidade 2 é composta por 6 resíduos, dentre eles resíduos da díade catalítica, Cys25 e His162, e resíduos Gly23 e Gly65, que coevoluem por questões estruturais. As comunidades 3, 4 e 5 possuem apenas 2 resíduos cada. O resíduo Cys155 da comunidade 3 possui importância estrutural, bem como os resíduos Val13 e Lys181 (comunidade 4) e Gln51 (comunidade 5). Os resíduos da comunidade 5, Leu48 e Gln51 possuem importância funcional e o resíduo Ser55 (comunidade 3) ainda não apresenta sua função descrita na literatura. Também foram encontrados 10 resíduos que são altamente conservados evolutivamente nessa família. Enquanto alguns resíduos, incluindo os que compõem a díade catalítica (Cys25 e His162), são conservados por questões funcionais, outros o são por questões estruturais (Ser49, Gly168, Gly189 e Gly192). Proteins evolve over time, due to the accumulation of mutations, insertions and deletions in the genes that encode them. Highly conserved positions generally refer to residues that are strictly necessary for the structure of the protein or for its function. The DRCN (Decomposition of Residue Coevolution Networks) methodology allows the study of coevolution of residues in protein families. To do this, it consists of five basic steps: filtering the alignment, calculating correlations, decomposing the network, generating auxiliary visualization files and annotating positions by automatic search in UniProt. Cysteine proteases are a large and diverse group of enzymes belonging to the CA clan and the C1 family (papain family) and can be found in all animals and plant kingdoms as well as in many viruses and prokaryotes. The present work aims to analyze and discuss the importance of correlated and conserved residues of the cysteine protease family from the multiple sequence alignment. For this, we use the PFSTATS software, which employs the DRCN methodology. There were 34 correlated residues found in 5 coevoluted communities. Community 1 is the largest community, consisting of 22 residues coevolved by structural issues, such as Cys22 and Cys63 and Cys56 and Cys101, which carry disulfide bridges; as well as for functional issues, such as Pro2, Phe28, Gly35 and Tyr89, which are located in predicted allosteric sites of cruzain. Community 2 is composed of 6 residues, among them residues of the catalytic dyad, Cys25 and His162, and residues Gly23 and Gly65, which coevolve for structural reasons. Communities 3, 4 and 5 have only 2 residues each. The Cys155 residue of community 3 has structural importance, as well as residues Val13 and Lys181 (community 4) and Gln51 (community 5). The residues of community 5, Leu48 and Gln51 have functional importance and the residue Ser55 (community 3) still does not present its function described in the literature. Also 10 residues were found that are highly conserved evolutionarily in this family. While some residues, including those that make up the catalytic dyad (Cys25 and His162), are conserved for functional reasons, others are for structural reasons (Ser49, Gly168, Gly189 and Gly192).