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.     

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.     

Insights into the interaction and binding mode of a set of antifungal azoles as inhibitors of potential fungal enzyme-based targets

Azole-containing compounds are a kind of chemical entities of natural and synthetic origin having a wide-range of activities. They are therefore considered as important moieties for fungicide development, mostly due to the possible action on several enzyme-based targets. As part of our research on fungicidal agents, the relationship between the ligand-enzyme affinities of several synthetic azole-containing compounds against a set of fungal enzyme-based targets was in silico evaluated through molecular docking. The affinity values of the test compounds were mostly higher than those of the respective test controls. Binding modes between enzymes and test compounds were firstly investigated through Vina scores and ligand–residue interactions. Furthermore, statistically relationships among docking scores were successfully found by multivariate analysis. They were mostly correlated with reported MIC₈₀ values, so it denoted an evident discrimination of the test compounds. Strong electron withdrawing groups on phenylacrylamide moiety were responsible for establishing stronger complexes with the enzyme targets, being trichodiene synthase and α-l-fucosidase the most important ones. Moreover, stability of a set of representative protein/ligand complexes was also analyzed by 10 ns molecular dynamics simulations (MD). Significant differences into the MD runs were detected and directly correlated to docking performances. Finally, docking affinity scores and HOMO–LUMO energy gaps resulted well predicted by comparative molecular field analysis (CoMFA) models, demonstrating the structure type is particularly associated with those calculated properties and these results were thus consistent with the respective validation parameters.     

Spectroscopic and docking studies of the binding of two stereoisomeric antioxidant catechins to serum albumins

The interactions of two stereoisomeric antioxidant flavonoids, catechin (C) and epicatechin (EC) with bovine serum albumin (BSA) and human serum albumin (HSA), have been investigated by steady state and time resolved fluorescence, phosphorescence, circular dichroism (CD), FTIR and protein–ligand docking studies. The steady-state fluorescence studies indicate a single binding site for both the ligands. FTIR spectra suggest that in both the albumins, C and EC stabilize the α-helix at the cost of a corresponding loss in the β-sheet structure. CD studies have been carried out using (±)C, and both the epimers (+)C and (?)C. The low temperature phosphorescence and protein–ligand [(+), (?) and (±) forms of C and EC] docking studies indicate that the ligands bind in the proximity of Trp 134 of BSA and Trp 214 of HSA, thereby changing their solvent accessible surface areas (ASA). Asn 158 and Glu 130 side chains are found to be within the hydrogen bonding distance from the phenolic –OH groups of C and EC in the case of BSA complex. C and EC are located within the binding pocket of sub-domain IIa of HSA.     

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.     

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.     

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.     

Thermally Stable Schiff Base and its Metal Complexes: Molecular Docking and Protein Binding Studies

In this paper, interaction of Schiff base and its metal complexes carrying naphthalene ring in the structure with bovine serum albumin (BSA) were investigated using UV-vis absorption, fluorescence spectroscopies and molecular docking methods. The effect on the binding mechanism and properties of these compounds containing metal-free, iron and copper ions were also investigated. The fluorescence spectroscopy results showed that fluorescence intensity of BSA in the presence of different concentration of ligands was decreased through a static quenching mechanism. Binding constants (KSV, Kbin and Ka) and thermodynamic parameters (ΔG, ΔH and ΔS) for the ligand-protein interactions were also determined. ΔG values of ligand-protein interaction were calculated in the range ? 6.3 to ?5.5 kcal/mol. These negative values showed that binding process is spontaneous and, hydrogen bonding and van der Waals force were main interaction of the protein and ligands. ΔH and ΔS value were also calculated in the range of 1.10 to 1.26 kJ/mol and 0.133 to 0.135 kJ/mol. K, respectively. These positive values indicated that the binding process between ligands and BSA are endothermic and electrostatic interaction, respectively.     

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.     

Studies on Comonomer Compositional Distribution of Poly(propylene carbonate-propylene oxide) Copolymer and Its Effect on the Thermal,Mechanical and Oxygen Barrier Properties of Fractions

Poly(propylene carbonate) (PPC) was synthesized by the alternating copolymerization of carbon dioxide and propylene oxide (PO). However, during the polymerization, two by-products tended to produce cyclic propylene carbonate (CPC) and a polyether (PE) segment. The excess PO repeat units (PE segment) can easily insert into the PPC backbone and eventually produce the PPC–PO copolymer. The production of CPC and PE segments affected the increase of polymer chain length. In order to investigate the effects of the existence of PE segments, CPC, and molecular weight of fractions on the physical properties of PPC–PO copolymer, a series of fractions with narrow molecular weight distribution were obtained by repeated fractionation. Based on a solvent/non-solvent (chloroform/n-heptane) mixture, an original PPC–PO sample was fractionated into nine fractions with number–average molecular weights (M_n) from 0.34 × 10⁵ to 5.56 × 10⁵ and PE content from 0.4 to 15 mol%. The M_n of PPC–PO fractions decreased with the increase of PE content in the PPC–PO backbone, and the thermal and mechanical properties of the PPC–PO copolymers were affected by their M_n and PE contents. Furthermore, the lower the PE content, the higher the M_n value. Higher M_n means better tensile and lower oxygen permeability of PPC–PO copolymer.     

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.     

