Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly² and Phe⁴ groups. The three-dimensional model of the δ-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the δ-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible δ-receptor binding conformation. Figure Schematic model of the extended antiparallel dimer structure of Met-enkephalin     

Study of the influence of temperature on the dynamics of the catalytic cleft in 1,3-1,4-β-glucanase by molecular dynamics simulations

The dependence of some molecular motions in the enzyme 1,3-1,4-β-glucanase from Bacillus licheniformis on temperature changes and the role of the calcium ion in them were explored. For this purpose, two molecular dynamics simulated trajectories along 4 ns at low (300 K) and high (325 K) temperatures were generated by the GROMOS96 package. Several structural and thermodynamic parameters were calculated, including entropy values, solvation energies, and essential dynamics (ED). In addition, thermoinactivation experiments to study the influence of the calcium ion and some residues on the activity were conducted. The results showed the release of the calcium ion, which, in turn, significantly affected the movements of loops 1, 2, and 3, as shown by essential dynamics. These movements differ at low and high temperatures and affect dramatically the activity of the enzyme, as observed by thermoinactivation studies. The first two authors contributed equally to this work     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid–liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. Figure Log–log plot of the self-diffusion coefficient D versus the chain length N. The error bars are the standard deviation measured in three repeated simulations     

Activation energies of selenoxide elimination from Se-substituted selenocysteine

Transition states for selenoxide elimination have been determined for a series of Se-substituted selenocysteine (RSeCys) derivatives that have potential use in the prevention and treatment of cancer, either directly or in conjunction with cisplatin (to reduce its nephrotoxic effects). Reduced activation barriers vs R=Me and R=Ph are found when the alkyl chain length is increased or when activating groups are para to the selenide. Ortho substitution of Lewis bases stabilizes the transition state by directly donating electron density to the selenoxide. The results suggest that RSeCys derivatives incorporating the properties of glutathione peroxidase mimics will, upon oxidation, rapidly eliminate selenenic acid, a precursor to chemopreventative selenols. Scheme Mechanism of selenoxide elimination from Se-substituted SeCys.     

CHIH-DFT determination of the molecular structure infrared spectra,UV spectra and chemical reactivity of three antitubercular compounds: Rifampicin,Isoniazid and Pyrazinamide

Three of the most frequent antitubercular agents employed against Mycobacterium tuberculosis are: Rifampicin, Isoniazid and Pyrazinamide. It has been proven that the use of these antitubercular agents together, shortens the treatment period from 12–18 months to 6 months [1]. In this work we use a new Density Functional Theory chemistry model called CHIH-DFT (Chihuahua-Heterocycles-Density Functional Theory) that reflects the mixture of Hartree Fock exchange and DFT exchange, according to a mixing parameter based on empirical rules suited for heterocyclic systems. This new chemistry model was used to calculate the molecular structure of these antitubercular compounds, as well as their infrared, UV spectra, chemical reactivity and electronic properties. The UV and infrared spectra were obtained by experimental techniques. The calculated molecular structure, UV and IR spectra values from CHIH-DFT were compared with experimentally obtained values and theoretical studies. These results are in good agreement with experimental and theoretical studies. We also predicted using the relative electrophilicity and relative nucleophilicity concepts as defined by Roy et al. [2] the chemical active sites for the three antitubercular compounds as well as their electronegativity, ionization potential, electron affinity, hardness, dipole moment, E_HOMO-E_LUMO gap energy, etc.      

SOMMER: self-organising maps for education and research

SOMMER is a publicly available, Java-based toolbox for training and visualizing two- and three-dimensional unsupervised self-organizing maps (SOMs). Various map topologies are implemented for planar rectangular, toroidal, cubic-surface and spherical projections. The software allows for visualization of the training process, which has been shown to be particularly valuable for teaching purposes. Spread of a spherical self-organizing map (SOM) in a three-dimensional data space     

Three-dimensional model of zeaxanthin binding PsbS protein associated with nonphotochemical quenching of excess quanta of light energy absorbed by the photosynthetic apparatus

A three-dimensional model of the PsbS protein was built with the help of homology-modeling methods. This protein is also known as CP22 and is associated with the protection of photosystem II of thylakoid from excess quanta of light energy absorbed by the photosynthetic apparatus. PsbS is reported to bind two molecules of zeaxanthin at low pH (<5.0) and is believed to be essential for rapid nonphotochemical quenching (qE) of chlorophyll a fluorescence in photosystem II. An attempt was made to explain the pH modulation of the conformation of protein through salt-bridges Glu⁻(122)-Lys⁺(113) and Glu⁻(226)-Lys⁺(217). Binding of two molecules of zeaxanthin in the three-dimensional model of PsbS is postulated. The molecular mechanism of photoprotection by PsbS is explained through the model. 1 Backbone structure of the PsbS protein with two molecules of all trans zeaxanthin (ZEX). Residues Glu 90, 122, 194, 226 and Lys 113, 217 are shown. The figure is drawn with RASMOL (Molecular Visualization Program, RasMol V2.6, Roger Sayle, Glaxo Wellcome Research and Development, Stevenage, Hertfordshire, UK) Electronic Supplementary Material Supplementary material is available for this article at     

