The Art and Science of Molecular Docking.

Molecular docking has become an essential part of a structural biologist's and medicinal chemist's toolkits. Given a chemical compound and the three-dimensional structure of a molecular target-for example, a protein-docking methods fit the compound into the target, predicting the compound's bound structure and binding energy. Docking can be used to discover novel ligands for a target by screening large virtual compound libraries. Docking can also provide a useful starting point for structure-based ligand optimization or for investigating a ligand's mechanism of action. Advances in computational methods, including both physics-based and machine learning approaches, as well as in complementary experimental techniques, are making docking an even more powerful tool. We review how docking works and how it can drive drug discovery and biological research. We also describe its current limitations and ongoing efforts to overcome them.

Structure-based discovery of CFTR potentiators and inhibitors.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, whereas its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here, we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify CFTR modulators. We docked ∼155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered mid-nanomolar potentiators, as well as inhibitors, that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

Comprehensive structure and functional adaptations of the yeast nuclear pore complex.

Nuclear pore complexes (NPCs) mediate the nucleocytoplasmic transport of macromolecules. Here we provide a structure of the isolated yeast NPC in which the inner ring is resolved by cryo-EM at sub-nanometer resolution to show how flexible connectors tie together different structural and functional layers. These connectors may be targets for phosphorylation and regulated disassembly in cells with an open mitosis. Moreover, some nucleoporin pairs and transport factors have similar interaction motifs, which suggests an evolutionary and mechanistic link between assembly and transport. We provide evidence for three major NPC variants that may foreshadow functional specializations at the nuclear periphery. Cryo-electron tomography extended these studies, providing a model of the in situ NPC with a radially expanded inner ring. Our comprehensive model reveals features of the nuclear basket and central transporter, suggests a role for the lumenal Pom152 ring in restricting dilation, and highlights structural plasticity that may be required for transport.

Ribonucleotide Reductases: Structure, Chemistry, and Metabolism Suggest New Therapeutic Targets.

Ribonucleotide reductases (RNRs) catalyze the de novo conversion of nucleotides to deoxynucleotides in all organisms, controlling their relative ratios and abundance. In doing so, they play an important role in fidelity of DNA replication and repair. RNRs' central role in nucleic acid metabolism has resulted in five therapeutics that inhibit human RNRs. In this review, we discuss the structural, dynamic, and mechanistic aspects of RNR activity and regulation, primarily for the human and Escherichia coli class Ia enzymes. The unusual radical-based organic chemistry of nucleotide reduction, the inorganic chemistry of the essential metallo-cofactor biosynthesis/maintenance, the transport of a radical over a long distance, and the dynamics of subunit interactions all present distinct entry points toward RNR inhibition that are relevant for drug discovery. We describe the current mechanistic understanding of small molecules that target different elements of RNR function, including downstream pathways that lead to cell cytotoxicity. We conclude by summarizing novel and emergent RNR targeting motifs for cancer and antibiotic therapeutics.

Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding.

The receptor binding domain (RBD) of the SARS-CoV-2 spike glycoprotein mediates viral attachment to ACE2 receptor and is a major determinant of host range and a dominant target of neutralizing antibodies. Here, we experimentally measure how all amino acid mutations to the RBD affect expression of folded protein and its affinity for ACE2. Most mutations are deleterious for RBD expression and ACE2 binding, and we identify constrained regions on the RBD's surface that may be desirable targets for vaccines and antibody-based therapeutics. But a substantial number of mutations are well tolerated or even enhance ACE2 binding, including at ACE2 interface residues that vary across SARS-related coronaviruses. However, we find no evidence that these ACE2-affinity-enhancing mutations have been selected in current SARS-CoV-2 pandemic isolates. We present an interactive visualization and open analysis pipeline to facilitate use of our dataset for vaccine design and functional annotation of mutations observed during viral surveillance.

Cryo-EM Structure of the Human Cannabinoid Receptor CB2-Gi Signaling Complex.

