Ab initio molecular dynamics study of proton transfer in a polyglycine analog of the ion channel gramicidin A.

Proton transfer in biological systems is thought to often proceed through hydrogen-bonded chains of water molecules. The ion channel, gramicidin A (gA), houses within its helical structure just such a chain. Using the density functional theory based ab initio molecular dynamics Car-Parrinello method, the structure and dynamics of proton diffusion through a polyglycine analog of the gA ion channel has been investigated. In the channel, a proton, which is initially present as hydronium (H3O+), rapidly forms a strong hydrogen bond with a nearest neighbor water, yielding a transient H5O2+ complex. As in bulk water, strong hydrogen bonding of this complex to a second neighbor solvation shell is required for proton transfer to occur. Within gA, this second neighbor shell included not only a channel water molecule but also a carbonyl of the channel backbone. The present calculations suggest a transport mechanism in which a priori carbonyl solvation is a requirement for proton transfer.     

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

The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.     

A molecular dynamics study of water chain formation in the proton-conducting K channel of cytochrome c oxidase

The formation of water chains in cytochrome c oxidase (CcO) is studied by molecular dynamics (MD). Focus is on water chains in the K channel that can supply a proton to the binuclear center (the heme a3 Fe/CuB region), the site of O2 reduction. By assessing the presence of chains of any length on a short time scale (0.1 ps), a view of the kinds of chains and their persistence is obtained. Chains from the entry of the channel on the inner membrane to Thr359 (Rhodobacter sphaeroides numbering) are often present but are blocked at that point until a rotation of the Thr359 side chain occurs, permitting formation of chains from Thr359 towards the binuclear center. No continuous hydrogen-bonded water chains are found connecting Thr359 and the binuclear center. Instead, waters hydrogen bond from Thr359 to the hydroxyl of the heme a3 farnesyl and then continue to the binuclear center via Tyr288, which has been identified as a source of a proton for O2 reduction. Three hydrogen-bonded waters are found to be present in the binuclear center after a sufficiently long simulation time. One is ligated to the CuB and could be associated with a water (or hydroxyl) identified in the crystal structure as the fourth ligand of CuB. The water hydrogen-bonded to the hydroxyl of Tyr288 is extremely persistent and well positioned to participate in O2 reduction. The third water is located where O2 is often suggested to reside in mechanistic studies of O2 reduction.     

Essential glycine in the proton channel of Escherichia coli transhydrogenase

The nicotinamide nucleotide transhydrogenases of mitochondria and bacteria are proton pumps that couple hydride ion transfer between NAD(H) and NADP(H) bound, respectively, to extramembranous domains I and III, to proton translocation by the membrane-intercalated domain II. Previous experiments have established the involvement of three conserved domain II residues in the proton pumping function of the enzyme: His91, Ser139, and Asn222, located on helices 9, 10, and 13, respectively. Eight highly conserved domain II glycines in helices 9, 10, 13, and 14 were mutated to alanine, and the mutant enzymes were assayed for hydride transfer between domains I and III and for proton translocation by domain II. One of the glycines on helix 14, Gly252, was further mutated to Cys, Ser, Thr, and Val, expression levels of the mutant enzymes were evaluated, and each was purified and assayed. The results show that Gly252 is essential for function and support a model for the proton channel composed of helices 9, 10, 13, and 14. Gly252 would allow spatial proximity of His91, Ser139, and Asn222 for proton conductance within the channel. Gly252 mutants are distinguished by high levels of cyclic transhydrogenation activity in the absence of added NADP(H) and by complete loss of proton pumping activity. The purified G252A mutant has <1% proton translocation and reverse transhydrogenation activity, retains 0.9 mol of NADP(H) per domain III, and has 96% intrinsic cyclic transhydrogenation activity, which does not exceed 100% upon the addition of NADP(H). These properties imply that Gly252 mutants exhibit a native-like domain II conformation while blocking proton translocation and coupled exchange of NADP(H) in domain III.     

Noncontact dipole effects on channel permeation. III. Anomalous proton conductance effects in gramicidin

Proton transport on water wires, of interest for many problems in membrane biology, is analyzed in side-chain analogs of gramicidin A channels. In symmetrical 0.1 N HCl solutions, fluorination of channel Trp(11), Trp-(13), or Trp(15) side chains is found to inhibit proton transport, and replacement of one or more Trps with Phe enhances proton transport, the opposite of the effects on K(+) transport in lecithin bilayers. The current-voltage relations are superlinear, indicating that some membrane field-dependent process is rate limiting. The interfacial dipole effects are usually assumed to affect the rate of cation translocation across the channel. For proton conductance, however, water reorientation after proton translocation is anticipated to be rate limiting. We propose that the findings reported here are most readily interpreted as the result of dipole-dipole interactions between channel waters and polar side chains or lipid headgroups. In particular, if reorientation of the water column begins with the water nearest the channel exit, this hypothesis explains the negative impact of fluorination and the positive impact of headgroup dipole on proton conductance.     

