K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study

Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K(+) compared with Na(+) ions. K(+) ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na(+) ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K(+) ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na(+) ion. Whereas K(+) ions prefer to sit within a cage of eight oxygen atoms of the filter, Na(+) ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K(+) ions (in addition to the differences in dehydration energy between K(+) and Na(+)) and the block of KcsA by internal Na(+) ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na(+) than by K(+) ions.     

Characterizing the binding of annexin V to a lipid bilayer using molecular dynamics simulations

Annexins play critical roles in membrane organization, membrane trafficking and vesicle transport. The family members share the ability to bind to membranes with high affinities, but the interactions between annexins and membranes remain unclear. Here, using long‐time molecular dynamics simulations, we provide detailed information for the binding of an annexin V trimer to a POPC/POPS lipid bilayer. Calcium ions function as bridges between several negatively charged residues of annexin V and the oxygen atoms of lipids. The preferred calcium‐bridges are those formed via the carboxyl oxygen atoms of POPS lipids. H‐bonds and hydrophobic interactions formed by several critical residues have also been observed in the annexin‐membrane interface. The annexin‐membrane binding causes small changes of annexin trimer structures, while has significant effects on lipid bilayer structures. The lipid bilayer shows a bent shape and forms a concave region in the annexin‐membrane interaction interface, which provides an atomic‐level evidence to support the view that annexins could disturb the stability of lipids and bend membranes. This study provides insights into the commonly occurring PS‐dependent and calcium‐dependent binding of proteins to membranes. Proteins 2014; 82:312–322. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulation of a dipalmitoylphosphatidylcholine bilayer with NaCl

Molecular dynamics simulations are performed on two hydrated dipalmitoylphosphatidylcholine bilayer systems: one with pure water and one with added NaCl. Due to the rugged nature of the membrane/electrolyte interface, ion binding to the membrane surface is characterized by the loss of ion hydration. Using this structural characterization, binding of Na(+) and Cl(-) ions to the membrane is observed, although the binding of Cl(-) is seen to be slightly weaker than that of Na(+). Dehydration is seen to occur to a different extent for each type of ion. In addition, the excess binding of Na(+) gives rise to a net positive surface charge density just outside the bilayer. The positive density produces a positive electrostatic potential in this region, whereas the system without salt shows an electrostatic potential of zero.     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

Molecular dynamics simulations of pentapeptides at interfaces: salt bridge and cation-pi interactions

Peptide-membrane interactions are important for understanding the binding, partitioning, and folding of membrane proteins; the activity of antimicrobial and fusion peptides; and a number of other processes. We describe molecular dynamics simulations (10-25 ns) of two pentapeptides Ace-WLXLL (with X = Arg or Lys side chain) (White, S. H., and Wimley, W.C. (1996) Nat. Struct. Biol. 3, 842-848) in water and three different membrane mimetic systems: (i) a water/cyclohexane interface, (ii) water-saturated octanol, and (iii) a solvated dioleoylphosphatidylcholine bilayer. A salt bridge is found between the protonated Arg or Lys side chains with the carboxyl terminus at the three interfaces. In water/cyclohexane, the salt bridge is most exposed to the water phase and least stable. In water/octanol and the lipid bilayer systems, the salt bridge once formed persists throughout the simulations. In the lipid bilayer, the salt bridge is more stable when the peptide penetrates deeper into the bilayer. In one of two peptides, a cation-pi interaction between the Arg and the Trp side chains is stable in the lipid bilayer for about 15 ns before breaking. In all cases, the conformations of the peptides are restricted by their presence at the interface and can be assigned to a few major conformational clusters. Side chains facing the water phase are most mobile. In the lipid bilayer, the peptides remain in the interface area, where they overlap with the carbonyl area of the lipid bilayer and perturb the local density profile of the bilayer. The tryptophan side chain remains in the water-lipid interface, where it interacts with the lipid choline group and forms hydrogen bonds with the ester carbonyl of the lipid and with water in the interface.     

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.     

Examining the contributions of lipid shape and headgroup charge on bilayer behavior

To better understand bilayer property dependency on lipid electrostatics and headgroup size, we use atomistic molecular dynamics simulations to study negatively charged and neutral lipid membranes. We compare the negatively charged phosphatidic acid (PA), which at physiological pH and salt concentration has a negative spontaneous curvature, with the negatively charged phosphatidylglycerol (PG) and neutrally charged phosphatidylcholine (PC), both of which have zero spontaneous curvature. The PA lipids are simulated using two different sets of partial charges for the headgroup and the varied charge distribution between the two PA systems results in significantly different locations for the Na⁺ ions relative to the water/membrane interface. For one PA system, the Na⁺ ions are localized around the phosphate group. In the second PA system, the Na⁺ ions are located near the ester carbonyl atoms, which coincides with the preferred location site for the PG Na⁺ ions. We find that the Na⁺ ion location has a larger effect on bilayer fluidity properties than lipid headgroup size, where the A_lipid and acyl chain order parameter values are more similar between the PA and PG bilayers that have Na⁺ ions located near the ester groups than between the two PA bilayers.     

