Alpha-helical hydrophobic polypeptides form proton-selective channels in lipid bilayers

Prebiotic Peptides: Molecular Hubs in the Origin of Life

The fundamental roles that peptides and proteins play in today's biology makes it almost indisputable that peptides were key players in the origin of life. Insofar as it is appropriate to extrapolate back from extant biology to the prebiotic world, one must acknowledge the critical importance that interconnected molecular networks, likely with peptides as key components, would have played in life's origin. In this review, we summarize chemical processes involving peptides that could have contributed to early chemical evolution, with an emphasis on molecular interactions between peptides and other classes of organic molecules. We first summarize mechanisms by which amino acids and similar building blocks could have been produced and elaborated into proto-peptides. Next, non-covalent interactions of peptides with other peptides as well as with nucleic acids, lipids, carbohydrates, metal ions, and aromatic molecules are discussed in relation to the possible roles of such interactions in chemical evolution of structure and function. Finally, we describe research involving structural alternatives to peptides and covalent adducts between amino acids/peptides and other classes of molecules. We propose that ample future breakthroughs in origin-of-life chemistry will stem from investigations of interconnected chemical systems in which synergistic interactions between different classes of molecules emerge.

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na⁺ and K⁺ ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the K_v1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine’s NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of −70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

Reconstituted human TPC1 is a proton-permeable ion channel and is activated by NAADP or Ca2+

NAADP potently triggers Ca2+ release from acidic lysosomal and endolysosomal Ca2+ stores. Human two-pore channels (TPC1 and TPC2), which are located on these stores, are involved in this process, but there is controversy over whether TPC1 and TPC2 constitute the Ca2+ release channels. We therefore examined the single-channel properties of human TPC1 after reconstitution into bilayers of controlled composition. We found that TPC1 was permeable not only to Ca2+ but also to monovalent cations and that permeability to protons was the highest (relative permeability sequence: H+ >> K+ > Na(+) ≥ Ca2+). NAADP or Ca2+ activated TPC1, and the presence of one of these ligands was required for channel activation. The endolysosome-located lipid phosphatidylinositol 3,5-bisphosphate [PI(3,5)P2] had no effect on TPC1 open probability but significantly increased the relative permeability of Na+ to Ca2+ and of H+ to Ca2+. Furthermore, our data showed that, although both TPC1 and TPC2 are stimulated by NAADP, these channels differ in ion selectivity and modulation by Ca2+ and pH. We propose that NAADP triggers H+ release from lysosomes and endolysomes through activation of TPC1, but that the Ca2+ -releasing ability of TPC1 will depend on the ionic composition of the acidic stores and may be influenced by other regulators that affect TPC1 ion permeation.

Toward understanding protocell mechanosensation

Mechanosensitive (MS) channels can prevent bacterial bursting during hypo-osmotic shocks by responding to increases in lateral tension at the membrane level through an integrated and coordinated opening mechanism. Mechanical regulation in protocells could have been one of the first mechanisms to evolve in order to preserve their integrity against changing environmental conditions. How has the rich functional diversity found in present cells been created throughout evolution, and what did the primordial MS channels look like? This review has been written with the aim of identifying which factors may have been important for the appearance of the first osmotic valve in a prebiotic context, and what this valve may have been like. It highlights the mechanical properties of lipid bilayers, the association of peptides as aggregates in membranes, and the conservation of sequence motifs as central aspects to understand the evolution of proteins that gate below the tension required for spontaneous pore formation and membrane rupture. The arguments developed here apply to both MscL and MscS homologs, but could be valid to mechano-susceptible proteins in general.

Membrane transport in primitive cells

Although model protocellular membranes consisting of monoacyl lipids are similar to membranes composed of contemporary diacyl lipids, they differ in at least one important aspect. Model protocellular membranes allow for the passage of polar solutes and thus can potentially support cell-to functions without the aid of transport machinery. The ability to transport polar molecules likely stems from increased lipid dynamics. Selectively permeable vesicle membranes composed of monoacyl lipids allow for many lifelike processes to emerge from a remarkably small set of molecules.

On the way towards ‘basic autonomous agents’: Stochastic simulations of minimal lipid–peptide cells