Structures of N-termini of helices in proteins.

We have surveyed 393 N-termini of alpha-helices and 156 N-termini of 3(10)-helices in 85 high resolution, non-homologous protein crystal structures for N-cap side-chain rotamer preferences, hydrogen bonding patterns, and solvent accessibilities. We find very strong rotamer preferences that are unique to N-cap sites. The following rules are generally observed for N-capping in alpha-helices: Thr and Ser N-cap side chains adopt the gauche - rotamer, hydrogen bond to the N3 NH and have psi restricted to 164 +/- 8 degrees. Asp and Asn N-cap side chains either adopt the gauche - rotamer and hydrogen bond to the N3 NH with psi = 172 +/- 10 degrees, or adopt the trans rotamer and hydrogen bond to both the N2 and N3 NH groups with psi = 1-7 +/- 19 degrees. With all other N-caps, the side chain is found in the gauche + rotamer so that the side chain does not interact unfavorably with the N-terminus by blocking solvation and psi is unrestricted. An i, i + 3 hydrogen bond from N3 NH to the N-cap backbone C = O in more likely to form at the N-terminus when an unfavorable N-cap is present. In the 3(10)-helix Asn and Asp remain favorable N-caps as they can hydrogen bond to the N2 NH while in the trans rotamer; in contrast, Ser and Thr are disfavored as their preferred hydrogen bonding partner (N3 NH) is inaccessible. This suggests that Ser is the optimum choice of N-cap when alpha-helix formation is to be encouraged while 3(10)-helix formation discouraged. The strong energetic and structural preferences found for N-caps, which differ greatly from positions within helix interiors, suggest that N-caps should be treated explicitly in any consideration of helical structure in peptides or proteins.     

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N(2)O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666-70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb(3) oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.     

Investigating the dielectric effects of channel pore water on the electrostatic barriers of the permeation ion by the finite difference Poisson-Boltzmann method

In this paper, the finite difference Poisson-Boltzmann (FDPB) method with four dielectric constants is developed to study the effect of dielectric saturation on the electrostatic barriers of the permeation ion. In this method, the inner shape of the channel pore is explicitly represented, and the fact that the dielectric constant inside the channel pore is different from that of bulk water is taken into account. A model channel system which is a right-handed twist bundle with four α-helical segments is provided for this study. From the FDPB calculations, it is found that the difference of the ionic electrostatic solvation energy for wider domains depends strongly on the pore radius in the vicinity of the ion when the pore dielectric constant is changed from 78 to 5. However, the electrostatic solvation energy of the permeation ion can not be significantly affected by the dielectric constant in regions with small pore radii. Our results indicate that the local electrostatic interactions inside the ion channel are of major importance for ion electrostatic solvation energies, and the effect of dielectric saturation on the electrostatic barriers is coupled to the interior channel dimensions. Received: 28 January 1997 / Accepted: 24 September 1997     

Plasticity of proton pathway structure and water coordination in cytochrome c oxidase

Cytochrome c oxidase (CytcO) is a redox-driven, membrane-bound proton pump. One of the proton transfer pathways of the enzyme, the D pathway, used for the transfer of both substrate and pumped protons, accommodates a network of hydrogen-bonded water molecules that span the distance between an aspartate (Asp(132)), near the protein surface, and glutamate Glu(286), which is an internal proton donor to the catalytic site. To investigate how changes in the environment around Glu(286) affect the mechanism of proton transfer through the pathway, we introduced a non-hydrogen-bonding (Ala) or an acidic residue (Asp) at position Ser(197) (S197A or S197D), located approximately 7 A from Glu(286). Although Ser(197) is hydrogen-bonded to a water molecule that is part of the D pathway "proton wire," replacement of the Ser by an Ala did not affect the proton transfer rate. In contrast, the S197D mutant CytcO displayed a turnover activity of approximately 35% of that of the wild-type CytcO, and the O(2) reduction reaction was not linked to proton pumping. Instead, a fraction of the substrate protons was taken from the positive ("incorrect") side of the membrane. Furthermore, the pH dependence of the proton transfer rate was altered in the mutant CytcO. The results indicate that there is plasticity in the water coordination of the proton pathway, but alteration of the electrostatic potential within the pathway results in uncoupling of the proton translocation machinery.     

Mechanism of enolase: the crystal structure of enolase-Mg2(+)-2-phosphoglycerate/phosphoenolpyruvate complex at 2.2-A resolution

