Signaling Modulation Mediated by Ligand Water Interactions with the Sodium Site at µOR.

The mu opioid receptor (µOR) is a target for clinically used analgesics. However, adverse effects, such as respiratory depression and physical dependence, necessitate the development of alternative treatments. Recently we reported a novel strategy to design functionally selective opioids by targeting the sodium binding allosteric site in µOR with a supraspinally active analgesic named C6guano. Presently, to improve systemic activity of this ligand, we used structure-based design, identifying a new ligand named RO76 where the flexible alkyl linker and polar guanidine guano group is swapped with a benzyl alcohol, and the sodium site is targeted indirectly through waters. A cryoEM structure of RO76 bound to the µOR-Gi complex confirmed that RO76 interacts with the sodium site residues through a water molecule, unlike C6guano which engages the sodium site directly. Signaling assays coupled with APEX based proximity labeling show binding in the sodium pocket modulates receptor efficacy and trafficking. In mice, RO76 was systemically active in tail withdrawal assays and showed reduced liabilities compared to those of morphine. In summary, we show that targeting water molecules in the sodium binding pocket may be an avenue to modulate signaling properties of opioids, and which may potentially be extended to other G-protein coupled receptors where this site is conserved.

CryoEM structures of the human CLC-2 voltage-gated chloride channel reveal a ball-and-chain gating mechanism.

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different tissues. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating among closely related homologs has been a long-standing mystery, in part because few CLC channel structures are available. Here, we report cryoEM structures of human CLC-2 at 2.46 - 2.76 Å, in the presence and absence of the selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl--permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct conformations involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl--permeation pathway. This peptide is highly conserved among species variants of CLC-2 but is not present in other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a "ball-and-chain" gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl--permeation pathway.

Xanomeline displays concomitant orthosteric and allosteric binding modes at the M4 mAChR.

The M4 muscarinic acetylcholine receptor (M4 mAChR) has emerged as a drug target of high therapeutic interest due to its expression in regions of the brain involved in the regulation of psychosis, cognition, and addiction. The mAChR agonist, xanomeline, has provided significant improvement in the Positive and Negative Symptom Scale (PANSS) scores in a Phase II clinical trial for the treatment of patients suffering from schizophrenia. Here we report the active state cryo-EM structure of xanomeline bound to the human M4 mAChR in complex with the heterotrimeric Gi1 transducer protein. Unexpectedly, two molecules of xanomeline were found to concomitantly bind to the monomeric M4 mAChR, with one molecule bound in the orthosteric (acetylcholine-binding) site and a second molecule in an extracellular vestibular allosteric site. Molecular dynamic simulations supports the structural findings, and pharmacological validation confirmed that xanomeline acts as a dual orthosteric and allosteric ligand at the human M4 mAChR. These findings provide a basis for further understanding xanomeline's complex pharmacology and highlight the myriad of ways through which clinically relevant ligands can bind to and regulate GPCRs.

A non-canonical mechanism of GPCR activation.

The goal of designing safer, more effective drugs has led to tremendous interest in molecular mechanisms through which ligands can precisely manipulate signaling of G-protein-coupled receptors (GPCRs), the largest class of drug targets. Decades of research have led to the widely accepted view that all agonists-ligands that trigger GPCR activation-function by causing rearrangement of the GPCR's transmembrane helices, opening an intracellular pocket for binding of transducer proteins. Here we demonstrate that certain agonists instead trigger activation of free fatty acid receptor 1 by directly rearranging an intracellular loop that interacts with transducers. We validate the predictions of our atomic-level simulations by targeted mutagenesis; specific mutations which disrupt interactions with the intracellular loop convert these agonists into inverse agonists. Further analysis suggests that allosteric ligands could regulate signaling of many other GPCRs via a similar mechanism, offering rich possibilities for precise control of pharmaceutically important targets.

CryoEM structures of the human CLC-2 voltage gated chloride channel reveal a ball and chain gating mechanism.

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different mammalian tissues and cell types. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating mechanisms among closely related CLC homologs has been a long-standing mystery, in part because few CLC channel structures are available, and those that exist exhibit high conformational similarity. Here, we report cryoEM structures of human CLC-2 at 2.46 - 2.76 Å, in the presence and absence of the potent and selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl--permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct apo conformations of CLC-2 involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl--permeation pathway from the intracellular side. This peptide is highly conserved among species variants of CLC-2 but is not present in any other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a "ball-and-chain" gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we show that loss of this short sequence increases the magnitude and decreases the rectification of CLC-2 currents expressed in mammalian cells. Furthermore, we show that with repetitive hyperpolarization WT CLC-2 currents increase in resemblance to the hairpin-deleted CLC-2 currents. These functional results combined with our structural data support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl--permeation pathway.

