Search

There are currently no drugs known to rescue the function of Kv1.1 voltage-gated potassium channels carrying loss-of-function sequence variants underlying the inherited movement disorder, Episodic Ataxia 1 (EA1). The Kwakwaka'wakw First Nations of the Pacific Northwest Coast used Fucus gardneri (bladderwrack kelp), Physocarpus capitatus (Pacific ninebark) and Urtica dioica (common nettle) to treat locomotor ataxia. Here, we show that extracts of these plants enhance wild-type Kv1.1 current, especially at subthreshold potentials. Screening of their constituents revealed that gallic acid and tannic acid similarly augment wild-type Kv1.1 current, with submicromolar potency. Crucially, the extracts and their constituents also enhance activity of Kv1.1 channels containing EA1-linked sequence variants. Molecular dynamics simulations reveal that gallic acid augments Kv1.1 activity via a small-molecule binding site in the extracellular S1-S2 linker. Thus, traditional Native American ataxia treatments utilize a molecular mechanistic foundation that can inform small-molecule approaches to therapeutically correcting EA1 and potentially other Kv1.1-linked channelopathies.

The human voltage-gated proton channel Hv1 is a drug target for cancer, ischemic stroke, and neuroinflammation. It resides on the plasma membrane and endocytic compartments of a variety of cell types, where it mediates outward proton movement and regulates the activity of NOX enzymes. Its voltage-sensing domain (VSD) contains a gated and proton-selective conduction pathway, which can be blocked by aromatic guanidine derivatives such as 2-guanidinobenzimidazole (2GBI). Mutation of Hv1 residue F150 to alanine (F150A) was previously found to increase 2GBI apparent binding affinity more than two orders of magnitude. Here, we explore the contribution of aromatic interactions between the inhibitor and the channel in the presence and absence of the F150A mutation, using a combination of electrophysiological recordings, classic mutagenesis, and site-specific incorporation of fluorinated phenylalanines via nonsense suppression methodology. Our data suggest that the increase in apparent binding affinity is due to a rearrangement of the binding site allowed by the smaller residue at position 150. We used this information to design new arginine mimics with improved affinity for the nonrearranged binding site of the wild-type channel. The new compounds, named "Hv1 Inhibitor Flexibles" (HIFs), consist of two "prongs," an aminoimidazole ring, and an aromatic group connected by extended flexible linkers. Some HIF compounds display inhibitory properties that are superior to those of 2GBI, thus providing a promising scaffold for further development of high-affinity Hv1 inhibitors.

Voltage-gated sodium, potassium, and calcium channels consist of four voltage-sensing domains (VSDs) that surround a central pore domain and transition from a down state to an up state in response to membrane depolarization. While many types of drugs bind pore domains, the number of organic molecules known to bind VSDs is limited. The Hv1 voltage-gated proton channel is made of two VSDs and does not contain a pore domain, providing a simplified model for studying how small ligands interact with VSDs. Here, we describe a ligand, named HIF, that interacts with the Hv1 VSD in the up and down states. We find that HIF rapidly inhibits proton conduction in the up state by blocking the open channel, as previously described for 2-guanidinobenzimidazole and its derivatives. HIF, however, interacts with a site slowly accessible in the down state. Functional studies and MD simulations suggest that this interaction traps the compound in a narrow pocket lined with charged residues within the VSD intracellular vestibule, which results in slow recovery from inhibition. Our findings point to a "wrench in gears" mechanism whereby side chains within the binding pocket trap the compound as the teeth of interlocking gears. We propose that the use of screening strategies designed to target binding sites with slow accessibility, similar to the one identified here, could lead to the discovery of new ligands capable of interacting with VSDs of other voltage-gated ion channels in the down state.

The voltage-gated proton channel Hv1 mediates efflux of protons from the cell. Hv1 integrally contributes to various physiological processes including pH homeostasis and the respiratory burst of phagocytes. Inhibition of Hv1 may provide therapeutic avenues for the treatment of inflammatory diseases, breast cancer, and ischemic brain damage. In this work, we investigate two prototypical Hv1 inhibitors, 2-guanidinobenzimidazole (2GBI), and 5-chloro-2-guanidinobenzimidazole (GBIC), from an experimentally screened class of guanidine derivatives. Both compounds block proton conduction by binding the same site located on the intracellular side of the channel. However, when added to the extracellular medium, the compounds strongly differ in their ability to inhibit proton conduction, suggesting substantial differences in membrane permeability. Here, we compute the potential of mean force for each compound to permeate through the membrane using atomistic molecular dynamics simulations with the adaptive biasing force method. Our results rationalize the putative distinction between these two blockers with respect to their abilities to permeate the cellular membrane.

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Hv1 is a voltage-gated proton channel whose main function is to facilitate extrusion of protons from the cell. The development of effective channel blockers for Hv1 can lead to new therapeutics for the treatment of maladies related to Hv1 dysfunction. Although the mechanism of proton permeation in Hv1 remains to be elucidated, a series of small molecules have been discovered to inhibit Hv1. Here, we computed relative binding free energies of a prototypical Hv1 blocker on a model of human Hv1 in an open state. We used alchemical free energy perturbation techniques based on atomistic molecular dynamics simulations. The results support our proposed open state model and shed light on the preferred tautomeric state of the channel blocker. This work lays the groundwork for future studies on adapting the blocker molecule for more effective inhibition of Hv1.

