Unraveling the mechanism of proton translocation in the extracellular half‐channel of bacteriorhodopsin

Bacteriorhodopsin, a light activated protein that creates a proton gradient in halobacteria, has long served as a simple model of proton pumps. Within bacteriorhodopsin, several key sites undergo protonation changes during the photocycle, moving protons from the higher pH cytoplasm to the lower pH extracellular side. The mechanism underlying the long‐range proton translocation between the central (the retinal Schiff base SB216, D85, and D212) and exit clusters (E194 and E204) remains elusive. To obtain a dynamic view of the key factors controlling proton translocation, a systematic study using molecular dynamics simulation was performed for eight bacteriorhodopsin models varying in retinal isomer and protonation states of the SB216, D85, D212, and E204. The side‐chain orientation of R82 is determined primarily by the protonation states of the residues in the EC. The side‐chain reorientation of R82 modulates the hydrogen‐bond network and consequently possible pathways of proton transfer. Quantum mechanical intrinsic reaction coordinate calculations of proton‐transfer in the methyl guanidinium‐hydronium‐hydroxide model system show that proton transfer via a guanidinium group requires an initial geometry permitting proton donation and acceptance by the same amine. In all the bacteriorhodopsin models, R82 can form proton wires with both the CC and the EC connected by the same amine. Alternatively, rare proton wires for proton transfer from the CC to the EC without involving R82 were found in an O′ state where the proton on D85 is transferred to D212. Proteins 2016; 84:639–654. © 2016 Wiley Periodicals, Inc.     

The mechanism of proton exclusion in aquaporin channels

The mechanism of proton exclusion in aquaporin channels is elucidated through free energy calculations of the pathway of proton transport. The second generation multistate empirical valence bond (MS-EVB2) model was applied to simulate the interaction of an excess proton with the channel environment. Jarzynski's equality was employed for rapid convergence of the free energy profile. A barrier sufficiently high to block proton transport is located near the channel center at the NPA motif-a site involved in bi-orientational ordering of the embedded water-wire in absence of the excess proton. A second and lower barrier is observed at the selectivity filter near the periplasmic outlet where the channel is narrowest. This secondary barrier may be essential in filtering other large solutes and cations.     

Enhancement of proton conductance by mutations of the selectivity filter of aquaporin-1

Prevention of cation permeation in wild-type aquaporin-1 (AQP1) is believed to be associated with the Asn-Pro-Ala (NPA) region and the aromatic/arginine selectivity filter (SF) domain. Previous work has suggested that the NPA region helps to impede proton permeation due to the protein backbone collective macrodipoles that create an environment favoring a directionally discontinuous channel hydrogen-bonded water chain and a large electrostatic barrier. The SF domain contributes to the proton permeation barrier by a spatial restriction mechanism and direct electrostatic interactions. To further explore these various effects, the free-energy barriers and the maximum cation conductance for the permeation of various cations through the AQP1-R195V and AQP1-R195S mutants are predicted computationally. The cations studied included the hydrated excess proton that utilizes the Grotthuss shuttling mechanism, a model “classical” charge localized hydronium cation that exhibits no Grotthuss shuttling, and a sodium cation. The hydrated excess proton was simulated using a specialized multi-state molecular dynamics method including a proper physical treatment of the proton shuttling and charge defect delocalization. Both AQP1 mutants exhibit a surprising cooperative effect leading to a reduction in the free-energy barrier for proton permeation around the NPA region due to altered water configurations in the SF region, with AQP1-R195S having a higher conductance than AQP1-R195V. The theoretical predictions are experimentally confirmed in wild-type AQP1 and the mutants expressed in Xenopus oocytes. The combined results suggest that the SF domain is a specialized structure that has evolved to impede proton permeation in aquaporins.     

The structural basis of water permeation and proton exclusion in aquaporins (Review)

Aquaporins (AQPs) represent a ubiquitous class of integral membrane proteins that play critical roles in cellular osmoregulations in microbes, plants and mammals. AQPs primarily function as water-conducting channels, whereas members of a sub-class of AQPs, termed aquaglyceroporins, are permeable to small neutral solutes such as glycerol. While AQPs facilitate transmembrane permeation of water and/or small neutral solutes, they preclude the conduction of protons. Consequently, openings of AQP channels allow rapid water diffusion down an osmotic gradient without dissipating electrochemical potentials. Molecular structures of AQPs portray unique features that define the two central functions of AQP channels: effective water permeation and strict proton exclusion. This review describes AQP structures known to date and discusses the mechanisms underlying water permeation, proton exclusion and water permeability regulation.     

