WHOLE-CELL CLAMP STUDY OF XENOPUS EMBRYONIC CHOLINERGIC NEURONS

Whole-cell clamped myoballs are placed into direct visible contact with the growth cones of isolated neurons in embryonic Xenopus culture to serve as probe of acetylchollne (AcCHo) release in order to determine whether these neurons are cholinergic or not. Using a GQ-seal, whole-cell recording method, the electrophysiological properties of these identified cholinergic neurons are studied. It is found that these embryonic neurons, like adult frog motor neurons, exhibit repetitive firings in a certain embryonic developing stage. A development of repetitive firings is observed simultaneously. Tracing the development of one neuron, we find that the development of repetitive firing is completed at the 48th h after fertilization. Tetrodotoxin (TTX) which blocks Na~+ channels can abolish all firings; and tetraethyl ammonium chloride (TEA), the blocker of K~+ channels, reverses this development, i. e. it makes the repetitive firings disappear again. These data show that the nature of the development of repetiti     

R-leaping: accelerating the stochastic simulation algorithm by reaction leaps

A novel algorithm is proposed for the acceleration of the exact stochastic simulation algorithm by a predefined number of reaction firings (R-leaping) that may occur across several reaction channels. In the present approach, the numbers of reaction firings are correlated binomial distributions and the sampling procedure is independent of any permutation of the reaction channels. This enables the algorithm to efficiently handle large systems with disparate rates, providing substantial computational savings in certain cases. Several mechanisms for controlling the accuracy and the appearance of negative species are described. The advantages and drawbacks of R-leaping are assessed by simulations on a number of benchmark problems and the results are discussed in comparison with established methods.     

Effects of Muscle Electrical Activity on the Transmission of Developing Neuromuscular Junction

Miniature endplate potentials (MEPPS) caused by the spontaneous release of ACh from the growth cone of cholinergic neurons, are recorded by the whole-cell patch-clamp technique on a large number of 1-day cultured myoballs which have contact neurites of co-cultured neurons. Both muscle cell and neuron are dissociated from the 1-day-old (about stage 20) Xenopus embryo. Frequency and/or amplitude of MEPPs can obviously increase after the repetitive high-level depolarization caused by the stimuli on muscle cells. No detectable changes of single ACh receptor channel property are observed by using the single-channel recording technique. These results suggest that the mechanism of the increase of MEPPs after electrical activity of postsynaptic muscle cells probably involve some alteration of presynaptic membrane.     

单向耦合Hindmarsh-Rose系统的同步

采用数值模拟方法, 研究了由两个具有不同初始条件单向耦合的Hindmarsh-Rose神经元所组成系统的动力学行为. 研究结果表明, 耦合强度对两个神经元的同步具有重要的影响. 当两个神经元的控制参数取值不同时, 即在驱动系统和响应系统均处于混沌态及驱动系统处于周期振荡态, 响应系统处于混沌态两种情况下, 随耦合强度的增加其放电活动都能从不同步达到相位同步, 最后实现近似完全同步.     

Artificial ion channel formed by cucurbit[n]uril derivatives with a carbonyl group fringed portal reminiscent of the selectivity filter of K+ channels

Novel artificial ion channels (1 and 2) based on CB[n] (n = 6 and 5, respectively) synthetic receptors with carbonyl-fringed portals (diameter 3.9 and 2.4 A, respectively) can transport proton and alkali metal ions across a lipid membrane with ion selectivity. Fluorometric experiments using large unilamellar vesicles showed that 1 mediates proton transport across the membranes, which can be blocked by a neurotransmitter, acetylcholine, reminiscent of the blocking of the K+ channels by polyamines. The alkali metal ion transport activity of 1 follows the order of Li+ > Cs+ approximately Rb+ > K+ > Na+, which is opposite to the binding affinity of CB[6] toward alkali metal ions. On the other hand, the transport activity of 2 follows the order of Li+ > Na+, which is also opposite to the binding affinity of 2 toward these metal ions, but virtually no transport was observed for K+, Rb+, and Cs+. It is presumably because the carbonyl-fringed portal size of 2 (diameter 2.4 A) is smaller than the diameters of these alkali metal ions. To determine the transport mechanism, voltage-clamp experiments on planar bilayer lipid membranes were carried out. The experiments showed that a single-channel current of 1 for Cs+ transport is approximately 5 pA, which corresponds to an ion flux of approximately 3 x 107 ions/s. These results are consistent with an ion channel mechanism. Not only the structural resemblance to the selectivity filter of K+ channels but also the remarkable ion selectivity makes this model system unique.     

