Brownian dynamics simulations of ion transport through the VDAC

It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.     

Ion permeation through the alpha-hemolysin channel: theoretical studies based on Brownian dynamics and Poisson-Nernst-Plank electrodiffusion theory

Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal alpha-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is approximately 10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30-50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions.     

Molecular dynamics simulations of the Cx26 hemichannel: evaluation of structural models with Brownian dynamics

The recently published crystal structure of the Cx26 gap junction channel provides a unique opportunity for elucidation of the structure of the conductive connexin pore and the molecular determinants of its ion permeation properties (conductance, current-voltage [I-V] relations, and charge selectivity). However, the crystal structure was incomplete, most notably lacking the coordinates of the N-terminal methionine residue, which resides within the pore, and also lacking two cytosolic domains. To allow computational studies for comparison with the known channel properties, we completed the structure. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) simulations of the completed and the published Cx26 hemichannel crystal structure indicate that the pore is too narrow to permit significant ion flux. The GCMC/BD simulations predict marked inward current rectification and almost perfect anion selectivity, both inconsistent with known channel properties. The completed structure was refined by all-atom molecular dynamics (MD) simulations (220 ns total) in an explicit solvent and POPC membrane system. These MD simulations produced an equilibrated structure with a larger minimal pore diameter, which decreased the height of the permeation barrier formed by the N terminus. GCMC/BD simulations of the MD-equilibrated structure yielded more appropriate single-channel conductance and less anion/cation selectivity. However, the simulations much more closely matched experimentally determined I-V relations when the charge effects of specific co- and posttranslational modifications of Cx26 previously identified by mass spectrometry were incorporated. We conclude that the average equilibrated structure obtained after MD simulations more closely represents the open Cx26 hemichannel structure than does the crystal structure, and that co- and posttranslational modifications of Cx26 hemichannels are likely to play an important physiological role by defining the conductance and ion selectivity of Cx26 channels. Furthermore, the simulations and data suggest that experimentally observed heterogeneity in Cx26 I-V relations can be accounted for by variation in co- and posttranslational modifications.     

The individual and mean activity coefficients of an electrolyte from the inverse GCMC simulation

The grand-canonical Monte Carlo (GCMC) technique, which is primarily used to calculate the concentration of a solute for a given activity coefficient, can be inverted and applied to calculate the activity coefficients corresponding to a given concentration of a solute. The changes needed to be introduced in the GCMC algorithm are discussed in the paper. The new method called the inverse GCMC (IGCMC) technique is applied to calculate the mean activity coefficient of 1:1, 2:2, 2:1, 3:1 salts with equal and unequal ion sizes and over a wide range of the electrolyte concentrations. The results for equal ion sizes are compared with the corresponding data obtained from the GCMC of Valleau and Cohen. It is shown that the IGCMC technique, after some extensions, can be used to compute the individual ionic activity coefficients. A comparison is made with the Sloth and Sørensen individual activity results and with the theoretical predictions for the 1:1 electrolyte with unequal ion sizes.     

Simulations of ion current in realistic models of ion channels: the KcsA potassium channel

Realistic studies of ion current in biologic channels present a major challenge for computer simulation approaches. All-atom molecular dynamics simulations involve serious time limitations that prevent their use in direct evaluation of ion current in channels with significant barriers. The alternative use of Brownian dynamics (BD) simulations can provide the current for simplified macroscopic models. However, the time needed for accurate calculations of electrostatic energies can make BD simulations of ion current expensive. The present work develops an approach that overcomes some of the above challenges and allows one to simulate ion currents in models of biologic channels. Our method provides a fast and reliable estimate of the energetics of the system by combining semimacroscopic calculations of the self-energy of each ion and an implicit treatment of the interactions between the ions, as well as the interactions between the ions and the protein-ionizable groups. This treatment involves the use of the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) model in its linear response approximation (LRA) implementation, which reduces the uncertainties about the value of the protein "dielectric constant." The resulting free energy surface is used to generate the forces for on-the-fly BD simulations of the corresponding ion currents. Our model is examined in a preliminary simulation of the ion current in the KcsA potassium channel. The complete free energy profile for a single ion transport reflects reasonable energetics and captures the effect of the protein-ionized groups. This calculated profile indicates that we are dealing with the channel in its closed state. Reducing the barrier at the gate region allows us to simulate the ion current in a reasonable computational time. Several limiting cases are examined, including those that reproduce the observed current, and the nature of the productive trajectories is considered. The ability to simulate the current in realistic models of ion channels should provide a powerful tool for studies of the biologic function of such systems, including the analysis of the effect of mutations, pH, and electric potentials.     

