RESUMO
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
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
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 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
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
Introduction: Nephroblastoma (Wilms tumor (WT)) is an embryonal tumor accounting for >90% of pediatric renal cancers. About 10% of WTs harbor pathogenic germline mutations. The REST gene, classified as a putative tumor suppressor, is affected in 2% of WTs. High-throughput molecular methods facilitate advanced diagnostics of cancer. In addition to this, germline mutations in REST are also associated with familial gingival fibromatosis (GFM). Reciprocally, none of the articles on RESTmut WT mentions GFM as a comorbid condition. This report provides unique evidence on the WT-GFM comorbidity in RESTmut carriers. Case presentation: Patient 1 (a 5-year-old boy with unilateral WT) is a proband, who has two healthy siblings. Patient 2 (a 4-year-old girl with bilateral WT) is a proband from in vitro fertilization (IVF) triplets, with a sister and brother without WT. We analyzed probands' DNA extracted from peripheral blood leucocytes with a custom-targeted next-generation sequencing (NGS)-198 gene panel. The detected variants were checked in family members by Sanger sequencing. Patient 1 had a pathogenic germline mutation in REST: c.1035_1036insTA, p.(E346*), as did his mother and both brothers. There were two other WT cases in this family (proband's maternal uncles). Patient 2 had a pathogenic germline variant in REST: c.2668_2671del, p.(E891Pfs*6), as well as her sister. The mutation was probably inherited from their deceased father, as he had gingival fibromatosis. Family members with REST mutations from both families had gingival fibromatosis. A somatic REST c.663C>A p.C221* mutation was identified in one patient with WT. At the moment both patients with WT are under dynamic observation without signs of the disease. Conclusion: Here, we describe two clinical cases of WT in nonrelated young children with germline-inactivating REST variants identified by next-generation sequencing. Both patients present with familial gingival fibromatosis, regarded as clinically useful comorbidity indicative of the tumor predisposition syndrome. The two cases illustrate Wilms tumor-gingival fibromatosis comorbidity in carriers of germline-inactivated REST alleles previously identified as a predisposition factor for both conditions.
RESUMO
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
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.
Assuntos
Canais de Cátion TRPM , Canais de Cátion TRPM/metabolismo , Diferenciação CelularRESUMO
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
Protein-ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)âa multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein-ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of â¼0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.
Assuntos
Simulação de Dinâmica Molecular , Proteínas , Ligantes , Ligação Proteica , Entropia , Conformação Molecular , Proteínas/química , TermodinâmicaRESUMO
N-methyl-D-aspartate receptors (NMDARs) are transmembrane proteins that are activated by the neurotransmitter glutamate and are found at most excitatory vertebrate synapses. NMDAR channel blockers, an antagonist class of broad pharmacological and clinical significance, inhibit by occluding the NMDAR ion channel. A vast literature demonstrates that NMDAR channel blockers, including MK-801, phencyclidine, ketamine, and the Alzheimer's disease drug memantine, can bind and unbind only when the NMDAR channel is open. Here we use electrophysiological recordings from transfected tsA201 cells and cultured neurons, NMDAR structural modeling, and custom-synthesized compounds to show that NMDAR channel blockers can enter the channel through two routes: the well-known hydrophilic path from extracellular solution to channel through the open channel gate, and also a hydrophobic path from plasma membrane to channel through a gated fenestration ("membrane-to-channel inhibition" (MCI)). Our demonstration that ligand-gated channels are subject to MCI, as are voltage-gated channels, highlights the broad expression of this inhibitory mechanism.
Assuntos
Receptores de N-Metil-D-Aspartato , Maleato de Dizocilpina , Canais Iônicos , Ketamina/farmacologia , Memantina/farmacologia , Receptores de N-Metil-D-Aspartato/antagonistas & inibidores , Receptores de N-Metil-D-Aspartato/metabolismoRESUMO
Currently, of the few accessible symptomatic therapies for Alzheimer's disease (AD), memantine is the only N-methyl-d-aspartate receptor (NMDAR) blocker approved by the FDA. This work further explores a series of memantine analogs featuring a benzohomoadamantane scaffold. Most of the newly synthesized compounds block NMDARs in the micromolar range, but with lower potency than previously reported hit IIc, results that were supported by molecular dynamics simulations. Subsequently, electrophysiological studies with the more potent compounds allowed classification of IIc, a low micromolar, uncompetitive, voltage-dependent, NMDAR blocker, as a memantine-like compound. The excellent in vitro DMPK properties of IIc made it a promising candidate for in vivo studies in Caenorhabditis elegans (C. elegans) and in the 5XFAD mouse model of AD. Administration of IIc or memantine improved locomotion and rescues chemotaxis behavior in C. elegans. Furthermore, both compounds enhanced working memory in 5XFAD mice and modified NMDAR and CREB signaling, which may prevent synaptic dysfunction and modulate neurodegenerative progression.
Assuntos
Doença de Alzheimer , Memantina , Doença de Alzheimer/tratamento farmacológico , Animais , Caenorhabditis elegans , Modelos Animais de Doenças , Memantina/farmacologia , Camundongos , Receptores de N-Metil-D-AspartatoRESUMO
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
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
The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools' computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.
RESUMO
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
BACKGROUND AND PURPOSE: AMPA receptors, which shape excitatory postsynaptic currents and are directly involved in overactivation of synaptic function during seizures, represent a well-accepted target for anti-epileptic drugs. Trans-4-butylcyclohexane carboxylic acid (4-BCCA) has emerged as a new promising anti-epileptic drug in several in vitro and in vivo seizure models, but the mechanism of its action remained unknown. The purpose of this study is to characterize structure and dynamics of 4-BCCA interaction with AMPA receptors. EXPERIMENTAL APPROACH: We studied the molecular mechanism of AMPA receptor inhibition by 4-BCCA using a combination of X-ray crystallography, mutagenesis, electrophysiological assays, and molecular dynamics simulations. KEY RESULTS: We identified 4-BCCA binding sites in the transmembrane domain (TMD) of AMPA receptor, at the lateral portals formed by transmembrane segments M1-M4. At this binding site, 4-BCCA is very dynamic, assumes multiple poses, and can enter the ion channel pore. CONCLUSION AND IMPLICATIONS: 4-BCCA represents a low-affinity inhibitor of AMPA receptors that acts at the TMD sites distinct from non-competitive inhibitors, such as the anti-epileptic drug perampanel and the ion channel blockers. Further studies might examine the possibsility of synergistic use of these inhibitors in treatment of epilepsy and a wide range of neurological disorders and gliomas. LINKED ARTICLES: This article is part of a themed issue on Structure Guided Pharmacology of Membrane Proteins (BJP 75th Anniversary). To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.14/issuetoc.
Assuntos
Ácidos Carboxílicos , Receptores de AMPA , Cicloexanos , Humanos , Receptores de AMPA/metabolismo , ConvulsõesRESUMO
Background. This study deals with a rare (orphan) monogenic connective tissue disorder - Ehlers-Danlos syndrome kyphoscoliotic type 2 (EDSKS2). Kyphoscoliotic type 2 Ehlers-Danlos syndrome is an autosomal recessive disorder caused by mutations in the FKBP14 gene (7p14.3), which encodes the FKBP22 protein. According to the 2017 classification, this type is in group seven - collagen spatial structure and cross-linking defects. We present results of clinical examination and molecular genetic analysis for five patients with age varying from two to fifteen years. Methods. Five patients were examined using clinical and laboratory methods. DNA samples used for the analysis were extracted from whole blood samples using a Wizard® Genomic DNA Purification Kit (Promega, USA) according to the manufacturer's protocol. Results. The major clinical findings were kyphoscoliosis, early motor development delay, muscular weakness, hypotonia and hearing loss. Molecular genetic analysis detected a homozygous c.362dupC duplication in exon 3 of the FKBP14 gene in all five patients. This mutation is common in various countries. Differential diagnostics were carried out to exclude other Ehlers-Danlos syndrome types and myopathies. Conclusions. Literature analysis and examination of five EDSKS2 patients demonstrated the involvement of major organs and systems, such as joints, spine, muscles, cardiovascular system, respiratory system, hearing, and vision, into the pathological process. Kidney mobility increases and nephroptosis seems to be secondary caused by muscular weakness. During molecular genetic analysis, to verify EDSKS2 it is recommended to initially search for the c.362dupC duplication, which appears to be common in European countries, including Russia.
Assuntos
Síndrome de Ehlers-Danlos , Adolescente , Criança , Pré-Escolar , Síndrome de Ehlers-Danlos/diagnóstico , Síndrome de Ehlers-Danlos/genética , Éxons , Humanos , Mutação , Peptidilprolil Isomerase/genética , FenótipoRESUMO
This work measures baseline sampling characteristics that highlight fundamental differences between sampling methods for assembly driven by short-ranged pair potentials. Such granular comparison is essential for fast, flexible, and accurate hybridization of complementary methods. Besides sampling speed, efficiency, and accuracy of uniform grid coverage, other sampling characteristics measured are (i) accuracy of covering narrow low energy regions that have low effective dimension (ii) ability to localize sampling to specific basins, and (iii) flexibility in sampling distributions. As a proof of concept, we compare a recently developed geometric methodology EASAL (Efficient Atlasing and Search of Assembly Landscapes) and the traditional Monte Carlo (MC) method for sampling the energy landscape of two assembling trans-membrane helices, driven by short-range pair potentials. By measuring the above-mentioned sampling characteristics, we demonstrate that EASAL provides localized and accurate coverage of crucial regions of the energy landscape of low effective dimension, under flexible sampling distributions, with much fewer samples and computational resources than MC sampling. EASAL's empirically validated theoretical guarantees permit credible extrapolation of these measurements and comparisons to arbitrary number and size of assembling units. Promising avenues for hybridizing the complementary advantages of the two methods are discussed.
RESUMO
Nature has developed complexity-generating reactions within natural product biosynthetic pathways. However, direct utilization of these pathways to prepare compound libraries remains challenging due to limited substrate scopes, involvement of multiple-step reactions, and moderate robustness of these sophisticated enzymatic transformations. Synthetic chemistry, on the other hand, offers an alternative approach to prepare natural product analogs. However, owing to complex and diverse functional groups appended on the targeted molecules, dedicated design and development of synthetic strategies are typically required. Herein, by leveraging the power of chemo-enzymatic synthesis, we report an approach to bridge the gap between biological and synthetic strategies in the preparation of quinolone alkaloid analogs. Leading by in silico analysis, the predicted substrate analogs were chemically synthesized. The AsqJ-catalyzed asymmetric epoxidation of these substrate analogues was followed by an Lewis Acid-triggered ring contraction to complete the viridicatin formation. We evaluated the robustness of this method in gram-scale reactions. Lastly, through chemoenzymatic cascades, a library of quinolone alkaloids is effectively prepared.