An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Single-molecule spectroscopy of the β2 adrenergic receptor: Observation of conformational substates in a membrane protein

Single-molecule studies of the conformations of the intact β2 adrenergic receptor were performed in solution. Photon bursts from the fluorescently tagged adrenergic receptor in a micelle were recorded. A photon-burst algorithm and a Poisson time filter were implemented to characterize single molecules diffusing across the probe volume of a confocal microscope. The effects of molecular diffusion and photon number fluctuations were deconvoluted by assuming that Poisson distributions characterize the molecular occupation and photon numbers. Photon-burst size histograms were constructed, from which the source intensity distributions were extracted. Different conformations of the β2 adrenergic receptor cause quenching of the bound fluorophore to different extents and hence produce different photon-burst sizes. An analysis of the photon-burst histograms shows that there are at least two distinct substates for the native adrenergic membrane receptor. This behavior is in contrast to one peak observed for the dye molecule, rhodamine 6G. We test the reliability and robustness of the substate number determination by investigating the application of different binning criteria. Conformational changes associated with agonist binding result in a marked change in the distribution of photon-burst sizes. These studies provide insight into the conformational heterogeneity of G protein-coupled receptors in the presence and absence of a bound agonist.

Imaging ROMK1 inwardly rectifying ATP-sensitive K+ channel protein using atomic force microscopy.

The inwardly rectifying K+ channel ROMK1 has been implicated as being significant in K+ secretion in the distal nephron. ROMK1 has been shown by immunocytochemistry to be expressed in relevant nephron segments. The development of the atomic force microscope has made possible the production of high resolution images of small particles, including a variety of biological macromolecules. Recently, a fusion protein of glutathione S-transferase (GST) and ROMK1 (ROMK1-GST) has been used to produce a polyclonal antibody for immunolocalization of ROMK1. We have used atomic force microscopy to examine ROMK1-GST and the native ROMK1 polypeptide cleaved from GST. Imaging was conducted with the proteins in physiological solutions attached to mica. ROMK1-GST appears in images as a particle composed of two units of similar size. Analyses of images indicate that the two units have volumes of approximately 118 nm3, which is close to the theoretical volume of a globular protein of approximately 65 kDa (the molecular mass of ROMK1-GST). Native GST exists as a dimer, and the images obtained here are consistent with the ROMK1-GST fusion protein's existence as a heterodimer. In experiments on ROMK1 in aqueous solution, single molecules appear to aggregate, but contact to the mica was maintained. Addition of ATP to the solution produced a change in height of the aggregates. This change (which was reversible) suggests that ATP induces a structural change in the ROMK1 protein. The data show that atomic force microscopy is a useful tool for examination of purified protein molecules under near-physiological conditions, and furthermore, that structural alterations in the proteins may be continuously investigated.

A single residue in the M2-M3 loop is a major determinant of coupling between binding and gating in neuronal nicotinic receptors.

Binding of agonists to nicotinic acetylcholine receptors generates a sequence of changes that activate a cation-selective conductance. By measuring electrophysiological responses in chimeric alpha7/alpha3 receptors expressed in Xenopus oocytes, we have showed the involvement of the M2-M3 loop in coupling agonist binding to the channel gate. An aspartate residue therein, Asp-266 in the alpha7 subunit, was identified by site-directed mutagenesis as crucial, since mutants at this position exhibited very poor functional responses to three different nicotinic agonists. We have extended this investigation to another neuronal nicotinic receptor (alpha3/beta4), and found that a homologous residue in the beta4 subunit, Asp-268, played a similar role in coupling. These findings are consistent with a hypothesis that the aspartate residue in the M2-M3 loop, which is conserved in all homomer-forming alpha-type subunits and all neuronal beta-type subunits that combine to form functional receptors, is a major determinant of information transmission from binding site to channel gate in all neuronal nicotinic receptors.