Efficient coupling of ligand binding to channel opening by the binding domain of a modulatory (beta) subunit of the olfactory cyclic nucleotide-gated channel.

CNG channels in vivo are heteromers of homologous alpha and beta subunits that each contain a six-transmembrane segment domain and a COOH-terminal cytoplasmic cyclic nucleotide binding domain (BD). In heterologous expression systems, heteromeric alphabeta channels activate with greater sensitivity to ligand than do homomeric alpha channels; however, ligand-gating of channels containing only beta subunit BDs has never been studied because beta subunits cannot form functional homomeric CNG channels. To characterize directly the contribution of the beta subunit BD to ligand-gating, we constructed a chimeric subunit, X-beta, whose BD sequence was that of the beta subunit CNG5 from rat, but whose sequence outside the BD was derived from alpha subunits. For comparison, we constructed another chimera, X-alpha, whose sequence outside the BD was identical to that of X-beta, but whose BD sequence was that of the alpha subunit CNG2 from catfish. When expressed in Xenopus oocytes, X-beta and X-alpha each formed functional homomeric channels activated by both cAMP and cGMP. This is the first demonstration that the beta subunit BD can couple ligand binding to activation in the absence of alpha subunit BD residues. Notably, both agonists activate X-beta more effectively than X-alpha (higher opening efficacy and lower K(1/2)). The BD is believed to comprise two functionally distinct subdomains: (1) the roll subdomain (beta-roll and flanking A- and B-helices) and (2) the C-helix subdomain. Opening efficacy was previously believed to be controlled primarily by the C-helix, but when we made additional chimeras by exchanging the subdomains between X-beta and X-alpha, we found that both subdomains contain significant determinants of efficacy and agonist selectivity. In particular, only channels containing the roll subdomain of the beta subunit had high efficacy. Thermodynamic linkage analysis shows that interaction between the two subdomains accounts for a significant portion of their contribution to activation energetics.     

Testing the applicability of Nernst-Planck theory in ion channels: comparisons with Brownian dynamics simulations

The macroscopic Nernst-Planck (NP) theory has often been used for predicting ion channel currents in recent years, but the validity of this theory at the microscopic scale has not been tested. In this study we systematically tested the ability of the NP theory to accurately predict channel currents by combining and comparing the results with those of Brownian dynamics (BD) simulations. To thoroughly test the theory in a range of situations, calculations were made in a series of simplified cylindrical channels with radii ranging from 3 to 15 Å, in a more complex ‘catenary’ channel, and in a realistic model of the mechanosensitive channel MscS. The extensive tests indicate that the NP equation is applicable in narrow ion channels provided that accurate concentrations and potentials can be input as the currents obtained from the combination of BD and NP match well with those obtained directly from BD simulations, although some discrepancies are seen when the ion concentrations are not radially uniform. This finding opens a door to utilising the results of microscopic simulations in continuum theory, something that is likely to be useful in the investigation of a range of biophysical and nano-scale applications and should stimulate further studies in this direction.     

Competitive substitution of hexammine cobalt(III) for Na+ and K+ ions in oriented DNA fibers

Competition of the trivalent cation, Co(NH3)(3+)(6), with K+ and Na+ ions in binding to DNA was studied by equilibrating oriented DNA fibers with ethanol/water solutions (65 and 52% v/v EtOH), containing different combinations and concentrations of KCl and NaCl and constant concentration (0.8 mM) of Co(NH3)(6)Cl(3). The degree of Co(NH3)(3+)(6) binding to DNA does not depend significantly on the ethanol concentration or on the kind of univalent cation (Na+ or K+). The ion exchange selectivity coefficient of monovalent-trivalent ion competition, D(1)(c3), increases with the concentration of Me+, C(o)(+), and the monotonic dependence of log D(1)(c3) vs log C(o)(+) has an inflection between 100 and 300 mM that is caused by a structural transformation of DNA from A- to B-form. The ion exchange experimental data are compared with results of grand canonical Monte Carlo (GCMC) simulations of systems of parallel and hexagonally ordered, discretely charged polyions with density and spatial distribution of the charged groups modeling B- and A-forms of DNA. The GCMC method for discretely charged models of the DNA polyion produces a quantitative agreement with experimental data on trivalent-monovalent ion competition in dependence on DNA structural state and salt concentration. Based on this and previous studies it is concluded that the affinity of DNA for the cations decreases in the order Co(NH3)(3+)(6) > Ca2+ > Mg2+ > Na+ approximately K+ > Li+. DNA does not exhibit selectivity for Na+ or K+ in ethanol/water solutions either in the absence or in the presence of Co(NH3)(3+)(6), Ca2+, and Mg2+.     

