Ion Selectivity in the KcsA Potassium Channel from the Perspective of the Ion Binding Site

To understand the thermodynamic exclusion of Na⁺ relative to K⁺ from the S₂ site of the selectivity filter, the distribution P_X(ɛ) (X = K⁺ or Na⁺) of the binding energy (ɛ) of the ion with the channel is analyzed using the potential distribution theorem. By expressing the excess chemical potential of the ion as a sum of mean-field 〈ɛ〉 and fluctuation μ^ex_flux,X contributions, we find that selectivity arises from a higher value of μflux,Na+ex relative to μflux,K+ex. To understand the role of site-site interactions on μ^ex_flux,X, we decompose P_X(ɛ) into n-dependent distributions, where n is the number of ion-coordinating ligands within a distance λ from the ion. For λ comparable to typical ion-oxygen bond distances, investigations building on this multistate model reveal an inverse correlation between favorable ion-site and site-site interactions: the ion-coordination states that most influence the thermodynamics of the ion are also those for which the binding site is energetically less strained and vice versa. This correlation motivates understanding entropic effects in ion binding to the site and leads to the finding that μ^ex_flux,X is directly proportional to the average site-site interaction energy, a quantity that is sensitive to the chemical type of the ligand coordinating the ion. Increasing the coordination number around Na⁺ only partially accounts for the observed magnitude of selectivity; acknowledging the chemical type of the ion-coordinating ligand is essential.     

Engineered Bi-Histidine Metal Chelation Sites Map the Structure of the Mechanical Unfolding Transition State of an Elastomeric Protein Domain GB1

Determining the structure of the transition state is critical for elucidating the mechanism behind how proteins fold and unfold. Due to its high free energy, however, the transition state generally cannot be trapped and studied directly using traditional structural biology methods. Thus, characterizing the structure of the transition state that occurs as proteins fold and unfold remains a major challenge. Here, we report a novel (to our knowledge) method that uses engineered bi-histidine (bi-His) metal-binding sites to directly map the structure of the mechanical unfolding transition state of proteins. This method is adapted from the traditional ψ-value analysis, which uses engineered bi-His metal chelation sites to probe chemical (un)folding transition-state structure. The ϕ^M2+_U-value is defined as ΔΔG_‡-N/ΔΔG_U-N, which is the energetic effects of metal chelation by the bi-His site on the unfolding energy barrier (ΔG_‡-N) relative to its thermodynamic stability (ΔG_U-N) and can be used to obtain information about the transition state in the mutational site. As a proof of principle, we used the small protein GB1 as a model system and set out to map its mechanical unfolding transition-state structure. Using single-molecule atomic force microscopy and spectrofluorimetry, we directly quantified the effect of divalent metal ion binding on the mechanical unfolding free energy and thermodynamic stability of GB1, which allowed us to quantify ϕ^M2+_U-values for different sites in GB1. Our results enabled us to map the structure of the mechanical unfolding transition state of GB1. Within GB1’s mechanical unfolding transition state, the interface between force-bearing β-strands 1 and 4 is largely disrupted, and the first β-hairpin is partially disordered while the second β-hairpin and the α-helix remain structured. Our results demonstrate the unique application of ψ-value analysis in elucidating the structure of the transition state that occurs during the mechanical unfolding process, offering a potentially powerful new method for investigating the design of novel elastomeric proteins.     

Prion formation and polyglutamine aggregation are controlled by two classes of genes

Prions are self-perpetuating aggregated proteins that are not limited to mammalian systems but also exist in lower eukaryotes including yeast. While much work has focused around chaperones involved in prion maintenance, including Hsp104, little is known about factors involved in the appearance of prions. De novo appearance of the [PSI ⁺] prion, which is the aggregated form of the Sup35 protein, is dramatically enhanced by transient overexpression of SUP35 in the presence of the prion form of the Rnq1 protein, [PIN ⁺]. When fused to GFP and overexpressed in [ps⁻] [PIN ⁺] cells, Sup35 forms fluorescent rings, and cells with these rings bud off [PSI ⁺] daughters. We investigated the effects of over 400 gene deletions on this de novo induction of [PSI ⁺]. Two classes of gene deletions were identified. Class I deletions (bug1Δ, bem1Δ, arf1Δ, and hog1Δ) reduced the efficiency of [PSI ⁺] induction, but formed rings normally. Class II deletions (las17Δ, vps5Δ, and sac6Δ) inhibited both [PSI⁺] induction and ring formation. Furthermore, class II deletions reduced, while class I deletions enhanced, toxicity associated with the expanded glutamine repeats of the huntingtin protein exon 1 that causes Huntington''s disease. This suggests that prion formation and polyglutamine aggregation involve a multi-phase process that can be inhibited at different steps.     

