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.

Structure-based discovery of conformationally selective inhibitors of the serotonin transporter.

The serotonin transporter (SERT) removes synaptic serotonin and is the target of anti-depressant drugs. SERT adopts three conformations: outward-open, occluded, and inward-open. All known inhibitors target the outward-open state except ibogaine, which has unusual anti-depressant and substance-withdrawal effects, and stabilizes the inward-open conformation. Unfortunately, ibogaine's promiscuity and cardiotoxicity limit the understanding of inward-open state ligands. We docked over 200 million small molecules against the inward-open state of the SERT. Thirty-six top-ranking compounds were synthesized, and thirteen inhibited; further structure-based optimization led to the selection of two potent (low nanomolar) inhibitors. These stabilized an outward-closed state of the SERT with little activity against common off-targets. A cryo-EM structure of one of these bound to the SERT confirmed the predicted geometry. In mouse behavioral assays, both compounds had anxiolytic- and anti-depressant-like activity, with potencies up to 200-fold better than fluoxetine (Prozac), and one substantially reversed morphine withdrawal effects.

An infectivity-enhancing site on the SARS-CoV-2 spike protein targeted by antibodies.

Antibodies against the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein prevent SARS-CoV-2 infection. However, the effects of antibodies against other spike protein domains are largely unknown. Here, we screened a series of anti-spike monoclonal antibodies from coronavirus disease 2019 (COVID-19) patients and found that some of antibodies against the N-terminal domain (NTD) induced the open conformation of RBD and thus enhanced the binding capacity of the spike protein to ACE2 and infectivity of SARS-CoV-2. Mutational analysis revealed that all of the infectivity-enhancing antibodies recognized a specific site on the NTD. Structural analysis demonstrated that all infectivity-enhancing antibodies bound to NTD in a similar manner. The antibodies against this infectivity-enhancing site were detected at high levels in severe patients. Moreover, we identified antibodies against the infectivity-enhancing site in uninfected donors, albeit at a lower frequency. These findings demonstrate that not only neutralizing antibodies but also enhancing antibodies are produced during SARS-CoV-2 infection.

ATG8-Binding UIM Proteins Define a New Class of Autophagy Adaptors and Receptors.

During autophagy, vesicle dynamics and cargo recruitment are driven by numerous adaptors and receptors that become tethered to the phagophore through interactions with lipidated ATG8/LC3 decorating the expanding membrane. Most currently described ATG8-binding proteins exploit a well-defined ATG8-interacting motif (AIM, or LC3-interacting region [LIR]) that contacts a hydrophobic patch on ATG8 known as the LIR/AIM docking site (LDS). Here we describe a new class of ATG8 interactors that exploit ubiquitin-interacting motif (UIM)-like sequences for high-affinity binding to an alternative ATG8 interaction site. Assays with candidate UIM-containing proteins together with unbiased screens identified a large collection of UIM-based ATG8 interactors in plants, yeast, and humans. Analysis of a subset also harboring ubiquitin regulatory X (UBX) domains revealed a role for UIM-directed autophagy in clearing non-functional CDC48/p97 complexes, including some impaired in human disease. With this new class of adaptors and receptors, we greatly extend the reach of selective autophagy and identify new factors regulating autophagic vesicle dynamics.

Irisin acts through its integrin receptor in a two-step process involving extracellular Hsp90α.

Exercise benefits the human body in many ways. Irisin is secreted by muscle, increased with exercise, and conveys physiological benefits, including improved cognition and resistance to neurodegeneration. Irisin acts via αV integrins; however, a mechanistic understanding of how small polypeptides like irisin can signal through integrins is poorly understood. Using mass spectrometry and cryo-EM, we demonstrate that the extracellular heat shock protein 90α (eHsp90α) is secreted by muscle with exercise and activates integrin αVß5. This allows for high-affinity irisin binding and signaling through an Hsp90α/αV/ß5 complex. By including hydrogen/deuterium exchange data, we generate and experimentally validate a 2.98 Å RMSD irisin/αVß5 complex docking model. Irisin binds very tightly to an alternative interface on αVß5 distinct from that used by known ligands. These data elucidate a non-canonical mechanism by which a small polypeptide hormone like irisin can function through an integrin receptor.

PEX5 translocation into and out of peroxisomes drives matrix protein import.

