Multibasin Quasi-Harmonic Approach for the Calculation of the Configurational Entropy of Small Molecules in Solution

CREST-A program for the exploration of low-energy molecular chemical space

Conformer-rotamer sampling tool (CREST) is an open-source program for the efficient and automated exploration of molecular chemical space. Originally developed in Pracht et al. [Phys. Chem. Chem. Phys. 22, 7169 (2020)] as an automated driver for calculations at the extended tight-binding level (xTB), it offers a variety of molecular- and metadynamics simulations, geometry optimization, and molecular structure analysis capabilities. Implemented algorithms include automated procedures for conformational sampling, explicit solvation studies, the calculation of absolute molecular entropy, and the identification of molecular protonation and deprotonation sites. Calculations are set up to run concurrently, providing efficient single-node parallelization. CREST is designed to require minimal user input and comes with an implementation of the GFNn-xTB Hamiltonians and the GFN-FF force-field. Furthermore, interfaces to any quantum chemistry and force-field software can easily be created. In this article, we present recent developments in the CREST code and show a selection of applications for the most important features of the program. An important novelty is the refactored calculation backend, which provides significant speed-up for sampling of small or medium-sized drug molecules and allows for more sophisticated setups, for example, quantum mechanics/molecular mechanics and minimum energy crossing point calculations.

ChemFlow_py: a flexible toolkit for docking and rescoring

The design of accurate virtual screening tools is an open challenge in drug discovery. Several structure-based methods have been developed at different levels of approximation. Among them, molecular docking is an established technique with high efficiency, but typically low accuracy. Moreover, docking performances are known to be target-dependent, which makes the choice of the docking program and corresponding scoring function critical when approaching a new protein target. To compare the performances of different docking protocols, we developed ChemFlow_py, an automated tool to perform docking and rescoring. Using four protein systems extracted from DUD-E with 100 known active compounds and 3000 decoys per target, we compared the performances of several rescoring strategies including consensus scoring. We found that the average docking results can be improved by consensus ranking, which emphasizes the relevance of consensus scoring when little or no chemical information is available for a given target. ChemFlow_py is a free toolkit to optimize the performances of virtual high-throughput screening (vHTS). The software is publicly available at https://github.com/IFMlab/ChemFlow_py .

Molecular insights into the inhibition mechanism of harringtonine against essential proteins associated with SARS-CoV-2 entry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently posed a serious threat to global public health. Harringtonine (HT), as a small-molecule antagonist, has antiviral activity against a variety of viruses. There is evidence that HT can inhibit the SARS-CoV-2 entry into host cells by blocking the Spike protein and transmembrane protease serine 2 (TMPRSS2). However, the molecular mechanism underlying the inhibition effect of HT is largely elusive. Here, docking and all-atom molecular dynamics simulations were used to investigate the mechanism of HT against the receptor binding domain (RBD) of Spike, TMPRSS2, as well as the complex of RBD and angiotensin-converting enzyme 2 complex (RBD-ACE2). The results reveal that HT binds to all proteins primarily through hydrogen bond and hydrophobic interactions. Binding with HT influences the structural stability and dynamic motility processes of each protein. The interactions of HT with residues N33, H34 and K353 of ACE2, and residue K417 and Y453 of RBD contribute to disrupting the binding affinity between RBD and ACE2, which may hinder the virus entry into host cells. Our research provides molecular insights into the inhibition mechanism of HT against SARS-CoV-2 associated proteins, which will help for the novel antiviral drugs development.

ChemFlow─From 2D Chemical Libraries to Protein-Ligand Binding Free Energies

The accurate prediction of protein-ligand binding affinities is a fundamental problem for the rational design of new drug entities. Current computational approaches are either too expensive or inaccurate to be effectively used in virtual high-throughput screening campaigns. In addition, the most sophisticated methods, e.g., those based on configurational sampling by molecular dynamics, require significant pre- and postprocessing to provide a final ranking, which hinders straightforward applications by nonexpert users. We present a novel computational platform named ChemFlow to bridge the gap between 2D chemical libraries and estimated protein-ligand binding affinities. The software is designed to prepare a library of compounds provided in SMILES or SDF format, dock them into the protein binding site, and rescore the poses by simplified free energy calculations. Using a data set of 626 protein-ligand complexes and GPU computing, we demonstrate that ChemFlow provides relative binding free energies with an RMSE < 2 kcal/mol at a rate of 1000 ligands per day on a midsize computer cluster. The software is publicly available at https://github.com/IFMlab/ChemFlow.

Capturing the Recognition Dynamics of para-Sulfonato-calix[4]arenes by Cytochrome c: Toward a Quantitative Free Energy Assessment

