Cell-specific models of hiPSC-CMs developed by the gradient-based parameter optimization method fitting two different action potential waveforms.

Parameter optimization (PO) methods to determine the ionic current composition of experimental cardiac action potential (AP) waveform have been developed using a computer model of cardiac membrane excitation. However, it was suggested that fitting a single AP record in the PO method was not always successful in providing a unique answer because of a shortage of information. We found that the PO method worked perfectly if the PO method was applied to a pair of a control AP and a model output AP in which a single ionic current out of six current species, such as IKr, ICaL, INa, IKs, IKur or IbNSC was partially blocked in silico. When the target was replaced by a pair of experimental control and IKr-blocked records of APs generated spontaneously in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), the simultaneous fitting of the two waveforms by the PO method was hampered to some extent by the irregular slow fluctuations in the Vm recording and/or sporadic alteration in AP configurations in the hiPSC-CMs. This technical problem was largely removed by selecting stable segments of the records for the PO method. Moreover, the PO method was made fail-proof by running iteratively in identifying the optimized parameter set to reconstruct both the control and the IKr-blocked AP waveforms. In the lead potential analysis, the quantitative ionic mechanisms deduced from the optimized parameter set were totally consistent with the qualitative view of ionic mechanisms of AP so far described in physiological literature.

Mathematical analysis of left ventricular elastance with respect to afterload change during ejection phase.

Since the left ventricle (LV) has pressure (Plv) and volume (Vlv), we can define LV elastance from the ratio between Plv and Vlv, termed as "instantaneous elastance." On the other hand, end-systolic elastance (Emax) is known to be a good index of LV contractility, which is measured by the slope of several end-systolic Plv-Vlv points obtained by using different loads. The word Emax originates from the assumption that LV elastance increases during the ejection phase and attains its maximum at the end-systole. From this concept, we can define another elastance determined by the slope of isochronous Plv-Vlv points, that is Plv-Vlv points at a certain time after the ejection onset time by using different loads. We refer to this elastance as "load-dependent elastance." To reveal the relation between these two elastances, we used a hemodynamic model that included a detailed ventricular myocyte contraction model. From the simulation results, we found that the isochronous Plv-Vlv points lay in one line and that the line slope corresponding to the load-dependent elastance slightly decreased during the ejection phase, which is quite different from the instantaneous elastance. Subsequently, we analyzed the mechanism determining these elastances from the model equations. We found that instantaneous elastance is directly related to contraction force generated by the ventricular myocyte, but the load-dependent elastance is determined by two factors: one is the transient characteristics of the cardiac cell, i.e., the velocity-dependent force drops characteristics in instantaneous shortening. The other is the force-velocity relation of the cardiac cell. We also found that the linear isochronous pressure-volume relation is based on the approximately linear relation between the time derivative of the cellular contraction force and the cellular shortening velocity that results from the combined characteristics of LV and aortic compliances.

Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model.

Although repolarization has been suggested to propagate in cardiac tissue both theoretically and experimentally, it has been challenging to estimate how and to what extent the propagation of repolarization contributes to relaxation because repolarization only occurs in the course of membrane excitation in normal hearts. We established a mathematical model of a 1D strand of 600 myocytes stabilized at an equilibrium potential near the plateau potential level by introducing a sustained component of the late sodium current (INaL). By applying a hyperpolarizing stimulus to a small part of the strand, we succeeded in inducing repolarization which propagated along the strand at a velocity of 1~2 cm/s. The ionic mechanisms responsible for repolarization at the myocyte level, i.e., the deactivation of both the INaL and the L-type calcium current (ICaL), and the activation of the rapid component of delayed rectifier potassium current (IKr) and the inward rectifier potassium channel (IK1), were found to be important for the propagation of repolarization in the myocyte strand. Using an analogy with progressive activation of the sodium current (INa) in the propagation of excitation, regenerative activation of the predominant magnitude of IK1 makes the myocytes at the wave front start repolarization in succession through the electrical coupling via gap junction channels.

Clarifying the composition of the ATP consumption factors required for maintaining ion homeostasis in mouse rod photoreceptors.

