A major part of chemical conversions is carried out in the fluid phase, where an accurate modeling of the involved reactions requires to also take into account solvation effects. Implicit solvation models often cover these effects with sufficient accuracy but can fail drastically when specific solvent–solute interactions are important. In those cases, microsolvation, i.e., the explicit inclusion of

Quantitative structure–activity relationship (QSAR) and quantitative structure–property relationship (QSPR) models predict biological activity and molecular property based on the numerical relationship between chemical structures and activity (property) values. Molecular representations are of importance in QSAR/QSPR analysis. Topological information of molecular structures is usually utilized (2D

The parallel artificial membrane permeability assay (PAMPA), a non-cellular lab-based assay, is extensively used to measure the permeability of pharmaceutical compounds. PAMPA experiments provide a working mimic of a molecule passing through cells and PAMPA values are widely used to estimate drug absorption parameters. There is an increased interest in developing computational methods to predict PAMPA

The prediction of acid dissociation constants (pKa) is a prerequisite for predicting many other properties of a small molecule, such as its protein–ligand binding affinity, distribution coefficient (log D), membrane permeability, and solubility. The prediction of each of these properties requires knowledge of the relevant protonation states and solution free energy penalties of each state. The SAMPL6

In the context of the SAMPL7 challenge, we computed, employing a non-equilibrium (NE) alchemical technique, the standard binding free energy of two series of host-guest systems, involving as a host the Isaac’s TrimerTrip, a Cucurbituril-like open cavitand, and the Gilson’s Cyclodextrin derivatives. The adopted NE alchemy combines enhanced sampling molecular dynamics simulations with driven fast out-of-equilibrium

Redox-sensitive cysteine (RSC) thiol contributes to many biological processes. The identification of RSC plays an important role in clarifying some mechanisms of redox-sensitive factors; nonetheless, experimental investigation of RSCs is expensive and time-consuming. The computational approaches that quickly and accurately identify candidate RSCs using the sequence information are urgently needed.

The SAMPL challenges focus on testing and driving progress of computational methods to help guide pharmaceutical drug discovery. However, assessment of methods for predicting binding affinities is often hampered by computational challenges such as conformational sampling, protonation state uncertainties, variation in test sets selected, and even lack of high quality experimental data. SAMPL blind challenges

The alternative oxidase (AOX) is a monotopic diiron carboxylate protein that catalyses the oxidation of ubiquinol and the reduction of oxygen to water. Although a number of AOX inhibitors have been discovered, little is still known about the ligand–protein interaction and essential chemical characteristics of compounds required for a potent inhibition. Furthermore, owing to the rapidly growing resistance

In the current work we report on our participation in the SAMPL7 challenge calculating absolute free energies of the host–guest systems, where 2 guest molecules were probed against 9 hosts-cyclodextrin and its derivatives. Our submission was based on the non-equilibrium free energy calculation protocol utilizing an averaged consensus result from two force fields (GAFF and CGenFF). The submitted prediction

Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) challenges provide routes to compare chemical quantities determined using computational chemistry approaches to experimental measurements that are shared after the competition. For this effort, several computational methods have been used to calculate the binding energies of Octa Acid (OA) and exo-Octa Acid (exoOA) host–guest systems

As part of the SAMPL7 host–guest binding challenge, the AMOEBA force field was applied to calculate the absolute binding free energy for 16 charged organic ammonium guests to the TrimerTrip host, a recently reported acyclic cucurbituril-derived clip host structure with triptycene moieties at its termini. Here we report binding free energy calculations for this system using the AMOEBA polarizable atomic

Conventional protein:ligand crystallographic refinement uses stereochemistry restraints coupled with a rudimentary energy functional to ensure the correct geometry of the model of the macromolecule—along with any bound ligand(s)—within the context of the experimental, X-ray density. These methods generally lack explicit terms for electrostatics, polarization, dispersion, hydrogen bonds, and other key

The recent outbreak of the respiratory syndrome-related coronavirus (SARS-CoV-2) is stimulating an unprecedented scientific campaign to alleviate the burden of the coronavirus disease (COVID-19). One line of research has focused on targeting SARS-CoV-2 proteins fundamental for its replication by repurposing drugs approved for other diseases. The first interaction between the virus and the host cell

Selective progesterone receptor modulators are promising therapeutic options for the treatment of uterine fibroids. Vilaprisan, a new chemical entity that was discovered at Bayer is currently in clinical development. In this study we provide a combined experimental and quantum chemical approach providing the data that allowed to present hydroxyestradienone as an acceptable starting material for drug