Reliable prediction of adsorption isotherms via genetic algorithm molecular simulation

Conventional molecular simulation techniques such as grand canonical Monte Carlo (GCMC) strictly rely on purely random search inside the simulation box for predicting the adsorption isotherms. This blind search is usually extremely time demanding for providing a faithful approximation of the real isotherm and in some cases may lead to non-optimal solutions. A novel approach is presented in this article which does not use any of the classical steps of the standard GCMC method, such as displacement, insertation, and removal. The new approach is based on the well-known genetic algorithm to find the optimal configuration for adsorption of any adsorbate on a structured adsorbent under prevailing pressure and temperature. The proposed approach considers the molecular simulation problem as a global optimization challenge. A detailed flow chart of our so-called genetic algorithm molecular simulation (GAMS) method is presented, which is entirely different from traditions molecular simulation approaches. Three real case studies (for adsorption of CO₂ and H₂ over various zeolites) are borrowed from literature to clearly illustrate the superior performances of the proposed method over the standard GCMC technique. For the present method, the average absolute values of percentage errors are around 11% (RHO-H₂), 5% (CHA-CO₂), and 16% (BEA-CO₂), while they were about 70%, 15%, and 40% for the standard GCMC technique, respectively.     

Estimating the dielectric constant of the channel protein and pore

When modelling biological ion channels using Brownian dynamics (BD) or Poisson–Nernst–Planck theory, the force encountered by permeant ions is calculated by solving Poisson’s equation. Two free parameters needed to solve this equation are the dielectric constant of water in the pore and the dielectric constant of the protein forming the channel. Although these values can in theory be deduced by various methods, they do not give a reliable answer when applied to channel-like geometries that contain charged particles. To determine the appropriate values of the dielectric constants, here we solve the inverse problem. Given the structure of the MthK channel, we attempt to determine the values of the protein and pore dielectric constants that minimize the discrepancies between the experimentally-determined current–voltage curve and the curve obtained from BD simulations. Two different methods have been applied to determine these values. First, we use all possible pairs of the pore dielectric constant of water, ranging from 20 to 80 in steps of 10, and the protein dielectric constant of 2–10 in steps of 2, and compare the simulated results with the experimental values. We find that the best agreement is obtained with experiment when a protein dielectric constant of 2 and a pore water dielectric constant of 60 is used. Second, we employ a learning-based stochastic optimization algorithm to pick out the optimum combination of the two dielectric constants. From the algorithm we obtain an optimum value of 2 for the protein dielectric constant and 64 for the pore dielectric constant.     

Spermine: an "invisible" component in the crystals of B-DNA. A grand canonical Monte Carlo and molecular dynamics simulation study

The association of spermine(4+) (Spm(4+)), Mg(2+) and monovalent (M(+)) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm(4+), Mg(2+), and M(+) between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm(4+), 32 Na(+) and four Cl(-) ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na(+) and Spm(4+) compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.     

Exploring the origin of the ion selectivity of the KcsA potassium channel

