Observation of magnetically induced transparency in a classical magnetized plasma.

We report the first demonstration of magnetically induced transmission in an opaque magnetized plasma. Magnetically induced transmission in a plasma is a classical analog to the electromagnetically induced transparency in atomic systems. The transmission of radiation through an axially magnetized plasma is obtained by applying an additional one dimensional transverse spatial periodic magnetic field. The transverse-periodic magnetic field uncouples the right-hand electromagnetic wave from interacting with plasma electrons, rendering the plasma band-stop transparent. This provides means to control the extent of absorption of electromagnetic radiation in magnetized plasma.

Resonant behavior of dielectric objects (electrostatic resonances).

Resonant behavior of dielectric objects occurs at certain frequencies for which the object permittivity is negative and the free-space wavelength is large in comparison with the object dimensions. Unique physical features of these resonances are studied and a novel technique for the calculation of resonance values of permittivity, and hence resonance frequencies, is proposed. Scale invariance of resonance frequencies, unusually strong orthogonality properties of resonance modes, and a two-dimensional phenomenon of "twin" spectra are reported. The paper concludes with brief discussions of optical controllability of these resonances in semiconductor nanoparticles and a plausible, electrostatic resonance based, mechanism for nucleation and formation of ball lightning.

Hydrodynamic flow of ions and atoms in partially ionized plasmas.

We have derived the hydrodynamic equations of motion for a partially ionized plasma, when the ionized component and the neutral components have different flow velocities and kinetic temperatures. Starting from the kinetic equations for a gas of ions and a gas of atoms we have considered various processes of encounters between the two species: self-collisions, interspecies collisions, ionization, recombination, and charge exchange. Our results were obtained by developing a general approach for the hydrodynamics of a gas in a binary mixture, in particular when the components drift with respect to each other. This was applied to a partially ionized plasma, when the neutral-species gas and the charged-species gas have separate velocities. We have further suggested a generalized version of the relaxation time approximation and obtained the contributions of the interspecies encounters to the transport equations.

Persistence: a new statistic for characterizing ion-channel activity.

We introduce and illustrate by examples a new statistical technique, the persistence function, for characterizing ion-channel activity in a single-channel patch-clamp recording. Persistence is a function of both current and time. It is the probability that the current is at a given level (conditional on it having been at that level at an earlier time). Viewed as a function of current it exhibits the prominent conductance levels present in the recording, and viewed as a function of time for a conductance level it portrays the kinetics at that level.

Ion channel formation by synthetic analogues of staphylococcal delta-toxin.

Ion channel formation by three analogues of staphylococcal delta-toxin, an amphipathic and alpha-helical channel-forming peptide, has been evaluated by measurement of ionic currents across planar lipid bilayers. Replacement of beta-branched, hydrophobic residues by leucine and movement of a tryptophan residue from the hydrophilic to the hydrophobic face of the helix does not significantly alter ion channel activity. Removal of the N-terminal blocking group combined with the substitution of glycine-10 by leucine changes the single channel properties of delta-toxin, without altering macroscopic conductance/voltage behaviour. Truncation of the N-terminus by three residues results in complete loss of channel-forming activity. These changes in channel-forming properties upon altering the peptide sequence do not mirror changes in haemolytic activity. The results lend support to the proposal that channel formation and haemolysis are distinct events. Channel properties are discussed in the context of a model in which the pore is formed by a bundle of approximately parallel transbilayer helices.

Electrostatic calculations of amino acid titration and electron transfer, Q-AQB-->QAQ-B, in the reaction center.

The titration of amino acids and the energetics of electron transfer from the primary electron acceptor (QA) to the secondary electron acceptor (QB) in the photosynthetic reaction center of Rhodobacter sphaeroides are calculated using a continuum electrostatic model. Strong electrostatic interactions between titrating sites give rise to complex titration curves. Glu L212 is calculated to have an anomalously broad titration curve, which explains the seemingly contradictory experimental results concerning its pKa. The electrostatic field following electron transfer shifts the average protonation of amino acids near the quinones. The pH dependence of the free energy between Q-AQB and QAQ-B calculated from these shifts is in good agreement with experiment. However, the calculated absolute free energy difference is in severe disagreement (by approximately 230 meV) with the observed experimental value, i.e., electron transfer from Q-A to QB is calculated to be unfavorable. The large stabilization energy of the Q-A state arises from the predominantly positively charged residues in the vicinity of QA in contrast to the predominantly negatively charged residues near QB. The discrepancy between calculated and experimental values for delta G(Q-AQB-->QAQ-B) points to limitations of the continuum electrostatic model. Inclusion of other contributions to the energetics (e.g., protein motion following quinone reduction) that may improve the agreement between theory and experiment are discussed.

Maximum likelihood estimation and identification directly from single-channel recordings.

We present a method for analysis of noisy sampled data from a single-channel patch clamp which bypasses restoration of an idealized quantal signal. We show that, even in the absence of a specific model, the conductance levels and mean dwell times within those levels can be estimated. Estimation of the rate constants of a hypothesized kinetic scheme is more difficult. We present examples in which the rate constants can be effectively estimated and examples in which they cannot.

