Helical water wires

A helitetrahedral model has been proposed to help explain reports of low-frequency oscillations in pure water following electromagnetic excitation at the hydronium ion cyclotron resonance frequency. The Lorentz force and the intrinsic structure constrain the motion of the H₃O⁺ ion so that it enjoys a unique form of proton-hopping, one whose path is helical. This model may also explain the numerous previously observed cyclotron resonance (ICR) biological couplings for cations other than hydronium by merely substituting hydrogen-bonded versions of these for hydronium in the tetrahedral structure. Thus the effectiveness of resonance stimulation in biological systems is explained in terms of the enhanced conductivity and reduced scattering associated with proton-hopping. It is further shown that the addition of charge-balancing hydroxyl ions act to enable oscillatory electric dipole moments that propagate along the helical axis, giving rise to weak power (≈ femtoWatts) radiation patterns. It is conceivable that the radiation associated with this process may play a role in the interactions at the interface between water and living matter.     

Weak-field H3O+ ion cyclotron resonance alters water refractive index

Heretofore only observed in living systems, we report that weak-field ion cyclotron resonance (ICR) also occurs in inanimate matter. Weak magnetic field (50 nT) hydronium ICR at the field combination (7.84 Hz, 7.5 µT) markedly changes water structure, as evidenced by finding an altered index of refraction exactly at this combined field. This observation utilizes a novel technique which measures the scattering of a He–Ne laser beam as the sample is exposed to a ramped magnetic field frequency. In addition to the hydronium resonance, we find evidence of ICR coupling to a more massive structure, possibly a tetrahedral combination of three waters and a single hydronium ion. To check our observations, we extended this technique to D₂O, successfully predicting the specific ICR charge-to-mass ratio for D₃O⁺ that alters the index of refraction.     

Lorentz force in water: evidence that hydronium cyclotron resonance enhances polymorphism

Abstract

There is an ongoing question regarding the structure forming capabilities of water at ambient temperatures. To probe for different structures, we studied effects in pure water following magnetic field exposures corresponding to the ion cyclotron resonance of H₃O⁺. Included were measurements of conductivity and pH. We find that under ion cyclotron resonance (ICR) stimulation, water undergoes a transition to a form that is hydroxonium-like, with the subsequent emission of a transient 48.5?Hz magnetic signal, in the absence of any other measurable field. Our results indicate that hydronium resonance stimulation alters the structure of water, enhancing the concentration of EZ-water. These results are not only consistent with Del Giudice's model of electromagnetically coherent domains, but they can also be interpreted to show that these domains exist in quantized spin states.     

Selectivity and conductance among the glycerol and water conducting aquaporin family of channels

The atomic structures of a transmembrane water plus glycerol conducting channel (GlpF), and now of aquaporin Z (AqpZ) from the same species, Escherichia coli, bring the total to three atomic resolution structures in the aquaporin (AQP) family. Members of the AQP family each assemble as tetramers of four channels. Common helical axes support a wider channel in the glycerol plus water channel paradigm, GlpF. Water molecules form a single hydrogen bonded file throughout the 28 A long channel in AqpZ. The basis for absolute exclusion of proton or hydronium ion conductance through the line of water is explored using simulations.     

Structure and dynamics of hydronium in the ion channel gramicidin A.

The effects of the hydronium ion, H(3)0+, on the structure of the ion channel gramicidin A and the hydrogen-bonded network of waters within the channel were studied to help elucidate a possible mechanism for proton transport through the channel. Several classical molecular dynamics studies were carried out with the hydronium in either the center of a gramicidin monomer or in the dimer junction. Structural reorganization of the channel backbone was observed for different hydronium positions, which were most apparent when the hydronium was within the monomer. In both cases the average O-O distance between the hydronium ion and its nearest neighbor water molecule was found to be approximately 2.55 A, indicating a rather strong hydrogen bond. Importantly, a subsequent break in the hydrogen-bonded network between the nearest neighbor and the next-nearest neighbor(approximately 2.7 -3.0 A) was repeatedly observed. Moreover, the carbonyl groups of gramicidin A were found to interact with the charge on the hydronium ion, helping in its stabilization. These facts may have significant implications for the proton hopping mechanism. The presence of the hydronium ion in the channel also inhibits to some degree the reorientational motions of the channel water molecules.     

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.     

A Cellular Automata Model of Proton Hopping Down a Channel

