Comparison of end-point continuum-solvation methods for the calculation of protein-ligand binding free energies

We have compared the predictions of ligand‐binding affinities from several methods based on end‐point molecular dynamics simulations and continuum solvation, that is, methods related to MM/PBSA (molecular mechanics combined with Poisson–Boltzmann and surface area solvation). Two continuum‐solvation models were considered, viz., the Poisson–Boltzmann (PB) and generalised Born (GB) approaches. The nonelectrostatic energies were also obtained in two different ways, viz., either from the sum of the bonded, van der Waals, nonpolar solvation energies, and entropy terms (as in MM/PBSA), or from the scaled protein–ligand van der Waals interaction energy (as in the linear interaction energy approach, LIE). Three different approaches to calculate electrostatic energies were tested, viz., the sum of electrostatic interaction energies and polar solvation energies, obtained either from a single simulation of the complex or from three independent simulations of the complex, the free protein, and the free ligand, or the linear‐response approximation (LRA). Moreover, we investigated the effect of scaling the electrostatic interactions by an effective internal dielectric constant of the protein (?_int). All these methods were tested on the binding of seven biotin analogues to avidin and nine 3‐amidinobenzyl‐1H‐indole‐2‐carboxamide inhibitors to factor Xa. For avidin, the best results were obtained with a combination of the LIE nonelectrostatic energies with the MM+GB electrostatic energies from a single simulation, using ?_int = 4. For fXa, standard MM/GBSA, based on one simulation and using ?_int = 4–10 gave the best result. The optimum internal dielectric constant seems to be slightly higher with PB than with GB solvation. © Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Cytokinins: Synthesis of 2-, 8-, and 2,8-substituted 6-(3-methyl-2-butenylamino)purines and their relative activities in promoting cell growth

We have synthesized 2- and 8-monosubstituted and 2,8-disubstituted derivatives of the cytokinin 6-(3-methyl-2-butenylamino)purine (N⁶-isopentenyladenine) and have shown the dependence of growth-promoting activity in the tobacco bioassay upon the position, number, and type of substituent. The representative substituent groups were MeS, Me, MeSO₂, C₆H₅CH₂S, HS and Cl. The 8-methyl derivative was exceptional in being more active than the unsubstituted parent compound. In general, substitution in the 8-position decreases activity less than substitution in the 2-position, with the exception of the electron-attracting methylsulfonyl. Substitution in both the 2- and 8-positions lowers the activity more than substitution at either single position on the adenine nucleus, with the exception of the 2,8-dimethyl derivative. The chloro and methylthio derivatives show activity in the same range as the methyl derivatives, and the mercapto compounds, which exist mainly as CS tautomers, show somewhat less activity than the corresponding methylthio compounds. Bulky (C₆H₅CH₂S and MeSO₂) and strongly electron-attracting (MeSO₂) substituents cause relatively great reduction in cytokinin activity.     

Design,synthesis and antifungal evaluation of borrelidin derivatives

Borrelidin, a nitrile containing 18-membered polyketide macrolide, display potent antifungal activity. In this study, a library of borrelidin derivatives were synthesized. Their structures were elucidated by detailed spectroscopic data analysis. The antifungal activity and cytotoxicity of these target compounds were evaluated by broth microdilution and 3-(4,5-dimethylthiazol-2-yl)-3,5-phenytetrazoliumromide (MTT) methods. Among forty-seven prepared analogues, compound 3b had the inhibitory effect on Candida albicans and Candida parapsilosis (MIC: 50 and 12.5?μg/mL, respectively). Furthermore, compounds 4n and 4r presented better antifungal activity against Aspergillus fumigatus with 12.5?μg/mL MIC value, which were insensitive to borrelidin. Preliminary structure-activity relationships (SAR) revealed that the ester analogues containing fragment -OCH₂CH₂N- had an important effect on the antifungal activity. Meanwhile, the molecular docking study indicated the carboxyl substituents in BN could provide extra interaction with pathogenic fungal threonyl-tRNA synthetase (ThrRS).     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Design,synthesis, and evaluation of curcumin analogues as potential inhibitors of bacterial sialidase

Sialidases are key virulence factors that remove sialic acid from the host cell surface glycan, unmasking receptors that facilitate bacterial adherence and colonisation. In this study, we developed potential agents for treating bacterial infections caused by Streptococcus pneumoniae Nan A that inhibit bacterial sialidase using Turmeric and curcumin analogues. Design, synthesis, and structure analysis relationship (SAR) studies have been also described. Evaluation of the synthesised derivatives demonstrated that compound 5e was the most potent inhibitor of S. pneumoniae sialidase (IC₅₀?=?0.2?±?0.1?µM). This compound exhibited a 3.0-fold improvement in inhibitory activity over that of curcumin and displayed competitive inhibition. These results warrant further studies confirming the antipneumococcal activity 5e and indicated that curcumin derivatives could be potentially used to treat sepsis by bacterial infections.     