Bayesian restoration of single-channel patch clamp recordings.

The technique of patch clamp recording makes possible the measurement of current flowing through a single ion channel in a cell membrane. Examination of such recordings suggests that the current is quantal in nature, alternating in a seemingly random manner between "on" and "off," but the recordings are corrupted by noise from a variety of sources. In this paper we propose and illustrate methods for restoring the underlying quantal signal from such noisy measurements. The methods use a Markov chain prior distribution for the underlying quantal process and base the restoration on the resulting posterior distribution.

On the superposition of currents from ion channels.

We derive a number of statistical properties of the superposition of several independent channels contributing to a patch-clamp recording. Failure of these properties indicates dependence of the channels and may suggest the nature of interactions. We show how properties such as dwell-time distributions of the individual channels may be determined from those of the superposition in the case that the channels are independent.

Protonation of interacting residues in a protein by a Monte Carlo method: application to lysozyme and the photosynthetic reaction center of Rhodobacter sphaeroides.

We used Monte Carlo methods to treat statistical problem of electrostatic interactions among many titrating amino acids and applied these methods to lysozyme and the photosynthetic reaction center of Rhodobacter sphaeroides, including all titrating sites. We computed the average protonation of residues as a function of pH from an equilibrium distribution of states generated by random sampling. Electrostatic energies were calculated from a finite difference solution to the linearized Poisson-Boltzmann equation using the coordinates from solved protein structures. For most calculations we used the Metropolis algorithm to sample protonation states; for strongly coupled sites, we substantially reduced sampling errors by using a modified algorithm that allows multiple site transitions. The Monte Carlo method agreed with calculations for a small test system, lysozyme, for which the complete partition function was calculated. We also calculated the pH dependence of the free energy change associated with electron transfer from the primary to the secondary quinone in the photosynthetic reaction center. The shape of the resulting curve agreed fairly well with experiment, but the proton uptake from which the free energy was calculated agreed only to within a factor of two with the observed values. We believe that this discrepancy resulted from errors in the individual electrostatic energy calculations rather than from errors in the Monte Carlo sampling.

Kinetic analysis of channel gating. Application to the cholinergic receptor channel and the chloride channel from Torpedo californica.

Identification of the minimum number of ways in which open and closed states communicate is a crucial step in defining the gating kinetics of multistate channels. We used certain correlation functions to extract information about the pathways connecting the open and closed states of the cation channel of the purified nicotinic acetylcholine receptor and of the chloride channel of Torpedo californica electroplax membranes. Single channel currents were recorded from planar lipid bilayers containing the membrane channel proteins under investigation. The correlation functions are conveniently computed from single channel current records and yield information on E, the minimum number of entry/exit states into the open or closed aggregates. E gives a lower limit on the numbers of transition pathways between open and closed states. For the acetylcholine receptor, the autocorrelation analysis shows that there are at least two entry/exit states through which the open and closed aggregates communicate. The chloride channel fluctuates between three conductance substates, here indentified as C, M, and H for closed, intermediate, and high conductance, respectively. Correlation analysis shows that E is greater than or equal to 2 for the M aggregate, indicating that there are at least two distinct entry/exit states in the M aggregate. In contrast, there is no evidence for the existence of more than one entry/exit state in the C or H aggregates. Thus, these correlation functions provide a simple and general strategy to extract information on channel gating kinetics.

The electronic structure of Fe2+ in reaction centers from Rhodopseudomonas sphaeroides. III. EPR measurements of the reduced acceptor complex.

Electron paramagnetic resonance (EPR) spectra of the reduced quinone-iron acceptor complex in reaction centers were measured in a variety of environments and compared with spectra calculated from a theoretical model. Spectra were obtained at microwave frequencies of 1, 9, and 35 GHz and at temperatures from 1.4 to 30 K. The spectra are characterized by a broad absorption peak centered at g = 1.8 with wings extending from g approximately equal to 5 to g less than 0.8. The peak is split with the low-field component increasing in amplitude with temperature. The theoretical model is based on a spin Hamiltonian, in which the reduced quinone, Q-, interacts magnetically with Fe2+. In this model the ground manifold of the interacting Q-Fe2+ system has two lowest doublets that are separated by approximately 3 K. Both perturbation analyses and exact numerical calculations were used to show how the observed spectrum arises from these two doublets. The following spin Hamiltonian parameters optimized the agreement between simulated and observed spectra: the electronic g tensor gFe, x = 2.16, gFe, y = 2.27, gFez = 2.04, the crystal field parameters D = 7.60 K and E/D = 0.25, and the antiferromagnetic magnetic interaction tensor, Jx = -0.13 K, Jy = -0.58 K, Jz = -0.58 K. The model accounts well for the g value (1.8) of the broad peak, the observed splitting of the peak, the high and low g value wings, and the observed temperature dependence of the shape of the spectra. The structural implications of the value of the magnetic interaction, J, and the influence of the environment on the spin Hamiltonian parameters are discussed. The similarity of spectra and relaxation times observed from the primary and secondary acceptor complexes Q-AFe2+ and Fe2+Q-B leads to the conclusion that the Fe2+ is approximately equidistant from QA and QB.

