Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative QSAR studies on peptide deformylase inhibitors

Comparative quantitative structure–activity relationship (QSAR) analyses of peptide deformylase (PDF) inhibitors were performed with a series of previously published (British Biotech Pharmaceuticals, Oxford, UK) reverse hydroxamate derivatives having antibacterial activity against Escherichia coli PDF, using 2D and 3D QSAR methods, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and hologram QSAR (HQSAR). Statistically reliable models with good predictive power were generated from all three methods (CoMFA r ² = 0.957, q ² = 0.569; CoMSIA r ² = 0.924, q ² = 0.520; HQSAR r ² = 0.860, q ² = 0.578). The predictive capability of these models was validated by a set of compounds that were not included in the training set. The models based on CoMFA and CoMSIA gave satisfactory predictive r ² values of 0.687 and 0.505, respectively. The model derived from the HQSAR method showed a low predictability of 0.178 for the test set. In this study, 3D prediction models showed better predictive power than 2D models for the test set. This might be because 3D information is more important in the case of datasets containing compounds with similar skeletons. Superimposition of CoMFA contour maps in the active site of the PDF crystal structure showed a meaningful correlation between receptor–ligand binding and biological activity. The final QSAR models, along with information gathered from 3D contour and 2D contribution maps, could be useful for the design of novel active inhibitors of PDF. Figure Superimposition of comparative molecular field analysis (CoMFA) contour plot in the active site of peptide deformylase (PDF)     

Theoretical study of [ XN5 ]− (X=O, S, Se, Te) systems

A series of [XN₅]⁻ (X=O, S, Se, Te) compounds has been examined with ab initio and Density Functional Theory (DFT) methods. The five-membered nitrogen ring series of structures are global minima and may exist or be characterized due to their significant dissociation barriers (29.7–32.7 kcal mol⁻¹). Nucleus-independent chemical shifts (NICS) criteria and the presence of (4n+2) π-electrons confirmed that the five-membered nitrogen ring in their structures exhibits characteristics of aromaticity. Thus, the strong stability of the five-membered nitrogen ring structures may be attributed partially to their aromaticity.      

