In-Silico Study on the Interaction of Saffron Ligands and Beta-Lactoglobulin by Molecular Dynamics and Molecular Docking Approach

Safranal, crocetin, and dimethylcrocetin are secondary metabolites found in saffron and have a wide range of biological activities. An investigation of their interaction with a transport protein, such as β-lactoglobulin (β-lg), at the atomic level could be a valuable factor in controlling their transport to biological sites. The interaction of these ligands and β-lg as a transport protein was investigated using molecular docking and molecular dynamics (MD) simulation methods. The molecular docking results showed that safranal and crocetin bind on the surface of β-lg. However, dimethylcrocetin binds in the internal cavity of β-lg. The β-lg affinity for binding saffron ligands decreases in the following order: crocetin > dimethylcrocetin > safranal. The analysis of MD simulation trajectories showed that the β-lg and β-lg–ligand complexes became stable at approximately 3000 ps and that there was little conformational change in the β-lg–safranal and β-lg–dimethylcrocetin complexes over a 10-ns timescale. In addition, the profiles of atomic fluctuations showed the rigidity of the ligand binding site during the simulation time.     

Deoxyribonucleic Acid and Bovine Serum Albumin Interaction with the Asymmetric Schiff Base Ligand and Its Molybdenum (VI) Complex: Multi Spectroscopic and Molecular Docking Studies

The interaction of an asymmetric Schiff base ligand derived from allylamine and 2,3-dihydroxybenzaldehyde and its molybdenum (VI) complex with deoxyribonucleic acid (DNA) and bovine serum albumin (BSA) were studied using spectroscopic and molecular docking methods. The spectroscopic results revealed that the DNA and BSA affinity for binding the Mo(VI) complex is greater than its ligand. Furthermore, the molecular docking calculations showed that H-bond, hydrophobic, π-π and π-cation interactions had the dominant roles in the stability of the compound-BSA complexes. The DNA interaction results suggested that the compounds interacted with DNA by the groove binding mechanism.     

Spectro Analytical,Computational and In Vitro Biological Studies of Novel Substituted Quinolone Hydrazone and it’s Metal Complexes

Some novel transition metal [Cu (II), Ni (II) and Co (II)] complexes of nalidixic acid hydrazone have been prepared and characterized by employing spectro-analytical techniques viz: elemental analysis, ¹H-NMR, Mass, UV–Vis, IR, TGA-DTA, SEM-EDX, ESR and Spectrophotometry studies. The HyperChem 7.5 software was used for geometry optimization of title compound in its molecular and ionic forms. Quantum mechanical parameters, contour maps of highest occupied molecular orbitals (HOMO) and lowest unoccupied molecular orbitals (LUMO) and corresponding binding energy values were computed using semi empirical single point PM3 method. The stoichiometric equilibrium studies of metal complexes carried out spectrophotometrically using Job’s continuous variation and mole ratio methods inferred formation of 1:2 (ML₂) metal complexes in respective systems. The title compound and its metal complexes screened for antibacterial and antifungal properties, exemplified improved activity in metal complexes. The studies of nuclease activity for the cleavage of CT- DNA and MTT assay for in vitro cytotoxic properties involving metal complexes exhibited high activity. In addition, the DNA binding properties of Cu (II), Ni (II) and Co (II) complexes investigated by electronic absorption and fluorescence measurements revealed their good binding ability and commended agreement of K_b values obtained from both the techniques. Molecular docking studies were also performed to find the binding affinity of synthesized compounds with DNA (PDB ID: 1N37) and “Thymidine phosphorylase from E.coli” (PDB ID: 4EAF) protein targets.     

Docking accuracy enhanced by QM-derived protein charges

