APROXIMACIÓN A LA ESTRUCTURA TRIDIMENSIONAL DE LA DESMINA BOVINA Y DETERMINACIÓN DEL EFECTO DE MUTACIONES PUNTUALES IN-SILICO

APPROACH TO TRIDIMENSIONAL STRUCTURE OF BOVINE DESMIN AND IN-SILICO DETERMINATION OF PUNCTUAL MUTATIONS EFFECT

J. D. Leal ^1*, L. M. Jiménez ¹

¹ Departamento de Producción Animal, Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional de Colombia, sede Bogotá. Carrera 45 nro. 26-85, Bogotá (Colombia).
^* Autor para correspondencia: jdlealg@unal.edu.co

Artículo recibido: 3 de junio de 2013 • Aprobado: 26 de septiembre de 2013

RESUMEN

La desmina es el mayor filamento proteico intermedio del miocito y desempeña un papel importante respecto de las características de calidad cárnica, dada su función intracelular de sostén y que es sustrato de los principales sistemas proteolíticos post-mortem. En la determinación de los parámetros fisicoquímicos de dicha proteína y de su ARNm, se consideraron siete mutaciones (polimorfismos de nucleótido simple, SNP) en bovino y cinco en porcino ubicadas en regiones exónicas. Mediante procedimientos computacionales se obtuvo un modelo tridimensional que incluyó desde el aminoácido 39 al 470 de la secuencia DAA32384.1. Se identificó que las mutaciones T49C y A45C del ARNm del bovino son responsables de una modificación en la estructura bidimensional del ARNm y de la disminución de su estabilidad in-silico, por lo que se les considera como las mutaciones más significativas para evaluar experimentalmente en bovinos.

Palabras clave: estructura celular, proteína, mutación, ARNm.

---

ABSTRACT

Desmin is the major proteic intermediate filament protein of muscular cell and it has an important effect on meat quality features because this protein is structural and during post-mortem conversion is substrate of proteolysis systems. Seven functional SNP (Single Nucleotide Polymorphisms) of bovine and five of porcine were evaluated in several physicochemical features of the molecules RNAm and protein of Desmin. A 3D model was obtained from aminoacid 39 to 470 of the sequence DAA32384.1. Mutations T49C and A45C of bovine RNAm modify the molecular structure and these have less in-silico stability. On this way, these should be tested to experimental level.

Key words: cell structure, protein, mutation, RNAm.

INTRODUCCIÓN

La desmina (DES) es el mayor filamento proteico intermedio de la célula muscular y su gen se considera altamente conservado desde el tiburón hasta el humano (Bar et al. 2004; Goldfarb et al. 2004). Se encuentra localizado principalmente en la periferia del disco-Z de los sarcómeros (Paulin y Li 2004) y los conecta con la membrana plasmática y la lámina nuclear en las células musculares estriadas, cardíacas y lisas (Goldfarb et al. 2004). Esta proteína juega un papel importante en el mantenimiento de la integridad de la estructura celular (Paulin y Li 2004) por ser uno de los filamentos intermedios principales involucrados en la conservación del alineamiento lateral de las miofibrillas adyacentes al unir la capa periférica de éstas con la membrana celular (Zhang et al. 2006).

Huff-Lonergan y Lonergan (2005), Ouali et al. (2006) y Zhang et al. (2006) determinaron que la degradación de proteínas estructurales, como la DES, durante el perído post-mortem temprano tiene varias consecuencias a nivel estructural que modifican atributos del músculo permitiendo su conversión a carne (Campo et al. 2000; Monsón et al. 2005, Monsón et al. 2004; Sañudo et al. 2004; Kolczak et al. 2005).Varias investigaciones han podido establecer que la desmina tiene un rol importante en varios parámetros de calidad cárnica, por ser uno de los sustratos del sistema proteolítico de la calpaína-calpastatina (Huff-Lonergan y Lonergan, 2005; Zhang et al. 2006).

