Application of the novel molecular alignment method using the Hopfield Neural Network to 3D-QSAR

Recently, we investigated and proposed the novel molecular alignment method with the Hopfield Neural Network (HNN). Molecules are represented by four kinds of chemical properties (hydrophobic group, hydrogen-bonding acceptor, hydrogen-bonding donor, and hydrogen-bonding donor/acceptor), and then those properties between two molecules correspond to each other using HNN. The 12 pairs of enzyme-inhibitors were used for validation, and our method could successfully reproduce the real molecular alignments obtained from X-ray crystallography. In this paper, we apply the molecular alignment method to three-dimensional quantitative structure-activity relationship (3D-QSAR) analysis. The two data sets (human epidermal growth factor receptor-2 inhibitors and cyclooxygenase-2 inhibitors) were investigated to validate our method. As a result, the robust and predictive 3D-QSAR models were successfully obtained in both data sets.     

Recoordination of a metal ion in the cavity of a crown compound: a theoretical study. 1. Conformers of arylazacrown ethers and their complexes with Ca2+ cation

Photoinduced recoordination of Ca²⁺ complexes of the photochromic azacrown ethers is studied by the density functional method. The study included model arylazacrown ethers containing various acceptor groups in the aromatic ring in the para position to the azacrown ether moiety and a real azacrown-containing styryl dye. It is found that both free azacrown ethers and their complexes can adopt two types of conformations: (1) axial conformations, in which the aromatic ring axis passing through the crown ether nitrogen N_cr and the opposite atom of the aromatic ring is perpendicular to the root-mean-square (RMS) plane of the crown ether (least-squares fitted plane for all the crown ether atoms), and (2) equatorial conformations, in which the aromatic ring axis only slightly deflects from the RMS plane of the crown ether. In the equatorial conformers, the metal cation is coordinated only to the O atoms of the azacrown ether cycle, the metal—nitrogen bond is broken, and N_cr is conjugated with the aromatic ring. In the axial conformers, the metal cation is additionally coordinated to N_cr. It is found that the presence of an acceptor group bearing a formal positive charge decreases the relative energy of the equatorial conformer and favors metal—nitrogen bond dissociation, which results in the recoordination of the metal cation. However, a long distance between the charged group and N_cr has the reverse effect. The photoinduced recoordination observed in the alkaline-earth metal complexes of the photochromic azacrown ethers is explained by the transitions between the axial and equatorial conformers facilitated by the charge transfer in the excited state of the complex.     

Hierarchical self-assembly of all-organic photovoltaic devices

Photovoltaic devices built by a hierarchical self-assembly process using hydrogen-bonding terminated self-assembled monolayers (SAMs) on gold and the combination of a hydrogen-bonding barbituric acid appended fullerene and a complementary melamine terminated π-conjugated thiophene-based oligomer are presented. The incorporation of these electron donor (oligomer) and electron acceptor (methanofullerene) assemblies into simple photovoltaic (PV) devices as thin films leads to a 2.5 fold-enhancement in photocurrent compared to analogous systems comprising non-hydrogen-bonding C₆₀-oligomer systems, which is ascribed to higher molecular-level ordering. The modification of the gold electrode surface with self-assembled monolayers bearing hydrogen-bonding molecular recognition endgroups was seen to further enhance the PV response of the corresponding functional supramolecular device. This superposition of two types of self-assembly facilitates the generation of binary supramolecular fullerene-containing architectures. Importantly, all functional materials are accessible in a direct fashion.     

An automated method for predicting the positions of hydrogen-bonding atoms in binding sites

Hydrogen bonds are the most specific, and therefore predictable of the intermolecular interactions involved in ligand–protein binding. Given the structure of a molecule, it is possible to estimate the positions at which complementary hydrogen-bonding atoms could be found. Crystal-survey data are used in the design of a program, HBMAP, that generates a hydrogen-bond map for any given ligand, which contains all the feasible positions at which a complementary atom could be found. On superposition of ligands, the overlapping regions of their maps represent positions of receptor atoms to which each molecule can bind. The certainty of these positions is increased by the incorporation of a larger number and diversity of molecules. In this work, superposition is achieved using the program HBMATCH, which uses simulated annealing to generate the correspondence between points from the hydrogen-bonding maps of the two molecules. Equivalent matches are distinguished on the basis of their steric similarity. The strategy is tested on a number of ligands for which ligand–protein complexes have been solved crystallographically, which allows validation of the techniques. The receptor atom positions of thermolysin are successfully predicted when the correct superposition is obtained.     

