Search

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. [ABSTRACT FROM AUTHOR]

Mass spectrometry is a powerful technique for the structural and functional characterization of biomolecules. However, it remains challenging to accurately gauge the gas-phase structure of biomolecular ions and assess to what extent native-like structures are maintained. Here we propose a synergistic approach which utilizes Förster resonance energy transfer and two types of ion mobility spectrometry (i.e., traveling wave and differential) to provide multiple constraints (i.e., shape and intramolecular distance) for structure-refinement of gas-phase ions. We add microsolvation calculations to assess the interaction sites and energies between the biomolecular ions and gaseous additives. This combined strategy is employed to distinguish conformers and understand the gas-phase structures of two isomeric α-helical peptides that might differ in helicity. Our work allows more stringent structural characterization of biologically relevant molecules (e.g., peptide drugs) and large biomolecular ions than using only a single structural methodology in the gas phase. It remains challenging to gauge the gas-phase structure of biomolecular ions and assess to what extent native-like structures are maintained. Here, the authors utilize Förster resonance energy transfer and ion mobility spectrometry for more stringent structural characterization of biomolecular ions. [ABSTRACT FROM AUTHOR]

This article documents enu, a freely-downloadable, open-source and stand-alone program written in C++ for the enumeration of the constitutional isomers and stereoisomers of a molecular formula. The program relies on graph theory to enumerate all the constitutional isomers of a given formula on the basis of their canonical adjacency matrix. The stereoisomers of a given constitutional isomer are enumerated as well, on the basis of the automorphism group of this matrix. The isomer list is then reported in the form of canonical SMILES strings within files in XML format. The specification of the molecule family of interest is very flexible and the code is optimized for computational efficiency. The algorithms and implementations underlying enu are described, and simple illustrative applications are presented. The enu code is freely available on GitHub at https://github.com/csms-ethz/CombiFF. [ABSTRACT FROM AUTHOR]

The calculation of relative binding free energies (RBFE) involves the choice of the end-state/system representation, of a sampling approach, and of a free-energy estimator. System representations are usually termed "single topology" or "dual topology". As the terminology is often used ambiguously in the literature, a systematic categorization of the system representations is proposed here. In the dual-topology approach, the molecules are simulated as separate molecules. Such an approach is relatively easy to automate for high-throughput RBFE calculations compared to the single-topology approach. Distance restraints are commonly applied to prevent the molecules from drifting apart, thereby improving the sampling efficiency. In this study, we introduce the program RestraintMaker, which relies on a greedy algorithm to find (locally) optimal distance restraints between pairs of atoms based on geometric measures. The algorithm is further extended for multi-state methods such as enveloping distribution sampling (EDS) or multi-site λ -dynamics. The performance of RestraintMaker is demonstrated for toy models and for the calculation of relative hydration free energies. The Python program can be used in script form or through an interactive GUI within PyMol. The selected distance restraints can be written out in GROMOS or GROMACS file formats. Additionally, the program provides a human-readable JSON format that can easily be parsed and processed further. The code of RestraintMaker is freely available on GitHub https://github.com/rinikerlab/restraintmaker. [ABSTRACT FROM AUTHOR]

The calculation of relative free-energy differences between different compounds plays an important role in drug design to identify potent binders for a given protein target. Most rigorous methods based on molecular dynamics simulations estimate the free-energy difference between pairs of ligands. Thus, the comparison of multiple ligands requires the construction of a "state graph", in which the compounds are connected by alchemical transformations. The computational cost can be optimized by reducing the state graph to a minimal set of transformations. However, this may require individual adaptation of the sampling strategy if a transformation process does not converge in a given simulation time. In contrast, path-free methods like replica-exchange enveloping distribution sampling (RE-EDS) allow the sampling of multiple states within a single simulation without the pre-definition of alchemical transition paths. To optimize sampling and convergence, a set of RE-EDS parameters needs to be estimated in a pre-processing step. Here, we present an automated procedure for this step that determines all required parameters, improving the robustness and ease of use of the methodology. To illustrate the performance, the relative binding free energies are calculated for a series of checkpoint kinase 1 inhibitors containing challenging transformations in ring size, opening/closing, and extension, which reflect changes observed in scaffold hopping. The simulation of such transformations with RE-EDS can be conducted with conventional force fields and, in particular, without soft bond-stretching terms. [ABSTRACT FROM AUTHOR]

