Hydroxide and proton migration in aquaporins

Hypothetical hydroxide and proton migration along the linear water chain in Aquaporin GlpF from Escherichia coli are studied by ab initio Car-Parrinello molecular dynamics simulations. It is found that the protein stabilizes a bipolar single file of water. The single file features a contiguous set of water-water hydrogen bonds in which polarization of the water molecules vary with position along the channel axis. Deprotonation of the water chain promotes the reorientation of water molecules while the hydroxide ion rapidly migrates by sequentially accepting protons from the neighboring water molecules. The hydroxide ion is not attracted by a conserved, channel-lining arginine residue, but is immobilized at two centrally located, conserved Asparagine-Proline-Alanine motifs where fourfold coordination stabilizes the ion. Hydroxide transition from the channel vestibules into the channel lumen is strongly influenced by electrostatic coupling to two conserved oppositely aligned macrodipoles. This suggests that the macrodipole's negative poles play a role in preventing hydroxide ions from entering into the channel's inner vestibules. Water protonation within the lumen facilitates water reorientation and subsequent proton expelling occurs. In the periplasmic half-channel, expelling occurs via the Grotthuss mechanism. Protonation within the cytoplasmic half-channel implies wire-breakage at the Asn-Pro-Ala motifs. The proton is here diffusively rejected as (H(5)O(2))(+).     

Permeation of water through the KcsA K+ channel

Previous studies have reported that the KcsA potassium channel has an osmotic permeability coefficient of 4.8 x 10(-12) cm3/s, giving it a significantly higher osmotic permeability coefficient than that of some membrane channels specialized in water transport. This high osmotic permeability is proposed to occur when the channel is depleted of potassium ions, the presence of which slow down the water permeation process. The atomic structure of the potassium-depleted KcsA channel and the mechanisms of water permeation have not been well characterized so far. Here, all-atom molecular dynamics simulations, in conjunction with an umbrella sampling strategy and a nonequilibrium approach to simulate pressure gradients are employed to illustrate the permeation of water in the absence of ions through the KcsA K+ channel. Equilibrium molecular dynamics simulations (95 ns combined total length) identified a possible structure of the potassium-depleted KcsA channel, and umbrella sampling calculations (160 ns combined total length) revealed that this structure is not permeable by water molecules moving along the channel axis. The simulation of a pressure gradient across the channel (30 ns combined total length) identified an alternative permeation pathway with a computed osmotic permeability of approximately (2.7 +/- 0.9) x 10(-13) cm3/s. Water fluxes along this pathway did not proceed through collective water motions or transitions to vapor state. All of the major results of this study were robust against variations in a wide set of simulation parameters (force field, water model, membrane model, and channel conformation).     

Cluster analysis of consensus water sites in thrombin and trypsin shows conservation between serine proteases and contributions to ligand specificity.

Cluster analysis is presented as a technique for analyzing the conservation and chemistry of water sites from independent protein structures, and applied to thrombin, trypsin, and bovine pancreatic trypsin inhibitor (BPTI) to locate shared water sites, as well as those contributing to specificity. When several protein structures are superimposed, complete linkage cluster analysis provides an objective technique for resolving the continuum of overlaps between water sites into a set of maximally dense microclusters of overlapping water molecules, and also avoids reliance on any one structure as a reference. Water sites were clustered for ten superimposed thrombin structures, three trypsin structures, and four BPTI structures. For thrombin, 19% of the 708 microclusters, representing unique water sites, contained water molecules from at least half of the structures, and 4% contained waters from all 10. For trypsin, 77% of the 106 microclusters contained water sites from at least half of the structures, and 57% contained waters from all three. Water site conservation correlated with several environmental features: highly conserved microclusters generally had more protein atom neighbors, were in a more hydrophilic environment, made more hydrogen bonds to the protein, and were less mobile. There were significant overlaps between thrombin and trypsin conserved water sites, which did not localize to their similar active sites, but were concentrated in buried regions including the solvent channel surrounding the Na+ site in thrombin, which is associated with ligand selectivity. Cluster analysis also identified water sites conserved in thrombin but not trypsin, and vice versa, providing a list of water sites that may contribute to ligand discrimination. Thus, in addition to facilitating the analysis of water sites from multiple structures, cluster analysis provides a useful tool for distinguishing between conserved features within a protein family and those conferring specificity.     

The binding site of sodium in the gramicidin A channel: comparison of molecular dynamics with solid-state NMR data.

