AlignMe--a membrane protein sequence alignment web server

Structural Determinants and Functional Significance of Dimerization for Osmosensing Transporter ProP in Escherichia coli

Osmosensing transporter ProP forestalls cellular dehydration by detecting environments with high osmotic pressure and mediating the accumulation of organic osmolytes by bacterial cells. It is composed of 12 transmembrane helices with cytoplasmic N- and C-termini. In Escherichia coli, dimers form when the C-terminal domains of ProP molecules form homodimeric, antiparallel, α-helical coiled coils. No dominant negative effect was detected when inactive and active ProP molecules formed heterodimers in vivo. Purification of ProP in detergent dodecylmaltoside yielded monomers, which were functional after reconstitution in proteoliposomes. With other evidence, this suggests that ProP monomers function independently whether in the monomeric or dimeric state. Amino acid replacements that disrupted or reversed the coiled coil did not prevent in vivo dimerization of ProP detected with a bacterial two-hybrid system. Maleimide labeling detected no osmolality-dependent variation in the reactivities of cysteine residues introduced to transmembrane helix (TM) XII. In contrast, coarse-grained molecular dynamic simulations detected deformation of the lipid around TMs III and VI, on the lipid-exposed protein surface opposite to TM XII. This suggests that the dimer interface of ProP includes the surfaces of TMs III and VI, not of TM XII as previously suggested by crosslinking data. Homology modeling suggested that coiled-coil formation and dimerization via such an interface are not mutually exclusive. In previous work, alterations to the C-terminal coiled coil blocked co-localization of ProP with phospholipid cardiolipin at E. coli cell poles. Thus, dimerization may contribute to ProP targeting, adjust its lipid environment, and hence indirectly modify its osmotic stress response.

Homology Modeling of Transporter Proteins

Membrane transporter proteins are divided into channels/pores and carriers and constitute protein families of physiological and pharmacological importance. Several presently used therapeutic compounds elucidate their effects by targeting membrane transporter proteins, including anti-arrhythmic, anesthetic, antidepressant, anxiolytic and diuretic drugs. The lack of three-dimensional structures of human transporters hampers experimental studies and drug discovery. In this chapter, the use of homology modeling for generating structural models of membrane transporter proteins is reviewed. The increasing number of atomic resolution structures available as templates, together with improvements in methods and algorithms for sequence alignments, secondary structure predictions, and model generation, in addition to the increase in computational power have increased the applicability of homology modeling for generating structural models of transporter proteins. Different pitfalls and hints for template selection, multiple-sequence alignments, generation and optimization, validation of the models, and the use of transporter homology models for structure-based virtual ligand screening are discussed.

AlignMe: an update of the web server for alignment of membrane protein sequences

The AlignMe web server is dedicated to accurately aligning sequences of membrane proteins, a particularly challenging task due to the strong evolutionary divergence and the low compositional complexity of hydrophobic membrane-spanning proteins. AlignMe can create pairwise alignments of either two primary amino acid sequences or two hydropathy profiles. The web server for AlignMe has been continuously available for >10 years, supporting 1000s of users per year. Recent improvements include anchoring, multiple submissions, and structure visualization. Anchoring is the ability to constrain a position in an alignment, which allows expert information about related residues in proteins to be incorporated into an alignment without manual modification. The original web interface to the server limited the user to one alignment per submission, hindering larger scale studies. Now, batches of alignments can be initiated with a single submission. Finally, to provide structural context for the relationship between proteins, sequence similarity can now be mapped onto one or more structures (or structural models) of the proteins being aligned, by links to MutationExplorer, a web-based visualization tool. Together with a refreshed user interface, these features further enhance an important resource in the membrane protein community. The AlignMe web server is freely available at https://www.bioinfo.mpg.de/AlignMe/.

Spontaneous Suppressors against Debilitating Transmembrane Mutants of CaMdr1 Disclose Novel Interdomain Communication via Signature Motifs of the Major Facilitator Superfamily

The Major Facilitator Superfamily (MFS) drug:H+ antiporter CaMdr1, from Candida albicans, is responsible for the efflux of structurally diverse antifungals. MFS members share a common fold of 12−14 transmembrane helices (TMHs) forming two N- and C-domains. Each domain is arranged in a pseudo-symmetric fold of two tandems of 3-TMHs that alternatively expose the drug-binding site towards the inside or the outside of the yeast to promote drug binding and release. MFS proteins show great diversity in primary structure and few conserved signature motifs, each thought to have a common function in the superfamily, although not yet clearly established. Here, we provide new information on these motifs by having screened a library of 64 drug transport-deficient mutants and their corresponding suppressors spontaneously addressing the deficiency. We found that five strains recovered the drug-resistance capacity by expressing CaMdr1 with a secondary mutation. The pairs of debilitating/rescuing residues are distributed either in the same TMH (T127ATMH1- > G140DTMH1) or 3-TMHs repeat (F216ATMH4- > G260ATMH5), at the hinge of 3-TMHs repeats tandems (R184ATMH3- > D235HTMH4, L480ATMH10- > A435TTMH9), and finally between the N- and C-domains (G230ATMH4- > P528HTMH12). Remarkably, most of these mutants belong to the different signature motifs, highlighting a mechanistic role and interplay thought to be conserved among MFS proteins. Results also point to the specific role of TMH11 in the interplay between the N- and C-domains in the inward- to outward-open conformational transition.

Disrupted Choline Clearance and Sustained Acetylcholine Release In Vivo by a Common Choline Transporter Coding Variant Associated with Poor Attentional Control in Humans

Transport of choline via the neuronal high-affinity choline transporter (CHT; SLC5A7) is essential for cholinergic terminals to synthesize and release acetylcholine (ACh). In humans, we previously demonstrated an association between a common CHT coding substitution (rs1013940; Ile89Val) and reduced attentional control as well as attenuated frontal cortex activation. Here, we used a CRISPR/Cas9 approach to generate mice expressing the I89V substitution and assessed, in vivo, CHT-mediated choline transport, and ACh release. Relative to wild-type (WT) mice, CHT-mediated clearance of choline in male and female mice expressing one or two Val89 alleles was reduced by over 80% in cortex and over 50% in striatum. Choline clearance in CHT Val89 mice was further reduced by neuronal inactivation. Deficits in ACh release, 5 and 10 min after repeated depolarization at a low, behaviorally relevant frequency, support an attenuated reloading capacity of cholinergic neurons in mutant mice. The density of CHTs in total synaptosomal lysates and neuronal plasma-membrane-enriched fractions was not impacted by the Val89 variant, indicating a selective impact on CHT function. When challenged with a visual disruptor to reveal attentional control mechanisms, Val89 mice failed to adopt a more conservative response bias. Structural modeling revealed that Val89 may attenuate choline transport by altering conformational changes of CHT that support normal transport rates. Our findings support the view that diminished sustained cholinergic signaling capacity underlies perturbed attentional performance in individuals expressing CHT Val89. The CHT Val89 mouse serves as a valuable model to study heritable risk for cognitive disorders arising from cholinergic dysfunction.SIGNIFICANCE STATEMENT Acetylcholine (ACh) signaling depends on the functional capacity of the neuronal choline transporter (CHT). Previous research demonstrated that humans expressing the common CHT coding variant Val89 exhibit attentional vulnerabilities and attenuated fronto-cortical activation during attention. Here, we find that mice engineered to express the Val89 variant exhibit reduced CHT-mediated choline clearance and a diminished capacity to sustain ACh release. Additionally, Val89 mice lack cognitive flexibility in response to an attentional challenge. These findings provide a mechanistic and cognitive framework for interpreting the attentional phenotype associated with the human Val89 variant and establish a model that permits a more invasive interrogation of CNS effects as well as the development of therapeutic strategies for those, including Val89 carriers, with presynaptic cholinergic perturbations.

