Predicting anti-HIV activity of 2,3-diaryl-1,3-thiazolidin-4-ones: computational approach using reformed eccentric connectivity index

The eccentric connectivity index, which has recently been employed successfully for the development of numerous mathematical models for the prediction of biological activities of diverse nature, has been reformed to overcome its limitations caused by degeneracy and insensitivity towards heteroatoms. The reformed eccentric connectivity index, termed the eccentric connectivity topochemical index, overcomes the limitations of the eccentric connectivity index by exhibiting very low degeneracy and displaying sensitivity to both the presence and relative position of heteroatoms without compromizing the discriminating power of the eccentric connectivity index. The relationship of the eccentric connectivity topochemical index, eccentric connectivity index and Wieners index with regard to the anti-HIV activity of 2, 3-diaryl-1, 3-thiazolidin-4-one derivatives was subsequently investigated. The values of the eccentric connectivity topochemical index, the eccentric connectivity index and Wieners index of each of 31 analogues comprizing the data set were computed using in-house computer program. Resultant data was analyzed and suitable models developed after identification of active ranges. Subsequently, each derivative was assigned a biological activity using these models, which was then compared with the reported anti-HIV activity. The accuracy of prediction using these models was found to vary from 81 to 90%. The proposed index offers a vast potential for virtual screening of combinatorial libraries, structure property/activity studies and drug design.Figure Basic structure of 2,3-diaryl-1, 3-thiazoidin-4-ones.     

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     

Modeling study of the influences of the aromatic transitions and the local environment on the far-UV rotational strengths in TEM-1 β-lactamase

Rotational strengths in the far-UV of TEM-1 β-lactamase have been investigated with two theoretical models based on the matrix method. The first model excludes, and a second includes, effects of the local electrostatic interactions on the chromophore energies. Special attention is given to the contributions of the aromatic side-chain chromophores, and the mechanisms of generation of rotational strengths are analyzed. The sensitivity of the computational models with respect to the structural changes of the protein are discussed. Figure Structure of TEM-1 β-lactamase. Both domains—α and αβ, the secondary structural elements and the aromatic and disulfide chromophores are shown     

Novel distance-based atom-type topological indices DAI for QSPR/QSAR studies of alcohols

In this work, we propose a distance-based atom-type topological index (DAI) for quantitative structure-property/activity relationship (QSPR/QSAR) studies. The newly constructed index, which codes the structural environment of each atom type in a molecule, can be calculated simply. These atom-type topological indices, along with our recently proposed Lu index, were used to construct QSPR/QSAR models for several representative physical properties and biological activities of several data sets of alcohols with a range of non-hydrogen atoms by using multiple linear regression (MLR) analysis. The efficiency of these indices is verified by high quality QSPR models. The results indicate that the combined use of Lu and DAI indices promises to be a useful method for QSPR/QSAR analysis of complex compounds.      

Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of various benzodiazepine analogues of γ-secretase inhibitors

A 3D QSAR analysis has been performed on a series of 67 benzodiazepine analogues reported as γ-secretase inhibitors using molecular field analysis (MFA), with G/PLS to predict steric and electrostatic molecular field interaction for the activity. The MFA study was carried out using a training set of 54 compounds. The predictive ability of model developed was assessed using a test set of 13 compounds ( as high as 0.729). The analyzed MFA model has demonstrated a good fit, having r² value of 0.858 and cross validated coefficient, value as 0.790. The analysis of the best MFA model provided insight into possible modification of the molecules for better activity.   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)     

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.     

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.      

Topological research on diamagnetic susceptibilities of organic compounds

A novel molecular connectivity index, , based on the adjacency matrix of molecular graphs and novel atomic valence connectivities, , for predicting the molar diamagnetic susceptibilities of organic compounds is proposed. The is defined as: , where and E_i are the atomic valence connectivity and the valence orbital energy of atom i, respectively. A good QSPR model for molar diamagnetic susceptibilities can be constructed from and using multivariate linear regression (MLR). The correlation coefficient r, standard error, and average absolute deviation of the MLR model are 0.9918, 5.56 cgs, and 4.26 cgs, respectively, for the 721 organic compounds tested (training set). Cross-validation using the leave-one-out method demonstrates that the MLR model is highly reliable statistically. Using the MLR model, the average absolute deviations of the predicted values of molar diamagnetic susceptibility of another 360 organic compounds (test set) is 4.34 cgs. The results show that the current method is more effective than literature methods for estimating the molar diamagnetic susceptibility of an organic compound. The MLR method thus provides an acceptable model for the prediction of molar diamagnetic susceptibilities of organic compounds. Figure Plot of calculated vs experimental values of molar diamagnetic susceptibilities using the multivariate linear regression (MLR) model (Eq. 8)     

Theoretical study of the intermolecular H-bonding and intermolecular proton transfer between isocytosine tautomeric forms and <Emphasis Type="Italic">R</Emphasis>,<Emphasis Type="Italic">S</Emphasis>-lactic acid