The availability of structural information about biological ion channels provides an opportunity to gain a detailed understanding of the control of ion selectivity by biological systems. However, accomplishing this task by computer simulation approaches is very challenging. First, although the activation barriers for ion transport can be evaluated by microscopic simulations, it is hard to obtain accurate results by such approaches. Second, the selectivity is related to the actual ion current and not directly to the individual activation barriers. Thus, it is essential to simulate the ion currents and this cannot be accomplished at present by microscopic MD approaches. In order to address this challenge, we developed and refined an approach capable of evaluating ion current while still reflecting the realistic features of the given channel. Our method involves generation of semimacroscopic free energy surfaces for the channel/ions system and Brownian dynamics (BD) simulations of the corresponding ion current. In contrast to most alternative macroscopic models, our approach is able to reproduce the difference between the free energy surfaces of different ions and thus to address the selectivity problem. Our method is used in a study of the selectivity of the KcsA channel toward the K+ and Na+ ions. The BD simulations with the calculated free energy profiles produce an appreciable selectivity. To the best of our knowledge, this is the first time that the trend in the selectivity in the ion current is produced by a computer simulation approach. Nevertheless, the calculated selectivity is still smaller than its experimental estimate. Recognizing that the calculated profiles are not perfect, we examine how changes in these profiles can account for the observed selectivity. It is found that the origin of the selectivity is more complex than generally assumed. The observed selectivity can be reproduced by increasing the barrier at the exit and the entrance of the selectivity filter, but the necessary changes in the barrier approach the limit of the error in the PDLD/S-LRA calculations. Other options that can increase the selectivity are also considered, including the difference between the Na+...Na+ and K+...K+ interaction. However, this interesting effect does not appear to lead to a major difference in selectivity since the Na+ ions at the limit of strong interaction tend to move in a less concerted way than the K+ ions. Changes in the relative binding energies at the different binding sites are also not so effective in changing the selectivity. Finally, it is pointed out that using the calculated profiles as a starting point and forcing the model to satisfy different experimentally based constraints, should eventually provide more detailed understanding of the different complex factors involved in ion selectivity of biological channels.     

The cognitive impact of the ANK3 risk variant for bipolar disorder: initial evidence of selectivity to signal detection during sustained attention

Background

Abnormalities in cognition have been reported in patients with Bipolar Disorder (BD) and their first degree relatives, suggesting that susceptibility genes for BD may impact on cognitive processes. Recent genome-wide genetic studies have reported a strong association with BD in a single nucleotide polymorphism (SNP) (rs10994336) within ANK3, which codes for Ankyrin 3. This protein is involved in facilitating the propagation of action potentials by regulating the assembly of sodium gated ion channels. Since ANK3 influences the efficiency of transmission of neuronal impulses, allelic variation in this gene may have widespread cognitive effects. Preclinical data suggest that this may principally apply to sequential signal detection, a core process of sustained attention.

Methodology/Principal Findings

One hundred and eighty-nine individuals of white British descent were genotyped for the ANK3 rs10994336 polymorphism and received diagnostic interviews and comprehensive neurocognitive assessment of their general intellectual ability, memory, decision making, response inhibition and sustained attention. Participants comprised euthymic BD patients (n = 47), their unaffected first-degree relatives (n = 75) and healthy controls (n = 67). The risk allele T was associated with reduced sensitivity in target detection (p = 0.0004) and increased errors of commission (p = 0.0018) during sustained attention regardless of diagnosis. We found no effect of the ANK3 genotype on general intellectual ability, memory, decision making and response inhibition.

Conclusions/Significance

Our results suggest that allelic variation in ANK3 impacts cognitive processes associated with signal detection and this mechanism may relate to risk for BD. However, our results require independent replication and confirmation that ANK3 (rs10994336) is a direct functional variant.     

The normalized slope conductance as a tool for quantitative analysis of current-voltage relations

The patch-clamp method, which was awarded the Nobel Prize in 1991, is a well-established and indispensable method to study ion channels in living cells and to biophysically characterize non-voltage-gated ion channels, which comprise about 70% of all ion channels in the human genome. To investigate the biophysical properties of non-voltage-gated ion channels, whole-cell measurements with application of continuous voltage ramps are routinely conducted to obtain current-voltage (IV) relationships. However, adequate tools for detailed and quantitative analysis of IV curves are still missing. We use the example of the transient receptor potential classical (TRPC) channel family to elucidate whether the normalized slope conductance (NSC) is an appropriate tool for reliable discrimination of the IV curves of diverse TRPC channels that differ in their individual curve progression. We provide a robust calculation method for the NSC, and, by applying this method, we find that TRPC channel activators and modulators can evoke different NSC progressions independent from their expression levels, which points to distinguishable active channel states. TRPC6 mutations in patients with focal segmental glomerulosclerosis resulted in distinct NSC progressions, suggesting that the NSC is suitable for investigating structure-function relations and might help unravel the unknown pathomechanisms leading to focal segmental glomerulosclerosis. The NSC is an effective algorithm for extended biophysical characterization of non-voltage-gated ion channels.     

