Assessing the local structural quality of transmembrane protein models using statistical potentials (QMEANBrane)

MyPMFs: a simple tool for creating statistical potentials to assess protein structural models

Evaluating the model quality of protein structures that evolve in environments with particular physicochemical properties requires scoring functions that are adapted to their specific residue compositions and/or structural characteristics. Thus, computational methods developed for structures from the cytosol cannot work properly on membrane or secreted proteins. Here, we present MyPMFs, an easy-to-use tool that allows users to train statistical potentials of mean force (PMFs) on the protein structures of their choice, with all parameters being adjustable. We demonstrate its use by creating an accurate statistical potential for transmembrane protein domains. We also show its usefulness to study the influence of the physical environment on residue interactions within protein structures. Our open-source software is freely available for download at https://github.com/bibip-impmc/mypmfs.

On the contributing role of the transmembrane domain for subunit-specific sensitivity of integrin activation

Integrins are α/β heterodimeric transmembrane adhesion receptors. Evidence exists that their transmembrane domain (TMD) separates upon activation. Subunit-specific differences in activation sensitivity of integrins were reported. However, whether sequence variations in the TMD lead to differential TMD association has remained elusive. Here, we show by molecular dynamics simulations and association free energy calculations on TMDs of integrin αIIbβ3, αvβ3, and α5β1 that αIIbβ3 TMD is most stably associated; this difference is related to interaction differences across the TMDs. The order of TMD association stability is paralleled by the basal activity of these integrins, which suggests that TMD differences can have a decisive effect on integrin conformational free energies. We also identified a specific order of clasp disintegration upon TMD dissociation, which suggests that the closed state of integrins may comprise several microstates. Our results provide unprecedented insights into a possibly contributing role of TMD towards subunit-specific sensitivity of integrin activation.

Disease-associated missense mutations in GluN2B subunit alter NMDA receptor ligand binding and ion channel properties

Genetic and bioinformatic analyses have identified missense mutations in GRIN2B encoding the NMDA receptor GluN2B subunit in autism, intellectual disability, Lennox Gastaut and West Syndromes. Here, we investigated several such mutations using a near-complete, hybrid 3D model of the human NMDAR and studied their consequences with kinetic modelling and electrophysiology. The mutants revealed reductions in glutamate potency; increased receptor desensitisation; and ablation of voltage-dependent Mg2+ block. In addition, we provide new views on Mg2+ and NMDA channel blocker binding sites. We demonstrate that these mutants have significant impact on excitatory transmission in developing neurons, revealing profound changes that could underlie their associated neurological disorders. Of note, the NMDAR channel mutant GluN2BV618G unusually allowed Mg2+ permeation, whereas nearby N615I reduced Ca2+ permeability. By identifying the binding site for an NMDAR antagonist that is used in the clinic to rescue gain-of-function phenotypes, we show that drug binding may be modified by some GluN2B disease-causing mutations.

A multidrug ABC transporter with a taste for GTP

During the evolution of cellular bioenergetics, many protein families have been fashioned to match the availability and replenishment in energy supply. Molecular motors and primary transporters essentially need ATP to function while proteins involved in cell signaling or translation consume GTP. ATP-Binding Cassette (ABC) transporters are one of the largest families of membrane proteins gathering several medically relevant members that are typically powered by ATP hydrolysis. Here, a Streptococcus pneumoniae ABC transporter responsible for fluoroquinolones resistance in clinical settings, PatA/PatB, is shown to challenge this concept. It clearly favors GTP as the energy supply to expel drugs. This preference is correlated to its ability to hydrolyze GTP more efficiently than ATP, as found with PatA/PatB reconstituted in proteoliposomes or nanodiscs. Importantly, the ATP and GTP concentrations are similar in S. pneumoniae supporting the physiological relevance of GTP as the energy source of this bacterial transporter.

Diversity of Nicotinic Acetylcholine Receptor Positive Allosteric Modulators Revealed by Mutagenesis and a Revised Structural Model

By combining electrophysiological and computational approaches we have examined a series of positive allosteric modulators (PAMs) acting on the human α7 nicotinic acetylcholine receptor (nAChR). Electrophysiological studies have focused on three α7-selective PAMs (A-867744, TBS-516, and TQS) that display similar effects on wild-type α7 nAChRs. In addition to potentiating agonist-evoked responses, all three compounds reduce receptor desensitization and, consequently, are classed as type II PAMs. Despite having similar effects on wild-type receptors, A-867744 was found to have profoundly differing effects on mutated receptors compared with TBS-516 and TQS, a finding that is consistent with previous studies indicating that A-867744 may have a different mechanism of action compare with other α7-selective type II PAMs. Due to evidence that these PAMs bind within the α7 nAChR transmembrane region, we generated and validated new structural models of α7. Importantly, we have corrected a previously identified error in the transmembrane region of the original cryo-electron microscopy Torpedo model; the only pentameric ligand-gated ion channel imaged in a native lipid membrane. Real-space refinement was used to generate closed and open conformations on which the α7 models were based. Consensus docking with an extended series of PAMs with chemical similarity to A-867744, TBS-516, and TQS suggests that all bind to a broadly similar intersubunit transmembrane site. However, differences in the predicted binding of A-867744, compared with TBS-516 and TQS, may help to explain the distinct functional effects of A-867744. Thus, our revised structural models may provide a useful tool for interpreting functional effects of PAMs.