Drugs selectively targeting CB2 hold promise for treating neurodegenerative disorders, inflammation, and pain while avoiding psychotropic side effects mediated by CB1. The mechanisms underlying CB2 activation and signaling are poorly understood but critical for drug design. Here we report the cryo-EM structure of the human CB2-Gi signaling complex bound to the agonist WIN 55,212-2. The 3D structure reveals the binding mode of WIN 55,212-2 and structural determinants for distinguishing CB2 agonists from antagonists, which are supported by a pair of rationally designed agonist and antagonist. Further structural analyses with computational docking results uncover the differences between CB2 and CB1 in receptor activation, ligand recognition, and Gi coupling. These findings are expected to facilitate rational structure-based discovery of drugs targeting the cannabinoid system.

Full-Length P2X7 Structures Reveal How Palmitoylation Prevents Channel Desensitization.

P2X receptors are trimeric, non-selective cation channels activated by extracellular ATP. The P2X7 receptor subtype is a pharmacological target because of involvement in apoptotic, inflammatory, and tumor progression pathways. It is the most structurally and functionally distinct P2X subtype, containing a unique cytoplasmic domain critical for the receptor to initiate apoptosis and not undergo desensitization. However, lack of structural information about the cytoplasmic domain has hindered understanding of the molecular mechanisms underlying these processes. We report cryoelectron microscopy structures of full-length rat P2X7 receptor in apo and ATP-bound states. These structures reveal how one cytoplasmic element, the C-cys anchor, prevents desensitization by anchoring the pore-lining helix to the membrane with palmitoyl groups. They show a second cytoplasmic element with a unique fold, the cytoplasmic ballast, which unexpectedly contains a zinc ion complex and a guanosine nucleotide binding site. Our structures provide first insights into the architecture and function of a P2X receptor cytoplasmic domain.

Crystal Structure of the Human Cannabinoid Receptor CB2.

The cannabinoid receptor CB2 is predominately expressed in the immune system, and selective modulation of CB2 without the psychoactivity of CB1 has therapeutic potential in inflammatory, fibrotic, and neurodegenerative diseases. Here, we report the crystal structure of human CB2 in complex with a rationally designed antagonist, AM10257, at 2.8 Å resolution. The CB2-AM10257 structure reveals a distinctly different binding pose compared with CB1. However, the extracellular portion of the antagonist-bound CB2 shares a high degree of conformational similarity with the agonist-bound CB1, which led to the discovery of AM10257's unexpected opposing functional profile of CB2 antagonism versus CB1 agonism. Further structural analysis using mutagenesis studies and molecular docking revealed the molecular basis of their function and selectivity for CB2 and CB1. Additional analyses of our designed antagonist and agonist pairs provide important insight into the activation mechanism of CB2. The present findings should facilitate rational drug design toward precise modulation of the endocannabinoid system.

Allosteric Coupling of Drug Binding and Intracellular Signaling in the A2A Adenosine Receptor.

Signaling across cellular membranes, the 826 human G protein-coupled receptors (GPCRs) govern a wide range of vital physiological processes, making GPCRs prominent drug targets. X-ray crystallography provided GPCR molecular architectures, which also revealed the need for additional structural dynamics data to support drug development. Here, nuclear magnetic resonance (NMR) spectroscopy with the wild-type-like A2A adenosine receptor (A2AAR) in solution provides a comprehensive characterization of signaling-related structural dynamics. All six tryptophan indole and eight glycine backbone 15N-1H NMR signals in A2AAR were individually assigned. These NMR probes provided insight into the role of Asp522.50 as an allosteric link between the orthosteric drug binding site and the intracellular signaling surface, revealing strong interactions with the toggle switch Trp 2466.48, and delineated the structural response to variable efficacy of bound drugs across A2AAR. The present data support GPCR signaling based on dynamic interactions between two semi-independent subdomains connected by an allosteric switch at Asp522.50.

Structure of the Nanobody-Stabilized Active State of the Kappa Opioid Receptor.