Structural model of an antistasin/notch-like fusion protein from the cocoon wall of the aquatic leech, Theromyzon tessulatum

The aquatic leech, Theromyzon tessulatum, secretes a proteinaceous cocoon with extraordinary physical properties (e.g., proteolytic, thermal resiliency). The deduced amino acid sequence of a major protein (Tcp—Theromyzon cocoon protein) from the T. tessulatum cocoon wall has been used to model the endogenous structure of the Tcp protein. The Tcp protein sequence comprises six internal repeats, each containing 12 ordered Cys residues. Amino acid alignments suggest that the region Cys1→6 is homologous to antistasin, a leech anticoagulant, and Cys7→12 is homologous to an epidermal growth factor-like domain found in notch-class proteins, which play critical roles in development, signaling, and adhesion throughout the Animalia. Modeling of individual domains (i.e., antistasin and notch) positions multiple hydrophobic and charged residues on the surface. When the antistasin and notch domains were fused, hydrophobic pockets appeared that may facilitate a polymerization mechanism. Collectively, the predicted features of our Tcp model are consistent with the physical properties of the leech cocoon wall.      

Influence of different transition metals in phthalocyanines on their interaction energies with volatile organic compounds: an experimental and computational study

Experimental partition coefficients were determined for a series of volatile organic compounds (VOCs) (acetonitrile, n-butylamine, n-octane tetrachloroethene, and toluene) for the interaction with 2,3,9,10,16,17,23,24-octakis(octyloxy)-phthalocyaninato complexes, PcM(OR)₈, with varying central metal atoms [M=H₂ (metal-free), Ni, Pd, Cu, Zn]. Large partition coefficients for toluene were observed in the case of the nickel and palladium phthalocyanines, whereas for the corresponding zinc-containing compound, interaction with n-butylamine resulted in a high value for the partition coefficient. Interaction energies for model coordination complexes were obtained at the ab initio LMP2/ LACVP* level of theory. The interaction of various small volatiles with the various PcM(OR)₈ compounds was studied using the PM3 semiempirical Hamiltonian. Large values for interaction energies correspond to particularly strong partition coefficients, suggesting that coordination of the volatiles to the central metal dominates over the often discussed π-system stacking at the PcM(OR)₈’s. Figure: Chemical structure of the phthalocyanines investigated.     

Molecular model for host–guest interaction of tetraamino-<Emphasis Type="Italic">tert</Emphasis>-butylthiacalix[4]arene and tetraamino-<Emphasis Type="Italic">tert</Emphasis>-butylcalix[4]arene receptors with carboxylate and dicarboxylate guests: an ONIOM study

Geometry optimizations of tetraamino-tert-butylthiacalix[4]arene (tatbtc4a) and tetraamino-tert-butylcalix[4]arene (tatbc4a) complexes with acetate, oxalate, malonate, succinate, glutarate, adipate, and pimelate were carried out using the integrated MO:MO method. Thermodynamic quantities, preorganization energies and complexation energies of these complexes were obtained at the ONIOM(B3LYP/6-31G(d):AM1) level of theory. The relative stabilities of the tatbtc4a and tatbc4a complexes with carboxylate guests are reported. The complexes tatbtc4a/malonate and tatbc4a/oxalate were found to be the most stable species. The selectivity of the tatbtc4a receptor toward to malonate with respect to oxalate, in terms of selectivity coefficient, is 9.90×10². Figure Atom labeling of tatbtc4a/oxalate complex as a representative of host-guest system.     

Alpha particle chemistry. On the formation of stable complexes between He2+ and other simple species: implications for atmospheric and interstellar chemistry

The possibility that stable complexes may be formed between alpha particles (He²⁺) and small molecules is investigated using QCISD quantum mechanical calculations. Implications for their presence in the terrestrial atmosphere and/or in interstellar space are discussed. Figure Optimized structure of a stable H₂OHe²⁺ complex     

CoMFA study of distamycin analogs binding to the minor-groove of DNA: a unified model for broad-spectrum activity

A 3D-QSAR analysis has been carried out by comparative molecular field analysis (CoMFA) on a series of distamycin analogs that bind to the DNA of drug-resistant bacterial strains MRSA, PRSP and VSEF. The structures of the molecules were derived from the X-ray structure of distamycin bound to DNA and were aligned using the Database alignment method in Sybyl. Statistically significant CoMFA models for each activity were generated. The CoMFA contours throw light on the structure activity relationship (SAR) and help to identify novel features that can be incorporated into the distamycin framework to improve the activity. Common contours have been gleaned from the three models to construct a unified model that explains the steric and electrostatic requirements for antimicrobial activity against the three resistant strains. Figure A unified CoMFA model for broad-spectrum DNA minor-groove binders     

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.      

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