Calcium channels, potassium channels and membrane potential of smooth muscle cells of human allantochorial placental vessels.

The membrane potential (Um), the main factor of the excitation-contraction coupling, of human allantochorial placental vascular smooth muscle cells (VSMCs) has been previously shown to depend on voltage-sensitive K+ channels. These channels were blocked by high external K+. To characterize other channels which regulated Um, various constrictor or/and vasodilators and channel blockers were used. Serotonin depolarized VSMCs, in normal medium, but induced a more marked depolarization in VSMCs predepolarized by high external K+. This depolarization was inhibited by nifedipine, a blocker of voltage-gated Ca2+ channels. Acetylcholine, sodium nitroprusside (without effect on Um in normal medium), hyperpolarized the predepolarized-high K+ medium VSMCs. This hyperpolarization was inhibited after addition of charybotoxin (a blocker of Ca2+-activated K+ channels) or/and glibenclamide (a blocker of ATP-sensitive K+ channels). A similar effect was obtained with isoproterenol. These results indicated that membrane potential of human placental allantochorial VSMCs was regulated by voltage-gated, Ca2+- and ATP-sensitive K+ channels and by voltage-dependent Ca2+ channels.     

N-甲硝胺异构化和分解反应机理和动力学的理论研究

The complex potential energy surface and reaction mechanisms for the unimolecular isomerization and decomposition of methyl-nitramine (CH3NHNO2) were theoretically probed at the QCISD(T)/6-311+G*//B3LYP/6-311+G* level of theory. The results demonstrated that there are four low-lying energy channels: (i) the NN bond fission pathway; (ii) a sequence of isomerization reactions via CH3NN(OH)O; (IS2a); (iii) the HONO elimination pathway; (iv) the isomerization and the dissociation reactions via CH3NHONO (IS3). The rate constants of each initial step (rate-determining step) for these channels were calculated using the canonical transition state theory. The Arrhenius expressions of the channels over the temperature range 298-2000 K are k6(T)=1014:8e-46:0=RT , k7(T)=1013:7e-42:1=RT , k8(T)=1013:6e-51:8=RT and k9(T)=1015:6e-54:3=RT s-1, respectively. The calculated overall rate constants is 6.9￡10-4 at 543 K, which is in good agreement with the experimental data. Based on the analysis of the rate constants, the dominant pathway is the isomerization reaction to form CH3NN(OH)O at low temperatures, while the NN bond fission and the isomerization reaction to produce CH3NHONO are expected to be competitive with the isomerization reaction to form CH3NN(OH)O at high temperatures.     

Theoretical study on the OH + CH3NHCOOCH3 reaction

The multiple-channel reactions OH + CH3NHC(O)OCH3 --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-311+G(d,p) level, and energetic information is further refined by the BMC-CCSD (single-point) method. The rate constants for every reaction channels, R1, R2, R3, and R4, are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-1000 K. The total rate constants are in good agreement with the available experimental data and the two-parameter expression k(T) = 3.95 x 10(-12) exp(15.41/T) cm3 molecule(-1) s(-1) over the temperature range 200-1000 K is given. Our calculations indicate that hydrogen abstraction channels R1 and R2 are the major channels due to the smaller barrier height among four channels considered, and the other two channels to yield CH3NC(O)OCH3 + H2O and CH3NHC(O)(OH)OCH3 + H2O are minor channels over the whole temperature range.     

Cell-based fluorescence screen for K+ channels and transporters using an extracellular triazacryptand-based K+ sensor

K+ channels and K+-coupled membrane transporters are important targets for drug discovery. We previously developed a triazacryptand (TAC)-based K+ sensor, TAC-Red, and demonstrated its utility to image K+ waves in mouse brain in vivo (Padmawar et al. Nat. Methods. 2005, 2, 825-827). Here, we synthesized a green-fluorescing dextran conjugate of TAC-bodipy ("TAC-Limedex") for use as an extracellular K+ sensor and demonstrated its utility in measuring K+ transport across cell membranes. TAC-Limedex fluorescence increased by 50% with increasing [K+] from 0 to 2 mM and was insensitive to [Na+], [Cl-], or pH. K+ efflux from cells was quantified from increasing extracellular TAC-Limedex fluorescence following cell immersion in K+-free buffer. In HT-29 cells, K+ efflux was 2.0 +/- 0.1 micromol/cm2/s, increasing 8-fold following K+ channel activation by ATP; the increase in K+ efflux was inhibited by a K+ channel blocker or by preventing cytoplasmic calcium elevation. Electroneutral K+/Cl- cotransport was demonstrated in SiHa cells, in which K+ efflux was increased 3-fold by hypotonic challenge; the increase in K+ efflux was fully inhibited by a K+/Cl- transport blocker. K+ efflux measurements were adapted to a commercial fluorescence platereader for automated screening. The fluorescence-based K+ transport assay largely replaces assays requiring radioactive rubidium and is suitable for high-throughput identification of K+ transport modulators.     

