Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Comprehensive strategies of machine-learning-based quantitative structure-activity relationship models

Early quantitative structure-activity relationship (QSAR) technologies have unsatisfactory versatility and accuracy in fields such as drug discovery because they are based on traditional machine learning and interpretive expert features. The development of Big Data and deep learning technologies significantly improve the processing of unstructured data and unleash the great potential of QSAR. Here we discuss the integration of wet experiments (which provide experimental data and reliable verification), molecular dynamics simulation (which provides mechanistic interpretation at the atomic/molecular levels), and machine learning (including deep learning) techniques to improve QSAR models. We first review the history of traditional QSAR and point out its problems. We then propose a better QSAR model characterized by a new iterative framework to integrate machine learning with disparate data input. Finally, we discuss the application of QSAR and machine learning to many practical research fields, including drug development and clinical trials.

Identification and analysis of structurally critical fragments in HopS2

Background

Among the diverse roles of the Type III secretion-system (T3SS), one of the notable functions is that it serves as unique nano machineries in gram-negative bacteria that facilitate the translocation of effector proteins from bacteria into their host. These effector proteins serve as potential targets to control the pathogenicity conferred to the bacteria. Despite being ideal choices to disrupt bacterial systems, it has been quite an ordeal in the recent times to experimentally reveal and establish a concrete sequence-structure-function relationship for these effector proteins. This work focuses on the disease-causing spectrum of an effector protein, HopS2 secreted by the phytopathogen Pseudomonas syringae pv. tomato DC3000.

Results

The study addresses the structural attributes of HopS2 via a bioinformatics approach to by-pass some of the experimental shortcomings resulting in mining some critical regions in the effector protein. We have elucidated the functionally important regions of HopS2 with the assistance of sequence and structural analyses. The sequence based data supports the presence of important regions in HopS2 that are present in the other functional parts of Hop family proteins. Furthermore, these regions have been validated by an ab-initio structure prediction of the protein followed by 100 ns long molecular dynamics (MD) simulation. The assessment of these secondary structural regions has revealed the stability and importance of these regions in the protein structure.

Conclusions

The analysis has provided insights on important functional regions that may be vital to the effector functioning. In dearth of ample experimental evidence, such a bioinformatics approach has helped in the revelation of a few structural regions which will aid in future experiments to attain and evaluate the structural and functional aspects of this protein family.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2551-1) contains supplementary material, which is available to authorized users.

KCNJ11, ABCC8 and TCF7L2 polymorphisms and the response to sulfonylurea treatment in patients with type 2 diabetes: a bioinformatics assessment

BACKGROUND

Type 2 diabetes (T2D) is a worldwide epidemic with considerable health and economic consequences. Sulfonylureas are widely used drugs for the treatment of patients with T2D. KCNJ11 and ABCC8 encode the Kir6.2 (pore-forming subunit) and SUR1 (regulatory subunit that binds to sulfonylurea) of pancreatic β cell KATP channel respectively with a critical role in insulin secretion and glucose homeostasis. TCF7L2 encodes a transcription factor expressed in pancreatic β cells that regulates insulin production and processing. Because mutations of these genes could affect insulin secretion stimulated by sulfonylureas, the aim of this study is to assess associations between molecular variants of KCNJ11, ABCC8 and TCF7L2 genes and response to sulfonylurea treatment and to predict their potential functional effects.

METHODS

Based on a comprehensive literature search, we found 13 pharmacogenetic studies showing that single nucleotide polymorphisms (SNPs) located in KCNJ11: rs5219 (E23K), ABCC8: rs757110 (A1369S), rs1799854 (intron 15, exon 16 -3C/T), rs1799859 (R1273R), and TCF7L2: rs7903146 (intron 4) were significantly associated with responses to sulfonylureas. For in silico bioinformatics analysis, SIFT, PolyPhen-2, PANTHER, MutPred, and SNPs3D were applied for functional predictions of 36 coding (KCNJ11: 10, ABCC8: 24, and TCF7L2: 2; all are missense), and HaploReg v4.1, RegulomeDB, and Ensembl's VEP were used to predict functions of 7 non-coding (KCNJ11: 1, ABCC8: 1, and TCF7L2: 5) SNPs, respectively.

RESULTS

Based on various in silico tools, 8 KCNJ11 missense SNPs, 23 ABCC8 missense SNPs, and 2 TCF7L2 missense SNPs could affect protein functions. Of them, previous studies showed that mutant alleles of 4 KCNJ11 missense SNPs and 5 ABCC8 missense SNPs can be successfully rescued by sulfonylurea treatments. Further, 3 TCF7L2 non-coding SNPs (rs7903146, rs11196205 and rs12255372), can change motif(s) based on HaploReg v4.1 and are predicted as risk factors by Ensembl's VEP.

CONCLUSIONS

Our study indicates that a personalized medicine approach by tailoring sulfonylurea therapy of T2D patients according to their genotypes of KCNJ11, ABCC8, and TCF7L2 could attain an optimal treatment efficacy.