Thermodynamic analysis of F1‐ATPase rotary catalysis using high‐speed imaging

F₁-ATPase (F₁) is a rotary motor protein fueled by ATP hydrolysis. Although the mechanism for coupling rotation and catalysis has been well studied, the molecular details of individual reaction steps remain elusive. In this study, we performed high-speed imaging of F₁ rotation at various temperatures using the total internal reflection dark-field (TIRDF) illumination system, which allows resolution of the F₁ catalytic reaction into elementary reaction steps with a high temporal resolution of 72 µs. At a high concentration of ATP, F₁ rotation comprised distinct 80° and 40° substeps. The 80° substep, which exhibited significant temperature dependence, is triggered by the temperature-sensitive reaction, whereas the 40° substep is triggered by ATP hydrolysis and the release of inorganic phosphate (P_i). Then, we conducted Arrhenius analysis of the reaction rates to obtain the thermodynamic parameters for individual reaction steps, that is, ATP binding, ATP hydrolysis, P_i release, and TS reaction. Although all reaction steps exhibited similar activation free energy values, ΔG^‡ = 53–56 kJ mol⁻¹, the contributions of the enthalpy (ΔH^‡), and entropy (ΔS^‡) terms were significantly different; the reaction steps that induce tight subunit packing, for example, ATP binding and TS reaction, showed high positive values of both ΔH^‡ and ΔS^‡. The results may reflect modulation of the excluded volume as a function of subunit packing tightness at individual reaction steps, leading to a gain or loss in water entropy.     

Long-pore Electrostatics in Inward-rectifier Potassium Channels