Spider venoms: a rich source of acylpolyamines and peptides as new leads for CNS drugs

Advances in NMR and mass spectrometry as well as in peptide biochemistry coupled to modern methods in electrophysiology have permitted the isolation and identification of numerous products from spider venoms, previously explored due to technical limitations. The chemical composition of spider venoms is diverse, ranging from low molecular weight organic compounds such as acylpolyamines to complex peptides. First, acylpolyamines (< 1000 Da) have an aromatic moiety linked to a hydrophilic lateral chain. They were characterized for the first time in spider venoms and are ligand-gated ion channel antagonists, which block mainly postsynaptic glutamate receptors in invertebrate and vertebrate nervous systems. Acylpolyamines represent the vast majority of organic components from the spider venom. Acylpolyamine analogues have proven to suppress hippocampal epileptic discharges. Moreover, acylpolyamines could suppress excitatory postsynaptic currents inducing Ca+ accumulation in neurons leading to protection against a brain ischemic insult. Second, short spider peptides (< 6000 Da) modulate ionic currents in Ca2+, Na+, or K+ voltage-gated ion channels. Such peptides may contain from three to four disulfide bridges. Some spider peptides act specifically to discriminate among Ca2+, Na+, or K+ ion channel subtypes. Their selective affinities for ion channel subfamilies are functional for mapping excitable cells. Furthermore, several of these peptides have proven to hyperpolarize peripheral neurons, which are associated with supplying sensation to the skin and skeletal muscles. Some spider N-type calcium ion channel blockers may be important for the treatment of chronic pain. A special group of spider peptides are the amphipathic and positively charged peptides. Their secondary structure is alpha-helical and they insert into the lipid cell membrane of eukaryotic or prokaryotic cells leading to the formation of pores and subsequently depolarizing the cell membrane. Acylpolyamines and peptides from spider venoms represent an interesting source of molecules for the design of novel pharmaceutical drugs.     

Patch-clamp recording of human retinal photoreceptors and bipolar cells

Photoreceptors and retinal bipolar cells are considered as nonspiking neurons; however, we recently showed that human rod photoreceptors can generate sodium action potentials in response to membrane depolarization from membrane potentials of -60 or -70 mV (Kawai et al., Neuron 30 [2001] 451). We performed patch-clamp recording of human cone photoreceptors and retinal bipolar cells to examine whether functional voltage-gated sodium channels are expressed in these cells as well as rod photoreceptors. Under current-clamp conditions, the injection of depolarizing current steps into a cone photoreceptor-induced marked action potentials. These action potentials were blocked by 1 microM tetrodotoxin, a voltage-gated sodium channel blocker. Under voltage-clamp conditions, depolarizing voltage steps-induced a fast transient inward current in several bipolar cells (n = 4/78). This current was activated from -70 to + 20 mV (maximal at -10 mV) and inactivated within 5 ms. The 10-90% rise time of this current was shorter than another inward current (less than one-hundredth). These results indicate that human cones and bipolar cells express voltage-gated sodium channels as rod photoreceptors. Sodium channels may serve to amplify the release of a neurotransmitter and to accelerate the light-dark change in photosignals.     

Theoretical study on the Br + CH3SCH3 reaction

The multiple-channel reactions Br + CH(3)SCH(3) --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-31+G(d,p) level, and energetic information is further refined by the G3(MP2) (single-point) theory. The rate constants for every reaction channels, Br + CH(3)SCH(3) --> CH(3)SCH(2) + HBr (R1), Br + CH(3)SCH(3) --> CH(3)SBr + CH(3) (R2), and Br + CH(3)SCH(3) -->CH(3)S + CH(3)Br (R3), are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-3000 K. The total rate constants are in good agreement with the available experimental data, and the two-parameter expression k(T) = 2.68 x 10(-12) exp(-1235.24/T) cm(3)/(molecule s) over the temperature range 200-3000 K is given. Our calculations indicate that hydrogen abstraction channel is the major channel due to the smallest barrier height among three channels considered, and the other two channels to yield CH(3)SBr + CH(3) and CH(3)S + CH(3)Br are minor channels over the whole temperature range.     