Unraveling the mechanism of selective ion transport in hydrophobic subnanometer channels

Recently reported synthetic organic nanopore (SONP) can mimic a key feature of natural ion channels, i.e., selective ion transport. However, the physical mechanism underlying the K(+)/Na(+) selectivity for the SONPs is dramatically different from that of natural ion channels. To achieve a better understanding of the selective ion transport in hydrophobic subnanometer channels in general and SONPs in particular, we perform a series of ab initio molecular dynamics simulations to investigate the diffusivity of aqua Na(+) and K(+) ions in two prototype hydrophobic nanochannels: (i) an SONP with radius of 3.2 Å, and (ii) single-walled carbon nanotubes (CNTs) with radii of 3-5 Å (these radii are comparable to those of the biological potassium K(+) channels). We find that the hydration shell of aqua Na(+) ion is smaller than that of aqua K(+) ion but notably more structured and less yielding. The aqua ions do not lower the diffusivity of water molecules in CNTs, but in SONP the diffusivity of aqua ions (Na(+) in particular) is strongly suppressed due to the rugged inner surface. Moreover, the aqua Na(+) ion requires higher formation energy than aqua K(+) ion in the hydrophobic nanochannels. As such, we find that the ion (K(+) vs. Na(+)) selectivity of the (8, 8) CNT is ∼20× higher than that of SONP. Hence, the (8, 8) CNT is likely the most efficient artificial K(+) channel due in part to its special interior environment in which Na(+) can be fully solvated, whereas K(+) cannot. This work provides deeper insights into the physical chemistry behind selective ion transport in nanochannels.

Pharmacological rescue of trafficking-impaired ATP-sensitive potassium channels

ATP-sensitive potassium (K_ATP) channels link cell metabolism to membrane excitability and are involved in a wide range of physiological processes including hormone secretion, control of vascular tone, and protection of cardiac and neuronal cells against ischemic injuries. In pancreatic β-cells, K_ATP channels play a key role in glucose-stimulated insulin secretion, and gain or loss of channel function results in neonatal diabetes or congenital hyperinsulinism, respectively. The β-cell K_ATP channel is formed by co-assembly of four Kir6.2 inwardly rectifying potassium channel subunits encoded by KCNJ11 and four sulfonylurea receptor 1 subunits encoded by ABCC8. Many mutations in ABCC8 or KCNJ11 cause loss of channel function, thus, congenital hyperinsulinism by hampering channel biogenesis and hence trafficking to the cell surface. The trafficking defects caused by a subset of these mutations can be corrected by sulfonylureas, K_ATP channel antagonists that have long been used to treat type 2 diabetes. More recently, carbamazepine, an anticonvulsant that is thought to target primarily voltage-gated sodium channels has been shown to correct K_ATP channel trafficking defects. This article reviews studies to date aimed at understanding the mechanisms by which mutations impair channel biogenesis and trafficking and the mechanisms by which pharmacological ligands overcome channel trafficking defects. Insight into channel structure-function relationships and therapeutic implications from these studies are discussed.

Homology modeling a fast tool for drug discovery: current perspectives

Major goal of structural biology involve formation of protein-ligand complexes; in which the protein molecules act energetically in the course of binding. Therefore, perceptive of protein-ligand interaction will be very important for structure based drug design. Lack of knowledge of 3D structures has hindered efforts to understand the binding specificities of ligands with protein. With increasing in modeling software and the growing number of known protein structures, homology modeling is rapidly becoming the method of choice for obtaining 3D coordinates of proteins. Homology modeling is a representation of the similarity of environmental residues at topologically corresponding positions in the reference proteins. In the absence of experimental data, model building on the basis of a known 3D structure of a homologous protein is at present the only reliable method to obtain the structural information. Knowledge of the 3D structures of proteins provides invaluable insights into the molecular basis of their functions. The recent advances in homology modeling, particularly in detecting and aligning sequences with template structures, distant homologues, modeling of loops and side chains as well as detecting errors in a model contributed to consistent prediction of protein structure, which was not possible even several years ago. This review focused on the features and a role of homology modeling in predicting protein structure and described current developments in this field with victorious applications at the different stages of the drug design and discovery.

Conductance properties of the inwardly rectifying channel, Kir3.2: molecular and Brownian dynamics study

Using the recently unveiled crystal structure, and molecular and Brownian dynamics simulations, we elucidate several conductance properties of the inwardly rectifying potassium channel, Kir3.2, which is implicated in cardiac and neurological disorders. We show that the pore is closed by a hydrophobic gating mechanism similar to that observed in Kv1.2. Once open, potassium ions move into, but not out of, the cell. The asymmetrical current-voltage relationship arises from the lack of negatively charged residues at the narrow intracellular mouth of the channel. When four phenylalanine residues guarding the intracellular gate are mutated to glutamate residues, the channel no longer shows inward rectification. Inward rectification is restored in the mutant Kir3.2 when it becomes blocked by intracellular Mg(2+). Tertiapin, a polypeptide toxin isolated from the honey bee, is known to block several subtypes of the inwardly rectifying channels with differing affinities. We identify critical residues in the toxin and Kir3.2 for the formation of the stable complex. A lysine residue of tertiapin protrudes into the selectivity filter of Kir3.2, while two other basic residues of the toxin form hydrogen bonds with acidic residues located just outside the channel entrance. The depth of the potential of mean force encountered by tertiapin is -16.1kT, thus indicating that the channel will be half-blocked by 0.4μM of the toxin.