To date, no effective treatment has been established for photoreceptor loss due to energy imbalances, but numerous therapeutic approaches have reported some success in slowing photoreceptor degeneration by downregulating energy demand. However, the detailed mechanisms remain unclear. This study aimed to clarify the composition of ATP consumption factors in photoreceptors in darkness and in light. We introduced mathematical formulas for ionic current activities combined with a phototransduction model to form a new mathematical model for estimating the energy expenditure of each ionic current. The proposed model included various ionic currents identified in mouse rods using a gene expression database incorporating an available electrophysiological recording of each specific gene. ATP was mainly consumed by Na+/K+-ATPase and plasma membrane Ca2+-ATPase pumps to remove excess Na+ and Ca2+. The rod consumed 7 [Formula: see text] 107 molecules of ATP s-1, where 65% was used to remove ions from the cyclic nucleotide-gated channel and 20% from the hyperpolarization-activated current in darkness. Increased light intensity raised the energy requirements of the complex phototransduction cascade mechanisms. Nevertheless, the overall energy consumption was less than that in darkness due to the significant reduction in ATPase activities, where the hyperpolarization-activated current proportion increased to 83%. A better understanding of energy demand/supply may provide an effective tool for investigating retinal pathophysiological changes and analyzing novel therapeutic treatments related to the energy consumption of photoreceptors.

Theoretical prediction of early afterdepolarization-evoked triggered activity formation initiating ventricular reentrant arrhythmias.

BACKGROUND AND OBJECTIVE: Excessive prolongation of QT interval on ECGs in patients with congenital/acquired long QT syndrome and heart failure is a sign suggesting the development of early afterdepolarization (EAD), an abnormal repolarization in the action potential of ventricular cardiomyocytes. The development of EAD has been believed to be a trigger for fatal tachyarrhythmia, which can be a risk for sudden cardiac death. The role of EAD in triggering ventricular tachycardia (VT) remains unclear. The aim of this study was to elucidate the mechanism of EAD-induced triggered activity formation that leads to the VT such as Torsades de Pointes. METHODS: We investigated the relationship between EAD and tachyarrhythmia initiation by constructing homogeneous myocardial sheet models consisting of the mid-myocardial cell version of a human ventricular myocyte model and performing simulations of excitation propagation. RESULTS: A solitary island-like (clustering) occurrence of EADs in the homogeneous myocardial sheet could induce a focal excitation wave. However, reentrant excitation, an entity of tachyarrhythmia, was not able to be triggered regardless of the EAD cluster size when the focal excitation wave formed a repolarization potential difference boundary consisting of only a convex surface. The discontinuous distribution of multiple EAD clusters in the ventricular tissue formed a specific repolarization heterogeneity due to the repolarization potential difference, the shape of which depended on EAD cluster size and placed intervals. We found that the triggered activity was formed in such a manner that the repolarization potential difference boundary included a concave surface. CONCLUSIONS: The formation of triggered activity that led to tachyarrhythmia required not only the occurrence of EAD onset-mediated focal excitation wave but also a repolarization heterogeneity-based specific repolarization potential difference boundary shape formed within the tissue.

Gradient-based parameter optimization method to determine membrane ionic current composition in human induced pluripotent stem cell-derived cardiomyocytes.

Premature cardiac myocytes derived from human induced pluripotent stem cells (hiPSC-CMs) show heterogeneous action potentials (APs), probably due to different expression patterns of membrane ionic currents. We developed a method for determining expression patterns of functional channels in terms of whole-cell ionic conductance (Gx) using individual spontaneous AP configurations. It has been suggested that apparently identical AP configurations can be obtained using different sets of ionic currents in mathematical models of cardiac membrane excitation. If so, the inverse problem of Gx estimation might not be solved. We computationally tested the feasibility of the gradient-based optimization method. For a realistic examination, conventional 'cell-specific models' were prepared by superimposing the model output of AP on each experimental AP recorded by conventional manual adjustment of Gxs of the baseline model. Gxs of 4-6 major ionic currents of the 'cell-specific models' were randomized within a range of ± 5-15% and used as an initial parameter set for the gradient-based automatic Gxs recovery by decreasing the mean square error (MSE) between the target and model output. Plotting all data points of the MSE-Gx relationship during optimization revealed progressive convergence of the randomized population of Gxs to the original value of the cell-specific model with decreasing MSE. The absence of any other local minimum in the global search space was confirmed by mapping the MSE by randomizing Gxs over a range of 0.1-10 times the control. No additional local minimum MSE was obvious in the whole parameter space, in addition to the global minimum of MSE at the default model parameter.

