Molecular basis for variations in the sensitivity of pathogenic rhodopsin variants to 9-cis-retinal.

Over 100 mutations in the rhodopsin gene have been linked to a spectrum of retinopathies that include retinitis pigmentosa and congenital stationary night blindness. Though most of these variants exhibit a loss of function, the molecular defects caused by these underlying mutations vary considerably. In this work, we utilize deep mutational scanning to quantitatively compare the plasma membrane expression of 123 known pathogenic rhodopsin variants in the presence and absence of the stabilizing cofactor 9-cis-retinal. We identify 69 retinopathy variants, including 20 previously uncharacterized variants, that exhibit diminished plasma membrane expression in HEK293T cells. Of these apparent class II variants, 67 exhibit a measurable increase in expression in the presence of 9-cis-retinal. However, the magnitude of the response to this molecule varies considerably across this spectrum of mutations. Evaluation of the observed shifts relative to thermodynamic estimates for the coupling between binding and folding suggests underlying differences in stability constrains the magnitude of their response to retinal. Nevertheless, estimates from computational modeling suggest that many of the least sensitive variants also directly compromise binding. Finally, we evaluate the functional properties of three previous uncharacterized, retinal-sensitive variants (ΔN73, S131P, and R135G) and show that two of these retain residual function in vitro. Together, our results provide a comprehensive experimental characterization of the proteostatic properties of retinopathy variants and their response to retinal.

Coordination of -1 programmed ribosomal frameshifting by transcript and nascent chain features revealed by deep mutational scanning.

Programmed ribosomal frameshifting (PRF) is a translational recoding mechanism that enables the synthesis of multiple polypeptides from a single transcript. During translation of the alphavirus structural polyprotein, the efficiency of -1PRF is coordinated by a 'slippery' sequence in the transcript, an adjacent RNA stem-loop, and a conformational transition in the nascent polypeptide chain. To characterize each of these effectors, we measured the effects of 4530 mutations on -1PRF by deep mutational scanning. While most mutations within the slip-site and stem-loop reduce the efficiency of -1PRF, the effects of mutations upstream of the slip-site are far more variable. We identify several regions where modifications of the amino acid sequence of the nascent polypeptide impact the efficiency of -1PRF. Molecular dynamics simulations of polyprotein biogenesis suggest the effects of these mutations primarily arise from their impacts on the mechanical forces that are generated by the translocon-mediated cotranslational folding of the nascent polypeptide chain. Finally, we provide evidence suggesting that the coupling between cotranslational folding and -1PRF depends on the translation kinetics upstream of the slip-site. These findings demonstrate how -1PRF is coordinated by features within both the transcript and nascent chain.

Towards generalizable predictions for G protein-coupled receptor variant expression.

Missense mutations that compromise the plasma membrane expression (PME) of integral membrane proteins are the root cause of numerous genetic diseases. Differentiation of this class of mutations from those that specifically modify the activity of the folded protein has proven useful for the development and targeting of precision therapeutics. Nevertheless, it remains challenging to predict the effects of mutations on the stability and/ or expression of membrane proteins. In this work, we utilize deep mutational scanning data to train a series of artificial neural networks to predict the PME of transmembrane domain variants of G protein-coupled receptors from structural and/ or evolutionary features. We show that our best-performing network, which we term the PME predictor, can recapitulate mutagenic trends within rhodopsin and can differentiate pathogenic transmembrane domain variants that cause it to misfold from those that compromise its signaling. This network also generates statistically significant predictions for the relative PME of transmembrane domain variants for another class A G protein-coupled receptor (ß2 adrenergic receptor) but not for an unrelated voltage-gated potassium channel (KCNQ1). Notably, our analyses of these networks suggest structural features alone are generally sufficient to recapitulate the observed mutagenic trends. Moreover, our findings imply that networks trained in this manner may be generalizable to proteins that share a common fold. Implications of our findings for the design of mechanistically specific genetic predictors are discussed.

Systematic profiling of temperature- and retinal-sensitive rhodopsin variants by deep mutational scanning.