Proton hopping is the process where a H‐atom on a hydronium ion forms a H‐bond with the O‐atom of a neighboring H₂O molecule. There is then an exchange of bonding forces when that covalent bond of the H‐atom in the hydronium ion changes to a H‐bond, and the previous H‐bond changes to a covalent bond with the neighboring O‐atom. The neighboring molecule now becomes a hydronium (H₃O⁺) ion. This process repeats itself very rapidly among neighboring hydronium and H₂O molecules. There is a flow of protonic character through bulk H₂O, referred to as proton hopping. This process carries information through living systems where H₂O is present. A cellular automata model of proton hopping down a channel has been created and studied. Variations in the rate of proton entry into the channel and the effects of the polar character of the channel walls was studied using the model. The behavior of the models corresponds to experimental results.     

Proton exchange of nucleic acids: the amino protons of 5'-AMP and poly A

Approximate values of the hydronium and hydroxyl second order rate constants for the exchange of the amino protons of 5′-AMP and Poly A were obtained from the observed variation with pH of the 100 MHz PMR resonance signal line width for the amino protons. These values are several orders of magnitude lower than corresponding rate constants for amino protons and purine hydrogens, which may resolve some important interpretative difficulties that are encountered both in the hydrogen exchange of helical polynucleotides and in recent data on the exchange of the non-hydrogen bonded amino protons of DNA.     

Charge delocalization in proton channels, I: the aquaporin channels and proton blockage

The explicit contribution to the free energy barrier and proton conductance from the delocalized nature of the excess proton is examined in aquaporin channels using an accurate all-atom molecular dynamics computer simulation model. In particular, the channel permeation free energy profiles are calculated and compared for both a delocalized (fully Grotthuss shuttling) proton and a classical (nonshuttling) hydronium ion along two aquaporin channels, Aqp1 and GlpF. To elucidate the effects of the bipolar field thought to arise from two alpha-helical macrodipoles on proton blockage, free energy profiles were also calculated for computational mutants of the two channels where the bipolar field was eliminated by artificially discharging the backbone atoms. Comparison of the free energy profiles between the proton and hydronium cases indicates that the magnitude of the free energy barrier and position of the barrier peak for the fully delocalized and shuttling proton are somewhat different from the case of the (localized) classical hydronium. The proton conductance through the two aquaporin channels is also estimated using Poisson-Nernst-Planck theory for both the Grotthuss shuttling excess proton and the classical hydronium cation.     

Many-body energies during proton transfer in an aqueous system

The energetics of the mechanism of proton transfer from a hydronium ion to one of the water molecules in its first solvation shell are studied using density functional theory and the Møller–Plesset perturbation (MP2) method. The potential energy surface of the proton transfer mechanism is obtained at the B3LYP and MP2 levels with the 6-311++G** basis set. Many-body analysis is applied to the proton transfer mechanism to obtain the change in relaxation energy, two-body, three-body and four-body energies when proton transfer occurs from the hydronium ion to one of the water molecules in its first solvation shell. It is observed that the binding energy (BE) of the complex decreases during the proton transfer process at both levels of theory. During the proton transfer process, the % contribution of the total two-body energy to the binding energy of the complex increases from 62.9 to 68.09% (39.9 to 45.95%), and that of the total three-body increases from 25.9 to 27.09% (24.16 to 26.17%) at the B3LYP/6-311++G** (MP2/ 6-311++G**) level. There is almost no change in the water–water–water three-body interaction energy during the proton transfer process at both levels of theory. The contribution of the relaxation energy and the total four-body energy to the binding energy of the complex is greater at the MP2 level than at the B3LYP level. Significant differences are found between the relaxation energies, the hydronium–water interaction energies and the four-body interaction energies at the B3LYP and MP2 levels.     

The Creation of Proton Hopping from a Drug?Receptor Encounter

We extend recent modeling studies of proton hopping, used to describe the functioning of membrane channels and axon nerve conduction, to offer an explanation of the initiation of the nerve impulse at an effector? ligand encounter. This encounter is proposed to create a hydronium ion in the vicinity of the effector and ligand, which leads to a continuous flow of protons, called proton hopping, through water adjacent to this encounter. This proton hopping is proposed to be the message carried from the encounter to the axon of a particular nerve system associated with that particular effector? ligand system.     

Splitting of the spectra of MHD waves and the structure of a satellite alfvén resonance in a cold plasma in a strong axial magnetic field and a weak helical field

The possibility is demonstrated of splitting the eigenfrequencies of MHD plasma waves in a stellarator with a weakly rippled helical confining magnetic field. The distribution of the fields of an Alfvén wave in the satellite Alfvén resonance region is investigated when the influence of the helical ripple in a confining magnetic field on the resonance structure is comparable with the effects of the finite ion Larmor radius, electron inertia, and collisions between plasma particles.     

