RESUMEN
We employ a combination of accelerated molecular dynamics and machine learning to unravel how the dynamic characteristics of CBL-B and C-CBL confer their binding affinity and selectivity for ligands from subtle structural disparities within their binding pockets and dissociation pathways. Our predictive model of dissociation rate constants (k off) demonstrates a moderate correlation between predicted k off and experimental IC50 values, which is consistent with experimental k off and τ-random accelerated molecular dynamics (τRAMD) results. By employing a linear regression of dissociation trajectories, we identified key amino acids in binding pockets and along the dissociation paths responsible for activity and selectivity. These amino acids are statistically significant in achieving activity and selectivity and contribute to the primary structural discrepancies between CBL-B and C-CBL. Moreover, the binding free energies calculated from molecular mechanics with generalized Born and surface area solvation (MM/GBSA) highlight the ΔG difference between CBL-B and C-CBL. The k off prediction, together with the key amino acids, provides important guides for designing drugs with high selectivity.
RESUMEN
Exclusively prioritizing the precision of energy prediction frequently proves inadequate in satisfying multifaceted requirements. A heightened focus is warranted on assessing the rationality of potential energy curves predicted by machine learning-based force fields (MLFFs), alongside evaluating the pragmatic utility of these MLFFs. This study introduces SWANI, an optimized neural network potential stemming from the ANI framework. Through the incorporation of supplementary physical constraints, SWANI aligns more cohesively with chemical expectations, yielding rational potential energy profiles. It also exhibits superior predictive precision compared with that of the ANI model. Additionally, a comprehensive comparison is conducted between SWANI and a prominent graph neural network-based model. The findings indicate that SWANI outperforms the latter, particularly for molecules exceeding the dimensions of the training set. This outcome underscores SWANI's exceptional capacity for generalization and its proficiency in handling larger molecular systems.
RESUMEN
Drug design based on kinetic properties is growing in application. Here, we applied retrosynthesis-based pre-trained molecular representation (RPM) in machine learning (ML) to train 501 inhibitors of 55 proteins and successfully predicted the dissociation rate constant (koff) values of 38 inhibitors from an independent dataset for the N-terminal domain of heat shock protein 90α (N-HSP90). Our RPM molecular representation outperforms other pre-trained molecular representations such as GEM, MPG, and general molecular descriptors from RDKit. Furthermore, we optimized the accelerated molecular dynamics to calculate the relative retention time (RT) for the 128 inhibitors of N-HSP90 and obtained the protein-ligand interaction fingerprints (IFPs) on their dissociation pathways and their influencing weights on the koff value. We observed a high correlation among the simulated, predicted, and experimental -log(koff) values. Combining ML, molecular dynamics (MD) simulation, and IFPs derived from accelerated MD helps design a drug for specific kinetic properties and selectivity profiles to the target of interest. To further validate our koff predictive ML model, we tested our model on two new N-HSP90 inhibitors, which have experimental koff values and are not in our ML training dataset. The predicted koff values are consistent with experimental data, and the mechanism of their kinetic properties can be explained by IFPs, which shed light on the nature of their selectivity against N-HSP90 protein. We believe that the ML model described here is transferable to predict koff of other proteins and will enhance the kinetics-based drug design endeavor.
RESUMEN
In recent years, machine learning (ML) models have been found to quickly predict various molecular properties with accuracy comparable to high-level quantum chemistry methods. One such example is the calculation of electrostatic potential (ESP). Different ESP prediction ML models were proposed to generate surface molecular charge distribution. Electrostatic complementarity (EC) can apply ESP data to quantify the complementarity between a ligand and its binding pocket, leading to the potential to increase the efficiency of drug design. However, there is not much research discussing EC score functions and their applicability domain. We propose a new EC score function modified from the one originally developed by Bauer and Mackey, and confirm its effectiveness against the available Pearson's R correlation coefficient. Additionally, the applicability domain of the EC score and two indices used to define the EC score application scope will be discussed.