Oscillations in KATP channel activity promote oscillations in cytoplasmic free Ca2+ concentration in the pancreatic beta cell.

Pancreatic beta cells exhibit oscillations in electrical activity, cytoplasmic free Ca2+ concentration ([Ca2+](i)), and insulin release upon glucose stimulation. The mechanism by which these oscillations are generated is not known. Here we demonstrate fluctuations in the activity of the ATP-dependent K+ channels (K(ATP) channels) in single beta cells subject to glucose stimulation or to stimulation with low concentrations of tolbutamide. During stimulation with glucose or low concentrations of tolbutamide, K(ATP) channel activity decreased and action potentials ensued. After 2-3 min, despite continuous stimulation, action potentials subsided and openings of K(ATP) channels could again be observed. Transient suppression of metabolism by azide in glucose-stimulated beta cells caused reversible termination of electrical activity, mimicking the spontaneous changes observed with continuous glucose stimulation. Thus, oscillations in K(ATP) channel activity during continuous glucose stimulation result in oscillations in electrical activity and [Ca2+](i).

An ATP-activated, ligand-gated ion channel on a cholinergic presynaptic nerve terminal.

ATP has recently been identified as a fast neurotransmitter in both the central and peripheral nervous systems. Several studies have suggested that ATP can also affect the release of classical neurotransmitters, including acetylcholine with which it is co-released. We have searched for ATP receptors on a cholinergic presynaptic nerve terminal using the calyx-type synapse of the chicken ciliary ganglion. ATP was pulsed onto the terminals under voltage clamp and induced a short latency cation current that exhibited inward rectification and marked desensitization. This current was not seen with adenosine but was mimicked by several sterically restricted ATP analogs and was blocked by suramin. ATP-activated single ion channels exhibited prominent flickering and had a conductance of approximately 17 pS. Our results demonstrate a ligand-gated P2X-like purinergic receptor on a cholinergic presynaptic nerve terminal.

A single amino acid in gamma-aminobutyric acid rho 1 receptors affects competitive and noncompetitive components of picrotoxin inhibition.

A class of bicuculline-insensitive gamma-aminobutyric acid (GABA) receptors, GABAC, has been identified in retina. Several lines of evidence indicate that GABAC receptors are formed partially or wholly of GABA rho subunits. These receptors generate a Cl- current in response to GABA but differ from GABAA receptors in a number of ways. Picrotoxin, widely accepted as a noncompetitive antagonist of GABAA receptors, displays competitive and noncompetitive antagonism of GABAC receptors in perch and bovine retina and GABA rho 1 receptors expressed in Xenopus oocytes. The aim of this study was to identify the molecular basis of the two components of picrotoxin inhibition of GABA rho 1 receptors. By using a domain-swapping and mutagenesis strategy, a difference in picrotoxin sensitivity between rho 1 and rho 2 receptors was localized to a single amino acid in the putative second transmembrane domain. Substitution of this amino acid with residues found in the analogous position in highly picrotoxin-sensitive glycine alpha and GABAA subunits increased the sensitivity of rho 1 mutants 10- to 500-fold. Importantly, the competitive component of picrotoxin inhibition of the rho 1 mutant receptors was almost eliminated. These findings demonstrate that an amino acid in the putative channel domain of GABA rho 1 receptors influences picrotoxin sensitivity and mediates agonist binding by an allosteric mechanism.

Cloning of a gamma-aminobutyric acid type C receptor subunit in rat retina with a methionine residue critical for picrotoxinin channel block.