Molecular models of protein kinase 6 from Plasmodium falciparum

Cyclin-dependent kinases (CDKs) have been identified as potential targets for development of drugs, mainly against cancer. These studies generated a vast library of chemical inhibitors of CDKs, and some of these molecules can also inhibit kinases identified in the Plasmodium falciparum genome. Here we describe structural models for Protein Kinase 6 from P. falciparum (PfPK6) complexed with Roscovitine and Olomoucine. These models show clear structural evidence for differences observed in the inhibition, and may help designing inhibitors for PfPK6 generating new potential drugs against malaria. Figure Ribbon diagram of PfPK6 complexed with a roscovitine and b olomoucine     

Experimental and theoretical investigations of spectroscopic properties of azobenzene derivatives in solution

The UV-Vis spectra of series of polymethylmethacrylate (PMMA) copolymers with attached trans-azobenzene derivatives were measured in 1,1,2-trichloroethane. In order to gain some insight into the recorded spectra, the quantum chemical calculations were performed for the substituted azobenzenes using both configuration interaction with single excitations method (CIS) as well as density functional theory (DFT) with B3LYP and PBE0 functionals. The calculations were performed in solvent. In particular, we found that the PBE0 excitation energies are in very good agreement with the experimental data. Figure The plots of orbital contour surfaces for molecule II. The molecular orbitals were calculated at the PBE0/6-311++G(d,p) level of theory. The upper plot presents contour surface of HOMO and the lower presents contour surface of LUMO. Shown are the contour surfaces of orbital amplitude 0.04 (red) and -0.04 (blue)     

Localized and delocalized perfluorosemibullvalenes

On the basis of the experimental Gibbs free-energy barrier of the degenerate Cope arrangement in semibullvalene, B3P86 shows the best agreement, while B3LYP and MP2 underestimate and CCSD(T) overestimates the barrier. The substituent effect proposal by Hoffmann has been verified. In contrast to semibullvalenes with either localized energy-minimum structures or delocalized transition-state structures, perfluorosemibullvalene has both localized and delocalized energy-minimum structures that are very close in energy. Localized and delocalized perfluorosemibullvalenes Electronic supplementary material Supplementary material is available for this article at and is accessible for authorized users.     

Insights into catalytic activity of industrial enzyme Co-nitrile hydratase. Docking studies of nitriles and amides

Nitrile hydratase (NHase) is an enzyme containing non-corrin Co³⁺ in the non-standard active site. NHases from Pseudonocardia thermophila JCM 3095 catalyse hydration of nitriles to corresponding amides. The efficiency of the enzyme is 100 times higher for aliphatic nitriles then aromatic ones. In order to understand better this selectivity dockings of a series of aliphatic and aromatic nitriles and related amides into a model protein based on an X-ray structure were performed. Substantial differences in binding modes were observed, showing better conformational freedom of aliphatic compounds. Distinct interactions with postranslationally modified cysteines present in the active site of the enzyme were observed. Modeling shows that water molecule activated by a metal ion may easily directly attack the docked acrylonitrile to transform this molecule into acryloamide. Thus docking studies provide support for one of the reaction mechanisms discussed in the literature. Figure Crystalographic structure of Pseudonocardia thermophila JCM 3095 nitrile hydratase (a) and the non-standard active site (b)     

Molecular geometry,vibrations and electrode potentials of 2-(4,5-dihydroxy-2-methylphenyl)-2-phenyl-2H-indene-1,3-dione; experimental and theoretical attempts

The electrode potential of 2-(4,5-dihydroxy-2-methylphenyl)-2-phenyl-2H-indene-1,3-dione (DMPID) in acetonitrile has been calculated. The calculations were performed using ab initio molecular orbital calculations (HF), and density functional theory (DFT) with the inclusion of entropic and thermochemical corrections to yield free energies of redox reactions. The electrode potential of DMPID was also obtained experimentally with the aid of an electrochemical technique (cyclic voltammetry). The values for geometric parameters and the vibrational frequencies of DMPID and 2-(6-methyl-3,4-dioxocyclohexa-1,5-dienyl)-2-phenyl-2H-indene-1,3-dione (MDPID) were also computed using the same levels with the basis set of 6-31G(d). The calculated IR spectrum of DMPID used for the assignment of IR frequencies was observed in the experimental FT-IR spectrum and the calculated IR and FT-IR observed spectra of DMPID were compared with correlation factor of 0.996. It should be mentioned that the present work is the first research on coagulant derivative molecules in which the electrode potential of a molecule is calculated. Optimized structures of 2-(6-methyl-3,4-dioxocyclohexa-1,5-dienyl)-2-phenyl-2H-indene-1,3-dione (MDPID)