Se han identificado varios polimorfismos en genes de este sistema proteolítico en bovinos y porcinos que se han asociado con variaciones en la calidad cárnica respecto de parámetros como la terneza (Curi et al. 2010; Iwanowska et al. 2011; Motter et al. 2009; Pinto et al. 2010) y la Capacidad de Retención de Agua (CRA) (Zhang et al. 2006; Bond y Warner 2007). En el estudio de Chang et al. (2003), que empleó un microsatélite ubicado en la región 3´UTR del gen de la desmina del cerdo, se logró determinar su asociación con el valor de pH a 45 minutos y 24 horas post-mortem, dada la co-segregación de este marcador neutro con el polimorfismo responsable. Sin embargo, es posible que la asociación pueda ser producida por una mutación en el gen PRKAG3, gen que yace muy próximo a la desmina, tanto en porcinos como en bovinos.

Es posible que las mutaciones en el gen de la desmina puedan contribuir a alteraciones fenotípicas en parámetros cárnicos asociados, como la terneza y la CRA en especies como el bovino y el cerdo, por lo que se le considera como uno de los genes-candidato asociados a parámetros de calidad cárnica (Williams et al. 2009). A esto contribuye la alta variabilidad reportada en otros mamíferos y en el humano (Goldfarb et al. 2004), en el cual se han detectado múltiples variaciones que van desde cambios nucleotídicos, hasta indels (del inglés insertion-deletion, o inserción/deleción) y splicings alternativos (del inglés alternative splicing, o empalme alternativo), lo que incrementa su valor en estudios de asociación fenotipo-genotipo.

El objetivo de este estudio fue determinar la estructura 3D de la proteína desmina en bovinos y el conjunto de SNP (del inglés Single Nucleotide Polymorphism, o polimorfismo en un solo nucleótido) con una mayor probabilidad de producir alteraciones funcionales en el ARNm y en la proteína mediante análisis computacional.

MATERIALES Y MÉTODOS

Se emplearon las secuencias de ARNm NM_001081575.1 y NM_001001535.1 (bovino y porcino, respectivamente) y para las proteínas DAA32384.1 y AAD46492.1 (bovina y porcina, respectivamente), las cuales serán denominadas en este artículo como 'alelo silvestre'. Las mutaciones puntuales de las regiones exónicas evaluadas se listan en la Tabla 1 y fueron ubicadas manualmente. También se emplearon las secuencias NM_001927.3 y NP_001918.3 correspondientes al ARNm y la proteína de la desmina humana, así como un total de 11 mutaciones missense que se ha reportado se encuentran asociadas con miopatías en humanos (Van Spaendonck et al. 2011).

Construcción de modelos de ARNm

Para el modelamiento del ARNm de la DES bovina fueron seleccionados siete SNPs (secuencias 1 a 7 de la Tabla 1) de los 25 reportados en la base de datos de NCBI (National Center for Biotechnology Information), debido a que se producen en un nucleótido de la región exónica y que, por consiguiente, alteran la secuencia del ARNm. Para el modelamiento del ARNm del gen de DES porcino se utilizaron las cinco mutaciones reportadas en la base de datos de NCBI que afectan regiones exónicas traducidas (secuencias 8 a 12 de la Tabla 1). También se realizó el modelamiento del ARNm de la DES humana normal el cual se designó como 'humano silvestre'.

Para determinar la estructura secundaria adoptada más probablemente por las secuencias del ARNm y establecer el valor de Energía Mínima Libre (ΔG) se empleó el software MFold (Zuker, 2003) a temperatura fisiológica.

Construcción de modelos de la proteína DES

Para el modelamiento de la proteína DES bovina fueron seleccionadas dos mutaciones (secuencias 3 y 4 de la Tabla 1) debido a que son missense. Para realizar una comparación con el modelo obtenido de la proteína DES bovina se ubicaron las mutaciones missense reportadas que alteran la proteína DES porcina (secuencias 10 y 12 de la Tabla 1) y las 11 mutaciones missense reportadas para la proteína DES humana (Van Spaendonck et al. 2011).

La estructura secundaria de las secuencias proteicas se determinó mediante el servidor Quick2d (toolkit.tuebingen.mpg.de/quick2_d) y las regiones desordenadas con Disopred (Ward et al. 2004). Adicionalmente se detectaron los dominios constituyentes de la proteína mediante el servidor CD-Search (Marchler y Bryant 2004) del NCBI.

