PLAS-20k: Extended Dataset of Protein-Ligand Affinities from MD Simulations for Machine Learning Applications.

Computing binding affinities is of great importance in drug discovery pipeline and its prediction using advanced machine learning methods still remains a major challenge as the existing datasets and models do not consider the dynamic features of protein-ligand interactions. To this end, we have developed PLAS-20k dataset, an extension of previously developed PLAS-5k, with 97,500 independent simulations on a total of 19,500 different protein-ligand complexes. Our results show good correlation with the available experimental values, performing better than docking scores. This holds true even for a subset of ligands that follows Lipinski's rule, and for diverse clusters of complex structures, thereby highlighting the importance of PLAS-20k dataset in developing new ML models. Along with this, our dataset is also beneficial in classifying strong and weak binders compared to docking. Further, OnionNet model has been retrained on PLAS-20k dataset and is provided as a baseline for the prediction of binding affinities. We believe that large-scale MD-based datasets along with trajectories will form new synergy, paving the way for accelerating drug discovery.

Revealing the drug resistance mechanism of saquinavir due to G48V and V82F mutations in subtype CRF01_AE HIV-1 protease: molecular dynamics simulation and binding free energy calculations.

Human immunodeficiency virus-1 (HIV-1) protease is one of the important targets in AIDS therapy. The majority of HIV infections are caused due to non-B subtypes in developing countries. The co-occurrence of mutations along with naturally occurring polymorphisms in HIV-1 protease cause resistance to the FDA approved drugs, thereby posing a major challenge in the treatment of antiretroviral therapy. In this work, the resistance mechanism against SQV due to active site mutations G48V and V82F in CRF01_AE (AE) protease was explored. The binding free energy calculations showed that the direct substitution of valine at position 48 introduces a bulkier side chain, directly impairing the interaction with SQV in the binding pocket. Also, the intramolecular hydrogen bonding network of the neighboring residues is altered, indirectly affecting the binding of SQV. Interestingly, the substitution of phenylalanine at position 82 induces conformational changes in the 80's loop and the flap region, thereby favoring the binding of SQV. The V82F mutant structure also maintains similar intramolecular hydrogen bond interactions as observed in AE-WT.Communicated by Ramaswamy H. Sarma.

Exploring the drug resistance mechanism of active site, non-active site mutations and their cooperative effects in CRF01_AE HIV-1 protease: molecular dynamics simulations and free energy calculations.

Human immunodeficiency virus type 1 protease is essential for virus replication and maturation and has been considered as one of the important drug target for the antiretroviral treatment of HIV infection. The majority of HIV infections are caused due to non-B subtypes in developing countries. Subtype AE is spreading rapidly and infecting huge population worldwide. Understanding the interdependence of active and non-active site mutations in conferring drug resistance is crucial for the development effective inhibitors in subtype AE protease. In this work, we have investigated the mechanism of resistance against indinavir (IDV) due to therapy selected active site mutation V82F, non-active site mutations PF82V and their cooperative effects PV82F in subtype AE-protease using molecular dynamics simulations and binding free energy calculations. The simulations suggested all the three complexes lead to decrease in binding affinity of IDV, whereas the PF82V complex resulted in an enhanced binding affinity compared to V82F and PV82F complexes. Large positional deviation of IDV was observed in V82F complex. The preservation of hydrogen bonds of IDV with active site Asp25/Asp25' and flap residue Ile50/50' via a water molecule is crucial for effective binding. Owing to the close contact of 80s loop with Ile50' and Asp25, the alteration between residues Thr80 and Val82, further induces conformational change thereby resulting in loss of interactions between IDV and the residues in the active site cavity, leading to drug resistance. Our present study shed light on the effect of active, non-active site mutations and their cooperative effects in AE protease. Communicated by Ramaswamy H. Sarma.

Evaluation of DNA/Protein interactions and cytotoxic studies of copper(II) complexes incorporated with N, N donor ligands and terpyridine ligand.

A series of four new copper(II) heteroleptic complexes, [Cu(2‴-pytpy) (L)] (NO3)2·2H2O (1-4), where 2‴-pytpy=4'-(2'''-Pyridyl)-2, 2':6', 2''-terpyridine, L=bipyridyl (bpy), 1, 10 phenanthroline(phen), dipyridoquinoxaline(dpq) and dipyridophenazine (dppz) were synthesized and characterized by spectroscopic techniques. Further, the molecular structure of the complex (2) was confirmed by single crystal X-ray diffraction technique and the data revealed a penta coordinated, distorted square-pyramidal geometry with triclinic system. The interactions of four complexes with calf thymus DNA and bovine serum albumin (BSA) were investigated by electronic absorption, fluorescence and circular dichroism spectroscopy techniques. Spectral studies substantiated an intercalative binding mode of metal complexes with ct-DNA. Significant binding interactions of the complexes with protein have been further revealed from fluorescence studies. Furthermore, all the four complexes show potential cytotoxicity towards the human liver carcinoma cell line (HepG-2).