Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

RNA polymerase SI3 domain modulates global transcriptional pausing and pause-site fluctuations

Transcriptional pausing aids gene regulation by cellular RNA polymerases (RNAPs). A surface-exposed domain inserted into the catalytic trigger loop (TL) of Escherichia coli RNAP, called SI3, modulates pausing and is essential for growth. Here we describe a viable E. coli strain lacking SI3 enabled by a suppressor TL substitution (β'Ala941→Thr; ΔSI3*). ΔSI3* increased transcription rate in vitro relative to ΔSI3, possibly explaining its viability, but retained both positive and negative effects of ΔSI3 on pausing. ΔSI3* inhibited pauses stabilized by nascent RNA structures (pause hairpins; PHs) but enhanced other pauses. Using NET-seq, we found that ΔSI3*-enhanced pauses resemble the consensus elemental pause sequence whereas sequences at ΔSI3*-suppressed pauses, which exhibited greater association with PHs, were more divergent. ΔSI3*-suppressed pauses also were associated with apparent pausing one nucleotide upstream from the consensus sequence, often generating tandem pause sites. These '-2 pauses' were stimulated by pyrophosphate in vitro and by addition of apyrase to degrade residual NTPs during NET-seq sample processing. We propose that some pauses are readily reversible by pyrophosphorolysis or single-nucleotide cleavage. Our results document multiple ways that SI3 modulates pausing in vivo and may explain discrepancies in consensus pause sequences in some NET-seq studies.

Unraveling Activation-Related Rearrangements and Intrinsic Divergence from Ligand-Specific Conformational Changes of the Dopamine D3 and D2 Receptors

Effective rational drug discovery hinges on understanding the functional states of the target protein and distinguishing it from homologues. However, for the G protein coupled receptors, both activation-related conformational changes (ACCs) and intrinsic divergence among receptors can be misled or obscured by ligand-specific conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics (MD) simulation results of the receptors bound with various ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors, including the extracellular portion of TM5 (TM5e) and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found more outward tilting of TM6e in the D2R compared to the D3R in both the experimental structures and simulations bound with ligands in different scaffolds. However, this difference was drastically reduced in the simulations bound with nonselective agonist quinpirole, suggesting a misleading effect of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may also obscure intrinsic divergence. Importantly, our MD simulations revealed divergence in the dynamics of these receptors. Specifically, the D2R exhibited heightened flexibility compared to the D3R in the extracellular loops and TMs 5e, 6e, and 7e, associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting the D2R and D3R with more precise pharmacological profiles.

An intracellular hydrophobic nexus critical for hERG1 channel slow deactivation

Slow deactivation is a critical property of voltage-gated K+ channels encoded by the human Ether-à-go-go-Related Gene 1 (hERG). hERG1 channel deactivation is modulated by interactions between intracellular N-terminal Per-Arnt-Sim (PAS) and C-terminal cyclic nucleotide-binding homology (CNBh) domains. The PAS domain is multipartite, comprising a globular domain (gPAS; residues 26-135) and an N-terminal PAS-cap that is further subdivided into an initial unstructured "tip" (residues 1-12) and an amphipathic α-helical region (residues 13-25). Although the PAS-cap tip has long been considered the effector of slow deactivation, how its position near the gating machinery is controlled has not been elucidated. Here, we show that a triad of hydrophobic interactions among the gPAS, PAS-cap α helix, and the CNBh domains is required to support slow deactivation in hERG1. The primary sequence of this "hydrophobic nexus" is highly conserved among mammalian ERG channels but shows key differences to fast-deactivating Ether-à-go-go 1 (EAG1) channels. Combining sequence analysis, structure-directed mutagenesis, electrophysiology, and molecular dynamics simulations, we demonstrate that polar serine substitutions uncover an intermediate deactivation mode that is also mimicked by deletion of the PAS-cap α helix. Molecular dynamics simulation analyses of the serine-substituted channels show an increase in distance among the residues of the hydrophobic nexus, a rotation of the intracellular gating ring, and a retraction of the PAS-cap tip from its receptor site near the voltage sensor domain and channel gate. These findings provide compelling evidence that the hydrophobic nexus coordinates the respective components of the intracellular gating ring and positions the PAS-cap tip to control hERG1 deactivation gating.

Unraveling activation-related rearrangements and intrinsic divergence from ligand-induced conformational changes of the dopamine D3 and D2 receptors

Effective rational drug discovery targeting a specific protein hinges on understanding their functional states and distinguishing it from homologs. However, for the G protein coupled receptors, both the activation-related conformational changes (ACCs) and the intrinsic divergence among receptors can be misled or obscured by ligand-induced conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics simulation results of the receptors bound with different ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors including TM5e and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found pronounced outward tilting of TM6e in the D2R compared to the D3R in both experimental structures and simulations with ligands in different scaffolds. This tilting was drastically reduced in the simulations of the receptors bound with nonselective full agonist quinpirole, suggesting a misleading impact of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may obscure intrinsic divergence. In addition, our analysis showed that the impact of the nonconserved TM1 propagated to conserved Trp7.40 and Glu2.65, both are ligand binding residues. We also found that the D2R exhibited heightened flexibility compared to the D3R in the extracellular portions of TMs 5, 6, and 7, potentially associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting D2R or D3R with more precise pharmacological profiles.

A Novel "Activation Switch" Motif Common to All Aminergic Receptors

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

Macrocyclic "tet a"-Derived Cobalt(III) Complex with a N,N'-Disubstituted Hexadentate Ligand: Crystal Structure, Photonuclease Activity, and as a Photosensitizer

A cobalt(III) complex, [Co(L)]Cl (complex 1, where L = 1,8-[N,N-bis{(3-formyl-2-hydroxy-5-methyl)benzyl}]-1,4,8,11-tetraaza-5,5,7,12,12,14-hexamethylcyclotetradecane) with distorted octahedral geometry has been synthesized and characterized using various spectroscopic techniques. The structure of the ligand has remarkably rich hydrogen intermolecular interactions such as H···H, H···C/C···H, and H···O/O···H that vary with the presence of the metal ion, and the structure of complex 1 has Cl···H interactions; this result has been proved by Hirshfeld surface and two-dimensional (2D) fingerprint maps analyses. The complex exhibits a quasi-reversible Co(III)/Co(II) redox couple with E _1/2 = -0.76 V. Calf thymus DNA (CT DNA) binding abilities of the ligand and complex 1 were confirmed by spectroscopic and electrochemical analyses. According to absorption studies, the ligand and complex 1 bind to CT DNA via intercalative binding mode, with intrinsic binding strengths of 1.41 × 10³ and 8.64 × 10³ M^-1, respectively. A gel electrophoresis assay shows that complex 1 promotes the pUC19 DNA cleavage under dark and light irradiation conditions. Complex 1 has superior antimicrobial activity than the ligand. The cytotoxicity of complex 1 was tested against MDA-MB-231 breast cancer cells with values of IC₅₀ of 1.369 μg mL^-1 in the dark and 0.9034 μg mL^-1 after light irradiation. Besides, cell morphological studies confirmed the morphological changes with AO/EB dual staining, reactive oxygen species (ROS) staining, mitochondria staining, and Hoechst staining on MDA-MB-231 cancer cells by fluorescence microscopy. Complex 1 was found to be a potent antiproliferative agent against MDA-MB-231 cells, and it can induce mitochondrial-mediated and caspase-dependent apoptosis with activation of downregulated caspases. The biotoxicity assay of complex 1 on the development of Artemia nauplii was evaluated at an IC₅₀ value of 200 μg mL^-1 and with excellent biocompatibility.

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context

Aims

The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.

Methods

In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein–protein interface.

Results

Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike “boat-shaped” receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the “stern” or “bow” regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the “hull” region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.

Conclusion

In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein

G-protein coupled receptors (GPCRs) are known for their low stability and large conformational changes upon transitions between multiple states. A widely used method for stabilizing these receptors is to make chimeric receptors by fusing soluble proteins (i.e., fusion partner proteins) into the intracellular loop 3 (ICL3) connecting the transmembrane helices 5 and 6 (TM5 and TM6). However, this fusion approach requires experimental trial and error to identify appropriate soluble proteins, residue positions, and linker lengths for making the fusion. Moreover, this approach has not provided state-targeting stabilization of GPCRs. Here, to rationally stabilize a class A GPCR, adenosine A2A receptor (A2AR) in a target state, we carried out the custom-made de novo design of α-helical fusion partner proteins, which can fix the conformation of TM5 and TM6 to that in an inactive state of A2AR through straight helical connections without any kinks or intervening loops. The chimeric A2AR fused with one of the designs (FiX1) exhibited increased thermal stability. Moreover, compared with the wild type, the binding affinity of the chimera against the agonist NECA was significantly decreased, whereas that against the inverse agonist ZM241385 was similar, indicating that the inactive state was selectively stabilized. Our strategy contributes to the rational state-targeting stabilization of GPCRs.

A GX2GX3G motif facilitates acyl chain sequestration by Saccharomyces cerevisiae acyl carrier protein

Saccharomyces cerevisiae acyl carrier protein (ScACP) is a component of the large fungal fatty acid synthase I (FAS I) complex. ScACP comprises two subdomains: a conserved ACP domain that shares extensive structural homology with other ACPs and a unique structural domain. Unlike the metazoan type I ACP that does not sequester the acyl chain, ScACP can partially sequester the growing acyl chain within its hydrophobic core by a mechanism that remains elusive. Our studies on the acyl-ScACP intermediates disclose a unique ¹⁸⁸GX₂GX₃G¹⁹⁵ sequence in helix II important for ACP function. Complete loss of sequestration was observed upon mutation of the three glycines in this sequence to valine (G188V/G191V/G195V), while G191V and G188V/G191V double mutants displayed a faster rate of acyl chain hydrolysis. Likewise, mutation of Thr216 to Ala altered the size of the hydrophobic cavity, resulting in loss of C₁₂- chain sequestration. Combining NMR studies with insights from the crystal structure, we show that three glycines in helix II and a threonine in helix IV favor conformational change, which in turn generate space for acyl chain sequestration. Furthermore, we identified the primary hydrophobic cavity of ScACP, present between the carboxyl end of helix II and IV. The opening of the cavity lies between the second and third turns of helix II and loop II. Overall, the study highlights a novel role of the GX₂GX₃G motif in regulating acyl chain sequestration, vital for ScACP function.

Cryo-EM structure of the mature and infective Mayaro virus at 4.4 Å resolution reveals features of arthritogenic alphaviruses

Mayaro virus (MAYV) is an emerging arbovirus of the Americas that may cause a debilitating arthritogenic disease. The biology of MAYV is not fully understood and largely inferred from related arthritogenic alphaviruses. Here, we present the structure of MAYV at 4.4 Å resolution, obtained from a preparation of mature, infective virions. MAYV presents typical alphavirus features and organization. Interactions between viral proteins that lead to particle formation are described together with a hydrophobic pocket formed between E1 and E2 spike proteins and conformational epitopes specific of MAYV. We also describe MAYV glycosylation residues in E1 and E2 that may affect MXRA8 host receptor binding, and a molecular “handshake” between MAYV spikes formed by N262 glycosylation in adjacent E2 proteins. The structure of MAYV is suggestive of structural and functional complexity among alphaviruses, which may be targeted for specificity or antiviral activity.

Mayaro virus (MAYV) is an emerging arbovirus in Central and South America that is transmitted by mosquitoes and causes arthritogenic disease. Here, the authors present the 4.4 Å resolution cryo-EM structure of MAYV and describe specific features of the virus, which could be exploited for the design of MAYV-specific diagnostics and therapeutics.

Mapping temperature-dependent conformational change in the voltage-sensing domain of an engineered heat-activated K+ channel

Temperature-dependent regulation of ion channel activity is critical for a variety of physiological processes ranging from immune response to perception of noxious stimuli. Our understanding of the structural mechanisms that underlie temperature sensing remains limited, in part due to the difficulty of combining high-resolution structural analysis with temperature stimulus. Here, we use NMR to compare the temperature-dependent behavior of Shaker potassium channel voltage sensor domain (WT-VSD) to its engineered temperature sensitive (TS-VSD) variant. Further insight into the molecular basis for temperature-dependent behavior is obtained by analyzing the experimental results together with molecular dynamics simulations. Our studies reveal that the overall secondary structure of the engineered TS-VSD is identical to the wild-type channels except for local changes in backbone torsion angles near the site of substitution (V369S and F370S). Remarkably however, these structural differences result in increased hydration of the voltage-sensing arginines and the S4-S5 linker helix in the TS-VSD at higher temperatures, in contrast to the WT-VSD. These findings highlight how subtle differences in the primary structure can result in large-scale changes in solvation and thereby confer increased temperature-dependent activity beyond that predicted by linear summation of solvation energies of individual substituents.

Analysis of Sequence Divergence in Mammalian ABCGs Predicts a Structural Network of Residues That Underlies Functional Divergence

The five members of the mammalian G subfamily of ATP-binding cassette transporters differ greatly in their substrate specificity. Four members of the subfamily are important in lipid transport and the wide substrate specificity of one of the members, ABCG2, is of significance due to its role in multidrug resistance. To explore the origin of substrate selectivity in members 1, 2, 4, 5 and 8 of this subfamily, we have analysed the differences in conservation between members in a multiple sequence alignment of ABCG sequences from mammals. Mapping sets of residues with similar patterns of conservation onto the resolved 3D structure of ABCG2 reveals possible explanations for differences in function, via a connected network of residues from the cytoplasmic to transmembrane domains. In ABCG2, this network of residues may confer extra conformational flexibility, enabling it to transport a wider array of substrates.

Ligand-Specific Factors Influencing GLP-1 Receptor Post-Endocytic Trafficking and Degradation in Pancreatic Beta Cells

The glucagon-like peptide-1 receptor (GLP-1R) is an important regulator of blood glucose homeostasis. Ligand-specific differences in membrane trafficking of the GLP-1R influence its signalling properties and therapeutic potential in type 2 diabetes. Here, we have evaluated how different factors combine to control the post-endocytic trafficking of GLP-1R to recycling versus degradative pathways. Experiments were performed in primary islet cells, INS-1 832/3 clonal beta cells and HEK293 cells, using biorthogonal labelling of GLP-1R to determine its localisation and degradation after treatment with GLP-1, exendin-4 and several further GLP-1R agonist peptides. We also characterised the effect of a rare GLP1R coding variant, T149M, and the role of endosomal peptidase endothelin-converting enzyme-1 (ECE-1), in GLP1R trafficking. Our data reveal how treatment with GLP-1 versus exendin-4 is associated with preferential GLP-1R targeting towards a recycling pathway. GLP-1, but not exendin-4, is a substrate for ECE-1, and the resultant propensity to intra-endosomal degradation, in conjunction with differences in binding affinity, contributes to alterations in GLP-1R trafficking behaviours and degradation. The T149M GLP-1R variant shows reduced signalling and internalisation responses, which is likely to be due to disruption of the cytoplasmic region that couples to intracellular effectors. These observations provide insights into how ligand- and genotype-specific factors can influence GLP-1R trafficking.

Three sodium channel mutations from Aedes albopictus confer resistance to Type I, but not Type II pyrethroids

Voltage-gated sodium channels are the major targets of several classes of insecticides, including pyrethroids. However, sensitivities of many insect pest species to pyrethroids have gradually decreased due to overuse in pest management programs. One major mechanism of pyrethroid resistance known as knockdown resistance (kdr) involves mutations in the sodium channel gene. Three new mutations in helix IIIS6 of sodium channel (I1532T and F1534S/L) are recently detected in several pyrethroid-resistant populations of Aedes albopictus. The roles of these mutations in pyrethroid resistance have not been functionally examined. We introduced mutations I1532T and F1534S/L alone or in combination into the pyrethroid-sensitive sodium channel AaNa_v1-1 from Aedes aegypti by site-directed mutagenesis and explored effects of these mutations on the channel gating and sensitivity to pyrethroids. No significant modifications in channel properties were detected, except for a slightly changed activation by F1534S and I1532T + F1534S. However, I1532T and F1534S/L substantially reduced the channel sensitivity to Type I pyrethroids, permethrin and bifenthrin, but not to two Type II pyrethroids, deltamethrin and cypermethrin. The double mutations did not increase the channel resistance to permethrin or bifenthrin. We have built a Na_v1.4-based homology model of the AaNa_v1-1 channel and docked pyrethroids in the model to explain different sensitivities of the mutants to Type I and Type II pyrethroids. The results will assist in developing molecular markers for monitoring pest resistance to pyrethroids. They also provide new insight in the molecular basis of different action of Type I and Type II pyrethroids on sodium channels.

GPR6 Structural Insights: Homology Model Construction and Docking Studies

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

Helminth Defense Molecules as Design Templates for Membrane Active Antibiotics

A design template for membrane active antibiotics against microbial and tumor cells is described. The template is an amino acid sequence that combines the properties of helminth defense molecules, which are not cytolytic, with the properties of host-defense peptides, which disrupt microbial membranes. Like helminth defense molecules, the template folds into an amphipathic helix in both mammalian host and microbial phospholipid membranes. Unlike these molecules, the template exhibits antimicrobial and anticancer properties that are comparable to those of antimicrobial and anticancer antibiotics. The selective antibiotic activity of the template builds upon a functional synergy between three distinctive faces of the helix, which is in contrast to two faces of membrane-disrupting amphipathic structures. This synergy enables the template to adapt pore formation mechanisms according to the nature of the target membrane, inducing the lysis of microbial and tumor cells.

Inter-residue interactions in alpha-helical transmembrane proteins

MOTIVATION

The number of available membrane protein structures has markedly increased in the last years and, in parallel, the reliability of the methods to detect transmembrane (TM) segments. In the present report, we characterized inter-residue interactions in α-helical membrane proteins using a dataset of 3462 TM helices from 430 proteins. This is by far the largest analysis published to date.

RESULTS

Our analysis of residue-residue interactions in TM segments of membrane proteins shows that almost all interactions involve aliphatic residues and Phe. There is lack of polar-polar, polar-charged and charged-charged interactions except for those between Thr or Ser sidechains and the backbone carbonyl of aliphatic and Phe residues. The results are discussed in the context of the preferences of amino acids to be in the protein core or exposed to the lipid bilayer and to occupy specific positions along the TM segment. Comparison to datasets of β-barrel membrane proteins and of α-helical globular proteins unveils the specific patterns of interactions and residue composition characteristic of α-helical membrane proteins that are the clue to understanding their structure.

AVAILABILITY AND IMPLEMENTATION

Results data and datasets used are available at http://lmc.uab.cat/TMalphaDB/interactions.php.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bactericidal Disruption of Magnesium Metallostasis in Mycobacterium tuberculosis Is Counteracted by Mutations in the Metal Ion Transporter CorA

Antimycobacterial agents might shorten the course of treatment by reducing the number of phenotypically tolerant bacteria if they could kill M. tuberculosis in diverse metabolic states. Here we report two chemically disparate classes of agents that kill M. tuberculosis both when it is replicating and when it is not. Under replicating conditions, the tricyclic 4-hydroxyquinolines and a barbituric acid analogue deplete intrabacterial magnesium as a mechanism of action, and for both compounds, mutations in CorA, a putative Mg²⁺/Co²⁺ transporter, conferred resistance to the compounds when M. tuberculosis was under replicating conditions but not under nonreplicating conditions, illustrating that a given compound can kill M. tuberculosis in different metabolic states by disparate mechanisms. Targeting magnesium metallostasis represents a previously undescribed antimycobacterial mode of action that might cripple M. tuberculosis in a Mg²⁺-deficient intraphagosomal environment of macrophages.

A defining characteristic of treating tuberculosis is the need for prolonged administration of multiple drugs. This may be due in part to subpopulations of slowly replicating or nonreplicating Mycobacterium tuberculosis bacilli exhibiting phenotypic tolerance to most antibiotics in the standard treatment regimen. Confounding this problem is the increasing incidence of heritable multidrug-resistant M. tuberculosis. A search for new antimycobacterial chemical scaffolds that can kill phenotypically drug-tolerant mycobacteria uncovered tricyclic 4-hydroxyquinolines and a barbituric acid derivative with mycobactericidal activity against both replicating and nonreplicating M. tuberculosis. Both families of compounds depleted M. tuberculosis of intrabacterial magnesium. Complete or partial resistance to both chemotypes arose from mutations in the putative mycobacterial Mg²⁺/Co²⁺ ion channel, CorA. Excess extracellular Mg²⁺, but not other divalent cations, diminished the compounds’ cidality against replicating M. tuberculosis. These findings establish depletion of intrabacterial magnesium as an antimicrobial mechanism of action and show that M. tuberculosis magnesium homeostasis is vulnerable to disruption by structurally diverse, nonchelating, drug-like compounds.

Towards A Molecular Understanding of The Cannabinoid Related Orphan Receptor GPR18: A Focus on Its Constitutive Activity

The orphan G-protein coupled receptor (GPCR), GPR18, has been recently proposed as a potential member of the cannabinoid family as it recognizes several endogenous, phytogenic, and synthetic cannabinoids. Potential therapeutic applications for GPR18 include intraocular pressure, metabolic disorders, and cancer. GPR18 has been reported to have high constitutive activity, i.e., activation/signaling occurs in the absence of an agonist. This activity can be reduced significantly by the A3.39N mutation. At the intracellular (IC) ends of (transmembrane helices) TMH3 and TMH6 in GPCRs, typically, a pair of oppositely charged amino acids form a salt bridge called the "ionic lock". Breaking of this salt bridge creates an IC opening for coupling with G protein. The GPR18 "ionic lock" residues (R3.50/S6.33) can form only a hydrogen bond. In this paper, we test the hypothesis that the high constitutive activity of GPR18 is due to the weakness of its "ionic lock" and that the A3.39N mutation strengthens this lock. To this end, we report molecular dynamics simulations of wild-type (WT) GPR18 and the A3.39N mutant in fully hydrated (POPC) phophatidylcholine lipid bilayers. Results suggest that in the A3.39N mutant, TMH6 rotates and brings R3.50 and S6.33 closer together, thus strengthening the GPR18 "ionic lock".

Functional and pharmacological characterization of an S5 domain hERG mutation associated with short QT syndrome

Congenital short QT syndrome (SQTS) is a repolarization disorder characterized by abbreviated QT intervals, atrial and ventricular arrhythmias and a risk of sudden death. This study characterized a missense mutation (I560T) in the S5 domain of the hERG K⁺ channel that has been associated with variant 1 of the SQTS. Whole cell patch clamp recordings of wild-type (WT) and I560T hERG current (I_hERG) were made at 37 °C from hERG expressing HEK 293 cells, and the structural context of the mutation was investigated using a recently reported cryo-EM structure of hERG. Under conventional voltage clamp, the I560T mutation increased I_hERG amplitude without altering the voltage-dependence of activation, although it accelerated activation time-course and also slowed deactivation time-course at some voltages. The voltage dependence of I_hERG inactivation was positively shifted (by ∼24 mV) and the time-course of inactivation was slowed by the I560T mutation. There was also a modest decrease in K⁺ over Na⁺ ion selectivity with the I560T mutation. Under action potential (AP) voltage clamp, the net charge carried by hERG was significantly increased during ventricular, Purkinje fibre and atrial APs, with maximal I_hERG also occurring earlier during the plateau phase of ventricular and Purkinje fibre APs. The I560T mutation exerted only a modest effect on quinidine sensitivity of I_hERG: the IC₅₀ for mutant I_hERG was 2.3 fold that for WT I_hERG under conventional voltage clamp. Under AP voltage clamp the inhibitory effect of 1 μM quinidine was largely retained for I560T hERG and the timing of peak I560T I_hERG was altered towards that of the WT channel. In both the open channel structure and a closed hERG channel model based on the closely-related EAG structure, I560T side-chains were oriented towards membrane lipid and away from adjacent domains of the channel, contrasting with previous predictions based on homology modelling. In summary, the I560T mutation produces multiple effects on hERG channel operation that result in a gain-of-function that is expected to abbreviate ventricular, atrial and Purkinje fibre repolarization. Quinidine is likely to be of value in offsetting the increase in I_hERG and altered I_hERG timing during ventricular APs in SQTS with this mutation.

Mechanisms Underlying Allosteric Molecular Switches of Metabotropic Glutamate Receptor 5

The metabotropic glutamate 5 (mGlu₅) receptor is a class C G protein-coupled receptor (GPCR) that is implicated in several CNS disorders making it a popular drug discovery target. Years of research have revealed allosteric mGlu₅ ligands showing an unexpected complete switch in functional activity despite only small changes in their chemical structure, resulting in positive allosteric modulators (PAM) or negative allosteric modulators (NAM) for the same scaffold. Up to now, the origins of this effect are not understood, causing difficulties in a drug discovery context. In this work, experimental data was gathered and analyzed alongside docking and Molecular Dynamics (MD) calculations for three sets of PAM and NAM pairs. The results consistently show the role of specific interactions formed between ligand substituents and amino acid side chains that block or promote local movements associated with receptor activation. The work provides an explanation for how such small structural changes lead to remarkable differences in functional activity. While this work can greatly help drug discovery programs avoid these switches, it also provides valuable insight into the mechanisms of class C GPCR allosteric activation. Furthermore, the approach shows the value of applying MD to understand functional activity in drug design programs, even for such close structural analogues.

Bipolar switching by HCN voltage sensor underlies hyperpolarization activation

Hyperpolarization-activated, cyclic nucleotide-gated (HCN) channels show an inverted voltage response compared with virtually all other voltage-gated channels, opening on hyperpolarization rather than depolarization. Although the structure of the HCN1 channel was recently solved, the structural element(s) responsible for the inverted gating polarity of HCN is not known. Here, we use a hierarchical approach, by first characterizing the functional contribution of each structural element to channel gating, and then identifying the critical interactions between these elements. Our studies reveal that the HCN voltage sensor can gate the same pore open on both depolarization and hyperpolarization, thereby acting as a bipolar switch. Elements in the pore domain shut off the depolarization-activation pathway in wild-type channels.

Despite sharing a common architecture with archetypal voltage-gated ion channels (VGICs), hyperpolarization- and cAMP-activated ion (HCN) channels open upon hyperpolarization rather than depolarization. The basic motions of the voltage sensor and pore gates are conserved, implying that these domains are inversely coupled in HCN channels. Using structure-guided protein engineering, we systematically assembled an array of mosaic channels that display the full complement of voltage-activation phenotypes observed in the VGIC superfamily. Our studies reveal that the voltage sensor of the HCN channel has an intrinsic ability to drive pore opening in either direction and that the extra length of the HCN S4 is not the primary determinant for hyperpolarization activation. Tight interactions at the HCN voltage sensor–pore interface drive the channel into an hERG-like inactivated state, thereby obscuring its opening upon depolarization. This structural element in synergy with the HCN cyclic nucleotide-binding domain and specific interactions near the pore gate biases the channel toward hyperpolarization-dependent opening. Our findings reveal an unexpected common principle underpinning voltage gating in the VGIC superfamily and identify the essential determinants of gating polarity.

Comparison of Peptide Ion Conformers Arising from Non-Helical and Helical Peptides Using Ion Mobility Spectrometry and Gas-Phase Hydrogen/Deuterium Exchange

The dominant gas-phase conformer of [M+3H]³⁺ ions of the model peptide acetyl-PSSSSKSSSSKSSSSKSSSSK has been examined with ion mobility spectrometry (IMS), gas-phase hydrogen deuterium exchange (HDX), and mass spectrometry (MS) techniques. The [M+3H]³⁺ peptide ions are observed predominantly as a relatively compact conformer type. Upon subjecting these ions to electron transfer dissociation (ETD), the level of protection for each amino acid residue in the peptide sequence is assessed. The overall per-residue deuterium uptake is observed to be relatively more efficient for the neutral residues than for the model peptide acetyl-PAAAAKAAAAKAAAAKAAAAK. In comparison, the N-terminal and C-terminal regions of the serine peptide show greater relative protection compared with interior residues. Molecular dynamics (MD) simulations have been used to generate candidate structures for collision cross section and HDX reactivity matching. Hydrogen accessibility scoring (HAS) for select structural candidates from MD simulations has been used to suggest conformer types that could contribute to the observed HDX patterns. The results are discussed with respect to recent studies employing extensive MD simulations of gas-phase structure establishment of a peptide system. Graphical Abstract ᅟ.

A cancer-associated polymorphism in ESCRT-III disrupts the abscission checkpoint and promotes genome instability

Cytokinetic abscission facilitates the irreversible separation of daughter cells. This process requires the endosomal-sorting complexes required for transport (ESCRT) machinery and is tightly regulated by charged multivesicular body protein 4C (CHMP4C), an ESCRT-III subunit that engages the abscission checkpoint (NoCut) in response to mitotic problems such as persisting chromatin bridges within the midbody. Importantly, a human polymorphism in CHMP4C (rs35094336, CHMP4CT232) increases cancer susceptibility. Here, we explain the structural and functional basis for this cancer association: The CHMP4CT232 allele unwinds the C-terminal helix of CHMP4C, impairs binding to the early-acting ESCRT factor ALIX, and disrupts the abscission checkpoint. Cells expressing CHMP4CT232 exhibit increased levels of DNA damage and are sensitized to several conditions that increase chromosome missegregation, including DNA replication stress, inhibition of the mitotic checkpoint, and loss of p53. Our data demonstrate the biological importance of the abscission checkpoint and suggest that dysregulation of abscission by CHMP4CT232 may synergize with oncogene-induced mitotic stress to promote genomic instability and tumorigenesis.

Transmembrane domain dependent inhibitory function of FcγRIIB

FcγRIIB, the only inhibitory IgG Fc receptor, functions to suppress the hyper-activation of immune cells. Numerous studies have illustrated its inhibitory function through the ITIM motif in the cytoplasmic tail of FcγRIIB. However, later studies revealed that in addition to the ITIM, the transmembrane (TM) domain of FcγRIIB is also indispensable for its inhibitory function. Indeed, recent epidemiological studies revealed that a non-synonymous single nucleotide polymorphism (rs1050501) within the TM domain of FcγRIIB, responsible for the I232T substitution, is associated with the susceptibility to systemic lupus erythematosus (SLE). In this review, we will summarize these epidemiological and functional studies of FcγRIIB-I232T in the past few years, and will further discuss the mechanisms accounting for the functional loss of FcγRIIB-I232T. Our review will help the reader gain a deeper understanding of the importance of the TM domain in mediating the inhibitory function of FcγRIIB and may provide insights to a new therapeutic target for the associated diseases.

Both positional and chemical variables control in vitro proteolytic cleavage of a presenilin ortholog

Mechanistic details of intramembrane aspartyl protease (IAP) chemistry, which is central to many biological and pathogenic processes, remain largely obscure. Here, we investigated the in vitro kinetics of a microbial intramembrane aspartyl protease (mIAP) fortuitously acting on the renin substrate angiotensinogen and the C-terminal transmembrane segment of amyloid precursor protein (C100), which is cleaved by the presenilin subunit of γ-secretase, an Alzheimer disease (AD)-associated IAP. mIAP variants with substitutions in active-site and putative substrate-gating residues generally exhibit impaired, but not abolished, activity toward angiotensinogen and retain the predominant cleavage site (His-Thr). The aromatic ring, but not the hydroxyl substituent, within Tyr of the catalytic Tyr-Asp (YD) motif plays a catalytic role, and the hydrolysis reaction incorporates bulk water as in soluble aspartyl proteases. mIAP hydrolyzes the transmembrane region of C100 at two major presenilin cleavage sites, one corresponding to the AD-associated Aβ42 peptide (Ala-Thr) and the other to the non-pathogenic Aβ48 (Thr-Leu). For the former site, we observed more favorable kinetics in lipid bilayer-mimicking bicelles than in detergent solution, indicating that substrate-lipid and substrate-enzyme interactions both contribute to catalytic rates. High-resolution MS analyses across four substrates support a preference for threonine at the scissile bond. However, results from threonine-scanning mutagenesis of angiotensinogen demonstrate a competing positional preference for cleavage. Our results indicate that IAP cleavage is controlled by both positional and chemical factors, opening up new avenues for selective IAP inhibition for therapeutic interventions.

Structure and function of yeast and fungal Na+ /H+ antiporters

Sodium proton antiporters (or sodium proton exchangers [NHEs]) are a critical family of membrane proteins that exchange sodium for protons across cell membranes. In yeast and plants, their primary function is to keep the sodium concentration low inside the cytoplasm. One class of NHE constitutively expressed in yeast is the plasma membrane Na⁺ /H⁺ antiporter, and another class is expressed on the endosomal/vacuolar membrane. At present, four bacterial plasma membrane antiporter structures are known and nuclear magnetic resonance structures are available for the membrane spanning transmembrane helices of mammalian and yeast NHEs. Additionally, a vast amount of mutational data are available on the role of individual amino acids and critical motifs involved in transport. We combine this information to obtain a more detailed picture of the yeast NHE plasma membrane protein and review mechanisms of transport, conserved motifs, unique residues important in function, and regulation of these proteins. The Na⁺ /H⁺ antiporter of Schizosaccharomyces pombe, SpNHE1, is an interesting model protein in an easy to study system and is representative of fungal Na⁺ /H⁺ antiporters. © IUBMB Life, 70(1):23-31, 2018.

Conformational states of HAMP domains interacting with sensory rhodopsin membrane systems: an integrated all-atom and coarse-grained molecular dynamics simulation approach

Understanding the downstream signaling mechanism of sensory rhodopsin and its cognate transducer complex (srII-htrII) has long been a challenge in the field of photoreceptor research. Here, an integration of all-atom and coarse-grained (CG) molecular dynamics (MD) simulations in different srII-htrII complex states is carried out. It is shown that the cytoplasmic four-helix HAMP dimer gives rise to a gear-box model interaction with discrete hydrophobic packing in Natronomonas pharaonis (Np). Structural analysis in all-atom and CG-MD reveals a stable conformational state in the physiological environment (323 K and 1.15 M salt). Comparative analysis in the ground and intermediate state conformations reveals substantial inter-HAMP interactions in the intermediate state with uniform clockwise (+10° to +30°) and counterclockwise (-20° to -40°) rotations in the α1 helix and the α2 helix of the monomer, respectively. Low temperature and low salt environments (283 K and 0.15 M) significantly affect srII-htrII binding affinity in both states with unusual helix bending. The distinguished control cable, knob-into-holes packing and piston-like movements in HAMP helices are found in the intermediate state complex. The N-terminal htrII (159 residues) coupled with srII yields a binding energy (ΔG_bind) of -309.22, -436.53 and -331.11 kJ mol^-1 in the MM/PBSA calculation for the NphtrII homodimer, the NpsrII-htrII ground state conformation and the NpsrII-htrII intermediate state conformation, respectively. Only the HAMP1 domain shows a very low ΔG_bind value (-21.03 kJ mol^-1) for the ground state in comparison to that for the intermediate state (-54.68 kJ mol^-1). The structural analysis highlights the key residues that include Y199^srII, T189^srII, E43^htrII, T86^htrII, M100^htrII, E116^htrII, E126^htrII and S130^htrII for complex stabilization and signal transduction.

Methods for the Development of In Silico GPCR Models

The Reggio group has constructed computer models of the inactive and G-protein-activated states of the cannabinoid CB1 and CB2 receptors, as well as, several orphan receptors that recognize a subset of cannabinoid compounds, including GPR55 and GPR18. These models have been used to design ligands, mutations, and covalent labeling studies. The resultant second-generation models have been used to design ligands with improved affinity, efficacy, and subtype selectivity. Herein, we provide a guide for the development of GPCR models using the most recent orphan receptor studied in our lab, GPR3. GPR3 is an orphan receptor that belongs to the Class A family of G-protein-coupled receptors. It shares high sequence similarity with GPR6, GPR12, the lysophospholipid receptors, and the cannabinoid receptors. GPR3 is predominantly expressed in mammalian brain and oocytes and it is known as a Gαs-coupled receptor activated constitutively in cells. GPR3 represents a possible target for the treatment of different pathological conditions such as Alzheimer's disease, oocyte maturation, or neuropathic pain. However, the lack of potent and selective GPR3 ligands is delaying the exploitation of this promising therapeutic target. In this context, we aim to develop a homology model that helps us to elucidate the structural determinants governing ligand-receptor interactions at GPR3. In this chapter, we detail the methods and rationale behind the construction of the GPR3 active-and inactive-state models. These homology models will enable the rational design of novel ligands, which may serve as research tools for further understanding of the biological role of GPR3.

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.

Adjacent channelrhodopsin-2 residues within transmembranes 2 and 7 regulate cation selectivity and distribution of the two open states

Channelrhodopsin-2 (ChR2) is a light-activated channel that can conduct cations of multiple valencies down the electrochemical gradient. Under continuous light exposure, ChR2 transitions from a high-conducting open state (O1) to a low-conducting open state (O2) with differing ion selectivity. The molecular basis for the O1 → O2 transition and how ChR2 modulates selectivity between states is currently unresolved. To this end, we used steered molecular dynamics, electrophysiology, and kinetic modeling to identify residues that contribute to gating and selectivity in discrete open states. Analysis of steered molecular dynamics experiments identified three transmembrane residues (Val-86, Lys-93, and Asn-258) that form a putative barrier to ion translocation. Kinetic modeling of photocurrents generated from ChR2 proteins with conservative mutations at these positions demonstrated that these residues contribute to cation selectivity (Val-86 and Asn-258), the transition between the two open states (Val-86), open channel stability, and the hydrogen-bonding network (K93I and K93N). These results suggest that this approach can be used to identify residues that contribute to the open-state transitions and the discrete ion selectivity within these states. With the rise of ChR2 use in optogenetics, it will be critical to identify residues that contribute to O1 or O2 selectivity and gating to minimize undesirable effects.

Impairment on the lateral mobility induced by structural changes underlies the functional deficiency of the lupus-associated polymorphism FcγRIIB-T232

Xu et al. show that the lupus-associated polymorphism FcγRIIB-T232 has structural changes of the TM domain that reduces lateral mobility and inhibitory functions.

FcγRIIB functions to suppress the activation of immune cells. A single-nucleotide polymorphism in the transmembrane (TM) domain of FcγRIIB, FcγRIIB-T232, is associated with lupus. In this study, we investigated the pathogenic mechanism of FcγRIIB-T232 at both functional and structural levels. Our results showed that FcγRIIB-T232 exhibited significantly reduced lateral mobility compared with FcγRIIB-I232 and was significantly less enriched into the microclusters of immune complexes (ICs) after stimulation. However, if sufficient responding time is given for FcγRIIB-T232 to diffuse and interact with the ICs, FcγRIIB-T232 can restore its inhibitory function. Moreover, substituting the FcγRIIB-T232 TM domain with that of a fast floating CD86 molecule restored both the rapid mobility and the inhibitory function, which further corroborated the importance of fast mobility for FcγRIIB to function. Mechanistically, the crippled lateral mobility of FcγRIIB-T232 can be explained by the structural changes of the TM domain. Both atomistic simulations and nuclear magnetic resonance measurement indicated that the TM helix of FcγRIIB-T232 exhibited a more inclined orientation than that of FcγRIIB-I232, thus resulting in a longer region embedded in the membrane. Therefore, we conclude that the single-residue polymorphism T232 enforces the inclination of the TM domain and thereby reduces the lateral mobility and inhibitory functions of FcγRIIB.

Structural Basis of the Disorder in the Tandem Zinc Finger Domain of the RNA-Binding Protein Tristetraprolin

Tristetraprolin (TTP) and TIS11d are two human RNA-binding proteins that belong to the CCCH-type tandem zinc finger family. In the RNA-free state, TIS11d coordinates a zinc ion in each of its two fingers, while TTP coordinates a single zinc ion with the N-terminal zinc finger. We have previously identified three residues, located in the C-terminal half of a short α-helix in the second zinc finger, that control how structured the RNA-binding domain is in these two proteins: Y151, L152, and Q153 in TTP and H201, T202, and I203 in TIS11d. Here, we have used molecular dynamics, NMR spectroscopy, and other biochemical methods to investigate the role of these three residues in the stability of the RNA-binding domain. We found that the intrahelical hydrogen bond formed by the T202 hydroxyl group in the C-terminal zinc finger of TIS11d is necessary to allow for π-π stacking between the side chains of a conserved phenylalanine and the zinc-coordinating histidine. We demonstrated that the lack of this hydrogen bond in TTP is responsible for the reduced zinc affinity of the C-terminal zinc finger.

Understanding Sodium Channel Function and Modulation Using Atomistic Simulations of Bacterial Channel Structures

Sodium channels are chief proteins involved in electrical signaling in the nervous system, enabling critical functions like heartbeat and brain activity. New high-resolution X-ray structures for bacterial sodium channels have created an opportunity to see how these proteins operate at the molecular level. An important challenge to overcome is establishing relationships between the structures and functions of mammalian and bacterial channels. Bacterial sodium channels are known to exhibit the main structural features of their mammalian counterparts, as well as several key functional characteristics, including selective ion conduction, voltage-dependent gating, pore-based inactivation and modulation by local anesthetic, antiarrhythmic and antiepileptic drugs. Simulations have begun to shed light on each of these features in the past few years. Despite deviations in selectivity signatures for bacterial and mammalian channels, simulations have uncovered the nature of the multiion conduction mechanism associated with Na(+) binding to a high-field strength site established by charged glutamate side chains. Simulations demonstrated a surprising level of flexibility of the protein, showing that these side chains are active participants in the permeation process. They have also uncovered changes in protein structure, leading to asymmetrical collapses of the activation gate that have been proposed to correspond to inactivated structures. These observations offer the potential to examine the mechanisms of state-dependent drug activity, focusing on pore-blocking and pore-based slow inactivation in bacterial channels, without the complexities of inactivation on multiple timescales seen in eukaryotic channels. Simulations have provided molecular views of the interactions of drugs, consistent with sites predicted in mammalian channels, as well as a wealth of other sites as potential new drug targets. In this chapter, we survey the new insights into sodium channel function that have emerged from studies of simpler bacterial channels, which provide an excellent learning platform, and promising avenues for mechanistic discovery and pharmacological development.

Inhibitor Discovery for the Human GLUT1 from Homology Modeling and Virtual Screening

The human Glucose Transporter 1 (hGLUT1 or SLC2A1) is a facilitative membrane transporter found in the liver, intestines, kidney, and brain, where it transports sugars such as d-glucose and d-galactose. Genetic variations in hGLUT1 are associated with a broad range of diseases and metabolic disorders. For example, hGLUT1 is upregulated in various cancer types (e.g., breast carcinoma) to support the increased anaerobic glycolysis and the Warburg effect. Thus, hGLUT1 is an emerging therapeutic target, which also transports commonly used cancer biomarkers (e.g., (18)F-DG). In this study, we use computational prediction followed by experimental testing, to characterize hGLUT1. We construct homology models of hGLUT1 in a partially occluded outward open ("occluded") conformation based on the X-ray structure of the E. coli xylose transporter, XylE. Comparison of the binding site of the occluded models to experimentally determined hGLUT structures revealed a hydrophobic pocket adjacent to the sugar-binding site, which was tested experimentally via site-directed mutagenesis. Virtual screening of various libraries of purchasable compounds against the occluded models, followed by experimental testing with cellular assays revealed seven previously unknown hGLUT1 ligands with IC50 values ranging from 0.45 μM to 59 μM. These ligands represent three unique chemotypes that are chemically different from any other known hGLUT1 ligands. The newly characterized hydrophobic pocket can potentially be utilized by the new ligands for increased affinity. Furthermore, the previously unknown hGLUT1 ligands can serve as chemical tools to further characterize hGLUT1 function or lead molecules for future drug development.

Interactions of a class IIb bacteriocin with a model lipid bilayer, investigated through molecular dynamics simulations

The emergence of antibiotic resistant microorganisms poses an alarming threat to global health. Antimicrobial peptides (AMPs) are considered a possible effective alternative to conventional antibiotic therapies. An understanding of the mechanism of action of AMPs is needed in order to better control and optimize their bactericidal activity. Plantaricin EF is a heterodimeric AMP, consisting of two peptides Plantaricin E (PlnE) and Plantaricin F (PlnF). We studied the behavior of these peptides on the surface of a model lipid bilayer. We identified the residues that facilitate peptide-peptide interactions. We also identified residues that mediate interactions of the dimer with the membrane. PlnE interacts with the membrane through amino acids at both its termini, while only the N terminus of PlnF approaches the membrane. By comparing the activity of single-site mutants of the two-peptide bacteriocin and the simulations of the bacteriocin on the surface of a model lipid bilayer, structure activity relationships are proposed. These studies allow us to generate hypotheses that relate biophysical interactions observed in simulations with the experimentally measured activity. We find that single-site amino acid substitutions result in markedly stronger antimicrobial activity when they strengthen the interactions between the two peptides, while, concomitantly, they weaken peptide-membrane association. This effect is more pronounced in the case of the PlnE mutant (G20A), which interacts the strongest with PlnF and the weakest with the membrane while displaying the highest activity.

Unusual H-Bond Topology and Bifurcated H-bonds in the 2-Fluoroethanol Trimer

By using a combination of rotational spectroscopy and ab initio calculations, an unusual H-bond topology was revealed for the 2-fluoroethanol trimer. The trimer exhibits a strong heterochiral preference and adopts an open OH⋅⋅⋅OH H-bond topology while utilizing two types of bifurcated H-bonds involving organic fluorine. This is in stark contrast to the cyclic OH⋅⋅⋅OH H-bond topology adopted by trimers of water and other simple alcohols. The strengths of different H-bonds in the trimer were analyzed by using the quantum theory of atoms in molecules. The study showcases a remarkable example of a chirality-induced switch in H-bond topology in a simple transient chiral fluoroalcohol. It provides important insight into the H-bond topologies of small fluoroalcohol aggregates, which are proposed to play a key role in protein folding and in enantioselective reactions and separations where fluoroalcohols serve as a (co)solvent.