Structural basis for activation of CB1 by an endocannabinoid analog.

Endocannabinoids (eCBs) are endogenous ligands of the cannabinoid receptor 1 (CB1), a G protein-coupled receptor that regulates a number of therapeutically relevant physiological responses. Hence, understanding the structural and functional consequences of eCB-CB1 interactions has important implications for designing effective drugs targeting this receptor. To characterize the molecular details of eCB interaction with CB1, we utilized AMG315, an analog of the eCB anandamide to determine the structure of the AMG315-bound CB1 signaling complex. Compared to previous structures, the ligand binding pocket shows some differences. Using docking, molecular dynamics simulations, and signaling assays we investigated the functional consequences of ligand interactions with the "toggle switch" residues F2003.36 and W3566.48. Further, we show that ligand-TM2 interactions drive changes to residues on the intracellular side of TM2 and are a determinant of efficacy in activating G protein. These intracellular TM2 rearrangements are unique to CB1 and are exploited by a CB1-specific allosteric modulator.

Structural basis of efficacy-driven ligand selectivity at GPCRs.

A drug's selectivity for target receptors is essential to its therapeutic utility, but achieving selectivity between similar receptors is challenging. The serendipitous discovery of ligands that stimulate target receptors more strongly than closely related receptors, despite binding with similar affinities, suggests a solution. The molecular mechanism of such 'efficacy-driven selectivity' has remained unclear, however, hindering design of such ligands. Here, using atomic-level simulations, we reveal the structural basis for the efficacy-driven selectivity of a long-studied clinical drug candidate, xanomeline, between closely related muscarinic acetylcholine receptors (mAChRs). Xanomeline's binding mode is similar across mAChRs in their inactive states but differs between mAChRs in their active states, with divergent effects on active-state stability. We validate this mechanism experimentally and use it to design ligands with altered efficacy-driven selectivity. Our results suggest strategies for the rational design of ligands that achieve efficacy-driven selectivity for many pharmaceutically important G-protein-coupled receptors.

Geometric Deep Learning for Structure-Based Ligand Design.

A pervasive challenge in drug design is determining how to expand a ligand-a small molecule that binds to a target biomolecule-in order to improve various properties of the ligand. Adding single chemical groups, known as fragments, is important for lead optimization tasks, and adding multiple fragments is critical for fragment-based drug design. We have developed a comprehensive framework that uses machine learning and three-dimensional protein-ligand structures to address this challenge. Our method, FRAME, iteratively determines where on a ligand to add fragments, selects fragments to add, and predicts the geometry of the added fragments. On a comprehensive benchmark, FRAME consistently improves predicted affinity and selectivity relative to the initial ligand, while generating molecules with more drug-like chemical properties than docking-based methods currently in widespread use. FRAME learns to accurately describe molecular interactions despite being given no prior information on such interactions. The resulting framework for quality molecular hypothesis generation can be easily incorporated into the workflows of medicinal chemists for diverse tasks, including lead optimization, fragment-based drug discovery, and de novo drug design.

Cryo-EM, Protein Engineering, and Simulation Enable the Development of Peptide Therapeutics against Acute Myeloid Leukemia.

Cryogenic electron microscopy (cryo-EM) has emerged as a viable structural tool for molecular therapeutics development against human diseases. However, it remains a challenge to determine structures of proteins that are flexible and smaller than 30 kDa. The 11 kDa KIX domain of CREB-binding protein (CBP), a potential therapeutic target for acute myeloid leukemia and other cancers, is a protein which has defied structure-based inhibitor design. Here, we develop an experimental approach to overcome the size limitation by engineering a protein double-shell to sandwich the KIX domain between apoferritin as the inner shell and maltose-binding protein as the outer shell. To assist homogeneous orientations of the target, disulfide bonds are introduced at the target-apoferritin interface, resulting in a cryo-EM structure at 2.6 Å resolution. We used molecular dynamics simulations to design peptides that block the interaction of the KIX domain of CBP with the intrinsically disordered pKID domain of CREB. The double-shell design allows for fluorescence polarization assays confirming the binding between the KIX domain in the double-shell and these interacting peptides. Further cryo-EM analysis reveals a helix-helix interaction between a single KIX helix and the best peptide, providing a possible strategy for developments of next-generation inhibitors.

