A nanosecond molecular dynamics study of antiparallel d(G)7 quadruplex structures: effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G7) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of approximately 1.5 A from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.     

The Molecular Mechanism of Ion-Dependent Gating in Secondary Transporters

LeuT-like fold Na-dependent secondary active transporters form a large family of integral membrane proteins that transport various substrates against their concentration gradient across lipid membranes, using the free energy stored in the downhill concentration gradient of sodium ions. These transporters play an active role in synaptic transmission, the delivery of key nutrients, and the maintenance of osmotic pressure inside the cell. It is generally believed that binding of an ion and/or a substrate drives the conformational dynamics of the transporter. However, the exact mechanism for converting ion binding into useful work has yet to be established. Using a multi-dimensional path sampling (string-method) followed by all-atom free energy simulations, we established the principal thermodynamic and kinetic components governing the ion-dependent conformational dynamics of a LeuT-like fold transporter, the sodium/benzyl-hydantoin symporter Mhp1, for an entire conformational cycle. We found that inward-facing and outward-facing states of Mhp1 display nearly the same free energies with an ion absent from the Na2 site conserved across the LeuT-like fold transporters. The barrier separating an apo-state from inward-facing or outward-facing states of the transporter is very low, suggesting stochastic gating in the absence of ion/substrate bound. In contrast, the binding of a Na2 ion shifts the free energy stabilizing the outward-facing state and promoting substrate binding. Our results indicate that ion binding to the Na2 site may also play a key role in the intracellular thin gate dynamics modulation by altering its interactions with the transmembrane helix 5 (TM5). The Potential of Mean Force (PMF) computations for a substrate entrance displays two energy minima that correspond to the locations of the main binding site S1 and proposed allosteric S2 binding site. However, it was found that substrate''s binds to the site S1 ∼5 kcal/mol more favorable than that to the site S2 for all studied bound combinations of ions and a substrate.     

A Nanosecond Molecular Dynamics Study of Antiparallel d(G)7 Quadruplex Structures: Effect of the Coordinated Cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ～ 1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.     

Computational study of the binding affinity and selectivity of the bacterial ammonium transporter AmtB

We report results from microscopic molecular dynamics and free energy perturbation simulations of substrate binding and selectivity for the Escherichia coli high-affinity ammonium transporter AmtB. The simulation system consists of the protein embedded in a model membrane/water surrounding. The calculated absolute binding free energies for the external NH(4)(+) ions are between -5.8 and -7.3 kcal/mol and are in close agreement with experimental data. The apparent pK(a) of the bound NH(4)(+) increases by more than 4 units, indicating a preference for binding ammonium ion and not neutral ammonia. The external binding site is also selective for NH(4)(+) toward monovalent metal cations by 2.4-4.4 kcal/mol. The externally bound NH(4)(+) shows strong electrostatic interactions with the proximal buried Asp160, stabilized in the anionic form, whereas the interactions with the aromatic rings of Phe107 and Trp148, lining the binding cavity, are less pronounced. Simulated mutation of the highly conserved Asp160 to Asn reduces the pK(a) of the bound ammonium ion by approximately 7 units and causes loss of its binding. The calculations further predict that the substrate affinity of E. coli AmtB depends on the ionization state of external histidines. The computed free energies of hypothetical intermediate states related to transfer of NH(3), NH(4)(+), or H(2)O from the external binding site to the first position inside the internal channel pore favor permeation of the neutral species through the channel interior. However, the predicted change in the apparent pK(a) of NH(4)(+) upon translocation from the external site, Am1, to the first internal site, Am2, indicates that ammonium ion becomes deprotonated only when it enters the channel interior.     

A compact RNA tertiary structure contains a buried backbone-K+ complex

The structure of a 58 nucleotide ribosomal RNA fragment buries several phosphate groups of a hairpin loop within a large tertiary core. During refinement of an X-ray crystal structure containing this RNA, a potassium ion was found to be contacted by six oxygen atoms from the buried phosphate groups; the ion is contained completely within the solvent-accessible surface of the RNA. The electrostatic potential at the ion chelation site is unusually large, and more than compensates for the substantial energetic penalties associated with partial dehydration of the ion and displacement of delocalized ions. The very large predicted binding free energy, approximately -30 kcal/mol, implies that the site must be occupied for the RNA to fold. These findings agree with previous studies of the ion-dependent folding of tertiary structure in this RNA, which concluded that a monovalent ion was bound in a partially dehydrated environment where Mg2+ could not easily compete for binding. By compensating the unfavorable free energy of buried phosphate groups with a chelated ion, the RNA is able to create a larger and more complex tertiary fold than would be possible otherwise.     