Functionalized Fullerene Targeting Human Voltage-Gated Sodium Channel, hNav1.7

Mutations of hNa_v1.7 that cause its activities to be enhanced contribute to severe neuropathic pain. Only a small number of hNa_v1.7 specific inhibitors have been identified, most of which interact with the voltage-sensing domain of the voltage-activated sodium ion channel. In our previous computational study, we demonstrated that a [Lys₆]-C₈₄ fullerene binds tightly (affinity of 46 nM) to Na_vAb, the voltage-gated sodium channel from the bacterium Arcobacter butzleri. Here, we extend this work and, using molecular dynamics simulations, demonstrate that the same [Lys₆]-C₈₄ fullerene binds strongly (2.7 nM) to the pore of a modeled human sodium ion channel hNa_v1.7. In contrast, the fullerene binds only weakly to a mutated model of hNa_v1.7 (I1399D) (14.5 mM) and a model of the skeletal muscle hNa_v1.4 (3.7 mM). Comparison of one representative sequence from each of the nine human sodium channel isoforms shows that only hNa_v1.7 possesses residues that are critical for binding the fullerene derivative and blocking the channel pore.

Molecular Modeling of Transporters: From Low Resolution Cryo-Electron Microscopy Map to Conformational Exploration. The Example of TSPO

This chapter describes a protocol to establish a three-dimensional (3D) model of a protein and to explore its conformational landscape. It combines predictions from up-to-date bioinformatics methods with low-resolution experimental data. It also proposes to examine rapidly the dynamics of the protein using molecular dynamics simulations with a coarse-grained force field. Tools for analyzing these trajectories are suggested as well as those for constructing all-atoms models. Thus, starting from a protein sequence and using free software, the user can get important conformational information, which might improve the knowledge about the protein function.

Conjecture Regarding Posttranslational Modifications to the Arabidopsis Type I Proton-Pumping Pyrophosphatase (AVP1)

Agbiotechnology uses genetic engineering to improve the output and value of crops. Altering the expression of the plant Type I Proton-pumping Pyrophosphatase (H⁺-PPase) has already proven to be a useful tool to enhance crop productivity. Despite the effective use of this gene in translational research, information regarding the intracellular localization and functional plasticity of the pump remain largely enigmatic. Using computer modeling several putative phosphorylation, ubiquitination and sumoylation target sites were identified that may regulate Arabidopsis H⁺-PPase (AVP1- Arabidopsis Vacuolar Proton-pump 1) subcellular trafficking and activity. These putative regulatory sites will direct future research that specifically addresses the partitioning and transport characteristics of this pump. We posit that fine-tuning H⁺-PPases activity and cellular distribution will facilitate rationale strategies for further genetic improvements in crop productivity.

Sorting protein decoys by machine-learning-to-rank

Much progress has been made in Protein structure prediction during the last few decades. As the predicted models can span a broad range of accuracy spectrum, the accuracy of quality estimation becomes one of the key elements of successful protein structure prediction. Over the past years, a number of methods have been developed to address this issue, and these methods could be roughly divided into three categories: the single-model methods, clustering-based methods and quasi single-model methods. In this study, we develop a single-model method MQAPRank based on the learning-to-rank algorithm firstly, and then implement a quasi single-model method Quasi-MQAPRank. The proposed methods are benchmarked on the 3DRobot and CASP11 dataset. The five-fold cross-validation on the 3DRobot dataset shows the proposed single model method outperforms other methods whose outputs are taken as features of the proposed method, and the quasi single-model method can further enhance the performance. On the CASP11 dataset, the proposed methods also perform well compared with other leading methods in corresponding categories. In particular, the Quasi-MQAPRank method achieves a considerable performance on the CASP11 Best150 dataset.

Increased efficiency of Campylobacter jejuni N-oligosaccharyltransferase PglB by structure-guided engineering

