The gates of ion channels and enzymes

Cryo-EM structures of the plant anion channel SLAC1 from Arabidopsis thaliana suggest a combined activation model

The anion channel SLAC1 functions as a crucial effector in the ABA signaling, leading to stomata closure. SLAC1 is activated by phosphorylation in its intracellular domains. Both a binding-activation model and an inhibition-release model for activation have been proposed based on only the closed structures of SLAC1, rendering the structure-based activation mechanism controversial. Here we report cryo-EM structures of Arabidopsis SLAC1 WT and its phosphomimetic mutants in open and closed states. Comparison of the open structure with the closed ones reveals the structural basis for opening of the conductance pore. Multiple phosphorylation of an intracellular domain (ICD) causes dissociation of ICD from the transmembrane domain. A conserved, positively-charged sequence motif in the intracellular loop 2 (ICL2) seems to be capable of sensing of the negatively charged phosphorylated ICD. Interactions between ICL2 and ICD drive drastic conformational changes, thereby widening the pore. From our results we propose that SLAC1 operates by a mechanism combining the binding-activation and inhibition-release models.

Rapid Determination of the Topology of Oligomeric α-Helical Membrane Proteins by Water- and Lipid-Edited Methyl NMR

Single-span oligomeric α-helical transmembrane proteins are common in virus ion channels, which are targets of antiviral drugs. Knowledge about the high-resolution structures of these oligomeric α-helical bundles is so far scarce. Structure determination of these membrane proteins by solid-state NMR traditionally requires resolving and assigning protein chemical shifts and measuring many interhelical distances, which are time-consuming. To accelerate experimental structure determination, here we introduce a simple solid-state NMR approach that uses magnetization transfer from water and lipid protons to the protein. By detecting the water- and lipid-transferred intensities of the high-sensitivity methyl 13C signals of Leu, Val, and Ile residues, which are highly enriched in these membrane proteins, we can derive models of the topology of these homo-oligomeric helical bundles. The topology is specified by the positions of amino acid residues in heptad repeats and the orientations of residues relative to the channel pore, lipids, and the helical interface. We demonstrate this water- and lipid-edited methyl NMR approach on the envelope (E) protein of SARS-CoV-2, the causative agent of the COVID-19 pandemic. We show that water-edited and lipid-edited 2D 13C-13C correlation spectra can be measured with sufficient sensitivity. Even without resolving multiple residues of the same type in the NMR spectra, we can obtain the helical bundle topology. We apply these experiments to the structurally unknown E proteins of the MERS coronavirus and the human coronavirus NL63. The resulting structural topologies show interesting differences in the positions of the aromatic residues in these three E proteins, suggesting that these viroporins may have different mechanisms of activation and ion conduction.

Understanding Life at High Temperatures: Relationships of Molecular Channels in Enzymes of Methanogenic Archaea and Their Growth Temperatures

Analyses of protein structures have shown the existence of molecular channels in enzymes from Prokaryotes. Those molecular channels suggest a critical role of spatial voids in proteins, above all, in those enzymes functioning under high temperature. It is expected that these spaces within the protein structure are required to access the active site and to maximize availability and thermal stability of their substrates and cofactors. Interestingly, numerous substrates and cofactors have been reported to be highly temperature-sensitive biomolecules. Methanogens represent a singular phylogenetic group of Archaea that performs anaerobic respiration producing methane during growth. Methanogens inhabit a variety of environments including the full range of temperatures for the known living forms. Herein, we carry out a dimensional analysis of molecular tunnels in key enzymes of the methanogenic pathway from methanogenic Archaea growing optimally over a broad temperature range. We aim to determine whether the dimensions of the molecular tunnels are critical for those enzymes from thermophiles. Results showed that at increasing growth temperature the dimensions of molecular tunnels in the enzymes methyl-coenzyme M reductase and heterodisulfide reductase become increasingly restrictive and present strict limits at the highest growth temperatures, i.e., for hyperthermophilic methanogens. However, growth at lower temperature allows a wide dimensional range for the molecular spaces in these enzymes. This is in agreement with previous suggestions on a potential major role of molecular tunnels to maintain biomolecule stability and activity of some enzymes in microorganisms growing at high temperatures. These results contribute to better understand archaeal growth at high temperatures. Furthermore, an optimization of the dimensions of molecular tunnels would represent an important adaptation required to maintain the activity of key enzymes of the methanogenic pathway for those methanogens growing optimally at high temperatures.

Ratio of Hydrophobic-Hydrophilic and Positive-Negative Residues at Lipid-Water-Interface Influences Surface Expression and Channel Gating of TRPV1

During evolution, TRPV1 has lost, retained or selected certain residues at Lipid-Water-Interface (LWI) and formed specific patterns there. The ratio of "hydrophobic-hydrophilic" and "positive-negative-charged" residues at the inner LWI remains conserved throughout vertebrate evolution and plays important role in regulating TRPV1 trafficking and localization. Arg575 is an important residue as Arg575Asp mutant has reduced surface expression, co-localization with lipid raft markers, cell area and increased cell lethality. This lethality is most likely due to the disruption of the ratio between positive-negative charges caused by the mutation. Such lethality can be rescued by either using TRPV1-specfic inhibitor 5'-IRTX or by restoring the positive-negative charge ratio at that position, i.e. by introducing Asp576Arg mutation in Arg575Asp backbone. We propose that Arg575Asp mutation confers TRPV1 in a "constitutive-open-like" condition. These findings have broader implication in understanding the molecular evolution of thermo-sensitive ion channels and the micro-environments involved in processes that goes erratic in different diseases. The segment of TRPV1 that is present at the inner lipid-water-interface (LWI) has a specific pattern of amino acid combinations. The overall ratio of +ve charge /-ve charge and the ratio of hydrophobicity/hydrophilicity remain constant throughout the vertebrate evolution (ca 450 million years). This specific pattern is not observed in the outer LWI region of TRPV1.

