Investigation of properties affecting controlled release of macromolecules from PLGA microspheres

Poly(lactide-co-glycolide), or PLGA, microspheres offer a widely-studied biodegradable option for controlled release of therapeutics. An array of fabrication methodologies have been developed to produce these microspheres with the capacity to encapsulate therapeutics of various types; and produce microspheres of a wide range of sizes for different methods of delivery. The encapsulation, stability, and release profiles of therapeutic release based on physical and thermodynamic properties has also been studied and modeled to an extent. Much research has been devoted to tailoring formulations for improved therapeutic encapsulation and stability as well as selective release profiles. Despite the breadth of available research on PLGA microspheres, further analysis of fundamental principles regarding the microsphere degradation, formation, and therapeutic encapsulation is necessary. This work aims to examine additional fundamental principles related to PLGA microsphere formation and degradation from solvent-evaporation of preformed polymer. In particular, mapping the development of the acidic microenvironment inside the microsphere during degradation and erosion is discussed. Also, the effect of macromolecule size and conformation is examined with respect to microsphere diameter and PLGA molecular weight. Lastly, the effects of mechanical shearing and protein exposure to aqueous media during microsphere formation are examined. In an effort to better understand the acidic microenvironment development across the microsphere diameter, pH sensitive dye conjugated to protein that undergoes conformational change at different acidic pH values was encapsulated in PLGA microspheres of diameters ranging from 40 µm to 80 µm, and used in conjunction with fluorescence resonance energy transfer to measure the radial pH change in the microspheres. Qualitative analysis of confocal micrographs was used to correlate fluorescence intensity with pH value, and obtain the radial pH across the center of the microsphere. Therapeutic encapsulation and release from polymeric microspheres is governed by an interconnected variety of factors, including the therapeutic itself. The globular protein bovine serum albumin, and the elongated and significantly smaller enzyme, lysozyme, were encapsulated in PLGA microspheres ranging from 40 µm to 80 µm in diameter. The initial surface morphology upon microsphere formation, release profiles, and microsphere erosion characteristics were explored in an effort to better understand the effect of protein size, conformation, and known PLGA interaction on the formation and degradation of PLGA microspheres and macromolecule release, with respect to PLGA molecular weight and microsphere diameter. In addition to PLGA behavior and macromolecule behavior, the effect of mechanical stresses during fabrication was examined. Two similar solvent extraction techniques were compared for the fabrication of albumin loaded microspheres. In particular, the homogeneity of the microspheres as well as capacity to retain encapsulated albumin were compared. This preliminary study paves the way for a more rigorous treatment of the effect of mechanical forces present in popular microsphere fabrication. Several factors affecting protein release from PLGA microspheres are examined herein. The technique explored for spatial resolution of the pH inside the microsphere proved mildly effective in producing a reliable method of mapping microsphere pH changes. However, notable trends with respect to microsphere size, PLGA molecular weight, and microsphere porosity were observed. Proposed methods of improving spatial resolution of the acidic microenvironment are also provided. With respect to microsphere formation, studies showed that albumin and lysozyme had little effect on the internal homogeneity of the microsphere. Rather, ionic interactions with PLGA played a more significant role in the encapsulation and release of each macromolecule. Studies also showed that higher instances of mechanical stress led to less homogeneous microspheres with lower protein encapsulation. This suggests that perhaps instead of or in addition to modifying the microsphere formation formulation, the fabrication technique itself should be more closely considered in achieving homogeneous microspheres with desired loading.

Three-dimensional, in-silico breast phantom for multimodality image simulations

We report on the construction of anatomically realistic three-dimensional in-silico breast phantoms with adjustable sizes, shapes and morphologic features. The concept of multiscale spatial resolution is implemented for generating breast tissue images from multiple modalities. Breast epidermal boundary and subcutaneous fat layer is generated by fitting an ellipsoid and 2nd degree polynomials to reconstructive surgical data and ultrasound imaging data. Intraglandular fat is simulated by randomly distributing and orienting adipose ellipsoids within a fibrous region immediately within the dermal layer. Cooper’s ligaments are simulated as fibrous ellipsoidal shells distributed within the subcutaneous fat layer. Individual ductal lobes are simulated following a random binary tree model which is generated based upon probabilistic branching conditions described by ramification matrices, as originally proposed by Bakic et al [3, 4]. The complete ductal structure of the breast is simulated from multiple lobes that extend from the base of the nipple and branch towards the chest wall. As lobe branching progresses, branches are reduced in height and radius and terminal branches are capped with spherical lobular clusters. Biophysical parameters are mapped onto the complete anatomical model and synthetic multimodal images (Mammography, Ultrasound, CT) are generated for phantoms of different adipose percentages (40%, 50%, 60%, and 70%) and are analytically compared with clinical examples. Results demonstrate that the in-silico breast phantom has applications in imaging performance evaluation and, specifically, great utility for solving image registration issues in multimodality imaging.