The voltage-gated Hv1 proton channel is a ubiquitous membrane protein that has roles in a variety of cellular processes, including proton extrusion, pH regulation, production of reactive oxygen species, proliferation of cancer cells, and increased brain damage during ischemic stroke. A crystal structure of an Hv1 construct in a putative closed state has been reported, and structural models for the channel open state have been proposed, but a complete characterization of the Hv1 conformational dynamics under an applied membrane potential has been elusive. We report structural models of the Hv1 voltage-sensing domain (VSD), both in a hyperpolarized state and a depolarized state resulting from voltage-dependent conformational changes during a 10-μs-timescale atomistic molecular dynamics simulation in an explicit membrane environment. In response to a depolarizing membrane potential, the S4 helix undergoes an outward displacement, leading to changes in the VSD internal salt-bridge network, resulting in a reshaping of the permeation pathway and a significant increase in hydrogen bond connectivity throughout the channel. The total gating charge displacement associated with this transition is consistent with experimental estimates. Molecular docking calculations confirm the proposed mechanism for the inhibitory action of 2-guanidinobenzimidazole (2GBI) derived from electrophysiological measurements and mutagenesis. The depolarized structural model is also consistent with the formation of a metal bridge between residues located in the core of the VSD. Taken together, our results suggest that these structural models are representative of the closed and open states of the Hv1 channel.

BackgroundThe eye lens crystallins are highly soluble proteins that are required to last the lifespan of an organism due to low protein turnover in the lens. Crystallin aggregation leads to formation of light-scattering aggregates known as cataract. The G18V mutation of human γS-crystallin (γS-G18V), which is associated with childhood-onset cataract, causes structural changes throughout the N-terminal domain and increases aggregation propensity. The holdase chaperone protein αB-crystallin does not interact with wild-type γS-crystallin, but does bind its G18V variant. The specific molecular determinants of αB-crystallin binding to client proteins is incompletely charcterized. Here, a new variant of γS, γS-G18A, was created to test the limits of αB-crystallin selectivity.MethodsMolecular dynamics simulations were used to investigate the structure and dynamics of γS-G18A. The overall fold of γS-G18A was assessed by circular dichroism (CD) spectroscopy and intrinsic tryptophan fluorescence. Its thermal unfolding temperature and aggregation propensity were characterized by CD and DLS, respectively. Solution-state NMR was used to characterize interactions between αB-crystallin and γS-G18A.ResultsγS-G18A exhibits minimal structural changes, but has compromised thermal stability relative to γS-WT. The placement of alanine, rather than valine, at this highly conserved glycine position produces minor changes in hydrophobic surface exposure. However, human αB-crystallin does not bind the G18A variant, in contrast to previous observations for γS-G18V, which aggregates at physiological temperature.ConclusionsαB-crystallin is capable of distinguishing between aggregation-prone and function-preserving variants, and recognizing the transient unfolding or minor conformers that lead to aggregation in the disease-related variant.General significanceHuman αB-crystallin distinguishes between highly similar variants of a structural crystallin, binding the cataract-related γS-G18V variant, but not the function-preserving γS-G18A variant, which is monomeric at physiological temperature.

Available experimental techniques cannot determine high-resolution three-dimensional structures of membrane proteins under a transmembrane voltage. Hence, the mechanism by which voltage-gated cation channels couple conformational changes within the four voltage sensor domains, in response to either depolarizing or polarizing transmembrane voltages, to opening or closing of the pore domain's ion channel remains unresolved. Single-membrane specimens, composed of a phospholipid bilayer containing a vectorially oriented voltage-gated K+ channel protein at high in-plane density tethered to the surface of an inorganic multilayer substrate, were developed to allow the application of transmembrane voltages in an electrochemical cell. Time-resolved neutron reflectivity experiments, enhanced by interferometry enabled by the multilayer substrate, were employed to provide directly the low-resolution profile structures of the membrane containing the vectorially oriented voltage-gated K+ channel for the activated, open and deactivated, closed states of the channel under depolarizing and hyperpolarizing transmembrane voltages applied cyclically. The profile structures of these single membranes were dominated by the voltage-gated K+ channel protein because of the high in-plane density. Importantly, the use of neutrons allowed the determination of the voltage-dependent changes in both the profile structure of the membrane and the distribution of water within the profile structure. These two key experimental results were then compared to those predicted by three computational modeling approaches for the activated, open and deactivated, closed states of three different voltage-gated K+ channels in hydrated phospholipid bilayer membrane environments. Of the three modeling approaches investigated, only one state-of-the-art molecular dynamics simulation that directly predicted the response of a voltage-gated K+ channel within a phospholipid bilayer membrane to applied transmembrane voltages by utilizing very long trajectories was found to be in agreement with the two key experimental results provided by the time-resolved neutron interferometry experiments.