High-resolution cryo-electron microscopy structure of photosystem II from the mesophilic cyanobacterium, Synechocystis sp. PCC 6803

Photosystem II (PSII) enables global-scale, light-driven water oxidation. Genetic manipulation of PSII from the mesophilic cyanobacterium Synechocystis sp. PCC 6803 has provided insights into the mechanism of water oxidation; however, the lack of a high-resolution structure of oxygen-evolving PSII from this organism has limited the interpretation of biophysical data to models based on structures of thermophilic cyanobacterial PSII. Here, we report the cryo-electron microscopy structure of PSII from Synechocystis sp. PCC 6803 at 1.93-Å resolution. A number of differences are observed relative to thermophilic PSII structures, including the following: the extrinsic subunit PsbQ is maintained, the C terminus of the D1 subunit is flexible, some waters near the active site are partially occupied, and differences in the PsbV subunit block the Large (O1) water channel. These features strongly influence the structural picture of PSII, especially as it pertains to the mechanism of water oxidation.

A biomimetic metal–organic framework with cuboid inner cavities for enantioselective separation

A biomimetic metal–organic framework with cuboid inner cavities and multiple recognition sites was constructed from a phenylalanine-derived ligand. It can enantioselectively separate various racemic alcohols, diols and epoxides with ee up to 99.5%.

Multiscale Approach for Computing Gated Ligand Binding from Molecular Dynamics and Brownian Dynamics Simulations

We develop an approach to characterize the effects of gating by a multiconformation protein consisting of macrostate conformations that are either accessible or inaccessible to ligand binding. We first construct a Markov state model of the apo-protein from atomistic molecular dynamics simulations from which we identify macrostates and their conformations, compute their relative macrostate populations and interchange kinetics, and structurally characterize them in terms of ligand accessibility. We insert the calculated first-order rate constants for conformational transitions into a multistate gating theory from which we derive a gating factor γ that quantifies the degree of conformational gating. Applied to HIV-1 protease, our approach yields a kinetic network of three accessible (semi-open, open, and wide-open) and two inaccessible (closed and a newly identified, "parted") macrostate conformations. The parted conformation sterically partitions the active site, suggesting a possible role in product release. We find that the binding kinetics of drugs and drug-like inhibitors to HIV-1 protease falls in the slow gating regime. However, because γ = 0.75, conformational gating only modestly slows ligand binding. Brownian dynamics simulations of the diffusional association of eight inhibitors to the protease─having a wide range of experimental association constants (∼10⁴-10¹⁰ M^-1 s^-1)─yields gated rate constants in the range of ∼0.5-5.7 × 10⁸ M^-1 s^-1. This indicates that, whereas the association rate of some inhibitors could be described by the model, for many inhibitors either subsequent conformational transitions or alternate binding mechanisms may be rate-limiting. For systems known to be modulated by conformational gating, the approach could be scaled computationally efficiently to screen association kinetics for a large number of ligands.

Fluorescence Approaches for Characterizing Ion Channels in Synthetic Bilayers

Ion channels are membrane proteins that play important roles in a wide range of fundamental cellular processes. Studying membrane proteins at a molecular level becomes challenging in complex cellular environments. Instead, many studies focus on the isolation and reconstitution of the membrane proteins into model lipid membranes. Such simpler, in vitro, systems offer the advantage of control over the membrane and protein composition and the lipid environment. Rhodopsin and rhodopsin-like ion channels are widely studied due to their light-interacting properties and are a natural candidate for investigation with fluorescence methods. Here we review techniques for synthesizing liposomes and for reconstituting membrane proteins into lipid bilayers. We then summarize fluorescence assays which can be used to verify the functionality of reconstituted membrane proteins in synthetic liposomes.

Tunnel engineering for modulating the substrate preference in cytochrome P450BsβHI

An active site is normally located inside enzymes, hence substrates should go through a tunnel to access the active site. Tunnel engineering is a powerful strategy for refining the catalytic properties of enzymes. Here, P450BsβHI (Q85H/V170I) derived from hydroxylase P450Bsβ from Bacillus subtilis was chosen as the study model, which is reported as a potential decarboxylase. However, this enzyme showed low decarboxylase activity towards long-chain fatty acids. Here, a tunnel engineering campaign was performed for modulating the substrate preference and improving the decarboxylation activity of P450BsβHI. The finally obtained BsβHI-F79A variant had a 15.2-fold improved conversion for palmitic acid; BsβHI-F173V variant had a 3.9-fold improved conversion for pentadecanoic acid. The study demonstrates how the substrate preference can be modulated by tunnel engineering strategy.

The structure of a potassium-selective ion channel reveals a hydrophobic gate regulating ion permeation

Protein dynamics are essential to function. One example of this is the various gating mechanisms within ion channels, which are transmembrane proteins that act as gateways into the cell. Typical ion channels switch between an open and closed state via a conformational transition which is often triggered by an external stimulus, such as ligand binding or pH and voltage differences. The atomic resolution structure of a potassium-selective ion channel named NaK2K has allowed us to observe that a hydro-phobic residue at the bottom of the selectivity filter, Phe92, appears in dual conformations. One of the two conformations of Phe92 restricts the diameter of the exit pore around the selectivity filter, limiting ion flow through the channel, while the other conformation of Phe92 provides a larger-diameter exit pore from the selectivity filter. Thus, it can be concluded that Phe92 acts as a hydro-phobic gate, regulating the flow of ions through the selectivity filter.

