Closed-loop investigation of the dynamical properties of cardiomyocyte electrophysiology by means of Dynamic clamp and Mathematical modelling

The research field of my PhD concerns mathematical modeling and numerical simulation, applied to the cardiac electrophysiology analysis at a single cell level. This is possible thanks to the development of mathematical descriptions of single cellular components, ionic channels, pumps, exchangers and subcellular compartments. Due to the difficulties of vivo experiments on human cells, most of the measurements are acquired in vitro using animal models (e.g. guinea pig, dog, rabbit). Moreover, to study the cardiac action potential and all its features, it is necessary to acquire more specific knowledge about single ionic currents that contribute to the cardiac activity. Electrophysiological models of the heart have become very accurate in recent years giving rise to extremely complicated systems of differential equations. Although describing the behavior of cardiac cells quite well, the models are computationally demanding for numerical simulations and are very difficult to analyze from a mathematical (dynamical-systems) viewpoint. Simplified mathematical models that capture the underlying dynamics to a certain extent are therefore frequently used. The results presented in this thesis have confirmed that a close integration of computational modeling and experimental recordings in real myocytes, as performed by dynamic clamp, is a useful tool in enhancing our understanding of various components of normal cardiac electrophysiology, but also arrhythmogenic mechanisms in a pathological condition, especially when fully integrated with experimental data.

γ-Aminobutyric acid type B receptor-dependent burst-firing in thalamic neurons: A dynamic clamp study

Synchronized network responses in thalamus depend on phasic inhibition originating in the thalamic reticular nucleus (nRt) and are mediated by the neurotransmitter γ-aminobutyric acid (GABA). A suggested role for intra-nRt connectivity in inhibitory phasing remains controversial. Recently, functional GABA type B (GABAB) receptors were demonstrated on nRt cells, and the slow time course of the GABAB synaptic response seems ideally suited to deinactivate low-threshold calcium channels. This promotes burst firing, a characteristic feature of synchronized responses. Here we investigate GABAB-mediated rebound burst firing in thalamic cells. Whole-cell current-clamp recordings were obtained from nRt cells and somatosensory thalamocortical relay cells in rat brain slices. Synthetic GABAB inhibitory postsynaptic potentials, generated by a hybrid computer–neuron synapse (dynamic clamp), triggered rebound low-threshold calcium spikes in both cell types when peak inhibitory postsynaptic potential hyperpolarization was greater than −92 mV. The threshold inhibitory postsynaptic potential conductance for rebound burst generation was comparable in nRt (7 nS) and thalamocortical (5 nS) cells. However, burst onset in nRt (1 s) was considerably delayed compared with thalamocortical (0.6 s) cells. Thus, GABAB inhibitory postsynaptic potentials can elicit low-threshold calcium spikes in both relay and nRt neurons, but the resultant oscillation frequency would be faster for thalamocortical–nRt networks (3 Hz) than for nRt–nRt networks (1–2 Hz). We conclude, therefore, that fast (>2 Hz) GABAB-dependent thalamic oscillations are maintained primarily by reciprocal connections between excitatory and inhibitory cells. These findings further indicate that when oscillatory neural networks contain both recurrent and reciprocal inhibition, then distinct population frequencies may result when one or the other type of inhibition is favored.

Single Synapse Information Coding in Intraburst Spike Patterns of Central Pattern Generator Motor Neurons

Burst firing is ubiquitous in nervous systems and has been intensively studied in central pattern generators (CPGs). Previous works have described subtle intraburst spike patterns (IBSPs) that, despite being traditionally neglected for their lack of relation to CPG motor function, were shown to be cell-type specific and sensitive to CPG connectivity. Here we address this matter by investigating how a bursting motor neuron expresses information about other neurons in the network. We performed experiments on the crustacean stomatogastric pyloric CPG, both in control conditions and interacting in real-time with computer model neurons. The sensitivity of postsynaptic to presynaptic IBSPs was inferred by computing their average mutual information along each neuron burst. We found that details of input patterns are nonlinearly and inhomogeneously coded through a single synapse into the fine IBSPs structure of the postsynaptic neuron following burst. In this way, motor neurons are able to use different time scales to convey two types of information simultaneously: muscle contraction (related to bursting rhythm) and the behavior of other CPG neurons (at a much shorter timescale by using IBSPs as information carriers). Moreover, the analysis revealed that the coding mechanism described takes part in a previously unsuspected information pathway from a CPG motor neuron to a nerve that projects to sensory brain areas, thus providing evidence of the general physiological role of information coding through IBSPs in the regulation of neuronal firing patterns in remote circuits by the CNS.