Aquaporin 0 (AQP0) is essential for eye lens homeostasis as is regulation of its water permeability by Ca2+, which occurs through interactions with calmodulin (CaM), but the underlying molecular mechanisms are not well understood. Here, we use molecular dynamics (MD) simulations on the microsecond timescale under an osmotic gradient to explicitly model water permeation through the AQP0 channel. To identify any structural features that are specific to water permeation through AQP0, we also performed simulations of aquaporin 1 (AQP1) and a pure mixed lipid bilayer under the same conditions. The relative single-channel water osmotic permeability coefficients (pf) calculated from all of our simulations are in reasonable agreement with experiment. Our simulations allowed us to characterize the dynamics of the key structural elements that modulate the diffusion of water single-files through the AQP0 and AQP1 pores. We find that CaM binding influences the collective dynamics of the whole AQP0 tetramer, promoting the closing of both the extracellular and intracellular gates by inducing cooperativity between neighboring subunits.

Indoor surfaces are often coated with organic compounds yet a molecular understanding of what drives these interactions is poorly understood. Herein, the adsorption and desorption of limonene, an organic compound found in indoor environments, on hydroxylated silica (SiO2) surfaces, used to mimic indoor glass surfaces, is investigated by combining vibrational spectroscopy, atomistic computer simulations and kinetic modeling. Infrared spectroscopy shows the interaction involves hydrogen-bonding between limonene and surface O-H groups. Atomistic molecular dynamics (MD) simulations confirm the existence of π-hydrogen bonding interactions, with one or two hydrogen bonds between the silica O-H groups and the carbon-carbon double bonds, roughly one third of the time. The concentration and temperature dependent adsorption/desorption kinetics as measured by infrared spectroscopy were reproduced with a kinetic model, yielding the adsorption enthalpy of ∼55 kJ mol-1, which is consistent with the value derived from the MD simulations. Importantly, this integrated experimental, theoretical and kinetic modeling study constitutes a conceptual framework for understanding the interaction of organic compounds with indoor relevant surfaces and thus provides important insights into our understanding of indoor air chemistry and indoor air quality.

In this simulation study, we investigate the influence of biomolecular confinement on dynamical processes in water. We compare water confined in a membrane protein nanopore at room temperature to pure liquid water at low temperatures with respect to structural relaxations, intermolecular vibrations, and the propagation of collective modes. We observe distinct potential energy landscapes experienced by water molecules in the two environments, which nevertheless result in comparable hydrogen bond lifetimes and sound propagation velocities. Hence, we show that a viscoelastic argument that links slow rearrangements of the water-hydrogen bond network to ice-like collective properties applies to both, the pure liquid and biologically confined water, irrespective of differences in the microscopic structure.

Dynamic disorder of the lipid bilayer presents a challenge for establishing structure-function relationships in membranous systems. The resulting structural heterogeneity is especially evident for peripheral and spontaneously inserting membrane proteins, which are not constrained by the well-defined transmembrane topology and exert their action in the context of intimate interaction with lipids. Here, we propose a concerted approach combining depth-dependent fluorescence quenching with Molecular Dynamics simulation to decipher dynamic interactions of membrane proteins with the lipid bilayers. We apply this approach to characterize membrane-mediated action of the diphtheria toxin translocation domain. First, we use a combination of the steady-state and time-resolved fluorescence spectroscopy to characterize bilayer penetration of the NBD probe selectively attached to different sites of the protein into membranes containing lipid-attached nitroxyl quenching groups. The constructed quenching profiles are analyzed with the Distribution Analysis methodology allowing for accurate determination of transverse distribution of the probe. The results obtained for 12 NBD-labeled single-Cys mutants are consistent with the so-called Open-Channel topology model. The experimentally determined quenching profiles for labeling sites corresponding to L350, N373, and P378 were used as initial constraints for positioning TH8-9 hairpin into the lipid bilayer for Molecular Dynamics simulation. Finally, we used alchemical free energy calculations to characterize protonation of E362 in soluble translocation domain and membrane-inserted conformation of its TH8-9 fragment. Our results indicate that membrane partitioning of the neutral E362 is more favorable energetically (by ~ 6 kcal/mol), but causes stronger perturbation of the bilayer, than the charged E362.