Cholesterol effects on the phosphatidylcholine bilayer polar region: a molecular simulation study

A molecular dynamics (MD) simulation of a fully hydrated, liquid-crystalline dimyristoylphosphatidylcholine (DMPC)-Chol bilayer membrane containing approximately 22 mol% Chol was carried out for 4.3 ns. The bilayer reached thermal equilibrium after 2.3 ns of MD simulation. A 2.0-ns trajectory generated during 2.3-4.3 ns of MD simulation was used for analyses to determine the effects of Chol on the membrane/water interfacial region. In this region, 70% of Chol molecules are linked to DMPC molecules via short-distance interactions, where the Chol hydroxyl group (OH-Chol) is 1) charge paired to methyl groups of the DMPC choline moiety ( approximately 34%), via the hydroxyl oxygen atom (Och); 2) water bridged to carbonyl ( approximately 19%) and nonester phosphate ( approximately 14%) oxygen atoms, via both Och and the hydroxyl hydrogen atom (Hch); and 3) directly hydrogen (H) bonded to carbonyl ( approximately 11%) and nonester phosphate ( approximately 5%) oxygen atoms, via Hch ( approximately 17% of DMPC-Chol links are multiple). DMPC's gamma-chain carbonyl oxygen atom is involved in 44% of water bridges and 51% of direct H bonds formed between DMPC and Chol. On average, a Chol molecule forms 0.9 links with DMPC molecules, while a DMPC molecule forms 2.2 and 0.3 links with DMPC and Chol molecules, respectively. OH-Chol makes hydrogen bonds with 1.1 water molecules, preferentially via Hch. The average number of water molecules H bonded to the DMPC headgroup is increased by 7% in the presence of Chol. These results indicate that inclusion of Chol decreases interlipid links and increases hydration in the polar region of the membrane.     

Structure, dynamics, and hydration of POPC/POPS bilayers suspended in NaCl, KCl, and CsCl solutions

Effects of alkali metal chlorides on the properties of mixed negatively charged lipid bilayers are experimentally measured and numerically simulated. Addition of 20mol% of negatively charged phosphatidylserine to zwitterionic phosphatidylcholine strengthens adsorption of monovalent cations revealing their specificity, in the following order: Cs(+)相似文献   

Molecular dynamics investigation of the structure of a fully hydrated gel-phase dipalmitoylphosphatidylcholine bilayer.

We report the results of a constant pressure and temperature molecular dynamics simulation of a gel-phase dipalmitoylphosphatidylcholine bilayer with nw = 11.8 water molecules/lipid at 19 degrees C. The results of the simulation were compared in detail with a variety of x-ray and neutron diffraction data. The average positions of specific carbon atoms along the bilayer normal and the interlamellar spacing and electron density profile were in very good agreement with neutron and x-ray diffraction results. The area per lipid and the details of the in-plane hydrocarbon chain structure were in excellent agreement with wide-angle x-ray diffraction results. The only significant deviation is that the chains met in a pleated arrangement at the bilayer center, although they should be parallel. Novel discoveries made in the present work include the observation of a bimodal headgroup orientational distribution. Furthermore, we found that there are a significant number of gauche conformations near the ends of the hydrocarbon chains and, in addition to verifying a previous suggestion that there is partial rotational ordering in the hydrocarbon chains, that the two chains in a given molecule are inequivalent with respect to rotations. Finally, we have investigated the lipid/water interface and found that the water penetrates beneath the headgroups, but not as far as the carbonyl groups, that the phosphates are strongly hydrated almost exclusively at the nonesterified oxygen atoms, and that the hydration of the ammonium groups is more diffuse, with some water molecules concentrated in the grooves between the methyl groups.     

Molecular dynamics simulations of the gramicidin A-dimyristoylphosphatidylcholine system with an ion in the channel pore region

To investigate the process of ion permeation in an ion channel systematically, we performed molecular dynamics (MD) simulations on a gramicidin A (GA)-phospholipid model system with an ion in the channel pore region. Each of the three types of ions (Ca2+, Na+ Cl-) was placed at five different positions along the channel axis by replacing a water molecule. MD simulations were performed on each system at constant pressure and constant temperature. The MD trajectories showed that the Ca2+ and Na+ ions could stably fluctuate in the pore region, but the Cl- ion was pushed out because of the unfavorable interaction with the channel. This result is consistent with experimental data. It was also found that the conformation of the GA channel underwent a significant change due to the presence of the ion, and the two ends of the GA monomer were more flexible than its middle region. In particular, the dramatic change of local pore radius near the ion indicated this kind of deformation. The strong interaction between the ion and carbonyl oxygen atoms of GA was the major contributor to this change. Furthermore, it was found that the ethanolamine group of the GA molecule was the most flexible group in the GA channel and often observed to block the entrance of GA. These results imply that the deformation of channel structure plays a very important factor in ion permeation, and the ethanolamine group may play a key role in regulating ion entry into the pore. In conclusion, our results indicate that the ion has a dominant influence on the structure of the GA channel and that the flexibility of the ion channel is a crucial factor in the ion permeation process.     