Ionotropic receptors for gamma-aminobutyric acid (GABA) are important to inhibitory neurotransmission in the mammalian retina, mediating GABAA and GABAC responses. In many species, these responses are blocked by the convulsant picrotoxinin (PTX), although the mechanism of block is not fully understood. In contrast, GABAC responses in the rat retina are extremely resistant to PTX. We hypothesized that this difference could be explained by molecular characterization of the receptors underlying the GABAC response. Here we report the cloning of two rat GABA receptor subunits, designated r rho 1 and r rho 2 after their previously identified human homologues. When coexpressed in Xenopus oocytes, r rho 1/r rho 2 heteromeric receptors mimicked PTX-resistant GABAC responses of the rat retina. PTX resistance is apparently conferred in native heteromeric receptors by r rho 2 subunits since homomeric r rho 1 receptors were sensitive to PTX; r rho 2 subunits alone were unable to form functional homomeric receptors. Site-directed mutagenesis confirmed that a single amino acid residue in the second membrane-spanning region (a methionine in r rho 2 in place of a threonine in r rho 1) is the predominant determinant of PTX resistance in the rat receptor. This study reveals not only the molecular mechanism underlying PTX blockade of GABA receptors but also the heteromeric nature of native receptors in the rat retina that underlie the PTX-resistant GABAC response.

Molecular determinants of drug access to the receptor site for antiarrhythmic drugs in the cardiac Na+ channel.

The clinical efficacy of local anesthetic and antiarrhythmic drugs is due to their voltage- and frequency-dependent block of Na+ channels. Quaternary local anesthetic analogs such as QX-314, which are permanently charged and membrane-impermeant, effectively block cardiac Na+ channels when applied from either side of the membrane but block neuronal Na+ channels only from the intracellular side. This difference in extracellular access to QX-314 is retained when rat brain rIIA Na+ channel alpha subunits and rat heart rH1 Na+ channel alpha subunits are expressed transiently in tsA-201 cells. Amino acid residues in transmembrane segment S6 of homologous domain IV (IVS6) of Na+ channel alpha subunits have important effects on block by local anesthetic drugs. Although five amino acid residues in IVS6 differ between brain rIIA and cardiac rH1, exchange of these amino acid residues by site-directed mutagenesis showed that only conversion of Thr-1755 in rH1 to Val as in rIIA was sufficient to reduce the rate and extent of block by extracellular QX-314 and slow the escape of drug from closed channels after use-dependent block. Tetrodotoxin also reduced the rate of block by extracellular QX-314 and slowed escape of bound QX-314 via the extracellular pathway in rH1, indicating that QX-314 must move through the pore to escape. QX-314 binding was inhibited by mutation of Phe-1762 in the local anesthetic receptor site of rH1 to Ala whether the drug was applied extracellularly or intracellularly. Thus, QX-314 binds to a single site in the rH1 Na+ channel alpha subunit that contains Phe-1762, whether it is applied from the extracellular or intracellular side of the membrane. Access to that site from the extracellular side of the pore is determined by the amino acid at position 1755 in the rH1 cardiac Na+ channel.

Multiple calcium channel transcripts in rat osteosarcoma cells: selective activation of alpha 1D isoform by parathyroid hormone.

Osteoblasts express calcium channels that are thought to be involved in the transduction of extracellular signals regulating bone metabolism. The molecular identity of the pore-forming subunit (alpha 1) of L-type calcium channel(s) was determined in rat osteosarcoma UMR-106 cells, which express an osteoblast phenotype. A homology-based reverse transcriptase-polymerase chain reaction cloning strategy was employed that used primers spanning the fourth domain. Three types of cDNAs were isolated, corresponding to the alpha 1S (skeletal), alpha 1C (cardiac), and alpha 1D (neuroendocrine) isoforms. In the transmembrane segment IVS3 and the extracellular loop formed by the IVS3-S4 linker, a single pattern of mRNA splicing was found that occurs in all three types of calcium channel transcripts. Northern blot analysis revealed an 8.6-kb mRNA that hybridized to the alpha 1C probe and 4.8- and 11.7-kb mRNAs that hybridized to the alpha 1S and alpha 1D probes. Antisense oligonucleotides directed to the calcium channel alpha 1D transcript, but not those directed to alpha 1S or alpha 1C transcripts, inhibited the rise of intracellular calcium induced by parathyroid hormone. However, alpha 1D antisense oligonucleotides had no effect on the accumulation of cAMP induced by parathyroid hormone. When L-type calcium channels were activated with Bay K 8644, antisense oligonucleotides to each of the three isoforms partially inhibited the rise of intracellular calcium. The present results provide evidence for the expression of three distinct calcium channel alpha 1-subunit isoforms in an osteoblast-like cell line. We conclude that the alpha 1D isoform is selectively activated by parathyroid hormone.

