Journal of Chemical Information and Modeling 最新文献

Machine learning (ML) models have become key in decision-making for many disciplines, including drug discovery and medicinal chemistry. ML models are generally evaluated prior to their usage in high-stakes decisions, such as compound synthesis or experimental testing. However, no ML model is robust or predictive in all real-world scenarios. Therefore, uncertainty quantification (UQ) in ML predictions has gained importance in recent years. Many investigations have focused on developing methodologies that provide accurate uncertainty estimates for ML-based predictions. Unfortunately, there is no UQ strategy that consistently provides robust estimates about model's applicability on new samples. Depending on the dataset, prediction task, and algorithm, accurate uncertainty estimations might be unfeasible to obtain. Moreover, the optimum UQ metric also varies across applications, and previous investigations have shown a lack of consistency across benchmarks. Herein, the UNIQUE (UNcertaInty QUantification bEnchmarking) framework is introduced to facilitate a comparison of UQ strategies in ML-based predictions. This Python library unifies the benchmarking of multiple UQ metrics, including the calculation of nonstandard UQ metrics (combining information from the dataset and model), and provides a comprehensive evaluation. In this framework, UQ metrics are evaluated for different application scenarios, e.g., eliminating the predictions with the lowest confidence or obtaining a reliable uncertainty estimate for an acquisition function. Taken together, this library will help to standardize UQ investigations and evaluate new methodologies.

In the realm of biomedical research, understanding the intricate structure of proteins is crucial, as these structures determine how proteins function within our bodies and interact with potential drugs. Traditionally, methods like X-ray crystallography and cryo-electron microscopy have been used to unravel these structures, but they are often challenging, time-consuming and costly. Recently, a breakthrough in computational biology has emerged with the development of deep learning algorithms capable of predicting protein structures based on their amino acid sequences (Jumper, J., et al. Nature 2021, 596, 583. Lane, T. J. Nature Methods 2023, 20, 170. Kryshtafovych, A., et al. Proteins: Structure, Function and Bioinformatics 2021, 89, 1607). This study focuses on predicting the dynamic changes that proteins undergo upon ligand binding, specifically when they bind to allosteric sites, i.e. a pocket different from the active site. Allosteric modulators are particularly important for drug discovery, as they open new avenues for designing drugs that can target proteins more effectively and with fewer side effects (Nussinov, R.; Tsai, C. J. Cell 2013, 153, 293). To study this, we curated a data set of 578 X-ray structures comprised of proteins displaying orthosteric and allosteric binding as well as a general framework to evaluate deep learning-based structure prediction methods. Our findings demonstrate the potential and current limitations of deep learning methods, such as AlphaFold2 (Jumper, J., et al. Nature 2021, 596, 583), NeuralPLexer (Qiao, Z., et al. Nat Mach Intell 2024, 6, 195), and RoseTTAFold All-Atom (Krishna, R., et al. Science 2024, 384, eadl2528) to predict not just static protein structures but also the dynamic conformational changes. Herein we show that predicting the allosteric induce-fit conformation still poses a challenge to deep learning methods as they more accurately predict the orthosteric bound conformation compared to the allosteric induce fit conformation. For AlphaFold2, we observed that conformational diversity, and sampling between the apo and holo state could be increased by modifying the MSA depth, but this did not enhance the ability to generate conformations close to the allosteric induced-fit conformation. To further support advancements in protein structure prediction field, the curated data set and evaluation framework are made publicly available.

Alchemical free energy campaigns can be planned using graph theory by building networks that contain nodes representing molecules that are connected by possible transformations as edges. We introduce Konnektor, an open-source Python package, for systematically planning, modifying, and analyzing free energy calculation networks. Konnektor is designed to aid in the drug discovery process by enabling users to easily setup free energy campaigns using complex graph manipulation methods. The package contains functions for network operations including concatenation of networks, deletion of transformations, and clustering of molecules along with a framework for combining these tools with existing network generation algorithms to enable the development of more complex methods for network generation. A comparison of the various network layout features offered is carried out using toy data sets. Additionally, Konnektor contains visualization and analysis tools, making the investigation of network features much simpler. Besides the content of the package, the paper also offers application examples, demonstrating how Konnektor can be used and how the different networks perform from a graph theory perspective. Konnektor is freely available via GitHub at https://github.com/OpenFreeEnergy/konnektor under the permissive MIT License.