Se utilizó el modelamiento por homología para la determinación de la estructura terciaria del Multidominio Tipo Filamento mediante el servidor Swiss-Model  y el software Deep Viewer (Guex y Peitsch 1997; www.expasy.org/spdbv/), mientras que para el modelamiento de la Superfamilia Cabeza de Filamento (aminoácidos 38 al 110) y los segmentos 240-271 y 407-470, se empleó el servidor I-Tasser (Zhang 2008; Roy et al. 2010), efectuando la minimización energética de cada segmento mediante la herramienta Gromos96 del software Deep Viewer (Guex y Peitsch 1997) para su posterior ensamblaje. El punto isoeléctrico (PI) y peso molecular se determinaron mediante el servidor Protparam (www.expasy.org/spdbv/).

RESULTADOS

En la Figura 1 se presenta la estructura y tamaño del gen de la DES bovina. Adicionalmente se presentan, los dominios predichos en la proteína lineal. De igual forma se ilustran las regiones modeladas con el servidor Swiss-Model . Para el proceso de modelamiento se emplearon los templetes 1gk7A con una identidad de 78,95% (entre los aminoácidos 106 y 143), 3s4rB con una identidad del 69,23% (entre los aminoácidos 139 y 194), 3swkA con una identidad del 68,61% (entre los aminoácidos 158 y 243), 3kltC con 76,06% de identidad (entre los aminoácidos 268 y 338) y 1gk4B con 77,22% de identidad modelando los aminoácidos 333 al 411.

Modelos de ARNm

En la Figura 2 se presenta la estructura secundaria adoptada por el ARNm silvestre del gen de la DES bovina (bovino silvestre) y las tres modificaciones estructurales que producían los mayores cambios visibles en el plegamiento del ARNm: en orden descendente se encuentran las secuencias 2, 1 y 5, lo que las convierte en buenas candidatas para contribuir a la variación fenotípica en el bovino, especialmente en cuanto a parámetros de calidad cárnica, puesto que generan polimorfismos conformacionales en el ARNm. Por su parte, en la Figura 3 se muestran las estructuras del ARNm del gen DES porcino silvestre y el modelo de la secuencia 8 debido a que fue la que tuvo el mayor impacto sobre el modelo bidimensional y, adicionalmente, la estructura del ARNm del gen DES silvestre del humano.

A nivel de estabilidad molecular, el ARNm del alelo silvestre de la DES bovina mostró un valor de Energía Mínima Libre (ΔG) de -919,81 kcal/mol. En orden descendente, las mutaciones de la DES bovina que desestabilizan en mayor medida a esta molécula son los SNP 2, 1 y 3, con valores de ΔG de -911,07, -911,83 y -914,41 kcal/mol respectivamente. Se destaca el SNP 4 cuyo valor de ΔG -926,91 kcal/mol le puede conferir un carácter estabilizante mayor a esta molécula. El alelo silvestre del cerdo presentó un valor de ΔG de -891,36, mientras que sus secuencias 8 y 11 mostraron un valor de -890,26 y -897,90 kcal/mol, respectivamente.

Modelo proteico

En la Figura 4 se presenta la estructura secundaria predicha de la proteína DES bovina. En la Tabla 2 se muestran los parámetros de punto isoeléctrico (PI) y peso molecular para las secuencias seleccionadas. Se debe tener en cuenta que, para obtener la predicción final del modelo, el segmento anterior al aminoácido 39 fue eliminado debido a que la predicción realizada por CD-Search no lo establecía como constituyente de la Superfamilia Cabeza de Filamento. Los ocho segmentos modelados de la proteína DES bovina (cinco resueltos mediante Swiss-Model  y tres con I-Tasser) tenían fragmentos solapados (≈5 aminoácidos) por lo que pudieron ser ensamblados con el programa Deep Viewer.

