Comment on "Free energy simulations of single and double ion occupancy in gramicidin A" [J. Chem. Phys. 126, 105103 (2007)

In a recent article published by Bastug and Kuyucak [J. Chem. Phys.126, 105103 (2007)] investigated the microscopic factors affecting double ion occupancy in the gramicidin channel. The analysis relied largely on the one-dimensional potential of mean force of ions along the axis of the channel (the so-called free energy profile of the ion along the channel axis), as well as on the calculation of the equilibrium association constant of the ions in the channel binding sites. It is the purpose of this communication to clarify this issue.     

短杆菌肽在十八烷基硫醇自组装单层膜上的通道行为

短杆菌肽是短杆芽孢菌在芽孢形成过程中合成的一组憎水十五肽,它能在脂双层膜上形成一价阳离子通道,允许一些一价阳离子通过．短杆菌肽离子通道研究在天然双层膜以及人工双层膜,如类脂双层膜（BLM）、脂质体上开展的较多,而基于支撑膜的研究则极少［1-4]．本文在金电极上制备硫醇自组装单层膜,研究了短杆菌肽的离子通道行为,并提出了可能的机理．     

Structural and functional consequences of an amide-to-ester substitution in the selectivity filter of a potassium channel

The selectivity filter of K(+) channels comprises four contiguous ion binding sites, S1 through S4. Structural and functional data indicate that the filter contains on average two K(+) ions at any given time and that these ions reside primarily in two configurations, namely to sites S1 and S3 or to sites S2 and S4. Maximum ion flux through the channel is expected to occur when the energy difference between these two binding configurations is zero. In this study, we have used protein semisynthesis to selectively perturb site 1 within the filter of the KcsA channel through use of an amide-to-ester substitution. The modification alters K(+) conduction properties. The structure of the selectivity filter is largely unperturbed by the modification, despite the loss of an ordered water molecule normally located just behind the filter. Introduction of the ester moiety was found to alter the distribution of K(+), Rb(+,) and Cs(+) within the filter, with the most dramatic change found for Rb(+). The redistribution of ions is associated with the appearance of a partially hydrated ion just external to the filter, at a position where no ion is observed in the wild-type channel. The appearance of this new ion-binding site creates a change in the distance between a pair of K(+) ions some fraction of the time, apparently leading to a reduction in the ion conduction rate. Importantly, this finding suggests that the selectivity filter of a potassium channel is optimized both in terms of absolute ion occupancy and in terms of the separation in distance between the conducting ions.     

Wavelength-selective fluorescence in ion channels formed by gramicidin A in membranes

Gramicidins are linear peptides that form ion channels that are specific for monovalent cations in membranes. The tryptophan residues in the gramicidin channel play a crucial role in the organization and function of the channel. The natural mixture of gramicidins, denoted as gramicidin A’, consists of mostly gramicidin A, but also contains gramicidins B, C and D as minor components. We have previously shown that the tryptophan residues in ion channels formed by the naturally occurring peptide, gramicidin A’, display wavelength-dependent fluorescence characteristics due to the motionally restricted environment in which they are localized. In order to check the influence of ground-state heterogeneity in the observed wavelength-selective fluorescence of gramicidin A’ in membranes, we performed similar experiments with pure gramicidin A in model membranes. Our results show that the observed wavelength-selective fluorescence characteristics of naturally occurring gramicidin A’ are not due to ground-state heterogeneity.     

(De)constructing the ryanodine receptor: modeling ion permeation and selectivity of the calcium release channel

Biological ion channels are proteins that passively conduct ions across membranes that are otherwise impermeable to ions. Here, we present a model of ion permeation and selectivity through a single, open ryanodine receptor (RyR) ion channel. Combining recent mutation data with electrodiffusion of finite-sized ions, the model reproduces the current/voltage curves of cardiac RyR (RyR2) in KCl, LiCl, NaCl, RbCl, CsCl, CaCl(2), MgCl(2), and their mixtures over large concentrations and applied voltage ranges. It also reproduces the reduced K(+) conductances and Ca(2+) selectivity of two skeletal muscle RyR (RyR1) mutants (D4899N and E4900Q). The model suggests that the selectivity filter of RyR contains the negatively charged residue D4899 that dominates the permeation and selectivity properties and gives RyR a DDDD locus similar to the EEEE locus of the L-type calcium channel. In contrast to previously applied barrier models, the current model describes RyR as a multi-ion channel with approximately three monovalent cations in the selectivity filter at all times. Reasons for the contradicting occupancy predictions are discussed. In addition, the model predicted an anomalous mole fraction effect for Na(+)/Cs(+) mixtures, which was later verified by experiment. Combining these results, the binding selectivity of RyR appears to be driven by the same charge/space competition mechanism of other highly charged channels.     