The location of the main binding site for sodium in the gramicidin A (GA) channel was investigated with molecular dynamics simulations, using an atomic model of the channel embedded in a fully hydrated dimyristoyl phosphatidycholine (DMPC) bilayer. Twenty-four separate simulations in which a sodium was restrained at different locations along the channel axis were generated. The results are compared with carbonyl 13C chemical shift anisotropy solid-state NMR experimental data previously obtained with oriented GA:DMPC samples. Predictions are made for other solid-state NMR properties that could be observed experimentally. The combined information from experiment and simulation strongly suggests that the main binding sites for sodium are near the channel's mouth, approximately 9.2 A from the center of the dimer channel. The 13C chemical shift anisotropy of Leu10 is the most affected by the presence of a sodium ion in the binding site. In the binding site, the sodium ion is lying off-axis, making contact with two carbonyl oxygens and two single-file water molecules. The main channel ligand is provided by the carbonyl group of the Leu10-Trp11 peptide linkage, which exhibits the largest deviation from the ion-free channel structure. Transient contacts with the carbonyl group of Val8 and Trp15 are also present. The influence of the tryptophan side chains on the channel conductance is examined based on the current information about the binding site.     

Lipoelectric modification of ion channel voltage gating by polyunsaturated fatty acids

Polyunsaturated fatty acids (PUFAs) have beneficial effects on epileptic seizures and cardiac arrhythmia. We report that ω-3 and ω-6 all-cis-PUFAs affected the voltage dependence of the Shaker K channel by shifting the conductance versus voltage and the gating charge versus voltage curves in negative direction along the voltage axis. Uncharged methyl esters of the PUFAs did not affect the voltage dependence, whereas changes of pH and charge mutations on the channel surface affected the size of the shifts. This suggests an electrostatic effect on the channel's voltage sensors. Monounsaturated and saturated fatty acids, as well as trans-PUFAs did not affect the voltage dependence. This suggests that fatty acid tails with two or more cis double bonds are required to place the negative carboxylate charge of the PUFA in a position to affect the channel's voltage dependence. We propose that charged lipophilic compounds could play a role in regulating neuronal excitability by electrostatically affecting the channel's voltage sensor. We believe this provides a new approach for pharmacological treatment that is voltage sensor pharmacology.     

Ion passage pathways and thermodynamics of the amphotericin B membrane channel

Amphotericin B is a polyene macrolide antibiotic used to treat systemic fungal infections. Amphotericin B's chemotherapeutic action requires the formation of transmembrane channels, which are known to transmit monovalent ions. We have investigated the ion passage pathways through the pore of a realistic model structure of the channel and computed the associated thermodynamic properties. Our calculations combined the free energy computations using the Poisson equation with a continuum solvent model and the molecular simulations in which solvent molecules were present explicitly. It was found that there are no substantial structural barriers to a single sodium or chloride ion passage. Thermodynamic free energy calculations showed that the path along which the ions prefer to move is off center from the channel's central axis. In accordance with experiments, Monte Carlo molecular simulations established that sodium ions can pass through the pore. When it encounters a chloride anion in the channel, the sodium cation prefers to form a solvent-bridged pair configuration with the anion.     

Time-correlation analysis of simulated water motion in flexible and rigid gramicidin channels.

Molecular dynamics simulations have been done on a system consisting of the polypeptide membrane channel former gramicidin, plus water molecules in the channel and caps of waters at the two ends of the channel. In the absence of explicit simulation of the surrounding membrane, the helical form of the channel was maintained by artificial restraints on the peptide motion. The characteristic time constant of the artificial restraint was varied to assess the effect of the restraints on the channel structure and water motions. Time-correlation analysis was done on the motions of individual channel waters and on the motions of the center of mass of the channel waters. It is found that individual water molecules confined in the channel execute higher frequency motions than bulk water, for all degrees of channel peptide restraint. The center-of-mass motion of the chain of channel waters (which is the motion that is critical for transmembrane transport, due to the mandatory single filing of water in the channel) does not exhibit these higher frequency motions. The mobility of the water chain is dramatically reduced by holding the channel rigid. Thus permeation through the channel is not like flow through a rigid pipe; rather permeation is facilitated by peptide motion. For the looser restraints we used, the mobility of the water chain was not very much affected by the degree of restraint. Depending on which set of experiments is considered, the computed mobility of our water chain in the flexible channel is four to twenty times too high to account for the experimentally measured resistance of the gramicidin channel to water flow. From this result it appears likely that the peptide motions of an actual gramicidin channel embedded in a lipid membrane may be more restrained than in our flexible channel model, and that these restraints may be a significant modulator of channel permeability. For the completely rigid channel model the "trapping" of the water molecules in preferred positions throughout the molecular dynamics run precludes a reasonable assessment of mobility, but it seems to be quite low.     