En la Figura 5a se presenta el modelo final de la proteína desmina bovina comprendida entre los aminoácidos 39 y 470 de la secuencia problema DAA32384.1. Se señala la presencia de los segmentos constituyentes de la secuencia Rod (PCD, 1A, L1, 1B, L12, 2A, L2 y 2B); dentro de este último se ubica el denominado 'Stutter' que corresponde a 10 aminoácidos (350 a 359) altamente conservados en las proteínas de los filamentos intermedios celulares (Bar et al. 2004).

En el modelo obtenido de la proteína DES bovina se indican las posiciones 38 y 50 correspondientes a las secuencias missense 3 y 4 (Figura 5b). Sobre este modelo se señalan también las posiciones homólogas a los aminoácidos 224 y 296 del cerdo de las secuencias 10 y 12 respectivamente (Figura 5c) y las posiciones 213, 337, 342, 345, 357, 360, 370, 385, 389, 406 y 449 correspondientes a las 11 mutaciones missense reportadas para la proteína DES humana (Figuras 5d y 5e).

DISCUSIÓN

El gen de la desmina se considera como uno de los genes candidato que se asocia con fenotipos relativos a parámetros de calidad cárnica. Por tanto, se puede asumir que la existencia de un bajo número de mutaciones reportadas en bases de datos (25 SNP en bovino y 5 en porcino), posiblemente se deriva de una reducida cantidad de estudios. lo que permite sugerir la necesidad de implementar análisis que contribuyan a la caracterización de su estructura y funcionamiento fisiológico.

Variaciones en el ARNm

Las configuraciones análogas adoptadas por la molécula de ARNm (Figura 2), ocasionadas por las mutaciones T49C, A45C, T177A y A213C (secuencias 2, 1, 5 y 4), pueden tener un efecto sobre la regulación en la expresión de la proteína (Andrade 2001). Adicionalmente a las variaciones en la estructura tridimensional adoptada por esta molécula, varias de las mutaciones mencionadas modifican su valor de Energía Mínima Libre, haciéndola in-silico más o menos estables. Las mutaciones T49C y A45C modifican, tanto la estructura conformacional como la Energía Mínima Libre, lo que hace menos estables a las moléculas de ARNm de la DES bovina.

Estudios previos utilizando estos mismos parámetros (configuración bidimensional y ΔG) han permitido confirmar que mutaciones sinónimas o mutaciones ubicadas en las regiones no traducidas del ARNm de genes candidato son causantes de variaciones fenotípicas en humanos (Duan et al. 2003; Johnson et al. 2011; Shen et al. 1999). Por su parte, Halvorsen et al. (2010), examinando SNP asociados con enfermedades del genoma humano que alteran significativamente la conformación de las regiones reguladoras UTR, definieron los RiboSNitches que inducen cambios conformacionales permanentes de la estructura del ARNm, actuando como interruptores moleculares (molecular switch) que apagan o encienden la expresión del gen. Ello se puede contrastar con los polimorfismos T49C y A45C, los cuales están ubicados en la región 5´UTR del ARNm de la DES y, por lo tanto, pueden tener un efecto fenotípico sobre parámetros de calidad cárnica, como la CRA y la terneza, entre otras.

Efecto de mutaciones missense en la proteína DES

De las mutaciones missense evaluadas en las tres especies analizadas, únicamente se presentó modificación del PI en tres de las 11 de la DES humana (N342D, K449T y R406W) (Van Spaendonck et al. 2011) siendo la R406W la que presenta el cuadro clínico más severo, que incluso puede conducir a la muerte (Park et al. 2000). El PI se modifica desde 5,21 para la DES humana silvestre hasta 5,17 en los tres casos citados anteriormente. En los casos de la DES de bovino y cerdo, el PI no se presentó variación alguna que permita postular a estas mutaciones como candidatas para una evaluación in-vivo.

Es interesante resaltar la ubicación de 10 de las 11 mutaciones humanas con efecto fenotípico en la región comprendida entre la región 2B del Rod (Figura 5e) y la parte terminal de la proteína, lo que podría señalar a este segmento como prometedor para llevar a cabo estudios de asociación fenotipo-genotipo en especies como el bovino y el porcino, con el fin de establecer su posible efecto sobre parámetros de calidad cárnica, dadas las funciones de sostén estructural, generación y transmisión de fuerza activa y pasiva desempeñadas por la desmina (Paulin y Li 2004).