Structural Insights into the Substrate Transport Mechanisms in GTR Transporters through Ensemble Docking

Glucosinolate transporters (GTRs) are part of the nitrate/peptide transporter (NPF) family, members of which also transport specialized secondary metabolites as substrates. Glucosinolates are defense compounds derived from amino acids. We selected 4-methylthiobutyl (4MTB) and indol-3-ylmethyl (I3M) glucosinolates to study how GTR1 from Arabidopsis thaliana transports these substrates in computational simulation approaches. The designed pipeline reported here includes massive docking of 4MTB and I3M in an ensemble of GTR1 conformations (in both inward and outward conformations) extracted from molecular dynamics simulations, followed by clustered and substrate–protein interactions profiling. The identified key residues were mutated, and their role in substrate transport was tested. We were able to identify key residues that integrate a major binding site of these substrates, which is critical for transport activity. In silico approaches employed here represent a breakthrough in the plant transportomics field, as the identification of key residues usually takes a long time if performed from a purely wet-lab experimental perspective. The inclusion of structural bioinformatics in the analyses of plant transporters significantly speeds up the knowledge-gaining process and optimizes valuable time and resources.

Structures of the TMC-1 complex illuminate mechanosensory transduction

The initial step in the sensory transduction pathway underpinning hearing and balance in mammals involves the conversion of force into the gating of a mechanosensory transduction channel¹. Despite the profound socioeconomic impacts of hearing disorders and the fundamental biological significance of understanding mechanosensory transduction, the composition, structure and mechanism of the mechanosensory transduction complex have remained poorly characterized. Here we report the single-particle cryo-electron microscopy structure of the native transmembrane channel-like protein 1 (TMC-1) mechanosensory transduction complex isolated from Caenorhabditis elegans. The two-fold symmetric complex is composed of two copies each of the pore-forming TMC-1 subunit, the calcium-binding protein CALM-1 and the transmembrane inner ear protein TMIE. CALM-1 makes extensive contacts with the cytoplasmic face of the TMC-1 subunits, whereas the single-pass TMIE subunits reside on the periphery of the complex, poised like the handles of an accordion. A subset of complexes additionally includes a single arrestin-like protein, arrestin domain protein (ARRD-6), bound to a CALM-1 subunit. Single-particle reconstructions and molecular dynamics simulations show how the mechanosensory transduction complex deforms the membrane bilayer and suggest crucial roles for lipid-protein interactions in the mechanism by which mechanical force is transduced to ion channel gating.

Uniport, Not Proton-Symport, in a Non-Mammalian SLC23 Transporter

SLC23 family members are transporters of either nucleobases or ascorbate. While the mammalian SLC23 ascorbate transporters are sodium-coupled, the non-mammalian nucleobase transporters have been proposed, but not formally shown, to be proton-coupled symporters. This assignment is exclusively based on in vivo transport assays using protonophores. Here, by establishing the first in vitro transport assay for this protein family, we demonstrate that a representative member of the SLC23 nucleobase transporters operates as a uniporter instead. We explain these conflicting assignments by identifying a critical role of uracil phosphoribosyltransferase, the enzyme converting uracil to UMP, in driving uracil uptake in vivo. Detailed characterization of uracil phosphoribosyltransferase reveals that the sharp reduction of uracil uptake in whole cells in presence of protonophores is caused by acidification-induced enzyme inactivation. The SLC23 family therefore consists of both uniporters and symporters in line with the structurally related SLC4 and SLC26 families that have previously been demonstrated to accommodate both transport modes as well.

Chloride-dependent conformational changes in the GlyT1 glycine transporter

The human GlyT1 glycine transporter requires chloride for its function. However, the mechanism by which Cl^- exerts its influence is unknown. To examine the role that Cl^- plays in the transport cycle, we measured the effect of Cl^- on both glycine binding and conformational changes. The ability of glycine to displace the high-affinity radioligand [³H]CHIBA-3007 required Na⁺ and was potentiated over 1,000-fold by Cl^- We generated GlyT1b mutants containing reactive cysteine residues in either the extracellular or cytoplasmic permeation pathways and measured changes in the reactivity of those cysteine residues as indicators of conformational changes in response to ions and substrate. Na⁺ increased accessibility in the extracellular pathway and decreased it in the cytoplasmic pathway, consistent with stabilizing an outward-open conformation as observed in other members of this transporter family. In the presence of Na⁺, both glycine and Cl^- independently shifted the conformation of GlyT1b toward an outward-closed conformation. Together, Na⁺, glycine, and Cl^- stabilized an inward-open conformation of GlyT1b. We then examined whether Cl^- acts by interacting with a conserved glutamine to allow formation of an ion pair that stabilizes the closed state of the extracellular pathway. Molecular dynamics simulations of a GlyT1 homolog indicated that this ion pair is formed more frequently as that pathway closes. Mutation of the glutamine blocked the effect of Cl^-, and substituting it with glutamate or lysine resulted in outward- or inward-facing transporter conformations, respectively. These results provide an unexpected insight into the role of Cl^- in this family of transporters.

Role of FAM134 paralogues in endoplasmic reticulum remodeling, ER-phagy, and Collagen quality control

Degradation of the endoplasmic reticulum (ER) via selective autophagy (ER‐phagy) is vital for cellular homeostasis. We identify FAM134A/RETREG2 and FAM134C/RETREG3 as ER‐phagy receptors, which predominantly exist in an inactive state under basal conditions. Upon autophagy induction and ER stress signal, they can induce significant ER fragmentation and subsequent lysosomal degradation. FAM134A, FAM134B/RETREG1, and FAM134C are essential for maintaining ER morphology in a LC3‐interacting region (LIR)‐dependent manner. Overexpression of any FAM134 paralogue has the capacity to significantly augment the general ER‐phagy flux upon starvation or ER‐stress. Global proteomic analysis of FAM134 overexpressing and knockout cell lines reveals several protein clusters that are distinctly regulated by each of the FAM134 paralogues as well as a cluster of commonly regulated ER‐resident proteins. Utilizing pro‐Collagen I, as a shared ER‐phagy substrate, we observe that FAM134A acts in a LIR‐independent manner and compensates for the loss of FAM134B and FAM134C, respectively. FAM134C instead is unable to compensate for the loss of its paralogues. Taken together, our data show that FAM134 paralogues contribute to common and unique ER‐phagy pathways.

Refining pairwise sequence alignments of membrane proteins by the incorporation of anchors

The alignment of primary sequences is a fundamental step in the analysis of protein structure, function, and evolution, and in the generation of homology-based models. Integral membrane proteins pose a significant challenge for such sequence alignment approaches, because their evolutionary relationships can be very remote, and because a high content of hydrophobic amino acids reduces their complexity. Frequently, biochemical or biophysical data is available that informs the optimum alignment, for example, indicating specific positions that share common functional or structural roles. Currently, if those positions are not correctly matched by a standard pairwise sequence alignment procedure, the incorporation of such information into the alignment is typically addressed in an ad hoc manner, with manual adjustments. However, such modifications are problematic because they reduce the robustness and reproducibility of the aligned regions either side of the newly matched positions. Previous studies have introduced restraints as a means to impose the matching of positions during sequence alignments, originally in the context of genome assembly. Here we introduce position restraints, or "anchors" as a feature in our alignment tool AlignMe, providing an aid to pairwise global sequence alignment of alpha-helical membrane proteins. Applying this approach to realistic scenarios involving distantly-related and low complexity sequences, we illustrate how the addition of anchors can be used to modify alignments, while still maintaining the reproducibility and rigor of the rest of the alignment. Anchored alignments can be generated using the online version of AlignMe available at www.bioinfo.mpg.de/AlignMe/.