Effects of Heat Stress on Heart Rate Variability in Free-Moving Sheep and Goats Assessed With Correction for Physical Activity.

Heart rate variability (HRV) is the heart beat-to-beat variation under control of the cardiovascular function of animals. Under stressed conditions, cardiac activity is generally regulated with an upregulated sympathetic tone and withdrawal of vagal tone; thus, HRV monitoring can be a non-invasive technique to assess stress level in animals especially related to animal welfare. Among several stress-induced factors, heat stress is one of the most serious causes of physiological damage to animals. The aim of this study was to assess the effects of heat stress on HRV in small ruminants under free-moving conditions. In three experimental periods (June, August, and October), inter-beat intervals in sheep and goats (three for each) in two consecutive days were measured. HRV parameters were calculated from the inter-beat interval data by three types of analyses: time domain, frequency domain, and non-linear analyses. The temperature-humidity index (THI) was used as an indicator of heat stress, and vectorial dynamic body acceleration (VeDBA) was calculated to quantify the physical activity of the animals tested. First, we investigated correlations of THI and VeDBA with HRV parameters; subsequently, THI was divided into five categories according to the values obtained (≤ 65, 65-70, 70-75, 75-80, and >80), and the effects of the THI categories on HRV parameters were investigated with and without correcting for the effects of physical activity based on the VeDBA. The results indicated that HRV significantly decreased with increasing THI and VeDBA. For non-linear HRV parameters that were corrected for the effects of physical activity, it was suggested that there would be a threshold of THI around 80 that strongly affected HRV; high heat stress can affect the autonomic balance of animals non-linearly by inducing the sympathetic nervous system. In conclusion, to assess psychophysiological conditions of unrestrained animals by HRV analysis, the confounding effect of physical activity on HRV should be minimized for a more precise interpretation of the results.

Mechanisms Underlying Spontaneous Action Potential Generation Induced by Catecholamine in Pulmonary Vein Cardiomyocytes: A Simulation Study.

Cardiomyocytes and myocardial sleeves dissociated from pulmonary veins (PVs) potentially generate ectopic automaticity in response to noradrenaline (NA), and thereby trigger atrial fibrillation. We developed a mathematical model of rat PV cardiomyocytes (PVC) based on experimental data that incorporates the microscopic framework of the local control theory of Ca2+ release from the sarcoplasmic reticulum (SR), which can generate rhythmic Ca2+ release (limit cycle revealed by the bifurcation analysis) when total Ca2+ within the cell increased. Ca2+ overload in SR increased resting Ca2+ efflux through the type II inositol 1,4,5-trisphosphate (IP3) receptors (InsP3R) as well as ryanodine receptors (RyRs), which finally triggered massive Ca2+ release through activation of RyRs via local Ca2+ accumulation in the vicinity of RyRs. The new PVC model exhibited a resting potential of -68 mV. Under NA effects, repetitive Ca2+ release from SR triggered spontaneous action potentials (APs) by evoking transient depolarizations (TDs) through Na+/Ca2+ exchanger (APTDs). Marked and variable latencies initiating APTDs could be explained by the time courses of the α1- and ß1-adrenergic influence on the regulation of intracellular Ca2+ content and random occurrences of spontaneous TD activating the first APTD. Positive and negative feedback relations were clarified under APTD generation.

Correcting the Activity-Specific Component of Heart Rate Variability Using Dynamic Body Acceleration Under Free-Moving Conditions.

