Contiguously hydrophobic sequences are functionally significant throughout the human exome.

Hydrophobic interactions have long been established as essential for stabilizing struc-tured proteins as well as drivers of aggregation, but the impact of hydrophobicity on thefunctional significance of sequence variants has rarely been considered in a genome-wide context. Here we test the role of hydrophobicity on functional impact across70,000 disease- and non­disease-associated single-nucleotide polymorphisms (SNPs),using enrichment of disease association as an indicator of functionality. We find thatfunctional impact is uncorrelated with hydrophobicity of the SNP itself and only weaklycorrelated with the average local hydrophobicity, but is strongly correlated with boththe size and minimum hydrophobicity of the contiguously hydrophobic sequence (or"blob") that contains the SNP. Disease association is found to vary by more than sixfoldas a function of contiguous hydrophobicity parameters, suggesting utility as a prior foridentifying causal variation. We further find signatures of differential selective constrainton hydrophobic blobs and that SNPs splitting a long hydrophobic blob or joiningtwo short hydrophobic blobs are particularly likely to be disease associated. Trends arepreserved for both aggregating and nonaggregating proteins, indicating that the role ofcontiguous hydrophobicity extends well beyond aggregation risk.

A lipid site shapes the agonist response of a pentameric ligand-gated ion channel.

Phospholipids are key components of cellular membranes and are emerging as important functional regulators of different membrane proteins, including pentameric ligand-gated ion channels (pLGICs). Here, we take advantage of the prokaryote channel ELIC (Erwinia ligand-gated ion channel) as a model to understand the determinants of phospholipid interactions in this family of receptors. A high-resolution structure of ELIC in a lipid-bound state reveals a phospholipid site at the lower half of pore-forming transmembrane helices M1 and M4 and at a nearby site for neurosteroids, cholesterol or general anesthetics. This site is shaped by an M4-helix kink and a Trp-Arg-Pro triad that is highly conserved in eukaryote GABAA/C and glycine receptors. A combined approach reveals that M4 is intrinsically flexible and that M4 deletions or disruptions of the lipid-binding site accelerate desensitization in ELIC, suggesting that lipid interactions shape the agonist response. Our data offer a structural context for understanding lipid modulation in pLGICs.

Spontaneous lipid binding to the nicotinic acetylcholine receptor in a native membrane.

The nicotinic acetylcholine receptor (nAChR) and other pentameric ligand-gated ion channels are native to neuronal membranes with an unusual lipid composition. While it is well-established that these receptors can be significantly modulated by lipids, the underlying mechanisms have been primarily studied in model membranes with few lipid species. Here, we use coarse-grained molecular dynamics simulation to probe specific binding of lipids in a complex quasi-neuronal membrane. We ran a total of 50 µs of simulations of a single nAChR in a membrane composed of 36 species of lipids. Competition between multiple lipid species produces a complex distribution. We find that overall, cholesterol selects for concave inter-subunit sites and polyunsaturated fatty acids select for convex M4 sites, while monounsaturated and saturated lipids are unenriched in the nAChR boundary. We propose the "density-threshold affinity" as a metric calculated from continuous density distributions, which reduces to a standard affinity in two-state binding. We find that the density-threshold affinity for M4 weakens with chain rigidity, which suggests that flexible chains may help relax packing defects caused by the conical protein shape. For any site, PE headgroups have the strongest affinity of all phospholipid headgroups, but anionic lipids still yield moderately high affinities for the M4 sites as expected. We observe cooperative effects between anionic headgroups and saturated chains at the M4 site in the inner leaflet. We also analyze affinities for individual anionic headgroups. When combined, these insights may reconcile several apparently contradictory experiments on the role of anionic phospholipids in modulating nAChR.

Sequence specificity despite intrinsic disorder: How a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein.

The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 µs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies "blobs" of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces ß-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form ß-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction.

Untangling Direct and Domain-Mediated Interactions Between Nicotinic Acetylcholine Receptors in DHA-Rich Membranes.

At the neuromuscular junction (NMJ), the nicotinic acetylcholine receptor (nAChR) self-associates to give rise to rapid muscle movement. While lipid domains have maintained nAChR aggregates in vitro, their specific roles in nAChR clustering are currently unknown. In the present study, we carried out coarse-grained molecular dynamics simulations (CG-MD) of 1-4 nAChR molecules in two membrane environments: one mixture containing domain-forming, homoacidic lipids, and a second mixture consisting of heteroacidic lipids. Spontaneous dimerization of nAChRs was up to ten times more likely in domain-forming membranes; however, the effect was not significant in four-protein systems, suggesting that lipid domains are less critical to nAChR oligomerization when protein concentration is higher. With regard to lipid preferences, nAChRs consistently partitioned into liquid-disordered domains occupied by the omega-3 ([Formula: see text]-3) fatty acid, docosahexaenoic acid (DHA); enrichment of DHA boundary lipids increased with protein concentration, particularly in homoacidic membranes. This result suggests dimer formation blocks access of saturated chains and cholesterol, but not polyunsaturated chains, to boundary lipid sites.