Inward-rectifier potassium (Kir) channels differ from the canonical K⁺ channel structure in that they possess a long extended pore (~85 Å) for ion conduction that reaches deeply into the cytoplasm. This unique structural feature is presumably involved in regulating functional properties specific to Kir channels, such as conductance, rectification block, and ligand-dependent gating. To elucidate the underpinnings of these functional roles, we examine the electrostatics of an ion along this extended pore. Homology models are constructed based on the open-state model of KirBac1.1 for four mammalian Kir channels: Kir1.1/ROMK, Kir2.1/IRK, Kir3.1/GIRK, and Kir6.2/KATP. By solving the Poisson-Boltzmann equation, the electrostatic free energy of a K⁺ ion is determined along each pore, revealing that mammalian Kir channels provide a favorable environment for cations and suggesting the existence of high-density regions in the cytoplasmic domain and cavity. The contribution from the reaction field (the self-energy arising from the dielectric polarization induced by the ion's charge in the complex geometry of the pore) is unfavorable inside the long pore. However, this is well compensated by the electrostatic interaction with the static field arising from the protein charges and shielded by the dielectric surrounding. Decomposition of the static field provides a list of residues that display remarkable correspondence with existing mutagenesis data identifying amino acids that affect conduction and rectification. Many of these residues demonstrate interactions with the ion over long distances, up to 40 Å, suggesting that mutations potentially affect ion or blocker energetics over the entire pore. These results provide a foundation for understanding ion interactions in Kir channels and extend to the study of ion permeation, block, and gating in long, cation-specific pores.     

On conduction in a bacterial sodium channel

Voltage-gated Na⁺-channels are transmembrane proteins that are responsible for the fast depolarizing phase of the action potential in nerve and muscular cells. Selective permeability of Na⁺ over Ca²⁺ or K⁺ ions is essential for the biological function of Na⁺-channels. After the emergence of the first high-resolution structure of a Na⁺-channel, an anionic coordination site was proposed to confer Na⁺ selectivity through partial dehydration of Na⁺ via its direct interaction with conserved glutamate side chains. By combining molecular dynamics simulations and free-energy calculations, a low-energy permeation pathway for Na⁺ ion translocation through the selectivity filter of the recently determined crystal structure of a prokaryotic sodium channel from Arcobacter butzleri is characterised. The picture that emerges is that of a pore preferentially occupied by two ions, which can switch between different configurations by crossing low free-energy barriers. In contrast to K⁺-channels, the movements of the ions appear to be weakly coupled in Na⁺-channels. When the free-energy maps for Na⁺ and K⁺ ions are compared, a selective site is characterised in the narrowest region of the filter, where a hydrated Na⁺ ion, and not a hydrated K⁺ ion, is energetically stable.     

Macoilin, a conserved nervous system-specific ER membrane protein that regulates neuronal excitability

Genome sequence comparisons have highlighted many novel gene families that are conserved across animal phyla but whose biological function is unknown. Here, we functionally characterize a member of one such family, the macoilins. Macoilins are characterized by several highly conserved predicted transmembrane domains towards the N-terminus and by coiled-coil regions C-terminally. They are found throughout Eumetazoa but not in other organisms. Mutants for the single Caenorhabditis elegans macoilin, maco-1, exhibit a constellation of behavioral phenotypes, including defects in aggregation, O₂ responses, and swimming. MACO-1 protein is expressed broadly and specifically in the nervous system and localizes to the rough endoplasmic reticulum; it is excluded from dendrites and axons. Apart from subtle synapse defects, nervous system development appears wild-type in maco-1 mutants. However, maco-1 animals are resistant to the cholinesterase inhibitor aldicarb and sensitive to levamisole, suggesting pre-synaptic defects. Using in vivo imaging, we show that macoilin is required to evoke Ca²⁺ transients, at least in some neurons: in maco-1 mutants the O₂-sensing neuron PQR is unable to generate a Ca²⁺ response to a rise in O₂. By genetically disrupting neurotransmission, we show that pre-synaptic input is not necessary for PQR to respond to O₂, indicating that the response is mediated by cell-intrinsic sensory transduction and amplification. Disrupting the sodium leak channels NCA-1/NCA-2, or the N-,P/Q,R-type voltage-gated Ca²⁺ channels, also fails to disrupt Ca²⁺ responses in the PQR cell body to O₂ stimuli. By contrast, mutations in egl-19, which encodes the only Caenorhabditis elegans L-type voltage-gated Ca²⁺ channel α1 subunit, recapitulate the Ca²⁺ response defect we see in maco-1 mutants, although we do not see defects in localization of EGL-19. Together, our data suggest that macoilin acts in the ER to regulate assembly or traffic of ion channels or ion channel regulators.     

THE COUPLED REDOX POTENTIAL OF THE LACTATE-ENZYME-PYRUVATE SYSTEM

1. The term "coupled redox potential" is defined. 2. The system lactic ion See PDF for Equation pyruvic ion + 2H⁺ + 2e is shown to be reversible (when the enzyme is lactic acid dehydrogenase) and its coupled redox potential between pH 5.2 and 7.2 at 32°C. is: See PDF for Equation 3. The free energy of the reaction: lactic ion (1m) → pyruvic ion (1m) = -ΔF = –14,572. 4. The standard free energy of formation (ΔF ₂₉₈) of pyruvic acid (l) is estimated at –108,127. This is merely an approximation as some necessary data are lacking. 5. The importance of coupled redox potentials as a factor in the regulation of the equilibrium of metabolites is indicated.     

Pseudomonas aeruginosa elastase provides an escape from phagocytosis by degrading the pulmonary surfactant protein-A

Pseudomonas aeruginosa is an opportunistic pathogen that causes both acute pneumonitis in immunocompromised patients and chronic lung infections in individuals with cystic fibrosis and other bronchiectasis. Over 75% of clinical isolates of P. aeruginosa secrete elastase B (LasB), an elastolytic metalloproteinase that is encoded by the lasB gene. Previously, in vitro studies have demonstrated that LasB degrades a number of components in both the innate and adaptive immune systems. These include surfactant proteins, antibacterial peptides, cytokines, chemokines and immunoglobulins. However, the contribution of LasB to lung infection by P. aeruginosa and to inactivation of pulmonary innate immunity in vivo needs more clarification. In this study, we examined the mechanisms underlying enhanced clearance of the ΔlasB mutant in mouse lungs. The ΔlasB mutant was attenuated in virulence when compared to the wild-type strain PAO1 during lung infection in SP-A^+/+ mice. However, the ΔlasB mutant was as virulent as PAO1 in the lungs of SP-A^-/- mice. Detailed analysis showed that the ΔlasB mutant was more susceptible to SP-A-mediated opsonization but not membrane permeabilization. In vitro and in vivo phagocytosis experiments revealed that SP-A augmented the phagocytosis of ΔlasB mutant bacteria more efficiently than the isogenic wild-type PAO1. The ΔlasB mutant was found to have a severely reduced ability to degrade SP-A, consequently making it unable to evade opsonization by the collectin during phagocytosis. These results suggest that P. aeruginosa LasB protects against SP-A-mediated opsonization by degrading the collectin.     

Collapse Dynamics of Single Proteins Extended by Force

Single-molecule force spectroscopy has opened up new approaches to the study of protein dynamics. For example, an extended protein folding after an abrupt quench in the pulling force was shown to follow variable collapse trajectories marked by well-defined stages that departed from the expected two-state folding behavior that is commonly observed in bulk. Here, we explain these observations by developing a simple approach that models the free energy of a mechanically extended protein as a combination of an entropic elasticity term and a short-range potential representing enthalpic hydrophobic interactions. The resulting free energy of the molecule shows a force-dependent energy barrier of magnitude, ΔE = ɛ(F − F_c)^3/2, separating the enthalpic and entropic minima that vanishes at a critical force F_c. By solving the Langevin equation under conditions of a force quench, we generate folding trajectories corresponding to the diffusional collapse of an extended polypeptide. The predicted trajectories reproduce the different stages of collapse, as well as the magnitude and time course of the collapse trajectories observed experimentally in ubiquitin and I27 protein monomers. Our observations validate the force-clamp technique as a powerful approach to determining the free-energy landscape of proteins collapsing and folding from extended states.     

A comparison of galvanic skin conductance and skin wettedness as indicators of thermal discomfort during moderate and high metabolic rates

The relationship between local thermal comfort, local skin wettedness (w_local) and local galvanic skin conductance (GSC) in four body segments during two different exercise intensities was compared in 10 males. In a balanced order, participants walked at 35% VO_2max for 45 min (WALK) (29.0±1.9°C, 29.8±3.6% RH, no wind) in one test and in a separate test ran at 70% VO_2max for 45 min (RUN) (26.2±2.1°C, 31.1±7.0% RH, no wind). During both tests, participants wore a loose fitting 100% polyester long sleeve top and trouser ensemble with a low resistance to heat and vapour transfer (total thermal resistance of 0.154 m² K W⁻¹ and total water vapour resistance of 35.9 m² Pa W⁻¹). w_local, change from baseline in GSC (ΔGSC) and local thermal comfort were recorded every 5 min. The results suggest that both w_local and ΔGSC are strong predictors of thermal comfort during the WALK when sweat production is low and thermal discomfort minimal (r²>0.78 and r²>0.71, respectively). Interestingly, during the RUN w_local plateaued at ~0.6 to 0.8 due to the high sweat production, whilst ΔGSC gradually increased throughout the experiment. ΔGSC had a similar relationship with thermal comfort to w_local during the RUN (r²>0.95 and r²>0.94, respectively). Despite the strength of these relationships, the ability of w_local to predict local thermal comfort accurately dramatically reduces in the exponential part of the curve. In a situation of uncompensated heat stress such as high metabolic rate in hot climate, where sweat production is high, ΔGSC shows to be a better predictor of local thermal comfort than w_local. The w_local data shows regional differences in the threshold which triggers local discomfort during the WALK than RUN; lower values are found for upper arms (0.22±0.03 and 0.28 ±0.22) and upper legs (0.22±0.11 and 0.22±0.10), higher values for upper back (0.30±0.12 and 0.36 ±0.10) and chest (0.27±0.10 and 0.39 ±0.32), respectively. However, no regional differences in the threshold of discomfort are found in the ?GSC data. Instead, the data suggests that the degree of discomfort experienced appears to be related to the amount of sweat within and around the skin (as indirectly measured by ΔGSC) at each body site.     

Ion channel inhibitors block caspase activation by mechanisms other than restoring intracellular potassium concentration

Ion fluxes at the plasma membrane have an important role in early stages of apoptosis. Accordingly, plasma membrane depolarization and gain of Na⁺ and loss of K⁺ are initial events in apoptosis. We have studied the effect of staurosporine (STS), a well-established apoptosis inducer, on the membrane potential of HeLa cells to determine the nature of STS-activated ion conductances and their role in the activation of different caspases. We observed that STS can activate tetraethylammonium (TEA⁺) and 4-aminopyridine-sensitive K⁺ channels and flufenamic-sensitive cation channels as an early response. The combination of these ion channel inhibitors significantly reduced cytochrome c (cyt c) release and activation of caspase-9, -3 and -8. STS also induced a large reduction in the intracellular [K⁺] that was not blocked by the ion channel inhibitors. Our data suggest that reduction in the [K⁺]_i is necessary but not sufficient and that ion channel inhibitors block activation of caspase-3 by two different mechanisms: the inhibitors of K⁺ channels by reducing cyt c release while flufenamic acid by a different, unrelated mechanism that does not involve cation channels at the plasma membrane. Our data also imply that these ion channels activated by STS are not responsible for the reduction in the [K⁺]_i associated with apoptosis.     

Anthrax lethal toxin suppresses murine cardiomyocyte contractile function and intracellular Ca2+ handling via a NADPH oxidase-dependent mechanism

Objectives

Anthrax infection is associated with devastating cardiovascular sequelae, suggesting unfavorable cardiovascular effects of toxins originated from Bacillus anthracis namely lethal and edema toxins. This study was designed to examine the direct effect of lethal toxins on cardiomyocyte contractile and intracellular Ca²⁺ properties.

Methods

Murine cardiomyocyte contractile function and intracellular Ca²⁺ handling were evaluated including peak shortening (PS), maximal velocity of shortening/ relengthening (± dL/dt), time-to-PS (TPS), time-to-90% relengthening (TR₉₀), intracellular Ca²⁺ rise measured as fura-2 fluorescent intensity (ΔFFI), and intracellular Ca²⁺ decay rate. Stress signaling and Ca²⁺ regulatory proteins were assessed using Western blot analysis.

Results

In vitro exposure to a lethal toxin (0.05 – 50 nM) elicited a concentration-dependent depression on cardiomyocyte contractile and intracellular Ca²⁺ properties (PS, ± dL/dt, ΔFFI), along with prolonged duration of contraction and intracellular Ca²⁺ decay, the effects of which were nullified by the NADPH oxidase inhibitor apocynin. The lethal toxin significantly enhanced superoxide production and cell death, which were reversed by apocynin. In vivo lethal toxin exposure exerted similar time-dependent cardiomyocyte mechanical and intracellular Ca²⁺ responses. Stress signaling cascades including MEK1/2, p38, ERK and JNK were unaffected by in vitro lethal toxins whereas they were significantly altered by in vivo lethal toxins. Ca²⁺ regulatory proteins SERCA2a and phospholamban were also differentially regulated by in vitro and in vivo lethal toxins. Autophagy was drastically triggered although ER stress was minimally affected following lethal toxin exposure.

Conclusions

Our findings indicate that lethal toxins directly compromised murine cardiomyocyte contractile function and intracellular Ca²⁺ through a NADPH oxidase-dependent mechanism.     

Efflux and Influx of Erythrocyte Water

Rabbit erythrocytes were washed in buffered NaCl solutions isotonic with rabbit serum (Δ^t -0.558°C.) and suspended in buffered NaCl solutions of tonicity equidistant from intracellular tonicity (Δ^t = -0.558°C. ± 0.112°C.) of varying pH and incubated at varying temperatures. After incubation, the freezing point depression (Δ^t) was measured on the supernatant. Change in the Δ^t measured change in the water content of the extracellular solutions—water being withdrawn by erythrocytes (W_I) from the hypotonic solutions and added (W_E) to the hypertonic solutions. W_E was always less than W_I and was inversely proportional to the pH in the range 6.5–8.0. W_E was significantly increased by lowering the temperature of the cell suspension to 4°C. W_I was increased by raising or lowering the pH or raising the temperature of the cell suspension. W_E x W_I ≠ k. W_E and W_I were affected differently by changes in pH and temperature. It was concluded that W_E and W_E were probably under different physicochemical control.     

Thermodynamic characterization of the multivalent binding of chartreusin to DNA

Characterization of the thermodynamics of DNA– drug interactions is a very useful part in rational drug design. Isothermal titration calorimetry (ITC), differential scanning calorimetry (DSC) and UV melting experiments have been used to analyze the multivalent (intercalation plus minor groove) binding of the antitumor antibiotic chartreusin to DNA. Using DNA UV melting studies in the presence of the ligand and the binding enthalpy determined by ITC, we determined that the binding constant for the interaction was 3.6 × 10⁵ M^–1 at 20°C, in a solution containing 18 mM Na⁺. The DNA–drug interaction was enthalpy driven, with a ΔH_b of –7.07 kcal/mol at 20°C. Binding enthalpies were determined by ITC in the 20–35°C range and used to calculate a binding-induced change in heat capacity (ΔCp) of –391 cal/mol K. We have obtained a detailed thermodynamic profile for the interaction of this multivalent drug, which makes possible a dissection of ΔG_obs into the component free energy terms. The hydrophobic transfer of the chartreusin chromophore from the solution to the DNA intercalating site is the main contributor to the free energy of binding.     

Calcium-release channels in paramecium. Genomic expansion, differential positioning and partial transcriptional elimination

The release of Ca²⁺ from internal stores is a major source of signal Ca²⁺ in almost all cell types. The internal Ca²⁺ pools are activated via two main families of intracellular Ca²⁺-release channels, the ryanodine and the inositol 1,4,5-trisphosphate (InsP₃) receptors. Among multicellular organisms these channel types are ubiquitous, whereas in most unicellular eukaryotes the identification of orthologs is impaired probably due to evolutionary sequence divergence. However, the ciliated protozoan Paramecium allowed us to prognosticate six groups, with a total of 34 genes, encoding proteins with characteristics typical of InsP₃ and ryanodine receptors by BLAST search of the Paramecium database. We here report that these Ca²⁺-release channels may display all or only some of the characteristics of canonical InsP₃ and ryanodine receptors. In all cases, prediction methods indicate the presence of six trans-membrane regions in the C-terminal domains, thus corresponding to canonical InsP₃ receptors, while a sequence homologous to the InsP₃-binding domain is present only in some types. Only two types have been analyzed in detail previously. We now show, by using antibodies and eventually by green fluorescent protein labeling, that the members of all six groups localize to distinct organelles known to participate in vesicle trafficking and, thus, may provide Ca²⁺ for local membrane-membrane interactions. Whole genome duplication can explain radiation within the six groups. Comparative and evolutionary evaluation suggests derivation from a common ancestor of canonical InsP₃ and ryanodine receptors. With one group we could ascertain, to our knowledge for the first time, aberrant splicing in one thoroughly analyzed Paramecium gene. This yields truncated forms and, thus, may indicate a way to pseudogene formation. No comparable analysis is available for any other, free-living or parasitic/pathogenic protozoan.     

A Steep Dependence of Inward-Rectifying Potassium Channels on Cytosolic Free Calcium Concentration Increase Evoked by Hyperpolarization in Guard Cells

Inactivation of inward-rectifying K⁺ channels (I_K,in) by a rise in cytosolic free [Ca²⁺] ([Ca²⁺]_i) is a key event leading to solute loss from guard cells and stomatal closure. However, [Ca²⁺]_i action on I_K,in has never been quantified, nor are its origins well understood. We used membrane voltage to manipulate [Ca²⁺]_i (A. Grabov and M.R. Blatt [1998] Proc Natl Acad Sci USA 95: 4778–4783) while recording I_K,in under a voltage clamp and [Ca²⁺]_i by Fura-2 fluorescence ratiophotometry. I_K,in inactivation correlated positively with [Ca²⁺]_i and indicated a K_i of 329 ± 31 nm with cooperative binding of four Ca²⁺ ions per channel. I_K,in was promoted by the Ca²⁺ channel antagonists Gd³⁺ and calcicludine, both of which suppressed the [Ca²⁺]_i rise, but the [Ca²⁺]_i rise was unaffected by the K⁺ channel blocker Cs⁺. We also found that ryanodine, an antagonist of intracellular Ca²⁺ channels that mediate Ca²⁺-induced Ca²⁺ release, blocked the [Ca²⁺]_i rise, and Mn²⁺ quenching of Fura-2 fluorescence showed that membrane hyperpolarization triggered divalent release from intracellular stores. These and additional results point to a high signal gain in [Ca²⁺]_i control of I_K,in and to roles for discrete Ca²⁺ flux pathways in feedback control of the K⁺ channels by membrane voltage.Ca²⁺ underlies many fundamental regulatory processes in plants, including adaptive responses to abiotic environmental stress (Knight et al., 1996; Russell et al., 1996; McAinsh et al., 1997) and programmed cell death evoked by pathogen attack (Low and Merida, 1996; Hammondkosack and Jones, 1997). Coordination of changes in [Ca²⁺]_i and its integration with downstream response elements are central in coupling stimulus input to cellular response in these processes.In stomatal guard cells, the best characterized higher-plant cell model, major downstream targets of [Ca²⁺]_i and their roles in stomatal function have been identified. Increasing [Ca²⁺]_i is known to inactivate I_K,in and to activate Cl⁻ channels, events that bias plasma membrane transport for net efflux of osmotically active solute and a loss of turgor, which drives stomatal closure (Blatt and Grabov, 1997). Furthermore, changes in [Ca²⁺]_i are associated with ABA, CO₂, and the growth hormone auxin (Blatt and Grabov, 1997; McAinsh et al., 1997). These [Ca²⁺]_i signals have been observed to oscillate (McAinsh et al., 1995; Webb et al., 1996), characteristics that may constitute “Ca²⁺ signatures” to encode specific downstream responses (Berridge, 1996). Yet, despite the evidence for [Ca²⁺]_i signaling in guard cells, surprisingly little detail is known about the link between [Ca²⁺]_i changes and ion channel activity at the plasma membrane or about the mechanisms mediating such [Ca²⁺]_i changes. To our knowledge, in no instance have the characteristics of ion channel regulation by Ca²⁺ been quantified directly in any higher-plant cell.We recently described the coupling of membrane voltage to [Ca²⁺]_i, demonstrating that hyperpolarization, whether under a voltage clamp or in the presence of low [K⁺]_o, evoked [Ca²⁺]_i increases in guard cells, and that the voltage threshold for [Ca²⁺]_i rise was profoundly altered by ABA (Grabov and Blatt, 1998). Our observations indicated a link to Ca²⁺ influx across the plasma membrane and raised questions about the efficacy of [Ca²⁺]_i in inactivating I_K,in and about the contributions of intracellular Ca²⁺ release to the [Ca²⁺]_i signal. We have used membrane voltage to experimentally manipulate [Ca²⁺]_i and report that I_K,in is strongly dependent on [Ca²⁺]_i, consistent with a cooperative binding of four Ca²⁺ ions to effect inactivation. Additional experiments indicate that voltage-evoked [Ca²⁺]_i increases depend both on Ca²⁺ influx and on release of Ca²⁺ from intracellular stores. These results underscore the role of [Ca²⁺]_i as a high-gain “switch” in the control of I_K,in, and implicate [Ca²⁺]_i in feedback control linking membrane voltage to the activity of the K⁺ channels.     

The power of hard-sphere models: explaining side-chain dihedral angle distributions of Thr and Val

The energy functions used to predict protein structures typically include both molecular-mechanics and knowledge-based terms. In contrast, our approach is to develop robust physics- and geometry-based methods. Here, we investigate to what extent simple hard-sphere models can be used to predict side-chain conformations. The distributions of the side-chain dihedral angle χ₁ of Val and Thr in proteins of known structure show distinctive features: Val side chains predominantly adopt χ₁ = 180°, whereas Thr side chains typically adopt χ₁ = 60° and 300° (i.e., χ₁ = ±60° or g− and g⁺ configurations). Several hypotheses have been proposed to explain these differences, including interresidue steric clashes and hydrogen-bonding interactions. In contrast, we show that the observed side-chain dihedral angle distributions for both Val and Thr can be explained using only local steric interactions in a dipeptide mimetic. Our results emphasize the power of simple physical approaches and their importance for future advances in protein engineering and design.