Heart rate variability (HRV) analysis is a widely used technique to assess sympatho-vagal regulation in response to various internal or external stressors. However, HRV measurements under free-moving conditions are highly susceptible to subjects' physical activity levels because physical activity alters energy metabolism, which inevitably modulates the cardiorespiratory system and thereby changes the sympatho-vagal balance, regardless of stressors. Thus, researchers must simultaneously quantify the effect of physical activity on HRV to reliably assess sympatho-vagal balance under free-moving conditions. In the present study, dynamic body acceleration (DBA), which was developed in the field of animal ecology as a quantitative proxy for activity-specific energy expenditure, was used as a factor to correct for physical activity when evaluating HRV in freely moving subjects. Body acceleration and heart inter-beat intervals were simultaneously measured in cattle and sheep, and the vectorial DBA and HRV parameters were evaluated at 5-min intervals. Next, the effects of DBA on the HRV parameters were statistically analyzed. The heart rate (HR) and most of the HRV parameters were affected by DBA in both animal species, and the inclusion of the effect of DBA in the HRV analysis greatly influenced the frequency domain and nonlinear HRV parameters. By removing the effect of physical activity quantified using DBA, we could fairly compare the stress levels of animals with different physical activity levels under different management conditions. Moreover, we analyzed and compared the HRV parameters before and after correcting for the mean HR, with and without inclusion of DBA. The results were somewhat unexpected, as the effect of DBA was a highly significant source of HRV also in parameters corrected for mean HR. In conclusion, the inclusion of DBA as a physical activity index is a simple and useful method for correcting the activity-specific component of HRV under free-moving conditions.

Regulation of the glucose supply from capillary to tissue examined by developing a capillary model.

A new glucose transport model relying upon diffusion and convection across the capillary membrane was developed, and supplemented with tissue space and lymph flow. The rate of glucose utilization (J util) in the tissue space was described as a saturation function of glucose concentration in the interstitial fluid (C glu,isf), and was varied by applying a scaling factor f to J max. With f = 0, the glucose diffusion ceased within ~20 min. While, with increasing f, the diffusion was accelerated through a decrease in C glu,isf, but the convective flux remained close to resting level. When the glucose supplying capacity of the capillary was measured with a criterion of J util /J max = 0.5, the capacity increased in proportion to the number of perfused capillaries. A consistent profile of declining C glu,isf along the capillary axis was observed at the criterion of 0.5 irrespective of the capillary number. Increasing blood flow scarcely improved the supplying capacity.

A simulation study on the constancy of cardiac energy metabolites during workload transition.

KEY POINTS: The cardiac energy metabolites such as ATP, phosphocreatine, ADP and NADH are kept relatively constant during physiological cardiac workload transition. How this is accomplished is not yet clarified, though Ca2+ has been suggested to be one of the possible mechanisms. We constructed a detailed mathematical model of cardiac mitochondria based on experimental data and studied whether known Ca2+ -dependent regulation mechanisms play roles in the metabolite constancy. Model simulations revealed that the Ca2+ -dependent regulation mechanisms have important roles under the in vitro condition of isolated mitochondria where malate and glutamate were mitochondrial substrates, while they have only a minor role and the composition of substrates has marked influence on the metabolite constancy during workload transition under the simulated in vivo condition where many substrates exist. These results help us understand the regulation mechanisms of cardiac energy metabolism during physiological cardiac workload transition. ABSTRACT: The cardiac energy metabolites such as ATP, phosphocreatine, ADP and NADH are kept relatively constant over a wide range of cardiac workload, though the mechanisms are not yet clarified. One possible regulator of mitochondrial metabolism is Ca2+ , because it activates several mitochondrial enzymes and transporters. Here we constructed a mathematical model of cardiac mitochondria, including oxidative phosphorylation, substrate metabolism and ion/substrate transporters, based on experimental data, and studied whether the Ca2+ -dependent activation mechanisms play roles in metabolite constancy. Under the in vitro condition of isolated mitochondria, where malate and glutamate were used as mitochondrial substrates, the model well reproduced the Ca2+ and inorganic phosphate (Pi ) dependences of oxygen consumption, NADH level and mitochondrial membrane potential. The Ca2+ -dependent activations of the aspartate/glutamate carrier and the F1 Fo -ATPase, and the Pi -dependent activation of Complex III were key factors in reproducing the experimental data. When the mitochondrial model was implemented in a simple cardiac cell model, simulation of workload transition revealed that cytoplasmic Ca2+ concentration ([Ca2+ ]cyt ) within the physiological range markedly increased NADH level. However, the addition of pyruvate or citrate attenuated the Ca2+ dependence of NADH during the workload transition. Under the simulated in vivo condition where malate, glutamate, pyruvate, citrate and 2-oxoglutarate were used as mitochondrial substrates, the energy metabolites were more stable during the workload transition and NADH level was almost insensitive to [Ca2+ ]cyt . It was revealed that mitochondrial substrates have a significant influence on metabolite constancy during cardiac workload transition, and Ca2+ has only a minor role under physiological conditions.

