De Novo Cloning and Functional Characterization of a Mechanosensitive Piezo-Like Ion Channel in the Crayfish.

BACKGROUND/AIMS: Mechanosensitive ion channels are the principal elements in the transduction of mechanical force to neural activity. To date, considerably fewer studies have been published about the molecular and structural properties of mechanosensitive channels. Piezo channels are the only ion channel family in eukaryotes which is selectively gated by the membrane tension. Piezo channels have been described in mammals and some other eukaryotes. However, not much information is available for the crustaceans. METHODS: Conventional cloning methods were used to clone the putative PIEZO channel mRNA in crayfish ganglia samples. HEK293T cells were transfected by the plasmid of the cloned gene for functional studies. The CDS of the mRNA translated into the protein sequence and three-dimensional structure of the channel has been calculated. RESULTS: An mRNA, 9378 bp, was firstly cloned from crayfish which codes a 2674 residues protein. The cloned sequence is similar to the piezo channel mRNAs reported in the other species. The sequence of the coded protein has been analyzed, and some functional domains have been identified. A three-dimensional structure of the coded protein was successfully calculated in reference to mouse piezo 1 channel protein data. A plasmid with a fluorescent protein indicator was synthesized for heterologous expression in HEK293T cells. The evoked calcium response to mechanical stimulation was not different from those observed in the control cells. However, the transfected cells were more sensitive to the gating modifier YODA-1. CONCLUSION: Based on the apparent similarity in sequence, structure and functional properties to other known piezo channels, it has been proposed that cloned mRNA may code a piezo-like ion channel in crayfish.

A novel homozygous missense variant in TBC1D31 in a consanguineous family with congenital anomalies of the kidney and urinary tract (CAKUT).

Congenital anomalies of the kidney and urinary tract (CAKUT) is the leading cause of chronic kidney disease in the first three decades of life. Until now, more than 180 monogenic causes of isolated or syndromic CAKUT have been described. In addition, copy number variants (CNV) have also been implicated, however, all of these causative factors only explain a small fraction of patients with CAKUT, suggesting that additional yet-to-be-discovered novel genes are present. Herein, we report three siblings (two of them are monozygotic twin) of a consanguineous family with CAKUT. Whole-exome sequencing identified a homozygous variant in TBC1D31. Three dimensional protein modeling as well as molecular dynamics simulations predicted it as pathogenic. We therefore showed for the first time an association between a homozygous TBC1D31 variant with CAKUT in humans, expanding its genetic spectrum.

Mechanism and energetics of ligand release in the aspartate transporter GltPh.

The bacterial aspartate transporter GltPh cotransports three Na(+) ions with the substrate. The mechanism and energetics of ligand binding have previously been studied using molecular dynamics simulations on the crystal structure of GltPh captured in the outward-facing state. Here we use the recent crystal structure of the inward-facing state of GltPh to study the reverse process of unbinding of ligands. Gating behavior is studied in the presence of different ligands. A detailed characterization of the intracellular gate is given, pointing out the differences from the extracellular gate. We then perform free energy simulations to calculate the binding affinities of all the ligands in different combinations, from which the unbinding order is determined as Na2, (gate opens), Asp, Na1, and Na3. The strong coupling between Asp and Na1 is quantified from several free energy calculations. Na3 has the largest affinity to GltPh, and therefore, its unbinding is proposed as the rate-limiting step in the transport cycle. The release time of Na3, estimated from Kramers' rate theory, is shown to be consistent with the experimental turnover rate of the transporter.

Position of the third Na+ site in the aspartate transporter GltPh and the human glutamate transporter, EAAT1.

Glutamate transport via the human excitatory amino acid transporters is coupled to the co-transport of three Na(+) ions, one H(+) and the counter-transport of one K(+) ion. Transport by an archaeal homologue of the human glutamate transporters, Glt(Ph), whose three dimensional structure is known is also coupled to three Na(+) ions but only two Na(+) ion binding sites have been observed in the crystal structure of Glt(Ph). In order to fully utilize the Glt(Ph) structure in functional studies of the human glutamate transporters, it is essential to understand the transport mechanism of Glt(Ph) and accurately determine the number and location of Na(+) ions coupled to transport. Several sites have been proposed for the binding of a third Na(+) ion from electrostatic calculations and molecular dynamics simulations. In this study, we have performed detailed free energy simulations for Glt(Ph) and reveal a new site for the third Na(+) ion involving the side chains of Threonine 92, Serine 93, Asparagine 310, Aspartate 312, and the backbone of Tyrosine 89. We have also studied the transport properties of alanine mutants of the coordinating residues Threonine 92 and Serine 93 in Glt(Ph), and the corresponding residues in a human glutamate transporter, EAAT1. The mutant transporters have reduced affinity for Na(+) compared to their wild type counterparts. These results confirm that Threonine 92 and Serine 93 are involved in the coordination of the third Na(+) ion in Glt(Ph) and EAAT1.