Simulations of substrate transport in the multidrug transporter EmrD

EmrD is a multidrug resistance (MDR) transporter from Escherichia coli, which is involved in the efflux of amphipathic compounds from the cytoplasm, and the first MDR member of the major facilitator superfamily to be crystallized. Molecular dynamics simulation of EmrD in a phospholipid bilayer was used to characterize the conformational dynamics of the protein. Motions that support a previously proposed lateral diffusion pathway for substrate from the cytoplasmic membrane leaflet into the EmrD central cavity were observed. In addition, the translocation pathway of meta-chloro carbonylcyanide phenylhydrazone (CCCP) was probed using both standard and steered molecular dynamics simulation. In particular, interactions of a few specific residues with CCCP have been identified. Finally, a large motion of two residues, Val 45 and Leu 233, was observed with the passage of CCCP into the periplasmic space, placing a lower bound on the extent of opening required at this end of the protein for substrate transport. Overall, our simulations probe details of the transport pathway, motions of EmrD at an atomic level of detail, and offer new insights into the functioning of MDR transporters.

Characterization of the Na⁺/H⁺ antiporter from Yersinia pestis

Yersinia pestis, the bacterium that historically accounts for the Black Death epidemics, has nowadays gained new attention as a possible biological warfare agent. In this study, its Na/H antiporter is investigated for the first time, by a combination of experimental and computational methodologies. We determined the protein's substrate specificity and pH dependence by fluorescence measurements in everted membrane vesicles. Subsequently, we constructed a model of the protein's structure and validated the model using molecular dynamics simulations. Taken together, better understanding of the Yersinia pestis Na/H antiporter's structure-function relationship may assist in studies on ion transport, mechanism of action and designing specific blockers of Na/H antiporter to help in fighting Yersinia pestis -associated infections. We hope that our model will prove useful both from mechanistic and pharmaceutical perspectives.

Diverse gating in K+ channels: differential role of the pore-helix glutamate in stabilizing the channel pore

The selectivity filter and adjacent regions in the bacterial KcsA and inwardly rectifying K(+) (Kir) channels reveal significant conformational changes that cause the channel pore to transition from an activated to inactive state (C-type inactivation) once the channel is open. The meshwork of residues stabilizing the pore of KcsA involves Glu71-Asp80 carboxyl-carboxylate interaction 'behind' the selectivity filter. Interestingly, the Kir channels do not have this exact interaction, but instead have a Glu-Arg salt bridge where the Glu is in the same position but the Arg is one position N-terminal compared to the Asp in KcsA. Also, the Kir channels lack the Trp that hydrogen bonds to Asp80 in KcsA. Here, the sequence and structural information are combined to understand the dissimilarity in the role of the pore-helix Glu in stabilizing the pore structure in KcsA and Kir channels. This review illustrates that although Glu is quite conserved among both types of channels, the network of interactions is not translatable from one channel to the other; thereby suggesting a unique phenomenon of diverse gating patterns in K(+) channels.

Comparative characterization, expression pattern and function analysis of the 12-oxo-phytodienoic acid reductase gene family in rice

The 12-oxo-phytodienoic acid reductases (OPRs) belong to the old yellow enzyme family of flavoenzymes and form multiple subfamilies in angiosperm plants. In our previous study, a comparative genomic analysis showed that five OPR subfamilies (subs. I-V) occur in monocots, and two subfamilies (subs. I and II) in dicots. Here, a comparative study of five OsOPR genes, representing five subfamilies (I-V) in rice, was performed to provide insights into OPR biochemical properties and physiological importance. Comparative analysis of the three-dimensional structure by homology modeling indicated all five OsOPR proteins contained a highly conserved backbone with (α/β)(8)-barrels, while two middle variable regions (MVR i and ii) were also detected and defined. Analysis of enzymatic characteristics revealed that all five OsOPR fusion proteins exhibit distinct substrate specificity. Different catalytic activity was observed using racemic OPDA and trans-2-hexen-1-al as substrates, suggesting OsOPR family genes participate in two main branches of the octadecanoid pathway, including the allene oxide synthase and hydroperoxide lyase pathways which regulate various developmental processes and/or defense responses. The transcript profiles of five OsOPR genes exhibited strong tissue-specific and inducible expression patterns under abiotic stress, hormones and plant wounding treatments. Furthermore, the transcriptions of OsOPR04-1 (OsOPR11) and OsOPR08-1 (OsOPR7), representing subs. I and II, respectively, were observed in all six selected tissues and with all above-stress treatments. This suggests that these two subfamilies play an important role during different developmental stages and in response to stresses; while the expressions of OsOPR06-1 (OsOPR6), OsOPR01-1 (OsOPR10) and OsOPR02-1 (OsOPR8), representing subs. III, IV and V respectively, were strongly up-regulated with abscisic acid (ABA) and indoleacetic acid (IAA) treatments in roots, suggesting these three subfamilies play an important role in responding to hormones especially ABA and IAA signals in roots.

