In Silico Prediction of the Dissociation Rate Constants of Small Chemical Ligands by 3D-Grid-Based VolSurf Method

Mixture-of-Experts Based Dissociation Kinetic Model for De Novo Design of HSP90 Inhibitors with Prolonged Residence Time

The dissociation rate constant (koff) significantly impacts the drug potency and dosing frequency. This work proposes a powerful optimization-based framework for de novo drug design guided by koff. First, a comprehensive database containing 2,773 unique koff values is created. Based on the database, a novel generic dissociation kinetic model is developed with a mixture-of-experts architecture, enabling high-throughput predictions of koff with high accuracy. The developed model is then integrated with an optimization-based mathematical programming approach to design drug candidates with low koff. Finally, the τ-RAMD method is utilized to rigorously verify the designed potential drug candidates. In a case study, the framework successfully identified numerous new potential HSP90 inhibitor candidates, achieving a maximum 45.7% improvement in residence time (τ = 1/koff) compared to that of a known exceptional HSP90 inhibitor. These findings demonstrate the feasibility and effectiveness of the kinetics-guided optimization-based de novo drug design framework in designing drug candidates with prolonged τ.

The translational value of ligand-receptor binding kinetics in drug discovery

The translation of in vitro potency of a candidate drug, as determined by traditional pharmacology metrics (such as EC50/IC50 and KD/Ki values), to in vivo efficacy and safety is challenging. Residence time, which represents the duration of drug-target interaction, can be part of a more comprehensive understanding of the dynamic nature of drug-target interactions in vivo, thereby enabling better prediction of drug efficacy and safety. As a consequence, a prolonged residence time may help in achieving sustained pharmacological activity, while transient interactions with shorter residence times may be favourable for targets associated with side effects. Therefore, integration of residence time into the early stages of drug discovery and development has yielded a number of clinical candidates with promising in vivo efficacy and safety profiles. Insights from residence time research thus contribute to the translation of in vitro potency to in vivo efficacy and safety. Further research and advances in measuring and optimizing residence time will bring a much-needed addition to the drug discovery process and the development of safer and more effective drugs. In this review, we summarize recent research progress on residence time, highlighting its importance from a translational perspective.

Unveiling Glucosinolate Diversity in Brassica Germplasm and In Silico Analysis for Determining Optimal Antioxidant Potential

This study explored the glucosinolate (GSL) content in Brassica plants and utilized in silico analysis approach to assess their antioxidant capabilities. GSLs, present abundantly in Brassica vegetables, offer potential health advantages, including antioxidant effects. Employing Ultra-Performance Liquid Chromatography (UPLC) coupled with tandem mass spectrometry (MS/MS), major GSLs were identified in 89 accessions from diverse species and subspecies. Statistical analysis and principal component analysis unveiled significant GSL variation and potential correlations among the Brassica germplasms. This study unveils the dominance of aliphatic GSLs over aromatic and indolyl compounds in all the accessions. Notably, Gluconapin (GNA) (33,049.23 µmol·kg-1 DW), Glucobrassicanapin (GBN) (9803.82 µmol·kg-1 DW), Progoitrin (PRO) (12,780.48 µmol·kg-1 DW) and Sinigrin (SIN) (14,872.93 µmol·kg-1 DW) were the most abundant compounds across the analyzed accessions. Moreover, in silico docking studies predicted promising antioxidant activity by evaluating the interactions of each GSL with antioxidant enzymes. Specifically, Sinigrin and Gluconapin exhibited a notably weaker influence on antioxidant enzymes. This provides key insights into the antioxidant potential of Brassica germplasm and highlights the importance of in silico analysis for evaluating bioactive properties. In general, the results of this study could be utilized in breeding programs to maximize GSL levels and antioxidant properties in Brassica crops and for developing functional foods with enhanced health benefits.

Structure-Kinetic Relationship for Drug Design Revealed by a PLS Model with Retrosynthesis-Based Pre-Trained Molecular Representation and Molecular Dynamics Simulation

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 (k_off) 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 k_off value. We observed a high correlation among the simulated, predicted, and experimental -log(k_off) 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 k_off predictive ML model, we tested our model on two new N-HSP90 inhibitors, which have experimental k_off values and are not in our ML training dataset. The predicted k_off 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 k_off of other proteins and will enhance the kinetics-based drug design endeavor.