Molecular dynamics simulations of membrane proteins.

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins-ion channels and transporters-which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.

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.

Comparative study of the energetics of ion permeation in Kv1.2 and KcsA potassium channels.

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K(+) ions in various configurations. For both channels, we find that a pair of K(+) ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K(+) ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics.

Importance of the peptide backbone description in modeling the selectivity filter in potassium channels.

A dihedral energy correction (CMAP) term has been recently included in the CHARMM force field to obtain a more accurate description of the peptide backbone. Its importance in improving dynamical properties of proteins and preserving their stability in long molecular-dynamics simulations has been established for several globular proteins. Here we investigate its role in maintaining the structure and function of two potassium channels, Shaker K(v)1.2 and KcsA, by performing molecular-dynamics simulations with and without the CMAP correction in otherwise identical systems. We show that without CMAP, it is not possible to maintain the experimentally observed orientations of the carbonyl groups in the selectivity filter in Shaker, and the channel loses its selectivity property. In the case of KcsA, the channel retains some selectivity even without CMAP because the carbonyl orientations are relatively better preserved compared to Shaker.

Potential of mean force calculations of ligand binding to ion channels from Jarzynski's equality and umbrella sampling.

Potential of mean force (PMF) calculations provide a reliable method for determination of the absolute binding free energies for protein-ligand systems. The common method used for this purpose -- umbrella sampling with weighted histogram analysis -- is computationally very laborious, which limits its applications. Recently, a much simpler alternative for PMF calculations has become available, namely, using Jarzynski's equality in steered molecular dynamics simulations. So far, there have been a few comparisons of the two methods and mostly in simple systems that do not reflect the complexities of protein-ligand systems. Here, we use both methods to calculate the PMF for ion permeation and ligand binding to ion channels. Comparison of results indicate that Jarzynski's method suffers from relaxation problems in complex systems and would require much longer simulation times to yield reliable PMFs for protein-ligand systems.

Binding of organic cations to gramicidin A channel studied with AutoDock and molecular dynamics simulations.

The accurate description of protein-ligand binding energies and configurations is an important problem in molecular biology with many applications in medicine and pharmacology. Molecular dynamics (MD) simulations provide the required accuracy but they are too slow for searching binding positions. Conversely, docking methods are much faster but have limited accuracy. An appropriate combination of the two methods could avoid the shortcomings associated with each, thus offering a better approach to the protein-ligand binding problem. Here we investigate the feasibility of such a combined docking-MD approach in a well-defined system: binding of organic cations to the gramicidin A channel. We use the AutoDock program to generate a set of protein-ligand binding configurations, which are then refined in MD simulations. For each system, we examine the binding configuration in detail and calculate the binding free energy by constructing the potential of mean force for the ligand. Our results show that AutoDock provides suitable initial configurations, which can be used profitably in MD simulations to obtain an accurate description of protein-ligand binding with a reasonable computational effort.

Free energy simulations of single and double ion occupancy in gramicidin A.

Simultaneous occupancy of the two binding sites in gramicidin A by monovalent cations is a well known property of this channel, but the energetic feasibility of this process in molecular dynamics simulations has not been established so far. Here the authors study the energetics of single and double ion occupancy in gramicidin A by constructing the potential of mean force for single and pair of cations. As representatives of small and large ions, they consider both Na+ and K+ ions in the calculations. Binding constants of ions are estimated from the free energy profiles. Comparisons with the experimental results indicate 3-4 kT discrepancy in the binding energies. They also study the coordination of the ions in their respective binding sites and the dynamic behavior of the channel water during the double ion binding process.

Molecular dynamics simulations of gramicidin A in a lipid bilayer: from structure-function relations to force fields.

Molecular dynamics simulations of membrane proteins have become a popular tool for studying their dynamic features, which are not easily accessible by experiments. Whether the force fields developed for globular proteins are adequate this purpose is an important question that is often glossed over. Here we determine the permeation properties of potassium ions in the gramicidin A channel in a lipid bilayer from free energy simulations, and compare the results to experimental data. In particular, we check the dependence of the free energy barriers ions face at the channel center on the membrane size. The results indicate that there is a serious problem with the current rigid force fields independent of the membrane size, and new, possibly polarizable, force fields need to be developed to resolve this problem.

Energetics of ion permeation, rejection, binding, and block in gramicidin A from free energy simulations.

