Developing imidazole analogues as potential inhibitor for Leishmania donovani trypanothione reductase: virtual screening,molecular docking,dynamics and ADMET approach

Visceral leishmaniasis (VL) affects Indian subcontinent, African and South American continent, and it covers 70 countries worldwide. Visceral form of leishmaniasis is caused by Leishmania donovani in Indian subcontinent which is lethal if left untreated. Extensive resistance to antileishmanial drugs such as sodium stibogluconate, pentamidine and miltefosine and their decreased efficacy has been reported in the endemic region. Amphotericin B drug has shown good antileishmanial activity with significant toxicity, but its cost of treatment has limited the outreach of this treatment to affected people living in endemic zone. So, there is an urgent need to identify new antileishmanial drugs with excellent activity and minimal toxicity issues. Trypanothione reductase, a component of antioxidant system, is necessary for parasite growth and survival to raise infection. To develop potential inhibitor, we docked nine hundred and eighty-four 5-nitroimidazole analogues along with clomipramine which is a well-known inhibitor for TR. Total one hundred and forty-seven 5-nitroimidazole analogues with better docking score than clomipramine were chosen for ADMET and QikProp studies. Among these imidazole analogues, total twenty-four imidazole analogues and clomipramine were chosen on the basis of their ADMET, QikProp, and prime MM-GBSA study. Later on, two analogues with best MM-GBSA dG bind were undergone molecular dynamic simulation to ensure protein–ligand interactions. Using above approach, we confirm that ethyl 2-acetyl-5-[4-butyl-2-(3-hydroxypentyl)-5-nitro-1H-imidazol-1-yl]pent-2-enoate can be a drug candidate against L. donovani for the treatment of VL in the Indian subcontinent.     

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800?K, for 25?ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500?K has several conformations reaching into the protein core, whereas at 600?K some of the alpha-helices were lost. At 800?K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.     

Evidence of salicylic acid pathway with EDS1 and PAD4 proteins by molecular dynamics simulation for grape improvement

Biotic stress is a major cause of heavy loss in grape productivity. In order to develop biotic stress-resistant grape varieties, the key defense genes along with its pathway have to be deciphered. In angiosperm plants, lipase-like protein phytoalexin deficient 4 (PAD4) is well known to be essential for systemic resistance against biotic stress. PAD4 functions together with its interacting partner protein enhanced disease susceptibility 1 (EDS1) to promote salicylic acid (SA)-dependent and SA-independent defense pathway. Existence and structure of key protein of systemic resistance EDS1 and PAD4 are not known in grapes. Before SA pathway studies are taken in grape, molecular evidence of EDS1: PAD4 complex is to be established. To establish this, EDS1 protein sequence was retrieved from NCBI and homologous PAD4 protein was generated using Arabidopsis thaliana as template and conserved domains were confirmed. In this study, computational methods were used to model EDS1 and PAD4 and simulated the interactions of EDS1 and PAD4. Since no structural details of the proteins were available, homology modeling was employed to construct three-dimensional structures. Further, molecular dynamic simulations were performed to study the dynamic behavior of the EDS1 and PAD4. The modeled proteins were validated and subjected to molecular docking analysis. Molecular evidence of stable complex of EDS1:PAD4 in grape supporting SA defense pathway in response to biotic stress is reported in this study. If SA defense pathway genes are explored, then markers of genes involved can play pivotal role in grape variety development especially against biotic stress leading to higher productivity.     

A comparison study of the interaction between β-lactoglobulin and retinol at two different conditions: spectroscopic and molecular modeling approaches

The interaction between bovin β-Lactoglobulin (β-LG) and retinol at two different pH values was investigated by multispectroscopic, zeta potential, molecular modeling, and conductometry measurements. The steady state and polarization fluorescence spectroscopy revealed that complex formation at two different pH values could occur through a remarkable static quenching. According to fluorescence quenching, one set of binding site at pH 2 and two sets of binding sites at pH 7 were introduced for binding of retinol to β-LG that show the enhancement of saturation score of β-LG to retinol in dimmer condition. The polarization fluorescence analysis represented that there is more affinity between β-LG and retinol at pH 7 rather than at pH 2. The effect of retinol on β-LG was studied by UV-visible, circular dichroism (CD), and synchronous fluorescence, which indicated that retinol induced more structural changes on β-LG at pH 7. β-LG–retinol complex formation at two different pH values was recorded via applying resonance light scattering (RLS) and zeta potential. Conductometry and RLS showed two different behaviors of interaction between β-LG and retinol at two different pH values; therefore, dimmer formation played important roles in different behaviors of interaction between β-LG and retinol. The zeta potential was the implied combination of electrostatic and hydrophobic forces which are involved in β-LG–retinol complex at two different pH values, and the hydrophobic interactions play a dominant role in complex formation. Molecular modeling was approved by all experimental results. The acquired results suggested that monomer and dimmer states of β-LG can be induced by retinol with different behaviors.     

3D-QSAR,molecular docking and molecular dynamics studies of a series of RORγt inhibitors

The discovery of clinically relevant inhibitors of retinoic acid receptor-related orphan receptor-gamma-t (RORγt) for autoimmune diseases therapy has proven to be a challenging task. In the present work, to find out the structural features required for the inhibitory activity, we show for the first time a three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for a series of novel thiazole/thiophene ketone amides with inhibitory activity at the RORγt receptor. The optimum CoMFA and CoMSIA models, derived from ligand-based superimposition I, exhibit leave-one-out cross-validated correlation coefficient (R²_cv) of .859 and .805, respectively. Furthermore, the external predictive abilities of the models were evaluated by a test set, producing the predicted correlation coefficient (R²_pred) of .7317 and .7097, respectively. In addition, molecular docking analysis was applied to explore the binding modes between the inhibitors and the receptor. MD simulation and MM/PBSA method were also employed to study the stability and rationality of the derived conformations, and the binding free energies in detail. The QSAR models and the results of molecular docking, MD simulation, binding free energies corroborate well with each other and further provide insights regarding the development of novel RORγt inhibitors with better activity.     