Advances in computational methods for ligand binding kinetics

Binding kinetic parameters can be correlated with drug efficacy, which in recent years led to the development of various computational methods for predicting binding kinetic rates and gaining insight into protein-drug binding paths and mechanisms. In this review, we introduce and compare computational methods recently developed and applied to two systems, trypsin-benzamidine and kinase-inhibitor complexes. Methods involving enhanced sampling in molecular dynamics simulations or machine learning can be used not only to predict kinetic rates, but also to reveal factors modulating the duration of residence times, selectivity, and drug resistance to mutations. Methods which require less computational time to make predictions are highlighted, and suggestions to reduce the error of computed kinetic rates are presented.

The Physical Chemistry and Chemical Physics (PCCP) Section of the International Journal of Molecular Sciences in Its Publications: The First 300 Thematic Articles in the First 3 Years

The Physical Chemistry and Chemical Physics Section (PCCP Section) is one of the youngest among the sections of the International Journal of Molecular Sciences (IJMS)-the year 2021 will only mark three years since its inception [...].

Prediction of the Drug-Target Binding Kinetics for Flexible Proteins by Comparative Binding Energy Analysis

There is growing consensus that the optimization of the kinetic parameters for drug-protein binding leads to improved drug efficacy. Therefore, computational methods have been developed to predict kinetic rates and to derive quantitative structure-kinetic relationships (QSKRs). Many of these methods are based on crystal structures of ligand-protein complexes. However, a drawback is that each ligand-protein complex is usually treated as having a single structure. Here, we present a modification of COMparative BINding Energy (COMBINE) analysis, which uses the structures of ligand-protein complexes to predict binding parameters. We introduce the option of using multiple structures to describe each ligand-protein complex in COMBINE analysis and apply this to study the effects of protein flexibility on the derivation of dissociation rate constants (k_off) for inhibitors of p38 mitogen-activated protein (MAP) kinase, which has a flexible binding site. Multiple structures were obtained for each ligand-protein complex by performing docking to an ensemble of protein configurations obtained from molecular dynamics simulations. Coefficients to scale ligand-protein interaction energies determined from energy-minimized structures of ligand-protein complexes were obtained by partial least squares regression, and they allowed for the computation of k_off values. The QSKR model obtained using single, energy-minimized crystal structures for each ligand-protein complex had higher predictive power than the QSKR model obtained with multiple structures from ensemble docking. However, incorporation of ligand-protein flexibility helped to highlight additional ligand-protein interactions that lead to longer residence times, such as interactions with residues Arg67 and Asp168, which are close to the ligand in many crystal structures. These results show that COMBINE analysis is a promising method to guide the design of compounds that bind to flexible proteins with improved binding kinetics.

Preventing the Interaction between Coronaviruses Spike Protein and Angiotensin I Converting Enzyme 2: An In Silico Mechanistic Case Study on Emodin as a Potential Model Compound

Emodin, a widespread natural anthraquinone, has many biological activities including health-protective and adverse effects. Amongst beneficial effects, potential antiviral activity against coronavirus responsible for the severe acute respiratory syndrome outbreak in 2002–2003 has been described associated with the inhibition of the host cells target receptors recognition by the viral Spike protein. However, the inhibition mechanisms have not been fully characterized, hindering the rational use of emodin as a model compound to develop more effective analogues. This work investigates emodin interaction with the Spike protein to provide a mechanistic explanation of such inhibition. A 3D molecular modeling approach consisting of docking simulations, pharmacophoric analysis and molecular dynamics was used. The plausible mechanism is described as an interaction of emodin at the protein–protein interface which destabilizes the viral protein-target receptor complex. This analysis has been extended to the Spike protein of the coronavirus responsible for the current pandemic hypothesizing emodin’s functional conservation. This solid knowledge-based foothold provides a possible mechanistic rationale of the antiviral activity of emodin as a future basis for the potential development of efficient antiviral cognate compounds. Data gaps and future work on emodin-related adverse effects in parallel to its antiviral pharmacology are explored.