A Scalable Molecular Force Field Parameterization Method Based on Density Functional Theory and Quantum-Level Machine Learning.

Fast and accurate molecular force field (FF) parameterization is still an unsolved problem. Accurate FF are not generally available for all molecules, like novel druglike molecules. While methods based on quantum mechanics (QM) exist to parameterize them with better accuracy, they are computationally expensive and slow, which limits applicability to a small number of molecules. Here, we present an automated FF parameterization method which can utilize either density functional theory (DFT) calculations or approximate QM energies produced by different neural network potentials (NNPs), to obtain improved parameters for molecules. We demonstrate that for the case of torchani-ANI-1x NNP, we can parameterize small molecules in a fraction of time compared with an equivalent parameterization using DFT QM calculations while producing more accurate parameters than FF (GAFF2). We expect our method to be of critical importance in computational structure-based drug discovery (SBDD). The current version is available at PlayMolecule ( www.playmolecule.org ) and implemented in HTMD, allowing to parameterize molecules with different QM and NNP options.

Molecular-Simulation-Driven Fragment Screening for the Discovery of New CXCL12 Inhibitors.

Fragment-based drug discovery (FBDD) has become a mainstream approach in drug design because it allows the reduction of the chemical space and screening libraries while identifying fragments with high protein-ligand efficiency interactions that can later be grown into drug-like leads. In this work, we leverage high-throughput molecular dynamics (MD) simulations to screen a library of 129 fragments for a total of 5.85 ms against the CXCL12 monomer, a chemokine involved in inflammation and diseases such as cancer. Our in silico binding assay was able to recover binding poses, affinities, and kinetics for the selected library and was able to predict 8 mM-affinity fragments with ligand efficiencies higher than 0.3. All of the fragment hits present a similar chemical structure, with a hydrophobic core and a positively charged group, and bind to either sY7 or H1S68 pockets, where they share pharmacophoric properties with experimentally resolved natural binders. This work presents a large-scale screening assay using an exclusive combination of thousands of short MD adaptive simulations analyzed with a Markov state model (MSM) framework.

Drug Discovery and Molecular Dynamics: Methods, Applications and Perspective Beyond the Second Timescale.

Bio-molecular dynamics (MD) simulations based on graphical processing units (GPUs) were first released to the public in the early 2009 with the code ACEMD. Almost 8 years after, applications now encompass a broad range of molecular studies, while throughput improvements have opened the way to millisecond sampling timescales. Based on an extrapolation of the amount of sampling in published literature, the second timescale will be reached by the year 2022, and therefore we predict that molecular dynamics is going to become one of the main tools in drug discovery in both academia and industry. Here, we review successful applications in the drug discovery domain developed over these recent years of GPU-based MD. We also retrospectively analyse limitations that have been overcome over the years and give a perspective on challenges that remain to be addressed.

SPECTRa-T: machine-based data extraction and semantic searching of chemistry e-theses.

The SPECTRa-T project has developed text-mining tools to extract named chemical entities (NCEs), such as chemical names and terms, and chemical objects (COs), e.g., experimental spectral assignments and physical chemistry properties, from electronic theses (e-theses). Although NCEs were readily identified within the two major document formats studied, only the use of structured documents enabled identification of chemical objects and their association with the relevant chemical entity (e.g., systematic chemical name). A corpus of theses was analyzed and it is shown that a high degree of semantic information can be extracted from structured documents. This integrated information has been deposited in a persistent Resource Description Framework (RDF) triple-store that allows users to conduct semantic searches. The strength and weaknesses of several document formats are reviewed.

SPECTRa: the deposition and validation of primary chemistry research data in digital repositories.

The SPECTRa (Submission, Preservation and Exposure of Chemistry Teaching and Research Data) project has investigated the practices of chemists in archiving and disseminating primary chemical data from academic research laboratories. To redress the loss of the large amount of data never archived or disseminated, we have developed software for data publication into departmental and institutional Open Access digital repositories (DSpace). Data adhering to standard formats in selected disciplines (crystallography, NMR, computational chemistry) is transformed to XML (CML, Chemical Markup Language) which provides added validation. Context-specific chemical metadata and persistent Handle identifiers are added to enable long-term data reuse. It was found essential to provide an embargo mechanism, and policies for operating this and other processes are presented.