Grand canonical Monte Carlo molecular and thermodynamic predictions of ion effects on binding of an oligocation (L8+) to the center of DNA oligomers.

Grand canonical Monte Carlo (GCMC) simulations are reported for aqueous solutions containing excess univalent salt (activities a +/- = 1.76-12.3 mM) and one of the following species: an octacationic rod-like ligand, L8+; a B-DNA oligomer with N phosphate charges (8 < or = N < or = 100); or a complex resulting from the binding of L8+ at the center of an N-mer (24 < or = N < or = 250). Simplified models of these multiply charged species are used in the GCMC simulations to predict the fundamental coulombic contributions to the following experimentally relevant properties: 1) the axial distance over which ligand binding affects local counterion concentrations at the surface of the N-mer; 2) the dependence on N of GCMC preferential interaction coefficients, gamma 32MC identical to delta C3/delta C2l a +/-, T, where C3 and C2 are, respectively, the molar concentrations of salt and the multiply charged species (ligand, N-mer or complex); and 3) the dependence on N of SaKobs identical to d in Kobs/d in a +/- = delta (magnitude of ZJ + 2 gamma 32J), where Kobs is the equilibrium concentration quotient for the binding of L8+ to the center of an N-mer and delta denotes the stoichiometric combination of terms, each of which pertains to a reactant or product J having magnitude of ZJ charges. The participation of electrolyte ions in the ligand binding interaction is quantified by the magnitude of SaKobs, which reflects the net (stoichiometrically weighted) difference in the extent of thermodynamic binding of salt ions to the products and reactants. Results obtained here from GCMC simulations yield a picture of the salient molecular consequences of binding a cationic ligand, as well as thermodynamic predictions whose applicability can be tested experimentally. Formation of the central complex is predicted to cause a dramatic reduction in the surface counterion (e.g., Na+) concentration over a region including but extending well beyond the location of the ligand binding site. For binding a cationic ligand, SaKobs is predicted to be negative, indicating net electrolyte ion release in the binding process. At small enough N, -SaKobs is predicted to decrease strongly toward zero with decreasing N. At intermediate N, -SaKobs appears to exceed its limiting value as N-->infinity.     

Bathocuproine-assisted reduction of copper(II) by human albumin

Human albumin (studied here as the recombinant protein rHA), a copper-binding protein in blood plasma, is shown to reduce Cu(II) to Cu(I) in the presence of a Cu(I) chelator, bathocuproinedisulfonate (BD). This reaction was accelerated at low pH, when there was little binding of Cu(II) to rHA. The addition of a competitive metal ion, Ni(II), or an increase in the concentration of BD, enhanced the reduction of Cu(II) to Cu(I). It was concluded that the oxidant was the Cu(II) complex of BD, which is likely to bind strongly to albumin. The free thiol at Cys34 was ruled out as the sole reducing agent, since Cys34-blocked albumin also gave rise to Cu(I) in the presence of BD. Reactions with amino acids and peptides suggested that Tyr and possibly His side-chains are potential reductants. BD and its homologues are frequently used as Cu(I)-specific chelators in biological experiments, but the strong oxidant activity of [Cu(II)(BD)2]2- and its ability to bind to biological macromolecules should not be overlooked, and may artificially trigger/accelerate Cu(II) reduction.     

Distinct structural determinants of efficacy and sensitivity in the ligand-binding domain of cyclic nucleotide-gated channels

Cyclic nucleotide-gated (CNG) channels open in response to direct binding of cyclic nucleotide messengers. Every subunit in a tetrameric CNG channel contains a cytoplasmic ligand-binding domain (BD) that includes a beta-roll (flanked by short helices) and a single C-terminal helix called the C-helix that was previously found to control efficacy (maximal open probability) and selectivity for cGMP versus cAMP. We constructed a series of chimeric CNG channel subunits, each containing a distinct BD sequence (chosen from among six phylogenetically divergent isoforms) fused to an invariant non-BD sequence. We assayed these "BD substitution" chimeras as homomeric CNG channels in Xenopus oo-cytes to compare their functions and found that the most efficient activation by both cAMP and cGMP derived from the BD of the catfish CNGA4 olfactory modulatory subunit (fCNGA4). We then tested the effects of replacing subregions of the bovine CNGA1 BD with corresponding fCNGA4 sequence and hence identified parts of the fCNGA4 BD producing efficient activation. For instance, replacing either the "hinge" that connects the roll and C-helix subdomains or the BD sequence N-terminal to the hinge greatly enhanced cAMP efficacy. Replacing the "loop-beta 8" region (the C-terminal end of the beta-roll) improved agonist sensitivity for cGMP selectively over cAMP. Our results thus identify multiple BD elements outside the C-helix that control selective ligand interaction and channel gating steps by distinct mechanisms. This suggests that the purine ring of the cyclic nucleotide may interact with both the beta-roll and the C-helix at different points in the mechanism.     