The platinum-olefin binding energy in series of (PH<Subscript>3</Subscript>)<Subscript>2</Subscript>Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R¹R²C = CR³R⁴)Pt(PH₃)₂ [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents. Figure Representation of (CH₂ = CHR)Pt(PPh₃)₂ and the stability chart     

Computational studies of the binding site of α<Subscript>1A</Subscript>-adrenoceptor antagonists

Aimed at achieving a good understanding of the 3-dimensional structures of human α_1A-adrenoceptor (α_1A-AR), we have successfully developed its homology model based on the crystal structure of β₂-AR. Subsequent structural refinements were performed to mimic the receptor’s natural membrane environment by using molecular mechanics (MM) and molecular dynamics (MD) simulations in the GBSW implicit membrane model. Through molecular docking and further simulations, possible binding modes of subtype-selective α_1A-AR antagonists, Silodosin, RWJ-69736 and (+)SNAP-7915, were examined. Results of the modeling and docking studies are qualitatively consistent with available experimental data from mutagenesis studies. The homology model built should be very useful for designing more potent subtype-selective α_1A-AR antagonists and for guiding further mutagenesis studies. Figure The superposition of β₂-AR crystal structure (gold ribbons) and α_1A-AR homology model (blue ribbons)     

Computer-based design of novel HIV-1 entry inhibitors: neomycin conjugated to arginine peptides at two specific sites

Aminoglycoside–arginine conjugates (AAC and APAC) are multi-target inhibitors of human immunodeficiency virus type-1 (HIV-1). Here, we predict new conjugates of neomycin with two arginine peptide chains binding at specific sites on neomycin [poly-arginine-neomycin-poly-arginine (PA-Neo-PA)]. The rationale for the design of such compounds is to separate two short arginine peptides with neomycin, which may extend the binding region of the CXC chemokine receptor type 4 (CXCR4). We used homology models of CXCR4 and unliganded envelope glycoprotein 120 (HIV-1_IIIB gp120) and docked PA-Neo-PAs and APACs to these using a multistep docking procedure. The results indicate that PA-Neo-PAs spread over two negatively charged patches of CXCR4. PA-Neo-PA–CXCR4 complexes are energetically more favorable than AACs/APAC–CXCR4 complexes. Notably, our CXCR4 model and docking procedure can be applied to predict new compounds that are either inhibitors of gp120–CXCR4 binding without affecting stromal cell-derived factor 1α (SDF-1α) chemotaxis activity, or inhibitors of SDF-1α–CXCR4 binding resulting in an anti-metastasis effect. We also predict that PA-Neo-PAs and APACs can interfere with CD4–gp120 binding in unliganded conformation. Figure The r5-Neo-r5-CXCR4 complex. CXCR4 is shown in CPK representation. The negatively charged residues are shown in red and positively charged residues in blue. The r5-Neo-r5 is shown in stick representation, neomycin core is colored yellow and arginine moieties are colored magenta. Two negatively charged patches separated by neutral and positively charged residues are visible.     

3D–QSAR studies of orvinol analogs as κ-opioid agonists

Orvinols are potent analgesics that target opioid receptors. However, their analgesic mechanism remains unclear and no significant preference for subtype opioid receptor has been achieved. In order to find new orvinols that target the κ-receptor, comparative 3D–QSAR studies were performed on 26 orvinol analogs using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best predictions for the κ-receptor were obtained with the CoMFA standard model (q ²=0.686, r ²=0.947) and CoMSIA model combined steric, electrostatic, hydrophobic, and hydrogen bond donor/acceptor fields (q ²=0.678, r ²=0.914). The models built were further validated by a test set made up of seven compounds, leading to predictive r ² values of 0.672 for CoMFA and 0.593 for CoMSIA. The study could be helpful for designing and prepare new category κ-agonists from orvinols.      

Theoretical study on the thermal decomposition of model compounds for Poly (dialkyl fumarate)

The thermal decomposition of model compounds for poly (dialkyl fumarate) was studied by using ab initio and density functional theory (DFT) calculations. To determine the most favorable reaction pathway of thermal decomposition, geometries, structures, and energies were evaluated for reactants, products, and transition states of the proposed pathways at the HF/6-31G(d) and B3LYP/6-31G(d) levels. Three possible paths (I, II and III) and subsequent reaction paths (IV and V) for the model compounds of poly (dialkyl fumarate) decomposition had been postulated. It has been found that the path (I) has the lowest activation energy 193.8 kJ mol⁻¹ at B3LYP/6-31G(d) level and the path (I) is considered as the main path for the thermal decomposition of model compounds for poly (dialkyl fumarate).      

Design of Benzene-1,2-diamines as selective inducible nitric oxide synthase inhibitors: a combined de novo design and docking analysis

Selective inhibition of inducible nitric oxide synthases (iNOS) has been a challenging problem for researchers pursuing work in finding methods to treat inflammatory disorders, shock, etc. Though many inhibitors have been studied to date, all are associated with selectivity or potency problems. Additionally, most of the reported compounds have several similarities and fewer number of novel structures are being tried. There is an increasing need to design novel molecules for this target. In this work, de novo design using LUDI, combined with docking analysis using FlexX has been employed in an attempt to identify novel scaffolds. Benzene-1,2-diamines were identified which could mimic the interactions of the substrate analogs and other inhibitors. Comparative docking scores in each of the isoforms of nitric oxide synthase were employed to recognize hits for iNOS selectivity. Figure Figure shows the docked poses of the ligand M226 along with that of the reference GW274150. (FlexX analysis)     

Vibrational analysis of a solvated green fluorescent protein chromophore

Resonance Raman (RR) spectra of green fluorescent protein (GFP) model chromophores in solution have been simulated with the CASSCF/MM methodology. Although several reports on vibrational analysis of GFP model chromophores have been recently published, the RR spectra were simulated for the first time in explicit solution with the inclusion of the counterion, as these effects are crucial for unambiguously reproducing the vibrational band assignment in the anionic form of the GFP chromophore. This strategy allows for a one-to-one correspondence of the calculated vibrational modes to the observed RR bands, concerning both the location and intensity pattern. In addition, these simulations were complemented with total energy distribution calculations to aid in the unambiguous assignment of the measured spectra. The current study helps to clarify some of the previous RR bands assignments as well as producing some new assignment for the anionic form of GFP chromophore. The explicit solvent simulations and PCM-based calculations are compared to the measured spectra, and these results demonstrate that explicit solvent simulations provide better agreement with experiment, both in terms of vibrational frequencies and intensity distribution. Figure a Correlation of explicit hydration calculations (CASSCF/6-31G*/MM) for the HBI model chromophore and experimental RR data [21]; slope = 0.982, intercept = 27.210 and regression coefficient = 0.997. b Correlation of implicit PCM calculations (CASSCF/6-31G*) for the HBI model chromophore and experimental RR data [21], slope = 1.017, intercept = −48.838 and regression coefficient = 0.984     

Structures of inclusion complexes of halogenbenzoic acids and α-cyclodextrin based on AM1 calculations

Semi-empirical quantum mechanics calculations using AM1 (Austin Method 1) were carried out for various host-guest combinations of α-cyclodextrin and mono-halogen benzoic acids. The energetically favorable inclusion structures were identified. The AM1 results show that α-cyclodextrin complexes with mono-halogen benzoic acid acids (where the halogen is chlorine, bromide, iodine) as guest compounds are more stable in the “head first” position than in the “tail-first” position for meta and para isomers while ortho mono-halogen benzoic acids complexes with α-cyclodextrin are more stable in “tail-first” position. The calculated structures were found to be in good agreement with those obtained from crystalographic databases.      

Synchronous fluorescence determination and molecular modeling of 5-Aminosalicylic acid (5-ASA) interacted with human serum albumin

In this paper, we proposed a new method for the determination of either human serum albumin (HSA) or 5-Aminosalicylic acid (5-ASA) by synchronous fluorescence spectra and examined the interaction between them using the molecular modeling method under simulative physiological conditions. The optimum conditions of synchronous fluorometric determination of HSA were investigated and the method was successfully applied to the determination of 5-ASA added to serum, urine, and saliva samples. The linear range of the determination of HSA and 5-ASA were 1.60 – 414 μg mL⁻¹ and 0.76 –22.95 μg mL⁻¹, the detection limits were 0.552 μg mL⁻¹ and 0.38 μg mL⁻¹, respectively. In addition, the effect of various common ions on the determination of HSA with 5-ASA was also discussed at room temperature. Figure The salicylic acid moiety is located within the binding pocket. The ring of 5-ASA was inserted in the hydrophobic cavity of site I, and it is important to note that the residue ARG-218 and the trptophan residue of HSA (Trp214) are in close proximity to the ring of 5-ASA suggesting the existence of hydrophobic interaction between them.     

The decomposition of α-aminophosphine oxides to phosphonic acid derivatives (P<Superscript>III</Superscript>)

Aminophosphine oxides and aminophosphonates are, in general, very stable compounds. However, following phosphorus–carbon bond cleavage in aqueous acidic media these compounds sometimes decompose to phosphonic acids derivatives (P^III). Despite some controversy in the literature, careful analysis supported by theoretical studies leads to the conclusion that decomposition to P^III derivatives proceeds via an elimination reaction. Figure The decomposition of α-aminophosphine oxides to phosphonic acid derivatives (P^III)     

3D-QSAR study of microsomal prostaglandin E<Subscript>2</Subscript> synthase(mPGES-1) inhibitors

Microsomal prostaglandin E₂ synthase (mPGES-1) has been identified recently as a novel target for treating pain and inflammation. The aim of this study is to understand the binding affinities of reported inhibitors for mPGES-1 and further to design potential new mPGES-1 inhibitors. 3D-QSAR-CoMFA (comparative molecular field analysis) and CoMSIA (comparative molecular similarity indices analysis) - techniques were employed on a series of indole derivatives that act as selective mPGES-1 inhibitors. The lowest energy conformer of the most active compound obtained from systematic conformational search was used as a template for the alignment of 32 compounds. The models obtained were used to predict the activities of the test set of eight compounds, and the predicted values were in good agreement with the experimental results. The 3D-QSAR models derived from the training set of 24 compounds were all statistically significant (CoMFA; q ² = 0.89, r ² = 0.95, , and CoMSIA; q ² = 0.84, r ² = 0.93, , ). Contour plots generated for the CoMFA and CoMSIA models reveal useful clues for improving the activity of mPGES-1 inhibitors. In particular, substitutions of an electronegative fluorine atom or a bulky hydrophilic phenoxy group at the meta or para positions of the biphenyl rings might improve inhibitory activity. A plausible binding mode between the ligands and mPGES-1 is also proposed.     

Docking studies suggest ligand-specific δ-opioid receptor conformations

An automated docking procedure was used to study binding of a series of δ-selective ligands to three models of the δ-opioid receptor. These models are thought to represent the three ligand-specific receptor conformations. Docking results are in agreement with point mutation studies and suggest that different ligands—agonists and antagonists—may bind to the same binding site under different receptor conformations. Docking to different receptor models (conformations) also suggests that by changing to a receptor-specific conformation, the receptor may open or close different binding sites to other ligands. Figure  Ligands 5 (green) and 6 (orange) in bindingpocket BP1 of the R1 δ-opioid receptor model     

Molecular modeling studies of poly lactic acid initiation mechanisms

The two possible routes to synthesize poly (lactic acid) are polycondensation of the lactic acid and ring opening polymerization (ROP) of the lactide. This work involves molecular modeling of the polymerization initiation mechanisms using different initiators a) H₂SO₄ for polycondensation b) aluminum isopropoxide for coordination-insertion ROP c)methyl triflate for cationic ROP, and d) potassium methoxide for anionic ROP. For molecular modeling of PLA, we have benchmarked our approach using Ryner’s work on ROP of L-lactide using stannous (II) 2-ethylhexanoate (Sn(Oct)₂) and methanol as initiators. Our values of -15.2 kcal mol^-1 and -14.1 kcal mol^-1 for enthalpy changes in the two steps of activated complex formation match with Ryner’s. Geometric and frequency optimizations have been done on Gaussian’03 using B3LYP density functional theory along with the basis sets LANL2DZ for metal atoms and 6–31G* and 6–31G** for non metal atoms. The kinetic rate constant for each mechanism has been calculated using the values of energy of activation, change in enthalpy, Gibbs free energy, entropy and the partition functions from the Gaussian’03 output. Our polycondensation rate constant value of 1.07 × 10^–4 se^-1 compares well with 1.51 × 10^–4 se^-1 as reported by Wang. However, ROP rate constants could not be validated due to lack of experimental data. Figure Cationic Ring Opening Polymerization of L-Lactide

---