RESUMEN
Stroke is one of the main causes of mortality and long-term disability worldwide. The pathophysiological mechanisms underlying this disease are not well understood, particularly in the chronic phase after the initial ischemic episode. In this study, a Macaca fascicularis stroke model consisting of two sample groups, as determined by MRI-quantified infarct volumes as a measure of the stroke severity 28 days after the ischemic episode, was evaluated using qualitative and quantitative proteomics analyses. By using multiple online multidimensional liquid chromatography platforms, 8790 nonredundant proteins were identified that condensed to 5223 protein groups at 1% global false discovery rate (FDR). After the application of a conservative criterion (5% local FDR), 4906 protein groups were identified from the analysis of cerebral cortex. Of the 2068 quantified proteins, differential proteomic analyses revealed that 31 and 23 were dysregulated in the elevated- and low-infarct-volume groups, respectively. Neurogenesis, synaptogenesis, and inflammation featured prominently as the cellular processes associated with these dysregulated proteins. Protein interaction network analysis revealed that the dysregulated proteins for inflammation and neurogenesis were highly connected, suggesting potential cross-talk between these processes in modulating the cytoskeletal structure and dynamics in the chronic phase poststroke. Elucidating the long-term consequences of brain tissue injuries from a cellular prospective, as well as the molecular mechanisms that are involved, would provide a basis for the development of new potentially neurorestorative therapies.
Asunto(s)
Corteza Cerebral/química , Regulación de la Expresión Génica , Proteómica/métodos , Accidente Cerebrovascular/metabolismo , Animales , Enfermedad Crónica , Modelos Animales de Enfermedad , Inflamación/genética , Macaca fascicularis , Imagen por Resonancia Magnética , Neurogénesis/genética , Mapas de Interacción de ProteínasRESUMEN
Metastatic renal cell carcinoma (RCC) is one of the most treatment-resistant malignancies, and patients have a dismal prognosis, with a <10% five-year survival rate. The identification of markers that can predict the potential for metastases will have a great effect in improving patient outcomes. In this study, we used differential proteomics with isobaric tags for relative and absolute quantitation (iTRAQ) labeling and LC-MS/MS analysis to identify proteins that are differentially expressed in metastatic and primary RCC. We identified 1256 non-redundant proteins, and 456 of these were quantified. Further analysis identified 29 proteins that were differentially expressed (12 overexpressed and 17 underexpressed) in metastatic and primary RCC. Dysregulated protein expressions of profilin-1 (Pfn1), 14-3-3 zeta/delta (14-3-3ζ), and galectin-1 (Gal-1) were verified on two independent sets of tissues by means of Western blot and immunohistochemical analysis. Hierarchical clustering analysis showed that the protein expression profile specific for metastatic RCC can distinguish between aggressive and non-aggressive RCC. Pathway analysis showed that dysregulated proteins are involved in cellular processes related to tumor progression and metastasis. Furthermore, preliminary analysis using a small set of tumors showed that increased expression of Pfn1 is associated with poor outcome and is a potential prognostic marker in RCC. In addition, 14-3-3ζ and Gal-1 also showed higher expression in tumors with poor prognosis than in those with good prognosis. Dysregulated proteins in metastatic RCC represent potential prognostic markers for kidney cancer patients, and a greater understanding of their involved biological pathways can serve as the foundation of the development of novel targeted therapies for metastatic RCC.
Asunto(s)
Carcinoma de Células Renales/genética , Neoplasias Renales/genética , Proteínas de Neoplasias/análisis , Proteoma/análisis , Proteínas 14-3-3/metabolismo , Biomarcadores de Tumor/análisis , Carcinoma de Células Renales/metabolismo , Carcinoma de Células Renales/mortalidad , Carcinoma de Células Renales/secundario , Cromatografía Liquida , Progresión de la Enfermedad , Galectina 1/metabolismo , Expresión Génica , Perfilación de la Expresión Génica , Regulación Neoplásica de la Expresión Génica , Humanos , Neoplasias Renales/metabolismo , Neoplasias Renales/mortalidad , Metástasis de la Neoplasia , Profilinas/metabolismo , Pronóstico , Proteómica , Espectrometría de Masas en TándemRESUMEN
BACKGROUND: A major barrier to effective treatment of glioblastoma multiforme (GBM) is the invasion of glioma cells into the brain parenchyma rendering local therapies such as surgery and radiation therapy ineffective. GBM patients with such highly invasive and infiltrative tumors have poor prognosis with a median survival time of only about a year. However, the mechanisms leading to increased cell migration, invasion and diffused behavior of glioma cells are still poorly understood. METHODS: In the current study, we applied quantitative proteomics for the identification of differentially expressed proteins in GBMs as compared to non-malignant brain tissues. RESULTS: Our study led to the identification of 23 proteins showing overexpression in GBM; these include membrane proteins, moesin and CD44. The results were verified using Western blotting and immunohistochemistry in independent set of GBM and non-malignant brain tissues. Both GBM tissues and glioma cell lines (U87 / U373) demonstrated membranous expression of moesin and CD44, as revealed by immunohistochemistry and immunofluorescence, respectively. Notably, glioma cells transfected with moesin siRNA displayed reduced migration and invasion on treatment with hyaluronan (HA), an important component of the extracellular matrix in GBM. CD44, a transmembrane glycoprotein, acts as a major receptor for hyaluronan (HA). Using co-immunoprecipitation assays, we further demonstrated that moesin interacts with CD44 in glioma cells only after treatment with HA; this implicates a novel role of moesin in HA-CD44 signaling in gliomas. CONCLUSIONS: Our results suggest that development of inhibitors which interfere with CD44-moesin interactions may open a new avenue in the future to mitigate cellular migration in gliomas.
Asunto(s)
Movimiento Celular/efectos de los fármacos , Glioblastoma/metabolismo , Ácido Hialurónico/farmacología , Proteínas de Microfilamentos/metabolismo , Ciclo Celular/efectos de los fármacos , Línea Celular Tumoral , Proliferación Celular/efectos de los fármacos , Regulación Neoplásica de la Expresión Génica/efectos de los fármacos , Glioblastoma/genética , Humanos , Receptores de Hialuranos/genética , Receptores de Hialuranos/metabolismo , Proteínas de Microfilamentos/genética , Unión Proteica/efectos de los fármacos , Proteoma , ProteómicaRESUMEN
Myogenesis is a well-characterized program of cellular differentiation that is exquisitely sensitive to the extracellular milieu. Systematic characterization of the myogenic secretome (i.e. the ensemble of secreted proteins) is, therefore, warranted for the identification of novel secretome components that regulate both the pluripotency of these progenitor mesenchymal cells, and also their commitment and passage through the differentiation program. Previously, we have successfully identified 26 secreted proteins in the mouse skeletal muscle cell line C2C12 (1). In an effort to attain a more comprehensive picture of the regulation of myogenesis by its extracellular milieu, quantitative profiling employing stable isotope labeling by amino acids in cell culture was implemented in conjunction with two parallel high throughput online reverse phase liquid chromatography-tandem mass spectrometry systems. In summary, 34 secreted proteins were quantified, 30 of which were shown to be differentially expressed during muscle development. Intriguingly, our analysis has revealed several novel up- and down-regulated secretome components that may have critical biological relevance for both the maintenance of pluripotency and the passage of cells through the differentiation program. In particular, the altered regulation of secretome components, including follistatin-like protein-1, osteoglycin, spondin-2, and cytokine-induced apoptosis inhibitor-1, along with constitutively expressed factors, such as fibulin-2, illustrate dynamic changes in the secretome that take place when differentiation to a specific lineage occurs.
Asunto(s)
Desarrollo de Músculos , Fibras Musculares Esqueléticas/metabolismo , Músculo Esquelético/fisiología , Mioblastos Esqueléticos/metabolismo , Proteoma/metabolismo , Secuencia de Aminoácidos , Animales , Isótopos de Carbono , Técnicas de Cultivo de Célula , Diferenciación Celular , Medios de Cultivo Condicionados/análisis , Regulación del Desarrollo de la Expresión Génica , Genes Reporteros , Marcaje Isotópico , Luciferasas/biosíntesis , Luciferasas/genética , Ratones , Datos de Secuencia Molecular , Fibras Musculares Esqueléticas/citología , Músculo Esquelético/citología , Músculo Esquelético/metabolismo , Mioblastos Esqueléticos/citología , Fragmentos de Péptidos/química , Regiones Promotoras Genéticas , Proteoma/química , Espectrometría de Masas en TándemRESUMEN
As the number of older adults increases, so does the pressure on health care systems due to age-related disorders. Attempts to reduce cognitive decline have focused on individual interventions such as exercise or diet, with limited success. This study adopted a different approach by investigating the impact of combined daily activities on memory decline. We used data from the National Institute of Aging's Health and Retirement Study to explore two new questions: does combining activities affect memory decline, and if yes, does this impact change across the lifespan? We created a new machine learning model using 33 daily activities and involving 3210 participants. Our results showed that the effect of combined activities on memory decline was stronger than any individual activity's impact. Moreover, this effect increased with age, whereas the importance of historical factors such as education, and baseline memory decreased. The present findings point out the importance of selecting multiple, diverse activities for older adults as they age. These results could have a significant impact on aging health policies promoting new programs such as social prescribing.
Asunto(s)
Cognición/fisiología , Disfunción Cognitiva/prevención & control , Envejecimiento Saludable/fisiología , Actividades Recreativas/psicología , Actividades Cotidianas , Anciano , Anciano de 80 o más Años , Relaciones Familiares/psicología , Femenino , Humanos , Aprendizaje Automático , MasculinoRESUMEN
Formalin-fixed paraffin-embedded (FFPE) tissues are the primary and preferred medium for archiving patients' samples. Here we demonstrate relative quantifications of protein biomarkers in extracts of laser microdissected epithelial cells from FFPE endometrial carcinoma tissues versus those from normal proliferative endometria by means of targeted proteomic analyses using LC-multiple reaction monitoring (MRM) MS with MRM Tags for Relative and Absolute Quantitation (mTRAQ) labeling. Comparable results of differential expressions for pyruvate kinase isoform M2 (PK-M2) and polymeric Ig receptor were observed between analyses on laser microdissected epithelial cells from FFPE tissues and corresponding homogenates from frozen tissues of the same individuals that had previously been analyzed and reported. We also identified PK-M2 in the normal proliferative phase of the endometrium. Other biomarkers in addition to PK-M2 and polymeric Ig receptor were also observed but not consistently and/or were at levels below the threshold for quantification.
Asunto(s)
Biomarcadores de Tumor/análisis , Neoplasias Endometriales/química , Marcaje Isotópico/métodos , Espectrometría de Masas/métodos , Secuencia de Aminoácidos , Biomarcadores de Tumor/metabolismo , Cromatografía por Intercambio Iónico , Neoplasias Endometriales/metabolismo , Neoplasias Endometriales/patología , Células Epiteliales/metabolismo , Femenino , Fase Folicular , Formaldehído , Humanos , Microdisección , Adhesión en Parafina , Fragmentos de Péptidos/análisis , Fragmentos de Péptidos/metabolismo , Piruvato Quinasa/análisis , Piruvato Quinasa/metabolismo , Receptores de Inmunoglobulina Polimérica/análisis , Receptores de Inmunoglobulina Polimérica/metabolismoRESUMEN
Using the notion of complexity and synchrony, this study presents a data-driven pipeline of nonlinear analysis of neuromagnetic sources reconstructed from human magnetoencephalographic (MEG) data collected in reaction to vibrostimulation of the right index finger. The dynamics of MEG source activity was reconstructed with synthetic aperture magnetometry (SAM) beam-forming technique. Considering brain as a complex system, we applied complexity-based tools to identify brain areas with dynamic patterns that remain regular across repeated stimulus presentations, and to characterize their synchronized behavior. Volumetric maps of brain activation were calculated using sample entropy as a measure of signal complexity. The complexity analysis identified activity in the primary somatosensory (SI) area contralateral to stimuli and bilaterally in the posterior parietal cortex (PPC) as regions with decreased complexity, consistently expressed in a group of subjects. Seeding an activated source with low complexity in the SI area, cross-sample entropy was used to generate synchrony maps. Cross-sample entropy analysis confirmed the synchronized dynamics of neuromagnetic activity between areas SI and PPC, robustly expressed across subjects. Our results extend the understanding of synchronization between co-activated brain regions, focusing on temporal coordination between events in terms of synchronized multidimensional signal patterns.
Asunto(s)
Mapeo Encefálico/métodos , Encéfalo/fisiología , Dedos/fisiología , Magnetoencefalografía/métodos , Procesamiento de Señales Asistido por Computador , Percepción del Tacto/fisiología , Algoritmos , Simulación por Computador , Lateralidad Funcional , Humanos , Imagen por Resonancia Magnética , Dinámicas no Lineales , Lóbulo Parietal/fisiología , Estimulación Física , Corteza Somatosensorial/fisiología , VibraciónRESUMEN
Large survey databases for aging-related analysis are often examined to discover key factors that affect a dependent variable of interest. Typically, this analysis is performed with methods assuming linear dependencies between variables. Such assumptions however do not hold in many cases, wherein data are linked by way of non-linear dependencies. This in turn requires applications of analytic methods, which are more accurate in identifying potentially non-linear dependencies. Here, we objectively compared the feature selection performance of several frequently-used linear selection methods and three non-linear selection methods in the context of large survey data. These methods were assessed using both synthetic and real-world datasets, wherein relationships between the features and dependent variables were known in advance. In contrast to linear methods, we found that the non-linear methods offered better overall feature selection performance than linear methods in all usage conditions. Moreover, the performance of the non-linear methods was more stable, being unaffected by the inclusion or exclusion of variables from the datasets. These properties make non-linear feature selection methods a potentially preferable tool for both hypothesis-driven and exploratory analyses for aging-related datasets.
Asunto(s)
Algoritmos , Bases de Datos Factuales , Procesamiento Automatizado de Datos , Modelos Teóricos , HumanosRESUMEN
There are no serum biomarkers for the accurate diagnosis of clear cell renal cell carcinoma (ccRCC). Diagnosis and decision of nephrectomy rely on imaging which is not always accurate. Non-invasive diagnostic biomarkers are urgently required. In this study, we preformed quantitative proteomics analysis on a total of 199 patients including 30 matched pairs of normal kidney and ccRCC using isobaric tags for relative and absolute quantitation (iTRAQ) labeling and LC-MS/MS analysis to identify differentially expressed proteins. We found 55 proteins significantly dysregulated in ccRCC compared to normal kidney tissue. 54 were previously reported to play a role in carcinogenesis, and 39 are secreted proteins. Dysregulation of alpha-enolase (ENO1), L-lactate dehydrogenase A chain (LDHA), heat shock protein beta-1 (HSPB1/Hsp27), and 10 kDa heat shock protein, mitochondrial (HSPE1) was confirmed in two independent sets of patients by western blot and immunohistochemistry. Pathway analysis, validated by PCR, showed glucose metabolism is altered in ccRCC compared to normal kidney tissue. In addition, we examined the utility of Hsp27 as biomarker in serum and urine. In ccRCC patients, Hsp27 was elevated in the urine and serum and high serum Hsp27 was associated with high grade (Grade 3-4) tumors. These data together identify potential diagnostic biomarkers for ccRCC and shed new light on the molecular mechanisms that are dysregulated and contribute to the pathogenesis of ccRCC. Hsp27 is a promising diagnostic marker for ccRCC although further large-scale studies are required. Also, molecular profiling may help pave the road to the discovery of new therapies.
Asunto(s)
Biomarcadores de Tumor/metabolismo , Carcinoma de Células Renales/metabolismo , Neoplasias Renales/metabolismo , Biomarcadores de Tumor/sangre , Biomarcadores de Tumor/genética , Biomarcadores de Tumor/orina , Carcinoma de Células Renales/sangre , Carcinoma de Células Renales/genética , Carcinoma de Células Renales/orina , Estudios de Casos y Controles , Femenino , Perfilación de la Expresión Génica , Proteínas de Choque Térmico HSP27/sangre , Proteínas de Choque Térmico HSP27/genética , Proteínas de Choque Térmico HSP27/orina , Proteínas de Choque Térmico , Humanos , Inmunohistoquímica , Neoplasias Renales/sangre , Neoplasias Renales/genética , Neoplasias Renales/orina , Masculino , Chaperonas Moleculares , Proteínas de Neoplasias/sangre , Proteínas de Neoplasias/genética , Proteínas de Neoplasias/metabolismo , Proteínas de Neoplasias/orina , Pronóstico , Proteómica/métodosRESUMEN
Linear and non-linear techniques for inferring causal relations between the brain signals representing the underlying neuronal systems have become a powerful tool to extract the connectivity patterns in the brain. Typically these tools employ the idea of Granger causality, which is ultimately based on the temporal precedence between the signals. At the same time, phase synchronization between coupled neural ensembles is considered a mechanism implemented in the brain to integrate relevant neuronal ensembles to perform a cognitive or perceptual task. Phase synchronization can be studied by analyzing the effects of phase-locking between the brain signals. However, we should expect that there is no one-to-one mapping between the observed phase lag and the time precedence as specified by physically interacting systems. Specifically, phase lag observed between two signals may interfere with inferring causal relations. This could be of critical importance for the coupled non-linear oscillating systems, with possible time delays in coupling, when classical linear cross-spectrum strategies for solving phase ambiguity are not efficient. To demonstrate this, we used a prototypical model of coupled non-linear systems, and compared three typical pipelines of inferring Granger causality, as established in the literature. Specifically, we compared the performance of the spectral and information-theoretic Granger pipelines as well as standard Granger causality in their relations to the observed phase differences for frequencies at which the signals become synchronized to each other. We found that an information-theoretic approach, which takes into account different time lags between the past of one signal and the future of another signal, was the most robust to phase effects.
Asunto(s)
Encéfalo/fisiología , Modelos Neurológicos , Animales , Causalidad , Haplorrinos , Neurofisiología , Dinámicas no LinealesRESUMEN
Dynamics of brain signals such as electroencephalogram (EEG) can be characterized as a sequence of quasi-stable patterns. Such patterns in the brain signals can be associated with coordinated neural oscillations, which can be modeled by non-linear systems. Further, these patterns can be quantified through dynamical non-stationarity based on detection of qualitative changes in the state of the systems underlying the observed brain signals. This study explored age-related changes in dynamical non-stationarity of the brain signals recorded at rest, longitudinally with 128-channel EEG during early adolescence (10 to 13 years of age, 56 participants). Dynamical non-stationarity was analyzed based on segmentation of the time series with subsequent grouping of the segments into clusters with similar dynamics. Age-related changes in dynamical non-stationarity were described in terms of the number of stationary states and the duration of the stationary segments. We found that the EEG signal became more non-stationary with age. Specifically, the number of states increased whereas the mean duration of the stationary segment decreased with age. These two effects had global and parieto-occipital distribution, respectively, with the later effect being most dominant in the alpha (around 10 Hz) frequency band.
Asunto(s)
Encéfalo/fisiología , Electroencefalografía , Adolescente , Factores de Edad , Ritmo alfa , Niño , Femenino , Humanos , Masculino , Estimulación LuminosaRESUMEN
Variability in source dynamics across the sources in an activated network may be indicative of how the information is processed within a network. Information-theoretic tools allow one not only to characterize local brain dynamics but also to describe interactions between distributed brain activity. This study follows such a framework and explores the relations between signal variability and asymmetry in mutual interdependencies in a data-driven pipeline of non-linear analysis of neuromagnetic sources reconstructed from human magnetoencephalographic (MEG) data collected as a reaction to a face recognition task. Asymmetry in non-linear interdependencies in the network was analyzed using transfer entropy, which quantifies predictive information transfer between the sources. Variability of the source activity was estimated using multi-scale entropy, quantifying the rate of which information is generated. The empirical results are supported by an analysis of synthetic data based on the dynamics of coupled systems with time delay in coupling. We found that the amount of information transferred from one source to another was correlated with the difference in variability between the dynamics of these two sources, with the directionality of net information transfer depending on the time scale at which the sample entropy was computed. The results based on synthetic data suggest that both time delay and strength of coupling can contribute to the relations between variability of brain signals and information transfer between them. Our findings support the previous attempts to characterize functional organization of the activated brain, based on a combination of non-linear dynamics and temporal features of brain connectivity, such as time delay.
RESUMEN
BACKGROUND: The number of patients with endometrial carcinoma (EmCa) with advanced stage or high histological grade is increasing and prognosis has not improved for over the last decade. There is an urgent need for the discovery of novel molecular targets for diagnosis, prognosis and treatment of EmCa, which will have the potential to improve the clinical strategy and outcome of this disease. METHODOLOGY AND RESULTS: We used a "drill-down" proteomics approach to facilitate the identification of novel molecular targets for diagnosis, prognosis and/or therapeutic intervention for EmCa. Based on peptide ions identified and their retention times in the first LC-MS/MS analysis, an exclusion list was generated for subsequent iterations. A total of 1529 proteins have been identified below the Proteinpilot® 5% error threshold from the seven sets of iTRAQ experiments performed. On average, the second iteration added 78% new peptides to those identified after the first run, while the third iteration added 36% additional peptides. Of the 1529 proteins identified, only 40 satisfied our criteria for significant differential expression in EmCa in comparison to normal proliferative tissues. These proteins included metabolic enzymes (pyruvate kinase M2 and lactate dehydrogenase A); calcium binding proteins (S100A6, calcyphosine and calumenin), and proteins involved in regulating inflammation, proliferation and invasion (annexin A1, interleukin enhancer-binding factor 3, alpha-1-antitrypsin, macrophage capping protein and cathepsin B). Network analyses revealed regulation of these molecular targets by c-myc, Her2/neu and TNF alpha, suggesting intervention with these pathways may be a promising strategy for the development of novel molecular targeted therapies for EmCa. CONCLUSIONS: Our analyses revealed the significance of drill-down proteomics approach in combination with iTRAQ to overcome some of the limitations of current proteomics strategies. This study led to the identification of a number of novel molecular targets having therapeutic potential for targeted molecular therapies for endometrial carcinoma.
Asunto(s)
Neoplasias Endometriales/química , Terapia Molecular Dirigida/métodos , Proteínas de Neoplasias/análisis , Proteómica/métodos , Cromatografía Liquida , Neoplasias Endometriales/diagnóstico , Neoplasias Endometriales/tratamiento farmacológico , Femenino , Humanos , Pronóstico , Espectrometría de Masas en TándemRESUMEN
The extent and effects of sequence scrambling in peptide ions during tandem mass spectrometry (MS/MS) have been examined using tryptic peptides from model proteins. Sequence-scrambled b ions appeared in about 35% of 43 tryptic peptides examined under MS/MS conditions. In general, these ions had relatively low abundances with averages of 8% and 16%, depending on the instrumentation used. A few tryptic peptides gave abundant scrambled b ions in MS/MS. However, peptide and protein identifications under proteomic conditions with Mascot were not affected, even for these peptides wherein scrambling was prominent. From the 43 tryptic peptides that have been investigated, the conclusion is that sequence scrambling is unlikely to impact negatively on the accuracy of automated peptide and protein identifications in proteomics.
Asunto(s)
Compuestos Macrocíclicos/química , Fragmentos de Péptidos/química , Proteómica/métodos , Análisis de Secuencia de Proteína/métodos , Espectrometría de Masas en Tándem/métodos , Secuencia de Aminoácidos , Animales , Bovinos , Caballos , Modelos Moleculares , Datos de Secuencia Molecular , Fragmentos de Péptidos/metabolismo , Proteínas/química , Proteínas/metabolismo , Proteómica/normas , Reproducibilidad de los Resultados , Análisis de Secuencia de Proteína/normas , Tripsina/metabolismoRESUMEN
Addressing the issue of effective connectivity, this study focuses on effects of indirect connections on inferring stable causal relations: partial transfer entropy. We introduce a Granger causality measure based on a multivariate version of transfer entropy. The statistic takes into account the influence of the rest of the network (environment) on observed coupling between two given nodes. This formalism allows us to quantify, for a specific pathway, the total amount of indirect coupling mediated by the environment. We show that partial transfer entropy is a more sensitive technique to identify robust causal relations than its bivariate equivalent. In addition, we demonstrate the confounding effects of the variation in indirect coupling on the detectability of robust causal links. Finally, we consider the problem of model misspecification and its effect on the robustness of the observed connectivity patterns, showing that misspecifying the model may be an issue even for model-free information-theoretic approach.
Asunto(s)
Causalidad , Entropía , Modelos Estadísticos , Algoritmos , Encéfalo/fisiología , Simulación por Computador , Electroencefalografía , Ambiente , Humanos , Memoria a Corto Plazo/fisiología , Pruebas Neuropsicológicas , Procesamiento de Señales Asistido por ComputadorRESUMEN
The blood oxygen level-dependent (BOLD) fMRI signal does not measure neuronal activity directly. This fact is a key concern for interpreting functional imaging data based on BOLD. Mathematical models describing the path from neural activity to the BOLD response allow us to numerically solve the inverse problem of estimating the timing and amplitude of the neuronal activity underlying the BOLD signal. In fact, these models can be viewed as an advanced substitute for the impulse response function. In this work, the issue of estimating the dynamics of neuronal activity from the observed BOLD signal is considered within the framework of optimization problems. The model is based on the extended "balloon" model and describes the conversion of neuronal signals into the BOLD response through the transitional dynamics of the blood flow-inducing signal, cerebral blood flow, cerebral blood volume and deoxyhemoglobin concentration. Global optimization techniques are applied to find a control input (the neuronal activity and/or the biophysical parameters in the model) that causes the system to follow an admissible solution to minimize discrepancy between model and experimental data. As an alternative to a local linearization (LL) filtering scheme, the optimization method escapes the linearization of the transition system and provides a possibility to search for the global optimum, avoiding spurious local minima. We have found that the dynamics of the neural signals and the physiological variables as well as the biophysical parameters can be robustly reconstructed from the BOLD responses. Furthermore, it is shown that spiking off/on dynamics of the neural activity is the natural mathematical solution of the model. Incorporating, in addition, the expansion of the neural input by smooth basis functions, representing a low-pass filtering, allows us to model local field potential (LFP) solutions instead of spiking solutions.