Peroxisomes are ubiquitous organelles whose dysfunction causes fatal human diseases. Most peroxisomal enzymes are imported from the cytosol by the receptor PEX5, which interacts with a docking complex in the peroxisomal membrane and then returns to the cytosol after monoubiquitination by a membrane-embedded ubiquitin ligase. The mechanism by which PEX5 shuttles between cytosol and peroxisomes and releases cargo inside the lumen is unclear. Here, we use Xenopus egg extract to demonstrate that PEX5 accompanies cargo completely into the lumen, utilizing WxxxF/Y motifs near its N terminus that bind a lumenal domain of the docking complex. PEX5 recycling is initiated by an amphipathic helix that binds to the lumenal side of the ubiquitin ligase. The N terminus then emerges in the cytosol for monoubiquitination. Finally, PEX5 is extracted from the lumen, resulting in the unfolding of the receptor and cargo release. Our results reveal the unique mechanism by which PEX5 ferries proteins into peroxisomes.

Recent advances in predicting and modeling protein-protein interactions.

Protein-protein interactions (PPIs) drive biological processes, and disruption of PPIs can cause disease. With recent breakthroughs in structure prediction and a deluge of genomic sequence data, computational methods to predict PPIs and model spatial structures of protein complexes are now approaching the accuracy of experimental approaches for permanent interactions and show promise for elucidating transient interactions. As we describe here, the key to this success is rich evolutionary information deciphered from thousands of homologous sequences that coevolve in interacting partners. This covariation signal, revealed by sophisticated statistical and machine learning (ML) algorithms, predicts physiological interactions. Accurate artificial intelligence (AI)-based modeling of protein structures promises to provide accurate 3D models of PPIs at a proteome-wide scale.

Localisation and function of key axonemal microtubule inner proteins and dynein docking complex members reveal extensive diversity among vertebrate motile cilia.

Vertebrate motile cilia are classified as (9+2) or (9+0), based on the presence or absence of the central pair apparatus, respectively. Cryogenic electron microscopy analyses of (9+2) cilia have uncovered an elaborate axonemal protein composition. The extent to which these features are conserved in (9+0) cilia remains unclear. CFAP53, a key axonemal filamentous microtubule inner protein (fMIP) and a centriolar satellites component, is essential for motility of (9+0), but not (9+2) cilia. Here, we show that in (9+2) cilia, CFAP53 functions redundantly with a paralogous fMIP, MNS1. MNS1 localises to ciliary axonemes, and combined loss of both proteins in zebrafish and mice caused severe outer dynein arm loss from (9+2) cilia, significantly affecting their motility. Using immunoprecipitation, we demonstrate that, whereas MNS1 can associate with itself and CFAP53, CFAP53 is unable to self-associate. We also show that additional axonemal dynein-interacting proteins, two outer dynein arm docking (ODAD) complex members, show differential localisation between types of motile cilia. Together, our findings clarify how paralogous fMIPs, CFAP53 and MNS1, function in regulating (9+2) versus (9+0) cilia motility, and further emphasise extensive structural diversity among these organelles.

Cyclin-Specific Docking Mechanisms Reveal the Complexity of M-CDK Function in the Cell Cycle.

Cyclin-dependent kinases (CDKs) coordinate hundreds of molecular events during the cell cycle. Multiple cyclins are involved, but the global role of cyclin-specific phosphorylation has remained unsolved. We uncovered a cyclin docking motif, LxF, that mediates binding of replication factor Cdc6 to mitotic cyclin. This interaction leads to phospho-adaptor Cks1-mediated inhibition of M-CDK to facilitate Cdc6 accumulation and sequestration in mitosis. The LxF motif and Cks1 also mediate the mutual inhibition between M-CDK and the tyrosine kinase Swe1. Additionally, the LxF motif is critical for targeting M-CDK to phosphorylate several mitotic regulators; for example, Spo12 is targeted via LxF to release the phosphatase Cdc14. The results complete the full set of G1, S, and M-CDK docking mechanisms and outline the unified role of cyclin specificity and CDK activity thresholds. Cooperation of cyclin and Cks1 docking creates a variety of CDK thresholds and switching orders, including combinations of last in, first out (LIFO) and first in, first out (FIFO) ordering.

Cyclin D-Cdk4,6 Drives Cell-Cycle Progression via the Retinoblastoma Protein's C-Terminal Helix.