It is now well established by numerous experimental and computational studies that the adsorption propensities of inorganic anions conform to the Hofmeister series. The adsorption propensities of inorganic cations, such as the alkali metal cations, have received relatively little attention. Here we use a combination of liquid-jet X-ray photoelectron experiments and molecular dynamics simulations to investigate the behavior of K+ and Li+ ions near the interfaces of their aqueous solutions with halide ions. Both the experiments and the simulations show that Li+ adsorbs to the aqueous solution-vapor interface, while K+ does not. Thus, we provide experimental validation of the "surfactant-like" behavior of Li+ predicted by previous simulation studies. Furthermore, we use our simulations to trace the difference in the adsorption of K+ and Li+ ions to a difference in the resilience of their hydration shells.

Hydrophobic amino acids are abundant in transmembrane (TM) helices of membrane proteins. Charged residues are sparse, apparently due to the unfavorable energetic cost of partitioning charges into nonpolar phases. Nevertheless, conserved arginine residues within TM helices regulate vital functions, such as ion channel voltage gating and integrin receptor inactivation. The energetic cost of arginine in various positions along hydrophobic helices has been controversial. Potential of mean force (PMF) calculations from atomistic molecular dynamics simulations predict very large energetic penalties, while in vitro experiments with Sec61 translocons indicate much smaller penalties, even for arginine in the center of hydrophobic TM helices. Resolution of this conflict has proved difficult, because the in vitro assay utilizes the complex Sec61 translocon, while the PMF calculations rely on the choice of simulation system and reaction coordinate. Here we present the results of computational and experimental studies that permit direct comparison with the Sec61 translocon results. We find that the Sec61 translocon mediates less efficient membrane insertion of Arg-containing TM helices compared with our computational and experimental bilayer-insertion results. In the simulations, a combination of arginine snorkeling, bilayer deformation, and peptide tilting is sufficient to lower the penalty of Arg insertion to an extent such that a hydrophobic TM helix with a central Arg residue readily inserts into a model membrane. Less favorable insertion by the translocon may be due to the decreased fluidity of the endoplasmic reticulum (ER) membrane compared with pure palmitoyloleoyl-phosphocholine (POPC). Nevertheless, our results provide an explanation for the differences between PMF- and experiment-based penalties for Arg burial.

The dimeric 44-residue E5 protein of bovine papillomavirus is the smallest known naturally occurring oncoprotein. This transmembrane protein binds to the transmembrane domain (TMD) of the platelet-derived growth factor β receptor (PDGFβR), causing dimerization and activation of the receptor. Here, we use Rosetta membrane modeling and all-atom molecular dynamics simulations in a membrane environment to develop a chemically detailed model of the E5 protein/PDGFβR complex. In this model, an active dimer of the PDGFβR TMD is sandwiched between two dimers of the E5 protein. Biochemical experiments showed that the major PDGFβR TMD complex in mouse cells contains two E5 dimers and that binding the PDGFβR TMD to the E5 protein is necessary and sufficient to recruit both E5 dimers into the complex. These results demonstrate how E5 binding induces receptor dimerization and define a molecular mechanism of receptor activation based on specific interactions between TMDs.

Voltage-sensing domains (VSDs) sense changes in the membrane electrostatic potential and, through conformational changes, regulate a specific function. The VSDs of wild-type voltage-dependent K+, Na+, and Ca2+ channels do not conduct ions, but they can become ion-permeable through pathological mutations in the VSD. Relatively little is known about the underlying mechanisms of conduction through VSDs. The most detailed studies have been performed on Shaker K+ channel variants in which ion conduction through the VSD is manifested in electrophysiology experiments as a voltage-dependent inward current, the so-called omega current, which appears when the VSDs are in their resting state conformation. Only monovalent cations appear to permeate the Shaker VSD via a pathway that is believed to be, at least in part, the same as that followed by the S4 basic side chains during voltage-dependent activation. We performed μs-time scale atomistic molecular dynamics simulations of a cation-conducting variant of the Shaker VSD under applied electric fields in an experimentally validated resting-state conformation, embedded in a lipid bilayer surrounded by solutions containing guanidinium chloride or potassium chloride. Our simulations provide insights into the Shaker VSD permeation pathway, the protein-ion interactions that control permeation kinetics, and the mechanism of voltage-dependent activation of voltage-gated ion channels.

Human voltage-gated proton (hHv1) channels are crucial for regulating essential biological processes such as immune cell respiratory burst, sperm capacitation, and cancer cell migration. Despite the significant concentration difference between protons and other ions in physiological conditions, hHv1 demonstrates remarkable proton selectivity. Our calculations of single-proton, cation, and anion permeation free energy profiles quantitatively demonstrate that the proton selectivity of the wild-type channel originates from its strong proton affinity via the titration of the key residues D112 and D174, although the channel imposes similar kinetic blocking effects for protons compared to other ions. A two-proton knock-on model is proposed to mathematically explain the electrophysiological measurements of the pH-dependent proton conductance in the conductive state. Moreover, it is shown that the anion selectivity of the D112N mutant channel is tied to impaired proton transport and substantial anion leakage. [ABSTRACT FROM AUTHOR]