Membrane protein variants with diminished conformational stability often exhibit enhanced cellular expression at reduced growth temperatures. The expression of "temperature-sensitive" variants is also typically sensitive to corrector molecules that bind and stabilize the native conformation. There are many examples of temperature-sensitive rhodopsin variants, the misfolding of which is associated with the molecular basis of retinitis pigmentosa. In this work, we employ deep mutational scanning to compare the effects of reduced growth temperature and 9-cis-retinal, an investigational corrector, on the plasma membrane expression of 700 rhodopsin variants in HEK293T cells. We find that the change in expression at reduced growth temperatures correlates with the response to 9-cis-retinal among variants bearing mutations within a hydrophobic transmembrane domain (TM2). The most sensitive variants appear to disrupt a native helical kink within this transmembrane domain. By comparison, mutants that alter the structure of a polar transmembrane domain (TM7) exhibit weaker responses to temperature and retinal that are poorly correlated. Statistical analyses suggest that this observed insensitivity cannot be attributed to a single variable, but likely arises from the composite effects of mutations on the energetics of membrane integration, the stability of the native conformation, and the integrity of the retinal-binding pocket. Finally, we show that the characteristics of purified temperature- and retinal-sensitive variants suggest that the proteostatic effects of retinal may be manifested during translation and cotranslational folding. Together, our findings highlight several biophysical constraints that appear to influence the sensitivity of genetic variants to temperature and small-molecule correctors.

Classification of the Molecular Defects Associated with Pathogenic Variants of the SLC6A8 Creatine Transporter.

More than 80 loss-of-function (LOF) mutations in the SLC6A8 creatine transporter (hCRT1) are responsible for cerebral creatine deficiency syndrome (CCDS), which gives rise to a spectrum of neurological defects, including intellectual disability, epilepsy, and autism spectrum disorder. To gain insight into the nature of the molecular defects caused by these mutations, we quantitatively profiled the cellular processing, trafficking, expression, and function of eight pathogenic CCDS variants in relation to the wild type (WT) and one neutral isoform. All eight CCDS variants exhibit measurable proteostatic deficiencies that likely contribute to the observed LOF. However, the magnitudes of their specific effects on the expression and trafficking of hCRT1 vary considerably, and we find that the LOF associated with two of these variants primarily arises from the disruption of the substrate-binding pocket. In conjunction with an analysis of structural models of the transporter, we use these data to suggest mechanistic classifications for these variants. To evaluate potential avenues for therapeutic intervention, we assessed the sensitivity of these variants to temperature and measured their response to the proteostasis regulator 4-phenylbutyrate (4-PBA). Only one of the tested variants (G132V) is sensitive to temperature, though its response to 4-PBA is negligible. Nevertheless, 4-PBA significantly enhances the activity of WT hCRT1 in HEK293T cells, which suggests it may be worth evaluating as a therapeutic for female intellectual disability patients carrying a single CCDS mutation. Together, these findings reveal that pathogenic SLC6A8 mutations cause a spectrum of molecular defects that should be taken into consideration in future efforts to develop CCDS therapeutics.

Divergent Folding-Mediated Epistasis Among Unstable Membrane Protein Variants.

Many membrane proteins are prone to misfolding, which compromises their functional expression at the plasma membrane. This is particularly true for the mammalian gonadotropin-releasing hormone receptor GPCRs (GnRHR). We recently demonstrated that evolutionary GnRHR modifications appear to have coincided with adaptive changes in cotranslational folding efficiency. Though protein stability is known to shape evolution, it is unclear how cotranslational folding constraints modulate the synergistic, epistatic interactions between mutations. We therefore compared the pairwise interactions formed by mutations that disrupt the membrane topology (V276T) or tertiary structure (W107A) of GnRHR. Using deep mutational scanning, we evaluated how the plasma membrane expression of these variants is modified by hundreds of secondary mutations. An analysis of 251 mutants in three genetic backgrounds reveals that V276T and W107A form distinct epistatic interactions that depend on both the severity and the mechanism of destabilization. V276T forms predominantly negative epistatic interactions with destabilizing mutations in soluble loops. In contrast, W107A forms positive interactions with mutations in both loops and transmembrane domains that reflect the diminishing impacts of the destabilizing mutations in variants that are already unstable. These findings reveal how epistasis is remodeled by conformational defects in membrane proteins and in unstable proteins more generally.

Structure-Activity Relationship Study of Cannabidiol-Based Analogs as Negative Allosteric Modulators of the µ-Opioid Receptor.

The US faces an unprecedented surge in fatal drug overdoses. Naloxone, the only antidote for opiate overdose, competes at the mu opioid receptor (µOR) orthosteric site. Naloxone struggles against fentanyl-class synthetic opioids that now cause ∼80% of deaths. Negative allosteric modulators (NAMs) targeting secondary sites may noncompetitively downregulate µOR activation. (-)-Cannabidiol ((-)-CBD) is a candidate µOR NAM. To explore its therapeutic potential, we evaluated the structure-activity relationships among CBD analogs to identify NAMs with increased potency. Using a cyclic AMP assay, we characterize reversal of µOR activation by 15 CBD analogs, several of which proved more potent than (-)-CBD. Comparative docking investigations suggest that potent compounds interact with a putative allosteric pocket to stabilize the inactive µOR conformation. Finally, these compounds enhance naloxone displacement of fentanyl from the orthosteric site. Our results suggest that CBD analogs offer considerable potential for the development of next-generation antidotes for opioid overdose.