Channel properties of the purified acetylcholine receptor from Torpedo californica reconstituted in planar lipid bilayer membranes.

The electrophysiological properties of the cation channel of the purified nicotinic acetylcholine receptor (AChR) reconstituted in planar lipid bilayers were characterized. Single-channel currents were activated by acetylcholine, carbamylcholine and suberyldicholine. The single channel conductance (28 pS in 0.3 M NaCl) was ohmic and independent of the agonist. Single channel currents increased with Na+ concentration to a maximum conductance of 95 pS and showed a half-saturation point of 395 mM. The apparent ion selectivity sequence, derived from single-channel current recordings, is: NH+4 greater than Cs+ greater than Rb+ greater than or equal to Na+ Cl-, F-, SO2-(4). The distribution of channel open times was fit by a sum of two exponentials, reflecting the existence of at least two distinct open states. The time constants depend on the choice of agonist, being consistently longer for suberyldicholine than for carbamylcholine. Similar channel properties were recorded in bilayers formed from monolayers at the tip of patch pipets . Single-channel currents occur in paroxysms of channel activity followed by quiescent periods. This pattern is more pronounced as the agonist concentration increases, and is reflected in histograms of channel-opening frequencies. Computer simulations with a three-state model, consisting of two closed (unliganded and liganded) and one open state, do not resemble the recorded pattern of channel activity, especially at high agonist concentration. Inclusion of a desensitized liganded state reproduces the qualitative features of channel recordings. The occurrence of paroxysms of channel activity thus seems to result from the transit of AChR through its active conformation, from which it can open several times before desensitizing.

Incorporation of acetylcholine receptors into liposomes. Vesicle structure and acetylcholine receptor function.

Functionally intact acetylcholine receptors can be solubilized from electric organ membranes of Torpedo californica and incorporated into liposomes by the cholate dialysis technique. Freezing and thawing of the reconstituted preparation appears to seal a population of initially leaky vesicles and leads to vesicle fusion. Inclusion of supplementary cholesterol at an optimal concentration of 20% (w/w) greatly enhances vesicle fusion during the freeze-thaw cycle. Size analysis by electron microscopy of negatively stained preparations indicates that fusion is accompanied by shifts in size and volume distributions of the vesicle population. Liposomes formed in the absence of acetylcholine receptors are distributed over a substantially smaller size range than liposomes containing receptors. Acetylcholine receptors appear in those liposomes as dimers of 80 A doughnut-shaped particles. Freeze-fracture replicas of reconstituted preparations reveal the presence of large vesicles containing particles which correspond in size to acetylcholine receptors and smaller liposomes devoid of particles. The distribution of particles in the reconstituted membranes is sparse compared to their dense packing in native electric organ membranes. The activation and desensitization of reconstituted acetylcholine receptors mediated by acetylcholine or carbamylcholine is dose dependent. The reconstituted receptors distinguish between these agonists in terms of binding affinity in a way similar to receptors in the native membrane. Correlation of the fractional occupancy of ligand binding sites by cobratoxin with inhibition of receptor function is used to demonstrate that in the reconstituted system the doubly liganded acetylcholine receptor prevails in controlling channel gating. The potential experimental advantages as well as limitations of this reconstituted system are discussed.

The electronic structure of Fe2+ in reaction centers from Rhodopseudomonas sphaeroides. I. Static magnetization measurements.

We have measured the static magnetization of unreduced and reduced reaction centers that vary in their quinone content. Measurements were performed in the temperature range 0.7 degrees K less than T less than 200 degrees K and magnetic fields of up to 10 kG. The electronic g-value, crystal field parameters D, E, and the exchange interaction, J, between the quinone spin and Fe2+ were determined using the spin Hamiltonian formalism. The effective moment mu eff/Fe2+ of both reduced and unreduced samples were determined to be 5.35 +/- 0.15 Bohr magnetons. This shows, in agreement with previous findings, that Fe2+ does not change its valence state when the reaction centers are reduced. Typical values of D congruent to +5 cm-1 and E/D congruent to 0.27 are consistent with Fe being in an octahedral environment with rhombic distortion. The values of D and E were approximately the same for reaction centers having one and two quinones. These findings imply that quinone is most likely not a ligand of Fe. The Fe2+ and the spin on the quinone in reduced reaction centers were found to be coupled with an exchange interaction 0 less than /J/ less than 1 cm-1. The validity of the spin Hamiltonian was checked by using an orbital Hamiltonian to calculate energy levels of the 25 states of the S = 2, L = 2 manifold and comparing the magnetization of the lowest five states with those obtained from the spin Hamiltonian. Using the orbital Hamiltonian, we calculated the position of the first excited quintet state to be 340 cm-1 above the ground state quintet. This is in good agreement with the temperature dependence of the quadrupole splitting as determined by Mossbauer spectroscopy.