Coronavirus disease 2019 (COVID-19) is caused by a new, highly pathogenic severe-acute-respiratory syndrome coronavirus 2 (SARS-CoV-2) that infects human cells through its transmembrane spike (S) glycoprotein. The receptor-binding domain (RBD) of the S protein interacts with the angiotensin-converting enzyme II (ACE2) receptor of the host cells. Therefore, pharmacological targeting of this interaction might prevent infection or spread of the virus. Here, we performed a virtual screening to identify small molecules that block S-ACE2 interaction. Large compound libraries were filtered for drug-like properties, promiscuity and protein-protein interaction-targeting ability based on their ADME-Tox descriptors and also to exclude pan-assay interfering compounds. A properly designed AI-based virtual screening pipeline was applied to the remaining compounds, comprising approximately 10% of the starting data sets, searching for molecules that could bind to the RBD of the S protein. All molecules were sorted according to their screening score, grouped based on their structure and postfiltered for possible interaction patterns with the ACE2 receptor, yielding 31 hits. These hit molecules were further tested for their inhibitory effect on Spike RBD/ACE2 (19-615) interaction. Six compounds inhibited the S-ACE2 interaction in a dose-dependent manner while two of them also prevented infection of human cells from a pseudotyped virus whose entry is mediated by the S protein of SARS-CoV-2. Of the two compounds, the benzimidazole derivative CKP-22 protected Vero E6 cells from infection with SARS-CoV-2, as well. Subsequent, hit-to-lead optimization of CKP-22 was effected through the synthesis of 29 new derivatives of which compound CKP-25 suppressed the Spike RBD/ACE2 (19-615) interaction, reduced the cytopathic effect of SARS-CoV-2 in Vero E6 cells (IC₅₀ = 3.5 μM) and reduced the viral load in cell culture supernatants. Early in vitro ADME-Tox studies showed that CKP-25 does not possess biodegradation or liver metabolism issues, while isozyme-specific CYP450 experiments revealed that CKP-25 was a weak inhibitor of the CYP450 system. Moreover, CKP-25 does not elicit mutagenic effect on Escherichia coli WP2 uvrA strain. Thus, CKP-25 is considered a lead compound against COVID-19 infection.

Natural product (NP) databases are crucial tools in computer-aided drug design (CADD). Over the past decade, there has been a worldwide effort to assemble information regarding natural products (NPs) isolated and characterized in certain geographical regions. In 2023, it was published LANaPDB, and to our knowledge, this is the first attempt to gather and standardize all the NP databases of Latin America. Herein, we present and analyze in detail the contents of an updated version of LANaPDB, which includes 619 newly added compounds from Colombia, Costa Rica, and Mexico. The present version of LANaPDB has a total of 13 578 compounds, coming from ten databases of seven Latin American countries. A chemoinformatic characterization of LANaPDB was carried out, which includes the structural classification of the compounds, calculation of six physicochemical properties of pharmaceutical interest, and visualization of the chemical space by employing and comparing two different fingerprints (MACCS keys (166-bit) and Morgan2 (2048-bit)). Furthermore, additional analyses were made, and valuable information not included in the first version of LANaPDB was added, which includes structural diversity, molecular complexity, synthetic feasibility, commercial availability, and reported and predicted biological activity. In addition, the LANaPDB compounds were cross-referenced to two of the largest public chemical compound databases annotated with biological activity: ChEMBL and PubChem.

MicroRNAs (miRs) are short, noncoding RNA strands that regulate the activity of mRNAs by affecting the repression of protein translation, and their dysregulation has been implicated in several pathologies. miR21 in particular has been implicated in tumorigenesis and anticancer drug resistance, making it a critical target for drug design. miR21 biogenesis involves precise biochemical pathways, including the cleavage of its precursor, pre-miR21, by the enzyme Dicer. The present work investigates the conformational dynamics of pre-miR21, focusing on the role of adenine29 in switching between Dicer-binding-prone and inactive states. We also investigated the effect of L50, a cyclic peptide binder of pre-miR21 and a weak inhibitor of its processing. Using time series data and our novel collective variable-based enhanced sampling technique, OneOPES, we simulated these conformational changes and assessed the effect of L50 on the conformational plasticity of pre-miR21. Our results provide insight into peptide-induced conformational changes and pave the way for the development of a computational platform for the screening of inhibitors of pre-miR21 processing that considers RNA flexibility, a stepping stone for effective structure-based drug design, with potentially broad applications in drug discovery.