Aquaporin 0 (AQP0), the major intrinsic protein of the eye lens, plays a vital role in maintaining lens clarity by facilitating the transport of water across lens fiber cell membranes. AQP0 reduces its osmotic water permeability constant (Pf) in response to increases in the external calcium concentration, an effect that is mediated by an interaction with the calcium-binding messenger protein, calmodulin (CaM), and phosphorylation of the CaM-binding site abolishes calcium sensitivity. Despite recent structural characterization of the AQP0-CaM complex, the mechanism by which CaM modulates AQP0 remains poorly understood. By combining atomistic molecular dynamics simulations and oocyte permeability assays, we conclude that serine phosphorylation of AQP0 does not inhibit CaM binding to the whole AQP0 protein. Instead, AQP0 phosphorylation alters calcium sensitivity by modifying the AQP0-CaM interaction interface, particularly at an arginine-rich loop that connects the fourth and fifth transmembrane helices. This previously unexplored loop, which sits outside of the canonical CaM-binding site on the AQP0 cytosolic face, mechanically couples CaM to the pore-gating residues of the second constriction site. We show that this allosteric loop is vital for CaM regulation of the channels, facilitating cooperativity between adjacent subunits and regulating factors such as serine phosphorylation. Similar allosteric interactions may also mediate CaM modulation of the properties of other CaM-regulated proteins.

Organized as bilayers, phospholipids are the fundamental building blocks of cellular membranes and determine many of their biological functions. Interactions between the two leaflets of the bilayer (interleaflet coupling) have been implicated in the passage of information through membranes. However, physically, the meaning of interleaflet coupling is ill defined and lacks a structural basis. Using all-atom molecular dynamics simulations of fluid phospholipid bilayers of five different lipids with differing degrees of acyl-chain asymmetry, we have examined interleaflet mixing to gain insights into coupling. Reasoning that the transbilayer distribution of terminal methyl groups is an appropriate measure of interleaflet mixing, we calculated the transbilayer distributions of the acyl chain terminal methyl groups for five lipids: dioleoylphosphatidylcholine (DOPC), palmitoyloleoylphosphatidylcholine (POPC), stearoyloleoylphosphatidylcholine (SOPC), oleoylmyristoylphosphatidylcholine (OMPC), and dimyristoylphosphatidylcholine (DMPC). We observed in all cases very strong mixing across the bilayer midplane that diminished somewhat with increasing acyl-chain ordering defined by methylene order parameters. A hallmark of the interleaflet coupling idea is complementarity, which postulates that lipids with short alkyl chains in one leaflet will preferentially associate with lipids with long alkyl chains in the other leaflet. Our results suggest a much more complicated picture for thermally disordered bilayers that we call distributed complementarity, as measured by the difference in the peak positions of the sn-1 and sn-2 methyl distributions in the same leaflet.

Leishmania major, the parasitic causative agent of leishmaniasis, produces a heme peroxidase (LmP), which catalyzes the peroxidation of mitochondrial cytochrome c (LmCytc) for protection from reactive oxygen species produced by the host. The association of LmP and LmCytc, which is known from kinetics measurements to be very fast (∼10(8) M(-1) s(-1)), does not involve major conformational changes and has been suggested to be dominated by electrostatic interactions. We used Brownian dynamics simulations to investigate the mechanism of formation of the LmP-LmCytc complex. Our simulations confirm the importance of electrostatic interactions involving the negatively charged D211 residue at the LmP active site, and reveal a previously unrecognized role in complex formation for negatively charged residues in helix A of LmP. The crystal structure of the D211N mutant of LmP reported herein is essentially identical to that of wild-type LmP, reinforcing the notion that it is the loss of charge at the active site, and not a change in structure, that reduces the association rate of the D211N variant of LmP. The Brownian dynamics simulations further show that complex formation occurs via a "bind and crawl" mechanism, in which LmCytc first docks to a location on helix A that is far from the active site, forming an initial encounter complex, and then moves along helix A to the active site. An atomistic molecular dynamics simulation confirms the helix A binding site, and steady state activity assays and stopped-flow kinetics measurements confirm the role of helix A charges in the association mechanism.

The heterotrimeric SecY translocon complex is required for the cotranslational assembly of membrane proteins in bacteria and archaea. The insertion of transmembrane (TM) segments during nascent-chain passage through the translocon is generally viewed as a simple partitioning process between the water-filled translocon and membrane lipid bilayer, suggesting that partitioning is driven by the hydrophobic effect. Indeed, the apparent free energy of partitioning of unnatural aliphatic amino acids on TM segments is proportional to accessible surface area, which is a hallmark of the hydrophobic effect [Öjemalm K, et al. (2011) Proc Natl Acad Sci USA 108(31):E359-E364]. However, the apparent partitioning solvation parameter is less than one-half the value expected for simple bulk partitioning, suggesting that the water in the translocon departs from bulk behavior. To examine the state of water in a SecY translocon complex embedded in a lipid bilayer, we carried out all-atom molecular-dynamics simulations of the Pyrococcus furiosus SecYE, which was determined to be in a "primed" open state [Egea PF, Stroud RM (2010) Proc Natl Acad Sci USA 107(40):17182-17187]. Remarkably, SecYE remained in this state throughout our 450-ns simulation. Water molecules within SecY exhibited anomalous diffusion, had highly retarded rotational dynamics, and aligned their dipoles along the SecY transmembrane axis. The translocon is therefore not a simple water-filled pore, which raises the question of how anomalous water behavior affects the mechanism of translocon function and, more generally, the partitioning of hydrophobic molecules. Because large water-filled cavities are found in many membrane proteins, our findings may have broader implications.