Ion channels in icosahedral virus: a comparative analysis of the structures and binding sites at their fivefold axes.

An analysis of the crystallographically determined structures of the icosahedral protein coats of Tomato Bushy Stunt Virus, Southern Bean Mosaic Virus, Satellite Tobacco Necrosis Virus, Human Rhinovirus 14 and Mengovirus around their fivefold axes is presented. Accessibilities surfaces, electrostatic energy profile calculations, ion-protein interaction energy calculations, free energy perturbation methods and comparisons with structures of chelating agents are used in this study. It is concluded that the structures built around the viral fivefold axes would be adequate for ion binding and transport. Relative ion preferences are derived for the binding sites, using free energy perturbation methods, which are consistent with the experimental data when available. In the cases where crystallographic studies determined the existence of ions on the fivefold axes, our results indicate that they would correspond to ions in crystallization or purification buffers. The environment of the fivefold axes are rich in polar residues in all icosahedral viral structures whose atomic coordinates are available, including some that are not being analyzed in detail in this work. The fivefold channel-like structures have most of the basic properties expected for real ion channels including a funnel at the entrance, a polar internal environment with frequent alternation of acidic and basic residues, ion binding sites, the capability to induce ion dehydration and ion transit from the external viral surface to the binding sites.     

Free energy simulations of ligand binding to the aspartate transporter Glt(Ph)

Glutamate/Aspartate transporters cotransport three Na⁺ and one H⁺ ions with the substrate and countertransport one K⁺ ion. The binding sites for the substrate and two Na⁺ ions have been observed in the crystal structure of the archeal homolog Glt_Ph, while the binding site for the third Na⁺ ion has been proposed from computational studies and confirmed by experiments. Here we perform detailed free energy simulations of Glt_Ph, giving a comprehensive characterization of the substrate and ion binding sites, and calculating their binding free energies in various configurations. Our results show unequivocally that the substrate binds after the binding of two Na⁺ ions. They also shed light into Asp/Glu selectivity of Glt_Ph, which is not observed in eukaryotic glutamate transporters.     

The structural basis of secondary active transport mechanisms

Secondary active transporters couple the free energy of the electrochemical potential of one solute to the transmembrane movement of another. As a basic mechanistic explanation for their transport function the model of alternating access was put forward more than 40 years ago, and has been supported by numerous kinetic, biochemical and biophysical studies. According to this model, the transporter exposes its substrate binding site(s) to one side of the membrane or the other during transport catalysis, requiring a substantial conformational change of the carrier protein. In the light of recent structural data for a number of secondary transport proteins, we analyze the model of alternating access in more detail, and correlate it with specific structural and chemical properties of the transporters, such as their assignment to different functional states in the catalytic cycle of the respective transporter, the definition of substrate binding sites, the type of movement of the central part of the carrier harboring the substrate binding site, as well as the impact of symmetry on fold-specific conformational changes. Besides mediating the transmembrane movement of solutes, the mechanism of secondary carriers inherently involves a mechanistic coupling of substrate flux to the electrochemical potential of co-substrate ions or solutes. Mainly because of limitations in resolution of available transporter structures, this important aspect of secondary transport cannot yet be substantiated by structural data to the same extent as the conformational change aspect. We summarize the concepts of coupling in secondary transport and discuss them in the context of the available evidence for ion binding to specific sites and the impact of the ions on the conformational state of the carrier protein, which together lead to mechanistic models for coupling.     

Possible existence of ion pairs at the mouths of ion channels

An electrostatic calculation suggests that when an ion is bound near the mouth of a channel penetrating a low-dielectric membrane, a counter ion may form an ion pair with this ion. The tendency towards ion-pair formation is remarkably enhanced at channel mouths by forces (image forces) arising from the charges induced on the boundaries between different dielectrics. The binding constant for the formation of ion-pairs of monovalent ions is estimated under the assumption that local interactions between the counter ion and the channel wall are negligibly small. It is of the order of 1–10 molal^?1 or more for the binding of a Cl^? (F^?) counter ion to an $\text{Na}{}^{\mathrm{+}}$ ($\text{Li}{}^{\mathrm{+}}$) ion if appropriate conditions are fulfilled. The binding constant depends on the position of the binding site, the dimensions and geometries of the channel and channel mouth, and the state of ion loading of the channel, as well as the ionic species. The present results also indicate that when cation (anion) channels have anionic (cationic) groups as integrant parts of their channel walls, interactions between these charged groups and permeant ions are markedly enhanced by the image forces.     