Role of solvent in determining conformational preferences of alanine dipeptide in water

Evidence from a variety of spectroscopic probes indicates that (phi, psi) values corresponding to the left-handed polyproline II helix (P(II)) are preferred for short alanine-based peptides in water. On the basis of results from theoretical studies, it is believed that the observed preference is dictated by favorable peptide-solvent interactions, which are realized through formation of optimal hydrogen-bonding water bridges between peptide donor and acceptor groups. In the present study, we address this issue explicitly by analyzing the hydration structure and thermodynamics of 16 low-energy conformers of the alanine dipeptide (N-acetylalanine-N'-methylamide) in liquid water. Monte Carlo simulations in the canonical ensemble were performed under ambient conditions with all-atom OPLS parameters for the alanine dipeptide and the TIP5P model for water. We find that the number of hydrogen-bonded water molecules connecting the peptide group donor and acceptor atoms has no effect on the solvation thermodynamics. Instead, the latter are determined by the work done to fully hydrate the peptide. This work is minimal for conformations that are characterized by a minimal overlap of the primary hydration shells around the peptide donor and acceptor atoms. As a result, peptide-solvent interactions favor "compact" conformations that do not include P(II)-like geometries. Our main conclusion is that the experimentally observed preference for P(II) does not arise due to favorable direct interactions between the peptide and water molecules. Instead, the latter act to unmask underlying conformational preferences that are a consequence of minimizing intrapeptide steric conflicts.     

The 2-naphthol-water2 cluster: two competing types of hydrogen-bonding arrangements

The potential energy surfaces of the S(0) and S(1)(pi(*)<--pi) states of the 2-naphthol(H(2)O)(n), n is an element of {1,2} clusters were explored at the level of coupled cluster (CC2) response theory. In the electronic ground state two different types of hydrogen-bonding networks coexist for n=2, (i) a cyclic one [similar to those of the water trimer and phenol(H(2)O)(2)] where the hydroxy group of the aryl alcohol acts simultaneously as H donor for the first, and as H acceptor for the second water molecule, and (ii) a hydrogen-bonding arrangement where the aromatic pi system is taking over the role as H acceptor. In the S(1) state, on the other hand, the cyclic conformers are unstable. Consequently, the first group of cyclic ground state conformers gives rise to broad unstructured band shapes in the absorption spectrum, whereas the second group of conformers involving the aromatic pi system gives rise to nicely structured band shapes. Based on these results the puzzling absorption spectrum of the n=2 cluster can properly be interpreted.     

Molecular simulations for the conformational assessment of a porphyrin-fullerene dyad in different environments

Conformational space of a porphyrin-fullerene dyad with the donor and acceptor connected by a relatively flexible linker is studied by molecular dynamics simulations in both non-polar and polar solvents, as well as in vacuum. The most probable conformations obtained from the vacuum MD simulations were optimized with semi-empirical (SE) and density functional theory (DFT) methods and the extent of the structural changes is assessed. The computational results indicate the co-existence of different conformers in both polar and nonpolar solvents showing agreement with experimental results. The most probable vacuum conformations at 300 K are similar to the ones at 0 K, while the structures most often observed in the solvents show less compact conformations. Optimization with SE and DFT calculations leads to structures, which represent relatively well the folded conformations in solvent, which validates the electronic structure calculations relevant to describing photoinduced electron-transfer in H2P-O34-C60.     

Conformers of diheteroaryl ketones and thioketones: a quantum chemical study of their properties and fundamental intramolecular energetic effects

The conformational behavior of ten diheteroaryl ketones and thioketones is investigated using various quantum chemical methods. These ketones and thioketones are formed by the disubstitution of formaldehyde and thioformaldehyde with such a heteroaryl group as 2-furanyl, 2-thiophenyl, 2-selenophenyl, 2-pyrrolyl or 1-methyl-2-pyrrolyl. For these compounds, their conformational preference and the energetic ordering of their conformers are determined at the MP2 and B3LYP levels of theory. Energetic barriers resulting from the interconversion between conformations are also estimated. The natural bond orbital (NBO) and interacting quantum atoms (IQA) methods are used to study fundamental intramolecular energetic effects influencing the stability of individual conformers. The results of the two methods indicate the great significance of the effect associated with electron delocalization (in the form of either NBO donor–acceptor interactions or the IQA interatomic exchange–correlation interaction energy) in governing the conformational behavior of the investigated diheteroaryl ketones and thioketones.     

Density functional theory and atoms-in-molecules investigation of intramolecular hydrogen bonding in derivatives of malonaldehyde and implications for resonance-assisted hydrogen bonding

A density functional theory (DFT) and atoms-in-molecules (AIM) analysis has been applied to the intramolecular hydrogen bonding in the enol conformers of malonaldehyde and its fluoro-, chloro-, cyano-, and nitro-substituted derivatives. With the B3LYP/6-311++G(2d,p) method, good agreement between the DFT geometries and published experimental structures has been found. The donor-acceptor distance was also varied in a series of constrained optimizations in order to determine if energetic, structural, and topological trends associated with intermolecular hydrogen bonding remain valid in the intramolecular case. At very short donor-acceptor distances (<2.24 A), the hydrogen is symmetrically located between donor and acceptor; at distances longer than this, the hydrogen bonding is no longer symmetric. The AIM methodology has been applied to explore the topology of the electron density in the intramolecular hydrogen bonds of the chosen model systems. Most AIM properties for intramolecular hydrogen bond distances longer than 2.24 A show smooth trends, consistent with intermolecular hydrogen bonds. Integrated AIM properties have also been used to explore the phenomenon of resonance-assisted hydrogen bonding (RAHB). It is shown that as the donor-acceptor distance is varied, pi-electron density is redistributed among the carbon atoms in the intramolecular hydrogen bond ring; however, contrary to prior studies, the integrated atomic charges on the donor-acceptor atoms were found to be insensitive to variation of hydrogen-bonding distance.     

QSAR of conformationally flexible molecules: Comparative Molecular Field Analysis of protein-tyrosine kinase inhibitors

Summary Comparative Molecular Field Analysis (CoMFA) has been applied to a study of quantitative structureactivity relationships (QSAR) of conformationally flexible molecules. The relationship between three-dimensional structure and activity of 20 styrene derivatives which inhibit protein-tyrosine kinase was determined. A technique was developed that allows accurate prediction of the inhibitory activity of these molecules and identification in each case of the active conformation. The problem of multiple energetically acceptable conformations was approached in an iterative procedure. Use was made of the varying degrees of symmetry among the molecules. First, CoMFA QSAR models were developed using only those compounds that possess a symmetrical substituent pattern on the phenyl ring. These CoMFA models were then used to select the active conformers of the less symmetrical compounds in the set. Allowing multiple conformers for each compound in the dataset yielded higher crossvalidated r² values and better predictivity of the QSAR models. Different probe atoms (C⁺, O^–, neutral C) were explored, the O^– probe atom exhibiting the highest selectivity in the conformer selection process.     

Long‐range charge transfer in donor‐peptide bridge‐acceptor model systems—A theoretical study

The quantum mechanical calculations were performed to study the effect of geometrical fluctuations of peptide on charge transfer in model oligopeptides linked between donor and acceptor molecules. The charge transfer parameters have been calculated based on the density functional theory method. Results show that the overall charge transfer in peptide mediated donor–acceptor complexes is determined by the conformations and chain length of the intermediate peptide bridge. The analysis of excess charge distribution show that the localization of an excess positive and negative charge are strongly depend on the conformations and chain length of the donor–bridge‐acceptor system. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Difference spatial distribution function analysis of aqueous solutions. IV. Hydration structure changes of ethylene glycol solutions throughout conformational change process from gGg' to tGg' conformers

Monte Carlo (MC) simulations were carried out for an infinitely dilute aqueous solution of two stable conformers (gGg' and tGg') and of three conformations between gGg' and tGg' conformers of ethylene glycol (EG) at 298K. Based on the spatial distribution function (SDF) goo(x,y,z), obtained from the MC simulation in the above conformations in liquid water, the high distribution of hydration water molecules could be divided into hydrogen acceptor (HA), hydrogen donor (HD), MIX (overlapped distribution of HA and HD), and hydrophobic hydration (HH) regions. The spatial orientations of hydrogen-bonded water molecules were found to be of a linear type with a triple-layer structure in the HA region and HA part (in the MIX region), and double-layer structures in the HD region and HD part (in the MIX region). In addition, it was apparent that the spatial orientations of these water molecules were of the linear type throughout the conformational change process from gGg' to tGg' conformers in liquid water. From the difference SDF (DSDF), deltagoo(x,y, z), between the SDFs of two conformations, we concluded that the distribution of hydration water molecules in the HA and HD parts of the MIX region are governed by the competition of internal hydrogen bonds between the hydrogen atom and two lone-pair electrons on the oxygen atom of an EG molecule.     

[4+3] Cycloaddition of Donor–Acceptor Cyclopropanes with Amphiphilic Benzodithioloimine as Surrogate for ortho‐Bisthioquinone

Donor–acceptor cyclopropanes were reacted with amphiphilic benzodithioloimine to give seven‐membered heterocycles with two sulfur atoms. Formally, this transformation can be regarded as a [4+3] cycloaddition reaction of the three‐membered ring and ortho‐bisthioquinone. The benzodithioloimine serves as a surrogate for this highly reactive diene. The structure of the products was confirmed by X‐ray crystallography. Broad signals in ¹³C NMR studies suggest that several conformers, slowly interconverting on the NMR timescale, are present at room temperature.     

Exploring the dynamics of calix[4]pyrrole: effect of solvent and fluorine substitution

Molecular dynamics simulations show that calix[4]pyrrole (CP) and octafluorocalix[4]pyrrole (8F-CP) are extremely flexible molecules. CP mainly adopts the 1,3-alternate conformation in all the solvents, although the percentage of alternative conformations increases in polar solvents, especially those with good hydrogen-bonding acceptor properties. However, in the case of 8F-CP, the cone conformation is the most populated in some solvents. Transitions between conformers are common and fast, and both CP and 8F-CP can adopt the cone conformation needed for optimum interaction with anions more easily than would be predicted on the basis of previous gas-phase calculations. Furthermore, the present studies show that when a fluoride anion is specifically placed initially in close proximity to CP and 8F-CP in their respective 1,3-alternate conformations, an extremely fast change to the cone conformation is observed in both cases. The results suggest that preorganization does not represent a major impediment to anion-binding for either CP or 8F-CP, and that ion-induced conformational changes can follow different mechanisms depending on the solvent and the chemical substituents present on the calix[4]pyrrole beta-pyrrolic positions.     

The presentation of an approach for estimating the intramolecular hydrogen bond strength in conformational study of β-Aminoacrolein

All the plausible conformations of β-aminoacrolein (AMAC) have been investigated by the Bekes-Lee-Yang-Parr (B3LYP) nonlocal density functional with extended 6-311++G** basis set for studying the stability order of conformers and the various possibilities of intramolecular hydrogen bonding formation. In general the ketoamine (KA) conformers of AMAC, by mean average, are more stable than the corresponding enolimine (EI) and ketoimine (KI) analogues and this stability is mainly due to the π-electron resonance in these conformers that established by NH2 functional group. The contribution of resonance to the stability of chelated KA conformers is about 75.6 kJ/mol, which is greater than that of the hydrogen bond energy (E_HB=35.0 kJ/mol). The relative decreasing order of the various hydrogen bond energies was found to be: O–HN_imine(strong)>N_amine-HO_keto (normal)>N_imine-HO_hydroxyl (weak) > N_imine-HO_keto (weak). Hydrogen bond energies for all systems were obtained from the method that we called related rotamers method (RRM). The topological properties of the electron density contributions for various type of intramolecular hydrogen bond have been analyzed in term of the Bader theory of atoms in molecules (AIM). The results of these calculations support the previous calculations, which obtained by the related rotamer methods.     

Conformation of 1,3,2,6-dioxaphosphazocines according to data obtained by molecular mechanics method

The molecular mechanics method has been used to determine the equilibrium geometric parameters of 39 conformations of the 1,3,2,6-dioxaphosphazocine molecule, among which the preferred form is the boat-boat (BB) conformer with C and O atoms at the vertices. The role of the anomeric effect in stabilizing the axial orientation of the P-S bond has been evaluated quantitatively in various conformers on the 1,3,2,6-dioxaphosphazocine molecule.Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 10, pp. 2283–2288, October, 1989.The authors wish to express their appreciation to R. R. Shagidullin for continuing interest in this work.     

Spectroscopic,X-ray Diffraction and Density Functional Theory Study of Intra- and Intermolecular Hydrogen Bonds in Ortho-(4-tolylsulfonamido)benzamides

The conformations of the title compounds were determined in solution (NMR and UV-Vis spectroscopy) and in the solid state (FT-IR and XRD), complemented with density functional theory (DFT) in the gas phase. The nonequivalence of the amide protons of these compounds due to the hindered rotation of the C(O)–NH₂ single bond resulted in two distinct resonances of different chemical shift values in the aromatic region of their ¹H-NMR spectra. Intramolecular hydrogen bonding interactions between the carbonyl oxygen and the sulfonamide hydrogen atom were observed in the solution phase and solid state. XRD confirmed the ability of the amide moiety of this class of compounds to function as a hydrogen bond acceptor to form a six-membered hydrogen bonded ring and a donor simultaneously to form intermolecular hydrogen bonded complexes of the type N–H···O=S. The distorted tetrahedral geometry of the sulfur atom resulted in a deviation of the sulfonamide moiety from co-planarity of the anthranilamide scaffold, and this geometry enabled oxygen atoms to form hydrogen bonds in higher dimensions.     

Crystal structure of the first SH-containing tetrahedral cobalt(II) complex, [Co(quinoline)2(SH)2]. Superoxide dismutase activity

The preparation, spectroscopic properties and crystal structure of the bis-quinoline–bis-mercaptocobalt(II) complex [Co(quinoline)₂(SH)₂] are reported. The complex is tetrahedral with nitrogen donor atoms from two quinoline ligands and sulphur donor atoms from two mercapto groups. The superoxide dismutase mimetic activity of the complex was investigated using the indirect xanthine–xanthine oxidase–nitroblue tetrazolium method and compared to that of the native enzyme.     

Quantum-chemical investigation of planar form of diethenediol ester of orthosilicic acid

Conclusions From this analysis of the electronic properties of a number of conformations of the Si (DOPh)₂ molecule and dimers of the model system Si (DOE)₂, it follows that isolated molecules of these esters have the form of a tetrahedron. The appreciable population of the p_z and d_z2 orbitals in these forms may tend to give a pseudo-octahedral configuration in solution as a result of coordination of solvent molecules, and this same configuration in the crystalline state as a result of formation of Si–Si bonds. The most stable structure should be that with a fourth-order screw axis, similar to the structure of the salts of Krogmann [6]. Such systems are of considerable interest in connection with the possibility of using them to form quasi-unidimensional conducting systems, this being achieved by doping with the appropriate donor or acceptor additives.Translated from Zhurnal Strukturnoi Khimii, Vol. 25, No. 3, pp. 29–33, May–June, 1984.     

A Method for Multiconformational Modeling of the Three‐Dimensional Shape of a Molecule

A method for multiconformational modeling of the threedimensional shape of a molecule is proposed that includes search for conformers, their optimum superposition, and analysis of spatial features of the resulting structure. The method allows one to determine features of various molecular conformations of compounds under study, to assess the contributions of conformers to particular properties of the substance, to evaluate the space occupied by the molecule, and to compare the average size of the multiconformational model of the molecule with the sizes of the most stable conformations. The potentials of the model are illustrated by density calculations for 137 organic liquids.