The organosulfur compound dimethylsulfoniopropionate (DMSP) is utilized as an osmoprotectant by Vibrio species

Dimethylsulfoniopropionate (DMSP), a key component of the global geochemical sulfur cycle, is a secondary metabolite produced in large quantities by marine phytoplankton and utilized as an osmoprotectant, thermoprotectant and antioxidant. Marine bacteria can use two pathways to degrade and catabolize DMSP, a demethylation pathway and a cleavage pathway that produces the climate active gas dimethylsulfide (DMS). Whether marine bacteria can also accumulate DMSP as an osmoprotectant to maintain the turgor pressure of the cell in response to changes in external osmolarity has received little attention. The marine halophile Vibrio parahaemolyticus, contains at least six osmolyte transporters, four betaine carnitine choline transport (BCCT) carriers BccT1-BccT4 and two ABC-family ProU transporters. In this study, we showed that DMSP is used as an osmoprotectant by V. parahaemolyticus and several other Vibrio species including V. cholerae and V. vulnificus Using a V. parahaemolyticus proU double mutant, we demonstrated that these ABC transporters are not required for DMSP uptake. However, a bccT null mutant lacking all four BCCTs had a growth defect compared to wild type in high salinity media supplemented with DMSP. Using mutants possessing only one functional BCCT in growth pattern assays, we identified two BCCT-family transporters, BccT1 and BccT2, which are carriers of DMSP. The only V. parahaemolyticus BccT homolog that V. cholerae and V. vulnificus possess is BccT3 and functional complementation in Escherichia coli MKH13 showed V. cholerae VcBccT3 could transport DMSP. In V. vulnificus strains, we identified and characterized an additional BCCT family transporter, which we named BccT5 that was also a carrier for DMSP.Importance DMSP is present in the marine environment, produced in large quantities by marine phytoplankton as an osmoprotectant, and is an important component of the global geochemical sulfur cycle. This algal osmolyte has not been previously investigated for its role in marine heterotrophic bacterial osmotic stress response. Vibrionaceae are marine species, many of which are halophiles exemplified by V. parahaemolyticus, a species that possesses at least six transporters for the uptake of osmolytes. Here, we demonstrated that V. parahaemolyticus and other Vibrio species can accumulate DMSP as an osmoprotectant and show that several BCCT family transporters uptake DMSP. These studies suggest that DMSP is a significant bacterial osmoprotectant, which may be important for understanding the fate of DMSP in the environment. DMSP is produced and present in coral mucus and Vibrio species form part of the microbial communities associated with them. The function of DMSP in these interactions is unclear, but could be an important driver for these associations allowing Vibrio proliferation. This work suggests that DMSP likely has an important role in heterotrophic bacteria ecology than previously appreciated.

Investigations of Dimethylglycine, Glycine Betaine, and Ectoine Uptake by a Betaine-Carnitine-Choline Transporter Family Transporter with Diverse Substrate Specificity in Vibrio Species

Fluctuations in osmolarity are one of the most prevalent stresses to which bacteria must adapt, both hypo- and hyperosmotic conditions. Most bacteria cope with high osmolarity by accumulating compatible solutes (osmolytes) in the cytoplasm to maintain the turgor pressure of the cell. Vibrio parahaemolyticus, a halophile, utilizes at least six compatible solute transporters for the uptake of osmolytes: two ABC family ProU transporters and four betaine-carnitine-choline transporter (BCCT) family transporters. The full range of compatible solutes transported by this species has yet to be determined. Using an osmolyte phenotypic microarray plate for growth analyses, we expanded the known osmolytes used by V. parahaemolyticus to include N,N-dimethylglycine (DMG), among others. Growth pattern analysis of four triple-bccT mutants, possessing only one functional BCCT, indicated that BccT1 (VP1456), BccT2 (VP1723), and BccT3 (VP1905) transported DMG. BccT1 was unusual in that it could take up both compounds with methylated head groups (glycine betaine [GB], choline, and DMG) and cyclic compounds (ectoine and proline). Bioinformatics analysis identified the four coordinating amino acid residues for GB in the BccT1 protein. In silico modeling analysis demonstrated that GB, DMG, and ectoine docked in the same binding pocket in BccT1. Using site-directed mutagenesis, we showed that a strain with all four residues mutated resulted in the loss of uptake of GB, DMG, and ectoine. We showed that three of the four residues were essential for ectoine uptake, whereas only one of the residues was important for GB uptake. Overall, we have demonstrated that DMG is a highly effective compatible solute for Vibrio species and have elucidated the amino acid residues in BccT1 that are important for the coordination of GB, DMG, and ectoine transport.IMPORTANCE Vibrio parahaemolyticus possesses at least six osmolyte transporters, which allow the bacterium to adapt to high-salinity conditions. In this study, we identified several additional osmolytes that were utilized by V. parahaemolyticus We demonstrated that the compound DMG, which is present in the marine environment, was a highly effective osmolyte for Vibrio species. We determined that DMG is transported via BCCT family carriers, which have not been shown previously to take up this compound. BccT1 was a carrier for GB, DMG, and ectoine, and we identified the amino acid residues essential for the coordination of these compounds. The data suggest that for BccT1, GB is more easily accommodated than ectoine in the transporter binding pocket.

Multiplexed measurement of variant abundance and activity reveals VKOR topology, active site and human variant impact

Vitamin K epoxide reductase (VKOR) drives the vitamin K cycle, activating vitamin K-dependent blood clotting factors. VKOR is also the target of the widely used anticoagulant drug, warfarin. Despite VKOR's pivotal role in coagulation, its structure and active site remain poorly understood. In addition, VKOR variants can cause vitamin K-dependent clotting factor deficiency or alter warfarin response. Here, we used multiplexed, sequencing-based assays to measure the effects of 2,695 VKOR missense variants on abundance and 697 variants on activity in cultured human cells. The large-scale functional data, along with an evolutionary coupling analysis, supports a four transmembrane domain topology, with variants in transmembrane domains exhibiting strongly deleterious effects on abundance and activity. Functionally constrained regions of the protein define the active site, and we find that, of four conserved cysteines putatively critical for function, only three are absolutely required. Finally, 25% of human VKOR missense variants show reduced abundance or activity, possibly conferring warfarin sensitivity or causing disease.

Structural Determinants of the Dopamine Transporter Regulation Mediated by G Proteins

Dopamine clearance in the brain is controlled by the dopamine transporter (DAT), a protein residing in the plasma membrane, which drives reuptake of extracellular dopamine into presynaptic neurons. Studies have revealed that the βγ subunits of heterotrimeric G proteins modulate DAT function through a physical association with the C-terminal region of the transporter. Regulation of neurotransmitter transporters by Gβγ subunits is unprecedented in the literature; therefore, it is interesting to investigate the structural details of this particular protein-protein interaction. Here, we refined the crystal structure of the Drosophila melanogaster DAT (dDAT), modeling de novo the N- and C-terminal domains; subsequently, we used the full-length dDAT structure to generate a comparative model of human DAT (hDAT). Both proteins were assembled with Gβ1γ2 subunits employing protein-protein docking, and subsequent molecular dynamics simulations were run to identify the specific interactions governing the formation of the hDAT:Gβγ and dDAT:Gβγ complexes. A [L/F]R[Q/E]R sequence motif containing the residues R588 in hDAT and R587 in dDAT was found as key to bind the Gβγ subunits through electrostatic interactions with a cluster of negatively charged residues located at the top face of the Gβ subunit. Alterations of DAT function have been associated with multiple devastating neuropathological conditions; therefore, this work represents a step toward better understanding DAT regulation by signaling proteins, allowing us to predict therapeutic target regions.