The cyclin-dependent kinases Cdk4 and Cdk6 form complexes with D-type cyclins to drive cell proliferation. A well-known target of cyclin D-Cdk4,6 is the retinoblastoma protein Rb, which inhibits cell-cycle progression until its inactivation by phosphorylation. However, the role of Rb phosphorylation by cyclin D-Cdk4,6 in cell-cycle progression is unclear because Rb can be phosphorylated by other cyclin-Cdks, and cyclin D-Cdk4,6 has other targets involved in cell division. Here, we show that cyclin D-Cdk4,6 docks one side of an alpha-helix in the Rb C terminus, which is not recognized by cyclins E, A, and B. This helix-based docking mechanism is shared by the p107 and p130 Rb-family members across metazoans. Mutation of the Rb C-terminal helix prevents its phosphorylation, promotes G1 arrest, and enhances Rb's tumor suppressive function. Our work conclusively demonstrates that the cyclin D-Rb interaction drives cell division and expands the diversity of known cyclin-based protein docking mechanisms.

ABAG-docking benchmark: a non-redundant structure benchmark dataset for antibody-antigen computational docking.

Accurate prediction of antibody-antigen complex structures is pivotal in drug discovery, vaccine design and disease treatment and can facilitate the development of more effective therapies and diagnostics. In this work, we first review the antibody-antigen docking (ABAG-docking) datasets. Then, we present the creation and characterization of a comprehensive benchmark dataset of antibody-antigen complexes. We categorize the dataset based on docking difficulty, interface properties and structural characteristics, to provide a diverse set of cases for rigorous evaluation. Compared with Docking Benchmark 5.5, we have added 112 cases, including 14 single-domain antibody (sdAb) cases and 98 monoclonal antibody (mAb) cases, and also increased the proportion of Difficult cases. Our dataset contains diverse cases, including human/humanized antibodies, sdAbs, rodent antibodies and other types, opening the door to better algorithm development. Furthermore, we provide details on the process of building the benchmark dataset and introduce a pipeline for periodic updates to keep it up to date. We also utilize multiple complex prediction methods including ZDOCK, ClusPro, HDOCK and AlphaFold-Multimer for testing and analyzing this dataset. This benchmark serves as a valuable resource for evaluating and advancing docking computational methods in the analysis of antibody-antigen interaction, enabling researchers to develop more accurate and effective tools for predicting and designing antibody-antigen complexes. The non-redundant ABAG-docking structure benchmark dataset is available at https://github.com/Zhaonan99/Antibody-antigen-complex-structure-benchmark-dataset.

GSScore: a novel Graphormer-based shell-like scoring method for protein-ligand docking.

Protein-ligand interactions (PLIs) are essential for cellular activities and drug discovery. But due to the complexity and high cost of experimental methods, there is a great demand for computational approaches to recognize PLI patterns, such as protein-ligand docking. In recent years, more and more models based on machine learning have been developed to directly predict the root mean square deviation (RMSD) of a ligand docking pose with reference to its native binding pose. However, new scoring methods are pressingly needed in methodology for more accurate RMSD prediction. We present a new deep learning-based scoring method for RMSD prediction of protein-ligand docking poses based on a Graphormer method and Shell-like graph architecture, named GSScore. To recognize near-native conformations from a set of poses, GSScore takes atoms as nodes and then establishes the docking interface of protein-ligand into multiple bipartite graphs within different shell ranges. Benefiting from the Graphormer and Shell-like graph architecture, GSScore can effectively capture the subtle differences between energetically favorable near-native conformations and unfavorable non-native poses without extra information. GSScore was extensively evaluated on diverse test sets including a subset of PDBBind version 2019, CASF2016 as well as DUD-E, and obtained significant improvements over existing methods in terms of RMSE, $R$ (Pearson correlation coefficient), Spearman correlation coefficient and Docking power.

RmsdXNA: RMSD prediction of nucleic acid-ligand docking poses using machine-learning method.