Calibrated measurement of gating-charge arginine displacement in the KvAP voltage-dependent K+ channel

Voltage-dependent ion channels open and conduct ions in response to changes in cell-membrane voltage. The voltage sensitivity of these channels arises from the motion of charged arginine residues located on the S4 helices of the channel's voltage sensors. In KvAP, a prokaryotic voltage-dependent K+ channel, the S4 helix forms part of a helical hairpin structure, the voltage-sensor paddle. We have measured the membrane depth of residues throughout the KvAP channel using avidin accessibility to different-length tethered biotin reagents. From these measurements, we have calibrated the tether lengths and derived the thickness of the membrane that forms a barrier to avidin penetration, allowing us to determine the magnitude of displacement of the voltage-sensor paddles during channel gating. Here we show that the voltage-sensor paddles are highly mobile compared to other regions of the channel and transfer the gating-charge arginines 15-20 A through the membrane to open the pore.     

Water and polypeptide conformations in the gramicidin channel. A molecular dynamics study.

Theoretical studies of ion channels address several important questions. The mechanism of ion transport, the role of water structure, the fluctuations of the protein channel itself, and the influence of structural changes are accessible from these studies. In this paper, we have carried out a 70-ps molecular dynamics simulation on a model structure of gramicidin A with channel waters. The backbone of the protein has been analyzed with respect to the orientation of the carbonyl and the amide groups. The results are in conformity with the experimental NMR data. The structure of water and the hydrogen bonding network are also investigated. It is found that the water molecules inside the channel act as a collective chain; whereas the conformation in which all the waters are oriented with the dipoles pointing along the axis of the channel is a preferred one, others are also accessed during the dynamics simulation. A collective coordinate involving the channel waters and some of the hydrogen bonding peptide partners is required to describe the transition of waters from one configuration to the other.     

Translocation and interactions of L-arabinose in OmpF porin: A molecular dynamics study

The passage of a natural substrate, L-arabinose (L-ARA) through Escherichia coli porin embedded in an artificial bilayer, is studied by equilibrium molecular dynamics simulations. We investigate the early stage of translocation process of L-ARA from intra-cellular to extra-cellular side (Int-to-Ext) across the bilayer. The average trajectory path over all L-ARA molecules along with quantum-mechanical configuration-optimizations at PM3 level predict the existence of at least three trapping zones. The common feature within all these zones is that L-ARA remains perpendicular to the channel axis. It is remarkable how the orientation and translational-rotational motion of L-ARA molecule play a role in its transport through OmpF channel. These simulations are important for better understanding of permeation process in OmpF channel. They also provide an insight into the chiral recognition of translocation process in protein nanochannels from substrate and protein prospects and help interpret experiments on permeation process of small dipolar molecules across biological membranes.     

Crystal and molecular structure of the A-DNA dodecamer d(CCGTACGTACGG). Choice of fragment helical axis.

The crystal structure of the dodecamer d(CCGTACGTACGG) has been determined at 2.5 A resolution. The crystals grow in the hexagonal space group P6(1)22, a = b = 46.2 A, c = 71.5 A with one strand as the asymmetric unit. Diffraction data were collected by the oscillation film method yielding 1664 unique reflections with an Rmerge of 0.04. The structure was solved by real-space rotational translational searches with idealized helical models of A, B and Z-DNA. The best agreement was given by an A-DNA model with its dyad axis along the diagonal crystallographic dyad axis, with an R-factor 0.43 and correlation coefficient of 0.59 for data between 10 and 5 A. Iterative map fitting and restrained least-squares refinement and addition of 40 solvent molecules brought the R-factor to 0.15 and the correlation coefficient to 0.97 for all data between 8.0 and 2.5 A. The stereochemistry of the atomic model is good, with a root-mean-square deviation in bond distances of 0.006 A. This is the first example of an A-DNA containing a full helical turn. The dodecamer displays a novel packing motif. In addition to the characteristic contacts between the terminal base-pairs and the minor grooves of symmetry-related molecules, there are also minor groove to minor groove interactions not previously observed. The packing leaves an approximately 25 A diameter solvent channel around the origin, along the c-axis. The presence of a prominent 3.4 A meridional reflection and other diffuse features in the diffraction pattern provided evidence for the presence of disordered B-DNA along the c-axis, which can be accommodated in these solvent channels. The molecular conformation of the dodecamer also displays novel features. The dyad-related halves of the molecule are bent at an angle of 20 degrees, and the helical parameters are affected by this bend. Unlike the shorter A-DNA octamers, the dimensions of the major groove can be directly measured. Novel correlations between local helical parameters and global conformational features are presented. Most of the solvent molecules are associated with the major groove and the sugar-phosphate backbone.     