A Novel Bifunctional Alkylphenol Anesthetic Allows Characterization of γ-Aminobutyric Acid, Type A (GABAA), Receptor Subunit Binding Selectivity in Synaptosomes.

Propofol, an intravenous anesthetic, is a positive modulator of the GABAA receptor, but the mechanistic details, including the relevant binding sites and alternative targets, remain disputed. Here we undertook an in-depth study of alkylphenol-based anesthetic binding to synaptic membranes. We designed, synthesized, and characterized a chemically active alkylphenol anesthetic (2-((prop-2-yn-1-yloxy)methyl)-5-(3-(trifluoromethyl)-3H-diazirin-3-yl)phenol, AziPm-click (1)), for affinity-based protein profiling (ABPP) of propofol-binding proteins in their native state within mouse synaptosomes. The ABPP strategy captured ∼4% of the synaptosomal proteome, including the unbiased capture of five α or ß GABAA receptor subunits. Lack of γ2 subunit capture was not due to low abundance. Consistent with this, independent molecular dynamics simulations with alchemical free energy perturbation calculations predicted selective propofol binding to interfacial sites, with higher affinities for α/ß than γ-containing interfaces. The simulations indicated hydrogen bonding is a key component leading to propofol-selective binding within GABAA receptor subunit interfaces, with stable hydrogen bonds observed between propofol and α/ß cavity residues but not γ cavity residues. We confirmed this by introducing a hydrogen bond-null propofol analogue as a protecting ligand for targeted-ABPP and observed a lack of GABAA receptor subunit protection. This investigation demonstrates striking interfacial GABAA receptor subunit selectivity in the native milieu, suggesting that asymmetric occupancy of heteropentameric ion channels by alkylphenol-based anesthetics is sufficient to induce modulation of activity.

Direct Interactions of Cholesterol With Pentameric Ligand-Gated Ion Channels: Testable Hypotheses From Computational Predictions.

Several essential ionotropic neurotransmitter receptors, including the nicotinic acetylcholine receptor (nAChR) and gamma-aminobutyric acid (GABA) type A receptor (GABAAr), belong to the family of pentameric ligand-gated ion channels (pLGICs). Function of these receptors is particularly sensitive to their lipid environment, including cholesterol and cholesterol-derived neurosteroids. Direct structural data investigating interactions between sterols and pLGICs, as well as their role in modulatory mechanisms, are largely unavailable. Physics-based computational approaches can serve a vital role in interpretation of more indirect data as well as hypothesis generation and experimental design. In this chapter, I report several examples in which computational approaches were used to predict direct binding interactions of steroids and pLGICs, evaluate the relative likelihood of possible interpretations of experimental data, and present rationally designed simple experiments. I conclude by offering several predictions that could be tested by future experiments.

Shedding Light on Anesthetic Mechanisms: Application of Photoaffinity Ligands.

Anesthetic photoaffinity ligands have had an increasing presence within anesthesiology research. These ligands mimic parent general anesthetics and allow investigators to study anesthetic interactions with receptors and enzymes; identify novel targets; and determine distribution within biological systems. To date, nearly all general anesthetics used in medicine have a corresponding photoaffinity ligand represented in the literature. In this review, we examine all aspects of the current methodologies, including ligand design, characterization, and deployment. Finally we offer points of consideration and highlight the future outlook as more photoaffinity ligands emerge within the field.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory.

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

A predicted binding site for cholesterol on the GABAA receptor.

Modulation of the GABA type A receptor (GABAAR) function by cholesterol and other steroids is documented at the functional level, yet its structural basis is largely unknown. Current data on structurally related modulators suggest that cholesterol binds to subunit interfaces between transmembrane domains of the GABAAR. We construct homology models of a human GABAAR based on the structure of the glutamate-gated chloride channel GluCl of Caenorhabditis elegans. The models show the possibility of previously unreported disulfide bridges linking the M1 and M3 transmembrane helices in the α and γ subunits. We discuss the biological relevance of such disulfide bridges. Using our models, we investigate cholesterol binding to intersubunit cavities of the GABAAR transmembrane domain. We find that very similar binding modes are predicted independently by three approaches: analogy with ivermectin in the GluCl crystal structure, automated docking by AutoDock, and spontaneous rebinding events in unbiased molecular dynamics simulations. Taken together, the models and atomistic simulations suggest a somewhat flexible binding mode, with several possible orientations. Finally, we explore the possibility that cholesterol promotes pore opening through a wedge mechanism.