Not all G-quadruplexes exhibit ion-channel-like properties: NMR study of ammonium ion (non)movement within the d(G(3)T(4)G(4))(2) quadruplex

A solution-state NMR study on 15NH4(+) ion movement within d(G(3)T(4)G(4))(2), a dimeric G-quadruplex consisting of three G-quartets and two T(4) loops, rather unexpectedly demonstrated the absence of 15NH4(+) ion movement between the binding sites U and L along the central axis of the G-quadruplex. Distinct temperature dependences of autocorrelation signals for U and L binding sites have been observed in 15N-1H NzExHSQC spectra which correlate with the local stiffness of the G-quadruplex. The volumes of the cross-peaks, which are the result of 15NH4(+) ion movement, have been interpreted in terms of rate constants, T(1) relaxation, and proton exchange. 15NH4(+) ion movements from the binding sites U and L into the bulk solution are characterized by lifetimes of 139 ms and 1.7 s at 298 K, respectively. The 12 times faster movement from the binding site U demonstrates that 15NH4(+) ion movement is controlled by the structure of T4 loop residues, which through diagonal- vs edge-type orientations impose distinct steric restraints for cations to leave or enter the G-quadruplex. Arrhenius-type analysis has afforded an activation energy of 66 kJ mol(-)1 for the UB process, while it could not be determined for the LB process due to slow rates at temperatures below 298 K. We further the use of the 15NH4(+) ion as an NMR probe to gain insight into the occupancy of binding sites by cations and kinetics of ion movement which are intrinsically correlated with the structural details, dynamic fluctuations, and local flexibility of the DNA structure.     

Molecular dynamics simulation exploration of cooperative migration mechanism of calcium ions in sarcoplasmic reticulum Ca2+‐ATPase

Calcium ATPase is a member of the P‐type ATPase, and it pumps calcium ions from the cytoplasm into the reticulum against a concentration gradient. Several X‐ray structures of different conformations have been solved in recent years, providing basis for elucidating the active transport mechanism of Ca²⁺ ions. In this work, molecular dynamics (MD) simulations were performed at atomic level to investigate the dynamical process of calcium ions moving from the outer mouth of the protein to their binding sites. Five initial locations of Ca²⁺ ions were considered, and the simulations lasted for 2 or 6 ns, respectively. Specific pathways leading to the binding sites and large structural rearrangements around binding sites caused by uptake of calcium ions were identified. A cooperative binding mechanism was observed from our simulation. Firstly, the first Ca²⁺ ion binds to site I , and then, the second Ca²⁺ ion approaches. The interactions between the second Ca²⁺ and the residues around site I disturb the binding state of site I and weaken its binding ability for the first bound Ca²⁺. Because of the electrostatic repulsion of the second Ca²⁺ and the electrostatic attraction of site II , the first bound Ca²⁺ shifts from site I to site II . Concertedly, the second Ca²⁺ binds to site I , forming a binding state with two Ca²⁺ ions, one at site I and the other at site II . Both of Glu908 and Asp800 coordinate with the two Ca²⁺ ions simultaneously during the concerted binding process, which is believed to be the hinge to achieve the concerted binding. In our simulations, four amino acid residues that serve as the channel to link the outer mouth and the binding sites during the binding process were recognized, namely Tyr837, Tyr763, Asn911, and Ser767. The analyses regarding the activity of the proteins via mutations of some key residues also supported our cooperative mechanism. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Matrix-free laser desorption/ionization of ions landed on plasma-treated metal surfaces

We report new experiments in which laser desorption/ionization time-of-flight mass spectrometry (LDI-TOF-MS) was applied to detection and characterization of gramicidin S and IgG pentapeptide (DSDPR) that were reactively landed on plasma-treated stainless steel surfaces. The distributions of [M + H](+), [M + Na](+) and [M + K](+) ion species in LDI-TOF for gramicidin S and IgG pentapeptide (DSDPR) were found to be markedly different from those in conventional MALDI-TOF spectra of the same samples. LDI-TOF mass spectra showed a strong preference for [M + K](+) adducts even in the presence of a large excess of sodium cations, or following surface treatment with trifluoroacetic acid. Alkali metal cations (K(+) and Cs(+)) can be exchanged in reactively landed peptide samples to provide the corresponding cationized peptide ions by LDI. Multiple charged trypsin cations were reactively landed into a layer of 2-(4-hydroxyphenylazo)benzoic acid and ionized by LDI. The ionization mechanisms for LDI of surface-deposited peptides are briefly discussed. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Direct NMR detection of alkali metal ions bound to G-quadruplex DNA