Potassium channel, ions, and water: simulation studies based on the high resolution X-ray structure of KcsA

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.     

Molecular ordering of interfacially localized tryptophan analogs in ester- and ether-lipid bilayers studied by 2H-NMR.

Perdeuterated indole-d6 and N-methylated indole-d6 were solubilized in lamellar liquid crystalline phases composed of either 1,2-diacyl-glycero-3-phosphocholine (14:0)/water or 1,2-dialkyl-glycero-3-phosphocholine(14:0/water. The molecular ordering of the tryptophan analogs was determined from deuteron quadrupole splittings observed in 2H-NMR spectra on macroscopically aligned lipid bilayers. NMR spectra were recorded with the bilayers oriented perpendicular to or parallel with the external magnetic field, and the values of the splittings differed by a factor of 2 between these distinct orientations, indicating fast rotational motion of the molecules about an axis parallel to the bilayer normal. In all cases the splittings were found to decrease with increasing temperature. Relatively large splittings were observed in all systems, demonstrating that the tryptophans partition into a highly anisotropic environment. Solubilization most likely occurs at the lipid/water interface, as indicated by 1H-NMR chemical shift studies. The 2H-NMR spectra obtained for each analog were found to be rather similar in ester and ether lipids, but with smaller splittings in the ether lipid under similar conditions. The difference was slightly less for the indole molecule. Furthermore, in both lipid systems the positions of the splittings from indole were different from those of N-methyl indole. The results suggest that 1) the tryptophan analogs are solubilized in the interfacial region of the lipid bilayer, 2) the behavior may be modulated by hydrogen bonding in the case of indole, and 3) hydrogen bonding with the lipid carbonyl groups is not likely to play a major role in the solubilization of single indole molecules in the ester lipid bilayer interface.     

The predominant role of coordination number in potassium channel selectivity

Potassium channels are exquisitely selective, allowing K+ to pass across cell membranes while blocking other ion types. Here we demonstrate that the number of carbonyl oxygen atoms that surround permeating ions is the most important factor in determining ion selectivity rather than the size of the pore or the strength of the coordinating dipoles. Although the electrostatic properties of the coordinating ligands can lead to Na+ or K+ selectivity at some values of the dipole moment, no significant selectivity arises at the specific value of the dipole moment for carbonyl groups found in potassium channels when the ligands have complete freedom. Rather, we show that the main contribution to selectivity arises from slight constraints on the conformational freedom of the channel protein that limit the number of carbonyl oxygen atoms to a value better suited to K+ than Na+, despite the pore being flexible. This mechanism provides an example of a general framework for explaining ion discrimination in a range of natural and synthetic macromolecules in which selectivity is controlled by the number of coordinating ligands in addition to their dipole moment.     

Cholesterol-lipid interactions: An infrared and raman spectroscopic study of the carbonyl stretching mode region of 1,2-dipalmitoyl phosphatidylcholine bilayers

Infrared and Raman spectra were obtained for the 1690–1770 cm^?1 carbonyl stretching mode region for 1,2-dipalmitoyl phosphatidylcholine (DPPC) bilayers in the anhydrous, partially hydrated and completely hydrated states. Spectral features at approx. 1740 and 1721 cm^?1 are assigned to CO stretching modes associated with the 1- and 2-chain carbonyl groups, respectively. Splittings of the primary transitions at 1743, 1738, ～1731 and ～1721 cm^?1 are attributed to rotational isomers involving the entire chain. Hydrogen bond formation between the fatty acid carbonyl and 3βOH cholesterol groups was investigated for anhydrous DPPC bilayers. Examination of frequencies, intensities and half-widths of the carbonyl bands indicates that no hydrogen bonding occurs at either of the two carbonyl sites. However, the addition of cholesterol to completely hydrated DPPC dispersions reduces the conformational inequivalence between the two fatty acid carbonyl groups by specifically perturbing the 2-chain. For cholesterol containing systems the carbonyl stretching mode transitions were also used to monitor lattice effects within the interface region as water binds to the bilayer head groups. Specifically, the addition of approx. 2 molecules of water per lipid molecule orders the lipid lattice and increases the bilayer packing density, while the subsequent addition of 4 molecules of water per lipid molecule releases the packing constraints within the interface region and thereby decreases the packing density.     

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.     

