RESUMO
Chemical structure optimization is a vital part of early drug discovery projects. Starting with compounds that show activity on the target of interest, the chemical structures are subsequently optimized toward a development candidate (DC) molecule with the best chances of clinical success. However, the DCs in the context of such optimization programs, as well as detailed characterization of major limiting factors, have not been investigated in detail so far. Here, we report an analysis of the historical DC molecules at Novartis since 2005 in the context of their optimization projects. Mapping the DCs into their respective chemical optimization series, we find that these tend to be synthesized rather early in a substantial number of cases. Further analysis of structural properties, ADMET, and potency-related readouts revealed that DC compounds tend to be generally significantly smaller, more permeable, and have higher ligand efficiency than other compounds sent to in vivo PK studies, which we also show for compounds from the same chemical series. Although this might seem obvious to most practitioners in medicinal chemistry, for all of these properties, we could show that they tend to evolve in an undesired direction during structure optimization. This highlights the difficulty of successfully translating our knowledge to medicinal chemistry optimizations.
RESUMO
Time-split cross-validation is broadly recognized as the gold standard for validating predictive models intended for use in medicinal chemistry projects. Unfortunately this type of data is not broadly available outside of large pharmaceutical research organizations. Here we introduce the SIMPD (simulated medicinal chemistry project data) algorithm to split public data sets into training and test sets that mimic the differences observed in real-world medicinal chemistry project data sets. SIMPD uses a multi-objective genetic algorithm with objectives derived from an extensive analysis of the differences between early and late compounds in more than 130 lead-optimization projects run within the Novartis Institutes for BioMedical Research. Applying SIMPD to the real-world data sets produced training/test splits which more accurately reflect the differences in properties and machine-learning performance observed for temporal splits than other standard approaches like random or neighbor splits. We applied the SIMPD algorithm to bioactivity data extracted from ChEMBL and created 99 public data sets which can be used for validating machine-learning models intended for use in the setting of a medicinal chemistry project. The SIMPD code and simulated data sets are available under open-source/open-data licenses at github.com/rinikerlab/molecular_time_series.
RESUMO
Early in silico assessment of the potential of a series of compounds to deliver a drug is one of the major challenges in computer-assisted drug design. The goal is to identify the right chemical series of compounds out of a large chemical space to then subsequently prioritize the molecules with the highest potential to become a drug. Although multiple approaches to assess compounds have been developed over decades, the quality of these predictors is often not good enough and compounds that agree with the respective estimates are not necessarily druglike. Here, we report a novel deep learning approach that leverages large-scale predictions of â¼100 ADMET assays to assess the potential of a compound to become a relevant drug candidate. The resulting score, which we termed bPK score, substantially outperforms previous approaches and showed strong discriminative performance on data sets where previous approaches did not.
Assuntos
Simulação por ComputadorRESUMO
In the drug development process, optimization of properties and biological activities of small molecules is an important task to obtain drug candidates with optimal efficacy when first applied in subsequent clinical studies. However, despite its importance, large-scale investigations of the optimization process in early drug discovery are lacking, likely due to the absence of historical records of different chemical series used in past projects. Here, we report a retrospective reconstruction of â¼3000 chemical series from the Novartis compound database, which allows us to characterize the general properties of chemical series as well as the time evolution of structural properties, ADMET properties, and target activities. Our data-driven approach allows us to substantiate common MedChem knowledge. We find that size, fraction of sp3-hybridized carbon atoms (Fsp3), and the density of stereocenters tend to increase during optimization, while the aromaticity of the compounds decreases. On the ADMET side, solubility tends to increase and permeability decreases, while safety-related properties tend to improve. Importantly, while ligand efficiency decreases due to molecular growth over time, target activities and lipophilic efficiency tend to improve. This emphasizes the heavy-atom count and logâ¯D as important parameters to monitor, especially as we further show that the decrease in permeability can be explained with the increase in molecular size. We highlight overlaps, shortcomings, and differences of the computationally reconstructed chemical series compared to the series used in recent internal drug discovery projects and investigate the relation to historical projects.
Assuntos
Descoberta de Drogas , Estudos Retrospectivos , Ligantes , Solubilidade , Bases de Dados FactuaisRESUMO
In electron cryomicroscopy (cryo-EM), molecular images of vitrified biological samples are obtained by conventional transmission microscopy (CTEM) using large underfocuses and subsequently computationally combined into a high-resolution three-dimensional structure. Here, we apply scanning transmission electron microscopy (STEM) using the integrated differential phase contrast mode also known as iDPC-STEM to two cryo-EM test specimens, keyhole limpet hemocyanin (KLH) and tobacco mosaic virus (TMV). The micrographs show complete contrast transfer to high resolution and enable the cryo-EM structure determination for KLH at 6.5 Å resolution, as well as for TMV at 3.5 Å resolution using single-particle reconstruction methods, which share identical features with maps obtained by CTEM of a previously acquired same-sized TMV data set. These data show that STEM imaging in general, and in particular the iDPC-STEM approach, can be applied to vitrified single-particle specimens to determine near-atomic resolution cryo-EM structures of biological macromolecules.
Assuntos
Microscopia Crioeletrônica , Microscopia Crioeletrônica/métodos , Microscopia Eletrônica de Transmissão e VarreduraRESUMO
The productivity of single-particle cryo-EM as a structure determination method has rapidly increased as many novel biological structures are being elucidated. The ultimate result of the cryo-EM experiment is an atomic model that should faithfully represent the computed image reconstruction. Although the principal approach of atomic model building and refinement from maps resembles that of the X-ray crystallographic methods, there are important differences due to the unique properties resulting from the 3D image reconstructions. In this review, we discuss the practiced work-flow from the cryo-EM image reconstruction to the atomic model. We give an overview of (i) resolution determination methods in cryo-EM including local and directional resolution variation, (ii) cryo-EM map contrast optimization including complementary map types that can help in identifying ambiguous density features, (iii) atomic model building and (iv) refinement in various resolution regimes including (v) their validation and (vi) discuss differences between X-ray and cryo-EM maps. Based on the methods originally developed for X-ray crystallography, the path from 3D image reconstruction to atomic coordinates has become an integral and important part of the cryo-EM structure determination work-flow that routinely delivers atomic models.
Assuntos
Microscopia Crioeletrônica/métodos , Glicosídeo Hidrolases/química , Cristalografia por Raios X , Processamento de Imagem Assistida por Computador , Imageamento Tridimensional , Modelos Moleculares , Peptídeos/química , Conformação Proteica , Dobramento de ProteínaRESUMO
Fourier shell correlation (FSC) has become a standard quantity for resolution estimation in electron cryo-microscopy. However, the resolution determination step is still subjective and not fully automated as it involves a series of map interventions before FSC computation and includes the selection of a common threshold. Here, we apply the statistical methods of permutation testing and false discovery rate (FDR) control to the resolution-dependent correlation measure. The approach allows fully automated and mask-free resolution determination based on statistical thresholding of FSC curves. We demonstrate the applicability for global, local and directional resolution estimation and show that the developed criterion termed FDR-FSC gives realistic resolution estimates based on statistical significance while eliminating the need of any map manipulations. The algorithms are implemented in a user-friendly GUI based software tool termed SPoC (https://github.com/MaximilianBeckers/SPOC).
Assuntos
Microscopia Crioeletrônica/métodos , Algoritmos , SoftwareRESUMO
Confidence maps provide complementary information for interpreting cryo-EM densities as they indicate statistical significance with respect to background noise. They can be thresholded by specifying the expected false-discovery rate (FDR), and the displayed volume shows the parts of the map that have the corresponding level of significance. Here, the basic statistical concepts of confidence maps are reviewed and practical guidance is provided for their interpretation and usage inside the CCP-EM suite. Limitations of the approach are discussed and extensions towards other error criteria such as the family-wise error rate are presented. The observed map features can be rendered at a common isosurface threshold, which is particularly beneficial for the interpretation of weak and noisy densities. In the current article, a practical guide is provided to the recommended usage of confidence maps.
Assuntos
Microscopia Crioeletrônica/métodos , ATPases Bacterianas Próton-Translocadoras/química , Gráficos por Computador , Modelos Moleculares , Conformação Proteica , Ribossomos/química , Eletricidade Estática , Estatística como Assunto , Vírus do Mosaico do Tabaco , Interface Usuário-ComputadorRESUMO
Stable capsid structures of viruses protect viral RNA while they also require controlled disassembly for releasing the viral genome in the host cell. A detailed understanding of viral disassembly processes and the involved structural switches is still lacking. This process has been extensively studied using tobacco mosaic virus (TMV), and carboxylate interactions are assumed to play a critical part in this process. Here, we present two cryo-EM structures of the helical TMV assembly at 2.0 and 1.9 Å resolution in conditions of high Ca2+ concentration at low pH and in water. Based on our atomic models, we identify the conformational details of the disassembly switch mechanism: In high Ca2+ /acidic pH environment, the virion is stabilized between neighboring subunits through carboxyl groups E95 and E97 in close proximity to a Ca2+ binding site that is shared between two subunits. Upon increase in pH and lower Ca2+ levels, mutual repulsion of the E95/E97 pair and Ca2+ removal destabilize the network of interactions between adjacent subunits at lower radius and release the switch for viral disassembly.
Assuntos
Modelos Moleculares , Vírus do Mosaico do Tabaco/fisiologia , Montagem de Vírus , Cálcio/química , Capsídeo/química , Concentração de Íons de Hidrogênio , Fenótipo , Reprodutibilidade dos Testes , VírionRESUMO
Cryo-EM now commonly generates close-to-atomic resolution as well as intermediate resolution maps from macromolecules observed in isolation and in situ. Interpreting these maps remains a challenging task owing to poor signal in the highest resolution shells and the necessity to select a threshold for density analysis. In order to facilitate this process, a statistical framework for the generation of confidence maps by multiple hypothesis testing and false discovery rate (FDR) control has been developed. In this way, three-dimensional confidence maps contain signal separated from background noise in the form of local detection rates of EM density values. It is demonstrated that confidence maps and FDR-based thresholding can be used for the interpretation of near-atomic resolution single-particle structures as well as lower resolution maps determined by subtomogram averaging. Confidence maps represent a conservative way of interpreting molecular structures owing to minimized noise. At the same time they provide a detection error with respect to background noise, which is associated with the density and is particularly beneficial for the interpretation of weaker cryo-EM densities in cases of conformational flexibility and lower occupancy of bound molecules and ions in the structure.
RESUMO
Compaction of DNA in chromatin is a hallmark of the eukaryotic cell and unravelling its structure is required for an understanding of DNA involving processes. Despite strong experimental efforts, many questions concerning the DNA packing are open. In particular, it is heavily debated whether an ordered structure referred to as the "30 nm fibre" exist in vivo. Scanning probe microscopy has become a cutting edge technology for the high-resolution imaging of DNA- protein complexes. Here, we perform high-resolution atomic force microscopy of non-cross-linked chromatin arrays in liquid, under different salt conditions. A statistical analysis of the data reveals that array compaction is salt dependent in a non-monotonic fashion. A simple physical model can qualitatively explain the observed findings due to the opposing effects of salt dependent stiffening of DNA, nucleosome stability and histone-histone interactions. While for different salt concentrations different compaction states are observed, our data do not provide support for the existence of regular chromatin fibres. Our studies add new insights into chromatin structure, and with that contribute to a further understanding of the DNA condensation.