Computational study of the Na+/H + antiporter from Vibrio parahaemolyticus

Sodium proton antiporters are ubiquitous membrane proteins that catalyze the exchange of Na(+) for protons throughout the biological world. The Escherichia coli NhaA is the archetypal Na(+)/H(+) antiporter and is absolutely essential for survival in high salt concentrations under alkaline conditions. Its crystal structure, accompanied by extensive molecular dynamics simulations, have provided an atomically detailed model of its mechanism. In this study, we utilized a combination of computational methodologies in order to construct a structural model for the Na(+)/H(+) antiporter from the gram-negative bacterium Vibrio parahaemolyticus. We explored its overall architecture by computational means and validated its stability and robustness. This protein belongs to a novel group of NhaA proteins that transports not only Na(+) and Li(+) as substrate ions, but K(+) as well, and was also found to miss a β-hairpin segment prevalent in other homologs of the Bacteria domain. We propose, for the first time, a structure of a prototype model of a β-hairpin-less NhaA that is selective to K(+). Better understanding of the Vibrio parahaemolyticus NhaA structure-function may assist in studies on ion transport, pH regulation and designing selective blockers.

Cooperative nature of gating transitions in K(+) channels as seen from dynamic importance sampling calculations

The growing dataset of K(+) channel x-ray structures provides an excellent opportunity to begin a detailed molecular understanding of voltage-dependent gating. These structures, while differing in sequence, represent either a stable open or closed state. However, an understanding of the molecular details of gating will require models for the transitions and experimentally testable predictions for the gating transition. To explore these ideas, we apply dynamic importance sampling to a set of homology models for the molecular conformations of K(+) channels for four different sets of sequences and eight different states. In our results, we highlight the importance of particular residues upstream from the Pro-Val-Pro (PVP) region to the gating transition. This supports growing evidence that the PVP region is important for influencing the flexibility of the S6 helix and thus the opening of the gating domain. The results further suggest how gating on the molecular level depends on intra-subunit motions to influence the cooperative behavior of all four subunits of the K(+) channel. We hypothesize that the gating process occurs in steps: first sidechain movement, then inter-S5-S6 subunit motions, and lastly the large-scale domain rearrangements.

PIK3CA somatic mutations in breast cancer: Mechanistic insights from Langevin dynamics simulations

Somatic mutations in PIK3CA (phosphatidylinositol-3 kinase, catalytic subunit, alpha isoform) are reported in breast and other human cancers to concentrate at hotspots within its kinase and helical domains. Most of these mutations cause kinase gain of function in vitro and are associated with oncogenicity in vivo. However, little is known about the mechanisms driving tumor development. We have performed computational structural studies on a homology model of wildtype PIK3CA plus recurrent H1047R, H1047L, and P539R mutations, located in the kinase and helical domains, respectively. The time evolution of the structures show that H1047R/L mutants exhibit a larger area of the catalytic cleft between the kinase N- and C-lobes compared with the wildtype that could facilitate the entrance of substrates. This larger area might yield enhanced substrate-to-product turnover associated with oncogenicity. In addition, the H1047R/L mutants display increased kinase activation loop mobility, compared with the wildtype. The P539R mutant forms more hydrogen bonds and salt-bridge interactions than the wildtype, properties that are associated with enhanced thermostability. Mutant-specific differences in the catalytic cleft and activation loop behavior suggest that structure-based mutant-specific inhibitors can be designed for PIK3CA-positive breast cancers.

Model development for the viral Kcv potassium channel

A computational model for the open state of the short viral Kcv potassium channel was created and tested based on homology modeling and extensive molecular-dynamics simulation in a membrane environment. Particular attention was paid to the structure of the highly flexible N-terminal region and to the protonation state of membrane-exposed lysine residues. Data from various experimental sources, NMR spectroscopy, and electrophysiology, as well as results from three-dimensional reference interaction site model integral equation theory were taken into account to select the most reasonable model among possible variants. The final model exhibits spontaneous ion transitions across the complete pore, with and without application of an external field. The nonequilibrium transport events could be induced reproducibly without abnormally large driving potential and without the need to place ions artificially at certain key positions along the transition path. The transport mechanism through the filter region corresponds to the classic view of single-file motion, which in our case is coupled to frequent exchange of ions between the innermost filter position and the cavity.

Molecular dynamics simulations of Na+/Cl(-)-dependent neurotransmitter transporters in a membrane-aqueous system