We describe a general multinuclear (1H, 23Na, 87Rb) NMR approach for direct detection of alkali metal ions bound to G-quadruplex DNA. This study is motivated by our recent discovery that alkali metal ions (Na+, K+, Rb+) tightly bound to G-quadruplex DNA are actually "NMR visible" in solution (Wong, A.; Ida, R.; Wu, G. Biochem. Biophys. Res. Commun. 2005, 337, 363). Here solution and solid-state NMR methods are developed for studying ion binding to the classic G-quadruplex structures formed by three DNA oligomers: d(TG4T), d(G4T3G4), and d(G4T4G4). The present study yields the following major findings. (1) Alkali metal ions tightly bound to G-quadruplex DNA can be directly observed by NMR in solution. (2) Competitive ion binding to the G-quadruplex channel site can be directly monitored by simultaneous NMR detection of the two competing ions. (3) Na+ ions are found to locate in the diagonal T4 loop region of the G-quadruplex formed by two strands of d(G4T4G4). This is the first time that direct NMR evidence has been found for alkali metal ion binding to the diagonal T4 loop in solution. We propose that the loop Na+ ion is located above the terminal G-quartet, coordinating to four guanine O6 atoms from the terminal G-quartet and one O2 atom from a loop thymine base and one water molecule. This Na+ ion coordination is supported by quantum chemical calculations on 23Na chemical shifts. Variable-temperature 23Na NMR results have revealed that the channel and loop Na+ ions in d(G4T4G4) exhibit very different ion mobilities. The loop Na+ ions have a residence lifetime of 220 micros at 15 degrees C, whereas the residence lifetime of Na+ ions residing inside the G-quadruplex channel is 2 orders of magnitude longer. (4) We have found direct 23Na NMR evidence that mixed K+ and Na+ ions occupy the d(G4T4G4) G-quadruplex channel when both Na+ and K+ ions are present in solution. (5) The high spectral resolution observed in this study is unprecedented in solution 23Na NMR studies of biological macromolecules. Our results strongly suggest that multinuclear NMR is a viable technique for studying ion binding to G-quadruplex DNA.     

High resolution and tandem fourier-transform mass spectrometry with californium-252 plasma desorption

For ions formed by plasma desorption (PD) in a Fourier-transform mass spectrometer, high resolution measurements are demonstrated, such as 65,000 (FWHH) for the protonated molecular ion of gramicidin S (MW 1140.7). Resolution is substantially improved by delaying measurements until a significant ion concentration has built up in the cell, and by collisionally deactivating the orbital kinetic energy of the ions. This also makes the ions available for subsequent dissociation steps, so that tandem mass spectrometry can be demonstrated for PD ions. With this for larger ions, collisionally activated dissociation (CAD) is affected with> 85% efficiency. The CAD spectra of (M + Na)⁺ and of fragment ions from the PD of gramicidin S provide structurally useful information.     

自组装杂化离子通道膜的构筑及阳离子传输机制

以3-异氰酸丙基三乙氧基硅烷和对甲氧基苯胺为原料合成了一种可以自组装形成有机-无机杂化材料的化合物--3-(脲基-4-甲氧基苯基)丙基三乙氧基硅烷. 采用FT-IR, ¹H NMR, DSC 和XRD 分析方法对该化合物的结构以及结晶性进行了表征. 将该化合物与聚乙烯醇(PVA)共混, 利用化合物的自组装性质构筑结构均一且致密无孔的离子通道杂化膜, 通过自制的膜运输实验装置测定膜对阳离子的传输性能并提出了相应的传输机制. SEM 照片显示, 自组装杂化膜致密无缺陷, 膜厚度为8 μm. 选择5 种阳离子进行运输实验测试, 结果表明, 自组装杂化离子通道膜对一价的碱金属离子Li⁺, Na⁺和K⁺有很好的传输功能, 这要归功于杂化材料中甲氧基苯基与碱金属阳离子形成的阳离子-π相互作用力. 碱金属阳离子在膜中的扩散过程可由溶解-扩散机制来解释, 结果显示, Li⁺, Na⁺和K⁺在杂化膜中传输的渗透率大小为: P_Na⁺ > P_K⁺ > P_Li⁺ , 说明本研究中的的自组装杂化离子通道膜对Na⁺有优先选择性. 杂化离子通道膜对二价的Ca²⁺和Mg²⁺没有传输作用, 此结果给一二价阳离子的分离带来很好的研究思路.     