Enolase in the presence of Mg2+ catalyzes the elimination of H2O from 2-phosphoglyceric acid (PGA) to form phosphoenolpyruvate (PEP) and the reverse reaction, the hydration of PEP to PGA. The structure of the ternary complex yeast enolase-Mg2(+)-PGA/PEP has been determined by X-ray diffraction and refined by crystallographic restrained least-squares to an R = 16.9% for those data with I/sigma (I) greater than or equal to 2 to 2.2-A resolution with a good geometry of the model. The structure indicates the substrate molecule in the active site has its hydroxyl group coordinated to the Mg2+ ion. The carboxylic group interacts with the side chains of His373 and Lys396. The phosphate group is H-bonded to the guanidinium group of Arg374. A water molecule H-bonded to the carboxylic groups of Glu168 and Glu211 is located at a 2.6-A distance from carbon-2 of the substrate in the direction of its proton. We propose that this cluster functions as the base abstracting the proton in the catalytic process. The proton is probably transferred, first to the water molecule, then to Glu168, and further to the substrate hydroxyl to form a water molecule. Some analogy is apparent between the initial stages of the enolase reverse reaction, the hydration of PEP, and the proteolytic mechanism of the metallohydrolases carboxypeptidase A and thermolysin. The substrate/product binding is accompanied by large movements of loops Ser36-His43 and Ser158-Gly162. The role of these conformational changes is not clear at this time.     

Molecular dynamics simulation of cation motion in water-filled gramicidinlike pores.

A model calculation is carried out to study the potential energy profile of a sodium ion with several water molecules inside a simplified model of the gramicidin ion channel. The sodium ion is treated as a Lennard-Jones sphere with a point charge at its center. The Barnes polarizable water model is used to mimic the water molecules. A polarizable and deformable gramicidinlike channel is constructed based on the model obtained by Koeppe and Kimura. Potential minima and saddle points are located and the static energy barriers are computed. The potential minima at the two mouths of the channel exhibit an aqueous solvation structure very different from that at any of the interior minima. These sites are approximately 23.6 and 24.4 A apart for binding of a sodium ion and a cesium ion, respectively. Ionic motion from these exterior sites to the first interior minimum requires substantial rearrangement of the waters of solvation; this rearrangement may be the hydration/dehydration step in ionic permeation through the channel. Based on these results, a mechanism by which the sodium ion moves from the exterior binding site to the interior of the channel is proposed. Our model channel accommodates about eight water molecules and the transport of the ion and water within the channel is found to be single file. Results of less extensive calculations for Cs+ and Li+ ions in a channel with or without water are also reported.     

Structural elements involved in proton translocation by cytochrome c oxidase as revealed by backbone amide hydrogen-deuterium exchange of the E286H mutant

Cytochrome c oxidase is the terminal electron acceptor in the respiratory chains of aerobic organisms and energetically couples the reduction of oxygen to water to proton pumping across the membrane. The mechanisms of proton uptake, gating, and pumping have yet to be completely elucidated at the molecular level for these enzymes. For Rhodobacter sphaeroides CytcO (cytochrome aa3), it appears as though the E286 side chain of subunit I is a branching point from which protons are shuttled either to the catalytic site for O2 reduction or to the acceptor site for pumped protons. Amide hydrogen-deuterium exchange mass spectrometry was used to investigate how mutation of this key branching residue to histidine (E286H) affects the structures and dynamics of four redox intermediate states. A functional characterization of this mutant reveals that E286H CytcO retains approximately 1% steady-state activity that is uncoupled from proton pumping and that proton transfer from H286 is significantly slowed. Backbone amide H-D exchange kinetics indicates that specific regions of CytcO, perturbed by the E286H mutation, are likely to be involved in proton gating and in the exit pathway for pumped protons. The results indicate that redox-dependent conformational changes around E286 are essential for internal proton transfer. E286H CytcO, however, is incapable of these specific conformational changes and therefore is insensitive to the redox state of the enzyme. These data support a model where the side chain conformation of E286 controls proton translocation in CytcO through its interactions with the proton gate, which directs the flow of protons either to the active site or to the exit pathway. In the E286H mutant, the proton gate does not function properly and the exit channel is unresponsive. These results provide new insight into the structure and mechanism of proton translocation by CytcO.     

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.     

Redox-driven proton pumping by heme-copper oxidases