Preventing the return of smallpox: molecular modeling studies on thymidylate kinase from Variola virus

Smallpox was one of the most devastating diseases in the human history and still represents a serious menace today due to its potential use by bioterrorists. Considering this threat and the non-existence of effective chemotherapy, we propose the enzyme thymidylate kinase from Variola virus (VarTMPK) as a potential target to the drug design against smallpox. We first built a homology model for VarTMPK and performed molecular docking studies on it in order to investigate the interactions with inhibitors of Vaccinia virus TMPK (VacTMPK). Subsequently, molecular dynamics (MD) simulations of these compounds inside VarTMPK and human TMPK (HssTMPK) were carried out in order to select the most promising and selective compounds as leads for the design of potential VarTMPK inhibitors. Results of the docking and MD simulations corroborated to each other, suggesting selectivity towards VarTMPK and, also, a good correlation with the experimental data.     

Memory enhancement by traditional Chinese medicine?

Abstract

Cognitive repair by insulin-like growth factor-I (IGF-I) through activation of insulin-like growth factor-I receptor (IGF-IR) is well established, but not used for clinical therapy due to its link to cancer. We hypothesize that IGF-IR activation rather than IGF-I per se may be essential for cognitive repair and attempted to identify ligands from traditional Chinese medicine (TCM) with drug-like potential towards IGF-IR. TCM ligands, 3-(2-carboxyphenyl)-4(3H)-quinazolinone from Isatisin digotica, (+)-N-methyllaurotetanine from Lindera aggregate, and (+)-1(R)-Coclaurine from Nelumbonucifera Gaertn, exhibited high binding affinities and good blood brain barrier (BBB) penetration crucial for accessing IGF-IR. Stable complex formation of the candidates was observed during molecular dynamics (MD) simulation. Interactions with Leu975 and Gly1055 or Asp1056 were important for ligand binding. Amino acid distance analysis revealed residues 974/975, 984–986, 996–1006, 1040–1056, and 1122–1135 as “hotspots” for ligand binding in IGF-IR. Versatile entry pathways for the TCM candidates suggest high accessibility to the binding site. Blockage of the binding site opening by the TCM candidates limits binding site access by other compounds. Multiple linear regression (R²?=?0.9715), support vector machine (R²?=?0.9084), Bayesian network (R²?=?0.8233) comparative molecular field analysis (CoMFA, R²?=?0.9941), and comparative molecular similarity indices analysis (CoMSIA, R²?=?0.9877) models consistently suggest that the TCM candidates might exert bioactivity on IGF-IR. Contour of representative MD conformations to CoMFA and CoMSIA maps exhibits similar results. Properties including BBB passage, evidence of ability to form stable complexes with IGF-IR by MD simulation, and predicted bioactivity suggest that the TCM candidates have drug-like properties and might have potential as cognitive-enhancing drugs.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:38     

Modeling of Peptides in Implicit Membrane-Mimetic Media

Abstract

Lipid bilayer plays a crucial role in folding of membrane peptides and their stabilization in the membrane-bound state. Correct treatment of the media effects is thus essential for realistic simulations of peptides in bilayers. Previously (Volynsky et al., 1999), we proposed an efficient solvation model which mimics heterogeneous membrane-water system. The model is based on combined employment of atomic solvation parameters for water and hydrocarbon, which approximate hydrated headgroups and acyl chains of lipids, respectively. In this study, the model is employed in non-restrained Monte Carlo simulations of several peptides: totally apolar 20-residue poly-L-Leu, hydrophobic peptide with polar edges, and strongly amphiphilic pep-tide. The principal goals are: to explore energy landscape of these peptides in membrane; to characterize the structures of low-energy states and their orientations with respect to the bilayer. Simulations were performed starting from different structures (unordered or helical) and orientations. It was found that the membrane environment significantly promotes an α-helical conformation for all the peptides, while their energetically favourable orientations are quite different. Thus, poly-Leu was immobilized inside the membrane, the hydrophobic peptide with polar termini adapted transbilayer orientation, whereas the amphiphilic peptide stayed on the lipid-water interface in peripherial orientation. Energy barriers between different states were characterized. The computational results were compared with the experimental structural data.     

9. LIST OF MEMBERS (as at date of publication)

We studied the breeding ecology of the Chestnut-backed Sparrowlark Eremopterix leucotis over three years between 2008 and 2010. The breeding season was bimodal with a main peak in laying in autumn (March–April) and another smaller peak in spring (September–October). Nest microhabitat analyses showed they prefer nesting in open areas with lots of bare ground (median 67.5%). Nest entrance directions were biased towards the south (mean vector (µ) = 186.44°). The majority of nests (78.2%) had an apron at the nest entrance. The mean clutch size was 1.88 but there was geographic variation in clutch size between northern and southern races of the species. The mean incubation and nestling periods were 10.33 d (range 10–11 d) and 9.20 d (range 8–10 d), respectively. The results suggest that parental contributions during incubation are almost equal, but females made significantly more food deliveries during the nestling period compared to males. The diet of nestlings comprised mainly of invertebrates (50.2%), seeds (34.4%) and unidentified food items (15.4%). Breeding success was low, averaging 16.1% (range 8.1–20.6%), and the average number of fledged young per pair was 0.36 ± 0.71. Replacement broods were common and we also recorded repeat brooding attempts.