In this paper, we apply a recently developed stochastic simulation platform to investigate the dynamic behaviour of minimal 'self-(re-)producing' cellular systems. In particular, we study a set of preliminary conditions for appearance of the simplest forms of autonomy in the context of lipid vesicles (more specifically, lipid-peptide vesicles) that enclose an autocatalytic/proto-metabolic reaction network. The problem is approached from a 'bottom-up' perspective, in the sense that we try to show how relatively simple cell components/processes could engage in a far-from-equilibrium dynamics, staying in those conditions thanks to a rudimentary but effective control of the matter-energy flow through it. In this general scenario, basic autonomy and, together with it, minimal agent systems would appear when (hypothetically pre-biological) cellular systems establish molecular trans-membrane mechanisms that allow them to couple internal chemical reactions with transport processes, in a way that they channel/transform external material-energetic resources into their own means and actively regulate boundary conditions (e.g., osmotic gradients, inflow/outflow of different compounds, ...) that are critical for their constitution and persistence as proto-metabolic cells. The results of our simulations indicate that, before that stage is reached, there are a number of relevant issues that have to be carefully analysed and clarified: especially the immediate effects that the insertion of peptide chains (channel precursors) in the lipid bilayer may have in the structural properties of the membrane (elasticity, permeability, ...) and in the overall dynamic behaviour of the cell.

Question 7: modelling minimal 'lipid-peptide' cells

This contribution is aimed to give support to 'bottom-up' approaches to the minimal or early cell research project. Even if, from this perspective, the most simple living cell still seems very far away, the analysis of less complex, infrabiological cellular systems (some of which could be relatively soon implemented in the lab) probably holds the key, or one of the fundamental keys, to the problem of origins of life. On these lines, we propose a simulation model to study the transition from proto-metabolic 'lipid' cells to 'lipid-peptide' cells, as a critical step in which self-reproducing vesicles could develop into more functionalized supramolecular systems.

Self-assembling of peptide/membrane complexes by atomistic molecular dynamics simulations

Model biological membranes consisting of peptide/lipid-bilayer complexes can nowadays be studied by classical molecular dynamics (MD) simulations at atomic detail. In most cases, the simulation starts with an assumed state of a peptide in a preformed bilayer, from which equilibrium configurations are difficult to obtain due to a relatively slow molecular diffusion. As an alternative, we propose an extension of reported work on the self-organization of unordered lipids into bilayers, consisting of including a peptide molecule in the initial random configuration to obtain a membrane-bound peptide simultaneous to the formation of the lipid bilayer. This strategy takes advantage of the fast reorganization of lipids, among themselves and around the peptide, in an aqueous environment. Model peptides of different hydrophobicity, CH3-CO-W2L18W2-NH2 (WL22) and CH3-CO-W2A18W2-NH2 (WA22), in dipalmitoyl-phosphatidylcholine (DPPC) are used as test cases. In the equilibrium states of the peptide/membrane complexes, achieved in time ranges of 50-100 ns, the two peptides behave as expected from experimental and theoretical studies. The strongly hydrophobic WL22 is inserted in a transmembrane configuration and the marginally apolar, alanine-based WA22 is found in two alternative states: transmembrane inserted or parallel to the membrane plane, embedded close to the bilayer interface, with similar stability. This shows that the spontaneous assembly of peptides and lipids is an unbiased and reliable strategy to produce and study models of equilibrated peptide/lipid complexes of unknown membrane-binding mode and topology.

Ion transport across model lipid membranes containing light-harvesting complex II: an effect of light

The effect of light on proton transport across lipid membranes of small unilamellar liposomes containing incorporated major light-harvesting pigment-protein complex of Photosystem II (LHCII) has been studied with the application of pH-sensitive dyes entrapped inside vesicles. Proton permeability coefficient for LHCII-modified membranes was found to be about twice as high as in the case of the control pure lipid vesicles. Illumination of the samples with light absorbed by the LHCII-bound photosynthetic pigments considerably affects the kinetics of proton transport: it increases the rate and decreases the steady-state level of proton gradient across the membranes. The effect was interpreted in terms of heat-induced conformational changes of LHCII molecular structures that affect proton buffering capacity of this protein. Both the control and the LHCII-modified lipid membranes have been found to be practically impermeable to Ca(++) ions, as demonstrated by fluorescence of liposome-entrapped calcium-sensitive probe calcium crimson. The slight differences in the proton transport across the LHCII-containing membranes under the presence of Ca(++) suggest calcium binding to this antenna protein.

