Estudiamos integrales fraccionarias de Mittag-Leffler (ML) involucradas en los procesos de solución de un sistema de ecuaciones diferenciales estocásticas fraccionarias acopladas. Introducimos el proceso estocástico fraccionario ML como una integral estocástica fraccionaria ML con respecto a un movimiento Browniano estándar. Proporcionamos algunas fórmulas de representación de procesos de solución en términos de integrales fraccionarias de Mittag-Leffler y procesos. Se dan expresiones computables de las funciones medias y de las covarianzas de dichos procesos de manera específica. Se proporciona la aplicación en modelado neuronal, y todas las funciones y procesos involucrados se determinan específicamente. Se realizan evaluaciones numéricas y se muestran y discuten algunos resultados.