La distrofia muscular (DM) es una enfermedad en la que los enfermos presentan mutaciones que les impiden producir total (DM de Duchenne) o parcialmente (DM de Becker) la proteína distrofina. Para estudiar la biomecánica de esta enfermedad, las propiedades mecánicas de los músculos en la ausencia de distrofina, es necesario usar un animal modelo. El nemátodo Caenorhabditis elegans ha sido el elegido por investigadores de la Universidad de Pennsylvania, Philadelphia, EEUU. Mediante vídeos de alta velocidad y modelos biomecánicos han estudiado cómo afecta la distrofia muscular a la natación de estos nemátodos en un fluido a bajo número de Reynolds. Un trabajo curioso e interesante que permitirá comprender mejor esta enfermedad y cómo las terapias genéticas actúan sobre la misma. El artículo técnico ha sido aceptado para publicación en el Biophysical Journal y está disponible gratis como J. Sznitman, Prashant K. Purohit, P. Krajacic, T. Lamitina, P.E. Arratia, “Material properties of Caenorhabditis elegans swimming at low Reynolds number,” ArXiv, 9 Nov. 2009.
El movimiento ondulatorio de los nemátodos les permite nadar en un fluido newtoniano a bajo número de Reynolds en el cual dominan las fuerzas viscosas lineales dominan a las fuerzas de inercia no lineales. Estos nemátodos nadan gracias a contracciones periódicas de sus músculos dorsales y ventrales, generando ondas que se propagan desde su cabeza a su cola. Los detalles de este movimiento ya se publicaron con anterioridad y se entienden bastante bien, lo que permite estudiar el efecto de ciertas mutaciones genéticas. En especial, se pueden inferir a partir de modelos matemáticos las propiedades mecánicas de los tejidos del animal (módulo de Young, viscosidad tisular) a partir de la grabación en vídeo y análisis de las imágenes de su movimiento natatorio.
El vídeo que acompaña esta entrada muestra el análisis del movimiento (motilidad) del nemátodo Caenorhabditis elegans, de sólo 1 mm. de longitud, cuyo movimiento ha sido grabado con una cámara de alta velocidad a 125 fotogramas por segundo en un fluido acuoso caracterizado por un número de Reynolds (
Gracias a un modelo (matemático) biomecánico muy sencillo los investigadores han logrado determinar los valores del módulo de Young (E) y la viscosidad tisular (eta) de los músculos tanto de los nemátodos sanos como de los mutantes (que no producen distrofina). Para los nemátodos sanos E = 3.77 +/- 0.62 kPa, y eta = -860.2 +/- 99.4 Pa s (valor negativo porque el nemátodo genera energía en su movimiento en lugar de disiparla). Estos valores se reducen hasta en un 40% dependiendo del tipo de mutante estudiado (han estudiado tres tipos) como muestra la figura de abajo.
En resumen, un interesante estudio biomecánico de interés para profesores de física, mecánica y biomecánica, quienes podrán presentar el modelo teórico a sus alumnos, así como para profesionales de la medicina interesados en estos temas (aunque las aplicaciones biomédicas de este tipo de estudios son todavía lejanas). Para los demás, resulta curioso que la biomecánica de la natación de un gusano nos permite entender la biomecánica más íntimo de un enfermo (humano) de distrofia muscular.
En mi opinión, merece la pena leer dicho artículo, aunque tiene algunas faltas que los autores podrían haber resuelto y que los revisores parecen no haber detectado. Lo más importante, hay que reinvindicar lo español. Investigadores españoles demostraron y publicaron en Nature que no es necesario la materia oscura para explicar las curvas de rotación de las galaxias espirales. Milgrom desarrolló su teoría MOND para explicar dichas curvas. Ni MOND ni materia oscura son necesarias para explicar el grueso de dichas curvas. El campo magnético de estas galaxias las explica perfectamente. Para los que no recuerden este trabajo (publicado en una época en lo que yo estaba subscrito a Nature en papel y me sorprendía al ver el nombre de la Universidad de Granada en dicha revista) les recuerdo la cita: E. Battaner, J. L. Garrido, M. Membrado, E. Florido, “
Los modelos más sencillos son del tipo Strogatz-Sprott y se basan en cuatro estados posibles de enamoramiento que se muestran en la figura de la izquierda: (I) deseo correspondido (eager beaver), saber que la otra persona nos ama refuerza nuestro propio amor hacia ella; (II) amor precavido (cautious lover), rechazamos nuestros propios sentimientos pero los de la otra persona refuerzan nuestro amor; (III) amor ermitaño (hermit), rechazamos nuestros propios sentimientos y los de la otra persona; y (IV) tímido narcisista (narcissistic nerd), nuestro amor es intenso pero nos hecha para atrás que la otra persona también nos ame. El modelo se basa en dos ecuaciones diferenciales acopladas para las variables 
que miden el amor hacia la persona amada, correspondiendo los valores positivos a sentimientos positivos (amistad, pasión, en función de la magnitud del valor) y valores negativos a sentimientos negativos (antagonismo, desdén). El modelo propuesto es el siguiente 
representan la atracción hacia al otro. Los parámetros 
indican el grado con que un individuo ha internalizado sus propios sentimientos y su propia autoestima. Los parámetros $\beta _{i}$ representan el efecto de refuerzo que los sentimiento de la otra persona provoca en nosotros. La constante 
introduce una función de retorno que, según los autores, modela el amor entre Steve Urkel y Laura Winslow en la teleserie “Cosas de casa”: Cuando Steve se desespera, el antagonismo de Laura se reduce por su sentimiento de compasión hacia él.
y 
es la función de retorno linealizada. En cualquier otro caso, el equilibrio es inestable.
Basándose en este modelo, Alhaji Cherif y Kamal Barley introducen un nuevo modelo de carácter estocástico que presenta una mayor diversidad de comportamientos dinámicos. Este modelo corresponde a una proceso de Markov continuo cuya tabla de transición aparece a la izquierda y que conduce a una ecuación diferencial estocástica en el sentido de Ito, de la forma
![dX_{1} =\left[-\alpha _{1} X_{1} +\beta _{1} X_{2} \left(1-\varepsilon X_{2} ^{2} \right)+A_{1} \right]dt-\sqrt{\alpha _{1} X_{1} } dW_{1} +\sqrt{\beta _{1} X_{2} \left(1-\varepsilon X_{2} ^{2} \right)+A_{1} } dW_{2},](http://s0.wp.com/latex.php?latex=dX_%7B1%7D+%3D%5Cleft%5B-%5Calpha+_%7B1%7D+X_%7B1%7D+%2B%5Cbeta+_%7B1%7D+X_%7B2%7D+%5Cleft%281-%5Cvarepsilon+X_%7B2%7D+%5E%7B2%7D+%5Cright%29%2BA_%7B1%7D+%5Cright%5Ddt-%5Csqrt%7B%5Calpha+_%7B1%7D+X_%7B1%7D+%7D+dW_%7B1%7D+%2B%5Csqrt%7B%5Cbeta+_%7B1%7D+X_%7B2%7D+%5Cleft%281-%5Cvarepsilon+X_%7B2%7D+%5E%7B2%7D+%5Cright%29%2BA_%7B1%7D+%7D+dW_%7B2%7D%2C&bg=ffffff&fg=000000&s=0 "dX_{1} =\left[-\alpha _{1} X_{1} +\beta _{1} X_{2} \left(1-\varepsilon X_{2} ^{2} \right)+A_{1} \right]dt-\sqrt{\alpha _{1} X_{1} } dW_{1} +\sqrt{\beta _{1} X_{2} \left(1-\varepsilon X_{2} ^{2} \right)+A_{1} } dW_{2},")
![dX_{2} =\left[-\alpha _{2} X_{2} +\beta _{2} X_{1} \left(1-\varepsilon X_{1} ^{2} \right)+A_{2} \right]dt+\sqrt{\beta _{2} X_{1} \left(1-\varepsilon X_{1} ^{2} \right)+A_{2} } dW_{3}-\sqrt{\alpha _{2} X_{2} } dW_{4}.](http://s0.wp.com/latex.php?latex=dX_%7B2%7D+%3D%5Cleft%5B-%5Calpha+_%7B2%7D+X_%7B2%7D+%2B%5Cbeta+_%7B2%7D+X_%7B1%7D+%5Cleft%281-%5Cvarepsilon+X_%7B1%7D+%5E%7B2%7D+%5Cright%29%2BA_%7B2%7D+%5Cright%5Ddt%2B%5Csqrt%7B%5Cbeta+_%7B2%7D+X_%7B1%7D+%5Cleft%281-%5Cvarepsilon+X_%7B1%7D+%5E%7B2%7D+%5Cright%29%2BA_%7B2%7D+%7D+dW_%7B3%7D-%5Csqrt%7B%5Calpha+_%7B2%7D+X_%7B2%7D+%7D+dW_%7B4%7D.&bg=ffffff&fg=000000&s=0 "dX_{2} =\left[-\alpha _{2} X_{2} +\beta _{2} X_{1} \left(1-\varepsilon X_{1} ^{2} \right)+A_{2} \right]dt+\sqrt{\beta _{2} X_{1} \left(1-\varepsilon X_{1} ^{2} \right)+A_{2} } dW_{3}-\sqrt{\alpha _{2} X_{2} } dW_{4}.")
Los cupratos están formados por capas alternas. ¿Cuántas capas planas son necesarias para observar la superconductividad? Sólo una. Un nuevo estudio experimental publicado en Science ha observado la superconductividad con una Tc de 32 K en una película “bicapa” con una sola capa metálica, dopada con zinc, LCZO, y una sola capa aislante, LCO. El dopado con zinc de toda la película de cuprato, elimina completamente la superconductividad. Cuando sólo se dopan ciertos planos, la temperatura crítica se reduce de 32 K a solo 18 K. Logvenov y sus colegas han dopado con zinc un solo plano de una “bicapa” y han observado que la Tc se mantiene en 32 K. Interpretan su experimento como que el origen de las superconductividad se encuentra en la capa monoatómica que hace de interface entre ambas capas de la bicapa, la metálica y la aislante. Han fabricado esta estructura utilizando la técnica de epitaxia por haces moleculares (MBE). El trabajo es un gran avance experimental que no sólo aporta gran información para los teóricos sino que además tendrá múltiples aplicaciones que requieren capas superconductoras ultradelgadas. El artículo técnico es G. Logvenov, A. Gozar, I. Bozovic, “
Jean-Pierre ha publicado poco en ArXiv desde entonces, sólo algunos artículos epistemológicos sobre la 
El libro de 
La Ciencia de la Mula Francis
(Centro Nacional de Física de Partículas, Astropartículas y Nuclear, Consolider-Ingenio 2010).