Quantitative characterization of physico-chemical properties of the active layers of reverse osmosis and nanofiltration membranes, and their relation to membrane performance

The main objectives of this dissertation were: (i) to develop experimental and analytical procedures to quantify different physico-chemical properties of the ultra-thin (~ 100 nm) active layers of reverse osmosis (RO) and nanofiltration (NF) membranes and their interactions with contaminants; (ii) to use such procedures to evaluate the similarities and differences between the active layers of different RO/NF membranes; and (iii) to relate characterization results to membrane performance. Such objectives were motivated by the current limited understanding of the physico-chemical properties of active layers as a result of traditional characterization techniques having limitations associated with the nanometer-scale spatial resolution required to study these ultra-thin films. Functional groups were chosen as the main active layer property of interest. Specific accomplishments of this study include the development of procedures to quantify in active layers as a function of pH: (1) the concentration of both negatively and positively ionized functional groups; (2) the stoichiometry of association between ions (i.e., barium) and ionized functional groups (i.e., carboxylate and sulfonate); and (3) the steric effects experienced by ions (i.e., barium). Conceptual and mathematical models were developed to describe experimental results. The depth heterogeneity of the active layer physico-chemical properties and interactions with contaminants studied in this dissertation was also characterized. Additionally, measured concentrations of ionized functional groups in the polyamide active layers of several commercial RO/NF membranes were used as input in a simplified RO/NF transport model to predict the rejection of a strong electrolyte (i.e., potassium iodide) and a weak acid (i.e., arsenious acid) at different pH values based on rejection results at one pH condition. The good agreement between predicted and experimental results showed that the characterization procedures developed in this study serve as useful tools in the advancement of the understanding of the properties and structure of the active layers of RO/NF membranes, and the mechanisms of contaminant transport through them.

In situ Lithium and Boron isotope determinations in mica, pyroxene, and serpentine by LA-MC-ICP-MS

Boron and Li are light, incompatible elements that preferentially partition into the liquid phase, whether melt or aqueous fluid, and thus are useful for tracking fluid-related processes in rocks. Most of the Li isotopic data presently available on subduction-related rocks are from whole-rock analyses; and the B isotopic analyses of subduction material have been carried out either on whole-rocks or in-situ on an accessory phase, such as tourmaline. The new method presented here couples an ESI New Wave UP-193-FX ArF* (193 nm) excimer laser-ablation microscope with a Neptune Plus (Thermo Scientific) MC-ICP-MS aiming to measure both Li and B isotopes in situ with good spatial resolution (metamorphic minerals are commonly chemically zoned, and whole-rock analyses lose this detail). The data thus obtained are compared with SIMS analyses on the same mineral samples for B, and with MC-ICP-MS analyses on whole-rock or mineral separates from the same sample for Li. Additionally, data acquired on tourmaline standards were compared to SIMS values. The results show that for B concentrations above 5 μg/g, the data obtained by LA-MC-ICP-MS and by SIMS are identical within error, for mica (phengitic muscovite), pyroxene (jadeite), serpentine (antigorite), and tourmaline. For Li concentrations above 10 μg/g, the data obtained by LA-MC-ICP-MS and by MC-ICP-MS are also identical, within error, for mica (phengitic muscovite), and pyroxene (jadeite). However, analyses of tourmaline standards have shown significant differences with reference values, so LA-MC-ICP-MS does not yet appear to be an appropriate method to analyze Li isotopes in tourmalines. Thus, LA-MC-ICP-MS is a suitable method to measure Li and B isotopes with good spatial resolution in major rock-forming silicates from subduction-related rocks where concentrations exceed 10 μg/g and 5 μg/g, respectively, with an error on individual measurements equal to or less than previously used methods, but obtainable in a significantly shorter amount of time. The external reproducibility is ± 2.88 to 3.31 ‰ for B and ± 1.50 to 1.75 for Li, which is lower than or equal to the variations encountered within a given chemically zoned sample (up to 10 ‰ of variation within a given natural sample).