Systematic and Comparative Evaluation of Software Programs for Template-Based Modeling of Protein Structures

Modeling protein structures is critical for understanding protein functions in various biological and biotechnological studies. Among representative protein structure modeling approaches, template-based modeling (TBM) is by far the most reliable and most widely used approach to model protein structures. However, it still remains as a challenge to select appropriate software programs for pairwise alignments and model building, two major steps of the TBM. In this paper, pairwise alignment methods for TBM are first compared with respect to the quality of structure models built using these methods. This comparative study is conducted using comprehensive datasets, which cover 6185 domain sequences from Structural Classification of Proteins extended for soluble proteins, and 259 Protein Data Bank entries (whole protein sequences) from Orientations of Proteins in Membranes database for membrane proteins. Overall, a profile-based method, especially PSI-BLAST, consistently shows high performance across the datasets and model evaluation metrics used. Next, use of two model building programs, MODELLER and SWISS-MODEL, does not seem to significantly affect the quality of protein structure models built except for the Hard group (a group of relatively less homologous proteins) of membrane proteins. The results presented in this study will be useful for more accurate implementation of TBM.

Structural and evolutionary analyses of the Plasmodium falciparum chloroquine resistance transporter

Mutations in the Plasmodium falciparum chloroquine resistance transporter (PfCRT) confer resistance to several antimalarial drugs such as chloroquine (CQ) or piperaquine (PPQ), a partner molecule in current artemisinin-based combination therapies. As a member of the Drug/Metabolite Transporter (DMT) superfamily, the vacuolar transporter PfCRT may translocate substrate molecule(s) across the membrane of the digestive vacuole (DV), a lysosome-like organelle. However, the physiological substrate(s), the transport mechanism and the functional regions of PfCRT remain to be fully characterized. Here, we hypothesized that identification of evolutionary conserved sites in a tertiary structural context could help locate putative functional regions of PfCRT. Hence, site-specific substitution rates were estimated over Plasmodium evolution at each amino acid sites, and the PfCRT tertiary structure was predicted in both inward-facing (open-to-vacuole) and occluded states through homology modeling using DMT template structures sharing <15% sequence identity with PfCRT. We found that the vacuolar-half and membrane-spanning domain (and especially the transmembrane helix 9) of PfCRT were more conserved, supporting that its physiological substrate is expelled out of the parasite DV. In the PfCRT occluded state, some evolutionary conserved sites, including positions related to drug resistance mutations, participate in a putative binding pocket located at the core of the PfCRT membrane-spanning domain. Through structural comparison with experimentally-characterized DMT transporters, we identified several conserved PfCRT amino acid sites located in this pocket as robust candidates for mediating substrate transport. Finally, in silico mutagenesis revealed that drug resistance mutations caused drastic changes in the electrostatic potential of the transporter vacuolar entry and pocket, facilitating the escape of protonated CQ and PPQ from the parasite DV.

Intestinal serotonin and fluoxetine exposure modulate bacterial colonization in the gut

The gut microbiota regulates levels of serotonin (5-hydroxytryptamine (5-HT)) in the intestinal epithelium and lumen^1-5. However, whether 5-HT plays a functional role in bacteria from the gut microbiota remains unknown. We demonstrate that elevating levels of intestinal lumenal 5-HT by oral supplementation or genetic deficiency in the host 5-HT transporter (SERT) increases the relative abundance of spore-forming members of the gut microbiota, which were previously reported to promote host 5-HT biosynthesis. Within this microbial community, we identify Turicibacter sanguinis as a gut bacterium that expresses a neurotransmitter sodium symporter-related protein with sequence and structural homology to mammalian SERT. T. sanguinis imports 5-HT through a mechanism that is inhibited by the selective 5-HT reuptake inhibitor fluoxetine. 5-HT reduces the expression of sporulation factors and membrane transporters in T. sanguinis, which is reversed by fluoxetine exposure. Treating T. sanguinis with 5-HT or fluoxetine modulates its competitive colonization in the gastrointestinal tract of antibiotic-treated mice. In addition, fluoxetine reduces the membership of T. sanguinis in the gut microbiota of conventionally colonized mice. Host association with T. sanguinis alters intestinal expression of multiple gene pathways, including those important for lipid and steroid metabolism, with corresponding reductions in host systemic triglyceride levels and inguinal adipocyte size. Together, these findings support the notion that select bacteria indigenous to the gut microbiota signal bidirectionally with the host serotonergic system to promote their fitness in the intestine.

A Chemosensory GPCR as a Potential Target to Control the Root-Knot Nematode Meloidogyne incognita Parasitism in Plants

Root-knot nematodes (RKN), from the Meloidogyne genus, have a worldwide distribution and cause severe economic damage to many life-sustaining crops. Because of their lack of specificity and danger to the environment, most chemical nematicides have been banned from use. Thus, there is a great need for new and safe compounds to control RKN. Such research involves identifying beforehand the nematode proteins essential to the invasion. Since G protein-coupled receptors GPCRs are the target of a large number of drugs, we have focused our research on the identification of putative nematode GPCRs such as those capable of controlling the movement of the parasite towards (or within) its host. A datamining procedure applied to the genome of Meloidogyne incognita allowed us to identify a GPCR, belonging to the neuropeptide GPCR family that can serve as a target to carry out a virtual screening campaign. We reconstructed a 3D model of this receptor by homology modeling and validated it through extensive molecular dynamics simulations. This model was used for large scale molecular dockings which produced a filtered limited set of putative antagonists for this GPCR. Preliminary experiments using these selected molecules allowed the identification of an active compound, namely C260-2124, from the ChemDiv provider, which can serve as a starting point for further investigations.

Topology prediction of insect olfactory receptors

Olfactory receptors are important transmembrane proteins that enable organisms to perceive odours and react to them. Structural understanding of insect olfactory receptors is scarce. In this review, we discuss different transmembrane helix prediction methods, consensus methods, topology prediction methods which can enable topology prediction of these proteins. We discuss the current success rates by applying the algorithms on few G-protein coupled receptors of known structure and olfactory receptor sequences and outstanding challenges. Finally, we discuss the impact of topology prediction on biology and modeling of ORs.

A Structural Model of the Inactivation Gate of Voltage-Activated Potassium Channels

After opening, the Shaker voltage-gated potassium (K_V) channel rapidly inactivates when one of its four N-termini enters and occludes the channel pore. Although it is known that the tip of the N-terminus reaches deep into the central cavity, the conformation adopted by this domain during inactivation and the nature of its interactions with the rest of the channel remain unclear. Here, we use molecular dynamics simulations coupled with electrophysiology experiments to reveal the atomic-scale mechanisms of inactivation. We find that the first six amino acids of the N-terminus spontaneously enter the central cavity in an extended conformation, establishing hydrophobic contacts with residues lining the pore. A second portion of the N-terminus, consisting of a long 24 amino acid α-helix, forms numerous polar contacts with residues in the intracellular entryway of the T1 domain. Double mutant cycle analysis revealed a strong relationship between predicted interatomic distances and empirically observed thermodynamic coupling, establishing a plausible model of the transition of K_V channels to the inactivated state.

Curvature induction and membrane remodeling by FAM134B reticulon homology domain assist selective ER-phagy