Concentration effects on capture rate and translocation configuration of nanopore-based DNA detection

Nanopore is a kind of powerful tool to detect single molecules and investigate fundamental biological processes. In biological cells or real detection systems, concentration of DNA molecules is various. Here, we report an experimental study of the effects of DNA concentration on capture rate and translocation configuration with different sized nanopores and applied voltages. Three classes of DNA translocation configurations have been observed including linear translocation, folded translocation, and cotranslocation. In the case of relatively large sized nanopore or high applied voltage, considerable cotranslocation events have been detected. The percentage of cotranslocation events also increases with DNA concentration, which leads to the relationship between capture rate and DNA concentration deviates from linearity. Therefore, in order to reflect the number of translocation molecules accurately, the capture rate should be corrected by double-counting cotranslocation events. These results will provide a valuable reference for the design of nanopore sensors.

Coarse-grain simulations on NMR conformational ensembles highlight functional residues in proteins

Dynamics are a key feature of protein function, and this is especially true of gating residues, which occupy cavity or tunnel lining positions in the protein structure, and will reversibly switch between open and closed conformations in order to control the diffusion of small molecules within a protein's internal matrix. Earlier work on globins and hydrogenases have shown that these gating residues can be detected using a multiscale scheme combining all-atom classic molecular dynamics simulations and coarse-grain calculations of the resulting conformational ensemble mechanical properties. Here, we show that the structural variations observed in the conformational ensembles produced by NMR spectroscopy experiments are sufficient to induce noticeable mechanical changes in a protein, which in turn can be used to identify residues important for function and forming a mechanical nucleus in the protein core. This new approach, which combines experimental data and rapid coarse-grain calculations and no longer needs to resort to time-consuming all-atom simulations, was successfully applied to five different protein families.

Structure and Dynamics of Stacking Interactions in an Antibody Binding Site

For years, antibodies (Abs) have been used as a paradigm for understanding how protein structure contributes to molecular recognition. However, with the ability to evolve Abs that recognize specific chromophores, they also have great potential as models for how protein dynamics contribute to molecular recognition. We previously raised murine Abs to different chromophores and, with the use of three-pulse photon echo peak shift spectroscopy, demonstrated that the immune system is capable of producing Abs with widely varying flexibility. We now report the characterization of the complexes formed between two Abs, 5D11 and 10A6, and the chromophoric ligand that they were evolved to recognize, 8-methoxypyrene-1,3,6-trisulfonic acid (MPTS). The sequences of the Ab genes indicate that they evolved from a common precursor. We also used a variety of spectroscopic methods to probe the photophysics and dynamics of the Ab-MPTS complexes and found that they are similar to each other but distinct from previously characterized anti-MPTS Abs. Structural studies revealed that this difference likely results from a unique mode of binding in which MPTS is sandwiched between the side chain of Phe^H98, which interacts with the chromophore via T-stacking, and the side chain of Trp^L91, which interacts with the chromophore via parallel stacking. The T-stacking interaction appears to mediate relaxation on the picosecond time scale, while the parallel stacking appears to mediate relaxation on an ultrafast, femtosecond time scale, which dominates the response. The anti-MPTS Abs thus not only demonstrate the simultaneous use of the two limiting modes of stacking for molecular recognition, but also provide a unique opportunity to characterize how dynamics might contribute to molecular recognition. Both types of stacking are common in proteins and protein complexes where they may similarly contribute to dynamics and molecular recognition.

An Energetically Favorable Ligand Entrance Gate of a Multidrug Transporter Revealed by Partial Nudged Elastic Band Simulations

P-glycoprotein (P-gp) is a multidrug transporter, which harnesses the chemical energy of ATP to power the efflux of diverse chemotherapeutics out of cells and thus contributes to the development of multidrug resistance (MDR) in cancer. It has been proved that the ligand-binding pocket of P-gp is located at the transmembrane domains (TMDs). However, the access of ligands into the binding pocket remains to be elucidated, which definitely hinder the development of P-gp inhibitors. Herein, the access pathways of a well-known substrate rhodamine-123 and a cyclopeptide inhibitor QZ-Leu were characterized by time-independent partial nudged elastic band (PNEB) simulations. The decreasing free energies along the PNEB-optimized access pathway indicated that TM4/6 cleft may be an energetically favorable entrance gate for ligand entry into the binding pocket of P-gp. The results can be reconciled with a range of experimental studies, further corroborating the reliability of the gate revealed by computational simulations. Our atomic level description of the ligand access pathway provides valuable insights into the gating mechanism for drug uptake and transport by P-gp and other multidrug transporters.

•

P-gp contributes to the development of multidrug resistance in cancer.

•

The entrance of drugs into P-gp binding pocket has yet to be elucidated.

•

An energetically favorable entrance gate was revealed by PNEB simulations.

•

The computational results were reconciled with the experimental data.

•

The atomic simulations provide insights into the gating mechanism of P-gp.

Lipid Headgroup Charge and Acyl Chain Composition Modulate Closure of Bacterial β-Barrel Channels