Role of fluctuations in a snug-fit mechanism of KcsA channel selectivity

The thermodynamic exclusion of Na+ relative to K+ in potassium channels is examined by calculating the distribution of binding energies for Na+ and K+ in a model of the selectivity filter of the KcsA potassium channel. These distributions are observed to take a surprisingly simple form: Gaussian with a slight positive skewness that is insignificant in the present context. Complications that might be anticipated from these distributions are not problematic here. Na+ occupies the filter with a mean binding energy substantially lower than that of K+. The difference is comparable to the difference in hydration free energies of Na+ and K+ in bulk aqueous solution. Thus, the average energies of binding to the filter do not discriminate Na+ from K+ when measured from a baseline of the difference in bulk hydration free energies. The strong binding of Na+ constricts the filter, consistent with a negative partial molar volume of Na+ in water in contrast with a positive partial molar volume of K+ in water. Discrimination in favor of K+)can be attributed to the scarcity of favorable binding configurations for Na+ compared to K+. That relative scarcity is quantified as enhanced binding energy fluctuations, which reflects both the energetically stronger binding of Na+ and the constriction of the filter induced by Na+ binding.     

Quantum chemical and theoretical kinetics study of the O(3P) + CS2 reaction

The triplet potential energy surface of the O((3)P) + CS(2) reaction is investigated by using various quantum chemical methods including CCSD(T), QCISD(T), CCSD, QCISD, G3B3, MPWB1K, BB1K, MP2, and B3LYP. The thermal rate coefficients for the formation of three major products, CS + SO ((3)Σ(-)), OCS + S ((3)P) and CO + S(2) ((3)Σ(-)(g)) were computed by using transition state and RRKM statistical rate theories over the temperature range of 200-2000 K. The computed k(SO + CS) by using high-level quantum chemical methods is in accordance with the available experimental data. The calculated rate coefficients for the formation of OCS + S ((3)P) and CO + S(2) ((3)Σ(-)(g)) are much lower than k(SO + CS); hence, it is predicted that these two product channels do not contribute significantly to the overall rate coefficient.     

Lead interferes with calcium entry through membrane pores

Calcium is an important intracellular messenger in all cells, represented here by nerve cells and osteoblast-like (OBL) cells. In neurons the intracellular calcium signal is related, e.g., to bioelectric phenomena. In OBL cells the intracellular calcium concentration ([Ca²⁺]_i) plays a role in the intercellular communication via gap junction channels. [Ca²⁺]_i might be affected by lead (Pb²⁺). In the nervous system even low Pb²⁺ concentrations impair learning and memory functions. Considering long-term potentiation (LTP) as a model for learning and memory it has been proven that the generation and maintenance of LTP is reduced by Pb²⁺ (1–10 μM). As the induction of LTP depends on a rise of [Ca²⁺]_i, we examined the effects of Pb²⁺ on [Ca²⁺]_i and on currents through calcium permeable membrane pores in dorsal-root ganglion (DRG) neurons, using calcium measurements (Fura-2/ AM) and whole cell patch clamp techniques. To study the effects of Pb²⁺ on intercellular communication via gap junctions we used rat OBL cells investigating interactions of Pb²⁺ with electric cell coupling. Furthermore, we examined calcium release activated channel currents (CRACCs) of these cells. Lead (1–10 μM) reduced the stimulated increase of [Ca²⁺]_i in a concentration dependent manner, by reducing both voltage-activated calcium channels (VACCs) and N-methyl-D-aspartate activated calcium channels (NACCs) in neurons. Voltage-activated calcium channel currents (VACCCs) were reduced by Pb²⁺ with an IC₅₀ of 0.46 μM. The effect was quite specific as voltage activated sodium and potassium channel currents were not significantly altered in the same concentration and voltage range. Furthermore, this effect was not voltage dependent and only partly reversible. A 100-fold higher concentration of Pb²⁺ (IC₅₀ of 46 μM) was found for the reduction of NACC currents. A small portion of this effect was not reversible. Other agonist activated channel currents (kainate and quisqualate) are not affected. In OBL cells, the calcium entry through calcium release activated channels (CRACs) was reduced in a concentration dependent manner by extracellular Pb²⁺, the concentrations were between 2 and 20 μM. Surprisingly the electric coupling through gap junction channels in OBL cells was not reduced by either extracellular or intracellular Pb²⁺ (5–25 μM). Received: 1 August 1997 / Revised: 24 October 1997 / Accepted: 30 October 1997     

