Identification and functional evaluation of GRIA1 missense and truncation variants in individuals with ID: An emerging neurodevelopmental syndrome.

GRIA1 encodes the GluA1 subunit of α-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) receptors, which are ligand-gated ion channels that act as excitatory receptors for the neurotransmitter L-glutamate (Glu). AMPA receptors (AMPARs) are homo- or heteromeric protein complexes with four subunits, each encoded by different genes, GRIA1 to GRIA4. Although GluA1-containing AMPARs have a crucial role in brain function, the human phenotype associated with deleterious GRIA1 sequence variants has not been established. Subjects with de novo missense and nonsense GRIA1 variants were identified through international collaboration. Detailed phenotypic and genetic assessments of the subjects were carried out and the pathogenicity of the variants was evaluated in vitro to characterize changes in AMPAR function and expression. In addition, two Xenopus gria1 CRISPR-Cas9 F0 models were established to characterize the in vivo consequences. Seven unrelated individuals with rare GRIA1 variants were identified. One individual carried a homozygous nonsense variant (p.Arg377Ter), and six had heterozygous missense variations (p.Arg345Gln, p.Ala636Thr, p.Ile627Thr, and p.Gly745Asp), of which the p.Ala636Thr variant was recurrent in three individuals. The cohort revealed subjects to have a recurrent neurodevelopmental disorder mostly affecting cognition and speech. Functional evaluation of major GluA1-containing AMPAR subtypes carrying the GRIA1 variant mutations showed that three of the four missense variants profoundly perturb receptor function. The homozygous stop-gain variant completely destroys the expression of GluA1-containing AMPARs. The Xenopus gria1 models show transient motor deficits, an intermittent seizure phenotype, and a significant impairment to working memory in mutants. These data support a developmental disorder caused by both heterozygous and homozygous variants in GRIA1 affecting AMPAR function.

A single historical substitution drives an increase in acetylcholine receptor complexity.

Human adult muscle-type acetylcholine receptors are heteropentameric ion channels formed from four different, but evolutionarily related, subunits. These subunits assemble with a precise stoichiometry and arrangement such that two chemically distinct agonist-binding sites are formed between specific subunit pairs. How this subunit complexity evolved and became entrenched is unclear. Here we show that a single historical amino acid substitution is able to constrain the subunit stoichiometry of functional acetylcholine receptors. Using a combination of ancestral sequence reconstruction, single-channel electrophysiology, and concatenated subunits, we reveal that an ancestral ß-subunit can not only replace the extant ß-subunit but can also supplant the neighboring δ-subunit. By forward evolving the ancestral ß-subunit with a single amino acid substitution, we restore the requirement for a δ-subunit for functional channels. These findings reveal that a single historical substitution necessitates an increase in acetylcholine receptor complexity and, more generally, that simple stepwise mutations can drive subunit entrenchment in this model heteromeric protein.

Coupling structure with function in acid-sensing ion channels: challenges in pursuit of proton sensors.

Acid-sensing ion channels (ASICs) are a class of trimeric cation-selective ion channels activated by changes in pH within the physiological range. They are widely expressed in the central and peripheral nervous systems where they participate in a range of physiological and pathophysiological situations such as learning and memory, pain sensation, fear and anxiety, substance abuse and cell death. ASICs are localized to cell bodies and dendrites, including the postsynaptic density, and within the last 5 years several examples of proton-evoked ASIC excitatory postsynaptic currents have emerged. Thus, ASICs have become bona fide neurotransmitter-gated ion channels, activated by the smallest neurotransmitter possible: protons. Here we review how protons are thought to drive the conformational changes associated with ASIC activation and desensitization. In particular, we weigh the evidence for and against the so-called 'acidic pocket' being a vital proton sensor and discuss the emerging role of the ß11-12 linker as a desensitization switch or 'molecular clutch'. We also examine how proton-induced conformational changes pose unique challenges to classical molecular dynamics simulations, as well as some possible solutions. Given the emergence of new methodologies and structures, the coming years will probably see many advances in the study of acid-sensing ion channels.