FAM134B/RETREG1 is a selective ER-phagy receptor that regulates the size and shape of the endoplasmic reticulum. The structure of its reticulon-homology domain (RHD), an element shared with other ER-shaping proteins, and the mechanism of membrane shaping remain poorly understood. Using molecular modeling and molecular dynamics (MD) simulations, we assemble a structural model for the RHD of FAM134B. Through MD simulations of FAM134B in flat and curved membranes, we relate the dynamic RHD structure with its two wedge-shaped transmembrane helical hairpins and two amphipathic helices to FAM134B functions in membrane-curvature induction and curvature-mediated protein sorting. FAM134B clustering, as expected to occur in autophagic puncta, amplifies the membrane-shaping effects. Electron microscopy of in vitro liposome remodeling experiments support the membrane remodeling functions of the different RHD structural elements. Disruption of the RHD structure affects selective autophagy flux and leads to disease states.

FAM134B/RETREG1 is a selective ER-phagy receptor that regulates the size and shape of the endoplasmic reticulum. Here authors use molecular modeling and molecular dynamics simulations to assemble a structural model for the reticulon-homology domain of FAM134B.

Entry Exclusion of Conjugative Plasmids of the IncA, IncC, and Related Untyped Incompatibility Groups

Conjugative plasmids of incompatibility group C (IncC), formerly known as A/C₂, disseminate antibiotic resistance genes globally in diverse pathogenic species of Gammaproteobacteria. Salmonella genomic island 1 (SGI1) can be mobilized by IncC plasmids and was recently shown to reshape the conjugative type IV secretion system (T4SS) encoded by these plasmids to evade entry exclusion. Entry exclusion blocks DNA translocation between cells containing identical or highly similar plasmids. Here, we report that the protein encoded by the entry exclusion gene of IncC plasmids (eexC) mediates entry exclusion in recipient cells through recognition of the IncC-encoded TraG_C protein in donor cells. Phylogenetic analyses based on EexC and TraG_C homologs predicted the existence of at least three different exclusion groups among IncC-related conjugative plasmids. Mating assays using Eex proteins encoded by representative IncC and IncA (former A/C₁) and related untyped plasmids confirmed these predictions and showed that the IncC and IncA plasmids belong to the C exclusion group, thereby explaining their apparent incompatibility despite their compatible replicons. Representatives of the two other exclusion groups (D and E) are untyped conjugative plasmids found in Aeromonas sp. Finally, we determined through domain swapping that the carboxyl terminus of the EexC and EexE proteins controls the specificity of these exclusion groups. Together, these results unravel the role of entry exclusion in the apparent incompatibility between IncA and IncC plasmids while shedding light on the importance of the TraG subunit substitution used by SGI1 to evade entry exclusion.IMPORTANCE IncA and IncC conjugative plasmids drive antibiotic resistance dissemination among several pathogenic species of Gammaproteobacteria due to the diversity of drug resistance genes that they carry and their ability to mobilize antibiotic resistance-conferring genomic islands such as SGI1 of Salmonella enterica While historically grouped as "IncA/C," IncA and IncC replicons were recently confirmed to be compatible and to abolish each other's entry into the cell in which they reside during conjugative transfer. The significance of our study is in identifying an entry exclusion system that is shared by IncA and IncC plasmids. It impedes DNA transfer to recipient cells bearing a plasmid of either incompatibility group. The entry exclusion protein of this system is unrelated to any other known entry exclusion proteins.

Structure of Outward-Facing PglK and Molecular Dynamics of Lipid-Linked Oligosaccharide Recognition and Translocation

PglK is a lipid-linked oligosaccharide (LLO) flippase essential for asparagine-linked protein glycosylation in Campylobacter jejuni. Previously we have proposed a non-alternating-access LLO translocation mechanism, where postulated outward-facing states play a primary role. To investigate this unusual mechanistic proposal, we have determined a high-resolution structure of PglK that displays an outward semi-occluded state with the two nucleotide binding domains forming an asymmetric closed dimer with two bound ATPγS molecules. Based on this structure, we performed extensive molecular dynamics simulations to investigate LLO recognition and flipping. Our results suggest that PglK may employ a "substrate-hunting" mechanism to locally increase the LLO concentration and facilitate its jump into the translocation pathway, for which sugars from the LLO head group are essential. We further conclude that the release of LLO to the outside occurs before ATP hydrolysis and is followed by the closing of the periplasmic cavity of PglK.

Unprecedented inhibition of P-gp activity by a novel ruthenium-cyclopentadienyl compound bearing a bipyridine-biotin ligand

Two new ruthenium complexes, [Ru(η⁵-Cp)(PPh₃)(2,2'-bipy-4,4'-R)]⁺ with R = -CH₂OH (Ru1) or dibiotin ester (Ru2) were synthesized and fully characterized. Both compounds were tested against two types of breast cancer cells (MCF7 and MDA-MB-231), showing better cytotoxicity than cisplatin in the same experimental conditions. Since multidrug resistance (MDR) is one of the main problems in cancer chemotherapy, we have assessed the potential of these compounds to overcome resistance to treatments. Ru2 showed exceptional selectivity as P-gp inhibitor, while Ru1 is possibly a substrate. In vivo studies in zebrafish showed that Ru2 is well tolerated up to 1.17 mg/L, presenting a LC₅₀ of 5.73 mg/L at 5 days post fertilization.

Molecular Modeling of Drug-Transporter Interactions-An International Transporter Consortium Perspective

Membrane transporters play diverse roles in the pharmacokinetics and pharmacodynamics of small-molecule drugs. Understanding the mechanisms of drug-transporter interactions at the molecular level is, therefore, essential for the design of drugs with optimal therapeutic effects. This white paper examines recent progress, applications, and challenges of molecular modeling of membrane transporters, including modeling techniques that are centered on the structures of transporter ligands, and those focusing on the structures of the transporters. The goals of this article are to illustrate current best practices and future opportunities in using molecular modeling techniques to understand and predict transporter-mediated effects on drug disposition and efficacy.Membrane transporters from the solute carrier (SLC) and ATP-binding cassette (ABC) superfamilies regulate the cellular uptake, efflux, and homeostasis of many essential nutrients and significantly impact the pharmacokinetics of drugs; further, they may provide targets for novel therapeutics as well as facilitate prodrug approaches. Because of their often broad substrate selectivity they are also implicated in many undesirable and sometimes life-threatening drug-drug interactions (DDIs).^5,6.

A dual substrate-accessing mechanism of a major facilitator superfamily protein facilitates lysophospholipid flipping across the cell membrane

Lysophospholipid transporter (LplT) is a member of the major facilitator superfamily present in many Gram-negative bacteria. LplT catalyzes flipping of lysophospholipids (LPLs) across the bacterial inner membrane, playing an important role in bacterial membrane homeostasis. We previously reported that LplT promotes both uptake of exogenous LPLs and intramembranous LPL flipping across the bilayer. To gain mechanistic insight into this dual LPL-flipping activity, here we implemented a combination of computational approaches and LPL transport analyses to study LPL binding of and translocation by LplT. Our results suggest that LplT translocates LPLs through an elongated cavity exhibiting an extremely asymmetric polarity. We found that two D(E)N motifs form a head group-binding site, in which the carboxylate group of Asp-30 is important for LPL head group recognition. Substitutions of residues in the head group-binding site disrupted both LPL uptake and flipping activities. However, alteration of hydrophobic residues on the interface between the N- and C-terminal domains impaired LPL flipping specifically, resulting in LPLs accumulation in the membrane, but LPL uptake remained active. These results suggest a dual substrate-accessing mechanism, in which LplT recruits LPLs to its substrate-binding site via two routes, either from its extracellular entry or through a membrane-embedded groove between transmembrane helices, and then moves them toward the inner membrane leaflet. This LPL-flipping mechanism is likely conserved in many bacterial species, and our findings illustrate how LplT adjusts the major facilitator superfamily translocation pathway to perform its versatile lipid homeostatic functions.