Modelo tridimensional de la DES bovina

En el modelo de estructura bidimensional de la DES bovina se predijo la existencia de una α-hélice continua desde el aminoácido 99 al 417; sin embargo Bar et al. (2004) y Van Spaendonck et al. (2011) mencionan en esta misma ubicación la existencia de una estructura de α-hélice interrumpida por tres enlaces o linkers (L1, L12 y L2), lo que permite la formación de cuatro segmentos separados en α-hélice denominados 1A, 1B, 2A y 2B (Figura 5a).

Cabe destacar que en humanos (Park et al. 2000; Goldfarb et al. 2004) se han encontrado mutaciones de novo en 45% de los casos clínicos asociados con fenotipos especiales, como en el caso de la presentación de cardiomiopatías causadas por la agregación proteica con un gran contenido de desmina y menores proporciones de alfa-β-cristalina, distrofina y miotilina, generando la inhabilidad de la desmina para interactuar con otras estructuras celulares (Van Spaendonck et al. 2011). Lo anterior evidencia una susceptibilidad del gen DES a sufrir mutaciones principalmente en el segmento 2B (Park et al. 2000; Van Spaendonck et al. 2011) estableciendo al exón 6 de este gen como un Hot-spot (Goldfarb et al. 2004).

El bajo número de mutaciones reportadas en el gen DES de bovinos y porcinos podría indicar la necesidad de realizar estudios que contribuyan al conocimiento de la estructura y funcionamiento de este gen, ya que se postula la existencia de una gran variabilidad con posibles efectos fenotípicos. Es posible que en especies de animales domésticos, como el bovino, no se encuentren mutaciones con un efecto tan evidente como en el humano, dados los mecanismos de selección artificial a los que vienen siendo sometidos; sin embargo, no se puede descartar la existencia de mutaciones que contribuyan a la variabilidad fenotípica (Goldfarb et al. 2004) y que permitan relacionar la estructura y funcionalidad del miocito con parámetros de calidad cárnica como la CRA y la terneza. Con relación a lo anterior, Goldfarb et al. (2004) destacan que los filamentos de desmina que poseen un plegamiento diferente al normal, presentan una mayor resistencia al remodelamiento in vitro producido por los sistemas enzimáticos responsables. Cabe resaltar que uno de esos sistemas de remodelamiento in vivo a nivel celular, el sistema calpaína-calpastatina, es el responsable del proceso de tenderización post-mortem en el paso de músculo a carne, por lo que una mutación en el gen de la desmina podría modificar este proceso, alterando su proteólisis post-mortem, lo que afectaría por consiguiente parámetros de calidad de la carne como la CRA y la terneza de la carne (Campo et al. 2000; Kolczak et al. 2005; Monsón et al. 2005, Monsón et al. 2004; Sañudo et al. 2004; Zhang et al. 2006).

CONCLUSIONES

Se logró modelar la proteína DES del bovino en su totalidad y se pudo confirmar la presencia de los segmentos característicos de los filamentos intermedios del miocito de los mamíferos, que habían sido reportados previamente. La desmina es uno de los genes candidato con un posible efecto sobre parámetros de calidad cárnica en bovinos y en porcinos; sin embargo, en la actualidad sólo existe un número reducido de polimorfismos identificados en ambas especies, en comparación con los reportados para la DES humana, proteína en la que se ha podido determinar una amplia variabilidad de la región 2B del Multidominio Tipo Filamento.

A través de la evaluación del polimorfismo conformacional del ARNm de la DES del bovino se logró determinar que las mutaciones T49C y A45C pueden ser reconocidas como las más importantes, ya que al incrementar el valor de Energía Mínima Libre de la molécula, la hacen menos estable.

REFERENCIAS

Citation:

Leal JD, Jiménez LM. 2013. Aproximación a la estructura tridimensional de la desmina bovina y determinación del efecto de mutaciones puntuales in-silico [Approach to tridimensional structure of bovine desmin and in-silico determination of punctual mutations effect]. Rev Fac Med Vet Zoot. 62(3): 157-168.