Effects of protein polari sation on docking accuracy were investigated using molecular docking programme AutoDock 4 in which topology-specific empirical Gasteiger charges were replaced with Polarised protein-specific charges (PPC) to represent quantum mechanics- polarised protein. Docking was successfully conducted for 50 diverse protein–ligand complexes. The docking with PPC charges shows a decrease in the root-mean-square deviation (RMSD) values of ligands compared to those from the docking with Gasteiger charges. Ligand binding orientations and their key interactions such as hydrogen bonding interactions in X-ray structures were substantially reproduced in complexes docked using PPC scheme with 98% of the RMSDs of the best docking poses less than 2 Å compared to 74% in the docking with Gasteiger charges. Considerable improvements in docking accuracy were observed by simply altering the atomic partial charges in the scoring function, which reflects the importance of protein polarisation in molecular docking. Further research can be carried out to (1) include polarisation of both ligands and proteins to account for polarisation effects within protein and between protein and ligand, and (2) develop a PPC-based scoring function to increase the docking accuracies for protein–ligand complexes from a larger variety of protein families.     

The Interaction of Polyphenol Flavonoids with β-lactoglobulin: Molecular Docking and Molecular Dynamics Simulation Studies

The interaction of quercetin, quercitrin, and rutin, as natural polyphenolic compounds, with β-lactoglobulin (BLG) using molecular docking and molecular dynamics simulation methods was examined. Molecular docking studies showed that quercetin and quercitrin were bounded to the internal cavity of protein, while rutin was bounded to the entrance of the cavity because of its large structural volume. It was found that there were one-, three-, and four-hydrogen bond interactions between BLG and quercetin, quercitrin, and rutin respectively. This showed that with an increase in the number of OH groups in the flavonoid structure, there was an increase in the number of hydrogen bond interactions. The binding constants for the binding of quercetin, quercitrin, and rutin to BLG were 1.2 × 10⁶, 1.9 × 10⁶, and 7.4 × 10⁴ M^?1 respectively. The results of molecular dynamics simulation showed that the root mean square deviation (RMSD) of non-liganded BLG and BLG–ligand complexes reached equilibration after 3500 ps. The study of the radius of gyration revealed that BLG and BLG–ligand complexes were stabilized around 2500 ps, and unlike the two other complexes, there was no conformational change for BLG–quercetin. Finally, analyzing the RMS fluctuations suggested that the structure of the ligand binding site remained approximately rigid during simulation.     

Computational Studies on the Interaction of Arctiin and Liquiritin With β-lactoglobulin

The study of the interaction of drugs purified from natural sources and a transport protein, such as β-lactoglobulin (BLG), at the atomic level could be a valuable factor to control their transport to biological sites. In the present study, molecular docking and molecular dynamics simulation methods were used to study the interaction of arctiin and liquiritin as natural drugs and BLG as the transport protein. The molecular docking results indicated that these drugs bind in the internal cavity of BLG and the BLG affinity for binding the liquiritin is greater than arctiin. The docking results also indicated that the hydrogen bond interactions have a dominant role in the BLG-drug complex stability. The analysis of MD simulation trajectories showed that the root mean square deviation (RMSD) of BLG-liquiritin, unliganded BLG, and BLG-arctiin reached equilibrium and fluctuated around the mean value at about 1000, 3500, and 4000 ps, respectively. The time evolution of the radius of gyration and total solvent accessible surface of the protein showed that BLG-arctiin and BLG-liquiritin complexes became stable around 2500 and 5000 ps, respectively. Also, the profiles of atomic fluctuations during the simulation showed the rigidity of the ligand binding sites.     

In Silico Study on the Interaction of Thiazolidinediones and β-Lactoglobulin by Molecular Dynamics and Docking Approach

Thiazolidinediones are widely used in the treatment of diabetes mellitus type 2. An investigation of their interaction with a transport protein, such as β-lactoglobulin (BLG), at the atomic level could be a valuable factor in controlling their transport to biological sites. The interaction of troglitazone, pioglitazone, and rosiglitazone, as representative thiazolidinediones, and BLG, as a transport protein, was investigated using molecular docking and molecular dynamics (MD) simulation methods. The molecular docking results showed that these thiazolidinediones bind to the internal cavity of BLG and the BLG affinity for binding the thiazolidinediones decreases in the following order: troglitazone > pioglitazone > rosiglitazone. The analysis of MD simulation trajectories showed that the BLG and BLG-thiazolidinedione complexes became stable at approximately 2500 ps and that there was little conformational change in the BLG-thiazolidinedione complexes over a 10 ns timescale. In addition, the profiles of atomic fluctuations showed the rigidity of the ligand-binding site during the simulation time.     