The in silico prediction of partition coefficients is an important task in computer-aided drug discovery. In particular the octanol–water partition coefficient is used as surrogate for lipophilicity. Various computational approaches have been proposed, ranging from simple group-contribution techniques based on the 2D topology of a molecule to rigorous methods based molecular dynamics (MD) or quantum chemistry. In order to balance accuracy and computational cost, we recently developed the MD fingerprints (MDFPs), where the information in MD simulations is encoded in a floating-point vector, which can be used as input for machine learning (ML). The MDFP-ML approach was shown to perform similarly to rigorous methods while being substantially more efficient. Here, we present the application of MDFP-ML for the prediction of octanol–water partition coefficients in the SAMPL6 blind challenge. The underlying computational pipeline is made freely available in form of the MDFPtools package. [ABSTRACT FROM AUTHOR]

Polytheonamide B (pTB), a highly cytotoxic peptide produced by a symbiotic bacterium of the marine sponge Theonella swinhoei, forms a transmembrane pore consisting of 49 residues. More than half of its residues are posttranslationally modified. Epimerizations result in alternating L- and D-amino acids that allow the peptide to adopt a $$\beta ^{6.3}$$ -helical conformation. Unusually, the wide $$\beta ^{6.3}$$ -helix of pTB is stable in a polar environment, which is in contrast to gramicidin A, an antibiotic with similar function and structure. The role of the other posttranslational modifications (PTMs) such as side chain hydroxylations, C- and N-methylations is not well understood. In this study, the importance of these PTMs for the stability of $$\beta ^{6.3}$$ -helix is investigated using computational tools. By reverting the modified residues to their precursors and monitoring the effect on the dominant structure, we show that the N-methylations are crucial for the stability of the $$\beta ^{6.3}$$ -helix in a polar environment. They are the driving force for the formation of stable side chain hydrogen-bond chains that act as an 'exoskeleton.' Such exoskeletons could present a general design strategy for helical peptides. [ABSTRACT FROM AUTHOR]

Similarity-search methods using molecular fingerprints are an important tool for ligand-based virtual screening. A huge variety of fingerprints exist and their performance, usually assessed in retrospective benchmarking studies using data sets with known actives and known or assumed inactives, depends largely on the validation data sets used and the similarity measure used. Comparing new methods to existing ones in any systematic way is rather difficult due to the lack of standard data sets and evaluation procedures. Here, we present a standard platform for the benchmarking of 2D fingerprints. The open-source platform contains all source code, structural data for the actives and inactives used (drawn from three publicly available collections of data sets), and lists of randomly selected query molecules to be used for statistically valid comparisons of methods. This allows the exact reproduction and comparison of results for future studies. The results for 12 standard fingerprints together with two simple baseline fingerprints assessed by seven evaluation methods are shown together with the correlations between methods. High correlations were found between the 12 fingerprints and a careful statistical analysis showed that only the two baseline fingerprints were different from the others in a statistically significant way. High correlations were also found between six of the seven evaluation methods, indicating that despite their seeming differences, many of these methods are similar to each other. [ABSTRACT FROM AUTHOR]