Mechanisms underlying the volume regulation of interstitial fluid by capillaries: a simulation study.

BACKGROUND: Control of the extracellular fluid volume is one of the most indispensable issues for homeostasis of the internal milieu. However, complex interdependence of the pressures involved in determination of fluid exchange makes it difficult to predict a steady-state tissue volume under various physiological conditions without mathematical approaches. METHODS: Here, we developed a capillary model based on the Starling's principle, which allowed us to clarify the mechanisms of the interstitial-fluid volume regulation. Three well known safety factors against edema: (1) low tissue compliance in negative pressure ranges; (2) lymphatic flow driven by the tissue pressure; and (3) protein washout by the lymph, were incorporated into the model in sequence. RESULTS: An increase in blood pressure at the venous end of the capillary induced an interstitial-fluid volume increase, which, in turn, reduced negative tissue pressure to prevent edema. The lymphatic flow alleviated the edema by both carrying fluid away from the tissue and decreasing the colloidal osmotic pressure. From the model incorporating all three factors, we found that the interstitial-fluid volume changed quickly after the blood pressure change, and that the protein movement towards a certain equilibrium point followed the volume change. CONCLUSION: Mathematical analyses revealed that the system of the capillary is stable near the equilibrium point at steady state and normal physiological capillary pressure. The time course of the tissue-volume change was determined by two kinetic mechanisms: rapid fluid exchange and slow protein fluxes.

A human ventricular myocyte model with a refined representation of excitation-contraction coupling.

Cardiac Ca(2+)-induced Ca(2+) release (CICR) occurs by a regenerative activation of ryanodine receptors (RyRs) within each Ca(2+)-releasing unit, triggered by the activation of L-type Ca(2+) channels (LCCs). CICR is then terminated, most probably by depletion of Ca(2+) in the junctional sarcoplasmic reticulum (SR). Hinch et al. previously developed a tightly coupled LCC-RyR mathematical model, known as the Hinch model, that enables simulations to deal with a variety of functional states of whole-cell populations of a Ca(2+)-releasing unit using a personal computer. In this study, we developed a membrane excitation-contraction model of the human ventricular myocyte, which we call the human ventricular cell (HuVEC) model. This model is a hybrid of the most recent HuVEC models and the Hinch model. We modified the Hinch model to reproduce the regenerative activation and termination of CICR. In particular, we removed the inactivated RyR state and separated the single step of RyR activation by LCCs into triggering and regenerative steps. More importantly, we included the experimental measurement of a transient rise in Ca(2+) concentrations ([Ca(2+)], 10-15 µM) during CICR in the vicinity of Ca(2+)-releasing sites, and thereby calculated the effects of the local Ca(2+) gradient on CICR as well as membrane excitation. This HuVEC model successfully reconstructed both membrane excitation and key properties of CICR. The time course of CICR evoked by an action potential was accounted for by autonomous changes in an instantaneous equilibrium open probability of couplons. This autonomous time course was driven by a core feedback loop including the pivotal local [Ca(2+)], influenced by a time-dependent decay in the SR Ca(2+) content during CICR.

Ionic mechanisms and Ca2+ dynamics underlying the glucose response of pancreatic ß cells: a simulation study.