In silico screening of neurokinin receptor antagonists as a therapeutic strategy for neuroinflammation in Alzheimer’s disease

Neuroinflammation is one of the detrimental factors leading to neurodegeneration in Alzheimer’s disease (AD) and other neurodegenerative disorders. The activation of microglial neurokinin 1 receptor (NK1R) by substance P (SP) enhances neuroinflammation which is mediated through pro-inflammatory pathways involving NFkB, ERK1/2, and P38 and thus projects the scope and importance of NK1R inhibitors. Emphasizing the inhibitory role of N Acetyl l Tryptophan (l-NAT) on NK1R, this is the first in silico screening of l-NAT mediated NK1R antagonism. In addition, FDA- approved ligands were screened for their potential NK1R antagonism. The l-NAT was docked in XP (Extra Precision) mode while FDA-approved ligands were screened in HTVS (High Throughput Virtual Screening), SP (Standard Precision), and XP mode onto NK1R (PDB:6HLO). The l-NAT and top 3 compounds FDA-approved ligands were subjected to molecular dynamics (MD) studies of 100 ns simulation time. The XP docking of l-NAT, indacaterol, modafinil and alosetron showed good docking scores. Their 100 ns MD showed brief protein–ligand interactions with an acceptable root mean square deviation. The protein–ligand contacts depicted pi-pi stacking, pi-cation, hydrogen bonds, and water bridges with the amino acids necessary for NK1R inhibition. The variable colour band intensities on the protein–ligand contact map indicated their binding strength with amino acids. The molecular mechanics/generalized born surface area (MM-GBSA) scores suggested favourable binding free energy of the complexes. Thus, our study predicted the ability of l-NAT, indacaterol, modafinil, and alosetron as capable NK1R inhibitors that can aid to curb neuroinflammation in conditions of AD which could be further ascertained in subsequent studies.

Graphic Abstract

     

Structure-based virtual screening,in silico docking,ADME properties prediction and molecular dynamics studies for the identification of potential inhibitors against SARS-CoV-2 Mpro

COVID-19 is a viral pandemic caused by SARS-CoV-2. Due to its highly contagious nature, millions of people are getting affected worldwide knocking down the delicate global socio-economic equilibrium. According to the World Health Organization, COVID-19 has affected over 186 million people with a mortality of around 4 million as of July 09, 2021. Currently, there are few therapeutic options available for COVID-19 control. The rapid mutations in SARS-CoV-2 genome and development of new virulent strains with increased infection and mortality among COVID-19 patients, there is a great need to discover more potential drugs for SARS-CoV-2 on a priority basis. One of the key viral enzymes responsible for the replication and maturation of SARS-CoV-2 is M^pro protein. In the current study, structure-based virtual screening was used to identify four potential ligands against SARS-CoV-2 M^pro from a set of 8,722 ASINEX library compounds. These four compounds were evaluated using ADME filter to check their ADME profile and druggability, and all the four compounds were found to be within the current pharmacological acceptable range. They were individually docked to SARS-CoV-2 M^pro protein to assess their molecular interactions. Further, molecular dynamics (MD) simulations was carried out on protein–ligand complex using Desmond at 100 ns to explore their binding conformational stability. Based on RMSD, RMSF and hydrogen bond interactions, it was found that the stability of protein–ligand complex was maintained throughout the entire 100 ns simulations for all the four compounds. Some of the key ligand amino acid residues participated in stabilizing the protein–ligand interactions includes GLN 189, SER 10, GLU 166, ASN 142 with PHE 66 and TRP 132 of SARS-CoV-2 M^pro. Further optimization of these compounds could lead to promising drug candidates for SARS-CoV-2 M^pro target.

     

Integrating sampling techniques and inverse virtual screening: toward the discovery of artificial peptide-based receptors for ligands