A "moving metal mechanism" for substrate cleavage by the DNA repair endonuclease APE-1

Apurinic/apyrimidinic endonuclease (APE-1) is essential for base excision repair (BER) of damaged DNA. Here molecular dynamics (MD) simulations of APE1 complexed with cleaved and uncleaved damaged DNA were used to determine the role and position of the metal ion(s) in the active site before and after DNA cleavage. The simulations started from an energy minimized wild-type structure of the metal-free APE1/damaged-DNA complex (1DE8). A grid search with one Mg2+ ion located two low energy clusters of Mg2+ consistent with the experimentally determined metal ion positions. At the start of the longer MD simulations, Mg2+ ions were placed at different positions as seen in the crystal structures and the movement of the ion was followed over the course of the trajectory. Our analysis suggests a "moving metal mechanism" in which one Mg2+ ion moves from the B- (more buried) to the A-site during substrate cleavage. The anticipated inversion of the phosphate oxygens occurs during the in-line cleavage reaction. Experimental results, which show competition between Ca2+ and Mg2+ for catalyzing the reaction, and high concentrations of Mg2+ are inhibitory, indicate that both sites cannot be simultaneously occupied for maximal activity.     

Properties of microsolvated ions: from the microenvironment of chromophore and alkali metal ions in proteins to negative ions in water clusters

Here we discuss the fascinating chemistry and physics of microsolvated ions that bridge the transition from bare ions in gas phase to ions in solution. Such ions occur in many situations in biochemistry and are crucial for several functions; metal ions, for example, must remove their water shell to pass through ion pumps in membranes. Furthermore, only a few water molecules are buried in the hydrophobic pockets of proteins where they are bound to charged amino acid residues or ionic chromophores. Another aspect is the reactivity of microsolvated ions and the importance in atmospheric, organic and inorganic chemistry. We close by a discussion of the stability of molecular dianions, and how hydration affects the electronic binding energy. There is a vast literature on microsolvated ions, and in this review we are far from being comprehensive, rather we mainly bring examples of our own work.     

Active transport of thallous ions by Streptococcus lactis.

Tl+ ions have been shown to mimic or compete with K+ in a number of membrane systems. We confirmed that in starved, valinomycin-treated cells of Streptococcus lactis 7962, Tl+ ions distributed themselves across the bacterial membrane in response to the potassium diffusion potential. In glucose-energized cells, however, Tl+ was taken up by a system specifically stimulated by sodium salts. The intracellular levels of Tl+ exceeded those attained by [3H]triphenylmethylphosphonium ion, a lipophilic cation which accumulates in response to the membrane potential. The uptake of Tl+ by (Na+ and glucose)-stimulated cells was strongly inhibited by potassium salts. These experiments suggest that metabolic energy is coupled to Tl+ transport by means of a high energy phosphate compound and that Tl+ ions are actively transported by a membrane carrier whose normal substrate is K+. The uptake of Tl+ is not a valid method for determining the streptococcal membrane potential.     

The effect of metal ions on the binding of ethanol to human alcohol dehydrogenase β<Subscript>2</Subscript>β<Subscript>2</Subscript>

Molecular docking simulations were performed in this study to investigate the importance of both structural and catalytic zinc ions in the human alcohol dehydrogenase beta(2)beta(2) on substrate binding. The structural zinc ion is not only important in maintaining the structural integrity of the enzyme, but also plays an important role in determining substrate binding. The replacement of the catalytic zinc ion or both catalytic and structural zinc ions with Cu(2+) results in better substrate binding affinity than with the wild-type enzyme. The width of the bottleneck formed by L116 and V294 in the substrate binding pocket plays an important role for substrate entrance. In addition, unfavorable contacts between the substrate and T48 and F93 prevent the substrate from moving too close to the metal ion. The optimal binding position occurs between 1.9 and 2.4 A from the catalytic metal ion.     

Identification of the Third Na Site and the Sequence of Extracellular Binding Events in the Glutamate Transporter