Features of MotA proton channel structure revealed by tryptophan-scanning mutagenesis.

The MotA protein of Escherichia coli is a component of the flagellar motors that functions in transmembrane proton conduction. Here, we report several features of MotA structure revealed by use of a mutagenesis-based approach. Single tryptophan residues were introduced at many positions within the four hydrophobic segments of MotA, and the effects on function were measured. Function was disrupted according to a periodic pattern that implies that the membrane-spanning segments are alpha-helices and that identifies the lipid-facing parts of each helix. The results support a hypothesis for MotA structure and mechanism in which water molecules form most of the proton-conducting pathway. The success of this approach in studying MotA suggests that it could be useful in structure-function studies of other integral membrane proteins.

Neural-network-based classification of cognitively normal, demented, Alzheimer disease and vascular dementia from single photon emission with computed tomography image data from brain.

Single photon emission with computed tomography (SPECT) hexamethylphenylethyleneamineoxime technetium-99 images were analyzed by an optimal interpolative neural network (OINN) algorithm to determine whether the network could discriminate among clinically diagnosed groups of elderly normal, Alzheimer disease (AD), and vascular dementia (VD) subjects. After initial image preprocessing and registration, image features were obtained that were representative of the mean regional tissue uptake. These features were extracted from a given image by averaging the intensities over various regions defined by suitable masks. After training, the network classified independent trials of patients whose clinical diagnoses conformed to published criteria for probable AD or probable/possible VD. For the SPECT data used in the current tests, the OINN agreement was 80 and 86% for probable AD and probable/possible VD, respectively. These results suggest that artificial neural network methods offer potential in diagnoses from brain images and possibly in other areas of scientific research where complex patterns of data may have scientifically meaningful groupings that are not easily identifiable by the researcher.

Optical probing of spin fluctuations of a single paramagnetic Mn atom in a semiconductor quantum dot

We analyzed the photoluminescence intermittency generated by a single paramagnetic spin localized in an individual semiconductor quantum dot. The statistics of the photons emitted by the quantum dot reflect the quantum fluctuations of the localized spin interacting with the injected carriers. Photon correlation measurements, which are reported here, reveal unique signatures of these fluctuations. A phenomenological model is proposed to quantitatively describe these observations, allowing a measurement of the spin dynamics of an individual magnetic atom at zero magnetic field. These results demonstrate the existence of an efficient spin-relaxation channel arising from a spin exchange with individual carriers surrounding the quantum dot. A theoretical description of a spin-flip mechanism involving spin exchange with surrounding carriers gives relaxation times in good agreement with the measured dynamics.

Logic hybrid simulation-optimization algorithm for distillation design

In this paper, we propose a novel algorithm for the rigorous design of distillation columns that integrates a process simulator in a generalized disjunctive programming formulation. The optimal distillation column, or column sequence, is obtained by selecting, for each column section, among a set of column sections with different number of theoretical trays. The selection of thermodynamic models, properties estimation etc., are all in the simulation environment. All the numerical issues related to the convergence of distillation columns (or column sections) are also maintained in the simulation environment. The model is formulated as a Generalized Disjunctive Programming (GDP) problem and solved using the logic based outer approximation algorithm without MINLP reformulation. Some examples involving from a single column to thermally coupled sequence or extractive distillation shows the performance of the new algorithm.