Leveraging nonstructural data to predict structures and affinities of protein-ligand complexes.

Over the past five decades, tremendous effort has been devoted to computational methods for predicting properties of ligands-i.e., molecules that bind macromolecular targets. Such methods, which are critical to rational drug design, fall into two categories: physics-based methods, which directly model ligand interactions with the target given the target's three-dimensional (3D) structure, and ligand-based methods, which predict ligand properties given experimental measurements for similar ligands. Here, we present a rigorous statistical framework to combine these two sources of information. We develop a method to predict a ligand's pose-the 3D structure of the ligand bound to its target-that leverages a widely available source of information: a list of other ligands that are known to bind the same target but for which no 3D structure is available. This combination of physics-based and ligand-based modeling improves pose prediction accuracy across all major families of drug targets. Using the same framework, we develop a method for virtual screening of drug candidates, which outperforms standard physics-based and ligand-based virtual screening methods. Our results suggest broad opportunities to improve prediction of various ligand properties by combining diverse sources of information through customized machine-learning approaches.

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray-derived conformational ensembles.

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme's reaction cycle. To examine positioning and motions and test catalytic proposals, we collected "room temperature" X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI's reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI's general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure-function relationships, providing the basis for a new era of "ensemble-function" interrogation of enzymes.

Structure of hepcidin-bound ferroportin reveals iron homeostatic mechanisms.

The serum level of iron in humans is tightly controlled by the action of the hormone hepcidin on the iron efflux transporter ferroportin. Hepcidin regulates iron absorption and recycling by inducing the internalization and degradation of ferroportin1. Aberrant ferroportin activity can lead to diseases of iron overload, such as haemochromatosis, or iron limitation anaemias2. Here we determine cryogenic electron microscopy structures of ferroportin in lipid nanodiscs, both in the apo state and in complex with hepcidin and the iron mimetic cobalt. These structures and accompanying molecular dynamics simulations identify two metal-binding sites within the N and C domains of ferroportin. Hepcidin binds ferroportin in an outward-open conformation and completely occludes the iron efflux pathway to inhibit transport. The carboxy terminus of hepcidin directly contacts the divalent metal in the ferroportin C domain. Hepcidin binding to ferroportin is coupled to iron binding, with an 80-fold increase in hepcidin affinity in the presence of iron. These results suggest a model for hepcidin regulation of ferroportin, in which only ferroportin molecules loaded with iron are targeted for degradation. More broadly, our structural and functional insights may enable more targeted manipulation of the hepcidin-ferroportin axis in disorders of iron homeostasis.

Azide-Alkyne Click Conjugation on Quantum Dots by Selective Copper Coordination.

Functionalization of nanocrystals is essential for their practical application, but synthesis on nanocrystal surfaces is limited by incompatibilities with certain key reagents. The copper-catalyzed azide-alkyne cycloaddition is among the most useful methods for ligating molecules to surfaces, but has been largely useless for semiconductor quantum dots (QDs) because Cu+ ions quickly and irreversibly quench QD fluorescence. To discover nonquenching synthetic conditions for Cu-catalyzed click reactions on QD surfaces, we developed a combinatorial fluorescence assay to screen >2000 reaction conditions to maximize cycloaddition efficiency while minimizing QD quenching. We identify conditions for complete coupling without significant quenching, which are compatible with common QD polymer surfaces and various azide/alkyne pairs. Based on insight from the combinatorial screen and mechanistic studies of Cu coordination and quenching, we find that superstoichiometric concentrations of Cu can promote full coupling if accompanied by ligands that selectively compete with the Cu from the QD surface but allow it to remain catalytically active. Applied to the conjugation of a K+ channel-specific peptidyl toxin to CdSe/ZnS QDs, we synthesize unquenched QD conjugates and image their specific and voltage-dependent affinity for K+ channels in live cells.

Growth and assembly of cobalt oxide nanoparticle rings at liquid nanodroplets with solid junction.