Examination of myocardial electrophysiology using novel panoramic optical mapping techniques

Optical mapping of voltage signals has revolutionised the field and study of cardiac electrophysiology by providing the means to visualise changes in electrical activity at a high temporal and spatial resolution from the cellular to the whole heart level under both normal and disease conditions. The aim of this thesis was to develop a novel method of panoramic optical mapping using a single camera and to study myocardial electrophysiology in isolated Langendorff-perfused rabbit hearts. First, proper procedures for selection, filtering and analysis of the optical data recorded from the panoramic optical mapping system were established. This work was followed by extensive characterisation of the electrical activity across the epicardial surface of the preparation investigating time and heart dependent effects. In an initial study, features of epicardial electrophysiology were examined as the temperature of the heart was reduced below physiological values. This manoeuvre was chosen to mimic the temperatures experienced during various levels of hypothermia in vivo, a condition known to promote arrhythmias. The facility for panoramic optical mapping allowed the extent of changes in conduction timing and pattern of ventricular activation and repolarisation to be assessed. In the main experimental section, changes in epicardial electrical activity were assessed under various pacing conditions in both normal hearts and in a rabbit model of chronic MI. In these experiments, there was significant changes in the pattern of electrical activation corresponding with the changes in pacing regime. These experiments demonstrated a negative correlation between activation time and APD, which was not maintained during ventricular pacing. This suggests that activation pattern is not the sole determinant of action potential duration in intact hearts. Lastly, a realistic 3D computational model of the rabbit left ventricle was developed to simulate the passive and active mechanical properties of the heart. The aim of this model was to infer further information from the experimental optical mapping studies. In future, it would be feasible to gain insight into the electrical and mechanical performance of the heart by simulating experimental pacing conditions in the model.

Sensoriamento óptico quase-distribuído aplicado ao monitoramento de trocadores de calor de usinas hidrelétricas

This work aims to develop optical sensors for temperature monitoring in hydroelectric power plant heat exchangers. The proposed sensors are based on the Fiber Bragg Gratings technology. First of all, a prototype with three sensors inscribed in a same fiber was developed. This fiber was then fixed to a conventional Pt100 sensor rod and inserted in a thermowell. The ensemble was then calibrated in a workbench, presenting a maximum combined uncertainty of 2,06 °C. The sensor was installed in one of the heat exchangers of the Salto Osório’s hydroelectric power plant. This power plant is situated in the Iguaçu river, at the Paraná state. Despite the satisfactory results, the sensor was improved to a second version. In this, fifteen optical Bragg sensors were inscribed in a same fiber. The fixation with a conventional sensor was no longer necessary, because the first version results comproved the efficiency and response time in comparison to a conventional sensor. For this reason, it was decided to position the fiber inside a stainless steel rod, due to his low thermal expansion coefficient and high corrosion immunity. The utilization of fifteen fiber Bragg gratings aims to improve the sensor spatial resolution. Therefore, measurements every ten centimeters with respect to the heat exchanger’s height are possible. This provides the generation of a thermal map of the heat exchanger’s surface, which can be used for determination of possible points of obstruction in the hydraulic circuit of the heat exchanger. The heat exchanger’s obstruction in hydroelectric power plants usually occur by bio-fouling, and has direct influence in the generator’s cooling system efficiency. The obtained results have demonstrated the feasibility in application of the optical sensors technology in hydroelectric power plants.

Caracterização in vitro de um dosímetro de fibra ótica para braquiterapia de HDR