Calix[n]arenes' selective recognition of protein surfaces covers a broad range of timely applications, from controlling protein assembly and crystallization to trapping partially disordered proteins. Here, the interaction of para-sulfonated calix-[4]-arenes with cytochrome c is investigated through all-atom, explicit water molecular dynamics simulations which allow characterization of two binding sites in quantitative agreement with experimental evidence. Free energy calculations based on the MM-PBSA and the attach-pull-release (APR) methods highlight key residues implicated in the recognition process and provide binding free energy results in quantitative agreement with isothermal titration calorimetry. Our study emphasizes the role of MD simulations to capture and describe the "walk" of sulfonated calix-[4]-arenes on the cytochrome c surface, with the arginine R13 as a pivotal interacting residue. Our MD investigation allows, through the quasi-harmonic multibasin (QHMB) method, probing an allosteric reinforcement of several per-residue interactions upon calixarene binding, which suggests a more complex mode of action of these supramolecular auxiliaries.

Toward Reliable and Insightful Entropy Calculations on Flexible Molecules

The absolute entropy of a flexible molecule can be approximated by the sum of a rigid-rotor-harmonic-oscillator (RRHO) entropy and a Gibbs-Shannon entropy associated to the Boltzmann distribution for the occupation of the conformational energy levels. Herein, we show that such partitioning, which has received renewed interest, leads to accurate entropies of single molecules of increasing size provided that the conformational part is estimated by means of a set of discretization and expansion techniques that are able to capture the significant correlation effects among the torsional motions. To ensure a reliable entropy estimation, we rely on extensive sampling as that produced by classical molecular dynamics simulations on the microsecond time scale, which is currently affordable for small- and medium-sized molecules. According to test calculations, the gas-phase entropy of simple organic molecules is predicted with a mean unsigned error of 0.9 cal/(mol K) when the RRHO entropies are computed at the B3LYP-D3/cc-pVTZ level. Remarkably, the same protocol gives small errors [<1 cal/(mol K)] for the extremely flexible linear alkane molecules (C_nH_2n+2, n = 14, 16, and 18). Similarly, we obtain well-converged entropies for a more challenging test of drug molecules, which exhibit more pronounced correlation effects. We also perform equivalent entropy calculations on a 76 amino acid protein, ubiquitin, by taking advantage of the cutoff-dependent formulation of an expansion technique (correlation-consistent multibody local approximation, CC-MLA), which incorporates genuine correlation effects among the neighboring dihedral angles. Moreover, we show that insightful descriptors of the coupled torsional motions can be obtained with the CC-MLA approach.

Towards understanding solvation effects on the conformational entropy of non-rigid molecules

The absolute molecular entropy is a fundamental quantity for the accurate description of thermodynamic properties. For non-rigid molecules, a substantial part of the entropy can be attributed to a conformational contribution. Systems and properties where this is relevant, e.g., protein-ligand binding affinities or pK_a values refer usually to the liquid phase. In this work, the influence of solvation on the conformational entropy is investigated. A recently introduced state-of-the-art and automated computational protocol for the computation of conformational entropies [Pracht et al., Chem. Sci., 2021, 12, 6551-6568.] is applied in combination with fast and accurate semiempirical quantum-chemical methods and implicit solvation models for a set of 25 commercially available drug molecules and five transition metal compounds. Computed gas-phase conformational entropies are compared with values obtained in implicit n-hexane and water. It is found that implicit solvation can have a substantial effect of several cal mol^-1 K^-1 on the entropy as a result of large conformational changes in the different phases. We conclude that for flexible molecules chemical accuracy for free energies in solution can only be achieved if solvation effects on the conformational ensemble are considered.

Calculation of absolute molecular entropies and heat capacities made simple

We propose a fully-automated composite scheme for the accurate and numerically stable calculation of molecular entropies by efficiently combining density-functional theory (DFT), semi-empirical methods (SQM), and force-field (FF) approximations. The scheme is systematically expandable and can be integrated seamlessly with continuum-solvation models. Anharmonic effects are included through the modified rigid-rotor-harmonic-oscillator (msRRHO) approximation and the Gibbs-Shannon formula for extensive conformer ensembles (CEs), which are generated by a metadynamics search algorithm and are extrapolated to completeness. For the first time, variations of the ro-vibrational entropy over the CE are consistently accounted-for through a Boltzmann-population average. Extensive tests of the protocol with the two standard DFT approaches B97-3c and B3LYP-D3 reveal an unprecedented accuracy with mean deviations <1 cal mol-1 K-1 (about <1-2%) for the total gas phase molecular entropy of medium-sized molecules. Even for the hardship case of extremely flexible linear alkanes (C14H30-C16H34), errors are only about 3 cal mol-1 K-1. Comprehensive tests indicate a relatively strong variation of the conformational entropy on the underlying level of theory for typical drug molecules, inferring the complex potential energy surfaces as the main source of error. Furthermore, we show some application examples for the calculation of free energy differences in typical chemical reactions.