A Comparative Study of Modern Homology Modeling Algorithms for Rhodopsin Structure Prediction

Rhodopsins are seven α-helical membrane proteins that are of great importance in chemistry, biology, and modern biotechnology. Any in silico study on rhodopsin properties and functioning requires a high-quality three-dimensional structure. Due to particular difficulties with obtaining membrane protein structures from the experiment, in silico prediction of the three-dimensional rhodopsin structure based only on its primary sequence is an especially important task. For the last few years, significant progress was made in the field of protein structure prediction, especially for methods based on comparative modeling. However, the majority of this progress was made for soluble proteins and further investigations are needed to achieve similar progress for membrane proteins. In this paper, we evaluate the performance of modern protein structure prediction methodologies (implemented in the Medeller, I-TASSER, and Rosetta packages) for their ability to predict rhodopsin structures. Three widely used methodologies were considered: two general methodologies that are commonly applied to soluble proteins and a methodology that uses constraints that are specific for membrane proteins. The test pool consisted of 36 target-template pairs with different sequence similarities that was constructed on the basis of 24 experimental rhodopsin structures taken from the RCSB database. As a result, we showed that all three considered methodologies allow obtaining rhodopsin structures with the quality that is close to the crystallographic one (root mean square deviation (RMSD) of the predicted structure from the corresponding X-ray structure up to 1.5 Å) if the target-template sequence identity is higher than 40%. Moreover, all considered methodologies provided structures of average quality (RMSD < 4.0 Å) if the target-template sequence identity is higher than 20%. Such structures can be subsequently used for further investigation of molecular mechanisms of protein functioning and for the development of modern protein-based biotechnologies.

Structural relationship between the putative hair cell mechanotransduction channel TMC1 and TMEM16 proteins

The hair cell mechanotransduction (MET) channel complex is essential for hearing, yet it's molecular identity and structure remain elusive. The transmembrane channel-like 1 (TMC1) protein localizes to the site of the MET channel, interacts with the tip-link responsible for mechanical gating, and genetic alterations in TMC1 alter MET channel properties and cause deafness, supporting the hypothesis that TMC1 forms the MET channel. We generated a model of TMC1 based on X-ray and cryo-EM structures of TMEM16 proteins, revealing the presence of a large cavity near the protein-lipid interface that also harbors the Beethoven mutation, suggesting that it could function as a permeation pathway. We also find that hair cells are permeable to 3 kDa dextrans, and that dextran permeation requires TMC1/2 proteins and functional MET channels, supporting the presence of a large permeation pathway and the hypothesis that TMC1 is a pore forming subunit of the MET channel complex.

A CLC-type F-/H+ antiporter in ion-swapped conformations

Fluoride/proton antiporters of the CLCF family combat F- toxicity in bacteria by exporting this halide from the cytoplasm. These transporters belong to the widespread CLC superfamily but display transport properties different from those of the well-studied Cl-/H+ antiporters. Here, we report a structural and functional investigation of these F--transport proteins. Crystal structures of a CLCF homolog from Enterococcus casseliflavus are captured in two conformations with simultaneous accessibility of F- and H+ ions via separate pathways on opposite sides of the membrane. Manipulation of a key glutamate residue critical for H+ and F- transport reverses the anion selectivity of transport; replacement of the glutamate with glutamine or alanine completely inhibits F- and H+ transport while allowing for rapid uncoupled flux of Cl-. The structural and functional results lead to a 'windmill' model of CLC antiport wherein F- and H+ simultaneously move through separate ion-specific pathways that switch sidedness during the transport cycle.

A direct interaction of cholesterol with the dopamine transporter prevents its out-to-inward transition

Monoamine transporters (MATs) carry out neurotransmitter reuptake from the synaptic cleft, a key step in neurotransmission, which is targeted in the treatment of neurological disorders. Cholesterol (CHOL), a major component of the synaptic plasma membrane, has been shown to exhibit a modulatory effect on MATs. Recent crystal structures of the dopamine transporter (DAT) revealed the presence of two conserved CHOL-like molecules, suggesting a functional protein-CHOL direct interaction. Here, we present extensive atomistic molecular dynamics (MD) simulations of DAT in an outward-facing conformation. In the absence of bound CHOL, DAT undergoes structural changes reflecting early events of dopamine transport: transition to an inward-facing conformation. In contrast, in the presence of bound CHOL, these conformational changes are inhibited, seemingly by an immobilization of the intracellular interface of transmembrane helix 1a and 5 by CHOL. We also provide evidence, from coarse grain MD simulations that the CHOL sites observed in the DAT crystal structures are preserved in all human monoamine transporters (dopamine, serotonin and norepinephrine), suggesting that our findings might extend to the entire family.

Iodide Binding in Sodium-Coupled Cotransporters

Several apical iodide translocation pathways have been proposed for iodide efflux out of thyroid follicular cells, including a pathway mediated by the sodium-coupled monocarboxylate transporter 1 (SMCT1), which remains controversial. Herein, we evaluate structural and functional similarities between SMCT1 and the well-studied sodium-iodide symporter (NIS) that mediates the first step of iodide entry into the thyroid. Free-energy calculations using a force field with electronic polarizability verify the presence of a conserved iodide-binding pocket between the TM2, TM3, and TM7 segments in hNIS, where iodide is coordinated by Phe67, Gln72, Cys91, and Gln94. We demonstrate the mutation of residue Gly93 of hNIS to a larger amino acid expels the side chain of a critical tryptophan residue (Trp255) into the interior of the binding pocket, partially occluding the iodide binding site and reducing iodide affinity, which is consistent with previous reports associating mutation of this residue with iodide uptake deficiency and hypothyroidism. Furthermore, we find that the position of Trp255 in this hNIS mutant mirrors that of Trp253 in wild-type hSMCT1, where a threonine (Thr91) occupies the position homologous to that occupied by glycine in wild-type hNIS (Gly93). Correspondingly, mutation of Thr91 to glycine in hSMCT1 makes the pocket structure more like that of wild-type hNIS, increasing its iodide affinity. These results suggest that wild-type hSMCT1 in the inward-facing conformation may bind iodide only very weakly, which may have implications for its ability to transport iodide.

Structural basis of nucleotide sugar transport across the Golgi membrane

Glycosylation is a fundamental cellular process that in eukaryotes occurs in the lumen of both the Golgi apparatus and endoplasmic reticulum 1. Nucleotide sugar transporters (NSTs) are an essential component of the glyscosylation pathway, providing the diverse range of substrates required for the glycosyltransferases 2,3. NSTs are linked to several developmental and immune disorders in humans and in pathogenic microbes play an important role in virulence 4–8. How NSTs recognise and transport activated monosaccharides however is currently unclear. Here we present the first crystal structure of an NST, the GDP-mannose transporter Vrg4, in both the substrate free and bound states. A hitherto unobserved requirement for short chain lipids in activating the transporter supports a model for regulation within the highly dynamic membranes of the Golgi apparatus. Our results provide a structural basis for understanding nucleotide sugar recognition and provide insights into the transport and regulatory mechanism for this family of intracellular transporters.

Molecular Recognition of Agonists and Antagonists by the Nucleotide-Activated G Protein-Coupled P2Y2 Receptor

A homology model of the nucleotide-activated P2Y₂R was created based on the X-ray structures of the P2Y₁ receptor. Docking studies were performed, and receptor mutants were created to probe the identified binding interactions. Mutation of residues predicted to interact with the ribose (Arg110) and the phosphates of the nucleotide agonists (Arg265, Arg292) or that contribute indirectly to binding (Tyr288) abolished activity. The Y114F, R194A, and F261A mutations led to inactivity of diadenosine tetraphosphate and to a reduced response of UTP. Significant reduction in agonist potency was observed for all other receptor mutants (Phe111, His184, Ser193, Phe261, Tyr268, Tyr269) predicted to be involved in agonist recognition. An ionic lock between Asp185 and Arg292 that is probably involved in receptor activation interacts with the phosphate groups. The antagonist AR-C118925 and anthraquinones likely bind to the orthosteric site. The updated homology models will be useful for virtual screening and drug design.