Ion permeation and selectivity of OmpF porin: a theoretical study based on molecular dynamics,Brownian dynamics,and continuum electrodiffusion theory

Three different theoretical approaches are used and compared to refine our understanding of ion permeation through the channel formed by OmpF porin from Escherichia coli. Those approaches are all-atom molecular dynamics (MD) in which ions, solvent, and lipids are represented explicitly, Brownian dynamics (BD) in which ions are represented explicitly, while solvent and lipids are represented as featureless dielectrics, and Poisson-Nernst-Planck (PNP) electrodiffusion theory in which both solvent and local ion concentrations are represented as a continuum. First, the ability of the different theoretical approaches in reproducing the equilibrium average ion density distribution in OmpF porin bathed by a 1M KCl symmetric salt solution is examined. Under those conditions the PNP theory is equivalent to the non-linear Poisson-Boltzmann (PB) theory. Analysis shows that all the three approaches are able to capture the important electrostatic interactions between ions and the charge distribution of the channel that govern ion permeation and selectivity in OmpF. The K(+) and Cl(-) density distributions obtained from the three approaches are very consistent with one another, which suggests that a treatment on the basis of a rigid protein and continuum dielectric solvent is valid in the case of OmpF. Interestingly, both BD and continuum electrostatics reproduce the distinct left-handed twisted ion pathways for K(+) and Cl(-) extending over the length of the pore which were observed previously in MD. Equilibrium BD simulations in the grand canonical ensemble indicate that the channel is very attractive for cations, particularly at low salt concentration. On an average there is 1.55 K(+) inside the pore in 10mM KCl. Remarkably, there is still 0.17 K(+) on average inside the pore even at a concentration as low as 1microM KCl. Secondly, non-equilibrium ion flow through OmpF is calculated using BD and PNP and compared with experimental data. The channel conductance in 0.2M and 1M KCl calculated using BD is in excellent accord with the experimental data. The calculations reproduce the experimentally well-known conductance-concentration relation and also reveal an asymmetry in the channel conductance (a larger conductance is observed under a positive transmembrane potential). Calculations of the channel conductance for three mutants (R168A, R132A, and K16A) in 1M KCl suggest that the asymmetry in the channel conductance arises mostly from the permanent charge distribution of the channel rather than the shape of the pore itself. Lastly, the calculated reversal potential in a tenfold salt gradient (0.1:1M KCl) is 27.4(+/-1.3)mV (BD) and 22.1(+/-0.6)mV (PNP), in excellent accord with the experimental value of 24.3mV. Although most of the results from PNP are qualitatively reasonable, the calculated channel conductance is about 50% higher than that calculated from BD probably because of a lack of some dynamical ion-ion correlations.     

Predicting ion channels and their types by the dipeptide mode of pseudo amino acid composition

Ion channels are integral membrane proteins that control movement of ions into or out of cells. They are key components in a wide range of biological processes. Different types of ion channels have different biological functions. With the appearance of vast proteomic data, it is highly desirable for both basic research and drug-target discovery to develop a computational method for the reliable prediction of ion channels and their types. In this study, we developed a support vector machine-based method to predict ion channels and their types using primary sequence information. A feature selection technique, analysis of variance (ANOVA), was introduced to remove feature redundancy and find out an optimized feature set for improving predictive performance. Jackknife cross-validated results show that the proposed method can discriminate ion channels from non-ion channels with an overall accuracy of 86.6%, classify voltage-gated ion channels and ligand-gated ion channels with an overall accuracy of 92.6% and predict four types (potassium, sodium, calcium and anion) of voltage-gated ion channels with an overall accuracy of 87.8%, respectively. These results indicate that the proposed method can correctly identify ion channels and provide important instructions for drug-target discovery. The predictor can be freely downloaded from http://cobi.uestc.edu.cn/people/hlin/tools/IonchanPred/.