Estrazione automatica di parametri fisiologici da registrazioni di potenziali d'azione cardiaci

Implementazione di uno strumento in Matlab in grado di analizzare i tracciati di potenziali d'azione cardiaci e di valutarne dei potenziali d'azione in risposta a variazioni di parametri del modello cardiaco della corrente di ripolarizzazione rapida di potassio.

Real-time closed-loop electrophysiology: towards new frontiers in in vitro investigations in the neurosciences

Reflected at any level of organization of the central nervous system, most of the processes ranging from ion channels to neuronal networks occur in a closed loop, where the input to the system depends on its output. In contrast, most in vitro preparations and experimental protocols operate autonomously, and do not depend on the output of the studied system. Thanks to the progress in digital signal processing and real-time computing, it is now possible to artificially close the loop and investigate biophysical processes and mechanisms under increased realism. In this contribution, we review some of the most relevant examples of a new trend in in vitro electrophysiology, ranging from the use of dynamic-clamp to multi-electrode distributed feedback stimulation. We are convinced these represents the beginning of new frontiers for the in vitro investigation of the brain, promising to open the still existing borders between theoretical and experimental approaches while taking advantage of cutting edge technologies.

Creating a dynamic picture of the sliding clamp during T4 DNA polymerase holoenzyme assembly by using fluorescence resonance energy transfer

The coordinated assembly of the DNA polymerase (gp43), the sliding clamp (gp45), and the clamp loader (gp44/62) to form the bacteriophage T4 DNA polymerase holoenzyme is a multistep process. A partially opened toroid-shaped gp45 is loaded around DNA by gp44/62 in an ATP-dependent manner. Gp43 binds to this complex to generate the holoenzyme in which gp45 acts to topologically link gp43 to DNA, effectively increasing the processivity of DNA replication. Stopped-flow fluorescence resonance energy transfer was used to investigate the opening and closing of the gp45 ring during holoenzyme assembly. By using two site-specific mutants of gp45 along with a previously characterized gp45 mutant, we tracked changes in distances across the gp45 subunit interface through seven conformational changes associated with holoenzyme assembly. Initially, gp45 is partially open within the plane of the ring at one of the three subunit interfaces. On addition of gp44/62 and ATP, this interface of gp45 opens further in-plane through the hydrolysis of ATP. Addition of DNA and hydrolysis of ATP close gp45 in an out-of-plane conformation. The final holoenzyme is formed by the addition of gp43, which causes gp45 to close further in plane, leaving the subunit interface open slightly. This open interface of gp45 in the final holoenzyme state is proposed to interact with the C-terminal tail of gp43, providing a point of contact between gp45 and gp43. This study further defines the dynamic process of bacteriophage T4 polymerase holoenzyme assembly.

Nouveaux indices de suppression de la lipolyse par l'insuline déterminés lors de l'hyperglycémie provoquée par voie orale : comparaisons avec le clamp euglycémique-hyperinsulinémique et les paramètres métaboliques chez les femmes