The outer membrane of Gram-negative bacteria contains β-barrel proteins that form high-conducting ion channels providing a path for hydrophilic molecules, including antibiotics. Traditionally, these proteins have been considered to exist only in an open state so that regulation of outer membrane permeability was accomplished via protein expression. However, electrophysiological recordings show that β-barrel channels respond to transmembrane voltages by characteristically switching from a high-conducting, open state, to a so-called 'closed' state, with reduced permeability and possibly exclusion of large metabolites. Here, we use the bacterial porin OmpF from E. coli as a model system to gain insight on the control of outer membrane permeability by bacterial porins through the modulation of their open state. Using planar bilayer electrophysiology, we perform an extensive study of the role of membrane lipids in the OmpF channel closure by voltage. We pay attention not only to the effects of charges in the hydrophilic lipid heads but also to the contribution of the hydrophobic tails in the lipid-protein interactions. Our results show that gating kinetics is governed by lipid characteristics so that each stage of a sequential closure is different from the previous one, probably because of intra- or intermonomeric rearrangements.

Molecular Dynamics Simulations of the SNARE Complex

Molecular dynamics (MD) simulations enable in silico investigations of the dynamic behavior of proteins and protein complexes. Here, we describe MD simulations of the SNARE complex and its interactions with the neuronal protein complexin. Complexin is an effector of neuronal secretion that inhibits spontaneous fusion and is thought to clamp the fusion process via the interactions with the SNARE complex. We describe MD simulations of the SNARE complex alone and bound to complexin. The MD simulations under external forces imitating the repulsion between lipid bilayers enabled us to investigate unraveling and assembly of the SNARE complex.

An odorant receptor from Anopheles sinensis in China is sensitive to oviposition attractants

Background

Anopheles sinensis is an important vector for the spread of malaria in China. Olfactory-related behaviours, particularly oviposition site seeking, offer opportunities for disrupting the disease-transmission process.

Results

This is the first report of the identification and characterization of AsinOrco and AsinOR10 in An. sinensis. AsinOrco and AsinOR10 share 97.49% and 90.37% amino acid sequence identity, respectively, with related sequences in Anopheles gambiae. A functional analysis demonstrated that AsinOrco- and AsinOR10-coexpressing HEK293 cells were highly sensitive to 3-methylindole, but showed no significant differences in response to other test odorants when compared to DMSO.

Conclusions

AsinOrco was characterized as a new member of the Orco ortholog subfamily. AsinOR10, which appears to be a member of the OR2-10 subfamily, is directly involved in identification of oviposition sites. This finding will help to elucidate the molecular mechanisms underlying olfactory signaling in An. sinensis and provide many more molecular targets for eco-friendly pest control.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2501-4) contains supplementary material, which is available to authorized users.

Selective Permeability of Carboxysome Shell Pores to Anionic Molecules

Carboxysomes are closed polyhedral cellular microcompartments that increase the efficiency of carbon fixation in autotrophic bacteria. Carboxysome shells consist of small proteins that form hexameric units with semipermeable central pores containing binding sites for anions. This feature is thought to selectively allow access to RuBisCO enzymes inside the carboxysome by HCO3- (the dominant form of CO2 in the aqueous solution at pH 7.4) but not O2, which leads to a nonproductive reaction. To test this hypothesis, here we use molecular dynamics simulations to characterize the energetics and permeability of CO2, O2, and HCO3- through the central pores of two different shell proteins, namely, CsoS1A of α-carboxysome and CcmK4 of β-carboxysome shells. We find that the central pores are in fact selectively permeable to anions such as HCO3-, as predicted by the model.

BALCONY: an R package for MSA and functional compartments of protein variability analysis

BACKGROUND

Here, we present an R package for entropy/variability analysis that facilitates prompt and convenient data extraction, manipulation and visualization of protein features from multiple sequence alignments. BALCONY can work with residues dispersed across a protein sequence and map them on the corresponding alignment of homologous protein sequences. Additionally, it provides several entropy and variability scores that indicate the conservation of each residue.

RESULTS

Our package allows the user to visualize evolutionary variability by locating the positions most likely to vary and to assess mutation candidates in protein engineering.

CONCLUSION

In comparison to other R packages BALCONY allows conservation/variability analysis in context of protein structure with linkage of the appropriate metrics with physicochemical features of user choice.

AVAILABILITY

CRAN project page: https://cran.r-project.org/package=BALCONY and our website: http://www.tunnelinggroup.pl/software/ for major platforms: Linux/Unix, Windows and Mac OS X.

Expressional and functional interactions of two Apis cerana cerana olfactory receptors

Apis cerana cerana relies on its sensitive olfactory system to perform foraging activities in the surrounding environment. Olfactory receptors (ORs) are a primary requirement for odorant recognition and coding. However, the molecular recognition of volatile compounds with ORs in A. cerana cerana is still not clear. Hence, in the present study, we achieved transient transfection and cell surface expression of A. cerana cerana ORs (AcerOr1 and AcerOr2; AcerOr2 is orthologous to the co-receptor) in Spodoptera frugiperda (Sf9) cells. AcerOr2 narrowly responded to N-(4-ethylphenyl)-2-((4-ethyl-5-(3-pyridinyl)-4H-1,2,4-triazol-3-yl) thio) acetamide (VUAA1), whereas AcerOr1 was sensitive to eugenol, lauric acid, ocimene, 1-nonanol, linolenic acid, hexyl acetate, undecanoic acid, 1-octyl alcohol, and nerol. Of the compounds tested, AcerOr1 showed the highest sensitivity to these odorants with EC₅₀ values of 10⁻⁷ and 10⁻⁸ M, and AcerOr2 recognized VUAA1 with higher sensitivity [EC₅₀ = (6.621 ± 0.26) × 10⁻⁸]. These results indicate that AcerOr2 is an essential gene for olfactory signaling, and AcerOr1 is a broadly tuned receptor. We discovered ligands that were useful for probing receptor activity during odor stimulation and validated three of them by electroantennography. The response increased with concentration of the odorant. The present study provides insight into the mechanism of olfactory discrimination in A. cerana cerana.