Recentemente foi desenvolvido um dosímetro baseado em fibras cintilantes (BCF-12 da companhia Saint Gobain Crystals com 1 e 0,5 mm de diâmetro e 5 mm de comprimento) para braquiterapia de baixa taxa de dose, em particular a braquiterapia direcionada para o tratamento do cancro da próstata. Este utiliza um novo fotomultiplicador de estado sólido dado pelo nome de MPPC - MultiPixel Photon Counter da companhia Hamamatsu Photonics (Japão). Nesta dissertação é estudado o mesmo dosímetro para a modalidade de braquiterapia de elevada taxa de dose (HDR). A informação sobre a dose neste tipo de dosímetros é obtida a partir de sinais óticos (em vez de sinais elétricos), que são imunes a interferências elétricas e eletromagnéticas. Adicionalmente as pequenas dimensões das fibras oferecem uma excelente resolução espacial e uma invasão mínima para uso em dosimetria in vivo, permitindo medir a dose diretamente ou próximo ao tumor e em tempo real. A sua utilização em braquiterapia para o cancro da próstata constitui-se assim como uma vantagem, uma vez que as fibras podem ser inseridas diretamente nos aplicadores utilizados neste tipo de tratamentos. Apesar de tudo, este tipo de dosímetros possui algumas desvantagens, como por exemplo a luz de Cherenkov e a fluorescência (forma de ruído dada pelo nome de stem effect) que, e a contrário da luz produzida pela fibra cintilante, não são diretamente proporcionais à energia depositada. Contudo, e para energias praticadas em braquiterapia de HDR, nesta dissertação, mostrou-se que este problema é pouco significativo dado que a percentagem de contribuição destes efeitos para o sinal medido é menor que 1% (ou 5% para distâncias menores que 25 mm). Ao longo desta dissertação é feita a caraterização do dosímetro (em modo corrente e impulso) e das suas várias partes em ambiente de laboratório e clínico. Nestes estudos o dosímetro, além de exibir uma boa reprodutibilidade (variação máxima de 3% entre medidas), mostrou uma alta linearidade para uma ampla gama de doses, assim como uma sensibilidade (µGy) semelhante à de uma câmara de ionização, tornando-o adequado para braquiterapia de HDR (tratamento que envolve altos gradientes de dose). Complementarmente, a sua grande versatilidade e simples utilização possibilita a sua aplicação prática em outras modalidades radioterapêuticas.

Analyse des signaux radars polarimétriques en bandes C et L pour le suivi de l'humidité du sol de sites forestiers

Résumé : Dans les couverts forestiers, le suivi de l’humidité du sol permet de prévenir plusieurs désastres tels que la paludification, les incendies et les inondations. Comme ce paramètre est très dynamique dans l’espace et dans le temps, son estimation à grande échelle présente un grand défi, d’où le recours à la télédétection radar. Le capteur radar à synthèse d’ouverture (RSO) est couramment utilisé grâce à sa vaste couverture et sa résolution spatiale élevée. Contrairement aux sols nus et aux zones agricoles, le suivi de l’humidité du sol en zone forestière est très peu étudié à cause de la complexité des processus de diffusion dans ce type de milieu. En effet, la forte atténuation de la contribution du sol par la végétation et la forte contribution de volume issue de la végétation réduisent énormément la sensibilité du signal radar à l’humidité du sol. Des études portées sur des couverts forestiers ont montré que le signal radar en bande C provient principalement de la couche supérieure et sature vite avec la densité de la végétation. Cependant, très peu d’études ont exploré le potentiel des paramètres polarimétriques, dérivés d’un capteur polarimétrique comme RADARSAT-2, pour suivre l’humidité du sol sur les couverts forestiers. L’effet du couvert végétal est moins important avec la bande L en raison de son importante profondeur de pénétration qui permet de mieux informer sur l’humidité du sol. L’objectif principal de ce projet est de suivre l’humidité du sol à partir de données radar entièrement polarimétriques en bandes C et L sur des sites forestiers. Les données utilisées sont celles de la campagne terrain Soil Moisture Active Passive Validation EXperiment 2012 (SMAPVEX12) tenue du 6 juin au 17 juillet 2012 au Manitoba (Canada). Quatre sites forestiers de feuillus ont été échantillonnés. L’espèce majoritaire présente est le peuplier faux-tremble. Les données utilisées incluent des mesures de l’humidité du sol, de la rugosité de surface du sol, des caractéristiques des sites forestiers (arbres, sous-bois, litières…) et des données radar entièrement polarimétriques aéroportées et satellitaires acquises respectivement, en bande L (UAVSAR) à 30˚ et 40˚ et en bande C (RADARSAT-2) entre 20˚ et 30˚. Plusieurs paramètres polarimétriques ont été dérivés des données UAVSAR et RADARSAT-2 : les coefficients de corrélation (ρHHVV, φHHVV, etc); la hauteur du socle; l’entropie (H), l’anisotropie (A) et l’angle alpha extraits de la décomposition de Cloude-Pottier; les puissances de diffusion de surface (Ps), de double bond (Pd) extraites de la décomposition de Freeman-Durden, etc. Des relations entre les données radar (coefficients de rétrodiffusion multifréquences et multipolarisations (linéaires et circulaires) et les paramètres polarimétriques) et l’humidité du sol ont été développées et analysées. Les résultats ont montré que 1) En bande L, plusieurs paramètres optimaux permettent le suivi de l’humidité du sol en zone forestière avec un coefficient de corrélation significatif (p-value < 0,05): σ[indice supérieur 0] linéaire et σ[indice supérieur 0] circulaire (le coefficient de corrélation, r, varie entre 0,60 et 0,96), Ps (r entre 0,59 et 0,84), Pd (r entre 0,6 et 0,82), ρHHHV_30˚, ρVVHV_30˚, φHHHV_30˚ and φHHVV_30˚ (r entre 0,56 et 0,81) alors qu’en bande C, ils sont réduits à φHHHV, φVVHV et φHHVV (r est autour de 0,90). 2) En bande L, les paramètres polarimétriques n’ont pas montré de valeur ajoutée par rapport aux signaux conventionnels multipolarisés d’amplitude, pour le suivi de l’humidité du sol sur les sites forestiers. En revanche, en bande C, certains paramètres polarimétriques ont montré de meilleures relations significatives avec l’humidité du sol que les signaux conventionnels multipolarisés d’amplitude.