Three ways in, one way out: water dynamics in the trans-membrane domains of the inner membrane translocase AcrB

Powered by proton-motive force, the inner membrane translocase AcrB is the engine of the AcrAB-TolC efflux pump in Escherichia coli. As proton conduction in proteins occurs along hydrogen-bonded networks of polar residues and water molecules, knowledge of the protein-internal water distribution and water-interacting residues allows drawing conclusions to possible pathways of proton conduction. Here, we report a series of 6× 50 ns independent molecular dynamics simulations of asymmetric AcrB embedded in a phospholipid/water environment. Simulating each monomer in its proposed protonation state, we calculated for each trans-membrane domain the average water distribution, identified residues interacting with these waters and quantified each residue's frequency of water hydrogen bond contact. Combining this information we find three possible routes of proton transfer connecting a continuously hydrated region of known key residues in the TMD interior to bulk water by one cytoplasmic and up to three periplasm water channels in monomer B and A. We find that water access of the trans-membrane domains is regulated by four groups of residues in a combination of side chain re-orientations and shifts of trans-membrane helices. Our findings support a proton release event via Arg971 during the C intermediate or in the transition to A, and proton uptake occurring in the A or B state or during a so far unknown intermediate in between B and C where cytoplasmic water access is still possible. Our simulations suggest experimentally testable hypotheses, which have not been investigated so far.     

Proton transfer along the external proton channel of bacteriorhodopsin

The process of proton transfer along a proton channel is considered using bacteriorhodopsin as a model system, for which a large body of experimental data is available. The possible amino acid composition of the external proton half-channel of bacteriorhodopsin and the stepwise scheme of proton transfer consistent with experimental data are proposed. The rate of proton transfer between fixed centers is assessed for certain regions of this channel for which spectroscopic data are available.     

Molecular dynamics simulations investigate the mechanism of Psalmotoxin 1 regulating gating process of an acid-sensing ion channel 1a at pH 5.5

Acid-sensing ion channel 1a (ASIC1a) is a cation channel activated by protons and causes neuronal death through central nervous system. Psalmotoxin1 (PcTx1) is a gating modifier for ASIC1a. The process of PcTx1 regulating the channel gating from the extracellular domain to the transmembrane domain is unclear. Here we used molecular dynamics (MD) simulations method to investigate how PcTx1 regulates the gating of the ASIC1a. Our results indicated that PcTx1can mainly regulate ASIC1a gating process through hydrogen bonds, which can affect their relative positions of several key domains in ASIC1a, further, a long-range conformational changes path was determined, which is composed of β1, β2, β10, α6, α7, β11, and β12 in ASIC1a.     

Trafficking mechanism of water channel aquaporin-2

Targeted positioning of the water channel AQP2 (aquaporin-2) strictly regulates body water homoeostasis. Trafficking of AQP2 to the apical membrane is critical for the reabsorption of water in renal collecting ducts. In addition to the cAMP-mediated effect of vasopressin on AQP2 trafficking to the apical membrane, other signalling cascades can also induce this sorting. Recently, AQP2-binding proteins which could regulate this trafficking have been discovered; SPA-1 (signal-induced proliferation-associated gene-1), a GAP (GTPase-activating protein) for Rap1, and the cytoskeletal protein actin. This review summarizes recent advances related to the trafficking mechanisms of AQP2.     

The mechanism of proton transfer between adjacent sites on the molecular surface

The surface of a protein, or a membrane, is spotted with a multitude of proton binding sites, some of which are only few Å apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Eventually, the released protons are dispersed in the bulk, but during the first few nanoseconds after the dissociation, the protons can be trapped by encounter with nearby acceptor sites. While the study of this reaction on the surface of a protein suffers from experimental and theoretical difficulties, it can be investigated with simple model compounds like derivatives of fluorescein. In the present study, we evaluate the mechanism of proton transfer reactions that proceed, preferentially, inside the Coulomb cage of the dye molecules. Kinetic analysis of the measured dynamics reveals the role of the dimension of the Coulomb cage on the efficiency of the reaction and how the ordering of the water molecules by the dye affects the kinetic isotope effect.     