Propofol directly binds and inhibits skeletal muscle ryanodine receptor 1 (RyR1).

As the primary Ca 2+ release channel in skeletal muscle sarcoplasmic reticulum (SR), mutations in the type 1 ryanodine receptor (RyR1) or its binding partners underlie a constellation of muscle disorders, including malignant hyperthermia (MH). In patients with MH mutations, exposure to triggering drugs such as the halogenated volatile anesthetics biases RyR1 to an open state, resulting in uncontrolled Ca 2+ release, sarcomere tension and heat production. Restoration of Ca 2+ into the SR also consumes ATP, generating a further untenable metabolic load. When anesthetizing patients with known MH mutations, the non-triggering intravenous general anesthetic propofol is commonly substituted for triggering anesthetics. Evidence of direct binding of anesthetic agents to RyR1 or its binding partners is scant, and the atomic-level interactions of propofol with RyR1 are entirely unknown. Here, we show that propofol decreases RyR1 opening in heavy SR vesicles and planar lipid bilayers, and that it inhibits activator-induced Ca 2+ release from SR in human skeletal muscle. In addition to confirming direct binding, photoaffinity labeling using m- azipropofol (AziP m ) revealed several putative propofol binding sites on RyR1. Prediction of binding affinity by molecular dynamics simulation suggests that propofol binds at least one of these sites at clinical concentrations. These findings invite the hypothesis that in addition to propofol not triggering MH, it may also be protective against MH by inhibiting induced Ca 2+ flux through RyR1.

The blobulator: a webtool for identification and visual exploration of hydrophobic modularity in protein sequences.

Motivation: Clusters of hydrophobic residues are known to promote structured protein stability and drive protein aggregation. Recent work has shown that identifying contiguous hydrophobic residue clusters (termed "blobs") has proven useful in both intrinsically disordered protein (IDP) simulation and human genome studies. However, a graphical interface was unavailable. Results: Here, we present the blobulator: an interactive and intuitive web interface to detect intrinsic modularity in any protein sequence based on hydrophobicity. We demonstrate three use cases of the blobulator and show how identifying blobs with biologically relevant parameters provides useful information about a globular protein, two orthologous membrane proteins, and an IDP. Other potential applications are discussed, including: predicting protein segments with critical roles in tertiary interactions, providing a definition of local order and disorder with clear edges, and aiding in predicting protein features from sequence. Availability: The blobulator GUI can be found at www.blobulator.branniganlab.org, and the source code with pip installable command line tool can be found on GitHub at www.GitHub.com/BranniganLab/blobulator.

Lipid nanodisc scaffold and size alter the structure of a pentameric ligand-gated ion channel.

Lipid nanodiscs have become a standard tool for studying membrane proteins, including using single particle cryo-electron microscopy (cryo-EM). We find that reconstituting the pentameric ligand-gated ion channel (pLGIC), Erwinia ligand-gated ion channel (ELIC), in different nanodiscs produces distinct structures by cryo-EM. The effect of the nanodisc on ELIC structure extends to the extracellular domain and agonist binding site. Additionally, molecular dynamic simulations indicate that nanodiscs of different size impact ELIC structure and that the nanodisc scaffold directly interacts with ELIC. These findings suggest that the nanodisc plays a crucial role in determining the structure of pLGICs, and that reconstitution of ion channels in larger nanodiscs may better approximate a lipid membrane environment.

General anesthetics predicted to block the GLIC pore with micromolar affinity.

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and 10 µM for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore.

Multiple binding sites for the general anesthetic isoflurane identified in the nicotinic acetylcholine receptor transmembrane domain.

An extensive search for isoflurane binding sites in the nicotinic acetylcholine receptor (nAChR) and the proton gated ion channel from Gloebacter violaceus (GLIC) has been carried out based on molecular dynamics (MD) simulations in fully hydrated lipid membrane environments. Isoflurane introduced into the aqueous phase readily partitions into the lipid membrane and the membrane-bound protein. Specifically, isoflurane binds persistently to three classes of sites in the nAChR transmembrane domain: (i) An isoflurane dimer occludes the pore, contacting residues identified by previous mutagenesis studies; analogous behavior is observed in GLIC. (ii) Several nAChR subunit interfaces are also occupied, in a site suggested by photoaffinity labeling and thought to positively modulate the receptor; these sites are not occupied in GLIC. (iii) Isoflurane binds to the subunit centers of both nAChR alpha chains and one of the GLIC chains, in a site that has had little experimental targeting. Interpreted in the context of existing structural and physiological data, the present MD results support a multisite model for the mechanism of receptor-channel modulation by anesthetics.