Design,synthesis and evaluation against Mycobacterium tuberculosis of azole piperazine derivatives as dicyclotyrosine (cYY) mimics

Three series of azole piperazine derivatives that mimic dicyclotyrosine (cYY), the natural substrate of the essential Mycobacterium tuberculosis cytochrome P450 CYP121A1, were prepared and evaluated for binding affinity and inhibitory activity (MIC) against M. tuberculosis. Series A replaces one phenol group of cYY with a C3-imidazole moiety, series B includes a keto group on the hydrocarbon chain preceding the series A imidazole, whilst series C explores replacing the keto group of the piperidone ring of cYY with a CH₂-imidazole or CH₂-triazole moiety to enhance binding interaction with the heme of CYP121A1. The series displayed moderate to weak type II binding affinity for CYP121A1, with the exception of series B 10a, which displayed mixed type I binding. Of the three series, series C imidazole derivatives showed the best, although modest, inhibitory activity against M. tuberculosis (17d MIC?=?12.5?μg/mL, 17a 50?μg/mL). Crystal structures were determined for CYP121A1 bound to series A compounds 6a and 6b that show the imidazole groups positioned directly above the haem iron with binding between the haem iron and imidazole nitrogen of both compounds at a distance of 2.2?Å. A model generated from a 1.5?Å crystal structure of CYP121A1 in complex with compound 10a showed different binding modes in agreement with the heterogeneous binding observed. Although the crystal structures of 6a and 6b would indicate binding with CYP121A1, the binding assays themselves did not allow confirmation of CYP121A1 as the target.     

Carbon nanotube functionalization with carboxylic derivatives: a DFT study

Chemical functionalization of a single-walled carbon nanotube (CNT) with different carboxylic derivatives including –COOX (X?=?H, CH₃, CH₂NH₂, CH₃Ph, CH₂NO₂, and CH₂CN) has been theoretically investigated in terms of geometric, energetic, and electronic properties. Reaction energies have been calculated to be in the range of ?0.23 to ?7.07 eV. The results reveal that the reaction energy is increased by increasing the electron withdrawing character of the functional groups so that the relative magnitude order is ?CH₂NO₂?>?CH₂CN?>?H?>?CH₂Ph?>?CH₃?>?CH₂NH₂. The chemical functionalization leads to an increase in HOMO/LUMO energy gap of CNT by about 0.32 to 0.35 eV (except for ?H). LUMO, HOMO, and Fermi level of the CNT are shifted to lower energies especially in the case of ?CH₂NO₂ and _?CH₂CN functional groups. Therefore, it leads to an increment in work function of the tube, impeding the field electron emission.     

Preparation of 4'-benzenesulfonyl-3'-deoxythymidine and its reaction with organoaluminum and organosilicon reagents

The 4′-benzenesulfonyl derivative of 3′-deoxythymidine was prepared from 3′-deoxythymidine-5′-aldehyde. The 4′-benzenesulfonyl leaving group undergoes a nucleophilic substitution with organoaluminum and organosilicon reagents to furnish a variety of 4′-substituted (Me, Et, i-Bu, trimethylsilylethynyl, CH₂CHCH₂, CN, N₃) analogues.     

A theoretical study on the coordination behavior of some phosphoryl,carbonyl and sulfoxide derivatives in lanthanide complexation

The selectivity of phosphoryl P(O)R₃, sulfoxide S(O)R₂, and carbonyl C(O)R₂ (R?=?NH₂, CH₃, OH, and F) derivatives with lanthanide cations (La³⁺, Eu³⁺, Lu³⁺) was studied by density functional theory calculations. Theoretical approaches were also used to investigate energy and the nature of metal–ligand interaction in the model complexes. Atoms in molecules and natural bond orbital (NBO) analyses were accomplished to understand the electronic structure of ligands, L, and the related complexes, L–Ln³⁺. NBO analysis demonstrated that the negative charge on phosphoryl, carbonyl, and sulfoxide oxygen (O_P, O_C, and O_S) has maximum and minimum values when the connected –R groups are –NH₂ and –F. The metal–ligand distance declines as, –F?>?–OH?>?–CH₃?>?–NH₂. Charge density at the bond critical point and on the lanthanide cation in the L–Ln³⁺ complexes varies in the order –F?<?–OH?<?–CH₃?<?–NH₂, due to greater ligand to metal charge transfer, which is well explained by energy decomposition analysis. It was also illustrated that E⁽²⁾ values of Lp(N)?→?σ*(Y–N) vary in the order P=O ? S=O ? C=O and the related values of Lp(N)?→?σ*(Y=O) change as C=O ? S=O ? P=O in (NH₂)_nYO ligands (Y?=?P, C, and S). Trends in the L–Ln³⁺ CP–corrected bond energies are in good accordance with the optimized O_Y?Ln distances. It seems that, comparing the three types of ligands studied, NH₂–substituted are the better coordination ligands.