The transport cycle in the glutamate transporter (GlT) is catalyzed by the cotransport of three Na⁺ ions. However, the positions of only two of these ions (Na1 and Na2 sites) along with the substrate have been captured in the crystal structures reported for both the outward-facing and the inward-facing states of Glt_ph. Characterizing the third ion binding site (Na3) is necessary for structure-function studies attempting to investigate the mechanism of transport in GlTs at an atomic level, particularly for the determination of the sequence of the binding events during the transport cycle. In this study, we report a series of molecular dynamics simulations performed on various bound states of Glt_ph (the apo state, as well as in the presence of Na⁺, the substrate, or both), which have been used to identify a putative Na3 site. The calculated trajectories have been used to determine the water accessibility of potential ion-binding residues in the protein, as a prerequisite for their ion binding. Combined with conformational analysis of the key regions in the protein in different bound states and several additional independent simulations in which a Na⁺ ion was randomly introduced to the interior of the transporter, we have been able to characterize a putative Na3 site and propose a plausible binding sequence for the substrate and the three Na⁺ ions to the transporter during the extracellular half of the transport cycle. The proposed Na3 site is formed by a set of highly conserved residues, namely, Asp³¹², Thr⁹², and Asn³¹⁰, along with a water molecule. Simulation of a fully bound state, including the substrate and the three Na⁺ ions, reveals a stable structure—showing closer agreement to the crystal structure when compared to previous models lacking an ion in the putative Na3 site. The proposed sequence of binding events is in agreement with recent experimental models suggesting that two Na⁺ ions bind before the substrate, and one after that. Our results, however, provide additional information about the sites involved in these binding events.     

Transitions between Closed and Open Conformations of TolC: The Effects of Ions in Simulations

Bacteria, such as Escherichia coli, use multidrug efflux pumps to export toxic substrates through their cell membranes. Upon formation of an efflux pump, the aperture of its outer membrane protein TolC opens and thereby enables the extrusion of substrate molecules. The specialty of TolC is its ability to dock to different transporters, making it a highly versatile export protein. Within this study, the transition between two conformations of TolC that are both available as crystal structures was investigated using all-atom molecular dynamics simulations. To create a partially open conformation from a closed one, the stability of the periplasmic aperture was weakened by a double point mutation at the constricting ring, which removes some salt bridges and hydrogen bonds. These mutants, which showed partial opening in previous experiments, did not spontaneously open during a 20-ns equilibration at physiological values of the KCl solution. Detailed analysis of the constricting ring revealed that the cations of the solvent were able to constitute ionic bonds in place of the removed salt bridges, which inhibited the opening of the aperture in simulations. To remove the ions from these binding positions within the available simulation time, an extra force was applied onto the ions. To keep the effect of this additional force rather flexible, it was applied in form of an artificial external electric field perpendicular to the membrane. Depending on the field direction and the ion concentration, these simulations led to a partial opening. In experiments, this energy barrier for the ions can be overcome by thermal fluctuations on a longer timescale.     

Crystal structure of the bromide-bound D85S mutant of bacteriorhodopsin: principles of ion pumping

We report the crystal structure of a bromide-bound form of the D85S mutant of bacteriorhodopsin, bR(D85S), a protein that uses light energy rather than ATP to pump halide ions across the cell membrane. Comparison of the structure of the halide-bound and halide-free states reveals that both displacements of individual side-chain positions and concerted helical movements occur on the extracellular side of the protein. Analysis of these structural changes reveals how this ion pump first facilitates ion uptake deep within the cell membrane and then prevents the backward escape of ions later in the pumping cycle. Together with the information provided by structures of intermediate states in the bacteriorhodopsin photocycle, this study also suggests the overall design principles that are necessary for ion pumping.     

Probing metal ion substrate-binding to the E. coli ZitB exporter in native membranes by solid state NMR

Metal ion homeostasis is important for healthy cell function and is regulated by metal ion transporters and chaperones. To explore metal ion binding to membrane transport proteins we have used cadmium-113 as a solid state NMR probe of the Escherichia coli zinc exporter ZitB present in native membrane preparations. Competition experiments with other metal ions indicated that nickel and copper are also able to bind to this protein. Metal ion uptake studies were also performed using ZitB-reconstituted into proteoliposomes for a well established fluorescence assay. The results of both the solid state NMR and the uptake studies demonstrate that ZitB is potentially capable of transporting not only zinc but also cadmium, nickel and copper. The solid state NMR approach therefore offers great potential for defining the substrate spectrum of metal ion transporter proteins in their native membrane environments. Further, it should be useful for functional dissection of transporter mechanisms by facilitating the identification of functional residues by mutational studies.     

Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters

The recent determination of high-resolution crystal structures of several transporters offers unprecedented insights into the structural mechanisms behind secondary transport. These proteins utilize the facilitated diffusion of the ions down their electrochemical gradients to transport the substrate against its concentration gradient. The structural studies revealed striking similarities in the structural organization of ion and solute binding sites and a well-conserved inverted-repeat topology between proteins from several gene families. In this paper we will overview recent atomistic simulations applied to study the mechanisms of selective binding of ion and substrate in LeuT, Glt, vSGLT and hSERT as well as its consequences for the transporter conformational dynamics. This article is part of a Special Issue entitled: Membrane protein structure and function.     

DNA cleavage by EcoRV endonuclease: two metal ions in three metal ion binding sites

Four crystal structures of EcoRV endonuclease mutants K92A and K38A provide new insight into the mechanism of DNA bending and the structural basis for metal-dependent phosphodiester bond cleavage. The removal of a key active site positive charge in the uncleaved K92A-DNA-M(2+) substrate complex results in binding of a sodium ion in the position of the amine nitrogen, suggesting a key role for a positive charge at this position in stabilizing the sharp DNA bend prior to cleavage. By contrast, two structures of K38A cocrystallized with DNA and Mn(2+) ions in different lattice environments reveal cleaved product complexes featuring a common, novel conformation of the scissile phosphate group as compared to all previous EcoRV structures. In these structures, the released 5'-phosphate and 3'-OH groups remain in close juxtaposition with each other and with two Mn(2+) ions that bridge the conserved active site carboxylates. The scissile phosphates are found midway between their positions in the prereactive substrate and postreactive product complexes of the wild-type enzyme. Mn(2+) ions occupy two of the three sites previously described in the prereactive complexes and are plausibly positioned to generate the nucleophilic hydroxide ion, to compensate for the incipient additional negative charge in the transition state, and to ionize a second water for protonation of the 3'-oxyanion. Reconciliation of these findings with earlier X-ray and fluorescence studies suggests a novel mechanism in which a single initially bound metal ion in a third distinct site undergoes a shift in position together with movement of the scissile phosphate deeper into the active site cleft. This reconfigures the local environment to permit binding of the second metal ion followed by movement toward the pentacovalent transition state. The new mechanism suggested here embodies key features of previously proposed two- and three-metal catalytic models, and offers a view of the stereochemical pathway that integrates much of the copious structural and functional data that are available from exhaustive studies in many laboratories.     

Mechanism of Cation Binding to the Glutamate Transporter EAAC1 Probed with Mutation of the Conserved Amino Acid Residue Thr101

The glutamate transporter excitatory amino acid carrier 1 (EAAC1) catalyzes the co-transport of three Na⁺ ions, one H⁺ ion, and one glutamate molecule into the cell, in exchange for one K⁺ ion. Na⁺ binding to the glutamate-free form of the transporter generates a high affinity binding site for glutamate and is thus required for transport. Moreover, sodium binding to the transporters induces a basal anion conductance, which is further activated by glutamate. Here, we used the [Na⁺] dependence of this conductance as a read-out of Na⁺ binding to the substrate-free transporter to study the impact of a highly conserved amino acid residue, Thr¹⁰¹, in transmembrane domain 3. The apparent affinity of substrate-free EAAC1 for Na⁺ was dramatically decreased by the T101A but not by the T101S mutation. Interestingly, in further contrast to EAAC1_WT, in the T101A mutant this [Na⁺] dependence was biphasic. This behavior can be explained by assuming that the binding of two Na⁺ ions prior to glutamate binding is required to generate a high affinity substrate binding site. In contrast to the dramatic effect of the T101A mutation on Na⁺ binding, other properties of the transporter, such as its ability to transport glutamate, were impaired but not eliminated. Our results are consistent with the existence of a cation binding site deeply buried in the membrane and involving interactions with the side chain oxygens of Thr¹⁰¹ and Asp³⁶⁷. A theoretical valence screening approach confirms that the predicted site of cation interaction has the potential to be a novel, so far undetected sodium binding site.