DopSCAT: A mission concept for simultaneous measurements of marine winds and surface currents

A radar scatterometer operates by transmitting a pulse of microwave energy toward the ocean's surface and measuring the normalized (per-unit-surface) radar backscatter coefficient (σ°). The primary application of scatterometry is the measurement of near-surface ocean winds. By combining σ° measurements from different azimuth angles, the 10 m vector wind can be determined through a Geophysical Model Function (GMF), which relates wind and backscatter. This paper proposes a mission concept for the measurement of both oceanic winds and surface currents, which makes full use of earlier C-band radar remote sensing experience. For the determination of ocean currents, in particular, the novel idea of using two chirps of opposite slope is introduced. The fundamental processing steps required to retrieve surface currents are given together with their associated accuracies. A detailed description of the mission proposal and comparisons between real and retrieved surface currents are presented. The proposed ocean Doppler scatterometer can be used to generate global surface ocean current maps with accuracies better than 0.2 m/s at a spatial resolution better than 25 km (i.e., 12.5 km spatial sampling) on a daily basis. These maps will allow gaining some insights on the upper ocean mesoscale dynamics. The work lies at a frontier, given that the present inability to measure ocean currents from space in a consistent and synoptic manner represents one of the greatest weaknesses in ocean remote sensing.

Extrapolating Satellite Winds to Turbine Operating Heights

Ocean wind retrievals from satellite sensors are typically performed for the standard level of 10 m. This restricts their full exploitation for wind energy planning, which requires wind information at much higher levels where wind turbines operate. A new method is presented for the vertical extrapolation of satellite-based wind maps. Winds near the sea surface are obtained from satellite data and used together with an adaptation of the Monin–Obukhov similarity theory to estimate the wind speed at higher levels. The thermal stratification of the atmosphere is taken into account through a long-term stability correction that is based on numerical weather prediction (NWP) model outputs. The effect of the long-term stability correction on the wind profile is significant. The method is applied to Envisat Advanced Synthetic Aperture Radar scenes acquired over the south Baltic Sea. This leads to maps of the long-term stability correction and wind speed at a height of 100 m with a spatial resolution of 0.02°. Calculations of the corresponding wind power density and Weibull parameters are shown. Comparisons with mast observations reveal that NWP model outputs can correct successfully for long-term stability effects and also, to some extent, for the limited number of satellite samples. The satellite-based and NWP-simulated wind profiles are almost equally accurate with respect to those from the mast. However, the satellite-based maps have a higher spatial resolution, which is particularly important in nearshore areas where most offshore wind farms are built.

Areal rainfall estimation using moving cars - computer experiments including hydrological modeling