The V2475F CPVT1 mutation yields distinct RyR2 channel populations that differ in their responses to cytosolic Ca2+ and Mg2.

Catecholaminergic polymorphic ventricular tachycardia type 1 (CPVT1) is a lethal genetic disease causing arrhythmias and sudden cardiac death in children and young adults and is linked to mutations in the cardiac ryanodine receptor (RyR2). The effects of CPVT1 mutations on RyR2 ion-channel function are often investigated using purified recombinant RyR2 channels homozygous for the mutation. However, CPVT1 patients are heterozygous for the disease, so this approach does not reveal the true changes to RyR2 function across the entire RyR2 population of channels in the heart. We therefore investigated the native cardiac RyR2 single-channel abnormalities in mice heterozygous for the CPVT1 mutation, V2475F(+/-)-RyR2, and applied molecular modelling techniques to investigate the possible structural changes that could initiate any altered function. We observed that increased sensitivity of cardiac V2475F(+/-)-RyR2 channels to both activating and inactivating levels of cytosolic Ca2+ , plus attenuation of Mg2+ inhibition, were the most marked changes. Severity of abnormality was not uniform across all channels, giving rise to multiple sub-populations with differing functional characteristics. For example, 46% of V2475F(+/-)-RyR2 channels exhibited reduced Mg2+ inhibition and 23% were actually activated by Mg2+ . Using homology modelling, we discovered that V2475 is situated at a hinge between two regions of the RyR2 helical domain 1 (HD1). Our model proposes that detrimental functional changes to RyR2 arise because mutation at this critical site reduces the angle between these regions. Our results demonstrate the necessity of characterising the total heterozygous population of CPVT1-mutated channels in order to understand CPVT1 phenotypes in patients. KEY POINTS: RyR2 mutations can cause type-1 catecholaminergic polymorphic ventricular tachycardia (CPVT1), a lethal, autosomal-dominant arrhythmic disease. However, the changes in RyR2 ion-channel function that result from the many different patient mutations are rarely investigated in detail and often only recombinant RyR2, homozygous for the mutation, is studied. As CPVT1 is a heterozygous disease and the tetrameric RyR2 channels expressed in the heart will contain varying numbers of mutated monomers, we have investigated the range of RyR2 single-channel abnormalities found in the hearts of mice heterozygous for the CPVT1 mutation, V2475F(+/-)-RyR2. Specific alterations to ligand regulation of V2475F(+/-)-RyR2 were observed. Multiple sub-populations of channels exhibited varying degrees of abnormality. In particular, an increased sensitivity to activating and inactivating cytosolic [Ca2+ ], and reduced sensitivity to Mg2+ inhibition were evident. Our results provide mechanistic insight into the changes to RyR2 gating that destabilise sarcoplasmic reticulum Ca2+ -release causing life-threatening arrhythmias in V2475F(+/-)-CPVT1 patients.

Identification of compounds that bind the centriolar protein SAS-6 and inhibit its oligomerization.

Centrioles are key eukaryotic organelles that are responsible for the formation of cilia and flagella, and for organizing the microtubule network and the mitotic spindle in animals. Centriole assembly requires oligomerization of the essential protein spindle assembly abnormal 6 (SAS-6), which forms a structural scaffold templating the organization of further organelle components. A dimerization interaction between SAS-6 N-terminal "head" domains was previously shown to be essential for protein oligomerization in vitro and for function in centriole assembly. Here, we developed a pharmacophore model allowing us to assemble a library of low-molecular-weight ligands predicted to bind the SAS-6 head domain and inhibit protein oligomerization. We demonstrate using NMR spectroscopy that a ligand from this family binds at the head domain dimerization site of algae, nematode, and human SAS-6 variants, but also that another ligand specifically recognizes human SAS-6. Atomistic molecular dynamics simulations starting from SAS-6 head domain crystallographic structures, including that of the human head domain which we now resolve, suggest that ligand specificity derives from favorable Van der Waals interactions with a hydrophobic cavity at the dimerization site.