We have performed molecular dynamics simulations of a homology model of the human serotonin transporter (hSERT) in a membrane environment and in complex with either the natural substrate 5-HT or the selective serotonin reuptake inhibitor escitalopram. We have also included a transporter homologue, the Aquifex aeolicus leucine transporter (LeuT), in our study to evaluate the applicability of a simple and computationally attractive membrane system. Fluctuations in LeuT extracted from simulations are in good agreement with crystallographic B factors. Furthermore, key interactions identified in the X-ray structure of LeuT are maintained throughout the simulations indicating that our simple membrane system is suitable for studying the transmembrane protein hSERT in complex with 5-HT or escitalopram. For these transporter complexes, only relatively small fluctuations are observed in the ligand-binding cleft. Specific interactions responsible for ligand recognition, are identified in the hSERT-5HT and hSERT-escitalopram complexes. Our findings are in good agreement with predictions from mutagenesis studies.

Time-resolved mechanism of extracellular gate opening and substrate binding in a glutamate transporter

Glutamate transporters, also referred to as excitatory amino acid transporters (EAATs), are membrane proteins that regulate glutamatergic signal transmission by clearing excess glutamate after its release at synapses. A structure-based understanding of their molecular mechanisms of function has been elusive until the recent determination of the x-ray structure of an archaeal transporter, Glt_Ph. Glt_Ph exists as a trimer, with each subunit containing a core region that mediates substrate translocation. In the present study a series of molecular dynamics simulations have been conducted and analyzed in light of new experimental data on substrate binding properties of EAATs. The simulations provide for the first time a full atomic description of the time-resolved events that drive the recognition and binding of substrate. The core region of each subunit exhibits an intrinsic tendency to open the helical hairpin HP2 loop, the extracellular gate, within tens of nanoseconds exposing conserved polar residues that serve as attractors for substrate binding. The NMDGT motif on the partially unwound part of the transmembrane helix TM7 and the residues Asp-390 and Asp-394 on TM8 are also distinguished by their important role in substrate binding and close interaction with mediating water molecules and/or sodium ions. The simulations reveal a Na⁺ binding site comprised in part of Leu-303 on TM7 and Asp-405 on TM8 and support a role for sodium ions in stabilizing substrate-bound conformers. The functional importance of Leu-303 or its counterpart Leu-391 in human EAAT1 (hEAAT1) is confirmed by site-directed mutagenesis and Na⁺ dependence assays conducted with hEAAT1 mutants L391C and L391A.

Molecular dynamics simulations of asymmetric NaCl and KCl solutions separated by phosphatidylcholine bilayers: potential drops and structural changes induced by strong Na+-lipid interactions and finite size effects

Differences of ionic concentrations across lipid bilayers are some of the primary energetic driving forces for cellular electrophysiology. While macroscopic models of asymmetric ionic solutions are well-developed, their connection to ion, water, and lipid interactions at the atomic scale are much more poorly understood. In this study, we used molecular dynamics to examine a system of two chambers of equal ionic strength, but differing amounts of NaCl and KCl, separated by a lipid bilayer. Our expectation was that the net electrostatic potential difference between the two chambers should be small or zero. Contrary to our expectation, a large potential difference (-70 mV) slowly evolved across the two water chambers over the course of our 172-ns simulation. This potential primarily originated from strong Na(+) binding to the carbonyls of the phosphatidylcholine lipids. This ion adsorption also led to significant structural and mechanical changes in the lipid bilayer. We discuss this surprising result in the context of indirect experimental evidence for Na(+) interaction with bilayers as well as potential caveats in current biomembrane simulation methodology, including force-field parameters and finite size effects.

Molecular dynamics simulations and membrane protein structure quality

Despite a growing repertoire of membrane protein structures (currently approximately 120 unique structures), considerations of low resolution and crystallization in the absence of a lipid bilayer require the development of techniques to assess the global quality of membrane protein folds. This is also the case for assessment of, e.g. homology models of human membrane proteins based on structures of (distant) bacterial homologues. Molecular dynamics (MD) simulations may be used to help evaluate the quality of a membrane protein structure or model. We have used a structure of the bacterial ABC transporter MsbA which has the correct transmembrane helices but an incorrect handedness and topology of their packing to test simulation methods of quality assessment. An MD simulation of the MsbA model in a lipid bilayer is compared to a simulation of another bacterial ABC transporter, BtuCD. The latter structure has demonstrated good conformational stability in the same bilayer environment and over the same timescale (20 ns) as for the MsbA model simulation. A number of comparative analyses of the two simulations were performed to assess changes in the structural integrity of each protein. The results show a significant difference between the two simulations, chiefly due to the dramatic structural deformations of MsbA. We therefore propose that MD could become a useful quality control tool for membrane protein structural biology. In particular, it provides a way in which to explore the global conformational stability of a model membrane protein fold.

Identification of the PIP2-binding site on Kir6.2 by molecular modelling and functional analysis