Functional dynamics of ion channels: modulation of proton movement by conformational switches

Detailed comparative studies of proton relay in native and chemically modified gramicidin channels provide a unique opportunity to uncover the structural basis of biological proton transport. The function of ion channels hinges on their ability to provide surrogate solvation in narrow pore filters so as to overcome the dielectric barrier presented by biological membranes. In the potassium channel KcsA and in the cation channel gramicidin, permeant selectivity and mobility are determined by the proteinaceous matrix via hydrogen bonding, charge-dipole, and dipole-dipole interactions. In particular, main-chain carbonyl groups in these pore interiors play an essential role in the solvation of alkali ions and of protons. In this study, molecular dynamics simulations reveal how the translocation of H(+) is controlled by nanosecond conformational transitions exchanging distorted states of the peptidic backbone in the single-file region of a dioxolane-linked analogue of the gramicidin dimer. These results underline the functional role of channel dynamics and provide a mechanism for the modulation of proton currents by fluctuating dipoles.     

How alkali metal ion binding alters the conformation preferences of gramicidin A: a molecular dynamics and ion mobility study

Here, we present a systematic study combing electrospray ionization-ion mobility experiments and an enhanced sampling molecular dynamics, specifically integrated tempering sampling molecular dynamics simulations (ITS-MDS), to explore the conformations of alkali metal ion (Na, K, and Cs) adducts of gramicidin A (GA) in vacuo. Folding simulation is performed to obtain inherent conformational preferences of neutral GA to provide insights about how the binding of metal ions influences the intrinsic conformations of GA. The comparison between conformations of neutral GA and alkali metal ion adducts reveals a high degree of structural similarity, especially between neutral GA and [GA + Na](+); however, the structural similarities decrease as ionic radius of the metal increases. Collision cross section (CCS) profiles for [GA + Na](+) and [GA + Cs](+) ions obtained from by ITS-MDS compare favorably with the experimental CCS, but there are significant differences from CCS profiles for [GA + K](+) ions. Such discrepancies between the calculated and measured CCS profiles for [GA + K](+) are discussed in terms of limitations in the simulation force field as well as possible size-dependent coordination of the [GA + K](+) ion complex.     

The effect of metallic cations on the electrophoretic properties of a pretreatedHevea latex

High-ammonia latex concentrate prepared from doubly-centrifuged fieldHevea latex was exhaustively dialysed to remove any residual water-soluble non-rubber constituents. The electrophoretic mobilities of the dialysed latex in the presence of various metallic cations were investigated as a function of electrolyte concentration. The mobility decreased with increasing concentration of the cations Na⁺, Ba²⁺, Mg²⁺ and Ca²⁺ in a manner consistent with the effect of simple electrolyte on compression of the electric double layer. Anomalous behaviour was noted for the divalent ion copper, in that it reversed the charge of the latex particles at a concentration even lower than that of uranyl ions. Multivalent cations (lanthanum, cerium and thorium) had a profound influence on the latex particles where very low charge reversal concentrations were observed. It is believed that strong adsorption of hydrolysed species from the metallic ions was responsible for reversing the charge of the originally negative latex particles. These experiments indicated that the efficiency with which the cations reversed the charge of the latex particle surface was in the order: lanthanum > cerium > copper > thorium > uranyl > calcium > magnesium > barium > sodium. The number of cation binding sites on the latex particle surface and the chemical free energies of cation adsorption were calculated. It was found that the interaction of the latex particle with the hydrolysable metallic cations was much more stronger than that with the simple divalent cations and that this intercation was comparable to that of biological surfaces.     

Steered molecular dynamics simulations of Na+ permeation across the gramicidin A channel

The potential of mean forces (PMF) governing Na+ permeation through gramicidin A (gA) channels with explicit water and membrane was characterized using steered molecular dynamics (SMD) simulations. Constant-force SMD with a steering force parallel to the channel axis revealed at least seven energy wells in each monomer of the channel dimer. Except at the channel dimer interface, each energy well is associated with at least three and at most four backbone carbonyl oxygens and two water oxygens in a pseudo-hexahedral or pseudo-octahedral coordination with the Na+ ion. Repeated constant-velocity SMD by dragging a Na+ ion from each energy well in opposite directions parallel to the channel axis allowed the computation of the PMF across the gA channel, revealing a global minimum corresponding to Na+ binding sites near the entrance of gA at +/-9.3 A from the geometric center of the channel. The effect of volatile anesthetics on the PMF was also analyzed in the presence of halothane molecules. Although the accuracy of the current PMF calculation from SMD simulations is not yet sufficient to quantify the PMF difference with and without anesthetics, the comparison of the overall PMF profiles nevertheless confirms that the anesthetics cause insignificant changes to the structural makeup of the free energy wells along the channel and the overall permeation barrier. On average, the PMF appears less rugged in the outer part of the channel in the presence of anesthetics, consistent with our earlier finding that halothane interaction with anchoring residues makes the gA channel more dynamic. A causal relationship was observed between the reorientation of the coordinating backbone carbonyl oxygen and Na+ transit from one energy well to another, suggesting the possibility that even minute changes in the conformation of pore-lining residues due to dynamic motion could be sufficient to trigger the ion permeation. Because some of the carbonyl oxygens contribute to Na+ coordination in two adjacent energy wells, our SMD results reveal that the atomic picture of ion "hopping" through a gA channel actually involves a Na+ ion being carried in a relay by the coordinating oxygens from one energy well to the next. Steered molecular dynamics complements other computational approaches as an attractive means for the atomistic interpretation of experimental permeation studies.     