Hydrophobic pore gates regulate ion permeation in polycystic kidney disease 2 and 2L1 channels

PKD2 and PKD1 genes are mutated in human autosomal dominant polycystic kidney disease. PKD2 can form either a homomeric cation channel or a heteromeric complex with the PKD1 receptor, presumed to respond to ligand(s) and/or mechanical stimuli. Here, we identify a two-residue hydrophobic gate in PKD2L1, and a single-residue hydrophobic gate in PKD2. We find that a PKD2 gain-of-function gate mutant effectively rescues PKD2 knockdown-induced phenotypes in embryonic zebrafish. The structure of a PKD2 activating mutant F604P by cryo-electron microscopy reveals a π- to α-helix transition within the pore-lining helix S6 that leads to repositioning of the gate residue and channel activation. Overall the results identify hydrophobic gates and a gating mechanism of PKD2 and PKD2L1.

Mutations in the cation channel PKD2 cause human autosomal dominant polycystic kidney disease but its channel function and gating mechanism are poorly understood. Here authors study PKD2 using electrophysiology and cryo-EM, which identifies hydrophobic gates and proposes a gating mechanism for PKD2.

Active site gate of M32 carboxypeptidases illuminated by crystal structure and molecular dynamics simulations

Enzyme gates are important dynamic features that regulate function. Study of these features is critical for understanding of enzyme mechanism. In this study, the active-site gate of M32 carboxypeptidases (M32CP) is illuminated. Only a handful of members of this family have been structurally and functionally characterized and various aspects of their activity and mechanism are yet not clarified. Here, crystal structure of putative M32CP from Deinococcus radiodurans (M32dr) was solved to 2.4Å resolution. Enzymatic assays confirmed its identity as a carboxypeptidase. Open and relatively closed conformations observed in the structure provided supporting evidence for previously hypothesized hinge motion in this family of enzymes. Molecular dynamics simulations of 1.5μs displayed distinct open and closed conformations revealing amplitude of the motion to be beyond what was observed in the crystal structure. Hinge region and anchoring region of this shell-type gate were identified. A small displacement of 3Å and a helical tilt of 9° propagated by the hinge region translates into a 10Å motion at the top of the gate. The dynamics of the gate was supported by our mutagenesis experiment involving formation of disulphide bond across helices of the gate. The nearly inactive mutant enzyme showed 65-fold increase in the enzymatic activity in presence of reducing agent. Further, while a previously proposed structural basis would have led to its classification in subfamily II, experimentally observed substrate length restriction places M32dr in subfamily I of M32CPs.

Identification of Intrahelical Bifurcated H-Bonds as a New Type of Gate in K+ Channels

Gating of ion channels is based on structural transitions between open and closed states. To uncover the chemical basis of individual gates, we performed a comparative experimental and computational analysis between two K⁺ channels, Kcv_S and Kcv_NTS. These small viral encoded K⁺ channel proteins, with a monomer size of only 82 amino acids, resemble the pore module of all complex K⁺ channels in terms of structure and function. Even though both proteins share about 90% amino acid sequence identity, they exhibit different open probabilities with ca. 90% in Kcv_NTS and 40% in Kcv_S. Single channel analysis, mutational studies and molecular dynamics simulations show that the difference in open probability is caused by one long closed state in Kcv_S. This state is structurally created in the tetrameric channel by a transient, Ser mediated, intrahelical hydrogen bond. The resulting kink in the inner transmembrane domain swings the aromatic rings from downstream Phes in the cavity of the channel, which blocks ion flux. The frequent occurrence of Ser or Thr based helical kinks in membrane proteins suggests that a similar mechanism could also occur in the gating of other ion channels.

A Gate Hinge Controls the Epithelial Calcium Channel TRPV5

TRPV5 is unique within the large TRP channel family for displaying a high Ca²⁺ selectivity together with Ca²⁺-dependent inactivation. Our study aims to uncover novel insights into channel gating through in-depth structure-function analysis. We identify an exceptional tryptophan (W583) at the terminus of the intracellular pore that is unique for TRPV5 (and TRPV6). A combination of site-directed mutagenesis, biochemical and electrophysiological analysis, together with homology modeling, demonstrates that W583 is part of the gate for Ca²⁺ permeation. The W583 mutants show increased cell death due to profoundly enhanced Ca²⁺ influx, resulting from altered channel function. A glycine residue above W583 might act as flexible linker to rearrange the tryptophan gate. Furthermore, we hypothesize functional crosstalk between the pore region and carboxy terminus, involved in Ca²⁺-calmodulin-mediated inactivation. This study proposes a unique channel gating mechanism and delivers detailed molecular insight into the Ca²⁺ permeation pathway that can be extrapolated to other Ca²⁺-selective channels.

Nanobodies to Study G Protein-Coupled Receptor Structure and Function

Ligand-induced activation of G protein-coupled receptors (GPCRs) is a key mechanism permitting communication between cells and organs. Enormous progress has recently elucidated the structural and dynamic features of GPCR transmembrane signaling. Nanobodies, the recombinant antigen-binding fragments of camelid heavy-chain-only antibodies, have emerged as important research tools to lock GPCRs in particular conformational states. Active-state stabilizing nanobodies have elucidated several agonist-bound structures of hormone-activated GPCRs and have provided insight into the dynamic character of receptors. Nanobodies have also been used to stabilize transient GPCR transmembrane signaling complexes, yielding the first structural insights into GPCR signal transduction across the cellular membrane. Beyond their in vitro uses, nanobodies have served as conformational biosensors in living systems and have provided novel ways to modulate GPCR function. Here, we highlight several examples of how nanobodies have enabled the study of GPCR function and give insights into potential future uses of these important tools.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Functional analysis of Orco and odorant receptors in odor recognition in Aedes albopictus