Hypertonic induction of aquaporin-1 water channel independent of transcellular osmotic gradient

Aquaporin-1 (AQP1) water channel plays a critical role for water reabsorption in the urinary concentrating mechanism. AQP1 expression in renal cells is upregulated by hypertonicity, but not urea, suggesting the requirement of an osmotic gradient. To investigate whether AQP1 expression is regulated by apical and/or basolateral hypertonicity, murine renal medullary mIMCD-K2 cells grown on permeable support were exposed to hypertonic medium. When the medium on the apical or basolateral membrane side was switched to hypertonic, the transcellular osmotic gradient was dissipated within 8h. Basolateral hypertonicity increased AQP1 expression more than apical hypertonicity. Comparable apical and basolateral hypertonicity without a transcellular hypertonic gradient, however, increased AQP1 expression. Cell surface biotinylation experiments revealed that hypertonicity promoted AQP1 trafficking to both plasma cell membranes. These results indicate that AQP1 expression is predominantly mediated by basolateral hypertonicity but a transcellular osmotic gradient is not necessary for its induction.     

Ion exclusion mechanism in aquaporin at an atomistic level

Although the mechanism of proton exclusion in aquaporin is investigated by many researchers, the detailed molecular mechanism for ion exclusion in aquaporin is still not completely understood. In the present work, a detailed mechanism for ion exclusion in aquaporin-1 (AQP1) at an atomistic level is investigated by calculating the free energy for transport of ions in AQP1 using an atomistic molecular dynamics simulation. For this purpose, sodium and chloride ions are chosen as representatives for nonprotonic ions. The simulation shows that the free energy barrier showing its maximum is located at the NPA region for sodium ion while it is located at both the front and the rear for chloride ion and that the barrier height is 18 and 9 kcal/mol, respectively, indicating that the ions are not able to pass through aquaporin. Analysis of the pair interaction energy between the permeating ion and its environment reveals that sodium ion is excluded by the positive charge generated by two alpha-helical macro-dipoles, while chloride ion is expelled by carbonyl oxygen atoms protruding from pore-making residues before it reaches the NPA motif. It is also found that the number of water molecules hydrating the ions is reduced as the ions enter the pore, implying that the energetic cost for detaching water molecules from a permeating ion also contributes to the free energy barriers of ion transport in AQP1.     

Molecular dynamics analysis of potent inhibitors of M2 proton channel against H1N1 swine influenza virus

The recent H1N1 (swine) influenza pandemic highlighted the urgent need of having effective anti‐viral strategies. In addition to neuraminidase inhibitors, there is another class of anti-viral drug known as M2 inhibitors that were, in the past, effective in treating seasonal influenza. However, due to the emergence of M2 inhibitor‐resistant influenza viruses, this class of drugs was not recommended for clinical usage in the latest influenza pandemic. In order to identify novel M2 inhibitors, we have performed molecular docking using a traditional Chinese medicine database (http://tcm.cmu.edu.tw/index.php). Docking and subsequent de novo designs gave 10 derivatives that have much better docking results than the control. Of these 10 derivatives, the top three, methyl isoferulate_1, genipin_1 and genipin_2, were selected for molecular dynamics simulation. During the simulation run, the top three derivatives all had stable interactions with M2 residues, Ser31 and Ala30. Methyl isoferulate_1 also has stable interaction to His37. Therefore, we recommend these three derivatives for further biomolecular experiments and clinical studies.     

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).     

The mechanism of proton translocation in respiratory complex I from molecular dynamics

Abstract

Respiratory complex I, the biggest enzyme of respiratory chain, plays a key role in energy production by the mitochondrial respiratory chain and has been implicated in many human neurodegenerative diseases. Recently, the crystal structure of respiratory complex I is reported. We perform 50?ns molecular dynamics simulations on the membrane domain of respiratory complex I under two hypothetical states (oxidized state and reduced state). We find that the density of water molecules in the trans-membrane domain under reduced state is bigger than that under oxidized state. The connecting elements (helix HL and β-hairpins-helix element) fluctuate stronger under reduced state than that under oxidized state, causing more internal water molecules and facilitating the proton conduction. The conformational changes of helix HL and the crucial charged residue Glu in TM5 play key roles in the mechanism of proton translocation. Our results illustrate the dynamic behavior and the potential mechanism of respiratory complex I, which provides the structural basis for drug design of respiratory complex I.     

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.     

