Rosace: a robust deep mutational scanning analysis framework employing position and mean-variance shrinkage.

Deep mutational scanning (DMS) measures the effects of thousands of genetic variants in a protein simultaneously. The small sample size renders classical statistical methods ineffective. For example, p-values cannot be correctly calibrated when treating variants independently. We propose Rosace, a Bayesian framework for analyzing growth-based DMS data. Rosace leverages amino acid position information to increase power and control the false discovery rate by sharing information across parameters via shrinkage. We also developed Rosette for simulating the distributional properties of DMS. We show that Rosace is robust to the violation of model assumptions and is more powerful than existing tools.

The full spectrum of SLC22 OCT1 mutations illuminates the bridge between drug transporter biophysics and pharmacogenomics.

Mutations in transporters can impact an individual's response to drugs and cause many diseases. Few variants in transporters have been evaluated for their functional impact. Here, we combine saturation mutagenesis and multi-phenotypic screening to dissect the impact of 11,213 missense single-amino-acid deletions, and synonymous variants across the 554 residues of OCT1, a key liver xenobiotic transporter. By quantifying in parallel expression and substrate uptake, we find that most variants exert their primary effect on protein abundance, a phenotype not commonly measured alongside function. Using our mutagenesis results combined with structure prediction and molecular dynamic simulations, we develop accurate structure-function models of the entire transport cycle, providing biophysical characterization of all known and possible human OCT1 polymorphisms. This work provides a complete functional map of OCT1 variants along with a framework for integrating functional genomics, biophysical modeling, and human genetics to predict variant effects on disease and drug efficacy.

Molecular basis of proton-sensing by G protein-coupled receptors.

Three proton-sensing G protein-coupled receptors (GPCRs), GPR4, GPR65, and GPR68, respond to changes in extracellular pH to regulate diverse physiology and are implicated in a wide range of diseases. A central challenge in determining how protons activate these receptors is identifying the set of residues that bind protons. Here, we determine structures of each receptor to understand the spatial arrangement of putative proton sensing residues in the active state. With a newly developed deep mutational scanning approach, we determined the functional importance of every residue in proton activation for GPR68 by generating ~9,500 mutants and measuring effects on signaling and surface expression. This unbiased screen revealed that, unlike other proton-sensitive cell surface channels and receptors, no single site is critical for proton recognition in GPR68. Instead, a network of titratable residues extend from the extracellular surface to the transmembrane region and converge on canonical class A GPCR activation motifs to activate proton-sensing GPCRs. More broadly, our approach integrating structure and unbiased functional interrogation defines a new framework for understanding the rich complexity of GPCR signaling.

Conserved regulatory motifs in the juxtamembrane domain and kinase N-lobe revealed through deep mutational scanning of the MET receptor tyrosine kinase domain.

MET is a receptor tyrosine kinase (RTK) responsible for initiating signaling pathways involved in development and wound repair. MET activation relies on ligand binding to the extracellular receptor, which prompts dimerization, intracellular phosphorylation, and recruitment of associated signaling proteins. Mutations, which are predominantly observed clinically in the intracellular juxtamembrane and kinase domains, can disrupt typical MET regulatory mechanisms. Understanding how juxtamembrane variants, such as exon 14 skipping (METΔEx14), and rare kinase domain mutations can increase signaling, often leading to cancer, remains a challenge. Here, we perform a parallel deep mutational scan (DMS) of MET intracellular kinase domain in two fusion protein backgrounds: wild type and METΔEx14. Our comparative approach has revealed a critical hydrophobic interaction between a juxtamembrane segment and the kinase αC helix, pointing to differences in regulatory mechanisms between MET and other RTKs. Additionally, we have uncovered a ß5 motif that acts as a structural pivot for kinase domain activation in MET and other TAM family of kinases. We also describe a number of previously unknown activating mutations, aiding the effort to annotate driver, passenger, and drug resistance mutations in the MET kinase domain.

The full spectrum of OCT1 (SLC22A1) mutations bridges transporter biophysics to drug pharmacogenomics.

Membrane transporters play a fundamental role in the tissue distribution of endogenous compounds and xenobiotics and are major determinants of efficacy and side effects profiles. Polymorphisms within these drug transporters result in inter-individual variation in drug response, with some patients not responding to the recommended dosage of drug whereas others experience catastrophic side effects. For example, variants within the major hepatic Human organic cation transporter OCT1 (SLC22A1) can change endogenous organic cations and many prescription drug levels. To understand how variants mechanistically impact drug uptake, we systematically study how all known and possible single missense and single amino acid deletion variants impact expression and substrate uptake of OCT1. We find that human variants primarily disrupt function via folding rather than substrate uptake. Our study revealed that the major determinants of folding reside in the first 300 amino acids, including the first 6 transmembrane domains and the extracellular domain (ECD) with a stabilizing and highly conserved stabilizing helical motif making key interactions between the ECD and transmembrane domains. Using the functional data combined with computational approaches, we determine and validate a structure-function model of OCT1s conformational ensemble without experimental structures. Using this model and molecular dynamic simulations of key mutants, we determine biophysical mechanisms for how specific human variants alter transport phenotypes. We identify differences in frequencies of reduced function alleles across populations with East Asians vs European populations having the lowest and highest frequency of reduced function variants, respectively. Mining human population databases reveals that reduced function alleles of OCT1 identified in this study associate significantly with high LDL cholesterol levels. Our general approach broadly applied could transform the landscape of precision medicine by producing a mechanistic basis for understanding the effects of human mutations on disease and drug response.

DIMPLE: deep insertion, deletion, and missense mutation libraries for exploring protein variation in evolution, disease, and biology.