Résumé : Une dysrégulation de la lipolyse des tissus adipeux peut conduire à une surexposition des tissus non-adipeux aux acides gras non-estérifiés (AGNE), qui peut mener à un certain degré de lipotoxicité dans ces tissus. La lipotoxicité constitue, par ailleurs, l’une des causes majeures du développement de la résistance à l’insuline et du diabète de type 2. En plus de ses fonctions glucorégulatrices, l’insuline a pour fonction d’inhiber la lipolyse et donc de diminuer les niveaux d’AGNE en circulation, prévenant ainsi la lipotoxicité. Il n’y a pas d’étalon d’or pour mesurer la sensibilité de la lipolyse à l’insuline. Le clamp euglycémique hyperinsulinémique constitue la méthode étalon d’or pour évaluer la sensibilité du glucose à l’insuline mais il est aussi utilisé pour mesurer la suppression de la lipolyse par l’insuline. Par contre, cette méthode est couteuse et laborieuse, et ne peut pas s’appliquer à de grandes populations. Il existe aussi des indices pour estimer la fonction antilipolytique de l’insuline dérivés de l’hyperglycémie provoquée par voie orale (HGPO), un test moins dispendieux et plus simple à effectuer à grande échelle. Cette étude vise donc à : 1) Étudier la relation entre les indices de suppressibilité des AGNE par l’insuline dérivés du clamp et ceux dérivés de l’HGPO; et 2) Déterminer laquelle de ces mesures corrèle le mieux avec les facteurs connus comme étant reliés à la dysfonction adipeuse : paramètres anthropométriques et indices de dysfonction métabolique. Les résultats montrent que dans le groupe de sujets étudiés (n=29 femmes, 15 témoins saines et 14 femmes avec résistance à l’insuline car atteintes du syndrome des ovaires polykystiques), certains indices de sensibilité à l’insuline pour la lipolyse dérivés de l’HGPO corrèlent bien avec ceux dérivés du clamp euglycémique hyperinsulinémique. Parmi ces indices, celui qui corrèle le mieux avec les indices du clamp et les paramètres anthropométriques et de dysfonction adipeuse est le T50[indice inférieur AGNE] (temps nécessaire pour diminuer de 50% le taux de base – à jeun – des AGNE). Nos résultats suggèrent donc que l’HGPO, facile à réaliser, peut être utilisée pour évaluer la sensibilité de la lipolyse à l’insuline. Nous pensons que la lipo-résistance à l’insuline peut être facilement quantifiée en clinique humaine.

Reference Values For High-density Lipoprotein Particle Size And Volume By Dynamic Light Scattering In A Brazilian Population Sample And Their Relationships With Metabolic Parameters.

Current data indicate that the size of high-density lipoprotein (HDL) may be considered an important marker for cardiovascular disease risk. We established reference values of mean HDL size and volume in an asymptomatic representative Brazilian population sample (n=590) and their associations with metabolic parameters by gender. Size and volume were determined in HDL isolated from plasma by polyethyleneglycol precipitation of apoB-containing lipoproteins and measured using the dynamic light scattering (DLS) technique. Although the gender and age distributions agreed with other studies, the mean HDL size reference value was slightly lower than in some other populations. Both HDL size and volume were influenced by gender and varied according to age. HDL size was associated with age and HDL-C (total population); non- white ethnicity and CETP inversely (females); HDL-C and PLTP mass (males). On the other hand, HDL volume was determined only by HDL-C (total population and in both genders) and by PLTP mass (males). The reference values for mean HDL size and volume using the DLS technique were established in an asymptomatic and representative Brazilian population sample, as well as their related metabolic factors. HDL-C was a major determinant of HDL size and volume, which were differently modulated in females and in males.

Interaction Between Fiscal And Monetary Policy In A Dynamic Nonlinear Model.

The objective of this study is to verify the dynamics between fiscal policy, measured by public debt, and monetary policy, measured by a reaction function of a central bank. Changes in monetary policies due to deviations from their targets always generate fiscal impacts. We examine two policy reaction functions: the first related to inflation targets and the second related to economic growth targets. We find that the condition for stable equilibrium is more restrictive in the first case than in the second. We then apply our simulation model to Brazil and United Kingdom and find that the equilibrium is unstable in the Brazilian case but stable in the UK case.

High Prevalence Of Insulin Resistance Assessed By The Glucose Clamp Technique In Hormonal And Non-hormonal Contraceptive Users.

Objective To assess the prevalence of insulin resistance (IR) and associated factors in contraceptive users. Methods A total of 47 women 18 to 40 years of age with a body mass index (kg/m(2)) < 30, fasting glucose levels < 100 mg/dl and 2-hour glucose level < 140 mg/dl after a 75-g oral glucose load were submitted to a hyperinsulinemic-euglycemic clamp. The women were distributed in tertiles regarding M-values. The analysed variables were use of combined hormonal/non-hormonal contraception, duration of use, body composition, lipid profile, glucose levels and blood pressure. Results IR was detected in 19% of the participants. The women with low M-values presented significantly higher body fat mass, waist-to-hip ratio, fasting insulin, HOMA-IR and were nulligravida, showed > 1 year of contraceptive use and higher triglyceride levels. IR was more frequent among combined oral contraceptive users, however no association was observed after regression analysis. Conclusions The prevalence of IR was high among healthy women attending a family planning clinic independent of the contraceptive method used with possible long-term negative consequences regarding their metabolic and cardiovascular health. Although an association between hormonal contraception and IR could not be found this needs further research. Family planning professionals should be proactive counselling healthy women about the importance of healthy habits.