Master equation methods for multiple well systems: application to the 1-,2-pentyl system

The master equation (ME) provides a powerful technique for modeling reactions that involve at least one potential energy well. It can be widely applied to reactions with several connected energy wells and multiple product channels. The application of the technique is reviewed by reference to the H + SO(2) reaction, where phenomenological rate constants for use, for example, in a combustion model can be extracted through an analysis of the eigenvalues and eigenvectors of the collision matrix, M, that describes formation of the adducts HSO(2) and HOSO from the source H + SO(2), collisional energy transfer in the adduct wells and reaction via the product channel (sink) OH + SO. The approach is extended to systems with more than one sink and it is demonstrated that macroscopic (phenomenological) rate coefficients derived from a ME obey detailed balance if the original ME is appropriately constructed. The method has been applied to the 1-, 2-pentyl radical system, that includes isomerisation and dissociation via two channels to form C(3)H(6) + C(2)H(5) and C(2)H(4) + C(3)H(7). The calculations clearly demonstrate the importance of indirect dissociation channels, in which an isomer can dissociate to form the product set to which it is not directly connected, e.g. formation of C(3)H(6) + C(2)H(5) from 1-pentyl, via the energized states of 2-pentyl. As in previous studies of pentyl dissociation, there is a convergence of the chemically significant eigenvalues and the internal energy relaxation eigenvalues above approximately 1000 K; the consequences of this convergence are discussed.     

K-leap method for accelerating stochastic simulation of coupled chemical reactions

Leap methods are very promising for accelerating stochastic simulation of a well stirred chemically reacting system, while providing acceptable simulation accuracy. In Gillespie's tau-leap method [D. Gillespie, J. Phys. Chem. 115, 1716 (2001)], the number of firings of each reaction channel during a leap is a Poisson random variable, whose sample values are unbounded. This may cause large changes in the populations of certain molecular species during a leap, thereby violating the leap condition. In this paper, we develop an alternative leap method called the K-leap method, in which we constrain the total number of reactions occurring during a leap to be a number K calculated from the leap condition. As the number of firings of each reaction channel during a leap is upper bounded by a properly chosen number, our K-leap method can better satisfy the leap condition, thereby improving simulation accuracy. Since the exact stochastic simulation algorithm (SSA) is a special case of our K-leap method when K=1, our K-leap method can naturally change from the exact SSA to an approximate leap method during simulation, whenever the leap condition allows to do so.     

High-temperature rate constants for CH3OH + Kr --> products, OH + CH3OH --> products, OH + (CH3)(2)CO --> CH2COCH3 + H2O, and OH + CH3 --> CH) + H2O

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the dissociation of methanol between 1591 and 2865 K. Rate constants for two product channels [CH3OH + Kr --> CH3 + OH + Kr (1) and CH3OH + Kr --> 1CH2 + H2O + Kr (2)] were determined. During the course of the study, it was necessary to determine several other rate constants that contributed to the profile fits. These include OH + CH3OH --> products, OH + (CH3)2CO --> CH2COCH3 + H2O, and OH + CH3 --> 1,3CH2 + H2O. The derived expressions, in units of cm(3) molecule(-1) s(-1), are k(1) = 9.33 x 10(-9) exp(-30857 K/T) for 1591-2287 K, k(2) = 3.27 x 10(-10) exp(-25946 K/T) for 1734-2287 K, kOH+CH3OH = 2.96 x 10-16T1.4434 exp(-57 K/T) for 210-1710 K, k(OH+(CH3)(2)CO) = (7.3 +/- 0.7) x 10(-12) for 1178-1299 K and k(OH+CH3) = (1.3 +/- 0.2) x 10(-11) for 1000-1200 K. With these values along with other well-established rate constants, a mechanism was used to obtain profile fits that agreed with experiment to within <+/-10%. The values obtained for reactions 1 and 2 are compared with earlier determinations and also with new theoretical calculations that are presented in the preceding article in this issue. These new calculations are in good agreement with the present data for both (1) and (2) and also for OH + CH3 --> products.