Potentials of mean force in acidic proton transfer reactions in constrained geometries

Free energy barriers associated with the transfer of an excess proton in water and related to the potentials of mean force in proton transfer episodes have been computed in a wide range of thermodynamic states, from low-density amorphous ices to high-temperature liquids under the critical point for unconstrained and constrained systems. The latter were represented by set-ups placed inside hydrophobic graphene slabs at the nanometric scale allocating a few water layers, namely one or two in the narrowest case. Water–proton and carbon–proton forces were modelled with a Multi-State Empirical Valence Bond method. As a general trend, a competition between the effects of confinement and temperature is observed on the local hydrogen-bonded structures around the lone proton and, consequently, on the mean force exerted by its environment on the water molecule carrying the proton. Free energy barriers estimated from the computed potentials of mean force tend to rise with the combined effect of increasing temperatures and the packing effect due to a larger extent of hydrophobic confinement. The main reason observed for such enhancement of the free energy barriers was the breaking of the second coordination shell around the lone proton.     

A QM/MM study of proton transport pathways in a [NiFe] hydrogenase

A theoretical QM/MM study of the [NiFe] hydrogenase from Desulfovibrio fructosovorans has been performed to investigate possible routes of proton transfer between the active site and the protein surface. We obtained the minimum energy paths, with a modified version of the nudged elastic band method, for a set of proposed pathways. The calculations were carried out for the crystallographic structure and for several structures of the protein obtained from a molecular dynamics simulation. The results show one of the studied pathways to be preferred for transport from the active site to the surface, but the preference is not so strong when transport occurs in the opposite direction.     

The proton pump/leak mechanism of unconsciousness

There are few explanations which account for the manner in which the catastrophic physiological consequences of anaesthesia, cold narcosis or, for that matter, a short, sharp upper-cut, come about. Most studies terminate with the presentation of ever-better correlations between an end-point in a model system (dough consistency, rubber elasticity, bacterial, protozoal or animal mobility, liposome permeability, luciferase activity, etc.) and oil/water partition co-efficients or with some arbitrary biological end-point. From what is currently known about (a) the permeating pathways of non-electrolytes, ions and protons across membranes e.g. liposomes, (b) the effect of anaesthetics on such pathways and (c) the effect of temperature and pressure on both liposomes and whole animals, it is possible to develop a testable hypothesis. It is called the ‘proton pump-leak’ hypothesis and involves a number of linked biophysical and biochemical processes. (d) It assumes that a living animal or plant is in a steady-state regarding all concentration gradients; passive leaks across membranes are balanced by temperature, pressure, and energy dependent ion/ion and/or proton/ion pumps (enzyme), working within an aqueous phase. (e) Consciousness is dependent upon inter-neuronal communication via release of transmitter substances. (f) Transmitter substances, characteristically either weak bases or weak acids e.g. catecholamines, accumulate passively in vesicles rich in acid-buffer, held to a low pH by the activity of H⁺/K⁺ energy-driven pumps. Interference with this finely-balanced system either by changing the chemical potential of the hydrophobic (membrane) phase at NTP (with anaesthetics), or by changing the chemical potential of both hydrophobic and aqueous (pump) phases by hyperbaric, hyphotermic, or anoxic conditions imposed (inevitably) on the whole animal, would result in the resetting of the steady-state parameters.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

The influences of protonation state of histidine on aromatic/arginine region of aquaporin-1 protein

Aquaporin-1 (AQP1) is widely distributed in the epithelial tissue for water absorption and secretion. The histidine (His182) residue in the aromatic/arginine (ar/R) constriction region plays an important role in transporting water through the membrane. In this study, we have performed a total of 46 ns equilibrium molecular dynamics (MD) simulations, and obtained the influence of His182 in two protonation states (Hsd is the proton at N_δ and Hse is the proton at N_?) on the ar/R region. Water permeation rate shows that it is easier for water molecules to permeate the ar/R region of the AQP1 with residue in the Hsd state than in the Hse state. The minimum radii of the pore in the ar/R region were calculated during the last 10 ns MD simulation. We have analysed the correlation among the state of the pore (open or close), the minimum radius of the ar/R region and the dihedral angles < C_β-C_γ-C_δ-N_? of Arg197. The results show that the minimum radius can be used to mark the state of the pore.