BAR-based multi-dimensional nonequilibrium pulling for indirect construction of QM/MM free energy landscapes: from semi-empirical to ab initio

Seeding the multi-dimensional nonequilibrium pulling for Hamiltonian variation: indirect nonequilibrium free energy simulations at QM levels

The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method could be used to obtain the QM thermodynamics by combining the SQM or MM results and the SQM-to-QM or MM-to-QM corrections. In our previous work, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations was introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (∼20 fold) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. Introducing very short ps-length equilibrium sampling to grab more initial seeds could also be helpful. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another deficiency of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Numerical tests are performed to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The presented selection-criterion-based multi-dimensional ASMD scheme follows the same perturbation network of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.

A General Picture of Cucurbit[8]uril Host-Guest Binding

Describing, understanding, and designing complex interaction networks within macromolecular systems remain challenging in modern chemical research. Host-guest systems, despite their relative simplicity in both the structural feature and interaction patterns, still pose problems in theoretical modeling. The barrel-shaped supramolecular container cucurbit[8]uril (CB8) shows promising functionalities in various areas, e.g., catalysis and molecular recognition. It can stably coordinate a series of structurally diverse guests with high affinities. In this work, we examine the binding of seven commonly abused drugs to the CB8 host, aiming at providing a general picture of CB8-guest binding. Extensive sampling of the configurational space of these host-guest systems is performed, and the binding pathway and interaction patterns of CB8-guest complexes are investigated. A thorough comparison of widely used fixed-charge models for drug-like molecules is presented. Iterative refitting of the atomic charges suggests significant conformation dependence of charge generation. The initial model generated at the original conformation could be inaccurate for new conformations explored during conformational search, and the newly fitted charge set improves the prediction-experiment correlation significantly. Our investigations of the configurational space of CB8-drug complexes suggest that the host-guest interactions are more complex than expected. Despite the structural simplicities of these molecules, the conformational fluctuations of the host and the guest molecules and orientations of functional groups lead to the existence of an ensemble of binding modes. The insights of the binding thermodynamics, performance of fixed-charge models, and binding patterns of the CB8-guest systems are useful for studying and elucidating the binding mechanism of other host-guest complexes.

Conformational Fluctuations in GTP-Bound K-Ras: A Metadynamics Perspective with Harmonic Linear Discriminant Analysis

Biomacromolecules often undergo significant conformational rearrangements during function. In proteins, these motions typically consist in nontrivial, concerted rearrangement of multiple flexible regions. Mechanistic, thermodynamics, and kinetic predictions can be obtained via molecular dynamics simulations, provided that the simulation time is at least comparable to the relevant time scale of the process of interest. Because of the substantial computational cost, however, plain MD simulations often have difficulty in obtaining sufficient statistics for converged estimates, requiring the use of more-advanced techniques. Central in many enhanced sampling methods is the definition of a small set of relevant degrees of freedom (collective variables) that are able to describe the transitions between different metastable states of the system. The harmonic linear discriminant analysis (HLDA) has been shown to be useful for constructing low-dimensional collective variables in various complex systems. Here, we apply HLDA to study the free-energy landscape of a monomeric protein around its native state. More precisely, we study the K-Ras protein bound to GTP, focusing on two flexible loops and on the region associated with oncogenic mutations. We perform microsecond-long biased simulations on the wild type and on G12C, G12D, G12 V mutants, describe the resulting free-energy landscapes, and compare our predictions with previous experimental and computational studies. The fast interconversion between open and closed macroscopic states and their similar thermodynamic stabilities are observed. The mutation-induced effects include the alternations of the relative stabilities of different conformational states and the introduction of many microscopic metastable states. Together, our results demonstrate the applicability of the HLDA-based protocol for the conformational sampling of multiple flexible regions in folded proteins.

SAMPL7 TrimerTrip host-guest binding affinities from extensive alchemical and end-point free energy calculations

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 of the host and the guests are provided, which indicates that the predictions of both the binding poses and the binding affinities are under assessment. In this work, starting from the binding poses obtained from our previous enhanced sampling simulations in the configurational space, we perform extensive alchemical and end-point free energy calculations to calculate the host-guest binding affinities retrospectively. The alchemical predictions with two widely accepted charge schemes (i.e. AM1-BCC and RESP) are in good agreement with the experimental reference, while the end-point estimates perform poorly in reproducing the experimental binding affinities. Aside from the absolute value of the binding affinity, the rank of binding free energies is also crucial in drug design. Surprisingly, the end-point MM/PBSA method seems very powerful in reproducing the experimental rank of binding affinities. Although the length of our simulations is long and the intermediate spacing is dense, the convergence behavior is not very good, which may arise from the flexibility of the host molecule. Enhanced sampling techniques in the configurational space may be required to obtain fully converged sampling. Further, as the length of sampling in alchemical free energy calculations already achieves several hundred ns, performing direct simulations of the binding/unbinding event in the physical space could be more useful and insightful. More details about the binding pathway and mechanism could be obtained in this way. The nonequilibrium method could also be a nice choice if one insists to use the alchemical method, as the intermediate sampling is avoided to some extent.

SAMPL7 TrimerTrip host-guest binding poses and binding affinities from spherical-coordinates-biased simulations

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 sampling technique metadynamics to enhance the sampling of the binding/unbinding event, search for possible binding poses and calculate the binding affinities in all three host-guest binding cases of the 6th SAMPL challenge. In this work, we report a retrospective study on the TrimerTrip host-guest systems by employing the same protocol to investigate the TrimerTrip host in the SAMPL7 challenge. As no binding pose is provided by the SAMPL7 host, our simulations initiate from randomly selected configurations and are proceeded long enough to obtain converged free energy estimates and search for possible binding poses. The calculated binding affinities are in good agreement with the experimental reference, and the obtained binding poses serve as a nice starting point for end-point or alchemical free energy calculations. Note that as the work is performed after the close of the SAMPL7 challenge, we do not participate in the challenge and the results are not formally submitted to the SAMPL7 challenge.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.

SAMPL6 host-guest binding affinities and binding poses from spherical-coordinates-biased simulations

Host-guest binding is a challenging problem in computer simulation. The prediction of binding affinities between hosts and guests is an important part of the statistical assessment of the modeling of proteins and ligands (SAMPL) challenges. In this work, the volume-based variant of well-tempered metadynamics is employed to calculate the binding affinities of the host-guest systems in the SAMPL6 challenge. By biasing the spherical coordinates describing the relative position of the host and the guest, the initial-configuration-induced bias vanishes and all possible binding poses are explored. The agreement between the predictions and the experimental results and the observation of new binding poses indicate that the volume-based technique serves as a nice candidate for the calculation of binding free energies and the search of the binding poses.