The Hydronium Tetrafluoroborate Dimer in Nonpolar Media and its Proton NMR Spectrum

Hydronium tetrafluoroborate ion pairs, H₃O⁺·BF₄ ^- have been shown computationally to be unstable toward decomposition, in the absence of solvation or electrostatic interactions existing in crystals. As the proton NMR spectrum of a hydronium salt with the octanesulfonate-antimony pentachloride complex anion was reported in freon solution, we investigated the hypothesis that larger ionic clusters were present in the nonpolar solvent. It was found that the dimer (H₃O⁺·BF₄ ^-)₂ was stable at the MP2/6-31G* level. GIAO-B3LYP chemical shift calculations with the same basis set and also with the 6-31G**, 6-31++G**, 6-311++G**, dzvp, tzp, tz2p, and qz2p basis sets conducted on the hydronium fluoroborate dimer reproduce the main features of the experimental spectrum: the existence of two signals with a two-to-one intensity ratio and the more intense resonance at higher frequency (more deshielded). The alternative structures, of hydronium tetrafluoroborate ion pairs with one and with two hydrogen bonds between anion and cation, give calculated chemical shifts which are farther from the experimental values.     

Molecular Disruption of the Power Stroke in the ATP-binding Cassette Transport Protein MsbA

ATP-binding cassette transporters affect drug pharmacokinetics and are associated with inherited human diseases and impaired chemotherapeutic treatment of cancers and microbial infections. Current alternating access models for ATP-binding cassette exporter activity suggest that ATP binding at the two cytosolic nucleotide-binding domains provides a power stroke for the conformational switch of the two membrane domains from the inward-facing conformation to the outward-facing conformation. In outward-facing crystal structures of the bacterial homodimeric ATP-binding cassette transporters MsbA from Gram-negative bacteria and Sav1866 from Staphylococcus aureus, two transmembrane helices (3 and 4) in the membrane domains have their cytoplasmic extensions in close proximity, forming a tetrahelix bundle interface. In biochemical experiments on MsbA from Escherichia coli, we show for the first time that a robust network of inter-monomer interactions in the tetrahelix bundle is crucial for the transmission of nucleotide-dependent conformational changes to the extracellular side of the membrane domains. Our observations are the first to suggest that modulation of tetrahelix bundle interactions in ATP-binding cassette exporters might offer a potent strategy to alter their transport activity.     

A phylogenetic study of U4 snRNA reveals the existence of an evolutionarily conserved secondary structure corresponding to 'free' U4 snRNA

The nucleotide sequence of Physarum polycephalum U4 snRNA*** was determined and compared to published U4 snRNA sequences. The primary structure of P polycephalum U4 snRNA is closer to that of plants and animals than to that of fungi. But, both fungi and P polycephalum U4 snRNAs are missing the 3' terminal hairpin and this may be a common feature of lower eucaryote U4 snRNAs. We found that the secondary structure model we previously proposed for 'free' U4 snRNA is compatible with the various U4 snRNA sequences published. The possibility to form this tetrahelix structure is preserved by several compensatory base substitutions and by compensatory nucleotide insertions and deletions. According to this finding, association between U4 and U6 snRNAs implies the disruption of 2 internal helical structures of U4 snRNA. One has a very low free energy, but the other, which represents one-half of the helical region of the 5' hairpin, requires 4 to 5 kcal to be open. The remaining part of the 5' hairpin is maintained in the U4/U6 complex and we observed the conservation, in all U4 snRNAs studied, of a U bulge residue at the limit between the helical region which has to be melted and that which is maintained. The 3' domain of U4 snRNA is less conserved in both size and primary structure than the 5' domain; its structure is also more compact in the RNA in solution. In this domain, only the Sm binding site and the presence of a bulge nucleotide in the hairpin on the 5' side of the Sm site are conserved throughout evolution.     

NMR studies of the aggregation of glucagon-like peptide-1: formation of a symmetric helical dimer

Nuclear magnetic resonance (NMR) spectroscopy reveals that higher-order aggregates of glucagon-like peptide-1-(7-36)-amide (GLP-1) in pure water at pH 2.5 are disrupted by 35% 2,2,2-trifluoroethanol (TFE), and form a stable and highly symmetric helical self-aggregate. NMR spectra show that the helical structure is identical to that formed by monomeric GLP-1 under the same experimental conditions [Chang et al., Magn. Reson. Chem. 37 (2001) 477-483; Protein Data Bank at RCSB code: 1D0R], while amide proton exchange rates reveal a dramatic increase of the stability of the helices of the self-aggregate. Pulsed-field gradient NMR diffusion experiments show that the TFE-induced helical self-aggregate is a dimer. The experimental data and model calculations indicate that the dimer is a parallel coiled coil, with a few hydrophobic residues on the surface that may cause aggregation in pure water. The results suggest that the coiled coil dimer is an intermediate state towards the formation of higher aggregates, e.g. fibrils.     