Basic autonomy as a fundamental step in the synthesis of life

In the search for the primary roots of autonomy (a pivotal concept in Varela's comprehensive understanding of living beings), the theory of autopoiesis provided an explicit criterion to define minimal life in universal terms, and was taken as a guideline in the research program for the artificial synthesis of biological systems. Acknowledging the invaluable contribution of the autopoietic school to present biological thinking, we offer an alternative way of conceiving the most basic forms of autonomy. We give a bottom-up account of the origins of "self-production" (or self-construction, as we propose to call it), pointing out which are the minimal material and energetic requirements for the constitution of basic autonomous systems. This account is, indeed, committed to the project of developing a general theory of biology, but well grounded in the universal laws of physics and chemistry. We consider that the autopoietic theory was formulated in highly abstract terms and, in order to advance in the implementation of minimal autonomous systems (and, at the same time, make major progress in exploring the origins of life), a more specific characterization of minimal autonomous systems is required. Such a characterization will not be drawn from a review of the autopoietic criteria and terminology (à la Fleischaker) but demands a whole reformulation of the question: a proper naturalization of the concept of autonomy. Finally, we also discuss why basic autonomy, according to our account, is necessary but not sufficient for life, in contrast with Varela's idea that autopoiesis was a necessary and sufficient condition for it.

A computer simulation study of the hydrated proton in a synthetic proton channel

Classical molecular dynamics simulations using the multistate empirical valence bond model for aqueous proton transport were performed to characterize the hydration structure of an excess proton inside a leucine-serine synthetic ion channel, LS2. For such a nonuniform pore size ion channel, it is found that the Zundel ion (H(5)O(2)(+)) solvation structure is generally more stable in narrow channel regions than in wider channel regions, which is in agreement with a recent study on idealized hydrophobic proton channels. However, considerable diversity in the relative stability of the Zundel to Eigen cation (H(9)O(4)(+)) was observed. Three of the five wide channel regions, one located at the channel's center and the other two located near the channel mouths, are found to show extraordinary preference for the Eigen solvation structure. This implies that proton hopping is inhibited in these regions and therefore suggests that these regions may behave as barriers in the proton conducting pathway inside the channel. The proton solvation is also greatly influenced by the local molecular environment of the protein. In particular, the polar side chains of the Ser residues, which are intimately involved in the solvation structure, can greatly influence proton solvation. However, no preference of the influence by the various Ser side chains was found; they can either promote or prevent the formation of certain solvation structures.

Voltage-gated proton channels and other proton transfer pathways

Proton channels exist in a wide variety of membrane proteins where they transport protons rapidly and efficiently. Usually the proton pathway is formed mainly by water molecules present in the protein, but its function is regulated by titratable groups on critical amino acid residues in the pathway. All proton channels conduct protons by a hydrogen-bonded chain mechanism in which the proton hops from one water or titratable group to the next. Voltage-gated proton channels represent a specific subset of proton channels that have voltage- and time-dependent gating like other ion channels. However, they differ from most ion channels in their extraordinarily high selectivity, tiny conductance, strong temperature and deuterium isotope effects on conductance and gating kinetics, and insensitivity to block by steric occlusion. Gating of H(+) channels is regulated tightly by pH and voltage, ensuring that they open only when the electrochemical gradient is outward. Thus they function to extrude acid from cells. H(+) channels are expressed in many cells. During the respiratory burst in phagocytes, H(+) current compensates for electron extrusion by NADPH oxidase. Most evidence indicates that the H(+) channel is not part of the NADPH oxidase complex, but rather is a distinct and as yet unidentified molecule.

How did cells get their size?

Cells exercise size homeostasis, and the origins of their ability to do so is the topic of this essay. Before there were cells, there were protocells. The most basic questions about protocells as objects are: What were they made of, and how big were they? Asking how big they were implies that the answer to the first part includes a boundary. The best candidate for that boundary is a self-assembling lipid bilayer. Therefore, protocells are defined here as Darwinian liposomes-bilayer vesicles with mutable on-board replicases linked to phenotypes. Because liposomes undergo spontaneous fission and fusion, and are subject to osmotic forces, size regulation in the earliest protocells would essentially have been liposome physics. For successful protocells, averting osmotic lysis would have been the first order of business. However, from the outset size mattered too, because of sex and reproduction (i.e., genome mixing and genome copying in entities with phenotypes). Protocell fission and fusion would have blended seamlessly into protocell sex and reproduction, making any gene product that furnished control over protocell size changes doubly adaptive. A recurrent theme is the feedback role of bilayer tension in protocell size control. Ways in which primitive peptides and their aggregates (e.g., channels) might have allowed liposomes to gain improved volume and surface area homeostasis are suggested. Domain-swapped proteins that polymerize as filaments are discussed as the origin of cytoskeleton structures that diversify and stabilize liposome shapes and sizes. Throughout, attention is paid to the question of set points for cell size.

Membrane self-assembly processes: steps toward the first cellular life

This review addresses the question of the origin of life, with emphasis on plausible boundary structures that may have initially provided cellular compartmentation. Some form of compartmentation is a necessary prerequisite for maintaining the integrity of interdependent molecular systems that are associated with metabolism, and for permitting variations required for speciation. The fact that lipid-bilayer membranes define boundaries of all contemporary living cells suggests that protocellular compartments were likely to have required similar, self-assembled boundaries. Amphiphiles such as short-chain fatty acids, which were presumably available on the early Earth, can self-assemble into stable vesicles that encapsulate hydrophilic solutes with catalytic activity. Their suspensions in aqueous media have therefore been used to investigate nutrient uptake across simple membranes and encapsulated catalyzed reactions, both of which would be essential processes in protocellular life forms.

Membrane insertion and orientation of polyalanine peptides: a (15)N solid-state NMR spectroscopy investigation

Polyalanine-based peptides were prepared by solid-phase peptide synthesis, labeled with (15)N at selected sites, reconstituted into oriented phosphatidylcholine bilayers, and investigated by proton-decoupled (15)N solid-state NMR spectroscopy. The anisotropic (15)N chemical shift is a direct indicator of helix alignment with respect to the membrane normal. The in-plane to transmembrane equilibrium is the focus of this study. Time- and solvent-dependent transmembrane alignments of K(3)A(18)K(3) have been obtained, and these are stabilized when a few alanine residues are replaced with leucine. The results are discussed in the context of a model where polyalanines adopt a variety of configurations, which are interconnected by multiple equilibria. The data indicate hydrophobicity values of alanine close to zero when studied in the context of helical polypeptides (> or =24 residues) and phospholipid bilayers.

Mitochondrial proton leak and the uncoupling protein 1 homologues

Mitochondrial proton leak is the largest single contributor to the standard metabolic rate (SMR) of a rat, accounting for about 20% of SMR. Yet the mechanisms by which proton leak occurs are incompletely understood. The available evidence suggests that both phospholipids and proteins in the mitochondrial inner membrane are important determinants of proton conductance. The uncoupling protein 1 homologues (e.g. UCP2, UCP3) may play a role in mediating proton leak, but it is unlikely they account for all of the observed proton conductance. Experimental data regarding the functions of these proteins include important ambiguities and contradictions which must be addressed before their function can be confirmed. The physiological role of the proton leak, and of the uncoupling protein 1 homologues, remains similarly unclear.

Transmembrane peptide NB of influenza B: a simulation, structure, and conductance study

The putative transmembrane segment of the ion channel forming peptide NB from influenza B was synthesized by standard solid-phase peptide synthesis. Insertion into the planar lipid bilayer revealed ion channel activity with conductance levels of 20, 61, 107, and 142 pS in a 0.5 M KCl buffer solution. In addition, levels at -100 mV show conductances of 251 and 413 pS. A linear current-voltage relation reveals a voltage-independent channel formation. In methanol and in vesicles the peptide appears to adopt an alpha-helical-like structure. Computational models of alpha-helix bundles using N = 4, 5, and 6 NB peptides per bundle revealed water-filled pores after 1 ns of MD simulation in a solvated lipid bilayer. Calculated conductance values [using HOLE (Smart et al. (1997) Biophys. J. 72, 1109-1126)] of ca. 20, 60, and 90 pS, respectively, suggested that the multiple conductance levels seen experimentally must correspond to different degrees of oligomerization of the peptide to form channels.

The conduction of protons in different stereoisomers of dioxolane-linked gramicidin A channels

Two different stereoisomers of the dioxolane-linked gramicidin A (gA) channels were individually synthesized (the SS and RR dimers;. Science. 244:813-817). The structural differences between these dimers arise from different chiralities within the dioxolane linker. The SS dimer mimics the helicity and the inter- and intramolecular hydrogen bonding of the monomer-monomer association of gA's. In contrast, there is a significant disruption of the helicity and hydrogen bonding pattern of the ion channel in the RR dimer. Single ion channels formed by the SS and RR dimers in planar lipid bilayers have different proton transport properties. The lipid environment in which the different dimers are reconstituted also has significant effects on single-channel proton conductance (g(H)). g(H) in the SS dimer is about 2-4 times as large as in the RR. In phospholipid bilayers with 1 M [H(+)](bulk), the current-voltage (I-V) relationship of the SS dimer is sublinear. Under identical experimental conditions, the I-V plot of the RR dimer is supralinear (S-shaped). In glycerylmonooleate bilayers with 1 M [H(+)](bulk), both the SS and RR dimers have a supralinear I-V plot. Consistent with results previously published (. Biophys. J. 73:2489-2502), the SS dimer is stable in lipid bilayers and has fast closures. In contrast, the open state of the RR channel has closed states that can last a few seconds, and the channel eventually inactivates into a closed state in either phospholipid or glycerylmonooleate bilayers. It is concluded that the water dynamics inside the pore as related to proton wire transfer is significantly different in the RR and SS dimers. Different physical mechanisms that could account for this hypothesis are discussed. The gating of the synthetic gA dimers seems to depend on the conformation of the dioxolane link between gA's. The experimental results provide an important framework for a detailed investigation at the atomic level of proton conduction in different and relatively simple ion channel structures.

Refsum disease diagnostic marker phytanic acid alters the physical state of membrane proteins of liver mitochondria

Phytanic acid (3,7,11,15-tetramethylhexadecanoic acid), a branched chain fatty acid accumulating in Refsum disease to high levels throughout the body, induces uncoupling of rat liver mitochondria similar to non-branched fatty acids (e.g. palmitic acid), but the contribution of the ADP/ATP carrier or the aspartate/glutamate carrier in phytanic acid-induced uncoupling is of minor importance. Possible deleterious effects of phytanic acid on membrane-linked energy coupling processes were studied by ESR spectroscopy using rat liver mitochondria and a membrane preparation labeled with the lipid-specific spin probe 5-doxylstearic acid (5-DSA) or the protein-specific spin probe MAL-TEMPO (4-maleimido-2,2,6, 6-tetramethyl-piperidine-1-oxyl). The effects of phytanic acid on phospholipid molecular dynamics and on the physical state of membrane proteins were quantified by estimation of the order parameter or the ratio of the amplitudes of the weakly to strongly immobilized MAL-TEMPO binding sites (W/S ratio), respectively. It was found, that phytanic acid (1) increased the mobility of phospholipid molecules (indicated by a decrease in the order parameter) and (2) altered the conformational state and/or the segmental mobility of membrane proteins (indicated by a drastic decrease in the W/S ratio). Unsaturated fatty acids with multiple cis-double bonds (e.g. linolenic or arachidonic acid), but not non-branched FFA (ranging from chain length C10:0 to C18:0), also decrease the W/S ratio. It is hypothesized that the interaction of phytanic acid with transmembrane proteins might stimulate the proton permeability through the mitochondrial inner membrane according to a mechanism, different to a protein-supported fatty acid cycling.

Structure-function study on a de novo synthetic hydrophobic ion channel

Ion conduction properties of a de novo synthesized channel, formed from cyclic octa-peptides consisting of four alternate L-alanine (Ala) and N'-acylated 3-aminobenzoic acid (Aba) moieties, were studied in bilayer membranes. The single-channel conductance was 9 pS in symmetrical 500 mM KCl. The channel favored permeation of cations over anions with a permeability ratio (PCl-/PK+) of 0.15. The selectivity sequence among monovalent cations based on permeability ratio (PX+/PK+) fell into an order: NH4+(1.4) > Cs+(1. 1) >/= K+(1.0) > Na+(0.4) >> Li+(0). The conductance-activity relationship of the channel in K+ solutions followed simple Michaelis-Menten kinetics with a half-maximal saturating activity of 8 mM and a maximal conductance of 9 pS. The permeability ratio PNa+/PK+ remained constant ( approximately 0.40) under biionic concentrations from 10 to 500 mM. These results suggests that the channel is a one-ion channel. The pore diameter probed by a set of organic cations was approximately 6 A. The single-channel current was blocked by Ca2+ in a dose-dependent manner that followed a single-site titration curve with a voltage-dependent dissociation constant of 0.6 mM at 100 mV. The electric distance of the binding site for Ca2+ was 0.07 from both entrances of the channel, indicating the presence of two symmetrical binding sites in each vicinity of the channel entrance. Correlations between conduction properties and structural aspects of the channel are discussed in terms of a three-barrier and two-binding-site (3B2S) model of Eyring rate theory. All available structural information supported an idea that the channel was formed from a tail-to-tail associated dimer of the molecule, the pore of which was lined with hydrophobic acyl chains. This is the first report to have made a systematic analysis of ion permeation through a hydrophobic pore.

The arachidonate-activatable, NADPH oxidase-associated H+ channel is contained within the multi-membrane-spanning N-terminal region of gp91-phox

The generation of superoxide by the NADPH oxidase of neutrophils is accompanied by the efflux of H+ ions through a H+ channel. gp91-phox, a protein component of the oxidase, has been shown previously to function as a H+ channel [Henderson, Banting and Chappell (1995) J. Biol. Chem. 270, 5909-5916]. We have constructed a CHO cell line (CHO-N) that expresses an N-terminal fragment of gp91-phox containing the predicted multiple transmembrane domains of the protein. These cells exhibit H+ fluxes in response to an imposed proton motive force and in the presence of arachidonate (to open the channel). The H+ fluxes were indistinguishable from those observed in cells expressing full-length gp91-phox. Therefore the N-terminal 230 amino acids of gp91-phox contain all that is required to function as the NADPH oxidase-associated H+ channel.

Structure and dynamics of a proton wire: a theoretical study of H+ translocation along the single-file water chain in the gramicidin A channel

The rapid translocation of H+ along a chain of hydrogen-bonded water molecules, or proton wire, is thought to be an important mechanism for proton permeation through transmembrane channels. Computer simulations are used to study the properties of the proton wire formed by the single-file waters in the gramicidin A channel. The model includes the polypeptidic dimer, with 22 water molecules and one excess proton. The dissociation of the water molecules is taken into account by the "polarization model" of Stillinger and co-workers. The importance of quantum effects due to the light mass of the hydrogen nuclei is examined with the use of discretized Feynman path integral molecular dynamics simulations. Results show that the presence of an excess proton in the pore orients the single-file water molecules and affects the geometry of water-water hydrogen bonding interactions. Rather than a well-defined hydronium ion OH3+ in the single-file region, the protonated species is characterized by a strong hydrogen bond resembling that of O2H5+. The quantum dispersion of protons has a small but significant effect on the equilibrium structure of the hydrogen-bonded water chain. During classical trajectories, proton transfer between consecutive water molecules is a very fast spontaneous process that takes place in the subpicosecond time scale. The translocation along extended regions of the chain takes place neither via a totally concerted mechanism in which the donor-acceptor pattern would flip over the entire chain in a single step, nor via a succession of incoherent hops between well-defined intermediates. Rather, proton transfer in the wire is a semicollective process that results from the subtle interplay of rapid hydrogen-bond length fluctuations along the water chain. These rapid structural fluctuations of the protonated single file of waters around an average position and the slow movements of the average position of the excess proton along the channel axis occur on two very different time scales. Ultimately, it is the slow reorganization of hydrogen bonds between single-file water molecules and channel backbone carbonyl groups that, by affecting the connectivity and the dynamics of the single-file water chain, also limits the translocation of the proton across the pore.

Antifreeze glycoproteins inhibit leakage from liposomes during thermotropic phase transitions

Antifreeze glycoproteins (AFGPs), found in the blood of polar fish at concentrations as high as 35 g/liter, are known to prevent ice crystal growth and depress the freezing temperature of the blood. Previously, Rubinsky et al. [Rubinsky, B., Mattioli, M., Arav, A., Barboni, B. & Fletcher, G. L. (1992) Am. J. Physiol. 262, R542-R545] provided evidence that AFGPs block ion fluxes across membranes during cooling, an effect that they ascribed to interactions with ion channels. We investigated the effects of AFGPs on the leakage of a trapped marker from liposomes during chilling. As these liposomes are cooled through the transition temperature, they leak approximately 50% of their contents. Addition of less than 1 mg/ml of AFGP prevents up to 100% of this leakage, both during chilling and warming through the phase transition. This is a general effect that we show here applies to liposomes composed of phospholipids with transition temperatures ranging from 12 degrees C to 41 degrees C. Because these results were obtained with liposomes composed of phospholipids alone, we conclude that the stabilizing effects of AFGPs on intact cells during chilling reported by Rubinsky et al. may be due to a nonspecific effect on the lipid components of native membranes. There are other proteins that prevent leakage, but only under specialized conditions. For instance, antifreeze proteins, bovine serum albumin, and ovomucoid all either have no effect or actually induce leakage. Following precipitation with acetone, all three proteins inhibited leakage, although not to the extent seen with AFGPs. Alternatively, there are proteins such as ovotransferrin that have no effect on leakage, either before or after acetone precipitation.

Molecular dynamics simulations of water within models of ion channels

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.

Packing interactions of Aib-containing helices: molecular modeling of parallel dimers of simple hydrophobic helices and of alamethicin

alpha-Aminoisobutyric acid (Aib) is a helicogenic alpha, alpha-dimethyl amino acid found in channel-forming peptaibols such as alamethicin. Possible effects of Aib on helix-helix packing are analyzed. Simulated annealing via restrained molecular dynamics is used to generate ensembles of approximately parallel helix dimers. Analysis of variations in geometrical and energetic parameters within ensembles defines how tightly a pair of helices interact. Simple hydrophobic helix dimers are compared: Ala20, Leu20, Aib20, and P20, the latter a simple channel-forming peptide [G. Menestrina, K.P. Voges, G. Jung, and G. Boheim (1986) Journal of Membrane Biology, Vol. 93, pp. 111-132]. Ala20 and Leu20 dimers exhibit well-defined ridges-in-grooves packing with helix crossing angles (omega) of the order of +20 degrees. Aib20 alpha-helix dimers are much more loosely packed, as evidenced by a wide range of omega values and small helix-helix interaction energies. However, when in a 3(10) conformation Aib20 helices pack in three well-defined parallel modes, with omega ca. -15 degrees, +5 degrees, and 10 degrees. Comparison of helix-helix interaction energies suggests that dimerization may favor the 3(10) conformation. P20, with 8 Aib residues, also shows looser packing of alpha-helices. The results of these studies of hydrophobic helix dimers are analyzed in the context of the ridges-in-grooves packing model. Simulations are extended to dimers of alamethicin, and of an alamethicin derivative in which all Aib residues are replaced by Leu. This substitution has little effect on helix-helix packing. Rather, such interactions appear to be sensitive to interactions between polar side chains. Overall, the results suggest that Aib may modulate the packing of simple hydrophobic helices, in favor of looser interactions. For more complex amphipathic helices, interactions between polar side chains may be more critical.

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.

Synthetic proteins. Designer proteins

Attempts to synthesize de novo proteins with particular properties should improve our understanding of the long-range cooperative interactions that are so crucial to the structure and function of proteins.