Insertions and deletions (indels) enable evolution and cause disease. Due to technical challenges, indels are left out of most mutational scans, limiting our understanding of them in disease, biology, and evolution. We develop a low cost and bias method, DIMPLE, for systematically generating deletions, insertions, and missense mutations in genes, which we test on a range of targets, including Kir2.1. We use DIMPLE to study how indels impact potassium channel structure, disease, and evolution. We find deletions are most disruptive overall, beta sheets are most sensitive to indels, and flexible loops are sensitive to deletions yet tolerate insertions.

Still rocking in the structural era: A molecular overview of the small multidrug resistance (SMR) transporter family.

The small multidrug resistance (SMR) family is composed of widespread microbial membrane proteins that fulfill different transport functions. Four functional SMR subtypes have been identified, which variously transport the small, charged metabolite guanidinium, bulky hydrophobic drugs and antiseptics, polyamines, and glycolipids across the membrane bilayer. The transporters possess a minimalist architecture, with ∼100-residue subunits that require assembly into homodimers or heterodimers for transport. In part because of their simple construction, the SMRs are a tractable system for biochemical and biophysical analysis. Studies of SMR transporters over the last 25 years have yielded deep insights for diverse fields, including membrane protein topology and evolution, mechanisms of membrane transport, and bacterial multidrug resistance. Here, we review recent advances in understanding the structures and functions of SMR transporters. New molecular structures of SMRs representing two of the four functional subtypes reveal the conserved structural features that have permitted the emergence of disparate substrate transport functions in the SMR family and illuminate structural similarities with a distantly related membrane transporter family, SLC35/DMT.

The structural basis of promiscuity in small multidrug resistance transporters.

By providing broad resistance to environmental biocides, transporters from the small multidrug resistance (SMR) family drive the spread of multidrug resistance cassettes among bacterial populations. A fundamental understanding of substrate selectivity by SMR transporters is needed to identify the types of selective pressures that contribute to this process. Using solid-supported membrane electrophysiology, we find that promiscuous transport of hydrophobic substituted cations is a general feature of SMR transporters. To understand the molecular basis for promiscuity, we solved X-ray crystal structures of a SMR transporter Gdx-Clo in complex with substrates to a maximum resolution of 2.3 Å. These structures confirm the family's extremely rare dual topology architecture and reveal a cleft between two helices that provides accommodation in the membrane for the hydrophobic substituents of transported drug-like cations.

Guanidinium export is the primal function of SMR family transporters.

The small multidrug resistance (SMR) family of membrane proteins is prominent because of its rare dual topology architecture, simplicity, and small size. Its best studied member, EmrE, is an important model system in several fields related to membrane protein biology, from evolution to mechanism. But despite decades of work on these multidrug transporters, the native function of the SMR family has remained a mystery, and many highly similar SMR homologs do not transport drugs at all. Here we establish that representative SMR proteins, selected from each of the major clades in the phylogeny, function as guanidinium ion exporters. Drug-exporting SMRs are all clustered in a single minority clade. Using membrane transport experiments, we show that these guanidinium exporters, which we term Gdx, are very selective for guanidinium and strictly and stoichiometrically couple its export with the import of two protons. These findings draw important mechanistic distinctions with the notably promiscuous and weakly coupled drug exporters like EmrE.

A topologically diverse family of fluoride channels.

Dual-topology proteins are likely evolutionary antecedents to a common motif in membrane protein structures, the inverted repeat. A family of fluoride channels, the Flucs, which protect microorganisms, fungi, and plants against cytoplasmic fluoride accumulation, has representatives of all topologies along this evolutionary trajectory, including dual-topology homodimers, antiparallel heterodimers, and, in eukaryotes, fused two-domain proteins with an inverted repeat motif. Recent high-resolution crystal structures of dual-topology homodimers, coupled with extensive functional information about both the homodimers and two-domain Flucs, provide a case study of the co-evolution of fold and function.

Antineoplastic agents. 599. Total synthesis of dolastatin 16.

The first 23-step total synthesis of the cyclodepsipeptide dolastatin 16 (1) has been achieved. Synthesis of the dolaphenvaline and dolamethylleuine amino acid units using simplified methods improved the overall efficiency. The formation of the 25-membered macrocycle employing lactonization with 2-methyl-6-nitrobenzoic anhydride completed a key step in the synthesis. Regrettably, the synthetic dolastatin 16 (1), while otherwise identical (by X-ray crystal structure and spectral analyses) with the natural product, did not reproduce the powerful (nanomolar) cancer cell growth inhibition displayed by the natural isolate. Presumably this result can be attributed to conformation(s) of the synthetic dolastatin 16 (1) or to a chemically undetected component isolated with the natural product.

Antineoplastic agents. 600. From the South Pacific Ocean to the silstatins.

The recent advances in the development of antibody and other drug conjugates for targeted cancer treatment have further increased the need for powerful cancer cell growth inhibitors. Toward that objective we have extended our earlier discovery of the remarkable anticancer bacillistatins 1 and 2 from Bacillus silvestris to SAR and other structural modifications such as availability of a free hydroxy group for antibody-drug conjugate (ADC) and other prodrug linkage. That direction has resulted in seven structural modifications designated silstatins 1-8 (7a, 8a, 8b, 14a, 15a, 15b, 18a, and 18b), where the exceptional cancer cell growth inhibition of some of them are in the range GI50 10(-3)-10(-4) µM/mL. Silstatin 7 (18a) was converted to a glucuronic conjugate (28) that displayed an impressive reduction in toxicity during transport.