Open image in new window

Graphical Abstract Density functional theory (B3LYP) calculations were used to compare structural, electronic and energy aspects of lanthanide (La, Eu, Lu) complexes of phosphine derivatives with those of carbonyls and sulfoxides in which the R– groups connected to the P=O, C=O and S=O are –NH₂, –CH₃, –OH and –F.

     

Synthesis,thymidine phosphorylase,angiogenic inhibition and molecular docking study of isoquinoline derivatives

Isoquinoline analogues (KA-1 to 16) have been synthesized and evaluated for their E. coli thymidine phosphorylase inhibitory activity. Except compound 11, all other analogs showed outstanding thymidine inhibitory potential ranging in between 4.40 ± 0.20 to 69.30 ± 1.80 µM when compared with standard drug 7-Deazaxanthine (IC₅₀ = 38.68 ± 4.42 µM). Structure Activity Relationships has been established for all compounds, mainly based on substitution pattern on phenyl ring. All analogs were characterized by various spectroscopic techniques such as ¹H NMR, ¹³C NMR and EI-MS. The binding interactions of isoquinoline analogues with the active site of TP enzyme, the molecular docking studies were performed. Furthermore, the angiogenic inhibitory potentials of isoquinoline analogues (KA-1-9, 14, 12 and 16) were determined in the presence of standard drug Dexamethasone based on percentage inhibitions at various concentrations. Herein this work analogue KA-12, 14 and 16 emerged with most potent angiogenic inhibitory potentials among the synthesized analogues.     

Synthesis and PGE2 inhibitory activity of novel diarylheptanoids

Prostaglandin E₂ (PGE₂) is a lipid mediator of inflammation and its inhibition has become a popular drug target due to its harmful physiological roles. Diarylheptanoids are one class of compounds that have shown successful inhibition of PGE₂. This paper reports the synthesis and PGE₂ inhibitory activity of a series of analogues of a naturally occurring diarylheptanoid. The most efficacious compounds were examined for dose-dependent PGE₂ inhibition. Among several promising compounds, the lead candidate exhibited an IC₅₀ value of 0.56?ng/µL or 1.7?µM with no detectable toxicity at the highest dose of 10?ng/µL.     

Evidence for regulation of Aerobacter aerogenes pH 6 acetolactate-forming enzyme by acetate ion.

Experiments have been carried out with acetic acid and monodeuteroacetic acid (CH₃COOD) at pH 5.5 in H₂O and D₂O respectively, and with acetate analogues at pH 6.5 with substitution in the methyl group. The results indicate that it is the CH₃COO^? which is the activator for the pH 6 acetolactate-forming enzyme.     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

Design and synthesis of pyrazolopyridine derivatives as sphingosine 1-phosphate receptor 2 ligands

Eleven new sphingosine 1-phosphate receptor 2 (S1PR2) ligands were synthesized by modifying lead compound N-(2,6-dichloropyridin-4-yl)-2-(4-isopropyl-1,3-dimethyl-1H-pyrazolo[3,4-b]pyridin-6-yl)hydrazine-1-carboxamide (JTE-013) and their binding affinities toward S1PRs were determined in vitro using [³²P]S1P and cell membranes expressing recombinant human S1PRs. Among these ligands, 35a (IC₅₀?=?29.1?±?2.6?nM) and 35b (IC₅₀?=?56.5?±?4.0?nM) exhibit binding potency toward S1PR2 comparable to JTE-013 (IC₅₀?=?58.4?±?7.4?nM) with good selectivity for S1PR2 over the other S1PRs (IC₅₀?>?1000?nM). Further optimization of these analogues may identify additional and more potent and selective compounds targeting S1PR2.     

Design and synthesis of cholecystokinin-4 dipeptide analogues with anxiolytic and anxiogenic activities

A new series of dipeptide analogues of the general formula Ph(CH₂) _n CO-NH(CH₂) _m CO-Trp-NH₂ (n = 1, 3–5; m = 1–3) was designed based on the structure of the endogenous tetrapeptide cholescystokinin-4 (CCK-4) and the topochemical Shemyakin-Ovchinnikov-Ivanov principle. The L-tryptophan derivatives exhibited anxiolytic properties and the D-tryptophan derivatives, anxiogenic properties. The dipeptide Ph(CH₂)₅CO-Gly-L-Trp-NH₂ (GB-115) with the activity in rats of 0.05–0.2 mg/kg after oral and intraperitoneal administration was chosen for further studies as a promising anxiolytic agent.     