Background

Aedes albopictus is a globally invasive mosquito and a major vector of arboviruses, including dengue, Zika and Chikungunya. Olfactory-related behaviors, particularly host-seeking, offer opportunities to disrupt the disease-transmission process. A better understanding of odorant receptors (ORs) may assist in explaining host selection and location, and contribute to novel strategy of vector control.

Methods

Based on previous prediction of 158 putative odorant receptors by Ae. albopictus genome analysis, 29 AalORs were selected for tissue-specific expression profiles in the present study. AalOrco (AalOR7), AalOR10 and AalOR88, highly expressed in female olfactory tissues, were chosen for further structure predictions as well as functional validation including calcium imaging assay in human embryonic kidney (HEK293) cells and RNA interference assay in Ae. albopictus. We also conducted electrophysiological and behavioral assays in mosquitoes after RNA interference of the three genes to determine their roles in host-seeking.

Results

The results support previous conclusions that individual conventional (ORXs) and Orco can form heteromeric complexes to recognize odorants and respond to components of human volatiles in HEK293 cells. The reduction of AalOrco transcript levels led to a significant decrease in host-seeking and confusion in host preference. In contrast, AalOR10 and AalOR88 knockdown mosquitoes showed no significant behavioral differences compared with controls. The functions of conventional ORs at least AalOR10 and AalOR88 are abolished with inhibited expression of the Orco gene orthologs, along with the concomitant relevant olfactory behavior.

Conclusions

Combining structural and functional data, we conclude that the product of the Orco gene in this mosquito is crucial for transmitting olfactory signaling and conventional ORs contribute directly to odorant recognition. Our results provide insight into the linkage between odorant receptors and host-seeking in this important vector species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1644-9) contains supplementary material, which is available to authorized users.

All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

Structural Changes Fundamental to Gating of the Cystic Fibrosis Transmembrane Conductance Regulator Anion Channel Pore

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), an epithelial cell anion channel. Potentiator drugs used in the treatment of cystic fibrosis act on the channel to increase overall channel function, by increasing the stability of its open state and/or decreasing the stability of its closed state. The structure of the channel in either the open state or the closed state is not currently known. However, changes in the conformation of the protein as it transitions between these two states have been studied using functional investigation and molecular modeling techniques. This review summarizes our current understanding of the architecture of the transmembrane channel pore that controls the movement of chloride and other small anions, both in the open state and in the closed state. Evidence for different kinds of changes in the conformation of the pore as it transitions between open and closed states is described, as well as the mechanisms by which these conformational changes might be controlled to regulate normal channel gating. The ways that key conformational changes might be targeted by small compounds to influence overall CFTR activity are also discussed. Understanding the changes in pore structure that might be manipulated by such small compounds is key to the development of novel therapeutic strategies for the treatment of cystic fibrosis.

Towards an understanding of the structural basis for insect olfaction by odorant receptors

Insects have co-opted a unique family of seven transmembrane proteins for odour sensing. Odorant receptors are believed to have evolved from gustatory receptors somewhere at the base of the Hexapoda and have expanded substantially to become the dominant class of odour recognition elements within the Insecta. These odorant receptors comprise an obligate co-receptor, Orco, and one of a family of highly divergent odorant "tuning" receptors. The two subunits are thought to come together at some as-yet unknown stoichiometry to form a functional complex that is capable of both ionotropic and metabotropic signalling. While there are still no 3D structures for these proteins, site-directed mutagenesis, resonance energy transfer, and structural modelling efforts, all mainly on Drosophila odorant receptors, are beginning to inform hypotheses of their structures and how such complexes function in odour detection. Some of the loops, especially the second extracellular loop that has been suggested to form a lid over the binding pocket, and the extracellular regions of some transmembrane helices, especially the third and to a less extent the sixth and seventh, have been implicated in ligand recognition in tuning receptors. The possible interaction between Orco and tuning receptor subunits through the final intracellular loop and the adjacent transmembrane helices is thought to be important for transducing ligand binding into receptor activation. Potential phosphorylation sites and a calmodulin binding site in the second intracellular loop of Orco are also thought to be involved in regulating channel gating. A number of new methods have recently been developed to express and purify insect odorant receptor subunits in recombinant expression systems. These approaches are enabling high throughput screening of receptors for agonists and antagonists in cell-based formats, as well as producing protein for the application of biophysical methods to resolve the 3D structure of the subunits and their complexes.

Temporal evolution of helix hydration in a light-gated ion channel correlates with ion conductance

The discovery of channelrhodopsins introduced a new class of light-gated ion channels, which when genetically encoded in host cells resulted in the development of optogenetics. Channelrhodopsin-2 from Chlamydomonas reinhardtii, CrChR2, is the most widely used optogenetic tool in neuroscience. To explore the connection between the gating mechanism and the influx and efflux of water molecules in CrChR2, we have integrated light-induced time-resolved infrared spectroscopy and electrophysiology. Cross-correlation analysis revealed that ion conductance tallies with peptide backbone amide I vibrational changes at 1,665(-) and 1,648(+) cm(-1). These two bands report on the hydration of transmembrane α-helices as concluded from vibrational coupling experiments. Lifetime distribution analysis shows that water influx proceeded in two temporally separated steps with time constants of 10 μs (30%) and 200 μs (70%), the latter phase concurrent with the start of ion conductance. Water efflux and the cessation of the ion conductance are synchronized as well, with a time constant of 10 ms. The temporal correlation between ion conductance and hydration of helices holds for fast (E123T) and slow (D156E) variants of CrChR2, strengthening its functional significance.