Joint academic–industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein–ligand interaction up to statistical analyses of biological data. This way, method

Anaplastic lymphoma kinase (ALK) has been thought to be a prospective target of anti-drug resistance design in treatment of tumors and specific neuron diseases. It is highly useful for the seeking of possible strategy alleviating drug resistance to probe the mutation-mediated effect on binding of inhibitors to ALK. In the current work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD)

Fatty acid binding proteins (FABPs) are small intracellular proteins that reversibly bind fatty acids and other hydrophobic ligands. In cestodes, due to their inability to synthesise fatty acids and cholesterol de novo, FABPs, together with other lipid binding proteins, have been proposed as essential, involved in the trafficking and delivery of such lipophilic metabolites. Pharmacological agents that

We investigate the binding of native β-cyclodextrin (β-CD) and eight novel β-CD derivatives with two different guest compounds, using isothermal calorimetry and 2D NOESY NMR. In all cases, the stoichiometry is 1:1 and binding is exothermic. Overall, modifications at the 3′ position of β-CD, which is at the secondary face, weaken binding by several kJ/mol relative to native β-CD, while modifications

The prediction of host–guest binding affinities with computational modelling is still a challenging task. In the 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge, a new host named TrimerTrip was synthesized and the thermodynamic parameters of 16 structurally diverse guests binding to the host were characterized. In the TrimerTrip-guest challenge, only structures

Computational protein–ligand docking is well-known to be prone to inaccuracies in input receptor structures, and it is challenging to obtain good docking results with computationally predicted receptor structures (e.g. through homology modeling). Here we introduce a fragment-based docking method and test if it reduces requirements on the accuracy of an input receptor structures relative to non-fragment

Atomistic simulations have become an invaluable tool for industrial applications ranging from the optimization of protein-ligand interactions for drug discovery to the design of new materials for energy applications. Here we review recent advances in the use of machine learning (ML) methods for accelerated simulations based on a quantum mechanical (QM) description of the system. We show how recent

The compound optimization monitor (COMO) approach was originally developed as a diagnostic approach to aid in evaluating development stages of analog series and progress made during lead optimization. COMO uses virtual analog populations for the assessment of chemical saturation of analog series and has been further developed to bridge between optimization diagnostics and compound design. Herein, we

Drug discovery is an expensive and time-consuming process. To make this process more efficient quantum chemistry methods can be employed. The electrophilicity index is one property that can be calculated by quantum chemistry methods, and if calculated correctly gives insight into the reactivity of covalent inhibitors. Herein we present the usage of the electrophilicity index on three common warheads

Human cathelicidin LL-37 has recently attracted interest as a potential therapeutic agent, mostly because of its ability to kill a wide variety of pathogens and cancer cells. In this study, we used molecular dynamics simulation aimed to get insights that help to correlate with the antibacterial activity of previously designed LL-37 anticancer derivative (i.e. GF-17). Two independent molecular dynamics

Water molecules can be found interacting with the surface and within cavities in proteins. However, water exchange between bulk and buried hydration sites can be slow compared to simulation timescales, thus leading to the inefficient sampling of the locations of water. This can pose problems for free energy calculations for computer-aided drug design. Here, we apply a hybrid method that combines nonequilibrium

A proinflammatory peptide (PIP) is a type of signaling molecules that are secreted from immune cells, which contributes to the first line of defense against invading pathogens. Numerous experiments have shown that PIPs play an important role in human physiology such as vaccines and immunotherapeutic drugs. Considering high-throughput laboratory methods that are time consuming and costly, effective

SARS-CoV-2 recently jumped species and rapidly spread via human-to-human transmission to cause a global outbreak of COVID-19. The lack of effective vaccine combined with the severity of the disease necessitates attempts to develop small molecule drugs to combat the virus. COVID19_GIST_HSA is a freely available online repository to provide solvation thermodynamic maps of COVID-19-related protein small

We propose a method to identify the correct binding mode of a ligand with a protein among multiple predicted docking poses. Our method consists of two steps. First, five independent MD simulations with different initial velocities are performed for each docking pose, in order to evaluate its stability. If the root-mean-square deviations (RMSDs) of heavy atoms from the docking pose are larger than a