ATP-sensitive potassium (K(ATP)) channels couple cell metabolism to electrical activity by regulating K(+) fluxes across the plasma membrane. Channel closure is facilitated by ATP, which binds to the pore-forming subunit (Kir6.2). Conversely, channel opening is potentiated by phosphoinositol bisphosphate (PIP(2)), which binds to Kir6.2 and reduces channel inhibition by ATP. Here, we use homology modelling and ligand docking to identify the PIP(2)-binding site on Kir6.2. The model is consistent with a large amount of functional data and was further tested by mutagenesis. The fatty acyl tails of PIP(2) lie within the membrane and the head group extends downwards to interact with residues in the N terminus (K39, N41, R54), transmembrane domains (K67) and C terminus (R176, R177, E179, R301) of Kir6.2. Our model suggests how PIP(2) increases channel opening and decreases ATP binding and channel inhibition. It is likely to be applicable to the PIP(2)-binding site of other Kir channels, as the residues identified are conserved and influence PIP(2) sensitivity in other Kir channel family members.

Molecular Dynamics Simulation of the Cytosolic Mouth in Kcv-Type Potassium Channels

The functional effect of mutations near the intracellular mouth of the short viral Kcv potassium channel was studied by molecular dynamics simulations. As a model system we used the analogously mutated and truncated KirBac1.1, a channel with known crystal structure that shares genuine local sequence motifs with Kcv. By a novel simulated annealing methodology for structural averaging, information about the structure and dynamics of the intracellular mouth was extracted and complemented by Poisson-Boltzmann and 3D-RISM (reference interaction site model) integral equation theory for the determination of the K+ free energy surface. Besides the wild-type analogue of Kcv with its experimental reference activity (truncated KirBac1.1), two variants were studied: a deletion mutant where the N-terminus is further truncated by eight amino acids, showing inactivity in the Kcv reference system, and a point mutant where the kink-forming proline at position 13 is substituted by alanine, resulting in hyperactivity. The computations reveal that the change of activity is closely related to a hydrophilic intracellular constriction formed by the C-terminal residues of the monomers. Hyperactivity of the point mutant is correlated with both sterical and electrostatic factors, while inactivity of the deletion mutant is related to a loss of specific salt bridge patterns between the C- and N-terminus at the constriction and to the consequences for ion passage barriers, as revealed by integral equation theory. The cytosolic gate, however, is probably formed by the N-terminal segment up to the proline kink and not by the constriction. The results are compared with design principles found for other channels.

Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study

Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.

Ion hydration in nanopores and the molecular basis of selectivity

Using a simple model, it is shown that the cost of constraining a hydrated potassium ion inside a narrow pore is smaller than the cost of constraining hydrated sodium or lithium ions in pores of radius around 1.5 A. The opposite is true for pores of radius around 2.5 A. The reason for the selectivity in the first region is that the potassium ion allows for a greater distortion of its hydration shell and can therefore maintain a better coordination, and the reason for the reverse selectivity in the second region is that the smaller ions retain their hydration shells in these pores. This is relevant to the molecular basis of ion selective channels, and since this mechanism does not depend on the molecular details of the pore, it could also operate in all sorts of nanotubes.

Gating of the ATP-sensitive K+ channel by a pore-lining phenylalanine residue

ATP-sensitive K(+) (K(ATP)) channels are gated by intracellular ATP, proton and phospholipids. The pore-forming Kir6.2 subunit has all essential machineries for channel gating by these ligands. It is known that channel gating involves the inner helix bundle of crossing in which a phenylalanine residue (Phe168) is found in the TM2 at the narrowest region of the ion-conduction pathway in the Kir6.2. Here we present evidence that Phe168-Kir6.2 functions as an ATP- and proton-activated gate via steric hindrance and hydrophobic interactions. Site-specific mutations of Phe168 to a small amino acid resulted in losses of the ATP- and proton-dependent gating, whereas the channel gating was well maintained after mutation to a bulky tryptophan, supporting the steric hindrance effect. The steric hindrance effect, though necessary, was insufficient for the gating, as mutating Phe168 to a bulky hydrophilic residue severely compromised the channel gating. Single-channel kinetics of the F168W mutant resembled the wild-type channel. Small residues increased P(open), and displayed long-lasting closures and long-lasting openings. Kinetic modeling showed that these resulted from stabilization of the channel to open and long-lived closed states, suggesting that a bulky and hydrophobic residue may lower the energy barrier for the switch between channel openings and closures. Thus, it is likely that the Phe168 acts as not only a steric hindrance gate but also potentially a facilitator of gating transitions in the Kir6.2 channel.

Modeling and simulations of a bacterial outer membrane protein: OprF from Pseudomonas aeruginosa