Mutational Analysis and Modeling of Negative Allosteric Modulator Binding Sites in AMPA Receptors.

The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) constitute a subclass of the ionotropic glutamate receptor superfamily, which functions as glutamate-gated cation channels to mediate the majority of excitatory neurotransmission in the central nervous system. AMPARs are therapeutic targets in a range of brain disorders associated with abnormal glutamate hyperactivity. Multiple classes of AMPAR inhibitors have been developed during the past decades, including competitive antagonists, ion channel blockers, and negative allosteric modulators (NAMs). At present, the NAM is the only class of AMPAR ligands that have been developed into safe and useful drugs in humans in the form of perampanel (Fycompa), which was recently approved for treatment of epilepsy. Compared with the detailed understanding of other AMPAR ligand classes, surprisingly little information has been available regarding the molecular mechanism of perampanel and other classes of NAMs at AMPARs; including the location and structure of NAM binding pockets in the receptor complex. However, structures of the AMPAR GluA2 in complex with NAMs were recently reported that unambiguously identified the NAM binding sites. In parallel with this work, our aim with the present study was to identify specific residues involved in the formation of the NAM binding site for three prototypical AMPAR NAMs. Hence, we have performed a mutational analysis of the AMPAR region that links the four extracellular ligand-binding domains to the central ion channel in the transmembrane domain region. Furthermore, we perform computational ligand docking of the NAMs into structural models of the homomeric GluA2 receptor and optimize side chain conformations around the NAMs to model how NAMs bind in this specific site. The new insights provide potentially valuable input for structure-based drug design of new NAMs. SIGNIFICANCE STATEMENT: The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors are glutamate-gated ion channels that mediate the majority of excitatory neurotransmission in the brain. Negative allosteric modulators of AMPA receptors are considered to have significant therapeutic potential in diseases linked to glutamate hyperactivity. The present work employs mutational analysis and molecular modeling of the binding site for prototypical NAMs to provide new molecular insight into how NAMs interact with the AMPA receptor, which is of potential use for future design of new types of NAMs.

Role of the Cys Loop and Transmembrane Domain in the Allosteric Modulation of α4ß2 Nicotinic Acetylcholine Receptors.

Allosteric modulators of pentameric ligand-gated ion channels are thought to act on elements of the pathways that couple agonist binding to channel gating. Using α4ß2 nicotinic acetylcholine receptors and the α4ß2-selective positive modulators 17ß-estradiol (ßEST) and desformylflustrabromine (dFBr), we have identified pathways that link the binding sites for these modulators to the Cys loop, a region that is critical for channel gating in all pentameric ligand-gated ion channels. Previous studies have shown that the binding site for potentiating ßEST is in the C-terminal (post-M4) region of the α4 subunit. Here, using homology modeling in combination with mutagenesis and electrophysiology, we identified the binding site for potentiating dFBr on the top half of a cavity between the third (M3) and fourth transmembrane (M4) α-helices of the α4 subunit. We found that the binding sites for ßEST and dFBr communicate with the Cys loop, through interactions between the last residue of post-M4 and Phe170 of the conserved FPF sequence of the Cys loop, and that these interactions affect potentiating efficacy. In addition, interactions between a residue in M3 (Tyr309) and Phe167, a residue adjacent to the Cys loop FPF motif, also affect dFBr potentiating efficacy. Thus, the Cys loop acts as a key control element in the allosteric transduction pathway for potentiating ßEST and dFBr. Overall, we propose that positive allosteric modulators that bind the M3-M4 cavity or post-M4 region increase the efficacy of channel gating through interactions with the Cys loop.