Calix[6]arene‐Based Cascade Complexes of Organic Ion Triplets Stable in a Protic Solvent

Herein we report a D_3h‐symmetric tail‐to‐tail bis‐calix[6]arene 3 featuring two divergent cavities triply connected by ureido linkages. This calix[6]tube was synthesized by a domino Staudinger/aza‐Wittig reaction followed by a macrocyclization reaction. This process also afforded a C_2h‐symmetric isomer that represents a rare example of a self‐threaded rotaxane based on calix[6]arene subunits. The binding properties of 3 have been evaluated by NMR studies. Thus, bis‐calix[6]arene 3 is able to bind simultaneously two neutral ureido guests through an induced‐fit process. The guests are located in the cavities and are recognized through multiple hydrogen‐bonding interactions with the ureido bridges. Host 3 can also simultaneously bind multiple ions and is especially efficient for the complexation of organic ion triplets. The anion is recognized through hydrogen‐bonding interactions at the ureido binding site and is thus located between the two ammonium ions accommodated in the cavities. The resulting [1+1+2] quaternary complexes represent rare examples of cascade complexes with organic cations. These complexes are unique: 1) They are stable even in a markedly protic solvent, 2) the recognition of the ion triplets proceeds in a cooperative way through an induced‐fit process and with a high selectivity, linear cations and doubly charged anions being particularly well recognized, 3) the ions are bound as contact ion triplets thanks to the closeness of the three binding sites, 4) the cationic guests can be exchanged and thus mixed [1+1+1+1] complexes can be obtained, 5) the ureido linkers wrapped around the anion adopt a helical shape and the resulting chirality is sensed by the cations. In other words, bis‐calix[6]arene 3 presents a selective inner tunnel in which multiple guests such as organic ion triplets can be aligned in a cooperative way through induced‐fit processes.     

Ag+ labeling: a convenient new tool for the characterization of hydrogen-bonded supramolecular assemblies by MALDI-TOF mass spectrometry

Herein we describe our results on the characterization of a wide variety of different hydrogen-bonded assemblies by means of a novel matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) technique with Ag+ labeling. The labeling technique with Ag+ ions is extremely mild and provides a nondestructive way to generate charged assemblies that can be detected by mass spectrometry. Up to now more than 25 different single (1(3).2(3)), double (3(3).2(6)), and tetrarosettes (4(3).2(12)) have been successfully characterized by the use of this method. The success of the method entirely depends on the presence of a suitable binding site for the Ag+ ion. A variety of functionalities has been identified that provide strong binding sites for Ag+, either acting in a cooperative way (pi-arene and pi-alkene donor functionalities) or individually (cyano and crown ether functionalities). The method works well for assemblies with molecular weights between 2,000 and 8,000 Da, and most likely far beyond this limit.     

The sodium ions inside a lipophilic G-quadruplex channel as probed by solid-state (23)Na NMR

We report solid-state 23Na NMR and X-ray crystallographic results for a self-assembled G-quadruplex channel formed by a guanine nucleoside, 5'-tert-butyl-dimethylsilyl-2',3'-O-isopropylidene guanosine (G 1). The study provides an unambiguous 23Na NMR identification for the Na+ ions inside a lipophilic G-quadruplex channel. The crystalline nature of the sample yields a remarkably high resolution in the 23Na multiple-quantum magic-angle spinning (MQMAS) spectrum, making it possible to extract very accurate 23Na NMR parameters for each of the three crystallographically distinct Na sites. The observation of a single Na+ ion from a 9-kDa system demonstrates the potential of solid-state 23Na NMR as a complementary technique to X-ray for detecting Na+ ions in biological structures.