Eight H-bonded complexes between isocytosine (isoC) tautomeric forms and R/S-lactic acid (LA) have been studied at the B3LYP and HF levels of theory using 6–31+G(d) basis set. The energy barriers of the intermolecular proton transfers were also estimated as the results showed that they are several times lower than those of the intramolecular proton transfers of isoC in the gas phase. Furthermore, the energy barriers of the tautomerizations in which the carboxylic H-atom takes part are several times lower than those in which the LA OH group assists the proton transfer. Figure     

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.      

Homology modeling of a novel epoxide hydrolase (EH) from Aspergillus niger SQ-6: structure-activity relationship in expoxides inhibiting EH activity

The 3D structure of a novel epoxide hydrolase from Aspergillus niger SQ-6 (sqEH) was constructed by using homology modeling and molecular dynamics simulations. Based on the 3D model, Asp191, His369 and Glu343 were predicted as catalytic triad. The putative active pocket is a hydrophobic environment and is rich in some important non—polar residues (Pro318, Trp282, Pro319, Pro317 and Phe242). Using three sets of epoxide inhibitors for docking study, the interaction energies of sqEH with each inhibitor are consistent with their inhibitory effects in previous experiments. Moreover, a critical water molecule which closes to the His369 was identified to be an ideal position for the hydrolysis step of the reaction. Two tyrosine residues (Tyr249 and Tyr312) are able to form hydrogen bonds with the epoxide oxygen atom to maintain the initial binding and positioning of the substrate in the active pocket. These docked complex models can well interpret the substrate specificity of sqEH, which could be relevant for the structural—based design of specific epoxide inhibitors. Figure       

Theoretical study of the hydrogen-bonded complexes serotonin–water/hydrogen peroxide

Eight H-bonded complexes between serotonin (5-hydroxy-tryptamine) and water/hydrogen peroxide were studied at the B3LYP and HF levels of theory, using the 6-31+G(d) basis set. A thermodynamic analysis was performed in order to find the most stable complex. The calculated bonding parameters showed that the most stable H-bonded complex is formed between serotonin and hydrogen peroxide by means of the intermolecular H-bond –H₂N...H–OOH. Fig. a Theoretical study of the hydrogen-bonded supersystems serotonin-water/hydrogen peroxide     

The theoretical investigation of one of the derivatives of 1, 2-dithienylcyclopentene as a molecular switch

The structural and electronic properties of a three-state molecular switch—an active device in a nano-electronic circuit—were studied using the B3LYP/6-31G* method. Due to its chemical stability, high conductivity upon doping, and non-linear optical properties, polythiophene is among the most widely studied conjugated organic polymers, both experimentally and theoretically. The aim of the present work was to theoretically study a very complex case: a three-state switch synthesized and experimentally investigated by Nishida et al. (Org Lett 6:2523–2526, 2004). An initial set of test calculations showed B3LYP level of theory and 6-31G* basis set to be the most appropriate for our purpose, i.e., the study of the structure, charge and spin distributions, as well as electrical characteristics such as electric polarizability, HOMO-LUMO gap (HLG) and electric dipole moment, for one of the 1,2-dithienylcyclopentene derivatives. Also, natural bond orbital analyses were performed to calculate local charges and charge transfers in order to study the capability of the molecule as a molecular switch. The results reported here are of general significance, and demonstrate that it is possible to use certain structural and electrical properties to understand and design electro-photochromic compounds showing a switching function in cases where stable forms can be exchanged by light or electron transfer. Figure Model of a thiophene wire incorporating a redox active unit     

Theoretical study of spin-spin coupling across the hydrogen (O-H⋯N) bond in adenosine derivatives

The study of spin-spin coupling constants across hydrogen bond provides useful information about configuration of complexes. The interesting case of such interactions was observed as a coupling across an intramolecular hydrogen bond in 8-bromo-2′,3′-O-isopropylideneadenosine between the -CH₂OH (at 5″ proton) group and the nitrogen atom of adenine. In this paper we report theoretical investigations on the ^4h J _NH coupling across the H″-C-O-H···N hydrogen bond in adenosine derivatives in various solvent models. Figure Coupling constants in 8-bromo-2′,3′-O-isopropylideneadenosine Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Interactions in pseudorotoxanes based on crown ether-secondary ammonium motifs. A theoretical study

A theoretical analysis of the nature of the interactions in dibenzo[24]crown-8 (DB24C8)-n-dibutylammonium (DBM)—pseudorotaxane complex at the MP2 and DFT levels shows that the main contribution to the binding energy is the electrostatic interaction with moderate (20–25%) correlation stabilization. The total binding energy in the DB24C8-DBM complex represents a sum of the binding energies of two NH–O and one CH–O hydrogen bonds and the latter constitutes about 25% of the total interaction energy, giving the total binding energy of −41.2 kcal mol⁻¹ at the BHandHLYP/6-311++G** level. Deprotonation of the DB24C8-DBM complex reduces the binding energy by some 50 kcal mol⁻¹, giving metastable complexes DB24C8-DBA-1 or DB24C8-DBA-2, which will dissociate to give free crown ether and n-dibutylamine because of the strong exchange repulsion that prevails in neutral complexes. Figure Formation of DB24C8-DBM pseudorotoxane complex