On the question of hydronium binding to ATP-synthase membrane rotors

A recently determined atomic structure of an H⁺-coupled ATP-synthase membrane rotor has revived the long-standing question of whether protons may be bound to these structures in the form of a hydronium ion. Using both classical and quantum-mechanical simulations, we show that this notion is implausible. Ab initio molecular dynamics simulations of the binding site demonstrate that the putative H₃O⁺ deprotonates within femtoseconds. The bound proton is thus transferred irreversibly to the carboxylate side chain found in the ion-binding sites of all ATP-synthase rotors. This result is consistent with classical simulations of the rotor in a phospholipid membrane, on the 100-nanosecond timescale. These simulations show that the hydrogen-bond network seen in the crystal structure is incompatible with a bound hydronium. The observed coordination geometry is shown to correspond instead to a protonated carboxylate and a bound water molecule. In conclusion, this study underscores the notion that binding and transient storage of protons in the membrane rotors of ATP synthases occur through a common chemical mechanism, namely carboxylate protonation.     

Alternative proton binding mode in ATP synthases

ATP synthases are rotary engines which use the energy stored in a transmembrane electrochemical gradient of protons or sodium ions to catalyze the formation of ATP by ADP and inorganic phosphate. Current models predict that protonation/deprotonation of specific amino acids of the rotating c-ring, extracting protons from one side and delivering them to the other side of the membrane, are at the core of the proton translocation mechanism of these enzymes. In this minireview, an alternative proton binding mechanism is presented, considering hydronium ion coordination as proposed earlier. Biochemical data and structural considerations provide evidence for two different proton binding modes in the c-ring of H⁺-translocating ATP synthases. Recent investigations in several other proton translocating membrane proteins suggest, that hydronium ion coordination by proteins might display a general principle which was so far underestimated in ATP synthases.     

Global structure of four-way RNA junctions studied using fluorescence resonance energy transfer.

Four-way helical junctions are found widely in natural RNA species. In this study, we have studied the conformation of two junctions by fluorescence resonance energy transfer. We show that the junctions are folded by pairwise coaxial helical stacking, forming one predominant stacking conformer in both examples studied. At low magnesium ion concentrations, the helical axes of both junctions are approximately perpendicular. One junction undergoes a rotation into a distorted antiparallel structure induced by the binding of a single magnesium ion. By contrast, the axes of the four-way junction of the U1 snRNA remain approximately perpendicular under all conditions examined, and we have determined the stacking conformer adopted.     

Involvement of the H3O+-Lys-164 -Gln-161-Glu-345 charge transfer pathway in proton transport of gastric H+,K+-ATPase

Gastric H(+),K(+)-ATPase is shown to transport 2 mol of H(+)/mol of ATP hydrolysis in isolated hog gastric vesicles. We studied whether the H(+) transport mechanism is due to charge transfer and/or transfer of hydronium ion (H(3)O(+)). From transport of [(18)O]H(2)O, 1.8 mol of water molecule/mol of ATP hydrolysis was found to be transported. We performed a molecular dynamics simulation of the three-dimensional structure model of the H(+),K(+)-ATPase alpha-subunit at E(1) conformation. It predicts the presence of a charge transfer pathway from hydronium ion in cytosolic medium to Glu-345 in cation binding site 2 (H(3)O(+)-Lys-164 -Gln-161-Glu-345). No charge transport pathway was formed in mutant Q161L, E345L, and E345D. Alternative pathways (H(3)O(+)-Gln-161-Glu-345) in mutant K164L and (H(3)O(+)-Arg-105-Gln-161-Gln-345) in mutant E345Q were formed. The H(+),K(+)-ATPase activity in these mutants reflected the presence and absence of charge transfer pathways. We also found charge transfer from sites 2 to 1 via a water wire and a charge transfer pathway (H(3)O(+)-Asn-794 -Glu-797). These results suggest that protons are charge-transferred from the cytosolic side to H(2)O in sites 2 and 1, the H(2)O comes from cytosolic medium, and H(3)O(+) in the sites are transported into lumen during the conformational transition from E(1)PtoE(2)P.