OprF is a major outer membrane protein from Pseudomonas aeruginosa, a homolog of OmpA from Escherichia coli. The N-terminal domains of both proteins have been demonstrated to form low conductance channels in lipid bilayers. Homology models, consisting of an eight-stranded beta-barrel, of the N-terminal domain OprF have been constructed based on the crystal structure of the corresponding domain from E. coli OmpA. OprF homology models have been evaluated via a set (6 x 10 ns) of simulations of the beta-barrel embedded within a solvated dimyristoyl-phosphatidylcholine (DMPC) bilayer. The conformational stability of the models is similar to that of the crystal structure of OmpA in comparable simulations. There is a degree of water penetration into the pore-like center of the OprF barrel. The presence of an acidic/basic (E8/K121) side-chain interaction within the OprF barrel may form a "gate" able to close/open a central pore. Lipid-protein interactions within the simulations were analyzed and revealed that aromatic side-chains (Trp, Tyr) of OprF interact with lipid headgroups. Overall, the behavior of the OprF model in simulations supports the suggestion that this molecule is comparable to OmpA. The simulations help to explain the mechanism of formation of low conductance pores within the outer membrane.

Understanding the energetics of helical peptide orientation in membranes

Understanding the energetic factors determining the positioning and orientation of single-helical peptides in membranes is of fundamental interest in structural biology. Here, a simple 5-slab continuum dielectric model for the membrane is examined that distinguishes between the solvent, headgroup, and core regions. An analytical solution for the electrostatic solvation of a single dipole and an all-atom model of N-methylacetamide are used to demonstrate the effect of the dielectric boundaries in the system on peptide dipole orientation. The dipole orientation energy is shown to dominate the electrostatic solvation energy of a polyalanine helix in the membrane. With an additional surface-area-dependent term to account for the cavity formation in the aqueous region, the continuum electrostatics description is used to examine several helical peptides, the atoms of which are explicitly represented with a molecular mechanics force field. The experimentally determined tilt angles of a number of peptides of alternating alanine and leucine residues, and of glycophorin and melittin, are accurately reproduced by the model. The factors determining the tilt angles and their fluctuations are analyzed. The tilt angles of the simpler peptides are found to increase approximately linearly with peptide length; this effect is also rationalized. The analysis and model presented here provide a step toward the prediction of helical membrane protein structure.

Common mechanism of pore opening shared by five different potassium channels

A fundamental question associated with the function of ion channels is the conformational changes that allow for reversibly opening/occluding the pore through which the cations permeate. The recently elucidated crystal structures of potassium channels reveal similar structural motifs at their pore-forming regions, suggesting that they share a common gating mechanism. The validity of this hypothesis is explored by analyzing the collective dynamics of five known K(+) channel structures. Normal-mode analysis using the Gaussian network model strikingly reveals that all five structures display the same intrinsic motions at their pore-forming region despite the differences in their sequences, structures, and activation mechanisms. Superposition of the most cooperative mode profiles shows that the identified common mechanism is a global corkscrew-like counterrotation of the extracellular and cytoplasmic (CP) regions, leading to the opening of the CP end of the pore. A second cooperative mode shared by all five K(+) channels is the extension of the extracellular and/or CP ends via alternating anticorrelated fluctuations of pairs of diagonally opposite monomers. Residues acting as hinges/anchors in both modes are highly conserved across the members of the family of K(+) channel proteins, consistent with their presently disclosed critical mechanical role in pore gating.

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Molecular modeling and molecular dynamics simulation of the human A2B adenosine receptor. The study of the possible binding modes of the A2B receptor antagonists

A molecular model of the human A(2B) adenosine receptor containing seven transmembrane alpha helices connected by three intracellular and three extracellular hydrophilic loops had been constructed. A molecular docking of seven structurally diverse xanthine antagonists of the A(2B) receptor was performed, and the differences in their binding modes were investigated. The 1 ns molecular dynamics (MD) simulations of several obtained ligand-receptor complexes inserted into the phospholipid bilayer were carried out. The conformational changes of the A(2B) receptor occurring during MD simulations were explored, and the stable binding modes of the studied antagonists were determined. According to the models presented in this work, the involvement of the His251, Asn282, Ser92, Thr89, and some aromatic residues in ligand recognition was determined. The obtained binding modes of the A(2B) antagonists demonstrate good agreement with the site-directed mutagenesis data.

Allosteric effects of external K+ ions mediated by the aspartate of the GYGD signature sequence in the Kv2.1 K+ channel

K+ channels achieve exquisite ion selectivity without jeopardizing efficient permeation by employing multiple, interacting K+-binding sites. Introduction ofa cadmium (Cd2+)-binding site in the external vestibule of Kv2.1 (drk1), allowed us to functionally characterize a binding site for external monovalent cations. Permeant ions displayed higher affinity for this site than non-permeant monovalent cations, although the selectivity profile was different from that of the channel. Point mutations identified the highly conserved aspartate residue immediately following the selectivity filter as a critical determinant of the antagonism between external K+ and Cd2+ ions. A conservative mutation at this position (D378E) significantly affected the open-state stability. Moreover, the mean open time was found to be modulated by external K+ concentration, suggesting a coupling between channel closing and the permeation process. Reducing the Rb+ conductance by mutating the selectivity filter to the sequence found inKv4.1, also significantly reduced the effectiveness ofRb+ ions to antagonize Cd2+ inhibition, thereby implicating the selectivity filter as the site at which K+ions exert their antagonistic effect on Cd2+ block. The equivalent of D378 in KcsA, D80, takes part in an inter-subunit hydrogen-bond network that allows D80to functionally interact with the selectivity filter. The results suggest that external K+ ions antagonize Cd2+inhibition (in I379C) and modulate the mean open time(in the wild-type Kv2.1) by altering the occupancy profile of the K+-binding sites in the selectivity filter.