To clarify the mechanisms underlying the pancreatic ß-cell response to varying glucose concentrations ([G]), electrophysiological findings were integrated into a mathematical cell model. The Ca(2+) dynamics of the endoplasmic reticulum (ER) were also improved. The model was validated by demonstrating quiescent potential, burst-interburst electrical events accompanied by Ca(2+) transients, and continuous firing of action potentials over [G] ranges of 0-6, 7-18, and >19 mM, respectively. These responses to glucose were completely reversible. The action potential, input impedance, and Ca(2+) transients were in good agreement with experimental measurements. The ionic mechanisms underlying the burst-interburst rhythm were investigated by lead potential analysis, which quantified the contributions of individual current components. This analysis demonstrated that slow potential changes during the interburst period were attributable to modifications of ion channels or transporters by intracellular ions and/or metabolites to different degrees depending on [G]. The predominant role of adenosine triphosphate-sensitive K(+) current in switching on and off the repetitive firing of action potentials at 8 mM [G] was taken over at a higher [G] by Ca(2+)- or Na(+)-dependent currents, which were generated by the plasma membrane Ca(2+) pump, Na(+)/K(+) pump, Na(+)/Ca(2+) exchanger, and TRPM channel. Accumulation and release of Ca(2+) by the ER also had a strong influence on the slow electrical rhythm. We conclude that the present mathematical model is useful for quantifying the role of individual functional components in the whole cell responses based on experimental findings.

Minor contribution of cytosolic Ca2+ transients to the pacemaker rhythm in guinea pig sinoatrial node cells.

The question of the extent to which cytosolic Ca(2+) affects sinoatrial node pacemaker activity has been discussed for decades. We examined this issue by analyzing two mathematical pacemaker models, based on the "Ca(2+) clock" (C) and "membrane clock" (M) hypotheses, together with patch-clamp experiments in isolated guinea pig sinoatrial node cells. By applying lead potential analysis to the models, the C mechanism, which is dependent on potentiation of Na(+)/Ca(2+) exchange current via spontaneous Ca(2+) release from the sarcoplasmic reticulum (SR) during diastole, was found to overlap M mechanisms in the C model. Rapid suppression of pacemaker rhythm was observed in the C model by chelating intracellular Ca(2+), whereas the M model was unaffected. Experimental rupturing of the perforated-patch membrane to allow rapid equilibration of the cytosol with 10 mM BAPTA pipette solution, however, failed to decrease the rate of spontaneous action potential within ∼30 s, whereas contraction ceased within ∼3 s. The spontaneous rhythm also remained intact within a few minutes when SR Ca(2+) dynamics were acutely disrupted using high doses of SR blockers. These experimental results suggested that rapid disruption of normal Ca(2+) dynamics would not markedly affect spontaneous activity. Experimental prolongation of the action potentials, as well as slowing of the Ca(2+)-mediated inactivation of the L-type Ca(2+) currents induced by BAPTA, were well explained by assuming Ca(2+) chelation, even in the proximity of the channel pore in addition to the bulk cytosol in the M model. Taken together, the experimental and model findings strongly suggest that the C mechanism explicitly described by the C model can hardly be applied to guinea pig sinoatrial node cells. The possible involvement of L-type Ca(2+) current rundown induced secondarily through inhibition of Ca(2+)/calmodulin kinase II and/or Ca(2+)-stimulated adenylyl cyclase was discussed as underlying the disruption of spontaneous activity after prolonged intracellular Ca(2+) concentration reduction for >5 min.

A novel method to quantify contribution of channels and transporters to membrane potential dynamics.

The action potential, once triggered in ventricular or atrial myocytes, automatically proceeds on its time course or is generated spontaneously in sinoatrial node pacemaker cells. It is induced by complex interactions among such cellular components as ion channels, transporters, intracellular ion concentrations, and signaling molecules. We have developed what is, to our knowledge, a new method using a mathematical model to quantify the contribution of each cellular component to the automatic time courses of the action potential. In this method, an equilibrium value, which the membrane potential is approaching at a given moment, is calculated along the time course of the membrane potential. The calculation itself is based on the time-varying conductance and the reversal potentials of individual ion channels and electrogenic ion transporters. Since the equilibrium potential moves in advance of the membrane potential change, we refer to it as the lead potential, V(L). The contribution of an individual current was successfully quantified by comparing dV(L)/dt before and after fixing the time-dependent change of a component of interest, such as the variations in the open probability of a channel or the turnover rate of an ion transporter. In addition to the action potential, the lead-potential analysis should also be applicable in all types of membrane excitation in many different kinds of cells.