The paired helical filaments (PHF) formed by the intrinsically disordered human protein tau are one of the pathological hallmarks of Alzheimer disease. PHF are fibers of amyloid nature that are composed of a rigid core and an unstructured fuzzy coat. The mechanisms of fiber formation, in particular the role that hydration water might play, remain poorly understood. We combined protein deuteration, neutron scattering, and all-atom molecular dynamics simulations to study the dynamics of hydration water at the surface of fibers formed by the full-length human protein htau40. In comparison with monomeric tau, hydration water on the surface of tau fibers is more mobile, as evidenced by an increased fraction of translationally diffusing water molecules, a higher diffusion coefficient, and increased mean-squared displacements in neutron scattering experiments. Fibers formed by the hexapeptide (306)VQIVYK(311) were taken as a model for the tau fiber core and studied by molecular dynamics simulations, revealing that hydration water dynamics around the core domain is significantly reduced after fiber formation. Thus, an increase in water dynamics around the fuzzy coat is proposed to be at the origin of the experimentally observed increase in hydration water dynamics around the entire tau fiber. The observed increase in hydration water dynamics is suggested to promote fiber formation through entropic effects. Detection of the enhanced hydration water mobility around tau fibers is conjectured to potentially contribute to the early diagnosis of Alzheimer patients by diffusion MRI.

Upon endoplasmic reticulum Ca(2+) store depletion, Orai channels in the plasma membrane are activated directly by endoplasmic reticulum-resident STIM proteins to generate the Ca(2+)-selective, Ca(2+) release-activated Ca(2+) (CRAC) current. After the molecular identification of Orai, a plethora of functional and biochemical studies sought to compare Orai homologs, determine their stoichiometry, identify structural domains responsible for the biophysical fingerprint of the CRAC current, identify the physiological functions, and investigate Orai homologs as potential therapeutic targets. Subsequently, the solved crystal structure of Drosophila Orai (dOrai) substantiated many findings from structure-function studies, but also revealed an unexpected hexameric structure. In this review, we explore Orai channels as elucidated by functional and biochemical studies, analyze the dOrai crystal structure and its implications for Orai channel function, and present newly available information from molecular dynamics simulations that shed light on Orai channel gating and permeation.

The voltage sensor domain (VSD) of voltage-gated cation (e.g., Na(+), K(+)) channels central to neurological signal transmission can function as a distinct module. When linked to an otherwise voltage-insensitive, ion-selective membrane pore, the VSD imparts voltage sensitivity to the channel. Proteins homologous with the VSD have recently been found to function themselves as voltage-gated proton channels or to impart voltage sensitivity to enzymes. Determining the conformational changes associated with voltage gating in the VSD itself in the absence of a pore domain thereby gains importance. We report the direct measurement of changes in the scattering-length density (SLD) profile of the VSD protein, vectorially oriented within a reconstituted phospholipid bilayer membrane, as a function of the transmembrane electric potential by time-resolved X-ray and neutron interferometry. The changes in the experimental SLD profiles for both polarizing and depolarizing potentials with respect to zero potential were found to extend over the entire length of the isolated VSD's profile structure. The characteristics of the changes observed were in qualitative agreement with molecular dynamics simulations of a related membrane system, suggesting an initial interpretation of these changes in terms of the VSD's atomic-level 3-D structure.

The crystal structure of the flavin mononucleotide (FMN)-containing redox partner to P450cin, cindoxin (Cdx), has been determined to 1.3 Å resolution. The overall structure is similar to that of the FMN domain of human cytochrome P450 reductase. A Brownian dynamics-molecular dynamics docking method was used to produce a model of Cdx with its redox partner, P450cin. This Cdx-P450cin model highlights the potential importance of Cdx Tyr96 in bridging the FMN and heme cofactors as well P450cin Arg102 and Arg346. Each of the single-site Ala mutants exhibits ~10% of the wild-type activity, thus demonstrating the importance of these residues for binding and/or electron transfer. In the well-studied P450cam system, redox partner binding stabilizes the open low-spin conformation of P450cam and greatly decreases the stability of the oxy complex. In sharp contrast, Cdx does not shift P450cin to a low-spin state, although the stability of oxy-P450cin is decreased 10-fold in the presence of Cdx. This indicates that Cdx may have a modest effect on the open-closed equilibrium in P450cin compared to that in P450cam. It has been postulated that part of the effector role of Pdx on P450cam is to promote a significant structural change that makes available a proton relay network involving Asp251 required for O2 activation. The structure around the corresponding Asp in P450cin, Asp241, provides a possible structural reason for why P450cin is less dependent on its redox partner for functionally important structural changes.