General trends in the effects of VX-661 and VX-445 on the plasma membrane expression of clinical CFTR variants.

Cystic fibrosis (CF) is caused by mutations that compromise the expression and/or function of the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel. Most people with CF harbor a common misfolded variant (ΔF508) that can be partially rescued by therapeutic "correctors" that restore its expression. Nevertheless, many other CF variants are insensitive to correctors. Using deep mutational scanning, we quantitatively compare the effects of two correctors on the plasma membrane expression of 129 CF variants. Though structural calculations suggest corrector binding provides similar stabilization to most variants, it's those with intermediate expression and mutations near corrector binding pockets that exhibit the greatest response. Deviations in sensitivity appear to depend on the degree of variant destabilization and the timing of misassembly. Combining correctors appears to rescue more variants by doubling the binding energy and stabilizing distinct cotranslational folding transitions. These results provide an overview of rare CF variant expression and establish new tools for precision pharmacology.

Structural analysis of the extracellular entrance to the serotonin transporter permeation pathway.

Neurotransmitter transporters are responsible for removal of biogenic amine neurotransmitters after release into the synapse. These transporters are the targets for many clinically relevant drugs, such as antidepressants and psychostimulants. A high resolution crystal structure for the monoamine transporters has yet to be solved. We have developed a homology model for the serotonin transporter (SERT) based on the crystal structure of the leucine transporter (LeuT(Aa)) from Aquifex aeolicus. The objective of the present studies is to identify the structural determinants forming the entrance to the substrate permeation pathway based on predictions from the SERT homology model. Using the substituted cysteine accessibility method, we identified residues predicted to reside at the entrance to the substrate permeation pathway that were reactive with methanethiosulfonate (MTS) reagents. Of these residues, Gln(332) in transmembrane helix (TMH) VI was protected against MTS inactivation in the presence of serotonin. Surprisingly, the reactivity of Gln(332) to MTS reagents was enhanced in the presence of cocaine. Bifunctional MTS cross-linkers also were used to examine the distances between helices predicted to form the entrance into the substrate and ion permeation pathway. Our studies suggest that substrate and ligand binding may induce conformational shifts in TMH I and/or VI, providing new opportunities to refine existing homology models of SERT and related monoamine transporters.

Molecular basis for the evolved instability of a human G-protein coupled receptor.

Membrane proteins are prone to misfolding and degradation. This is particularly true for mammalian forms of the gonadotropin-releasing hormone receptor (GnRHR). Although they function at the plasma membrane, mammalian GnRHRs accumulate within the secretory pathway. Their apparent instability is believed to have evolved through selection for attenuated GnRHR activity. Nevertheless, the molecular basis of this adaptation remains unclear. We show that adaptation coincides with a C-terminal truncation that compromises the translocon-mediated membrane integration of its seventh transmembrane domain (TM7). We also identify a series of polar residues in mammalian GnRHRs that compromise the membrane integration of TM2 and TM6. Reverting a lipid-exposed polar residue in TM6 to an ancestral hydrophobic residue restores expression with no impact on function. Evolutionary trends suggest variations in the polarity of this residue track with reproductive phenotypes. Our findings suggest that the marginal energetics of cotranslational folding can be exploited to tune membrane protein fitness.

Probing biophysical sequence constraints within the transmembrane domains of rhodopsin by deep mutational scanning.

Membrane proteins must balance the sequence constraints associated with folding and function against the hydrophobicity required for solvation within the bilayer. We recently found the expression and maturation of rhodopsin are limited by the hydrophobicity of its seventh transmembrane domain (TM7), which contains polar residues that are essential for function. On the basis of these observations, we hypothesized that rhodopsin's expression should be less tolerant of mutations in TM7 relative to those within hydrophobic TM domains. To test this hypothesis, we used deep mutational scanning to compare the effects of 808 missense mutations on the plasma membrane expression of rhodopsin in HEK293T cells. Our results confirm that a higher proportion of mutations within TM7 (37%) decrease rhodopsin's plasma membrane expression relative to those within a hydrophobic TM domain (TM2, 25%). These results in conjunction with an evolutionary analysis suggest solvation energetics likely restricts the evolutionary sequence space of polar TM domains.