Design,synthesis, cholinesterase inhibition and molecular modelling study of novel tacrine hybrids with carbohydrate derivatives

A series of hybrids containing tacrine linked to carbohydrate-based moieties, such as d-xylose, d-ribose, and d-galactose derivatives, were synthesized by the nucleophilic substitution between 9-aminoalkylamino-1,2,3,4-tetrahydroacridines and the corresponding sugar-based tosylates. All compounds were found to be potent inhibitors of both acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) in the nanomolar IC₅₀ scale. Most of the d-xylose derivatives (6a-e) were selective for AChE and the compound 6e (IC₅₀?=?2.2?nM for AChE and 4.93?nM for BuChE) was the most active compound for both enzymes. The d-galactose derivative 8a was the most selective for AChE exhibiting an IC₅₀ ratio of 7.6 for AChE over BuChE. Only two compounds showed a preference for BuChE, namely 7a (d-ribose derivative) and 6b (d-xylose derivative). Molecular docking studies indicated that the inhibitors are capable of interacting with the entire binding cavity and the main contribution of the linker is to enable the most favorable positioning of the two moieties with CAS, PAS, and hydrophobic pocket to provide optimal interactions with the binding cavity. This finding is reinforced by the fact that there is no linear correlation between the linker size and the observed binding affinities. The majority of the new hybrids synthesized in this work do not violate the Lipinski's rule-of-five according to FAF-Drugs4, and do not demonstrated predicted hepatotoxicity according ProTox-II.     

Insight into the structural requirements of thiophene-3-carbonitriles-based MurF inhibitors by 3D-QSAR,molecular docking and molecular dynamics study

The discovery of clinically relevant inhibitors against MurF enzyme has proven to be a challenging task. In order to get further insight into the structural features required for the MurF inhibitory activity, we performed pharmacophore and atom-based three-dimensional quantitative structure–activity relationship studies for novel thiophene-3-carbonitriles based MurF inhibitors. The five-feature pharmacophore model was generated using 48 inhibitors having IC₅₀ values ranging from 0.18 to 663?μm. The best-fitted model showed a higher coefficient of determination (R²?=?0.978), cross-validation coefficient (Q²?=?0.8835) and Pearson coefficient (0.9406) at four component partial least-squares factor. The model was validated with external data set and enrichment study. The effectiveness of the docking protocol was validated by docking the co-crystallized ligand into the catalytic pocket of MurF enzyme. Further, binding free energy calculated by the molecular mechanics generalized Born surface area approach showed that van der Waals and non-polar solvation energy terms are the main contributors to ligand binding in the active site of MurF enzyme. A 10-ns molecular dynamic simulation was performed to confirm the stability of the 3ZM6-ligand complex. Four new molecules are also designed as potent MurF inhibitors. These results provide insights regarding the development of novel MurF inhibitors with better binding affinity.     

The interaction of poly(N5-(3-hydroxypropyl)-L-glutamine) with solvent components in water/dioxane mixtures.

The preferential binding of solvent components with a nonionic homopolypeptide, poly(N⁵-(3-hydroxypropyl)-L -glutamine), ([Gln((CH₂)₃OH)]_n), has been determined in water/dioxane mixtures using differential refractometry. The degree of preferential binding was calculated from the difference between the refractive index increments of [Gln((CH₂)₃OH)]_n obtained from experiments carried out under two conditions: experiments where the molality of dioxane was kept identical in both compartments of the differential cell, and experiments where the chemical potential was kept identical. The polypeptide was preferentially hydrated between 10 and 70 wt % of dioxane; the amount of preferential hydration per gram of the mixed solvent increases monotonically (with a plateau region between 40 and 60 wt %) with the dioxane concentration. A monotonic increase was also observed in the degree of helicity of the polypeptide. The absolute amounts of water and dioxane bound by [Gln((CH₂)₃OH)]_n were investigated in the frozen state by the method of nuclear magnetic resonance. Hydration was measured using a mixed solvent, water/dioxane-d₈; dioxane solvation was measured using a mixed solvent, dioxane/D₂O. The polypeptide binds about 0.35 g of water per g of the polymer in aqueous solution, and hydration decreases gradually with an increase in dioxane concentration. On the other hand, the amount of dioxane solvation increases to 0.04 g per g of the polymer in the dioxane concentration range between 0 and 20 wt %, and then levels off. The rapid increase in solvation is observed before the conformational transition from random coil to α-helix occurs in [Gln((CH₂)₃OH)]_n. The dependence of the preferential and absolute binding of solvent components to [Gln((CH₂)₃OH)]_n on dioxane concentration and the conformational change in the homopolypeptide suggest that addition of dioxane to aqueous solutions induces lowering of water activity and that the helical structure of the polypeptide is enhanced by the formation of intrachain hydrogen bonds. The validity of the frozen method is also discussed.