The κ-opioid receptor (KOP) mediates the actions of opioids with hallucinogenic, dysphoric, and analgesic activities. The design of KOP analgesics devoid of hallucinatory and dysphoric effects has been hindered by an incomplete structural and mechanistic understanding of KOP agonist actions. Here, we provide a crystal structure of human KOP in complex with the potent epoxymorphinan opioid agonist MP1104 and an active-state-stabilizing nanobody. Comparisons between inactive- and active-state opioid receptor structures reveal substantial conformational changes in the binding pocket and intracellular and extracellular regions. Extensive structural analysis and experimental validation illuminate key residues that propagate larger-scale structural rearrangements and transducer binding that, collectively, elucidate the structural determinants of KOP pharmacology, function, and biased signaling. These molecular insights promise to accelerate the structure-guided design of safer and more effective κ-opioid receptor therapeutics.

Structures of DPAGT1 Explain Glycosylation Disease Mechanisms and Advance TB Antibiotic Design.

Protein N-glycosylation is a widespread post-translational modification. The first committed step in this process is catalysed by dolichyl-phosphate N-acetylglucosamine-phosphotransferase DPAGT1 (GPT/E.C. 2.7.8.15). Missense DPAGT1 variants cause congenital myasthenic syndrome and disorders of glycosylation. In addition, naturally-occurring bactericidal nucleoside analogues such as tunicamycin are toxic to eukaryotes due to DPAGT1 inhibition, preventing their clinical use. Our structures of DPAGT1 with the substrate UDP-GlcNAc and tunicamycin reveal substrate binding modes, suggest a mechanism of catalysis, provide an understanding of how mutations modulate activity (thus causing disease) and allow design of non-toxic "lipid-altered" tunicamycins. The structure-tuned activity of these analogues against several bacterial targets allowed the design of potent antibiotics for Mycobacterium tuberculosis, enabling treatment in vitro, in cellulo and in vivo, providing a promising new class of antimicrobial drug.

Structure of the Human cGAS-DNA Complex Reveals Enhanced Control of Immune Surveillance.

Cyclic GMP-AMP synthase (cGAS) recognition of cytosolic DNA is critical for immune responses to pathogen replication, cellular stress, and cancer. Existing structures of the mouse cGAS-DNA complex provide a model for enzyme activation but do not explain why human cGAS exhibits severely reduced levels of cyclic GMP-AMP (cGAMP) synthesis compared to other mammals. Here, we discover that enhanced DNA-length specificity restrains human cGAS activation. Using reconstitution of cGAMP signaling in bacteria, we mapped the determinant of human cGAS regulation to two amino acid substitutions in the DNA-binding surface. Human-specific substitutions are necessary and sufficient to direct preferential detection of long DNA. Crystal structures reveal why removal of human substitutions relaxes DNA-length specificity and explain how human-specific DNA interactions favor cGAS oligomerization. These results define how DNA-sensing in humans adapted for enhanced specificity and provide a model of the active human cGAS-DNA complex to enable structure-guided design of cGAS therapeutics.

Targeting KRAS Mutant Cancers with a Covalent G12C-Specific Inhibitor.

KRASG12C was recently identified to be potentially druggable by allele-specific covalent targeting of Cys-12 in vicinity to an inducible allosteric switch II pocket (S-IIP). Success of this approach requires active cycling of KRASG12C between its active-GTP and inactive-GDP conformations as accessibility of the S-IIP is restricted only to the GDP-bound state. This strategy proved feasible for inhibiting mutant KRAS in vitro; however, it is uncertain whether this approach would translate to in vivo. Here, we describe structure-based design and identification of ARS-1620, a covalent compound with high potency and selectivity for KRASG12C. ARS-1620 achieves rapid and sustained in vivo target occupancy to induce tumor regression. We use ARS-1620 to dissect oncogenic KRAS dependency and demonstrate that monolayer culture formats significantly underestimate KRAS dependency in vivo. This study provides in vivo evidence that mutant KRAS can be selectively targeted and reveals ARS-1620 as representing a new generation of KRASG12C-specific inhibitors with promising therapeutic potential.