Functional Validation of Heteromeric Kainate Receptor Models.

Kainate receptors require the presence of external ions for gating. Most work thus far has been performed on homomeric GluK2 but, in vivo, kainate receptors are likely heterotetramers. Agonists bind to the ligand-binding domain (LBD) which is arranged as a dimer of dimers as exemplified in homomeric structures, but no high-resolution structure currently exists of heteromeric kainate receptors. In a full-length heterotetramer, the LBDs could potentially be arranged either as a GluK2 homomer alongside a GluK5 homomer or as two GluK2/K5 heterodimers. We have constructed models of the LBD dimers based on the GluK2 LBD crystal structures and investigated their stability with molecular dynamics simulations. We have then used the models to make predictions about the functional behavior of the full-length GluK2/K5 receptor, which we confirmed via electrophysiological recordings. A key prediction and observation is that lithium ions bind to the dimer interface of GluK2/K5 heteromers and slow their desensitization.

Kainate receptor pore-forming and auxiliary subunits regulate channel block by a novel mechanism.

KEY POINTS: Kainate receptor heteromerization and auxiliary subunits, Neto1 and Neto2, attenuate polyamine ion-channel block by facilitating blocker permeation. Relief of polyamine block in GluK2/GluK5 heteromers results from a key proline residue that produces architectural changes in the channel pore α-helical region. Auxiliary subunits exert an additive effect to heteromerization, and thus relief of polyamine block is due to a different mechanism. Our findings have broad implications for work on polyamine block of other cation-selective ion channels. ABSTRACT: Channel block and permeation by cytoplasmic polyamines is a common feature of many cation-selective ion channels. Although the channel block mechanism has been studied extensively, polyamine permeation has been considered less significant as it occurs at extreme positive membrane potentials. Here, we show that kainate receptor (KAR) heteromerization and association with auxiliary proteins, Neto1 and Neto2, attenuate polyamine block by enhancing blocker permeation. Consequently, polyamine permeation and unblock occur at more negative and physiologically relevant membrane potentials. In GluK2/GluK5 heteromers, enhanced permeation is due to a single proline residue in GluK5 that alters the dynamics of the α-helical region of the selectivity filter. The effect of auxiliary proteins is additive, and therefore the structural basis of polyamine permeation and unblock is through a different mechanism. As native receptors are thought to assemble as heteromers in complex with auxiliary proteins, our data identify an unappreciated impact of polyamine permeation in shaping the signalling properties of neuronal KARs and point to a structural mechanism that may be shared amongst other cation-selective ion channels.

Steered Molecular Dynamics Simulations Predict Conformational Stability of Glutamate Receptors.

The stability of protein-protein interfaces can be essential for protein function. For ionotropic glutamate receptors, a family of ligand-gated ion channels vital for normal function of the central nervous system, such an interface exists between the extracellular ligand binding domains (LBDs). In the full-length protein, the LBDs are arranged as a dimer of dimers. Agonist binding to the LBDs opens the ion channel, and briefly after activation the receptor desensitizes. Several residues at the LBD dimer interface are known to modulate desensitization, and conformational changes around these residues are believed to be involved in the state transition. The general hypothesis is that the interface is disrupted upon desensitization, and structural evidence suggests that the disruption might be substantial. However, when cross-linking the central part of this interface, functional data suggest that the receptor can still undergo desensitization, contradicting the hypothesis of major interface disruption. Here, we illustrate how opening the dimer interface using steered molecular dynamics (SMD) simulations, and analyzing the work values required, provides a quantitative measure for interface stability. For one subtype of glutamate receptors, which is regulated by ion binding to the dimer interface, we show that opening the interface without ions bound requires less work than with ions present, suggesting that ion binding indeed stabilizes the interface. Likewise, for interface mutants with longer-lived active states, the interface is more stable, while the work required to open the interface is reduced for less active mutants. Moreover, a cross-linked mutant can still undergo initial interface opening motions similar to the native receptor and at similar energetic cost. Thus, our results support that interface opening is involved in desensitization. Furthermore, they provide reconciliation of apparently opposing data and demonstrate that SMD simulations can give relevant biological insight into longer time scale processes without the need for expensive calculations.