Crystal structure of the low-humidity form of aspartame sweetener.

The low-humidity IB crystal form of aspartame (L-alphaaspartyl-L-phenylalanine methyl ester) is prepared via humidity-induced transition from the highly hydrated IA crystal form and is used widely as a sweetener. The crystal structure of the low-humidity IB form is determined at 1.05 A resolution (0.476 A(-1) in maximum sintheta/lambda) from an extremely fine fibrous crystal using synchrotron radiation. There are three aspartame molecules and two water molecules in the asymmetric unit of the monoclinic space group P2(1). Each aspartame molecule adopts an almost identical extended conformation which is commonly observed in other crystal forms of aspartame. Three aspartame molecules are assembled into a triangular trimer, and trimer units are stacked along the b-axis via hydrogen-bonding and electrostatic interactions in the main chains and also via hydrophobic contacts in the phenyl side-chains. Six trimer units are related by pseudo 6(1)-screw axis symmetry and form a hydrophilic channel at their center. The hydrophilic channel in the IB form contains only four water molecules in the unit cell, compared with 16 in the IA form. Although the IB form exhibits a trimer structure similar to that of the IA form, one aspartame molecule is rotated by approximately equals 20 degrees from the orientation in the IA form. This arrangement of the molecule implies that the humidity-induced transition is accompanied by a flapping motion of its methyl ester group. These structural differences may imply the stepwise transition from the IA to the IB forms.     

Side-chain dynamics are critical for water permeation through aquaporin-1

Molecular dynamics simulations of aquaporin-1 embedded in a solvated lipid bilayer were carried out to investigate the mechanism of water permeation. The 2.2 Å resolution crystal structure of the bovine protein was used for five independent trajectories. During the equilibration and preparatory steps in which the protein was held fixed, water molecules inside the water channel adopted the same positions as observed in the crystal structure but they did not pass through the channel, suggesting that the dynamic motion of the protein is critical for water permeation. When the protein atoms were allowed to move, the side chains of the two asparagines in the two conserved Asn-Pro-Ala motifs near the center of the channel formed hydrogen bonds with water and helped water molecules move through the channel by actively aligning them for transport. The main-chain oxygen atoms, which were exposed to the pore surface in the crystal structure, also contributed to water transfer. Besides the constriction region observed in the crystal structure (Arg¹⁹⁷, Phe⁵⁸, His¹⁸², and Cys¹⁹¹), we found that His⁷⁶ and Val¹⁵⁵ act as a valve by dynamically blocking water permeation and helping control flow.     

Batrachotoxin-activated Na+ channels in planar lipid bilayers. Competition of tetrodotoxin block by Na+

Single Na+ channels from rat skeletal muscle plasma membrane vesicles were inserted into planar lipid bilayers formed from neutral phospholipids and were observed in the presence of batrachotoxin. The batrachotoxin-modified channel activates in the voltage range -120 to - 80 mV and remains open almost all the time at voltages positive to -60 mV. Low levels of tetrodotoxin (TTX) induce slow fluctuations of channel current, which represent the binding and dissociation of single TTX molecules to single channels. The rates of association and dissociation of TTX are both voltage dependent, and the association rate is competitively inhibited by Na+. This inhibition is observed only when Na+ is increased on the TTX binding side of the channel. The results suggest that the TTX receptor site is located at the channel's outer mouth, and that the Na+ competition site is not located deeply within the channel's conduction pathway.     

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.     

A MECHANISTIC VIEW OF THE NON-IDEAL OSMOTIC AND MOTIONAL BEHAVIOR OF INTRACELLULAR WATER

It is commonly assumed that essentially all of the water in cells has the same ideal motional and colligative properties as does water in bulk liquid state. This assumption is used in studies of volume regulation, transmembrane movement of solutes and electrical potentials, solute and solution motion, solute solubility and other phenomena. To get at the extent and the source of non-ideally behaved water (an operational term dependent on the measurement method), we studied the motional and colligative properties of water in cells, in solutions of amino acids and glycine peptides whose surface characteristics are known, and in solution of bovine serum albumin, hemoglobin and some synthetic polypeptides. Solutions of individual amino acids with progressively larger hydrophobic side chains showed one perturbed water molecule (structured-slowed in motion) per nine square angstroms of hydrophobic surface area. Water molecules adjacent to hydrophobic surfaces form pentagonal structural arrays, as shown by X-ray diffraction studies, that are reported to be disrupted by heat, electric field, hydrostatic pressure and phosphorylation state. Hydrophilic amino acids demonstrated water destructuring (increased motion) that was attributed to dielectric realignment of dipolar water molecules in the electric field between charge groups. In solutions of proteins, several methods indicate the equivalent of 2–8 layers of structured water molecules extending beyond the protein surface, and we have recently demonstrated that induced protein conformational change modifies the extent of non-ideally behaved water. Water self-diffusion rate as measured in three different cell types was about half that of bulk water, indicating that most of the water in these cells was slower in motion than bulk water. In different cell types the extent of osmotically perturbed water ranged from less than half to almost all of the intracellular water. The assumption that essentially all intracellular water has ideal osmotic and motional behavior is not supported by the experimental findings. The non-ideality of cell water is an operational term. Therefore, the amount of non-ideally behaving water is dependent on the characteristics of water targeted, i.e. the measurement method, and a large fraction of it is explainable in mechanistic terms at a molecular level based on solute—solvent interactions.     

