RESUMO
Kainate receptors, a subclass of ionotropic glutamate receptors, are tetrameric ligand-gated ion channels that mediate excitatory neurotransmission1-4. Kainate receptors modulate neuronal circuits and synaptic plasticity during the development and function of the central nervous system and are implicated in various neurological and psychiatric diseases, including epilepsy, depression, schizophrenia, anxiety and autism5-11. Although structures of kainate receptor domains and subunit assemblies are available12-18, the mechanism of kainate receptor gating remains poorly understood. Here we present cryo-electron microscopy structures of the kainate receptor GluK2 in the presence of the agonist glutamate and the positive allosteric modulators lectin concanavalin A and BPAM344. Concanavalin A and BPAM344 inhibit kainate receptor desensitization and prolong activation by acting as a spacer between the amino-terminal and ligand-binding domains and a stabilizer of the ligand-binding domain dimer interface, respectively. Channel opening involves the kinking of all four pore-forming M3 helices. Our structures reveal the molecular basis of kainate receptor gating, which could guide the development of drugs for treatment of neurological disorders.
Assuntos
Concanavalina A , Microscopia Crioeletrônica , Receptor de GluK2 Cainato , Ácido Glutâmico , Ativação do Canal Iônico , Modelos Moleculares , Domínios Proteicos , Receptores de Ácido Caínico , Receptores de Ácido Caínico/química , Receptores de Ácido Caínico/metabolismo , Receptores de Ácido Caínico/ultraestrutura , Humanos , Ácido Glutâmico/metabolismo , Ácido Glutâmico/química , Animais , Concanavalina A/química , Concanavalina A/metabolismo , Concanavalina A/farmacologia , Ligantes , Regulação Alostérica , Sítios de LigaçãoRESUMO
Ionotropic glutamate receptors (iGluRs) are tetrameric ligand-gated ion channels that open their pores in response to binding of the agonist glutamate1-3. An ionic current through a single iGluR channel shows up to four discrete conductance levels (O1-O4)4-6. Higher conductance levels have been associated with an increased number of agonist molecules bound to four individual ligand-binding domains (LBDs)6-10. Here we determine structures of a synaptic complex of AMPA-subtype iGluR and the auxiliary subunit γ2 in non-desensitizing conditions with various occupancy of the LBDs by glutamate. We show that glutamate binds to LBDs of subunits B and D only after it is already bound to at least the same number of LBDs that belong to subunits A and C. Our structures combined with single-channel recordings, molecular dynamics simulations and machine-learning analysis suggest that channel opening requires agonist binding to at least two LBDs. Conversely, agonist binding to all four LBDs does not guarantee maximal channel conductance and favours subconductance states O1 and O2, with O3 and O4 being rare and not captured structurally. The lack of subunit independence and low efficiency coupling of glutamate binding to channel opening underlie the gating of synaptic complexes to submaximal conductance levels, which provide a potential for upregulation of synaptic activity.
Assuntos
Receptores de Glutamato , Receptores Ionotrópicos de Glutamato , Ácido Glutâmico/metabolismo , Simulação de Dinâmica Molecular , Domínios Proteicos , Receptores de Glutamato/metabolismo , Receptores Ionotrópicos de Glutamato/metabolismoRESUMO
BACKGROUND: Deficiency of adenosine deaminase 2 (DADA2) results in heterogeneous manifestations including systemic vasculitis and red cell aplasia. The basis of different disease phenotypes remains incompletely defined. OBJECTIVE: We sought to further delineate disease phenotypes in DADA2 and define the mechanistic basis of ADA2 variants. METHODS: We analyzed the clinical features and ADA2 variants in 33 patients with DADA2. We compared the transcriptomic profile of 14 patients by bulk RNA sequencing. ADA2 variants were expressed experimentally to determine impact on protein production, trafficking, release, and enzymatic function. RESULTS: Transcriptomic analysis of PBMCs from DADA2 patients with the vasculitis phenotype or pure red cell aplasia phenotype exhibited similar upregulation of TNF, type I interferon, and type II interferon signaling pathways compared with healthy controls. These pathways were also activated in 3 asymptomatic individuals with DADA2. Analysis of ADA2 variants, including 7 novel variants, showed different mechanisms of functional disruption including (1) unstable transcript leading to RNA degradation; (2) impairment of ADA2 secretion because of retention in the endoplasmic reticulum; (3) normal expression and secretion of ADA2 that lacks enzymatic function; and (4) disruption of the N-terminal signal peptide leading to cytoplasmic localization of unglycosylated protein. CONCLUSIONS: Transcriptomic signatures of inflammation are observed in patients with different disease phenotypes, including some asymptomatic individuals. Disease-associated ADA2 variants affect protein function by multiple mechanisms, which may contribute to the clinical heterogeneity of DADA2.
Assuntos
Adenosina Desaminase , Vasculite , Humanos , Adenosina Desaminase/genética , Peptídeos e Proteínas de Sinalização Intercelular/genética , Fenótipo , MutaçãoRESUMO
A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.
RESUMO
In silico identification of potent protein inhibitors commonly requires prediction of a ligand binding free energy (BFE). Thermodynamics integration (TI) based on molecular dynamics (MD) simulations is a BFE calculation method capable of acquiring accurate BFE, but it is computationally expensive and time-consuming. In this work, we have developed an efficient automated workflow for identifying compounds with the lowest BFE among thousands of congeneric ligands, which requires only hundreds of TI calculations. Automated machine learning (AutoML) orchestrated by active learning (AL) in an AL-AutoML workflow allows unbiased and efficient search for a small set of best-performing molecules. We have applied this workflow to select inhibitors of the SARS-CoV-2 papain-like protease and were able to find 133 compounds with improved binding affinity, including 16 compounds with better than 100-fold binding affinity improvement. We obtained a hit rate that outperforms that expected of traditional expert medicinal chemist-guided campaigns. Thus, we demonstrate that the combination of AL and AutoML with free energy simulations provides at least 20× speedup relative to the naïve brute force approaches.
Assuntos
COVID-19 , Humanos , SARS-CoV-2/metabolismo , Desenho de Fármacos , Proteínas/química , Termodinâmica , Simulação de Dinâmica Molecular , Ligação Proteica , LigantesRESUMO
Zn2+, Mg2+ and Ca2+ are essential divalent cations implicated in many metabolic processes and signalling pathways. An emerging new paradigm is that the organismal balance of these cations predominantly depends on a common gatekeeper, the channel-kinase TRPM7. Despite extensive electrophysiological studies and recent cryo-EM analysis, an open question is how the channel activity of TRPM7 is activated. Here, we performed site-directed mutagenesis of mouse TRPM7 in conjunction with patch-clamp assessment of whole-cell and single-channel activity and molecular dynamics (MD) simulations to show that the side chains of conserved N1097 form an inter-subunit Mg2+ regulatory site located in the lower channel gate of TRPM7. Our results suggest that intracellular Mg2+ binds to this site and stabilizes the TRPM7 channel in the closed state, whereas the removal of Mg2+ favours the opening of TRPM7. Hence, our study identifies the structural underpinnings through which the TRPM7 channel is controlled by cytosolic Mg2+, representing a new structure-function relationship not yet explored among TRPM channels.
Assuntos
Canais de Cátion TRPM , Animais , Cátions Bivalentes/metabolismo , Magnésio/metabolismo , Camundongos , Fosfotransferases/metabolismo , Canais de Cátion TRPM/genética , Canais de Cátion TRPM/metabolismoRESUMO
Diamond-Blackfan anemia (DBA) is an inherited bone marrow failure syndrome, associated with mutations in ribosomal protein (RP) genes. Growing data on mutations in non-RP genes in patients with DBA-like phenotype became available over recent years. We describe two patients with the phenotype of DBA (onset of macrocytic anemia within the first year of life, paucity of erythroid precursors in bone marrow) and germline de novo variants in the TP53 gene. Both patients became transfusion independent, probably due to L-leucine therapy. The possible role of TP53 variants should be considered in patients with DBA-like phenotype and no mutations in RP genes.
Assuntos
Anemia de Diamond-Blackfan , Anemia de Diamond-Blackfan/genética , Anemia de Diamond-Blackfan/terapia , Células Germinativas , Humanos , Mutação , Fenótipo , Proteínas Ribossômicas/genética , Proteína Supressora de Tumor p53/genéticaRESUMO
Epithelial calcium channel TRPV6 is a member of the vanilloid subfamily of TRP channels that is permeable to cations and highly selective to Ca2+ ; it shows constitutive activity regulated negatively by Ca2+ and positively by phosphoinositol and cholesterol lipids. In this review, we describe the molecular structure of TRPV6 and discuss how its structural elements define its unique functional properties. High Ca2+ selectivity of TRPV6 originates from the narrow selectivity filter, where Ca2+ ions are directly coordinated by a ring of anionic aspartate side chains. Divalent cations Ca2+ and Ba2+ permeate TRPV6 pore according to the knock-off mechanism, while tight binding of Gd3+ to the aspartate ring blocks the channel and prevents Na+ from permeating the pore. The iris-like channel opening is accompanied by an α-to-π helical transition in the pore-lining transmembrane helix S6. As a result of this transition, the intracellular halves of the S6 helices bend and rotate by about 100 deg, exposing different residues to the channel pore in the open and closed states. Channel opening is also associated with changes in occupancy of the transmembrane domain lipid binding sites. The inhibitor 2-aminoethoxydiphenyl borate (2-APB) binds to TRPV6 in a pocket formed by the cytoplasmic half of the S1-S4 transmembrane helical bundle and shifts open-closed channel equilibrium towards the closed state by outcompeting lipids critical for activation. Ca2+ inhibits TRPV6 via binding to calmodulin (CaM), which mediates Ca2+ -dependent inactivation. The TRPV6-CaM complex exhibits 1:1 stoichiometry; one TRPV6 tetramer binds both CaM lobes, which adopt a distinct head-to-tail arrangement. The CaM C-terminal lobe plugs the channel through a unique cation-π interaction by inserting the side chain of lysine K115 into a tetra-tryptophan cage at the ion channel pore intracellular entrance. Recent studies of TRPV6 structure and function described in this review advance our understanding of the role of this channel in physiology and pathophysiology and inform new therapeutic design.
Assuntos
Canais de Cálcio , Cálcio , Sítios de Ligação , Cálcio/metabolismo , Canais de Cálcio/metabolismo , Calmodulina/metabolismo , Canais de Cátion TRPV/metabolismoRESUMO
Mechanisms of enzymatic epoxidation via oxygen atom transfer (OAT) to an olefin moiety is mainly derived from the studies on thiolate-heme containing epoxidases, such as cytochrome P450 epoxidases. The molecular basis of epoxidation catalyzed by nonheme-iron enzymes is much less explored. Herein, we present a detailed study on epoxidation catalyzed by the nonheme iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase, AsqJ. The native substrate and analogues with different para substituents ranging from electron-donating groups (e.g., methoxy) to electron-withdrawing groups (e.g., trifluoromethyl) were used to probe the mechanism. The results derived from transient-state enzyme kinetics, Mössbauer spectroscopy, reaction product analysis, X-ray crystallography, density functional theory calculations, and molecular dynamic simulations collectively revealed the following mechanistic insights: (1) The rapid O2 addition to the AsqJ Fe(II) center occurs with the iron-bound 2OG adopting an online-binding mode in which the C1 carboxylate group of 2OG is trans to the proximal histidine (His134) of the 2-His-1-carboxylate facial triad, instead of assuming the offline-binding mode with the C1 carboxylate group trans to the distal histidine (His211); (2) The decay rate constant of the ferryl intermediate is not strongly affected by the nature of the para substituents of the substrate during the OAT step, a reactivity behavior that is drastically different from nonheme Fe(IV)-oxo synthetic model complexes; (3) The OAT step most likely proceeds through a stepwise process with the initial formation of a C(benzylic)-O bond to generate an Fe-alkoxide species, which is observed in the AsqJ crystal structure. The subsequent C3-O bond formation completes the epoxide installation.
Assuntos
Aspergillus nidulans/metabolismo , Compostos de Epóxi/metabolismo , Proteínas Fúngicas/metabolismo , Ácidos Cetoglutáricos/metabolismo , Oxigênio/metabolismo , Oxigenases/metabolismo , Aspergillus nidulans/química , Aspergillus nidulans/enzimologia , Cristalografia por Raios X , Compostos de Epóxi/química , Proteínas Fúngicas/química , Ferro/química , Ferro/metabolismo , Modelos Moleculares , Oxigênio/química , Oxigenases/químicaRESUMO
The dependency of current-voltage characteristics of the α-hemolysin channel on the channel position within the membrane was studied using Poisson-Nernst-Planck theory of ion conductivity with soft repulsion between mobile ions and protein atoms (SP-PNP). The presence of the membrane environment also influences the protonation state of the residues at the boundary of the water-lipid interface. In this work, we predict that Asp and Lys residues at the protein rim change their protonation state upon penetration to the lipid environment. Free energies of protein insertion in the membrane for different penetration depths were estimated using the Poisson-Boltzmann/solvent-accessible surface area (PB/SASA) model. The results show that rectification and reversal potentials are very sensitive to the relative position of channel in the membrane, which in turn contributes to alternative protonation states of lipid-penetrating ionizable groups. The prediction of channel position based on the matching of calculated rectification with experimentally determined rectification is in good agreement with recent neutron reflection experiments. Based on the results, we conclude that α-hemolysin membrane position is determined by a combination of factors and not only by the pattern of the surface hydrophobicity as is typically assumed.
Assuntos
Proteínas Hemolisinas/metabolismo , Canais Iônicos/metabolismo , Proteínas Hemolisinas/química , Interações Hidrofóbicas e Hidrofílicas , Canais Iônicos/química , Potenciais da Membrana/fisiologia , Modelos Moleculares , Modelos TeóricosRESUMO
AsqJ, an iron(II)- and 2-oxoglutarate-dependent enzyme found in viridicatin-type alkaloid biosynthetic pathways, catalyzes sequential desaturation and epoxidation to produce cyclopenins. Crystal structures of AsqJ bound to cyclopeptin and its C3â epimer are reported. Meanwhile, a detailed mechanistic study was carried out to decipher the desaturation mechanism. These findings suggest that a pathway involving hydrogen atom abstraction at the C10 position of the substrate by a short-lived FeIV -oxo species and the subsequent formation of a carbocation or a hydroxylated intermediate is preferred during AsqJ-catalyzed desaturation.
Assuntos
Compostos de Epóxi/metabolismo , Proteínas Fúngicas/metabolismo , Peptídeos/metabolismo , Aspergillus nidulans/enzimologia , Biocatálise , Domínio Catalítico , Sistema Enzimático do Citocromo P-450/metabolismo , Compostos de Epóxi/química , Compostos Férricos/química , Proteínas Fúngicas/química , Ácidos Cetoglutáricos/química , Ácidos Cetoglutáricos/metabolismo , Simulação de Dinâmica Molecular , Peptídeos/química , Teoria Quântica , EstereoisomerismoRESUMO
Ionotropic glutamate receptors are a family of tetrameric ion channels with functional states consisting of nonconducting, conducting, and desensitized states that are starting to become well characterized by electrophysiological and biophysical studies. However, the structure and relative energetics of these states beyond the general structure of the receptor are still not well understood. It is known that the interface between monomeric subunits of the tetramer plays a major role in distinguishing these functional states. We have used umbrella sampling and multimicrosecond molecular dynamics simulations of the GluA2 AMPA subtype glutamate receptor ligand-binding domain (LBD) dimers to characterize a natural propensity of the LBD dimers for various configurational states. Our results show a proposed desensitized conformation of the LBD dimer is a highly preferable conformation of the LBD dimer without the influence of other receptor domains or crystallographic conditions. This has been demonstrated by both free protein simulations of 5 µs duration, as well as by computed free energy difference between the active and desensitized states. At the same time, the simulations performed using the same protocols revealed that for the LBD mutant L483Y, known to lack desensitization, the postulated active state of the LBD dimer is indeed the preferred configurational state, which remained stable in the simulations. Our findings pave the path for developing more detailed hypotheses of the full receptor activation mechanism. Combined with the energetics of glutamate binding to the LBD and the energy required to open the transmembrane pore helices, our results strongly support a hypothesis that the low absolute free-energy state is the desensitized state of the intact AMPA receptor.
Assuntos
Receptores de AMPA/metabolismo , Ácido Glutâmico/química , Ácido Glutâmico/metabolismo , Simulação de Dinâmica Molecular , Mutação , Ligação Proteica , Domínios Proteicos , Multimerização Proteica , Estabilidade Proteica , Receptores de AMPA/química , Receptores de AMPA/genética , Soluções , Água/químicaRESUMO
In the present studies, we analyzed the influence of temperature on the stability and dynamics of the α subunit of tryptophan synthase (TRPS) from hyperthermophilic, mesophilic, and psychrophilic homologues at different temperatures by molecular dynamics simulations. Employing different indicators such as root-mean-square deviations, root-mean-square fluctuations, principal component analysis, and free energy landscapes, this study manifests the diverse behavior of these homologues with changes in temperature. Especially, an enhancement in the collective motions, classified as representative motions, is observed at high temperature. Similarly, the criterion for the selection of electrostatic interactions in terms of their life span (duty cycle) has indeed helped in identifying the short- and long-lived electrostatic interactions and how they affect the protein's overall stability at different temperatures. Rigidity and flexibility patterns of the homologous proteins are examined using FIRST software along with the calculation of duty cycles with various threshold limits at different temperatures. Rigid cluster decomposition in TRPS of psychrophilic, mesophilic, and hyperthermophilic origin identifies the flexible and rigid regions in the protein. Early loss of rigidity is observed in mesophilic TRPS via loss of contact between the major fragments of the protein compared with the other homologues. In spite of the high similarity of their three-dimensional structures, the overall responses of the three proteins to varying temperatures are significantly different.
Assuntos
Simulação de Dinâmica Molecular , Homologia de Sequência de Aminoácidos , Temperatura , Triptofano Sintase/química , Estabilidade Enzimática , Ligação de Hidrogênio , Interações Hidrofóbicas e Hidrofílicas , Conformação Proteica , Subunidades Proteicas/químicaRESUMO
Diphtheria toxin translocation (T) domain inserts in lipid bilayers upon acidification of the environment. Computational and experimental studies have suggested that low pH triggers a conformational change of the T-domain in solution preceding membrane binding. The refolded membrane-competent state was modeled to be compact and mostly retain globular structure. In the present work, we investigate how this refolded state interacts with membrane interfaces in the early steps of T-domain's membrane association. Coarse-grained molecular dynamics simulations suggest two distinct membrane-bound conformations of the T-domain in the presence of bilayers composed of a mixture of zwitteronic and anionic phospholipids (POPC:POPG with a 1:3 molar ratio). Both membrane-bound conformations show a common near parallel orientation of hydrophobic helices TH8-TH9 relative to the membrane plane. The most frequently observed membrane-bound conformation is stabilized by electrostatic interactions between the N-terminal segment of the protein and the membrane interface. The second membrane-bound conformation is stabilized by hydrophobic interactions between protein residues and lipid acyl chains, which facilitate deeper protein insertion in the membrane interface. A theoretical estimate of a free energy of binding of a membrane-competent T-domain to the membrane is provided.
Assuntos
Toxina Diftérica/química , Bicamadas Lipídicas/química , Histidina/química , Concentração de Íons de Hidrogênio , Interações Hidrofóbicas e Hidrofílicas , Simulação de Dinâmica Molecular , Fosfatidilcolinas/química , Estrutura Secundária de ProteínaRESUMO
BACKGROUND: Diamond Blackfan anemia (DBA) is a genetically and clinically heterogeneous ribosomopathy and inherited bone marrow failure syndrome characterized by anemia, reticulocytopenia, and decreased erythroid precursors in the bone marrow with an increased risk of malignancy and, in approximately 50%, physical abnormalities. METHODS: We retrospectively analyzed clinical data from 77 patients with DBA born in the Russian Federation from 1993 to 2014. In 74 families there was one clinically affected individual; in only three instances a multiplex family was identified. Genomic DNA from 57 DBA patients and their first-degree relatives was sequenced for mutations in RPS19, RPS10, RPS24, RPS26, RPS7, RPS17, RPL5, RPL11, RPL35a, and GATA1. RESULTS: Severe anemia presented before 8 months of age in all 77 patients; before 2 months in 61 (78.2%); before 4 months in 71 (92.2%). Corticosteroid therapy was initiated after 1 year of age in the majority of patients. Most responded initially to steroids, while 5 responses were transient. Mutations in RP genes were detected in 35 of 57 patients studied: 15 in RPS19, 6 in RPL5, 3 in RPS7, 3 each in RPS10, RPS26, and RPL11 and 1 each in RPS24 and RPL35a; 24 of these mutations have not been previously reported. One patient had a balanced chromosomal translocation involving RPS19. No mutations in GATA1 were found. CONCLUSION: In our cohort from an ethnically diverse population the distribution of mutations among RP genes was approximately the same as was reported by others, although within genotypes most of the mutations had not been previously reported.
Assuntos
Anemia de Diamond-Blackfan/genética , Fator de Transcrição GATA1/genética , Mutação , Proteínas Ribossômicas/genética , Anormalidades Múltiplas/genética , Adolescente , Anemia de Diamond-Blackfan/epidemiologia , Criança , Pré-Escolar , Anormalidades Craniofaciais/genética , Análise Mutacional de DNA , Feminino , Fator de Transcrição GATA1/deficiência , Heterogeneidade Genética , Genótipo , Cardiopatias Congênitas/genética , Humanos , Lactente , Masculino , Fenótipo , Estudos Retrospectivos , Proteínas Ribossômicas/deficiência , Federação Russa/epidemiologia , Análise de Sequência de DNA , Adulto JovemRESUMO
Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile-mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.
Assuntos
Conformação Proteica , Ligação de Hidrogênio , Modelos Moleculares , Simulação de Dinâmica MolecularRESUMO
TRP channels are implicated in various diseases, but high structural similarity between them makes selective pharmacological modulation challenging. Here, we study the molecular mechanism underlying specific inhibition of the TRPM7 channel, which is essential for cancer cell proliferation, by the anticancer agent CCT128930 (CCT). Using cryo-EM, functional analysis, and MD simulations, we show that CCT binds to a vanilloid-like (VL) site, stabilizing TRPM7 in the closed non-conducting state. Similar to other allosteric inhibitors of TRPM7, NS8593 and VER155008, binding of CCT is accompanied by displacement of a lipid that resides in the VL site in the apo condition. Moreover, we demonstrate the principal role of several residues in the VL site enabling CCT to inhibit TRPM7 without impacting the homologous TRPM6 channel. Hence, our results uncover the central role of the VL site for the selective interaction of TRPM7 with small molecules that can be explored in future drug design.
Assuntos
1-Naftilamina/análogos & derivados , Antineoplásicos , Canais de Cátion TRPM , Canais de Cátion TRPM/metabolismo , Canais de Cátion TRPM/antagonistas & inibidores , Humanos , Antineoplásicos/farmacologia , Antineoplásicos/química , Proteínas Serina-Treonina Quinases/metabolismo , Proteínas Serina-Treonina Quinases/antagonistas & inibidores , Células HEK293 , Simulação de Dinâmica Molecular , Sítios de Ligação , Ligação Proteica , Microscopia CrioeletrônicaRESUMO
The Critical Assessment of Computational Hit-Finding Experiments (CACHE) Challenge series is focused on identifying small molecule inhibitors of protein targets using computational methods. Each challenge contains two phases, hit-finding and follow-up optimization, each of which is followed by experimental validation of the computational predictions. For the CACHE Challenge #1, the Leucine-Rich Repeat Kinase 2 (LRRK2) WD40 Repeat (WDR) domain was selected as the target for in silico hit-finding and optimization. Mutations in LRRK2 are the most common genetic cause of the familial form of Parkinson's disease. The LRRK2 WDR domain is an understudied drug target with no known molecular inhibitors. Herein we detail the first phase of our winning submission to the CACHE Challenge #1. We developed a framework for the high-throughput structure-based virtual screening of a chemically diverse small molecule space. Hit identification was performed using the large-scale Deep Docking (DD) protocol followed by absolute binding free energy (ABFE) simulations. ABFEs were computed using an automated molecular dynamics (MD)-based thermodynamic integration (TI) approach. 4.1 billion ligands from Enamine REAL were screened with DD followed by ABFEs computed by MD TI for 793 ligands. 76 ligands were prioritized for experimental validation, with 59 compounds successfully synthesized and 5 compounds identified as hits, yielding a 8.5% hit rate. Our results demonstrate the efficacy of the combined DD and ABFE approaches for hit identification for a target with no previously known hits. This approach is widely applicable for the efficient screening of ultra-large chemical libraries as well as rigorous protein-ligand binding affinity estimation leveraging modern computational resources.
RESUMO
The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 molecules, which were subsequently ranked to find 'consensus compounds'. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for biological activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (only the Nsp12 domain), and (alpha) spike protein S. Overall, 27 compounds with weak inhibition/binding were experimentally identified by binding-, cleavage-, and/or viral suppression assays and are presented here. Open science approaches such as the one presented here contribute to the knowledge base of future drug discovery efforts in finding better SARS-CoV-2 treatments.
Assuntos
COVID-19 , SARS-CoV-2 , Humanos , Pandemias , Bioensaio , Descoberta de DrogasRESUMO
The transient receptor potential channel TRPM7 is a master regulator of the organismal balance of divalent cations that plays an essential role in embryonic development, immune responses, cell mobility, proliferation, and differentiation. TRPM7 is implicated in neuronal and cardiovascular disorders, tumor progression and has emerged as a new drug target. Here we use cryo-EM, functional analysis, and molecular dynamics simulations to uncover two distinct structural mechanisms of TRPM7 activation by a gain-of-function mutation and by the agonist naltriben, which show different conformational dynamics and domain involvement. We identify a binding site for highly potent and selective inhibitors and show that they act by stabilizing the TRPM7 closed state. The discovered structural mechanisms provide foundations for understanding the molecular basis of TRPM7 channelopathies and drug development.