Using liquid cell TEM, we imaged the formation of CoO nanoparticle rings. Nanoparticles nucleated and grew tracing the perimeter of droplets sitting on the SiNx solid substrate, and finally formed necklace-like rings. By tracking single nanoparticle trajectories during the ring formation and an estimation of the forces between droplets and nanoparticles using a simplified model, we found the junction of liquid nanodroplets with a solid substrate is the attractive site for CoO nanoparticles. Coalescing droplets were capable of pushing nanoparticles to the perimeter of the new droplet and nanoparticles on top of the droplets rolled off toward the perimeter. We propose that the curved surface morphology of the droplets created a force gradient that contributed to the assembly of nanoparticles at the droplet perimeter. Revealing the dynamics of nanoparticle movements and the interactions of nanoparticles with the liquid nanodroplet provides insights on developing novel self-assembly strategies for building precisely defined nanostructures on solid substrates.

Covalent Protein Labeling and Improved Single-Molecule Optical Properties of Aqueous CdSe/CdS Quantum Dots.

Semiconductor quantum dots (QDs) have proven to be superior probes for single-molecule imaging compared to organic or genetically encoded fluorophores, but they are limited by difficulties in protein targeting, their larger size, and on-off blinking. Here, we report compact aqueous CdSe/CdS QDs with significantly improved bioconjugation efficiency and superior single-molecule optical properties. We have synthesized covalent protein labeling ligands (i.e., SNAP tags) that are optimized for nanoparticle use, and QDs functionalized with these ligands label SNAP-tagged proteins ∼10-fold more efficiently than existing SNAP ligands. Single-molecule analysis of these QDs shows 99% of time spent in the fluorescent on-state, ∼4-fold higher quantum efficiency than standard CdSe/ZnS QDs, and 350 million photons detected before photobleaching. Bright signals of these QDs enable us to track the stepping movement of a kinesin motor in vitro, and the improved labeling efficiency enables tracking of single kinesins in live cells.

Characterizing Photon Reabsorption in Quantum Dot-Polymer Composites for Use as Displacement Sensors.

The reabsorption of photoluminescence within a medium, an effect known as the inner filter effect (IFE), has been well studied in solutions, but has garnered less attention in regards to solid-state nanocomposites. Photoluminescence from a quantum dot (QD) can selectively excite larger QDs around it resulting in a net red-shift in the reemitted photon. In CdSe/CdS core/shell QD-polymer nanocomposites, we observe a large spectral red-shift of over a third of the line width of the photoluminescence of the nanocomposites over a distance of 100 µm resulting from the IFE. Unlike fluorescent dyes, which do not show a large IFE red-shift, QDs have a component of inhomogeneous broadening that originates from their size distribution and quantum confinement. By controlling the photoluminescence broadening as well as the sample dispersion and concentration, we show that the magnitude of the IFE within the nanocomposite can be tuned. We further demonstrate that this shift can be exploited in order to spectroscopically monitor the vertical displacement of a nanocomposite in a fluorescence microscope. Large energetic shifts in the measured emission with displacement can be maximized, resulting in a displacement sensor with submicrometer resolution. We further show that the composite can be easily attached to biological samples and is able to measure deformations with high temporal and spatial precision.

Tracking Nanoparticle Diffusion and Interaction during Self-Assembly in a Liquid Cell.

Nanoparticle self-assembly has been well studied theoretically, but it remains challenging to directly observe and quantify individual nanoparticle interactions. With our custom image analysis method, we track the trajectories of nanoparticle movement with high precision from a stack of relatively noisy images obtained using liquid cell transmission electron microscopy. In a time frame of minutes, Pt-Fe nanoparticles self-assembled into a loosely packed hcp lattice. The energetics and stability of the dynamic assembly were studied quantitatively. From velocity and diffusion measurements, we experimentally determined the magnitude of forces between single particles and the related physical properties. The results illustrate that long-range anisotropic forces drive the formation of chains, which then clump and fold to maximize close range van der Waals interactions.

Single-particle mapping of nonequilibrium nanocrystal transformations.

Chemists have developed mechanistic insight into numerous chemical reactions by thoroughly characterizing nonequilibrium species. Although methods to probe these processes are well established for molecules, analogous techniques for understanding intermediate structures in nanomaterials have been lacking. We monitor the shape evolution of individual anisotropic gold nanostructures as they are oxidatively etched in a graphene liquid cell with a controlled redox environment. Short-lived, nonequilibrium nanocrystals are observed, structurally analyzed, and rationalized through Monte Carlo simulations. Understanding these reaction trajectories provides important fundamental insight connecting high-energy nanocrystal morphologies to the development of kinetically stabilized surface features and demonstrates the importance of developing tools capable of probing short-lived nanoscale species at the single-particle level.