One of the key problems of molecular bioenergetics is the understanding of the function of redox-driven proton pumps on a molecular level. One such class of proton pumps are the heme-copper oxidases. These enzymes are integral membrane proteins in which proton translocation across the membrane is driven by electron transfer from a low-potential donor, such as, e.g. cytochrome c, to a high-potential acceptor, O(2). Proton pumping is associated with distinct exergonic reaction steps that involve gradual reduction of oxygen to water. During the process of O(2) reduction, unprotonated high pK(a) proton acceptors are created at the catalytic site. Initially, these proton acceptors become protonated as a result of intramolecular proton transfer from a residue(s) located in the membrane-spanning part of the enzyme, but removed from the catalytic site. This residue is then reprotonated from the bulk solution. In cytochrome c oxidase from Rhodobacter sphaeroides, the proton is initially transferred from a glutamate, E(I-286), which has an apparent pK(a) of 9.4. According to a recently published structure of the enzyme, the deprotonation of E(I-286) is likely to result in minor structural changes that propagate to protonatable groups on the proton output (positive) side of the protein. We propose that in this way, the free energy available from the O(2) reduction is conserved during the proton transfer. On the basis of the observation of these structural changes, a possible proton-pumping model is presented in this paper. Initially, the structural changes associated with deprotonation of E(I-286) result in the transfer of a proton to an acceptor for pumped protons from the input (negative) side of the membrane. After reprotonation of E(I-286) this acceptor releases a proton to the output side of the membrane.     

The H channel is not a proton transfer path in yeast cytochrome c oxidase

Cytochrome c oxidases (CcOs) in the respiratory chains of mitochondria and bacteria are primary consumers of molecular oxygen, converting it to water with the concomitant pumping of protons across the membrane to establish a proton electrochemical gradient. Despite a relatively well understood proton pumping mechanism of bacterial CcOs, the role of the H channel in mitochondrial forms of CcO remains debated. Here, we used site-directed mutagenesis to modify a central residue of the lower span of the H channel, Q413, in the genetically tractable yeast Saccharomyces cerevisiae. Exchange of Q413 to several different amino acids showed no effect on rates and efficiencies of respiratory cell growth, and redox potential measurements indicated minimal electrostatic interaction between the 413 locus and the nearest redox active component heme a. These findings clearly exclude a primary role of this section of the H channel in proton pumping in yeast CcO. In agreement with the experimental data, atomistic molecular dynamics simulations and continuum electrostatic calculations on wildtype and mutant yeast CcOs highlight potential bottlenecks in proton transfer through this route. Our data highlight the preference for neutral residues in the 413 locus, precluding sufficient hydration for formation of a proton conducting wire.     

Mechanism of proton transport by plant plasma membrane proton ATPases

The mechanism of proton translocation by P-type proton ATPases is poorly defined. Asp684 in transmembrane segment M6 of the Arabidopsis thaliana AHA2 plasma membrane P-type proton pump is suggested to act as an essential proton acceptor during proton translocation. Arg655 in transmembrane segment M5 seems to be involved in this proton translocation too, but in contrast to Asp684, is not essential for transport. Asp684 may participate in defining the E₁ proton-binding site, which could possibly exist as a hydronium ion coordination center. A model of proton translocation of AHA2 involving the side chains of amino acids Asp684 and Arg655 is discussed.     

Electrostatic modeling of dipole-ion interactions in gramicidinlike channels

Using an electrostatic model for the pore and membrane region in a gramicidinlike channel, the effect of dipoles located inside the membrane on the ion transport are analyzed. Calculated energy profiles for different orientations of dipoles show a predominant influence of their radial components. The results qualitatively agree with experimental measurements of conductance on different modified gramicidins and allow to understand the important role of polar side chains on ion permeation.     

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.     

Ab initio molecular orbital study on the effects of zinc ion,its ligands and ionic amino acid residues for proton transfer energetics between glu 270 and Zn co-ordinated water molecule in carboxypeptidase A

Thezinc-water-Glu 270 system was reported from the X-ray crystallographic study of native carboxypeptidase A(CPA) (Lipscomb et al., 1968). General base catalysis by the γ-carboxylate of Glu 270 was proposed for peptidase activity of CPA. The effects of zinc ion and its ligands (Glu 72, His 69-Asp 142, His 196) for proton transfer between Glu 270 and Zn co-ordinated water molecule in CPA were studied by the ab initio SCFLCAO-MO method. The results show that the proton transfer from the Zn co-ordinated water molecule to the γ-carboxylate of Glu 270 is greatly promoted by the Zn ion and, conversely, is greatly inhibited by its ligands. The facilitation effect of Zn ion and the inhibition effect of its ligands for the proton transfer were analysed by using the energy decomposition analysis. Moreover, calculations including all side chains of ionic amino acid residues and main chain residues in CPA as point fractional charges were performed. The results show that the proton transfer is affected by the ionic amino acid residues and is not affected by the main chain residues.