Alchemical Free Energy Calculations on Membrane-Associated Proteins.

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Identification of a fluorescent general anesthetic, 1-aminoanthracene.

We identified a fluorophore, 1-aminoanthracene (1-AMA), that is anesthetic, potentiates GABAergic transmission, and gives an appropriate dissociation constant, K(d) approximately 0.1 mM, for binding to the general anesthetic site in horse spleen apoferritin (HSAF). 1-AMA fluorescence is enhanced when bound to HSAF. Thus, displacement of 1-AMA from HSAF by other anesthetics attenuates the fluorescence signal and allows determination of K(d), as validated by isothermal titration calorimetry. This provides a unique fluorescence assay for compound screening and anesthetic discovery. Additional electrophysiology experiments in isolated cells indicate that 1-AMA potentiates chloride currents elicited by GABA, similar to many general anesthetics. Furthermore, 1-AMA reversibly immobilizes stage 45-50 Xenopus laevis tadpoles (EC(50) = 16 microM) and fluorescence micrographs show 1-AMA localized to brain and olfactory regions. Thus, 1-AMA provides an unprecedented opportunity for studying general anesthetic distribution in vivo at the cellular and subcellular levels.

Polyunsaturated fatty acids inhibit a pentameric ligand-gated ion channel through one of two binding sites.

Polyunsaturated fatty acids (PUFAs) inhibit pentameric ligand-gated ion channels (pLGICs) but the mechanism of inhibition is not well understood. The PUFA, docosahexaenoic acid (DHA), inhibits agonist responses of the pLGIC, ELIC, more effectively than palmitic acid, similar to the effects observed in the GABAA receptor and nicotinic acetylcholine receptor. Using photo-affinity labeling and coarse-grained molecular dynamics simulations, we identified two fatty acid binding sites in the outer transmembrane domain (TMD) of ELIC. Fatty acid binding to the photolabeled sites is selective for DHA over palmitic acid, and specific for an agonist-bound state. Hexadecyl-methanethiosulfonate modification of one of the two fatty acid binding sites in the outer TMD recapitulates the inhibitory effect of PUFAs in ELIC. The results demonstrate that DHA selectively binds to multiple sites in the outer TMD of ELIC, but that state-dependent binding to a single intrasubunit site mediates DHA inhibition of ELIC.

Open-channel structure of a pentameric ligand-gated ion channel reveals a mechanism of leaflet-specific phospholipid modulation.

Pentameric ligand-gated ion channels (pLGICs) mediate synaptic transmission and are sensitive to their lipid environment. The mechanism of phospholipid modulation of any pLGIC is not well understood. We demonstrate that the model pLGIC, ELIC (Erwinia ligand-gated ion channel), is positively modulated by the anionic phospholipid, phosphatidylglycerol, from the outer leaflet of the membrane. To explore the mechanism of phosphatidylglycerol modulation, we determine a structure of ELIC in an open-channel conformation. The structure shows a bound phospholipid in an outer leaflet site, and structural changes in the phospholipid binding site unique to the open-channel. In combination with streamlined alchemical free energy perturbation calculations and functional measurements in asymmetric liposomes, the data support a mechanism by which an anionic phospholipid stabilizes the activated, open-channel state of a pLGIC by specific, state-dependent binding to this site.

Embedded cholesterol in the nicotinic acetylcholine receptor.

The nicotinic acetylcholine receptor (nAChR) is a cation-selective channel central to both neuronal and muscular processes and is considered the prototype for ligand-gated ion channels, motivating a structural determination effort that spanned several decades [Unwin N (2005) Refined structure of the nicotinic acetylcholine receptor at 4 A resolution. J Mol Biol 346:967-989]. Purified nAChR must be reconstituted in a mixture containing cholesterol to function. Proposed modes of interaction between cholesterol and the protein range from specific binding to indirect membrane-mediated mechanisms. However, the underlying cause of nAChR sensitivity to cholesterol remains controversial, in part because the vast majority of functional studies were conducted before a medium resolution structure was reported. We show that the nAChR contains internal sites capable of containing cholesterol, whose occupation stabilizes the protein structure. We detect sites at the protein-lipid interface as conventionally predicted from functional data, as well as deeply buried sites that are not usually considered. Molecular dynamics simulations reveal that occupation of both superficial and deeply buried sites most effectively preserves the experimental structure; the structure collapses in the absence of bound cholesterol. In particular, we find that bound cholesterol directly supports contacts between the agonist-binding domain and the pore that are thought to be essential for activation of the receptor. These results likely apply to those other ion channels within the Cys-loop superfamily that depend on cholesterol, such as the GABA receptor.