Ion pathways in the sarcoplasmic reticulum Ca2+-ATPase.

The sarco/endoplasmic reticulum Ca(2+)-ATPase (SERCA) is a transmembrane ion transporter belonging to the P(II)-type ATPase family. It performs the vital task of re-sequestering cytoplasmic Ca(2+) to the sarco/endoplasmic reticulum store, thereby also terminating Ca(2+)-induced signaling such as in muscle contraction. This minireview focuses on the transport pathways of Ca(2+) and H(+) ions across the lipid bilayer through SERCA. The ion-binding sites of SERCA are accessible from either the cytoplasm or the sarco/endoplasmic reticulum lumen, and the Ca(2+) entry and exit channels are both formed mainly by rearrangements of four N-terminal transmembrane α-helices. Recent improvements in the resolution of the crystal structures of rabbit SERCA1a have revealed a hydrated pathway in the C-terminal transmembrane region leading from the ion-binding sites to the cytosol. A comparison of different SERCA conformations reveals that this C-terminal pathway is exclusive to Ca(2+)-free E2 states, suggesting that it may play a functional role in proton release from the ion-binding sites. This is in agreement with molecular dynamics simulations and mutational studies and is in striking analogy to a similar pathway recently described for the related sodium pump. We therefore suggest a model for the ion exchange mechanism in P(II)-ATPases including not one, but two cytoplasmic pathways working in concert.

Small-molecule positive allosteric modulation of homomeric kainate receptors GluK1-3: development of screening assays and insight into GluK3 structure.

The kainate receptors GluK1-3 (glutamate receptor ionotropic, kainate receptors 1-3) belong to the family of ionotropic glutamate receptors and are essential for fast excitatory neurotransmission in the brain, and are associated with neurological and psychiatric diseases. How these receptors can be modulated by small-molecule agents is not well understood, especially for GluK3. We show that the positive allosteric modulator BPAM344 can be used to establish robust calcium-sensitive fluorescence-based assays to test agonists, antagonists, and positive allosteric modulators of GluK1-3. The half-maximal effective concentration (EC50) of BPAM344 for potentiating the response of 100 µm kainate was determined to be 26.3 µm for GluK1, 75.4 µm for GluK2, and 639 µm for GluK3. Domoate was found to be a potent agonist for GluK1 and GluK2, with an EC50 of 0.77 and 1.33 µm, respectively, upon co-application of 150 µm BPAM344. At GluK3, domoate acts as a very weak agonist or antagonist with a half-maximal inhibitory concentration (IC50) of 14.5 µm, in presence of 500 µm BPAM344 and 100 µm kainate for competition binding. Using H523A-mutated GluK3, we determined the first dimeric structure of the ligand-binding domain by X-ray crystallography, allowing location of BPAM344, as well as zinc-, sodium-, and chloride-ion binding sites at the dimer interface. Molecular dynamics simulations support the stability of the ion sites as well as the involvement of Asp761, Asp790, and Glu797 in the binding of zinc ions. Using electron microscopy, we show that, in presence of glutamate and BPAM344, full-length GluK3 adopts a dimer-of-dimers arrangement.

Multiscale molecular dynamics simulations predict arachidonic acid binding sites in human ASIC1a and ASIC3 transmembrane domains.