Structural Insights into Yeast Telomerase Recruitment to Telomeres.

Telomerase maintains chromosome ends from humans to yeasts. Recruitment of yeast telomerase to telomeres occurs through its Ku and Est1 subunits via independent interactions with telomerase RNA (TLC1) and telomeric proteins Sir4 and Cdc13, respectively. However, the structures of the molecules comprising these telomerase-recruiting pathways remain unknown. Here, we report crystal structures of the Ku heterodimer and Est1 complexed with their key binding partners. Two major findings are as follows: (1) Ku specifically binds to telomerase RNA in a distinct, yet related, manner to how it binds DNA; and (2) Est1 employs two separate pockets to bind distinct motifs of Cdc13. The N-terminal Cdc13-binding site of Est1 cooperates with the TLC1-Ku-Sir4 pathway for telomerase recruitment, whereas the C-terminal interface is dispensable for binding Est1 in vitro yet is nevertheless essential for telomere maintenance in vivo. Overall, our results integrate previous models and provide fundamentally valuable structural information regarding telomere biology.

The p53 Pathway: Origins, Inactivation in Cancer, and Emerging Therapeutic Approaches.

Inactivation of the transcription factor p53, through either direct mutation or aberrations in one of its many regulatory pathways, is a hallmark of virtually every tumor. In recent years, screening for p53 activators and a better understanding of the molecular mechanisms of oncogenic perturbations of p53 function have opened up a host of novel avenues for therapeutic intervention in cancer: from the structure-guided design of chemical chaperones to restore the function of conformationally unstable p53 cancer mutants, to the development of potent antagonists of the negative regulators MDM2 and MDMX and other modulators of the p53 pathway for the treatment of cancers with wild-type p53. Some of these compounds have now moved from proof-of-concept studies into clinical trials, with prospects for further, personalized anticancer medicines. We trace the structural evolution of the p53 pathway, from germ-line surveillance in simple multicellular organisms to its pluripotential role in humans.

Accurate structure prediction of biomolecular interactions with AlphaFold 3.

The introduction of AlphaFold 21 has spurred a revolution in modelling the structure of proteins and their interactions, enabling a huge range of applications in protein modelling and design2-6. Here we describe our AlphaFold 3 model with a substantially updated diffusion-based architecture that is capable of predicting the joint structure of complexes including proteins, nucleic acids, small molecules, ions and modified residues. The new AlphaFold model demonstrates substantially improved accuracy over many previous specialized tools: far greater accuracy for protein-ligand interactions compared with state-of-the-art docking tools, much higher accuracy for protein-nucleic acid interactions compared with nucleic-acid-specific predictors and substantially higher antibody-antigen prediction accuracy compared with AlphaFold-Multimer v.2.37,8. Together, these results show that high-accuracy modelling across biomolecular space is possible within a single unified deep-learning framework.

The molecular basis of sugar detection by an insect taste receptor.

Animals crave sugars because of their energy potential and the pleasurable sensation of tasting sweetness. Yet all sugars are not metabolically equivalent, requiring mechanisms to detect and differentiate between chemically similar sweet substances. Insects use a family of ionotropic gustatory receptors to discriminate sugars1, each of which is selectively activated by specific sweet molecules2-6. Here, to gain insight into the molecular basis of sugar selectivity, we determined structures of Gr9, a gustatory receptor from the silkworm Bombyx mori (BmGr9), in the absence and presence of its sole activating ligand, D-fructose. These structures, along with structure-guided mutagenesis and functional assays, illustrate how D-fructose is enveloped by a ligand-binding pocket that precisely matches the overall shape and pattern of chemical groups in D-fructose. However, our computational docking and experimental binding assays revealed that other sugars also bind BmGr9, yet they are unable to activate the receptor. We determined the structure of BmGr9 in complex with one such non-activating sugar, L-sorbose. Although both sugars bind a similar position, only D-fructose is capable of engaging a bridge of two conserved aromatic residues that connects the pocket to the pore helix, inducing a conformational change that allows the ion-conducting pore to open. Thus, chemical specificity does not depend solely on the selectivity of the ligand-binding pocket, but it is an emergent property arising from a combination of receptor-ligand interactions and allosteric coupling. Our results support a model whereby coarse receptor tuning is derived from the size and chemical characteristics of the pocket, whereas fine-tuning of receptor activation is achieved through the selective engagement of an allosteric pathway that regulates ion conduction.