Water transport in human aquaporin-4: molecular dynamics (MD) simulations

Aquaporin-4 (AQP4) is the predominant water channel in the central nervous system, where it has been reported to be involved in many pathophysiological roles including water transport. In this paper, the AQP4 tetramer was modeled from its PDB structure file, embedded in a palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer, solvated in water, then minimized and equilibrated by means of molecular dynamics simulations. Analysis of the equilibrated structure showed that the central pore along the fourfold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-195, Leu-191 and Leu-75, form the narrowest part of the pore. Therefore water molecules are not expected to transport through the central pore, which was confirmed by MD simulations. Each monomer of the AQP4 tetramers forms a channel whose walls consist mostly of hydrophilic residues. There are eight water molecules in single file observed in each of the four channels, transporting through the selectivity filter containing Arg-216, His-201, Phe-77, Ala-210, and the two conserved Asn-Pro-Ala (NPA) motifs containing Asn-213 and Asn-97. By using Brownian dynamics fluctuation–dissipation-theorem (BD-FDT), the overall free-energy profile was obtained for water transporting through AQP4 for the first time, which gives a complete map of the entire channel of water permeation.     

The open state gating mechanism of gramicidin a requires relative opposed monomer rotation and simultaneous lateral displacement

The gating mechanism of the open state of the gramicidin A (gA) channel is studied by using a new Monte Carlo Normal Mode Following (MC-NMF) technique, one applicable even without a target structure. The results demonstrate that the lowest-frequency normal mode (NM) at approximately 6.5 cm(-1) is the crucial mode that initiates dissociation. Perturbing the gA dimer in either direction along this NM leads to opposed, nearly rigid-body rotations of the gA monomers around the central pore axis. Tracking this NM by using the eigenvector-following technique reveals the channel's gating mechanism: dissociation via relative opposed monomer rotation and simultaneous lateral displacement. System evolution along the lowest-frequency eigenvector shows that the large-amplitude motions required for gating (dissociation) are not simple relative rigid-body motions of the monomers. Gating involves coupling intermonomer hydrogen bond breaking, backbone realignment, and relative monomer tilt with complex side chain reorganization at the intermonomer junction.     

The location of a closed channel gate in the GABAA receptor channel

Considerable controversy surrounds the location of the closed channel gate in members of the Cys-loop receptor family of neurotransmitter-gated ion channels that includes the GABAA, glycine, acetylcholine, and 5-HT3 receptors. Cysteine-accessibility studies concluded that the gate is near the cytoplasmic end of the channel in acetylcholine and GABAA receptors but in the middle of the 5-HT3A receptor channel. Zn2+ accessibility studies in a chimeric 5-HT3-ACh receptor suggested the gate is near the channel's cytoplasmic end. In the 4-A resolution structure of the acetylcholine receptor closed state determined by cryoelectron microscopy, the narrowest region, inferred to be the gate, is in the channel's midsection from 9' to 14' but the M1-M2 loop residues at the channel's cytoplasmic end were not resolved in that structure. We used blocker trapping experiments with picrotoxin, a GABAA receptor open channel blocker, to determine whether a gate exists at a position more extracellular than the picrotoxin binding site, which is in the vicinity of alpha1Val257 (2') near the channel's cytoplasmic end. We show that picrotoxin can be trapped in the channel after removal of GABA. By using the state-dependent accessibility of engineered cysteines as reporters for the channel's structural state we infer that after GABA washout, with picrotoxin trapped in the channel, the channel appears to be in the closed state. We infer that a gate exists between the picrotoxin binding site and the channel's extracellular end, consistent with a closed channel gate in the middle of the channel. Given the homology with acetylcholine and 5-HT3 receptors there is probably a similar gate in those channels as well. This does not preclude the existence of an additional gate at a more cytoplasmic location.