Acid-sensing ion channels (ASICs) play important roles in inflammatory pathways by conducting ions across the neuronal membrane in response to proton binding under acidic conditions. Recent studies have shown that ASICs can be modulated by arachidonic acid (AA), and, in the case of the ASIC3 subtype, even activated by AA at physiological pH. However, the mechanism by which these fatty acids act on the channel is still unknown. Here, we have used multiscale molecular dynamics simulations to predict a putative, general binding region of AA to models of the human ASIC protein. We have identified, in agreement with recent studies, residues in the outer leaflet transmembrane region which interact with AA. In addition, despite their similar modulation, we observe subtle differences in the AA interaction pattern between human ASIC1a and human ASIC3, which can be reversed by mutating three key residues at the outer leaflet portion of TM1. We further probed interactions with these residues in hASIC3 using atomistic simulations and identified possible AA coordinating interactions; salt bridge interactions of AA with R65hASIC3 and R68hASIC3 and AA tail interactions with the Y58hASIC3 aromatic ring. We have shown that longer fatty acid tails with more double bonds have increased relative occupancy in this region of the channel, a finding supported by recent functional studies. We further proposed that the modulatory effect of AA on ASIC does not result from changes in local membrane curvature. Rather, we speculate that it may occur through structural changes to the ion channel upon AA binding.

Tracing cytoplasmic Ca(2+) ion and water access points in the Ca(2+)-ATPase.

Sarco(endo)plasmic reticulum Ca(2+)-ATPase (SERCA) transports two Ca(2+) ions across the membrane of the sarco(endo)plasmic reticulum against the concentration gradient, harvesting the required energy by hydrolyzing one ATP molecule during each transport cycle. Although SERCA is one of the best structurally characterized membrane transporters, it is still largely unknown how the transported Ca(2+) ions reach their transmembrane binding sites in SERCA from the cytoplasmic side. Here, we performed extended all-atom molecular dynamics simulations of SERCA. The calculated electrostatic potential of the protein reveals a putative mechanism by which cations may be attracted to and bind to the Ca(2+)-free state of the transporter. Additional molecular dynamics simulations performed on a Ca(2+)-bound state of SERCA reveal a water-filled pathway that may be used by the Ca(2+) ions to reach their buried binding sites from the cytoplasm. Finally, several residues that are involved in attracting and guiding the cations toward the possible entry channel are identified. The results point to a single Ca(2+) entry site close to the kinked part of the first transmembrane helix, in a region loaded with negatively charged residues. From this point, a water pathway outlines a putative Ca(2+) translocation pathway toward the transmembrane ion-binding sites.

Statin activation of skeletal ryanodine receptors (RyR1) is a class effect but separable from HMG-CoA reductase inhibition.

BACKGROUND AND PURPOSE: Statins, inhibitors of HMG-CoA reductase, are mainstay treatment for hypercholesterolaemia. However, muscle pain and weakness prevent many patients from benefiting from their cardioprotective effects. We previously demonstrated that simvastatin activates skeletal ryanodine receptors (RyR1), an effect that could be important in initiating myopathy. Using a range of structurally diverse statin analogues, we examined structural features associated with RyR1 activation, aiming to identify statins lacking this property. EXPERIMENTAL APPROACH: Compounds were screened for RyR1 activity utilising [3 H]ryanodine binding. Mechanistic insight into RyR1 activity was studied by incorporating RyR1 channels from sheep, mouse or rabbit skeletal muscle into bilayers. KEY RESULTS: All UK-prescribed statins activated RyR1 at nanomolar concentrations. Cerivastatin, withdrawn from the market due to life-threatening muscle-related side effects, was more effective than currently-prescribed statins and possessed the unique ability to open RyR1 channels independently of cytosolic Ca2+ . We synthesised the one essential structural moiety that all statins must possess for HMG-CoA reductase inhibition, the R-3,5-dihydroxypentanoic acid unit, and it did not activate RyR1. We also identified five analogues retaining potent HMG-CoA reductase inhibition that inhibited RyR1 and four that lacked the ability to modulate RyR1. CONCLUSION AND IMPLICATIONS: That cerivastatin activates RyR1 most strongly supports the hypothesis that RyR1 activation is implicated in statin-induced myopathy. Demonstrating that statin regulation of RyR1 and HMG-CoA reductase are separable effects will allow the role of RyR1 in statin-induced myopathy to be further elucidated by the tool compounds we have identified, allowing development of effective cardioprotective statins with improved patient tolerance.