The external gate of the human and Drosophila serotonin transporters requires a basic/acidic amino acid pair for 3,4-methylenedioxymethamphetamine (MDMA) translocation and the induction of substrate efflux.

The substituted amphetamine, 3,4-methylenedioxy-methamphetamine (MDMA, ecstasy), is a widely used drug of abuse that induces non-exocytotic release of serotonin, dopamine, and norepinephrine through their cognate transporters as well as blocking the reuptake of neurotransmitter by the same transporters. The resulting dramatic increase in volume transmission and signal duration of neurotransmitters leads to psychotropic, stimulant, and entactogenic effects. The mechanism by which amphetamines drive reverse transport of the monoamines remains largely enigmatic, however, promising outcomes for the therapeutic utility of MDMA for post-traumatic stress disorder and the long-time use of the dopaminergic and noradrenergic-directed amphetamines in treatment of attention-deficit hyperactivity disorder and narcolepsy increases the importance of understanding this phenomenon. Previously, we identified functional differences between the human and Drosophila melanogaster serotonin transporters (hSERT and dSERT, respectively) revealing that MDMA is an effective substrate for hSERT but not dSERT even though serotonin is a potent substrate for both transporters. Chimeric dSERT/hSERT transporters revealed that the molecular components necessary for recognition of MDMA as a substrate was linked to regions of the protein flanking transmembrane domains (TM) V through IX. Here, we performed species-scanning mutagenesis of hSERT, dSERT and C. elegans SERT (ceSERT) along with biochemical and electrophysiological analysis and identified a single amino acid in TM10 (Glu394, hSERT; Asn484, dSERT, Asp517, ceSERT) that is primarily responsible for the differences in MDMA recognition. Our findings reveal that an acidic residue is necessary at this position for MDMA recognition as a substrate and serotonin releaser.

Conformational flexibility of transmembrane helix VII of the human serotonin transporter impacts ion dependence and transport.

The serotonin transporter (SERT) regulates the serotonin concentration in the synapse and is a target of several antidepressant and psychostimulant drugs. Previous work suggested that the middle transmembrane helices (TMHs) of the biogenic amine transporters (TMHs) play a role in substrate and ion recognition. We focused our present studies on exploring the role of TMH VII in transporter function and ion recognition. Residues divergent between human SERT and Drosophila SERT (hSERT and dSERT, respectively) were identified and mutated in hSERT to the corresponding identity in dSERT. hSERT mutants V366S, M370L, S375A, and T381S exhibited a decrease in transport capacity. To further explore the role of these residues in the transport process, we generated cysteine mutants at multiple positions. Pretreatment with [2-(trimethylammonium)ethyl] methanethiosulfonate (MTSET) caused a decrease in transport of [(3)H]5-HT in the V366C and M370C mutants. The hSERT V366S, M370L, and M370C mutations also altered the sodium and chloride dependence for substrate transport. Interpretation of our results in the context of a homology model of SERT based on the crystal structure of the Aquifex aeolicus leucine transporter suggests flexibility in the conformation of TMH VII that impacts ion dependence and substrate transport.

Mouse embryonic fibroblast cells from transgenic mice overexpressing tNOX exhibit an altered growth and drug response phenotype.

Mouse embryonic fibroblast (MEF) cells prepared from transgenic mice overexpressing a cancer-specific and growth-related cell surface NADH oxidase with protein disulfide-thiol interchange activity grew at rates approximately twice those of wild-type embryonic fibroblast cells. Growth of transgenic MEF cells overexpressing tNOX was inhibited by low concentrations of the green tea catechin (-)-epigallocatechin-3-gallate (EGCg) or the synthetic isoflavene phenoxodiol. Both are putative tNOX-targeted inhibitors with anti-cancer activity. With both EGCg and phenoxodiol, growth inhibition was followed after about 48 h by apoptosis. Growth of wild-type mouse fibroblast cells from the same strain was unaffected by EGCg and phenoxodiol and neither compound induced apoptosis even at concentrations 100-1,000-fold higher than those that resulted in apoptotic death in the transgenic MEF cells. The findings validate earlier reports of evidence for tNOX presence as contributing to unregulated growth of cancer cells as well as the previous identification of the tNOX protein as the molecular target for the anti-cancer activities attributed to both EGCg and phenoxodiol. The expression of tNOX emerges as both necessary and sufficient to account for the cancer cell-specific growth inhibitions by both EGCg and phenoxodiol.