Gated molecular baskets

In this review, we describe the construction of gated molecular baskets, discuss their mechanism of action in regulating the exchange of guests and illustrate the potential of these concave hosts to act as catalysts for controlling chemical reactions. Importantly, a number of computational and experimental studies have suggested that gated baskets ought to unfold their gates at the rim for permitting the passage of guests to/from their inner space. These dynamic hosts are therefore offered as useful models for investigating the process of gating in artificial systems. Furthermore, gated baskets should permit examining the benefit of controlling the rate by which reactants access a gated catalyst for promoting chemical reactions occurring in its confined space.

CHEXVIS: a tool for molecular channel extraction and visualization

Background

Understanding channel structures that lead to active sites or traverse the molecule is important in the study of molecular functions such as ion, ligand, and small molecule transport. Efficient methods for extracting, storing, and analyzing protein channels are required to support such studies. Further, there is a need for an integrated framework that supports computation of the channels, interactive exploration of their structure, and detailed visual analysis of their properties.

Results

We describe a method for molecular channel extraction based on the alpha complex representation. The method computes geometrically feasible channels, stores both the volume occupied by the channel and its centerline in a unified representation, and reports significant channels. The representation also supports efficient computation of channel profiles that help understand channel properties. We describe methods for effective visualization of the channels and their profiles. These methods and the visual analysis framework are implemented in a software tool, ChExVis. We apply the method on a number of known channel containing proteins to extract pore features. Results from these experiments on several proteins show that ChExVis performance is comparable to, and in some cases, better than existing channel extraction techniques. Using several case studies, we demonstrate how ChExVis can be used to study channels, extract their properties and gain insights into molecular function.

Conclusion

ChExVis supports the visual exploration of multiple channels together with their geometric and physico-chemical properties thereby enabling the understanding of the basic biology of transport through protein channels. The ChExVis web-server is freely available at http://vgl.serc.iisc.ernet.in/chexvis/. The web-server is supported on all modern browsers with latest Java plug-in.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0545-9) contains supplementary material, which is available to authorized users.

Crystal structure of snake venom acetylcholinesterase in complex with inhibitory antibody fragment Fab410 bound at the peripheral site: evidence for open and closed states of a back door channel

The acetylcholinesterase found in the venom of Bungarus fasciatus (BfAChE) is produced as a soluble, non-amphiphilic monomer with a canonical catalytic domain but a distinct C terminus compared with the other vertebrate enzymes. Moreover, the peripheral anionic site of BfAChE, a surface site located at the active site gorge entrance, bears two substitutions altering sensitivity to cationic inhibitors. Antibody Elec410, generated against Electrophorus electricus acetylcholinesterase (EeAChE), inhibits EeAChE and BfAChE by binding to their peripheral sites. However, both complexes retain significant residual catalytic activity, suggesting incomplete gorge occlusion by bound antibody and/or high frequency back door opening. To explore a novel acetylcholinesterase species, ascertain the molecular bases of inhibition by Elec410, and document the determinants and mechanisms for back door opening, we solved a 2.7-Å resolution crystal structure of natural BfAChE in complex with antibody fragment Fab410. Crystalline BfAChE forms the canonical dimer found in all acetylcholinesterase structures. Equally represented open and closed states of a back door channel, associated with alternate positions of a tyrosine phenol ring at the active site base, coexist in each subunit. At the BfAChE molecular surface, Fab410 is seated on the long Ω-loop between two N-glycan chains and partially occludes the gorge entrance, a position that fully reflects the available mutagenesis and biochemical data. Experimentally based flexible molecular docking supports a similar Fab410 binding mode onto the EeAChE antigen. These data document the molecular and dynamic peculiarities of BfAChE with high frequency back door opening, and the mode of action of Elec410 as one of the largest peptidic inhibitors targeting the acetylcholinesterase peripheral site.

Probing protein multidimensional conformational fluctuations by single-molecule multiparameter photon stamping spectroscopy

Conformational motions of proteins are highly dynamic and intrinsically complex. To capture the temporal and spatial complexity of conformational motions and further to understand their roles in protein functions, an attempt is made to probe multidimensional conformational dynamics of proteins besides the typical one-dimensional FRET coordinate or the projected conformational motions on the one-dimensional FRET coordinate. T4 lysozyme hinge-bending motions between two domains along α-helix have been probed by single-molecule FRET. Nevertheless, the domain motions of T4 lysozyme are rather complex involving multiple coupled nuclear coordinates and most likely contain motions besides hinge-bending. It is highly likely that the multiple dimensional protein conformational motions beyond the typical enzymatic hinged-bending motions have profound impact on overall enzymatic functions. In this report, we have developed a single-molecule multiparameter photon stamping spectroscopy integrating fluorescence anisotropy, FRET, and fluorescence lifetime. This spectroscopic approach enables simultaneous observations of both FRET-related site-to-site conformational dynamics and molecular rotational (or orientational) motions of individual Cy3-Cy5 labeled T4 lysozyme molecules. We have further observed wide-distributed rotational flexibility along orientation coordinates by recording fluorescence anisotropy and simultaneously identified multiple intermediate conformational states along FRET coordinate by monitoring time-dependent donor lifetime, presenting a whole picture of multidimensional conformational dynamics in the process of T4 lysozyme open-close hinge-bending enzymatic turnover motions under enzymatic reaction conditions. By analyzing the autocorrelation functions of both lifetime and anisotropy trajectories, we have also observed the dynamic and static inhomogeneity of T4 lysozyme multidimensional conformational fluctuation dynamics, providing a fundamental understanding of the enzymatic reaction turnover dynamics associated with overall enzyme as well as the specific active-site conformational fluctuations that are not identifiable and resolvable in the conventional ensemble-averaged experiment.