Ovariectomy increases neuronal amyloid-beta binding alcohol dehydrogenase level in the mouse hippocampus.

Ovarian hormone decline after menopause may influence cognitive performance and increase the risk for Alzheimer's disease (AD) in women. Amyloid-beta peptide (Abeta) has been proposed to be the primary cause of AD. In this study, we examined whether ovariectomy (OVX) could affect the levels of cofactors Abeta-binding alcohol dehydrogenase (ABAD) and receptor for advanced glycation endproducts (RAGE), which have been reported to potentiate Abeta-mediated neuronal perturbation, in mouse hippocampus, correlating with estrogen and Abeta levels. Female ICR mice were randomly divided into ovariectomized or sham-operated groups, and biochemical analyses were carried out at 5 weeks after the operation. OVX for 5 weeks significantly decreased hippocampal 17beta-estradiol level, while it tended to reduce the hormone level in serum, compared with the sham-operated control. In contrast, OVX did not affect hippocampal Abeta(1-40) level, although it significantly increased serum Abeta(1-40) level. Furthermore, we demonstrated that OVX increased hippocampal ABAD level in neurons, but not astrocytes, while it did not affect RAGE level. These findings suggest that the expression of neuronal ABAD depends on estrogen level in the hippocampus and the increase in serum Abeta and hippocampal ABAD induced by ovarian hormone decline may be associated with pre-stage of memory deficit in postmenopausal women and Abeta-mediated AD pathology.

Enhanced activity of hippocampal BACE1 in a mouse model of postmenopausal memory deficits.

Ovarian hormone decline after menopause may influence cognitive performance and increase the risk for Alzheimer's disease (AD) in women. We have recently demonstrated that a combination of ovariectomy and chronic stress (OVX/stress) causes hippocampus-associated cognitive dysfunction in mice. In this study, we examined whether OVX/stress could affect the levels of AD-related molecules in the mouse hippocampus. Female ICR mice were ovariectomized or sham-operated, and then randomly divided into a daily restraint stress (21 days, 6 h/day) or non-stress group. Although OVX or stress alone did not affect beta-site amyloid precursor protein (APP)-cleaving enzyme-1 (BACE1) activity, OVX/stress increased activity in hippocampal CA1 and CA3 regions, compared with other groups. In contrast, OVX/stress did not affect gamma-secretase activity, Abeta(1-40), and phosphorylated-tau levels in the hippocampus. These findings suggest that a stressful life after menopause can influence the levels of AD-related molecules and that BACE1 is the most sensitive molecule for such a situation.

Ionic mechanisms underlying the positive chronotropy induced by beta1-adrenergic stimulation in guinea pig sinoatrial node cells: a simulation study.

Positive chronotropy induced by beta1-adrenergic stimulation is achieved by multiple interactions of ion channels and transporters in sinoatrial node pacemaker cells (SANs). To investigate the ionic mechanisms, we updated our SAN model developed in 2003 and incorporated the beta1-adrenergic signaling cascade developed by Kuzumoto et al. (2007). Since the slow component of the delayed rectifier K+ current (IKs) is one of the major targets of the beta1-adrenergic cascade, we developed a guinea pig model with a large IKs. The new model provided a good representation of the experimental characteristics of SANs. A comparison of individual current during diastole recorded before and after beta1-adrenergic stimulation clearly showed the negative shift of the L-type Ca2+ current (ICaL) takeoff potential, enlargement of the sustained inward current (I st), and the hyperpolarization-activated nonselective cation current (Iha) played major roles in increasing the firing frequency. Deactivation of IKs during diastole scarcely contributed to the time-dependent decrease in membrane K+ conductance, which was the major mechanism for slow diastolic depolarization, as indicated by calculating the instantaneous equilibrium potential (lead potential). This was because the activation of IKs during the preceding action potential was negligibly small. However, IKs was important in counterbalancing the increase in ICaL and the Na+/Ca2+ exchange current (INaCa), which otherwise compromised the positive chronotropic effect by elongating the action potential duration. Enhanced Ca2+ release from the sarcoplasmic reticulum failed to induce an obvious chronotropic effect in our model.