EPR Study of Ligand Effects in Spin Triads of Bis(o-Semiquinonato)Copper(II) Complexes

We report the X-band (9 GHz) electron paramagnetic resonance (EPR) study of series of bis(o-semiquinonato)copper(II) complexes with different ligands. It was found previously, that exchange interactions in spin triads of these compounds are very sensitive to the structure of the ligand coordinated to the central copper(II) ion. Ligand moderates the copper–radical and radical–radical exchange interactions and strongly changes the magnetic properties of the compound. Depending on a ligand, ferromagnetic or antiferromagnetic character of exchange dominates in the system. The EPR study of these complexes allowed us to obtain information on zero-field splitting parameters and their distributions in the studied compounds. The EPR results compliment previously obtained spectroscopic data on these compounds and suggest the pronounced plasticity of the clusters manifested in the broad distributions of their rhombicity parameters.     

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

     

CORCEMA refinement of the bound ligand conformation within the protein binding pocket in reversibly forming weak complexes using STD-NMR intensities

We describe an intensity-restrained optimization procedure for refining approximate structures of ligands within the protein binding pockets using STD-NMR intensity data on reversibly forming weak complexes. In this approach, the global minimum for the bound-ligand conformation is obtained by a hybrid structure refinement method involving CORCEMA calculation of intensities and simulated annealing optimization of torsion angles of the bound ligand using STD-NMR intensities as experimental constraints and the NOE R-factor as the pseudo-energy function to be minimized. This method is illustrated using simulated STD data sets for typical carbohydrate and peptide ligands. Our procedure also allows for the optimization of side chain torsion angles of protein residues within the binding pocket. This procedure is useful in refining and improving initial models based on crystallography or computer docking or other algorithms to generate models for the bound ligand (e.g., a lead compound) within the protein binding pocket compatible with solution STD-NMR data. This method may facilitate structure-based drug design efforts.     

Organic light‐emitting materials based on iridium(III) complexes bearing phenanthroimidazole ligands

We present the elegant synthesis and the photophysical and electroluminescent properties of a series of cyclometalated iridium(III) complexes [Ir(PPI)₂(pic), PPI: 1,2‐diphenyl‐1H‐phenanthro[9,10‐d]imidazole; pic: picolinic acid]. The Ir(PPI)₂(pic) complexes showed characteristic phosphorescence with an emission range of 556–579 nm and a high quantum efficiency with microsecond lifetimes. The strongly allowed phosphorescence in these complexes is the result of significant spin–orbit coupling of the Ir center. All bis(PPI) derivatives exhibit intense triplet metal‐to‐ligand charge transfer (MLCT) photoluminescence in the fluid solutions at room temperature. The impact of different solvents, substituents on the phenanthroimidazole ligands and complex concentrations upon their emissive behavior have been examined and demonstrate that their emission energies can be systematically modified. Weak bands located at longer wavelength have been assigned to the ¹MLCT ← S₀ and ³MLCT ← S₀ transitions of iridium complexes. Application of the ³MLCT excited state of the [Ir(PPI)₂(pic)] materials in organic light‐emitting devices are described. Copyright © 2014 John Wiley & Sons, Ltd.     

Computational predictive models for P-glycoprotein inhibition of in-house chalcone derivatives and drug-bank compounds

The human P-glycoprotein (P-gp) efflux pump is of great interest for medicinal chemists because of its important role in multidrug resistance (MDR). Because of the high polyspecificity as well as the unavailability of high-resolution X-ray crystal structures of this transmembrane protein, ligand-based, and structure-based approaches which were machine learning, homology modeling, and molecular docking were combined for this study. In ligand-based approach, individual two-dimensional quantitative structure–activity relationship models were developed using different machine learning algorithms and subsequently combined into the Ensemble model which showed good performance on both the diverse training set and the validation sets. The applicability domain and the prediction quality of the developed models were also judged using the state-of-the-art methods and tools. In our structure-based approach, the P-gp structure and its binding region were predicted for a docking study to determine possible interactions between the ligands and the receptor. Based on these in silico tools, hit compounds for reversing MDR were discovered from the in-house and DrugBank databases through virtual screening using prediction models and molecular docking in an attempt to restore cancer cell sensitivity to cytotoxic drugs.     

Origin of accidental degeneracy in ligand-field splittings of substituted octahedral complexes

Accidental degeneracy seems to be the rule rather than an exception amongst thed orbital energies of substituted octahedral complexes ofd ¹ configuration. By using symmetry and physical arguments, in conjunction with first-order and second-order degenerate perturbation theory, it is shown that such accidental degeneracies arise in crystal-field theory due to the choice of an inflexible basis set of metal orbitals which neglects the polarisation of metal orbitals by the ligand charges.