Key Participants Of The Tumor Microenvironment Of The Prostate: An Approach Of The Structural Dynamic Of Cellular Elements And Extracellular Matrix Components During Epithelial-stromal Transition.

Cancer is a multistep process that begins with the transformation of normal epithelial cells and continues with tumor growth, stromal invasion and metastasis. The remodeling of the peritumoral environment is decisive for the onset of tumor invasiveness. This event is dependent on epithelial-stromal interactions, degradation of extracellular matrix components and reorganization of fibrillar components. Our research group has studied in a new proposed rodent model the participation of cellular and molecular components in the prostate microenvironment that contributes to cancer progression. Our group adopted the gerbil Meriones unguiculatus as an alternative experimental model for prostate cancer study. This model has presented significant responses to hormonal treatments and to development of spontaneous and induced neoplasias. The data obtained indicate reorganization of type I collagen fibers and reticular fibers, synthesis of new components such as tenascin and proteoglycans, degradation of basement membrane components and elastic fibers and increased expression of metalloproteinases. Fibroblasts that border the region, apparently participate in the stromal reaction. The roles of each of these events, as well as some signaling molecules, participants of neoplastic progression and factors that promote genetic reprogramming during epithelial-stromal transition are also discussed.

Dynamic and steady: shear rheological properties of xanthan and guar gums dispersed in yellow passion fruit pulp (Passiflora edulis f. flavicarpa)

Yellow passion fruit pulp is unstable, presenting phase separation that can be avoided by the addition of hydrocolloids. For this purpose, xanthan and guar gum [0.3, 0.7 and 1.0% (w/w)] were added to yellow passion fruit pulp and the changes in the dynamic and steady - shear rheological behavior evaluated. Xanthan dispersions showed a more pronounced pseudoplasticity and the presence of yield stress, which was not observed in the guar gum dispersions. Cross model fitting to flow curves showed that the xanthan suspensions also had higher zero shear viscosity than the guar suspensions, and, for both gums, an increase in temperature led to lower values for this parameter. The gums showed different behavior as a function of temperature in the range of 5 - 35ºC. The activation energy of the apparent viscosity was dependent on the shear rate and gum concentration for guar, whereas for xanthan these values only varied with the concentration. The mechanical spectra were well described by the generalized Maxwell model and the xanthan dispersions showed a more elastic character than the guar dispersions, with higher values for the relaxation time. Xanthan was characterized as a weak gel, while guar presented a concentrated solution behavior. The simultaneous evaluation of temperature and concentration showed a stronger influence of the polysaccharide concentration on the apparent viscosity and the G' and G" moduli than the variation in temperature.

Spectroscopic characterization and investigation of the dynamic of charge compensation process of supramolecular films derived from tetra-2-pyridyl-1,4-pyrazine ligand

This work describes the infrared spectroscopy characterization and the charge compensation dynamics in supramolecular film FeTPPZFeCN derived from tetra-2-pyridyl-1,4-pyrazine (TPPZ) with hexacyanoferrate, as well as the hybrid film formed by FeTPPZFeCN and polypyrrole (PPy). For supramolecular film, it was found that anion flux is greater in a K+ containing solution than in Li+ solution, which seems to be due to the larger crystalline ionic radius of K+. The electroneutralization process is discussed in terms of electrostatic interactions between cations and metallic centers in the hosting matrix. The nature of the charge compensation process differs from others modified electrodes based on Prussian blue films, where only cations such as K+ participate in the electroneutralization process. In the case of FeTPPZFeCN/PPy hybrid film, the magnitude of the anions’s flux is also dependent on the identity of the anion of the supporting electrolyte.

Multidisciplinary optimization of a catamaran's hydrofoil design with dynamic response surface and distributive processing

This paper presents a rational approach to the design of a catamaran's hydrofoil applied within a modern context of multidisciplinary optimization. The approach used includes the use of response surfaces represented by neural networks and a distributed programming environment that increases the optimization speed. A rational approach to the problem simplifies the complex optimization model; when combined with the distributed dynamic training used for the response surfaces, this model increases the efficiency of the process. The results achieved using this approach have justified this publication.