Protonation states of important acidic residues in the central Ca²âº ion binding sites of the Ca²âº-ATPase: a molecular modeling study.

The P-type ATPases are responsible for the transport of cations across cell membranes. The sarco(endo)plasmic reticulum Ca²âº-ATPase (SERCA) transports two Ca²âº ions from the cytoplasm to the lumen of the sarco(endo)plasmic reticulum and countertransports two or three protons per catalytic cycle. Two binding sites for Ca²âº ions have been located via protein crystallography, including four acidic amino acid residues that are essential to the ion coordination. In this study, we present molecular dynamics (MD) simulations examining the protonation states of these amino acid residues in a Ca²âº-free conformation of SERCA. Such knowledge will be important for an improved understanding of atomistic details of the transport mechanism of protons and Ca²âº ions. Eight combinations of the protonation of four central acidic residues, Glu309, Glu771, Asp800, and Glu908, are tested from 10 ns MD simulations with respect to protein stability and ability to maintain a structure similar to the crystal structure. The trajectories for the most prospective combinations of protonation states were elongated to 50 ns and subjected to more detailed analysis, including prediction of pK(a) values of the four acidic residues over the trajectories. From the simulations we find that the combination leaving only Asp800 as charged is most likely. The results are compared to available experimental data and explain the observed destabilization upon full deprotonation, resulting in the entry of cytoplasmic K⁺ ions into the Ca²âº binding sites during the simulation in which Ca²âº ions are absent. Furthermore, a hypothesis for the exchange of protons from the central binding cavity is proposed.

Binding and orientation of tricyclic antidepressants within the central substrate site of the human serotonin transporter.

Tricyclic antidepressants (TCAs) have been used for decades, but their orientation within and molecular interactions with their primary target is yet unsettled. The recent finding of a TCA binding site in the extracellular vestibule of the bacterial leucine transporter 11 A above the central site has prompted debate about whether this vestibular site in the bacterial transporter is applicable to binding of antidepressants to their relevant physiological target, the human serotonin transporter (hSERT). We present an experimentally validated structural model of imipramine and analogous TCAs in the central substrate binding site of hSERT. Two possible binding modes were observed from induced fit docking calculations. We experimentally validated a single binding mode by combining mutagenesis of hSERT with uptake inhibition studies of different TCA analogs according to the paired mutation ligand analog complementation paradigm. Using this experimental method, we identify a salt bridge between the tertiary aliphatic amine and Asp(98). Furthermore, the 7-position of the imipramine ring is found vicinal to Phe(335), and the pocket lined by Ala(173) and Thr(439) is utilized by 3-substituents. These protein-ligand contact points unambiguously orient the TCA within the central binding site and reveal differences between substrate binding and inhibitor binding, giving important clues to the inhibition mechanism. Consonant with the well established competitive inhibition of uptake by TCAs, the resulting binding site for TCAs in hSERT is fully overlapping with the serotonin binding site in hSERT and dissimilar to the low affinity noncompetitive TCA site reported in the leucine transporter (LeuT).

Molecular Investigation of Chicken Acid-Sensing Ion Channel 1 ß11-12 Linker Isomerization and Channel Kinetics.