NudC guides client transfer between the Hsp40/70 and Hsp90 chaperone systems.

In the eukaryotic cytosol, the Hsp70 and the Hsp90 chaperone machines work in tandem with the maturation of a diverse array of client proteins. The transfer of nonnative clients between these systems is essential to the chaperoning process, but how it is regulated is still not clear. We discovered that NudC is an essential transfer factor with an unprecedented mode of action: NudC interacts with Hsp40 in Hsp40-Hsp70-client complexes and displaces Hsp70. Then, the interaction of NudC with Hsp90 allows the direct transfer of Hsp40-bound clients to Hsp90 for further processing. Consistent with this mechanism, NudC increases client activation in vitro as well as in cells and is essential for cellular viability. Together, our results show the complexity of the cooperation between the major chaperone machineries in the eukaryotic cytosol.

Mechanism of PINK1 activation by autophosphorylation and insights into assembly on the TOM complex.

Mutations in PINK1 cause autosomal-recessive Parkinson's disease. Mitochondrial damage results in PINK1 import arrest on the translocase of the outer mitochondrial membrane (TOM) complex, resulting in the activation of its ubiquitin kinase activity by autophosphorylation and initiation of Parkin-dependent mitochondrial clearance. Herein, we report crystal structures of the entire cytosolic domain of insect PINK1. Our structures reveal a dimeric autophosphorylation complex targeting phosphorylation at the invariant Ser205 (human Ser228). The dimer interface requires insert 2, which is unique to PINK1. The structures also reveal how an N-terminal helix binds to the C-terminal extension and provide insights into stabilization of PINK1 on the core TOM complex.

Pre-T cell receptor self-MHC sampling restricts thymocyte dedifferentiation.

Programming T cells to distinguish self from non-self is a vital, multi-step process that occurs in the thymus1-4. Signalling through the pre-T cell receptor (preTCR), a CD3-associated heterodimer comprising an invariant pTα chain and a clone-specific ß chain, is a critical early checkpoint in thymocyte development within the αß T cell lineage5,6. PreTCRs arrayed on CD4-CD8- double-negative thymocytes ligate peptides bound to major histocompatibility complex molecules (pMHC) on thymic stroma, similar to αß T cell receptors that appear on CD4+CD8+ double-positive thymocytes, but via a different molecular docking strategy7-10. Here we show the consequences of these distinct interactions for thymocyte progression using synchronized fetal thymic progenitor cultures that differ in the presence or absence of pMHC on support stroma, and single-cell transcriptomes at key thymocyte developmental transitions. Although major histocompatibility complex (MHC)-negative stroma fosters αß T cell differentiation, the absence of preTCR-pMHC interactions leads to deviant thymocyte transcriptional programming associated with dedifferentiation. Highly proliferative double-negative and double-positive thymocyte subsets emerge, with antecedent characteristics of T cell lymphoblastic and myeloid malignancies. Compensatory upregulation of diverse MHC class Ib proteins in B2m/H2-Ab1 MHC-knockout mice partially safeguards in vivo thymocyte progression, although disseminated double-positive thymic tumours may develop with ageing. Thus, as well as promoting ß chain repertoire broadening for subsequent αß T cell receptor utilization, preTCR-pMHC interactions limit cellular plasticity to facilitate normal thymocyte differentiation and proliferation that, if absent, introduce developmental vulnerabilities.