The major removal processes for gaseous nitric acid (HNO3) in the atmosphere are dry and wet deposition onto various surfaces. The surface in the boundary layer is often covered with organic films, but the interaction of gaseous HNO3 with them is not well understood. To better understand the factors controlling the uptake of gaseous nitric acid and its dissociation in organic films, studies were carried out using single component and mixtures of C8 and C18 alkyl self-assembled monolayers (SAMs) attached to a germanium (Ge) attenuated total reflectance (ATR) crystal upon which a thin layer of SiOx had been deposited. For comparison, diffuse reflectance infrared Fourier transform spectrometry (DRIFTS) studies were also carried out using a C18 SAM attached to the native oxide layer on the surface of silicon powder. These studies show that the alkyl chain length and order/disorder of the SAMs does not significantly affect the uptake or dissociation/recombination of molecular HNO3. Thus, independent of the nature of the SAM, molecular HNO3 is observed up to 70-90% relative humidity. After dissociation, molecular HNO3 is regenerated on all SAM surfaces when water is removed. Results of molecular dynamics simulations are consistent with experiments and show that defects and pores on the surfaces control the uptake, dissociation and recombination of molecular HNO3. Organic films on surfaces in the boundary layer will certainly be more irregular and less ordered than SAMs studied here, therefore undissociated HNO3 may be present on surfaces in the boundary layer to a greater extent than previously thought. The combination of this observation with the results of recent studies showing enhanced photolysis of nitric acid on surfaces suggests that renoxification of deposited nitric acid may need to be taken into account in atmospheric models.

Calmodulin (CaM) is a universal regulatory protein that communicates the presence of calcium to its molecular targets and correspondingly modulates their function. This key signaling protein is important for controlling the activity of hundreds of membrane channels and transporters. However, understanding of the structural mechanisms driving CaM regulation of full-length membrane proteins has remained elusive. In this study, we determined the pseudoatomic structure of full-length mammalian aquaporin-0 (AQP0, Bos taurus) in complex with CaM, using EM to elucidate how this signaling protein modulates water-channel function. Molecular dynamics and functional mutation studies reveal how CaM binding inhibits AQP0 water permeability by allosterically closing the cytoplasmic gate of AQP0. Our mechanistic model provides new insight, only possible in the context of the fully assembled channel, into how CaM regulates multimeric channels by facilitating cooperativity between adjacent subunits.

The Symposium ‘Frontiers in membrane and membrane protein biophysics: experiments and theory’, held this year at the University of California, Irvine (August 19–20), celebrated the 70th Birthday of Stephen H. White by bringing together distinguished experimentalists and theoreticians to discuss the state of the art and future challenges in the field of membrane and membrane protein biophysics. The meeting and this special issue highlight the highly interdisciplinary nature of membrane and membrane protein biophysics, and the tremendous contributions that S. H. White and his lab have brought to the field.

Several laboratories have carried out molecular dynamics (MD) simulations of arginine interactions with lipid bilayers and found that the energetic cost of placing arginine in lipid bilayers is an order of magnitude greater than observed in molecular biology experiments in which Arg-containing transmembrane helices are inserted across the endoplasmic reticulum membrane by the Sec61 translocon. We attempt here to reconcile the results of the two approaches. We first present MD simulations of guanidinium groups alone in lipid bilayers, and then, to mimic the molecular biology experiments, we present simulations of hydrophobic helices containing single Arg residues at different positions along the helix. We discuss the simulation results in the context of molecular biology results and show that the energetic discrepancy is reduced, but not eliminated, by considering free energy differences between Arg at the interface and at the center of the model helices. The reduction occurs because Arg snorkeling to the interface prevents Arg from residing in the bilayer center where the energetic cost of desolvation is highest. We then show that the problem with MD simulations is that they measure water-to-bilayer free energies, whereas the molecular biology experiments measure the energetics of partitioning from translocon to bilayer, which raises the fundamental question of the relationship between water-to-bilayer and water-to-translocon partitioning. We present two thermodynamic scenarios as a foundation for reconciliation of the simulation and molecular biology results. The simplest scenario is that translocon-to-bilayer partitioning is independent of water-to-bilayer partitioning; there is no thermodynamic cycle connecting the two paths.