Structures of the trimeric acid-sensing ion channel have been solved in the resting, toxin-bound open and desensitized states. Within the extracellular domain, there is little difference between the toxin-bound open state and the desensitized state. The main exception is that a loop connecting the 11th and 12th ß-strand, just two amino acid residues long, undergoes a significant and functionally critical re-orientation or flipping between the open and desensitized conformations. Here we investigate how specific interactions within the surrounding area influence linker stability in the "flipped" desensitized state using all-atom molecular dynamics simulations. An inherent challenge is bringing the relatively slow channel desensitization and recovery processes (in the milliseconds to seconds) within the time window of all-atom simulations (hundreds of nanoseconds). To accelerate channel behavior, we first identified the channel mutations at either the Leu414 or Asn415 position with the fastest recovery kinetics followed by molecular dynamics simulations of these mutants in a deprotonated state, accelerating recovery. By mutating one residue in the loop and examining the evolution of interactions in the neighbor, we identified a novel electrostatic interaction and validated prior important interactions. Subsequent functional analysis corroborates these findings, shedding light on the molecular factors controlling proton-mediated transitions between functional states of the channel. Together, these data suggest that the flipped loop in the desensitized state is stabilized by interactions from surrounding regions keeping both L414 and N415 in place. Interestingly, very few mutations in the loop allow for equivalent channel kinetics and desensitized state stability. The high degree of sequence conservation in this region therefore indicates that the stability of the ASIC desensitized state is under strong selective pressure and underlines the physiological importance of desensitization.

Mutation of a conserved glutamine residue does not abolish desensitization of acid-sensing ion channel 1.

Desensitization is a common feature of ligand-gated ion channels, although the molecular cause varies widely between channel types. Mutations that greatly reduce or nearly abolish desensitization have been described for many ligand-gated ion channels, including glutamate, GABA, glycine, and nicotinic receptors, but not for acid-sensing ion channels (ASICs) until recently. Mutating Gln276 to a glycine (Q276G) in human ASIC1a was reported to mostly abolish desensitization at both the macroscopic and the single channel levels, potentially providing a valuable tool for subsequent studies. However, we find that in both human and chicken ASIC1, the effect of Q276G is modest. In chicken ASIC1, the equivalent Q277G slightly reduces desensitization when using pH 6.5 as a stimulus but desensitizes, essentially like wild-type, when using more acidic pH values. In addition, steady-state desensitization is intact, albeit right-shifted, and recovery from desensitization is accelerated. Molecular dynamics simulations indicate that the Gln277 side chain participates in a hydrogen bond network that might stabilize the desensitized conformation. Consistent with this, destabilizing this network with the Q277N or Q277L mutations largely mimics the Q277G phenotype. In human ASIC1a, the Q276G mutation also reduces desensitization, but not to the extent reported previously. Interestingly, the kinetic consequences of Q276G depend on the human variant used. In the common G212 variant, Q276G slows desensitization, while in the rare D212 variant desensitization accelerates. Our data reveal that while the Q/G mutation does not abolish or substantially impair desensitization as previously reported, it does point to unexpected differences between chicken and human ASICs and the need for careful scrutiny before using this mutation in future studies.

P-type ATPases as drug targets: tools for medicine and science.

P-type ATPases catalyze the selective active transport of ions like H+, Na+, K+, Ca2+, Zn2+, and Cu2+ across diverse biological membrane systems. Many members of the P-type ATPase protein family, such as the Na+,K+-, H+,K+-, Ca2+-, and H+-ATPases, are involved in the development of pathophysiological conditions or provide critical function to pathogens. Therefore, they seem to be promising targets for future drugs and novel antifungal agents and herbicides. Here, we review the current knowledge about P-type ATPase inhibitors and their present use as tools in science, medicine, and biotechnology. Recent structural information on a variety of P-type ATPase family members signifies that all P-type ATPases can be expected to share a similar basic structure and a similar basic machinery of ion transport. The ion transport pathway crossing the membrane lipid bilayer is constructed of two access channels leading from either side of the membrane to the ion binding sites at a central cavity. The selective opening and closure of the access channels allows vectorial access/release of ions from the binding sites. Recent structural information along with new homology modeling of diverse P-type ATPases in complex with known ligands demonstrate that the most proficient way for the development of efficient and selective drugs is to target their ion transport pathway.