A novel heuristic using an iterative select-and-purge strategy is proposed. It combines statistical techniques for sampling and classification by rigid molecular docking through an inverse virtual screening scheme. This approach aims to the de novo discovery of short peptides that may act as docking receptors for small target molecules when there are no data available about known association complexes between them. The algorithm performs an unbiased stochastic exploration of the sample space, acting as a binary classifier when analyzing the entire peptides population. It uses a novel and effective criterion for weighting the likelihood of a given peptide to form an association complex with a particular ligand molecule based on amino acid sequences. The exploratory analysis relies on chemical information of peptides composition, sequence patterns, and association free energies (docking scores) in order to converge to those peptides forming the association complexes with higher affinities. Statistical estimations support these results providing an association probability by improving predictions accuracy even in cases where only a fraction of all possible combinations are sampled. False positives/false negatives ratio was also improved with this method. A simple rigid-body docking approach together with the proper information about amino acid sequences was used. The methodology was applied in a retrospective docking study to all 8000 possible tripeptide combinations using the 20 natural amino acids, screened against a training set of 77 different ligands with diverse functional groups. Afterward, all tripeptides were screened against a test set of 82 ligands, also containing different functional groups. Results show that our integrated methodology is capable of finding a representative group of the top-scoring tripeptides. The associated probability of identifying the best receptor or a group of the top-ranked receptors is more than double and about 10 times higher, respectively, when compared to classical random sampling methods.     

The Interaction of a New Schiff Base Ligand with Human Serum Albumin: Molecular Docking and Molecular Dynamics Simulation Studies

The interaction of a new heterocyclic Schiff base bearing pyridine and pyrimidine cycles, with human serum albumin (HSA) using molecular docking and molecular dynamics simulation methods was examined. Molecular docking studies showed that the ligand was bonded to the IB domain of the protein. It was found that there was one hydrogen bond interaction between HSA and the ligand. The standard Gibbs free energy for binding of the ligand to HSA was calculated as ?9.63 kcal.mol^?1. The results of the molecular dynamics simulation showed that the root mean square deviation (RMSD) of the non-liganded HSA and the HSA–ligand complex reached equilibration after 1000 ps. The study of the radius of gyration revealed that there was a conformational change when the HSA–ligand complex was formed. Finally, analyzing the RMS fluctuations (RMSF) suggested that the structure of the ligand binding site remained approximately rigid during the simulation.     

2D-SAR and 3D-QSAR analyses for acetylcholinesterase inhibitors

Alzheimer’s disease (AD) accounts for almost three quarters of dementia patients and interferes people’s normal life. Great progress has been made recently in the study of Acetylcholinesterase (AChE), known as one of AD’s biomarkers. In this study, acetylcholinesterase inhibitors (AChEI) were collected to build a two-dimensional structure–activity relationship (2D-SAR) model and three-dimensional quantitative structure–activity relationship (3D-QSAR) model based on feature selection method combined with random forest. After calculation, the prediction accuracy of the 2D-SAR model was 89.63% by using the tenfold cross-validation test and 87.27% for the independent test set. Three cutting ways were employed to build 3D-QSAR models. A model with the highest \({q}^{2}\) (cross-validated correlation coefficient) and \({r}^{2 }\)(non-cross-validated correlation coefficient) was obtained to predict AChEI activity. The mean absolute error (MAE) of the training set and the test set was 0.0689 and 0.5273, respectively. In addition, molecular docking was also employed to reveal that the ionization state of the compounds had an impact upon their interaction with AChE. Molecular docking results indicate that Ser124 might be one of the active site residues.     

Synthesis,Characterization, DNA Binding Properties and Antioxidant Activity of Ln(III) Complexes with Schiff Base Ligand Derived from 3-Carbaldehyde Chromone and Aminophenazone

A novel Schiff base ligand, chromone-3-carbaldehyde-aminophenazone (L) and its Ln(III) (Ln = La, Yb) complexes were synthesized and characterized by physicochemical methods. The interaction between the ligand, Ln(III) complexes and calf thymus DNA in physiological buffer (pH = 7.10) was investigated by using UV–vis spectroscopy, fluorescence spectra, ethidium bromide experiments and viscosity measurements, indicating that the studied compounds can all bind to DNA via an intercalation binding mode and the complexes have stronger binding affinity than the free ligand alone. Furthermore, antioxidant activity of the ligand and its complexes was determined by superoxide and hydroxyl radical scavenging methods in vitro, suggesting that Ln(III) complexes inhibit stronger antioxidant activity than the ligand alone and some standard antioxidants, such as mannitol and vitamin C.     

苯并咪唑/苯并噻唑类T315I突变型Abl抑制剂的结合机理和分子设计