Emerging evidence suggests that physiological C-terminal truncation of α-synuclein (αS) plays a critical role in regulating liquid-liquid phase separation and promoting amyloid aggregation, processes implicated in neurodegenerative diseases such as Parkinson's disease (PD). However, the molecular mechanisms through which C-terminal truncation influences αS conformation and modulates its aggregation remain poorly understood. In this study, we investigated the impact of C-terminal truncation on αS conformational dynamics by comparing full-length αS_1-140 with truncated αS_1-103 monomers using atomistic discrete molecular dynamics simulations. Our findings revealed that both αS_1-140 and αS_1-103 primarily adopted helical conformations around residues 7-32, while residues 36-95, located in the second half of the N-terminal and NAC domains, predominantly formed a dynamic β-sheet core. The C-terminus of αS_1-140 was largely unstructured and dynamically wrapped around the β-sheet core. While residues 1-95 exhibited similar secondary structure propensities in both αS_1-140 and αS_1-103, the dynamic capping by the C-terminus in αS_1-140 slightly enhanced β-sheet formation around residues 36-95. In contrast, key aggregation-driving regions (residues 2-9, 36-42, 45-57, and 68-78) were dynamically shielded by the C-terminus in αS_1-140, reducing their exposure and potentially preventing interpeptide interactions that drive aggregation. C-terminal truncation, on the other hand, increased the exposed surface area of these aggregation-prone regions, thereby enhancing interpeptide interactions, phase separation, and amyloid aggregation. Overall, our simulations provide valuable insights into the conformational effects of C-terminal truncation on αS and its role in promoting pathological aggregation.

Graphs are one of the most natural and powerful representations available for molecules; natural because they have an intuitive correspondence to skeletal formulas, the language used by chemists worldwide, and powerful, because they are highly expressive both globally (molecular topology) and locally (atom and bond properties). Graph kernels are used to transform molecular graphs into fixed-length vectors, which, based on their capacity of measuring similarity, can be used as fingerprints for machine learning (ML). To date, graph kernels have mostly focused on the atomic nodes of the graph. In this work, we developed a graph kernel based on atom-atom, bond-bond, and bond-atom (AABBA) autocorrelations. The resulting vector representations were tested on regression ML tasks on a data set of transition metal complexes; a benchmark motivated by the higher complexity of these compounds relative to organic molecules. In particular, we tested different flavors of the AABBA kernel in the prediction of the energy barriers and bond distances of the Vaska’s complex data set (Friederich et al., Chem. Sci., 2020, 11, 4584). For a variety of ML models, including neural networks, gradient boosting machines, and Gaussian processes, we showed that AABBA outperforms the baseline including only atom-atom autocorrelations. Dimensionality reduction studies also showed that the bond-bond and bond-atom autocorrelations yield many of the most relevant features. We believe that the AABBA graph kernel can accelerate the exploration of large chemical spaces and inspire novel molecular representations in which both atomic and bond properties play an important role.

To propose new mechanism-based therapeutics for Alzheimer’s disease (AD), it is crucial to study the kinetics and oligomerization/aggregation mechanisms of the hallmark tau proteins, which have various isoforms and are intrinsically disordered. In this study, multiple all-atom (AA) and coarse-grained (CG) force fields (FFs) have been benchmarked on molecular dynamics (MD) simulations of K18 tau (M243–E372), which is a truncated form (130 residues) of full-length tau (441 residues). FF19SB is first excluded because the dynamics are too slow, and the conformations are too stable. All other benchmarked AAFFs (Charmm36m, FF14SB, Gromos54A7, and OPLS-AA) and CGFFs (Martini3 and Sirah2.0) exhibit a trend of shrinking K18 tau into compact structures with the radius of gyration (ROG) around 2.0 nm, which is much smaller than the experimental value of 3.8 nm, within 200 ns of AA-MD or 2000 ns of CG-MD. Gromos54A7, OPLS-AA, and Martini3 shrink much faster than the other FFs. To perform meaningful postanalysis of various properties, we propose a strategy of selecting snapshots with 2.5 < ROG < 4.5 nm, instead of using all sampled snapshots. The calculated chemical shifts of all C, CA, and CB atoms have very good and close root-mean-square error (RMSE) values, while Charmm36m and Sirah2.0 exhibit better chemical shifts of N than other FFs. Comparing the calculated distributions of the distance between the CA atoms of CYS291 and CYS322 with the results of the FRET experiment demonstrates that Charmm36m is a perfect match with the experiment while other FFs exhibit limitations. In summary, Charmm36m is recommended as the best AAFF, and Sirah2.0 is recommended as an excellent CGFF for simulating tau K18.

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. To accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO₃)₃ dissolved in acetonitrile with a terpyridine ligand (terpyNO₂), and Nd(NO₃)₃ dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO₃)₃(terpyNO₂) and Nd(NO₃)(acetonitrile)₃, with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.