Small molecule drugs can be used to target nucleic acids (NA) to regulate biological processes. Computational modeling methods, such as molecular docking or scoring functions, are commonly employed to facilitate drug design. However, the accuracy of the scoring function in predicting the closest-to-native docking pose is often suboptimal. To overcome this problem, a machine learning model, RmsdXNA, was developed to predict the root-mean-square-deviation (RMSD) of ligand docking poses in NA complexes. The versatility of RmsdXNA has been demonstrated by its successful application to various complexes involving different types of NA receptors and ligands, including metal complexes and short peptides. The predicted RMSD by RmsdXNA was strongly correlated with the actual RMSD of the docked poses. RmsdXNA also outperformed the rDock scoring function in ranking and identifying closest-to-native docking poses across different structural groups and on the testing dataset. Using experimental validated results conducted on polyadenylated nuclear element for nuclear expression triplex, RmsdXNA demonstrated better screening power for the RNA-small molecule complex compared to rDock. Molecular dynamics simulations were subsequently employed to validate the binding of top-scoring ligand candidates selected by RmsdXNA and rDock on MALAT1. The results showed that RmsdXNA has a higher success rate in identifying promising ligands that can bind well to the receptor. The development of an accurate docking score for a NA-ligand complex can aid in drug discovery and development advancements. The code to use RmsdXNA is available at the GitHub repository https://github.com/laiheng001/RmsdXNA.

Enhancing cryo-EM structure prediction with DeepTracer and AlphaFold2 integration.

Understanding the protein structures is invaluable in various biomedical applications, such as vaccine development. Protein structure model building from experimental electron density maps is a time-consuming and labor-intensive task. To address the challenge, machine learning approaches have been proposed to automate this process. Currently, the majority of the experimental maps in the database lack atomic resolution features, making it challenging for machine learning-based methods to precisely determine protein structures from cryogenic electron microscopy density maps. On the other hand, protein structure prediction methods, such as AlphaFold2, leverage evolutionary information from protein sequences and have recently achieved groundbreaking accuracy. However, these methods often require manual refinement, which is labor intensive and time consuming. In this study, we present DeepTracer-Refine, an automated method that refines AlphaFold predicted structures by aligning them to DeepTracers modeled structure. Our method was evaluated on 39 multi-domain proteins and we improved the average residue coverage from 78.2 to 90.0% and average local Distance Difference Test score from 0.67 to 0.71. We also compared DeepTracer-Refine with Phenixs AlphaFold refinement and demonstrated that our method not only performs better when the initial AlphaFold model is less precise but also surpasses Phenix in run-time performance.

Linear motif specificity in signaling through p38α and ERK2 mitogen-activated protein kinases.

Mitogen-activated protein kinase (MAPK) cascades are essential for eukaryotic cells to integrate and respond to diverse stimuli. Maintaining specificity in signaling through MAPK networks is key to coupling distinct inputs to appropriate cellular responses. Docking sites-short linear motifs found in MAPK substrates, regulators, and scaffolds-can promote signaling specificity through selective interactions, but how they do so remains unresolved. Here, we screened a proteomic library for sequences interacting with the MAPKs extracellular signal-regulated kinase 2 (ERK2) and p38α, identifying selective and promiscuous docking motifs. Sequences specific for p38α had high net charge and lysine content, and selective binding depended on a pair of acidic residues unique to the p38α docking interface. Finally, we validated a set of full-length proteins harboring docking sites selected in our screens to be authentic MAPK interactors and substrates. This study identifies features that help define MAPK signaling networks and explains how specific docking motifs promote signaling integrity.

Protein Amyloid Cofactors: Charged Side-Chain Arrays Meet Their Match?

Recent advances in high-resolution structural studies of protein amyloids have revealed parallel in-register cross-ß-sheets with periodic arrays of closely spaced identical residues. What do these structures tell us about the mechanisms of action of common amyloid-promoting factors, such as heparan sulfate (HS), nucleic acids, polyphosphates, anionic phospholipids, and acidic pH?

Selective inhibition of overactive warmth-sensitive Ca2+-permeable TRPV3 channels by antispasmodic agent flopropione for alleviation of skin inflammation.

The temperature-sensitive Ca2+-permeable TRPV3 ion channel is robustly expressed in the skin keratinocytes, and its gain-of-function mutations are involved in the pathology of skin lesions. Here, we report the identification of an antispasmodic agent flopropione that alleviates skin inflammation by selective inhibition of TRPV3. In whole-cell patch clamp recordings, flopropione selectively inhibits macroscopic TRPV3 currents in a concentration-dependent manner with an IC50 value of 17.8 ± 3.5 µM. At the single-channel level, flopropione inhibits TRPV3 channel open probability without alteration of its unitary conductance. In an in vivo mouse model of skin inflammation induced by the skin sensitizer DNFB, flopropione also alleviates dorsal skin lesions and ear skin swelling. Further molecular docking combined with site-directed mutagenesis reveals that two residues E501 and I505 in the channel S2-helix are critical for flopropione-mediated inhibition of TRPV3. Taken together, our findings demonstrate that the spasmolytic drug flopropione as a selective inhibitor of TRPV3 channel not only provides a valuable tool molecule for understanding of TRPV3 channel pharmacology but also holds repurposing potential for therapy of skin disorders, such as dermatitis and pruritus.