Despite the efficacy of imatinib therapy in chronic myelogenous leukemia, the development of drug-resistant Abl mutants, especially the most difficult overcoming T315I mutant, makes the search for new Abl T315I inhibitors a very interesting challenge in medicinal chemistry. In this work, a multistep computational framework combining the three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, molecular dynamics (MD) simulation and binding free energy calculation, was performed to explore the structural requirements for the Abl T315I activities of benzimidazole/benzothiazole derivatives and the binding mechanism between the inhibitors and Abl T315I. The established 3D-QSAR models exhibited satisfactory internal and external predictability. Docking study elucidated the comformations of compounds and the key amino acid residues at the binding pocket, which were confirmed by MD simulation. The binding free energies correlated well with the experimental activities. The MM-GBSA energy decomposition revealed that the van der Waals interaction was the major driving force for the interaction between the ligands and Abl T315I. The hydrogen bond interactions between the inhibitors and Met318 also played an important role in stablizing the binding of compounds to Abl T315I. Finally, four new compounds with rather high Abl T315I activities were designed and presented to experimenters for reference.     

A Study on Spectro-Analytical Aspects,DNA – Interaction,Photo-Cleavage,Radical Scavenging,Cytotoxic Activities,Antibacterial and Docking Properties of 3 – (1 – (6 – methoxybenzo [d] thiazol – 2 – ylimino) ethyl) – 6 – methyl – 3H – pyran - 2, 4 – dione and its Metal Complexes

The focus of the present work is on the design, synthesis, characterization, DNA-interaction, photo-cleavage, radical scavenging, in-vitro cytotoxicity, antimicrobial, docking and kinetic studies of Cu (II), Cd (II), Ce (IV) and Zr (IV) metal complexes of an imine derivative, 3 – (1 – (6 – methoxybenzo [d] thiazol – 2 – ylimino) ethyl) – 6 – methyl – 3H – pyran – 2, 4 – dione. The investigation of metal ligand interactions for the determination of composition of metal complexes, corresponding kinetic studies and antioxidant activity in solution was carried out by spectrophotometric methods. The synthesized metal complexes were characterized by EDX analysis, Mass, IR, ¹H-NMR, ¹³C-NMR and UV–Visible spectra. DNA binding studies of metal complexes with Calf thymus (CT) DNA were carried out at room temperature by employing UV-Vis electron absorption, fluorescence emission and viscosity measurement techniques. The results revealed that these complexes interact with DNA through intercalation. The results of in vitro antibacterial studies showed the enhanced activity of chelating agent in metal chelated form and thus inferring scope for further development of new therapeutic drugs. Cell viability experiments indicated that all complexes showed significant dose dependent cytotoxicity in selected cell lines. The molecular modeling and docking studies were carried out with energy minimized structures of metal complexes to identify the receptor to metal interactions.     

High-throughput screening against $$\sim $$6.1 million structurally diverse,lead-like compounds to discover novel ROCK inhibitors for cerebral injury recovery

Rho-associated protein kinase (ROCK) has been recognized as an attractive therapeutic target to promote neurogenesis, neuroregeneration, and neurorecovery after cerebral injury. Here, a high-throughput screening protocol was described to discover novel ROCK inhibitors from a large chemical library containing \(\sim \)6.1 million structurally diverse, lead-like compounds. The protocol employed empirical rules such as ADMET evaluation and chemical similarity analysis to exclude those of drug-unlike candidates, and then molecular docking and binding affinity predictions were performed to suggest few promising candidates with high scores. Consequently, five compounds were successfully identified to have satisfactory activity profile with \(\hbox {IC}_{50}\) values at nanomolar level. In order to elucidate the molecular mechanism of inhibitor binding to target, the complex structures of ROCK kinase domain with the five identified compounds were modeled and examined in detail. A number of polar chemical forces such as hydrogen bonds and cation-\(\pi \) interactions as well as nonpolar contacts such as \(\pi \)–\(\pi \) stacking and hydrophobic forces were revealed at the complex interface, conferring high affinity and strong specificity to inhibitor binding. In addition, several key residues around the kinase active site, including Val90, Lys105, Asn203, and Phe368, were found to play an important role in binding.