The rigid force fields currently used in molecular dynamics (MD) simulations of biomolecules are optimized for globular proteins. Whether they can also be used in MD simulations of membrane proteins is an important issue that needs to be resolved. Here we address this issue using the gramicidin A channel, which provides an ideal test case because of the simplicity of its structure and the availability of a wealth of functional data. Permeation properties of gramicidin A can be summarized as "it conducts monovalent cations, rejects anions, and binds divalent cations." Hence, a comprehensive test should consider the energetics of permeation for all three types of ions. To that end, we construct the potential of mean force for K(+), Cl(-), and Ca(2+) ions along the channel axis. For an independent check of the potential-of-mean-force results, we also calculate the free energy differences for these ions at the channel center and binding sites relative to bulk. We find that "rejection of anions" is satisfied but there are difficulties in accommodating the other two properties using the current MD force fields.

Role of protein flexibility in ion permeation: a case study in gramicidin A.

Proteins have a flexible structure, and their atoms exhibit considerable fluctuations under normal operating conditions. However, apart from some enzyme reactions involving ligand binding, our understanding of the role of flexibility in protein function remains mostly incomplete. Here we investigate this question in the realm of membrane proteins that form ion channels. Specifically, we consider ion permeation in the gramicidin A channel, and study how the energetics of ion conduction changes as the channel structure is progressively changed from completely flexible to a fixed one. For each channel structure, the potential of mean force for a permeating potassium ion is determined from molecular dynamics (MD) simulations. Using the same molecular dynamics data for completely flexible gramicidin A, we also calculate the average densities and fluctuations of the peptide atoms and investigate the correlations between these fluctuations and the motion of a permeating ion. Our results show conclusively that peptide flexibility plays an important role in ion permeation in the gramicidin A channel, thus providing another reason--besides the well-known problem with the description of single file pore water--why this channel cannot be modeled using continuum electrostatics with a fixed structure. The new method developed here for studying the role of protein flexibility on its function clarifies the contributions of the fluctuations to energy and entropy, and places limits on the level of detail required in a coarse-grained model.

Test of molecular dynamics force fields in gramicidin A.

The force fields commonly used in molecular dynamics simulations of proteins are optimized under bulk conditions. Whether the same force fields can be used in simulations of membrane proteins is not well established, although they are increasingly being used for such purposes. Here we consider ion permeation in the gramicidin A channel as a test of the AMBER force field in a membrane environment. The potentials of mean force for potassium ions are calculated along the channel axis and compared with the one deduced from the experimental conductance data. The calculated result indicates a rather large central barrier similar to those obtained from other force fields, which are incompatible with the conductance data. We suggest that lack of polarizability is the most likely cause of this problem, and, therefore, urge development of polarizable force fields for simulations of membrane proteins.

Gramicidin A channel as a test ground for molecular dynamics force fields.

We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyristoylphosphatidylcholine bilayer, and the potential of mean force of a K(+) ion is calculated along the channel axis using the umbrella sampling method. Calculations are performed using two of the most common force fields in MD simulations: CHARMM and GROMACS. Both force fields lead to large central barriers for K(+) ion permeation, that are substantially higher than those deduced from the physiological data by inverse methods. In long MD simulations lasting over 60 ns, several ions are observed to enter the binding site but none of them crossed the channel despite the presence of a large driving field. The present results, taken together with many earlier studies, highlights the shortcomings of the standard force fields used in MD simulations of ion channels and calls for construction of more appropriate force fields for this purpose.

Role of the dielectric constants of membrane proteins and channel water in ion permeation.

Using both analytical solutions obtained from simplified systems and numerical results from more realistic cases, we investigate the role played by the dielectric constant of membrane proteins epsilon(p) and pore water epsilon(w) in permeation of ions across channels. We show that the boundary and its curvature are the crucial factors in determining how an ion's potential energy depends on the dielectric constants near an interface. The potential energy of an ion outside a globular protein has a dominant 1/epsilon(w) dependence, but this becomes 1/epsilon(p) for an ion inside a cavity. For channels, where the boundaries are in between these two extremes, the situation is more complex. In general, we find that variations in epsilon(w) have a much larger impact on the potential energy of an ion compared to those in epsilon(p). Therefore a better understanding of the effective epsilon(w) values employed in channel models is desirable. Although the precise value of epsilon(p) is not a crucial determinant of ion permeation properties, it still needs to be chosen carefully when quantitative comparisons with data are made.

Physics of ion channels.

We review the basic physics involved in transport of ions across membrane channels in cells. Electrochemical forces that control the diffusion of ions are discussed both from microscopic and macroscopic perspectives. A case is made for use of Brownian dynamics as the minimal phenomenological model that provides a bridge between experiments and more fundamental theoretical approaches. Application of Brownian and molecular dynamics methods to channels with known molecular structures is discussed.