The need for high temporal and spatial resolution precipitation data for hydrological analyses has been discussed in several studies. Although rain gauges provide valuable information, a very dense rain gauge network is costly. As a result, several new ideas have been emerged to help estimating areal rainfall with higher temporal and spatial resolution. Rabiei et al. (2013) observed that moving cars, called RainCars (RCs), can potentially be a new source of data for measuring rainfall amounts. The optical sensors used in that study are designed for operating the windscreen wipers and showed promising results for rainfall measurement purposes. Their measurement accuracy has been quantified in laboratory experiments. Considering explicitly those errors, the main objective of this study is to investigate the benefit of using RCs for estimating areal rainfall. For that, computer experiments are carried out, where radar rainfall is considered as the reference and the other sources of data, i.e. RCs and rain gauges, are extracted from radar data. Comparing the quality of areal rainfall estimation by RCs with rain gauges and reference data helps to investigate the benefit of the RCs. The value of this additional source of data is not only assessed for areal rainfall estimation performance, but also for use in hydrological modeling. The results show that the RCs considering measurement errors derived from laboratory experiments provide useful additional information for areal rainfall estimation as well as for hydrological modeling. Even assuming higher uncertainties for RCs as obtained from the laboratory up to a certain level is observed practical.

Método de elaboração de mapas dos cenários climáticos futuros para o Brasil.

Human activities are altering greenhouse gas concentrations in the atmosphere and causing global climate change. The issue of impacts of human-induced climate change has become increasingly important in recent years. The objective of this work was to develop a database of climate information of the future scenarios using a Geographic Information System (GIS) tools. Future scenarios focused on the decades of the 2020?s, 2050?s, and 2080?s (scenarios A2 and B2) were obtained from the General Circulation Models (GCM) available on Data Distribution Centre from the Third Assessment Report (TAR) of Intergovernmental Panel on Climate Change (IPCC). The TAR is compounded by six GCM with different spatial resolutions (ECHAM4:2.8125×2.8125º, HadCM3: 3.75×2.5º, CGCM2: 3.75×3.75º, CSIROMk2b: 5.625×3.214º, and CCSR/NIES: 5.625×5.625º). The mean monthly of the climate variables was obtained by the average from the available models using the GIS spatial analysis tools (arithmetic operation). Maps of mean monthly variables of mean temperature, minimum temperature, maximum temperature, rainfall, relative humidity, and solar radiation were elaborated adopting the spatial resolution of 0.5° X 0.5° latitude and longitude. The method of elaborating maps using GIS tools allowed to evaluate the spatial and distribution of future climate assessments. Nowadays, this database is being used in studies of impacts of climate change on plant disease of Embrapa projects.

Microelectrical characterisation of diamond films: an attempt to understand the structural influence on electrical transport phenomena

The properties of microelectrical conduction in microwave plasma assisted chemical vapour deposition (MPCVD) diamond films were investigated using an atomic force microscopy probe, giving a morphological map of the electrical conduction with a spatial resolution better than 500 nm. Also, a cathodoluminescence map with a spatial resolution of about 1 mu m was obtained, giving the possibility of correlating the defects involved in the different carrier transport phenomena. Using micro-Raman analysis several bands could be identified. It is found that the defects responsible for the cathodoluminescence (CL) blue band are responsible for the major part of the electrical conduction in diamond films, while the defects localised in < 111 > surfaces, responsible for the green CL emission, could be involved in a less conductive process. (C) 2000 Elsevier Science S.A. All rights reserved.

Desenvolvimento e avaliação de um sistema de cálculo de dose independente para controle de qualidade de IMRT do tipo jaws-only

Dissertação (mestrado)—Universidade de Brasília, Faculdade Gama, Programa de Pós-Graduação em Engenharia Biomédica, 2015.

Whole-cell-analysis of live cardiomyocytes using wide-field interferometric phase microscopy.

We apply wide-field interferometric microscopy techniques to acquire quantitative phase profiles of ventricular cardiomyocytes in vitro during their rapid contraction with high temporal and spatial resolution. The whole-cell phase profiles are analyzed to yield valuable quantitative parameters characterizing the cell dynamics, without the need to decouple thickness from refractive index differences. Our experimental results verify that these new parameters can be used with wide field interferometric microscopy to discriminate the modulation of cardiomyocyte contraction dynamics due to temperature variation. To demonstrate the necessity of the proposed numerical analysis for cardiomyocytes, we present confocal dual-fluorescence-channel microscopy results which show that the rapid motion of the cell organelles during contraction preclude assuming a homogenous refractive index over the entire cell contents, or using multiple-exposure or scanning microscopy.