Despite the growing number of atomic-resolution membrane protein structures, direct structural information about proteins in their native membrane environment is scarce. This problem is particularly relevant in the case of the highly charged S1-S4 voltage-sensing domains responsible for nerve impulses, where interactions with the lipid bilayer are critical for the function of voltage-activated ion channels. Here we use neutron diffraction, solid-state nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics simulations to investigate the structure and hydration of bilayer membranes containing S1-S4 voltage-sensing domains. Our results show that voltage sensors adopt transmembrane orientations and cause a modest reshaping of the surrounding lipid bilayer, and that water molecules intimately interact with the protein within the membrane. These structural findings indicate that voltage sensors have evolved to interact with the lipid membrane while keeping energetic and structural perturbations to a minimum, and that water penetrates the membrane, to hydrate charged residues and shape the transmembrane electric field.

We have investigated the sensitivity of two-dimensional infrared (2D IR) spectroscopy to peptide helicity with an experimental and theoretical study of Z-[l-(alphaMe)Val](8)-OtBu in CDCl(3). 2D IR experiments were carried out in the amide-I region under the parallel and the double-crossed polarization configurations. In the latter polarization configuration, the 2D spectra taken with the rephasing and nonrephasing pulse sequences exhibit a doublet feature and a single peak, respectively. These cross-peak patterns are highly sensitive to the underlying peptide structure. Spectral calculations were performed on the basis of a vibrational exciton model, with the local mode frequencies and couplings calculated from snapshots of molecular dynamics (MD) simulation trajectories using six different models for the Hamiltonian. Conformationally variant segments of the MD trajectory, while reproducing the main features of the experimental spectra, are characterized by extraneous features, suggesting that the structural ensembles sampled by the simulation are too broad. By imposing periodic restraints on the peptide dihedral angles with the crystal structure as a reference, much better agreement between the measured and the calculated spectra was achieved. The result indicates that the structure of Z-[l-(alphaMe)Val](8)-OtBu in CDCl(3) is a fully developed 3(10)-helix with only a small fraction of alpha-helical or nonhelical conformations in the middle of the peptide. Of the four different combinations of pulse sequences and polarization configurations, the nonrephasing double-crossed polarization 2D IR spectrum exhibits the highest sensitivity in detecting conformational variation. Of the six local mode frequency models tested, the electrostatic maps of Mukamel and Cho perform the best. Our results show that the high sensitivity of 2D IR spectroscopy can provide a useful basis for developing methods to improve the sampling accuracy of force fields and for characterizing the relative merits of the different protocols for the Hamiltonian calculation.

Heterogeneous reactions of sea salt aerosol with various oxides of nitrogen lead to replacement of chloride ion by nitrate ion. Studies of the photochemistry of a model system were carried out using deliquesced mixtures of NaCl and NaNO3 on a Teflon substrate. Varying molar ratios of NaCl to NaNO3 (1 : 9 Cl- : NO3-, 1 : 1 Cl- : NO3-, 3 : 1 Cl- : NO3-, 9 : 1 Cl- : NO3-) and NaNO3 at the same total concentration were irradiated in air at 299 +/- 3 K and at a relative humidity of 75 +/- 8% using broadband UVB light (270-380 nm). Gaseous NO2 production was measured as a function of time using a chemiluminescence NO(y) detector. Surprisingly, an enhanced yield of NO2 was observed as the chloride to nitrate ratio increased. Molecular dynamics (MD) simulations show that as the Cl- : NO3- ratio increases, the nitrate ions are drawn closer to the interface due to the existence of a double layer of interfacial Cl- and subsurface Na+. This leads to a decreased solvent cage effect when the nitrate ion photodissociates to NO2+O*-, increasing the effective quantum yield and hence the production of gaseous NO2. The implications of enhanced NO2 and likely OH production as sea salt aerosols become processed in the atmosphere are discussed.

The photooxidation of 0.6-0.9 ppm alpha-pinene in the presence of a deliquesced thin film of NaNO(3), and for comparison increasing concentrations of NO(2), was studied in a 100 L Teflon(R) chamber at relative humidities from 72-88% and temperatures from 296-304 K. The loss of alpha-pinene and the formation of gaseous products were followed with time using proton transfer mass spectrometry. The yields of gas phase products were smaller in the NaNO(3) experiments than in NO(2) experiments. In addition, pinonic acid, pinic acid, trans-sobrerol and other unidentified products were detected in the extracts of the wall washings only for the NaNO(3) photolysis. These data indicate enhanced loss of alpha-pinene at the NaNO(3) thin film during photolysis. Supporting the experimental results are molecular dynamics simulations which predict that alpha-pinene has an affinity for the surface of the deliquesced nitrate thin film, enhancing the opportunity for oxidation of the impinging organic gas during the nitrate photolysis. This new mechanism of oxidation of organics may be partially responsible for the correlation between nitrate and the organic component of particles observed in many field studies, and may also contribute to the missing source of SOA needed to reconcile model predictions and field measurements. In addition, photolysis of nitrate on surfaces in the boundary layer may lead to oxidation of co-adsorbed organics.