Changes in phosphatidylcholine headgroup tilt and water order induced by monovalent salts: molecular dynamics simulations

The association between monovalent salts and neutral lipid bilayers is known to influence global bilayer structural properties such as headgroup conformational fluctuations and the dipole potential. The local influence of the ions, however, has been unknown due to limited structural resolution of experimental methods. Molecular dynamics simulations are used here to elucidate local structural rearrangements upon association of a series of monovalent Na(+) salts to a palmitoyl-oleoyl-phosphatidylcholine bilayer. We observe association of all ion types in the interfacial region. Larger anions, which are meant to rationalize data regarding a Hofmeister series of anions, bind more deeply within the bilayer than either Cl(-) or Na(+). Although the simulations are able to reproduce experimentally measured quantities, the analysis is focused on local properties currently invisible to experiments, which may be critical to biological systems. As such, for all ion types, including Cl(-), we show local ion-induced perturbations to headgroup tilt, the extent and direction of which is sensitive to ion charge and size. Additionally, we report salt-induced ordering of the water well beyond the interfacial region, which may be significant in terms of hydration repulsion between stacked bilayers.     

Gramicidin tryptophans mediate formamidinium-induced channel stabilization.

Compared with alkali metal cations, formamidinium ions stabilize the gramicidin A channel molecule in monoolein bilayers (Seoh and Busath, 1993a). A similar effect is observed with N-acetyl gramicidin channel molecules in spite of the modified forces at the dimeric junction (Seoh and Busath, 1993b). Here we use electrophysiological measurements with tryptophan-to-phenylalanine-substituted gramicidin analogs to show that the formamidinium-induced channel molecule stabilization is eliminated when the four gramicidin tryptophans are replaced with phenylalanines in gramicidin M-. This suggests that the stabilization is mediated by the tryptophan side chains. Tryptophan residues 9, 13, and 15 must cooperate to produce the effect because replacement of any one of the three with phenylalanine significantly reduces stabilization; replacement of Trp-11 with phenylalanine causes negligible decrease in stabilization. In addition, formamidinium-related current-voltage supralinearity and open-channel noise are absent with gramicidin M-. When the lipid bilayer was formed with monoolein ether rather than monoolein ester, the channel lifetimes were reduced markedly and, at low voltage and relative to those in KCl solution, were decreased by a factor of 2, whereas the open-channel noise was unaffected and the current-voltage relation was only modestly affected. These results suggest that formamidinium modifies the state of the tryptophan side chains, which, in turn, affects channel lifetime, current-voltage supralinearity, and open-channel noise through interactions with water or lipid headgroup atoms including the lipid ester carbonyl.     

Ion leakage through transient water pores in protein-free lipid membranes driven by transmembrane ionic charge imbalance

We have employed atomic-scale molecular dynamics simulations to address ion leakage through transient water pores in protein-free phospholipid membranes. Our results for phospholipid membranes in aqueous solution with NaCl and KCl salts show that the formation of transient water pores and the consequent ion leakage can be induced and be driven by a transmembrane ionic charge imbalance, an inherent feature in living cells. These processes take place if the gradient is large enough to develop a sufficiently significant potential difference across the membrane. The transport of cations and anions through the water pores is then seen; it discharges the transmembrane potential, considerably reduces the size of a water pore, and makes the water pore metastable, leading eventually to its sealing. The ion transport is found to be sensitive to the type of ions. It turns out that Na(+) and Cl(-) ions leak through a membrane at approximately the same ratio despite the fact that Na(+) ions are expected to experience a lower potential barrier for the permeation through the pore. This is because of strong interactions of sodium ions with the carbonyl region of a phospholipid membrane as well as with lipid headgroups forming pore "walls," considerably slowing down the permeation of sodium ions. In contrast, we observed a pronounced selectivity of a phospholipid membrane to the permeation of potassium ions as compared to chloride ions: Potassium ions, being larger than sodium ions, interact only weakly with phospholipid headgroups, so that these interactions are not able to compensate for a large difference in free-energy barriers for permeation of K(+) and Cl(-) ions. These findings are found to be robust to a choice of force-field parameters for ions (tested by Gromacs and Charmm force-fields for ions). What is more, a potassium ion is found to be able to permeate a membrane along an alternate, "water-defect-mediated" pathway without actual formation of a pore. The "water-defect-mediated" leakage involves formation of a single water defect only and is found to be at least one order of magnitude faster than the pore-mediated ion leakage.     

Effect of Coordinated Ions on Structure and Flexibility of Parallel G-quadruplexes: A Molecular Dynamics Study

Abstract Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na(+) or K(+) ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 ? from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na(+) counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na(+) ions, within the quadruplex cavity, are more mobile than coordinated K(+) ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.