High-throughput virtual search of small molecules for controlling the mechanical stability of human CD4.

Protein mechanical stability determines the function of a myriad of proteins, especially proteins from the extracellular matrix. Failure to maintain protein mechanical stability may result in diseases and disorders such as cancer, cardiomyopathies, or muscular dystrophy. Thus, developing mutation-free approaches to enhance and control the mechanical stability of proteins using pharmacology-based methods may have important implications in drug development and discovery. Here, we present the first approach that employs computational high-throughput virtual screening and molecular docking to search for small molecules in chemical libraries that function as mechano-regulators of the stability of human cluster of differentiation 4, receptor of HIV-1. Using single-molecule force spectroscopy, we prove that these small molecules can increase the mechanical stability of CD4D1D2 domains over 4-fold in addition to modifying the mechanical unfolding pathways. Our experiments demonstrate that chemical libraries are a source of mechanoactive molecules and that drug discovery approaches provide the foundation of a new type of molecular function, that is, mechano-regulation, paving the way toward mechanopharmacology.

The small molecule activator S3969 stimulates the epithelial sodium channel by interacting with a specific binding pocket in the channel's ß-subunit.

The epithelial sodium channel (ENaC) is essential for mediating sodium absorption in several epithelia. Its impaired function leads to severe disorders, including pseudohypoaldosteronism type 1 and respiratory distress. Therefore, pharmacological ENaC activators have potential therapeutic implications. Previously, a small molecule ENaC activator (S3969) was developed. So far, little is known about molecular mechanisms involved in S3969-mediated ENaC stimulation. Here, we identified an S3969-binding site in human ENaC by combining structure-based simulations with molecular biological methods and electrophysiological measurements of ENaC heterologously expressed in Xenopus laevis oocytes. We confirmed a previous observation that the extracellular loop of ß-ENaC is essential for ENaC stimulation by S3969. Molecular dynamics simulations predicted critical residues in the thumb domain of ß-ENaC (Arg388, Phe391, and Tyr406) that coordinate S3969 within a binding site localized at the ß-γ-subunit interface. Importantly, mutating each of these residues reduced (R388H; R388A) or nearly abolished (F391G; Y406A) the S3969-mediated ENaC activation. Molecular dynamics simulations also suggested that S3969-mediated ENaC stimulation involved a movement of the α5 helix of the thumb domain of ß-ENaC away from the palm domain of γ-ENaC. Consistent with this, the introduction of two cysteine residues (ßR437C - γS298C) to form a disulfide bridge connecting these two domains prevented ENaC stimulation by S3969 unless the disulfide bond was reduced by DTT. Finally, we demonstrated that S3969 stimulated ENaC endogenously expressed in cultured human airway epithelial cells (H441). These new findings may lead to novel (patho-)physiological and therapeutic concepts for disorders associated with altered ENaC function.

A new linear cyclin docking motif that mediates exclusively S-phase CDK-specific signaling.

Cyclin-dependent kinases (CDKs), the master regulators of cell division, are activated by different cyclins at different cell cycle stages. In addition to being activators of CDKs, cyclins recognize various linear motifs to target CDK activity to specific proteins. We uncovered a cyclin docking motif, NLxxxL, that contributes to phosphorylation-dependent degradation of the CDK inhibitor Far1 at the G1/S stage in the yeast Saccharomyces cerevisiae. This motif is recognized exclusively by S-phase CDK (S-CDK) Clb5/6-Cdc28 and is considerably more potent than the conventional RxL docking motif. The NLxxxL and RxL motifs were found to overlap in some target proteins, suggesting that cyclin docking motifs can evolve to switch from one to another for fine-tuning of cell cycle events. Using time-lapse fluorescence microscopy, we show how different docking connections temporally control phosphorylation-driven target degradation. This also revealed a differential function of the phosphoadaptor protein Cks1, as Cks1 docking potentiated degron phosphorylation of RxL-containing but not of NLxxxL-containing substrates. The NLxxxL motif was found to govern S-cyclin-specificity in multiple yeast CDK targets including Fin1, Lif1, and Slx4, suggesting its wider importance.