Water molecules in the binding pocket of a protein and their role in ligand binding have increasingly raised interest in recent years. Displacement of such water molecules by ligand atoms can be either favourable or unfavourable for ligand binding depending on the change in free enthalpy. In this study, we investigate the displacement of water molecules by an apolar probe in the binding pocket of two proteins, cyclin-dependent kinase 2 and tRNA-guanine transglycosylase, using the method of enveloping distribution sampling (EDS) to obtain free enthalpy differences. In both cases, a ligand core is placed inside the respective pocket and the remaining water molecules are converted to apolar probes, both individually and in pairs. The free enthalpy difference between a water molecule and a CH group at the same location in the pocket in comparison to their presence in bulk solution calculated from EDS molecular dynamics simulations corresponds to the binding free enthalpy of CH at this location. From the free enthalpy difference and the enthalpy difference, the entropic contribution of the displacement can be obtained too. The overlay of the resulting occupancy volumes of the water molecules with crystal structures of analogous ligands shows qualitative correlation between experimentally measured inhibition constants and the calculated free enthalpy differences. Thus, such an EDS analysis of the water molecules in the binding pocket may give valuable insight for potency optimization in drug design. [ABSTRACT FROM AUTHOR]

Simulation of the dynamics of a protein in aqueous solution using an atomic model for both the protein and the many water molecules is still computationally extremely demanding considering the time scale of protein motions. The use of supra-atomic or supra-molecular coarse-grained (CG) models may enhance the computational efficiency, but inevitably at the cost of reduced accuracy. Coarse-graining solvent degrees of freedom is likely to yield a favourable balance between reduced accuracy and enhanced computational speed. Here, the use of a supra-molecular coarse-grained water model that largely preserves the thermodynamic and dielectric properties of atomic level fine-grained (FG) water in molecular dynamics simulations of an atomic model for four proteins is investigated. The results of using an FG, a CG, an implicit, or a vacuum solvent environment of the four proteins are compared, and for hen egg-white lysozyme a comparison to NMR data is made. The mixed-grained simulations do not show large differences compared to the FG atomic level simulations, apart from an increased tendency to form hydrogen bonds between long side chains, which is due to the reduced ability of the supra-molecular CG beads that represent five FG water molecules to make solvent-protein hydrogen bonds. But, the mixed-grained simulations are at least an order of magnitude faster than the atomic level ones. [ABSTRACT FROM AUTHOR]

Docking algorithms for computer-aided drug discovery and design often ignore or restrain the flexibility of the receptor, which may lead to a loss of accuracy of the relative free enthalpies of binding. In order to evaluate the contribution of receptor flexibility to relative binding free enthalpies, two host-guest systems have been examined: inclusion complexes of α-cyclodextrin (αCD) with 1-chlorobenzene (ClBn), 1-bromobenzene (BrBn) and toluene (MeBn), and complexes of DNA with the minor-groove binding ligands netropsin (Net) and distamycin (Dist). Molecular dynamics simulations and free energy calculations reveal that restraining of the flexibility of the receptor can have a significant influence on the estimated relative ligand-receptor binding affinities as well as on the predicted structures of the biomolecular complexes. The influence is particularly pronounced in the case of flexible receptors such as DNA, where a 50% contribution of DNA flexibility towards the relative ligand-DNA binding affinities is observed. The differences in the free enthalpy of binding do not arise only from the changes in ligand-DNA interactions but also from changes in ligand-solvent interactions as well as from the loss of DNA configurational entropy upon restraining. [ABSTRACT FROM AUTHOR]

New parameter sets of the GROMOS biomolecular force field, 54A7 and 54B7, are introduced. These parameter sets summarise some previously published force field modifications: The 53A6 helical propensities are corrected through new φ/ ψ torsional angle terms and a modification of the N-H, C=O repulsion, a new atom type for a charged −CH in the choline moiety is added, the Na and Cl ions are modified to reproduce the free energy of hydration, and additional improper torsional angle types for free energy calculations involving a chirality change are introduced. The new helical propensity modification is tested using the benchmark proteins hen egg-white lysozyme, fox1 RNA binding domain, chorismate mutase and the GCN4-p1 peptide. The stability of the proteins is improved in comparison with the 53A6 force field, and good agreement with a range of primary experimental data is obtained. [ABSTRACT FROM AUTHOR]