The free fatty acid receptor 1 (FFAR1, formerly GPR40), is a potential G protein-coupled receptor (GPCR) target for the treatment of type 2 diabetes mellitus (T2DM), as it enhances glucose-dependent insulin secretion upon activation by endogenous long-chain free fatty acids. The presence of two allosterically communicating binding sites and the lack of the conserved GPCR structural motifs challenge

Recent breakthroughs in G protein-coupled receptor (GPCR) crystallography and the subsequent increase in number of solved GPCR structures has allowed for the unprecedented opportunity to utilize their experimental structures for structure-based drug discovery applications. As virtual screening represents one of the primary computational methods used for the discovery of novel leads, the GPCR-Bench

There is a pressing need to improve the efficiency of drug development, and nowhere is that need more clear than in the case of neglected diseases like malaria. The peculiarities of pyrimidine metabolism in Plasmodium species make inhibition of dihydroorotate dehydrogenase (DHODH) an attractive target for antimalarial drug design. By applying a pair of complementary quantitative structure–activity

In contrast to the computational generation of conventional tautomers, the analogous operation that would produce ring-chain tautomers is rarely available in cheminformatics codes. This is partly due to the perceived unimportance of ring-chain tautomerism and partly because specialized algorithms are required to realize the non-local proton transfers that occur during ring-chain rearrangement. Nevertheless

Previously a fingerprint based on monomer composition (MCFP) of nonribosomal peptides (NRPs) has been introduced. MCFP is a novel method for obtaining a representative description of NRP structures from their monomer composition in a fingerprint form. An effective screening and prediction of biological activities has been obtained from Norine NRPs database. In this paper, we present an extension of

Host–guest binding remains a major challenge in modern computational modelling. The newest 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge contains a new series of host–guest systems. The TrimerTrip host binds to 16 structurally diverse guests. Previously, we have successfully employed the spherical coordinates as the collective variables coupled with the enhanced

G protein-coupled receptors (GPCR) comprise the largest family of membrane proteins and are of considerable interest as targets for drug development. However, many GPCR structures remain unsolved. To address the structural ambiguity of these receptors, computational tools such as homology modeling and loop modeling are often employed to generate predictive receptor structures. Here we combined both

The exchange proteins activated by cAMP (EPAC) are implicated in a large variety of physiological processes and they are considered as promising targets for a wide range of therapeutic applications. Several recent reports provided evidence for the therapeutic effectiveness of the inhibiting EPAC1 activity cardiac diseases. In that context, we recently characterized a selective EPAC1 antagonist named

Computer-aided strategies are useful for reducing the costs and increasing the success-rate in drug discovery. Among these strategies, methods based on pharmacophores (an ensemble of electronic and steric features representing the target active site) are efficient to implement over large compound libraries. However, traditional pharmacophore-based methods require knowledge of active compounds or ligand–receptor

Experimental 3D structures of calcium channels with phenylalkylamines (PAAs) provide basis for further analysis of atomic mechanisms of these important cardiovascular drugs. In the crystal structure of the engineered calcium channel CavAb with Br-verapamil and in the cryo-EM structure of the Cav1.1 channel with verapamil, the ligands bind in the inner pore. However, there are significant differences

In order to assess safety and efficacy of small molecule drugs as well as agrochemicals, it is key to understanding the nature of protein–ligand interaction on an atomistic level. Prothioconazole (PTZ), although commonly considered to be an azole-like inhibitor of sterol 14-α demethylase (CYP51), differs from classical azoles with respect to how it binds its target. The available evidence is only indirect

Nowadays, the importance of computational methods in the design of therapeutic agents in a more efficient way is indisputable. Particularly, these methods have been important in the design of novel acetylcholinesterase enzyme inhibitors related to Alzheimer’s disease. In this sense, in this report a computational model of linear prediction of acetylcholinesterase inhibitory activity of steroids and

We investigate the minimum-energy path for the rotation of formal C=N double bonds in molecules with guanidine-like substructures as present in the chemical class of neonicotinoids. The transitions between the E- and Z-isomers of several neonicotinoids using scans of the torsional potential energy hypersurfaces are quantified at the DFT-level of theory. The validity of using this ansatz is checked

Herein, the LASSBio Chemical Library is presented as a valuable source of compounds for screening to identify hits suitable for subsequent hit-to-lead optimization stages. A feature of the LASSBio Chemical Library worth highlighting is the fact that it is a smart library designed by medicinal chemists with pharmacological activity as the main priority. The great majority of the compounds part of this

Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H1 antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation

Phage virion protein (PVP) perforate the host cell membrane and eventually culminates in cell rupture thereby releasing replicated phages. The accurate identification of PVP is thus a crucial step towards improving our understanding of the biological function and mechanisms of PVPs. Therefore, it is desirable to develop a computational method that is capable of fast and accurate identification of PVPs

This study is directed toward assessing the predictive potential of eigenvalue-based topological molecular descriptors. The graph energy, Estrada index, resolvent energy, and the Laplacian energy were tested as parameters for the prediction of boiling points, heats of formation, and octanol/water partition coefficients of alkanes. It was shown that an eigenvalue-based molecular descriptor cannot be

We present a study based on density functional theory calculations to explore the rate limiting steps of product formation for oxidation by Flavin-containing Monooxygenase (FMO) and glucuronidation by the UDP-glucuronosyltransferase (UGT) family of enzymes. FMOs are responsible for the modification phase of metabolism of a wide diversity of drugs, working in conjunction with Cytochrome P450 (CYP) family

Classification of chemical compounds of plants as a source of medicaments for neurodegenerative diseases through computer screening is an efficient process in drug discovery, in advance of laboratory testing and clinical trials. The onset of neurodegenerative disorders incarcerates both sufferers and their families mentally and financially. This investigation emphasises the search for potent compounds

Optimization in medicinal chemistry often involves designing replacements for a section of a molecule which aim to retain potency while improving other properties of the compound. In this study, we perform a retrospective analysis using a number of computational methods to identify active side chains amongst a pool of random decoy side chains, mimicking a similar procedure that might be undertaken

Activity cliffs (ACs) consist of structurally similar compounds with a large difference in potency against their target. Accordingly, ACs introduce discontinuity in structure-activity relationships (SARs) and are a prime source of SAR information. In compound data sets, the vast majority of ACs are formed by differently sized groups of structurally similar compounds with large potency variations. As

The cholinesterases are essential targets implicated in the pathogenesis of Alzheimer’s disease (AD). In the present study, virtual screening and molecular docking are performed to identify the potential hits. Docking-post processing (DPP) and pose filtration protocols against AChE and BChE resulted in three hits (AW00308, HTS04089, and JFD03947). Molecular Mechanics-Generalized Born Surface Area (MM-GBSA)

CREB-binding protein (CBP) is a multi-subunit scaffold protein complex in transcription regulation process, binding and interacting with ligands such as mixed-lineage leukemia (MLL) and c-Myb allosterically. Here in this study, we have revisited the concept of allostery in CBP via residue-based interaction energy calculation based on molecular dynamics (MD) simulations. To this end, we conducted MD

The activity cliff (AC) concept is of comparable relevance for medicinal chemistry and chemoinformatics. An AC is defined as a pair of structurally similar compounds with a large potency difference against a given target. In medicinal chemistry, ACs are of interest because they reveal small chemical changes with large potency effects, a concept referred to as structure-activity relationship (SAR) discontinuity

Difficulties in interpreting machine learning (ML) models and their predictions limit the practical applicability of and confidence in ML in pharmaceutical research. There is a need for agnostic approaches aiding in the interpretation of ML models regardless of their complexity that is also applicable to deep neural network (DNN) architectures and model ensembles. To these ends, the SHapley Additive

Theoretical approaches for predicting physicochemical properties are valuable tools for accelerating the drug discovery process. In this work, quantum chemical methods are used to predict water-octanol partition coefficients as a part of the SAMPL6 blind challenge. The SMD continuum solvent model was employed with MP2 and eight DFT functionals in conjunction with correlation consistent basis sets to

Blind predictions of octanol/water partition coefficients at 298 K for 11 kinase inhibitor fragment like compounds were made for the SAMPL6 challenge. We used the conventional, "untrained", free energy based approach wherein the octanol/water partition coefficient was computed directly as the difference in solvation free energy in water and 1-octanol. We additionally proposed and used two different

All-atom molecular dynamics simulations with stratified alchemical free energy calculations were used to predict the octanol-water partition coefficient [Formula: see text] of eleven small molecules as part of the SAMPL6-[Formula: see text] blind prediction challenge using four different force field parametrizations: standard OPLS-AA with transferable charges, OPLS-AA with non-transferable CM1A charges

Two different types of approaches: (a) approaches that combine quantitative structure activity relationships, quantum mechanical electronic structure methods, and machine-learning and, (b) electronic structure vertical solvation approaches, were used to predict the logP coefficients of 11 molecules as part of the SAMPL6 logP blind prediction challenge. Using electronic structures optimized with density