General rules for the arrangements and gating motions of pore-lining helices in homomeric ion channels

The pore-lining helix (PLH) bundles are central to the function of all ion channels, as their conformational rearrangements dictate channel gating. Here we explore all plausible oligomeric arrangements of the PLH bundles within homomeric ion channels by building models using generic restraints. In particular, the distance between two neighbouring PLHs was bounded both below and above in order to avoid steric clash and allow proper packing. The resulting models provide a theoretical representation of the accessible space for oligomeric arrangements. While the represented space is confined, it encompasses nearly all the ion channel PLH bundles for which the structures are currently known. For a multitude of channels, gating models suggested by paths within the confined accessible space are in qualitative agreement with those established in previous structural and computational studies.

Rearrangements of the pore-lining helix (PLH) bundles of ion channels are central to their gating mechanisms. Here, Dai et al. use a modelling approach to define the general rules that govern the arrangements and gating motions of the PLHs in homomeric ion channels.

Channel-facilitated diffusion boosted by particle binding at the channel entrance

We investigate single-file diffusion of Brownian particles in arrays of closely confining microchannels permeated by a variety of attractive optical potentials and connecting two baths with equal particle concentration. We simultaneously test free diffusion in the channel, diffusion in optical traps coupled in the center of the channel, and diffusion in traps extending into the baths. We found that both classes of attractive optical potentials enhance the translocation rate through the channel with respect to free diffusion. Surprisingly, for the latter class of potentials we measure a 40-fold enhancement in the translocation rate with respect to free diffusion and find a sublinear power law dependence of the translocation rate on the average number of particles in the channel. Our results reveal the function of particle binding at the channel entrances for diffusive transport and open the way to a better understanding of membrane transport and design of synthetic membranes with enhanced diffusion rate.

Insights into water coordination associated with the CuII/CuI electron transfer at a biomimetic Cu centre

The coordination chemistry of an aqua Cu complex was investigated in non-coordinating solvents and in ionic liquids.

Homotropic cooperativity from the activation pathway of the allosteric ligand-responsive regulatory trp RNA-binding attenuation protein

The trp RNA-binding attenuation protein (TRAP) assembles into an 11-fold symmetric ring that regulates transcription and translation of trp-mRNA in bacilli via heterotropic allosteric activation by the amino acid tryptophan (Trp). Whereas nuclear magnetic resonance studies have revealed that Trp-induced activation coincides with both microsecond to millisecond rigidification and local structural changes in TRAP, the pathway of binding of the 11 Trp ligands to the TRAP ring remains unclear. Moreover, because each of 11 bound Trp molecules is completely surrounded by protein, its release requires flexibility of Trp-bound (holo) TRAP. Here, we used stopped-flow fluorescence to study the kinetics of Trp binding by Bacillus stearothermophilus TRAP over a range of temperatures and observed well-separated kinetic steps. These data were analyzed using nonlinear least-squares fitting of several two- and three-step models. We found that a model with two binding steps best describes the data, although the structural equivalence of the binding sites in TRAP implies a fundamental change in the time-dependent structure of the TRAP rings upon Trp binding. Application of the two-binding step model reveals that Trp binding is much slower than the diffusion limit, suggesting a gating mechanism that depends on the dynamics of apo TRAP. These data also reveal that dissociation of Trp from the second binding mode is much slower than after the first Trp binding mode, revealing insight into the mechanism for positive homotropic allostery, or cooperativity. Temperature-dependent analyses reveal that both binding modes imbue increases in bondedness and order toward a more compressed active state. These results provide insight into mechanisms of cooperative TRAP activation and underscore the importance of protein dynamics for ligand binding, ligand release, protein activation, and allostery.

Uncovering a dynamically formed substrate access tunnel in carbon monoxide dehydrogenase/acetyl-CoA synthase

The transport of small ligands to active sites of proteins is the basis of vital processes in biology such as enzymatic catalysis and cell signaling, but also of more destructive ones including enzyme inhibition and oxidative damage. Here, we show how a diffusion-reaction model solved by means of molecular dynamics and density functional theory calculations provides novel insight into the transport of small ligands in proteins. In particular, we unravel the existence of an elusive, dynamically formed gas channel, which CO2 takes to diffuse from the solvent to the active site (C-cluster) of the bifunctional multisubunit enzyme complex carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS). Two cavities forming this channel are temporarily created by protein fluctuations and are not apparent in the X-ray structures. The ligand transport is controlled by two residues at the end of this tunnel, His113 and His116, and occurs on the same time scale on which chemical binding to the active site takes place (0.1-1 ms), resulting in an overall binding rate on the second time scale. We find that upon reduction of CO2 to CO, the newly formed Fe-hydroxy ligand greatly strengthens the hydrogen-bond network, preventing CO from exiting the protein through the same way that CO2 takes to enter the protein. This is the basis for directional transport of CO from the production site (C-cluster of CODH subunit) to the utilization site (A-cluster of ACS subunit). In view of these results, a general picture emerges of how large proteins guide small ligands toward their active sites.