Refined topology model of the STT3/Stt3 protein subunit of the oligosaccharyltransferase complex

The oligosaccharyltransferase complex, localized in the endoplasmic reticulum (ER) of eukaryotic cells, is responsible for the N-linked glycosylation of numerous protein substrates. The membrane protein STT3 is a highly conserved part of the oligosaccharyltransferase and likely contains the active site of the complex. However, understanding the catalytic determinants of this system has been challenging, in part because of a discrepancy in the structural topology of the bacterial versus eukaryotic proteins and incomplete information about the mechanism of membrane integration. Here, we use a glycosylation mapping approach to investigate these questions. We measured the membrane integration efficiency of the mouse STT3-A and yeast Stt3p transmembrane domains (TMDs) and report a refined topology of the N-terminal half of the mouse STT3-A. Our results show that most of the STT3 TMDs are well inserted into the ER membrane on their own or in the presence of the natural flanking residues. However, for the mouse STT3-A hydrophobic domains 4 and 6 and yeast Stt3p domains 2, 3a, 3c, and 6 we measured reduced insertion efficiency into the ER membrane. Furthermore, we mapped the first half of the STT3-A protein, finding two extra hydrophobic domains between the third and the fourth TMD. This result indicates that the eukaryotic STT3 has 13 transmembrane domains, consistent with the structure of the bacterial homolog of STT3 and setting the stage for future combined efforts to interrogate this fascinating system.

Helical jackknives control the gates of the double-pore K+ uptake system KtrAB

Ion channel gating is essential for cellular homeostasis and is tightly controlled. In some eukaryotic and most bacterial ligand-gated K⁺ channels, RCK domains regulate ion fluxes. Until now, a single regulatory mechanism has been proposed for all RCK-regulated channels, involving signal transduction from the RCK domain to the gating area. Here, we present an inactive ADP-bound structure of KtrAB from Vibrio alginolyticus, determined by cryo-electron microscopy, which, combined with EPR spectroscopy and molecular dynamics simulations, uncovers a novel regulatory mechanism for ligand-induced action at a distance. Exchange of activating ATP to inactivating ADP triggers short helical segments in the K⁺-translocating KtrB dimer to organize into two long helices that penetrate deeply into the regulatory RCK domains, thus connecting nucleotide-binding sites and ion gates. As KtrAB and its homolog TrkAH have been implicated as bacterial pathogenicity factors, the discovery of this functionally relevant inactive conformation may advance structure-guided drug development.

DOI:http://dx.doi.org/10.7554/eLife.24303.001

Development of Potent and Selective Antagonists for the UTP-Activated P2Y4 Receptor

P2Y₄ is a G_q protein-coupled receptor activated by uridine-5'-triphosphate (UTP), which is widely expressed in the body, e.g., in intestine, heart, and brain. No selective P2Y₄ receptor antagonist has been described so far. Therefore, we developed and optimized P2Y₄ receptor antagonists based on an anthraquinone scaffold. Potency was assessed by a fluorescence-based assay measuring inhibition of UTP-induced intracellular calcium release in 1321N1 astrocytoma cells stably transfected with the human P2Y₄ receptor. The most potent compound of the present series, sodium 1-amino-4-[4-(2,4-dimethylphenylthio)phenylamino]-9,10-dioxo-9,10-dihydroanthracene-2-sulfonate (PSB-16133, 61) exhibited an IC₅₀ value of 233 nM, selectivity versus other P2Y receptor subtypes, and is thought to act as an allosteric antagonist. A receptor homology model was built and docking studies were performed to analyze ligand-receptor interactions. Compound 64 (PSB-1699, sodium 1-amino-4-[4-(3-pyridin-3-ylmethylthio)phenylamino]-9,10-dioxo-9,10-dihydroanthracene-2-sulfonate) represents the most selective P2Y₄ receptor antagonist known to date. Compounds 61 and 64 are therefore anticipated to become useful tools for studying this scarcely investigated receptor.

An Extra Amino Acid Residue in Transmembrane Domain 10 of the γ-Aminobutyric Acid (GABA) Transporter GAT-1 Is Required for Efficient Ion-coupled Transport

The GABA transporter GAT-1 mediates electrogenic transport of its substrate together with sodium and chloride. It is a member of the neurotransmitter:sodium:symporters, which are crucial for synaptic transmission. Compared with all other neurotransmitter:sodium:symporters, GAT-1 and other members of the GABA transporter subfamily all contain an extra amino acid residue at or near a conserved glycine in transmembrane segment 10. Therefore, we studied the functional impact of deletion and replacement mutants of Gly-457 and its two adjacent residues in GAT-1. The glycine replacement mutants were devoid of transport activity, but remarkably the deletion mutant was active, as were mutants obtained by deleting positions on either side of Gly-457. However, the inward rectification of GABA-induced transport currents by all three deletion mutants was diminished, and the charge-to-flux ratio was increased by more than 2.5-fold, both of which indicate substantial uncoupled transport. These observations suggest that the deletions render the transporters less tightly packed. Consistent with this interpretation, the inactive G457A mutant was partially rescued by removing the adjacent serine residue. Moreover, the activity of several gating mutants was also partially rescued upon deletion of Gly-457. Structural modeling showed that the stretch surrounding Gly-457 is likely to form a π-helix. Our data indicate that the "extra" residue in transmembrane domain 10 of the GABA transporter GAT-1 provides extra bulk, probably in the form of a π-helix, which is required for stringent gating and tight coupling of ion and substrate fluxes in the GABA transporter family.

Computational Approaches for Revealing the Structure of Membrane Transporters: Case Study on Bilitranslocase

The structural and functional details of transmembrane proteins are vastly underexplored, mostly due to experimental difficulties regarding their solubility and stability. Currently, the majority of transmembrane protein structures are still unknown and this present a huge experimental and computational challenge. Nowadays, thanks to X-ray crystallography or NMR spectroscopy over 3000 structures of membrane proteins have been solved, among them only a few hundred unique ones. Due to the vast biological and pharmaceutical interest in the elucidation of the structure and the functional mechanisms of transmembrane proteins, several computational methods have been developed to overcome the experimental gap. If combined with experimental data the computational information enables rapid, low cost and successful predictions of the molecular structure of unsolved proteins. The reliability of the predictions depends on the availability and accuracy of experimental data associated with structural information. In this review, the following methods are proposed for in silico structure elucidation: sequence-dependent predictions of transmembrane regions, predictions of transmembrane helix–helix interactions, helix arrangements in membrane models, and testing their stability with molecular dynamics simulations. We also demonstrate the usage of the computational methods listed above by proposing a model for the molecular structure of the transmembrane protein bilitranslocase. Bilitranslocase is bilirubin membrane transporter, which shares similar tissue distribution and functional properties with some of the members of the Organic Anion Transporter family and is the only member classified in the Bilirubin Transporter Family. Regarding its unique properties, bilitranslocase is a potentially interesting drug target.

A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling

It was demonstrated in recent studies that some rhodopsins can be used in optogenetics as fluorescent indicators of membrane voltage. One of the promising candidates for these applications is archaerhodopsin-3. While it has already shown encouraging results, there is still a large room for improvement. One of possible directions is increasing the intensity of the protein's fluorescent signal. Rational design of mutants with an improved signal is an important task, which requires both experimental and theoretical studies. Herein, we used a homology-based computational approach to predict the three-dimensional structure of archaerhodopsin-3, and a Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid approach with high-level multireference ab initio methodology (SORCI+Q/AMBER) to model optical properties of this protein. We demonstrated that this methodology allows for reliable prediction of structure and spectral properties of archaerhodopsin-3. The results of this study can be utilized for computational molecular design of efficient fluorescent indicators of membrane voltage for modern optogenetics on the basis of archaerhodopsin-3.

Transmembrane region prediction by using sequence-derived features and machine learning methods

Membrane proteins are central to carrying out impressive biological functions.

Emulating proton-induced conformational changes in the vesicular monoamine transporter VMAT2 by mutagenesis

Neurotransporters located in synaptic vesicles are essential for communication between nerve cells in a process mediated by neurotransmitters. Vesicular monoamine transporter (VMAT), a member of the largest superfamily of transporters, mediates transport of monoamines to synaptic vesicles and storage organelles in a process that involves exchange of two H⁺ per substrate. VMAT transport is inhibited by the competitive inhibitor reserpine, a second-line agent to treat hypertension, and by the noncompetitive inhibitor tetrabenazine, presently in use for symptomatic treatment of hyperkinetic disorders. During the transport cycle, VMAT is expected to occupy at least three different conformations: cytoplasm-facing, occluded, and lumen-facing. The lumen- to cytoplasm-facing transition, facilitated by protonation of at least one of the essential membrane-embedded carboxyls, generates a binding site for reserpine. Here we have identified residues in the cytoplasmic gate and show that mutations that disrupt the interactions in this gate also shift the equilibrium toward the cytoplasm-facing conformation, emulating the effect of protonation. These experiments provide significant insight into the role of proton translocation in the conformational dynamics of a mammalian H⁺-coupled antiporter, and also identify key aspects of the mode of action and binding of two potent inhibitors of VMAT2: reserpine binds the cytoplasm-facing conformation, and tetrabenazine binds the lumen-facing conformation.

SLC Transporters: Structure, Function, and Drug Discovery

The human Solute Carrier (SLC) transporters are important targets for drug development. Structure-based drug discovery for SLC transporters requires the description of their structure, dynamics, and mechanism of interaction with small molecule ligands and ions. The recent determination of atomic structures of human SLC transporters and their homologs, combined with improved computational power and prediction methods have led to an increased applicability of structure-based drug design methods for human SLC members. In this review, we provide an overview of the SLC transporters' structures and transport mechanisms. We then describe computational techniques, such as homology modeling and virtual screening that are emerging as key tools to discover chemical probes for human SLC members. We illustrate the utility of these methods by presenting case studies in which rational integration of computation and experiment was used to characterize SLC members that transport key nutrients and metabolites, including the amino acid transporters LAT-1 and ASCT2, the SLC13 family of citric acid cycle intermediate transporters, and the glucose transporter GLUT1. We conclude with a brief discussion about future directions in structure-based drug discovery for the human SLC superfamily, one of the most structurally and functionally diverse protein families in human.

Non-Amyloid-β Component of Human α-Synuclein Oligomers Induces Formation of New Aβ Oligomers: Insight into the Mechanisms That Link Parkinson's and Alzheimer's Diseases

Parkinson's disease (PD) is characterized by the formation of Lewy bodies (LBs), of which their major component is the non-amyloid-β component (NAC) of α-synuclein (AS). Clinical studies have identified a link between PD and Alzheimer's disease (AD), but the question of why PD patients are at risk to develop various types of dementia, such as AD, is still elusive. In vivo studies have shown that Aβ can act as a seed for NAC/AS aggregation, promoting NAC/AS aggregation and thus contributing to the etiology of PD. However, the mechanisms by which NAC/AS oligomers interact with Aβ oligomers are still elusive. This work presents the interactions between NAC oligomers and Aβ oligomers at atomic resolution by applying extensive molecular dynamics simulations for an ensemble of cross-seeded NAC-Aβ(1-42) oligomers. The main conclusions of this study are as follows: first, the cross-seeded NAC-Aβ(1-42) oligomers represent polymorphic states, yet NAC oligomers prefer to interact with Aβ(1-42) oligomers to form double-layer over single-layer conformations due to electrostatic/hydrophobic interactions; second, among the single-layer conformations, the NAC oligomers induce formation of new β-strands in Aβ(1-42) oligomers, thus leading to new Aβ oligomer structures; and third, NAC oligomers stabilize the cross-β structure of Aβ oligomers, i.e., yielding compact Aβ fibril-like structures.

Binding of the multimodal antidepressant drug vortioxetine to the human serotonin transporter

Selective inhibitors of the human serotonin transporter (hSERT) have been first-line treatment against depression for several decades. Recently, vortioxetine was approved as a new therapeutic option for the treatment of depression. Vortioxetine represents a new class of antidepressant drugs with a multimodal pharmacological profile that in addition to potent inhibition of hSERT include agonistic or antagonistic effects at different serotonin receptors. We used a combination of computational, chemical, and biological methods to decipher the molecular basis for high affinity binding of vortioxetine in hSERT. X-ray crystal structures of the bacterial amino acid transporter LeuT and the Drosophila melanogaster dopamine transporter were used to build homology models of hSERT. Comparative modeling and ligand docking suggest that vortioxetine can adopt several distinct binding modes within the central binding site of hSERT. To distinguish between the identified binding modes, we determined the effect of 57 functional hSERT point mutants on vortioxetine potency and characterized seven structurally related analogs of vortioxetine in a subset of the point mutants. This allowed us to determine the orientation of vortioxetine within the central binding site and showed that only one of the proposed binding modes is functionally relevant. The findings provide important new insight about the molecular basis for high affinity recognition of vortioxetine in hSERT, which is essential for future structure-based drug discovery of novel multimodal drugs with fine-tuned selectivity across different transporter and receptor proteins in the human brain.

Computational approaches to study the effects of small genomic variations

Advances in DNA sequencing technologies have led to an avalanche-like increase in the number of gene sequences deposited in public databases over the last decade as well as the detection of an enormous number of previously unseen nucleotide variants therein. Given the size and complex nature of the genome-wide sequence variation data, as well as the rate of data generation, experimental characterization of the disease association of each of these variations or their effects on protein structure/function would be costly, laborious, time-consuming, and essentially impossible. Thus, in silico methods to predict the functional effects of sequence variations are constantly being developed. In this review, we summarize the major computational approaches and tools that are aimed at the prediction of the functional effect of mutations, and describe the state-of-the-art databases that can be used to obtain information about mutation significance. We also discuss future directions in this highly competitive field.

A Proposed Atomic Structure of the Self-Assembly of the Non-Amyloid-β Component of Human α-Synuclein As Derived by Computational Tools

α-Synuclein (AS) fibrils are the major hallmarks of Parkinson's disease (PD). It is known that the central domain of the 140-residue AS protein, known as the non-amyloid-β component (NAC), plays a crucial role in aggregation. The secondary structure of AS fibrils (including the NAC domain) has been proposed on the basis of solid-state nuclear magnetic resonance studies, but the atomic structure of the self-assembly of NAC (or AS itself) is still elusive. This is the first study that presents a detailed three-dimensional structure of NAC at atomic resolution. The proposed self-assembled structure of NAC consists of three β-strands connected by two turn regions. Our study shows that calculated structural parameter values of the simulated fibril-like cross-β structure of NAC are in excellent agreement with the experimental values. Moreover, the diameter dimensions of the proposed fibril-like structure are also in agreement with experimental measurements. The proposed fibril-like structure of NAC may assist in future work aimed at understanding the formation of aggregates in PD and developing compounds to modulate aggregation.