Conjugate vaccines belong to the most efficient preventive measures against life-threatening bacterial infections. Functional expression of N-oligosaccharyltransferase (N-OST) PglB of Campylobacter jejuni in Escherichia coli enables a simplified production of glycoconjugate vaccines in prokaryotic cells. Polysaccharide antigens of pathogenic bacteria can be covalently coupled to immunogenic acceptor proteins bearing engineered glycosylation sites. Transfer efficiency of PglB_Cj is low for certain heterologous polysaccharide substrates. In this study, we increased glycosylation rates for Salmonella enterica sv. Typhimurium LT2 O antigen (which lacks N-acetyl sugars) and Staphylococcus aureus CP5 polysaccharides by structure-guided engineering of PglB. A three-dimensional homology model of membrane-associated PglB_Cj, docked to the natural C. jejuni N-glycan attached to the acceptor peptide, was used to identify potential sugar-interacting residues as targets for mutagenesis. Saturation mutagenesis of an active site residue yielded the enhancing mutation N311V, which facilitated fivefold to 11-fold increased in vivo glycosylation rates as determined by glycoprotein-specific ELISA. Further rounds of in vitro evolution led to a triple mutant S80R-Q287P-N311V enabling a yield improvement of S. enterica LT2 glycoconjugates by a factor of 16. Our results demonstrate that bacterial N-OST can be tailored to specific polysaccharide substrates by structure-guided protein engineering.

SIB Swiss Institute of Bioinformatics Members. The SIB Swiss Institute of Bioinformatics' resources: focus on curated databases

The SIB Swiss Institute of Bioinformatics (www.isb-sib.ch) provides world-class bioinformatics databases, software tools, services and training to the international life science community in academia and industry. These solutions allow life scientists to turn the exponentially growing amount of data into knowledge. Here, we provide an overview of SIB's resources and competence areas, with a strong focus on curated databases and SIB's most popular and widely used resources. In particular, SIB's Bioinformatics resource portal ExPASy features over 150 resources, including UniProtKB/Swiss-Prot, ENZYME, PROSITE, neXtProt, STRING, UniCarbKB, SugarBindDB, SwissRegulon, EPD, arrayMap, Bgee, SWISS-MODEL Repository, OMA, OrthoDB and other databases, which are briefly described in this article.

Characterization of the differences in the cyclopiazonic acid binding mode to mammalian and P. Falciparum Ca2+ pumps: a computational study

Despite the investments in malaria research, an effective vaccine has not yet been developed and the causative parasites are becoming increasingly resistant to most of the available drugs. PfATP6, the sarco/endoplasmic reticulum Ca2+ pump (SERCA) of P. falciparum, has been recently genetically validated as a potential antimalarial target and cyclopiazonic acid (CPA) has been found to be a potent inhibitor of SERCAs in several organisms, including P. falciparum. In position 263, PfATP6 displays a leucine residue, whilst the corresponding position in the mammalian SERCA is occupied by a glutamic acid. The PfATP6 L263E mutation has been studied in relation to the artemisinin inhibitory effect on P. falciparum and recent studies have provided evidence that the parasite with this mutation is more susceptible to CPA. Here, we characterized, for the first time, the interaction of CPA with PfATP6 and its mammalian counterpart to understand similarities and differences in the mode of binding of the inhibitor to the two Ca2+ pumps. We found that, even though CPA does not directly interact with the residue in position 263, the presence of a hydrophobic residue in this position in PfATP6 rather than a negatively charged one, as in the mammalian SERCA, entails a conformational arrangement of the binding pocket which, in turn, determines a relaxation of CPA leading to a different binding mode of the compound. Our findings highlight differences between the plasmodial and human SERCA CPA-binding pockets that may be exploited to design CPA derivatives more selective toward PfATP6. Proteins 2015; 83:564–574. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

TMFoldRec: a statistical potential-based transmembrane protein fold recognition tool

Background

Transmembrane proteins (TMPs) are the key components of signal transduction, cell-cell adhesion and energy and material transport into and out from the cells. For the deep understanding of these processes, structure determination of transmembrane proteins is indispensable. However, due to technical difficulties, only a few transmembrane protein structures have been determined experimentally. Large-scale genomic sequencing provides increasing amounts of sequence information on the proteins and whole proteomes of living organisms resulting in the challenge of bioinformatics; how the structural information should be gained from a sequence.

Results

Here, we present a novel method, TMFoldRec, for fold prediction of membrane segments in transmembrane proteins. TMFoldRec based on statistical potentials was tested on a benchmark set containing 124 TMP chains from the PDBTM database. Using a 10-fold jackknife method, the native folds were correctly identified in 77 % of the cases. This accuracy overcomes the state-of-the-art methods. In addition, a key feature of TMFoldRec algorithm is the ability to estimate the reliability of the prediction and to decide with an accuracy of 70 %, whether the obtained, lowest energy structure is the native one.

Conclusion

These results imply that the membrane embedded parts of TMPs dictate the TM structures rather than the soluble parts. Moreover, predictions with reliability scores make in this way our algorithm applicable for proteome-wide analyses.

Availability

The program is available upon request for academic use.

Electronic supplementary material

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