A strategy for integrative computational physiology

Organ function (the heart beat for example) can only be understood through knowledge of molecular and cellular processes within the constraints of structure-function relations at the tissue level. A quantitative modeling framework that can deal with these multiscale issues is described here under the banner of the International Union of Physiological Sciences Physiome Project.

Ion conduction and selectivity in K(+) channels

Potassium (K(+)) channels are tetrameric membrane-spanning proteins that provide a selective pore for the conductance of K(+) across the cell membranes. These channels are most remarkable in their ability to discriminate K(+) from Na(+) by more than a thousandfold and conduct at a throughput rate near diffusion limit. The recent progress in the structural characterization of K(+) channel provides us with a unique opportunity to understand their function at the atomic level. With their ability to go beyond static structures, molecular dynamics simulations based on atomic models can play an important role in shaping our view of how ion channels carry out their function. The purpose of this review is to summarize the most important findings from experiments and computations and to highlight a number of fundamental mechanistic questions about ion conduction and selectivity that will require further work.

Molecular dynamics simulation approaches to K channels: conformational flexibility and physiological function

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

A model of voltage gating developed using the KvAP channel crystal structure

Having inspected the crystal structure of the complete KvAP channel protein, we suspect that the voltage-sensing domain is too distorted to provide reliable information about its native tertiary structure or its interactions with the central pore-forming domain. On the other hand, a second crystal structure of the isolated voltage-sensing domain may well correspond to a native open conformation. We also observe that the paddle model of gating developed from these two structures is inconsistent with many experimental results, and suspect it to be energetically unrealistic. Here we show that the isolated voltage-sensing domain crystal structure can be docked onto the pore domain portion of the full-length KvAP crystal structure in an energetically favorable way to create a model of the open conformation. Using this as a starting point, we have developed rather conventional models of resting and transition conformations based on the helical screw mechanism for the transition from the open to the resting conformation. Our models are consistent with both theoretical considerations and experimental results.

Homology Models of the Tetramerization Domain of Six Eukaryotic Voltage-gated Potassium Channels Kv1.1-Kv1.6

The homology models of the tetramerization (T1) domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were constructed based on the crystal structure of the Shaker T1 domain. The results of amino acid sequence alignment indicate that the T1 domains of these K+ channels are highly conserved, with the similarities varying from 77% between Shaker and Kv1.6 to 93% between Kv1.2 and Kv1.3. The homology models reveal that the T1 domains of these Kv channels exhibit similar folds as those of Shaker K+ channel. These models also show that each T1 monomer consists of three distinct layers, with N-terminal layer 1 and C-terminal layer 3 facing the cytoplasm and the membrane, respectively. Layer 2 exhibits the highest structural conservation because it is located around the central hydrophobic core. For each Kv channel, four identical subunits assemble into the homotetramer architecture around a four-fold axis through the hydrogen bonds and salt bridges formed by 15 highly conserved polar residues. The narrowest opening of the pore is formed by the four conserved residues corresponding to R115 of the Shaker T1 domain. The homology models of these Kv T1 domains provide particularly attractive targets for further structure-based studies.

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Concerted gating mechanism underlying KATP channel inhibition by ATP

KATP channels assemble from four regulatory SUR1 and four pore-forming Kir6.2 subunits. At the single-channel current level, ATP-dependent gating transitions between the active burst and the inactive interburst conformations underlie inhibition of the KATP channel by intracellular ATP. Previously, we identified a slow gating mutation, T171A in the Kir6.2 subunit, which dramatically reduces rates of burst to interburst transitions in Kir6.2DeltaC26 channels without SUR1 in the absence of ATP. Here, we constructed all possible mutations at position 171 in Kir6.2DeltaC26 channels without SUR1. Only four substitutions, 171A, 171F, 171H, and 171S, gave rise to functional channels, each increasing Ki,ATP for ATP inhibition by >55-fold and slowing gating to the interburst by >35-fold. Moreover, we investigated the role of individual Kir6.2 subunits in the gating by comparing burst to interburst transition rates of channels constructed from different combinations of slow 171A and fast T171 "wild-type" subunits. The relationship between gating transition rate and number of slow subunits is exponential, which excludes independent gating models where any one subunit is sufficient for inhibition gating. Rather, our results support mechanisms where four ATP sites independently can control a single gate formed by the concerted action of all four Kir6.2 subunit inner helices of the KATP channel.