Aqueous Extracts of Organic Mulch Materials Have Nematicide and Repellent Effect on Meloidogyne incognita Infective Juveniles: A Laboratory Study.

While the nematicidal effectiveness of mulching against root-knot nematodes (Meloidogyne spp.) is calculated within organic crop protection, underlying mechanisms are not yet fully explored. Laboratory experiments were set up to determine whether mulch-derived substances cause mortality directly, or repel Meloidogyne juveniles from crop rhizosphere. Mortality and area choice tests were conducted with mulch-derived extracts, supported by the measurements on tannic acid content and the pH values of extracts as supplementary examinations. In our study, leaf litter and straw extracts were generally found lethal to the juveniles, which is in line with the results from area preference tests. However, compost extract had no effect on Meloidogyne incognita juveniles. Tannic acid content showed positive correlation with mortality only in the case of straw and sycamore leaf litter extracts. Tannic acid and pH weakly correlated with repellent effect of the applied extracts generally. Our results have inspired further experiments to explore nematicidal components of leaf litters, contributing to the development of a new approach in crop protection based on the repellent effect of these materials.

Replacement of Oxygen by Sulfur in Small Organic Molecules. 3. Theoretical Studies on the Tautomeric Equilibria of the 2OH and 4OH-Substituted Oxazole and Thiazole and the 3OH and 4OH-Substituted Isoxazole and Isothiazole in the Isolated State and in Solution.

This follow-up paper completes the author's investigations to explore the in-solution structural preferences and relative free energies of all OH-substituted oxazole, thiazole, isoxazole, and isothiazole systems. The polarizable continuum dielectric solvent method calculations in the integral-equation formalism (IEF-PCM) were performed at the DFT/B97D/aug-cc-pv(q+(d))z level for the stable neutral tautomers with geometries optimized in dichloromethane and aqueous solution. With the exception of the predictions for the predominant tautomers of the 3OH isoxazole and isothiazole, the results of the IEF-PCM calculations for identifying the most stable tautomer of the given species in the two selected solvents agreed with those from experimental investigations. The calculations predict that the hydroxy proton, with the exception for the 4OH isoxazole and 4OH isothiazole, moves preferentially to the ring nitrogen or to a ring carbon atom in parallel with the development of a C=O group. The remaining, low-fraction OH tautomers will not be observable in the equilibrium compositions. Relative solvation free energies obtained by the free energy perturbation method implemented in Monte Carlo simulations are in moderate accord with the IEF-PCM results, but consideration of the ΔGsolv/MC values in calculating ΔG(s)tot maintains the tautomeric preferences. It was revealed from the Monte Carlo solution structure analyses that the S atom is not a hydrogen-bond acceptor in any OH-substituted thiazole or isothiazole, and the OH-substituted isoxazole and oxazole ring oxygens may act as a weak hydrogen-bond acceptor at most. The molecules form 1.0-3.4 solute-water hydrogen bonds in generally unexplored numbers at some specific solute sites. Nonetheless, hydrogen-bond formation is favorable with the NH, C=O and OH groups.

Theoretical in-Solution Conformational/Tautomeric Analyses for Chain Systems with Conjugated Double Bonds Involving Nitrogen(s).

Conformational/tautomeric transformations for X=CH-CH=Y structures (X = CH2, O, NH and Y = NH) have been studied in the gas phase, in dichloromethane and in aqueous solutions. The paper is a continuation of a former study where s-cis/s-trans conformational equilibria were predicted for analogues. The s-trans conformation is preferred for the present molecules in the gas phase on the basis of its lowest internal free energy as calculated at the B97D/aug-cc-pvqz and CCSD(T)CBS (coupled-cluster singles and doubles with non-iterative triples extrapolated to the complete basis set) levels. Transition state barriers are of 29-36 kJ/mol for rotations about the central C-C bonds. In solution, an s-trans form is still favored on the basis of its considerably lower internal free energy compared with the s-cis forms as calculated by IEF-PCM (integral-equation formalism of the polarizable continuum dielectric solvent model) at the theoretical levels indicated. A tetrahydrate model in the supermolecule/continuum approach helped explore the 2solute-solvent hydrogen bond pattern. The calculated transition state barrier for rotation about the C-C bond decreased to 27 kJ/mol for the tetrahydrate. Considering explicit solvent models, relative solvation free energies were calculated by means of the free energy perturbation method through Monte Carlo simulations. These calculated values differ remarkably from those by the PCM approach in aqueous solution, nonetheless the same prevalent conformation was predicted by the two methods. Aqueous solution structure-characteristics were determined by Monte Carlo. Equilibration of conformers/tautomers through water-assisted double proton-relay is discussed. This mechanism is not viable, however, in non-protic solvents where the calculated potential of mean force curve does not predict remarkable solute dimerization and subsequent favorable orientation.

Competing intramolecular vs. intermolecular hydrogen bonds in solution.

A hydrogen bond for a local-minimum-energy structure can be identified according to the definition of the International Union of Pure and Applied Chemistry (IUPAC recommendation 2011) or by finding a special bond critical point on the density map of the structure in the framework of the atoms-in-molecules theory. Nonetheless, a given structural conformation may be simply favored by electrostatic interactions. The present review surveys the in-solution competition of the conformations with intramolecular vs. intermolecular hydrogen bonds for different types of small organic molecules. In their most stable gas-phase structure, an intramolecular hydrogen bond is possible. In a protic solution, the intramolecular hydrogen bond may disrupt in favor of two solute-solvent intermolecular hydrogen bonds. The balance of the increased internal energy and the stabilizing effect of the solute-solvent interactions regulates the new conformer composition in the liquid phase. The review additionally considers the solvent effects on the stability of simple dimeric systems as revealed from molecular dynamics simulations or on the basis of the calculated potential of mean force curves. Finally, studies of the solvent effects on the type of the intermolecular hydrogen bond (neutral or ionic) in acid-base complexes have been surveyed.

Are the intramolecular O-H···F and O-H···Cl hydrogen bonds maintained in solution? A theoretical study.

The present case study aims at calculating the equilibrium conformer compositions for 2X-ethanol and 2X-phenol (X = F, Cl) in solution, and exploring the effect of the applied theoretical method and basis set on the obtained results, as well as considering the usefulness of the continuum solvent approach in comparison with the explicit solvent Monte Carlo model utilizing the free energy perturbation method. Gas-phase optimizations at the DFT/B97D/aug-cc-pvtz and ab initio MP2/aug-cc-pvtz levels predicted structures in good agreement with the available experimental data for three test molecules. Because in-solution geometries change only slightly according to the IEF-PCM continuum solvent calculations in carbon tetrachloride and water, the two theoretical levels were applied further on, and complete basis set (CBS) relative internal free energies were estimated for the conformers under study. The predicted OCCF gauche/trans ratio for 2F-ethanol was well reproduced in comparison with available experimental compositions. The predominant gauche structure maintains an intramolecular hydrogen bond in carbon tetrachloride (HB structure), whereas HB and NoHB gauche conformers appear in nearly the same fraction in aqueous solution. The internally hydrogen-bonded conformer is predominant also for 2X-phenol species, as calculated on the basis of relative CBS internal free energies and IEF-PCM and FEP/MC solvation free energies. Use of a trihydrate supermolecule model for 2F-ethanol conformers leads to the prediction of the aqueous-solution composition in contrast to the experiment. Solution structure modeling predicts weak hydrogen-bond formation capacity for both the covalently bound F and Cl atoms, even in conformations where they are fully exposed to hydration.

In-solution conformational analysis of the XCYCH3 moiety for small esters and ethers with all combinations of X, Y = O, S.

Favorable steric and electrostatic fit of a ligand to a receptor is of central interest in theoretical drug design. This paper considers the effects of non-protic solvents, in comparison with the gas phase, on the preferred conformation of the XCYCH3 moiety of simple aliphatic esters and heterocyclic methyl ethers with all combinations of the X and Y atoms as oxygen and sulfur. An IEF-PCM/B97D/aug-cc-pv(t+d)z continuum dielectric solvent study in chloroform and acetonitrile explores the through-space polarization effect of the environment on the conformational preference, not affected by possible solute-solvent hydrogen bond formation. The inherently favored structure for the present molecules is important, since the hypothetical oxygen and sulfur lone-pairs point approximately in opposite directions in the cis conformation of esters, whereas the trans and gauche conformations for the methyl group in ethers define nearly parallel or perpendicular directionality for the lone pairs of the ring heteroatoms and the O or S atoms connecting to the ring. These different preferences for the studied two families of compounds allow for designing formation of hydrogen bonds with a protein in fairly different regions of the latter still within the ligand-binding cavity. For a fine-tuning of these hydrogen bonds, a replacement of an oxygen atom of the ligand by a sulfur atom could be a straightforward possibility.

Theoretical study of the gauche-trans equilibrium with and without an intramolecular hydrogen bond for +H3NCH2CH2X systems (X = OH, NH2, COO-) in solution.

The gauche-trans conformational equilibrium has been studied for (+)H(3)NCH(2)CH(2)X systems (X = OH, NH(2), COO(-)) and for the neutral model of the simplest ß-amino acid in aqueous solution and chloroform. Each structure exhibits an intramolecular hydrogen bond in the NCCX gauche conformation, which is necessarily disrupted in the local energy minimum trans form. Geometries were optimized at the IEF-PCM/B97D/aug-cc-pvtz level of theory and indicated that the solvent effect of chloroform vs. water is moderate when the geometries for the corresponding gauche and trans conformers are compared. The only remarkable difference was found for ß-alanine, which can exist in gauche, zwitterionic form in aqueous solution but not in chloroform. Total relative free energies were obtained as the sum of the relative IEF-PCM/B97D/aug-cc-pvtz internal free energy and ΔG(solv). The latter was calculated both by means of IEF-PCM and the explicit solvent FEP/Monte Carlo methods. The resulting ΔG(tot) values could differ by 1-2 kcal mol(-1) depending on the accepted ΔG(solv) value, but any calculation indicated that the internally bound gauche conformers are far more populated than the corresponding trans species. A difference of 2.6 kcal mol(-1) for ΔG(solv) by the two methods resulted in the preference of the zwitterionic gauche ß-alanine structure by FEP/MC in contrast to IEF-PCM favoring the neutral form in aqueous solution. Monte Carlo characterization of the solution structure in the first solvation shell of polar groups indicates that solvation of a N-H(+) bond is sensitive to its involvement in an intramolecular hydrogen bond. For the zwitterionic, gauche ß-alanine, almost no water oxygen can reach that N-H(+) bond, which is in strong interaction with the -COO(-) group by forming a six-membered ring and favorable local geometry for the hydrogen bond.

Structure simulations for the 0.22 and 1 molar aqueous dimethylammonium chloride solutions.

Monte Carlo simulations have been performed for characterizing the 0.22 and 1 molar aqueous dimethylammonium chloride solutions at p = 1 atm and T = 310 K. On the basis of potential of mean force curves for the two systems with increasing concentrations, the N···N separations of about 8.7 and 7.5 Å correspond to a vague and a more pronounced minimum, respectively. Nitrogen separations at the minima are considerably smaller than those what the cations would take if the solutes comprised uniform local solute density on a microscale. The derived N by N coordination numbers predict non-negligible cation association and concomitant inhomogeneity for the studied systems. The calculated N···N distance distribution in the molar solution indicates that about 12% of the N···N separations are shorter than 8.5 Å compared with R(N···N) = 11.84 Å corresponding to the closest distance in a uniform cation local density. Despite a global minimum of -1.79 ± 0.63 kcal mol(-1) at N···Cl separation of 3.24 Å, obtained from the pmf for the 0.22 molar model, the N by Cl coordination number is only 0.14 in the first coordination shell, suggesting frequent disruption of an N-H(+)···Cl(+) hydrogen bond in a relatively dilute solution. The expression for the chemical potential of a solute includes a concentration-dependent activity coefficient, whose varying values are expected to reflect the different degrees of solute association in the 0.2-1 molar range. Theoretical follow-up of the changes in the activity coefficient values is difficult, thus calculation of the K(c) equilibrium constant has been proposed by considering 1 molar solutions as the standard state.

Theoretical studies of the solvent effect on the conformation of the HO-C-C-X (X = F, NH2, NO2) moiety with competing intra- and intermolecular hydrogen bonds.

Theoretical calculations up to the ab initio IEF-PCM/CCSD(T)/CBS//IEF-PCM/B3LYP/6-311++G** and IEF-PCM/B97D/aug-cc-pvtz levels have been performed for 2X-ethanol and 2X-phenol systems with X = F, NH(2), NO(2) in chloroform and aqueous solution. The calculated relative free energies by means of the IEF-PCM continuum dielectric method do not differ very much at the DFT and ab initio levels. Application of explicit solvent models and the FEP/MC method for determining relative solvation free energies causes, however, large deviations in the predicted equilibrium compositions, although the predominant conformation for the solute is generally in agreement with that from the corresponding IEF-PCM calculations. Existence of an intramolecular hydrogen bond (HB structure) for species with the HO-C-C-X moiety is preferred compared with a conformation when the hydrogen bond is disrupted (NoHB) for the considered F- and NO(2)-substituted molecules both in chloroform and aqueous solution. For 2NH(2)-ethanol, the HB structure is predominant in chloroform, whereas the 93:7 ratio for the OCCN trans/gauche species was obtained in aqueous solution. 2NH(2)-phenol exhibits a subtle equilibrium of the HB and NoHB conformations in both solvents. Potential of mean force calculations predict about a 10% solute association for the trans 2NH(2)-ethanol solute even in the fairly dilute 0.22 molar solution, whereas direct MC simulations do not support the maintenance of a doubly hydrogen-bonded dimer. Aqueous solution characteristics, as coordination numbers and numbers of strongly bound water molecules to the solute at T = 298 K and p = 1 atm, correspond reasonably to the derived molecular structures.

Molecular dynamics and Monte Carlo simulations for the structure of the aqueous trimethylammonium chloride solution in the 0.2-1 molar range.

Molecular dynamics and Monte Carlo simulations have been performed for characterizing the structure of the 0.2 and 1 molar aqueous trimethylammonium chloride solutions. Atomic charges were derived through the CHELPG and RESP fits to the molecular electrostatic potentials calculated for the cation in water at the IEF-PCM/B3LYP level using the 6-31G* and 6-311++G** basis sets. Maxima and minima of the calculated radial distribution functions were not sensitive to the four atomic charge sets. Simulated solution structures suggest non-negligible solute-solute interactions and remarkable inhomogeneity at both concentrations. This means that equilibrium ratios, derived for conformers/tautomers by means of ab initio calculations with the IEF-PCM continuum dielectric solvent model, should be corrected for free energy changes following solute association when compared to experimental data obtained for the 0.1-1 molar aqueous solutions.

The catalytic effect of water on the keto-enol tautomerism. Pyruvate and acetylacetone: a computational challenge.

The catalytic effect of explicit water molecules on the keto-enol tautomerism in a system of biological interest (enolpyruvate) has been investigated at the B3LYP/6-31++G** level by exploring the potential energy surface in the presence of 1 or 2 water molecules and comparing the energy profiles to the direct tautomerisation one. The water-assisted mechanisms turned out to be more favourable than the direct ones, in agreement with what occurred for the acetylacetone tautomerism in the presence of a single water molecule. To compare the behaviour of water in pyruvate and in acetylacetone, the two-water-assisted mechanism has been also examined for the latter at the B3LYP/6-31G* level. A number of stationary points for both compounds in vacuo and in the presence of explicit water molecules have been computed with DFT coupled to larger basis sets and at the MP2 level. Two diketo forms more stable than any of the keto-enol tautomers have been located, while for pyruvate the keto form is always more favourable than the enol one. The equilibrium constant for acetylacetone tautomerisation has been computed with high accuracy, performing a complete basis set extrapolation for the relative internal energy, to determine whether the quality of the method/basis set is responsible for the earlier modest agreement with the experimental value. Monte Carlo and molecular dynamics simulations helped characterize the solution structure and association features in the 0.03-0.22 molar concentration range.

Monte Carlo structure simulations for aqueous 1,4-dioxane solutions.

Monte Carlo simulations in the NpT ensembles have been performed for the structure exploration of aqueous 1,4-dioxane solutions. Three different systems with all-atom dioxane:TIP4P water molar compositions of 2:500 (code:D2), 8:465 (D8), and 17:425 (D17) modeled solutions of 0.22, 0.88, and 1.86 mol/dm3 concentrations, respectively, at T = 298 K and p = 1 atm. The calculated solution densities increase from 0.992 to 1.002 g/cm3 with increasing dioxane concentration and approach the experimentally determined densities within 1%. This close agreement was achieved by utilizing RESP charges fitted to the in-solution IEF-PCM/B3LYP/6-31G* electrostatic potential of dioxane taken in its chair conformation and recently developed C, H steric parameters for ethers for calculations with a 12-6-1 all-atom potential. Solution structure analyses pointed out that the dioxane molecules arrange in the solutions with favorable distances of 4-8 angstroms for the ring symmetry centers. Within this range not only pairs of rings but triangular triads and tetrads have also been observed with center-center distances <8 angstroms. For the D8 system, about 25% of the sampled configurations included such a triad. In the case of the D17 model, two simulations starting from different solution configuration predicted different degrees for the dioxane aggregation in aqueous solution. In the more aggregated structure 3-21 triads are consistently maintained and 1-2 tetrads are formed in 58% of the configurations. Each dioxane oxygen forms about one hydrogen bond, on average, to a water molecule in the 0.22-1.86 molar range. The most likely O(dioxane)...H(water) hydrogen bond distance is 1.75-1.80 angstroms compared to the optimal distance of 1.72 angstroms in the isolated dimer. The optimal dioxane-water interaction energy of -5.65 kcal/mol indicates a remarkable hydrogen-bond acceptor character for dioxane.

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains.

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.

Design and synthesis of novel derivatives of the muscarinic agonist tetra(ethylene glycol)(3-methoxy-1,2,5-thiadiazol-4-yl) [3-(1-methyl-1,2,5,6-tetrahydropyrid-3-yl)-1,2,5-thiadiazol-4-yl] ether (CDD-0304): effects of structural modifications on the binding and activity at muscarinic receptor subtypes and chimeras.

As part of a continuing effort to design and synthesize highly selective muscarinic agonists for different muscarinic receptor subtypes, several tetra(ethylene glycol)(3-methoxy-1,2,5-thiadiazol-4-yl) [3-(1-methyl-1,2,5,6-tetrahydropyrid-3-yl)-1,2,5-thiadiazol-4-yl] ether (1) analogues were prepared and characterized. Different analogues were synthesized having hydrophilic spacers of di-, tri-, tetra-, penta(ethylene glycol) and tri(propylene glycol) separating the 1,2,5,6-tetrahydropyridine ring from the terminal heterocycle, which was either a 1,2,5-thiadiazole or 1,2,4-thiadiazole ring. Chimeric receptor and molecular modeling studies also were conducted to determine how the ligands interact with muscarinic receptors. The studies revealed that varying the distance of the terminal thiadiazole and the positioning of the methoxy group can increase binding affinity for certain muscarinic receptor subtypes (at M(2) for 13d and M(4) for 1) and enhance functional efficacy at M(4) receptors for 13e and 18b. Moreover, compound 1 exhibited antipsychotic activity as assessed by reversal of apomorphine-induced sensory motor gating deficits, suggesting potential utility in the treatment of schizophrenia.

Theoretical study of the enol imine <--> enaminone tautomeric equilibrium in organic solvents.

The tautomeric enol imine <--> enaminone (phenol <--> quinone) equilibrium of the 1-hydroxy-2-naphthaldehyde Schiff base (2-phenyliminomethyl-naphthalen-1-ol) was investigated by density functional theory (B3LYP) and ab initio (MP2) methods in the IEF-PCM polarizable continuum dielectric solvent approximation and by a combined ab initio + FEP/MC study by considering an explicit solvent model. Special emphasis was put on the effect of solvation on this equilibrium by using an apolar (CCl4), polar aprotic (CH3CN), and polar protic (CH3OH) solvent. Compared with experimental tautomerization Gibbs free energies, the IEF-PCM/B3LYP calculations apparently overestimate the stability of the quinone form both when the 6-31G(d,p) and the 6-311++G(d,p) basis sets are applied. IEF-PCM/MP2 studies with the above basis sets predict the preference of the aromatic phenol tautomer, in contrast to the experiment in methanol and acetonitrile solvent. Calculation of the total relative free energy as DeltaG(tot) = DeltaE(int)(IEF-PCM/QCISD(T)/6-31G(d)) + DeltaG(solv, FEP/MC) + DeltaG(thermal) provided agreement with the experimental values up to 0.6 kcal/mol in the three solvents, and the predominant tautomer was always correctly predicted. In-solution relevant atomic charges, derived by a fit to the molecular electrostatic potential generated by the IEF-PCM/B3LYP/6-31G(d,p) wave function, show strong dependence on the fitting procedure (CHELPG or RESP) and are fairly insensitive to the chemical nature of the actual solvent. Use of the CHELPG charges in FEP/MC simulations revealed to be superior in comparison with the use of the RESP charge set.

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures.

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.

Theoretical studies of the tautomeric equilibria for five-member N-heterocycles in the gas phase and in solution.

Tautomeric equilibria have been studied for five-member N-heterocycles and their methyl derivatives in the gas phase and in different solvents with dielectric constants of epsilon = 4.7-78.4. The free energy changes differently for tautomers upon solvation as compared to the gas phase, resulting in a shift of the equilibrium constant in solution. Solvents with increasing dielectric constant produce more negative solute-solvent interaction energies and increasing internal energies. The methyl-substituted imidazole and pyrrazole form delicate equilibria between two tautomeric forms. Depending on the solvent, the methyl-substituted triazoles and tetrazole have one or two major tautomers in solution. When estimating the relative solvation free energies by means of an explicit solvent model and using the FEP/MC method, one observes that the preferred tautomers differ in several cases from those predicted by the continuum solvent model. The 1,2-prototropic shift, as an intramolecular tautomerization path, requires about 50 kcal/mol activation energy for imidazole in the gas phase, and this route is also disfavored in a solution. The calculated activation free energy along the intramolecular path is 48-50 kcal/mol in chloroform and water as compared to a literature value of 13.6 kcal/mol for pyrrazole in DMSO. A molecular dynamics computer experiment favors the formation of an imidazole chain in chloroform, making the 1,3-tautomerization feasible along an intermolecular path in nonprotic solvents. In aqueous solution, one strong N-H...Ow hydrogen bond is formed for each species, whereas all other nitrogens in the ring form weaker, N...HwOw type hydrogen bonds. The tetrahydrofuran solvent acts as a hydrogen bond acceptor and forms N-H... Oether bonds. Molecules of the dichloromethane solvent are in favorable dipole-dipole interactions with the solute. The results obtained are useful in the design of N-heterocyclic ligands forming specified hydrogen bonds with protein side chains.

Synthesis and biological characterization of 1-methyl-1,2,5,6-tetrahydropyridyl-1,2,5-thiadiazole derivatives as muscarinic agonists for the treatment of neurological disorders.

Muscarinic agonists might be useful in the treatment of neurological disorders, including Alzheimer's disease, schizophrenia, chronic pain, and drug abuse. Previous studies identified a series of bis-1,2,5-thiadiazole derivatives of 1,2,5,6-tetrahydropyridine with high activity and selectivity for muscarinic receptors. To develop compounds with improved central nervous system penetration, several new derivatives were synthesized and characterized for muscarinic receptor binding and activity. One ligand (11) exhibited agonist activity at M(1), M(2), and M(4) receptors, a selectivity profile suggesting potential utility in the treatment of schizophrenia.

[EFfect of quinazolone-alkyl-carboxylic acid derivatives on the transmembrane Ca2+ ion flux mediated by AMPA receptors]. / Kinazolon-alkil-karbonsav származékok hatása az alpha-amino-3-hidroxi-5-metil-4-izoxazol propionsav (AMPA) receptor által szabályozott transzmembrán Ca2+ ion fluxusra.

The excitatory neurotransmitter, Glu, plays a crucial role in many sensory and motor functions as well as in brain development, learning and memory and it is also involved in the pathogenesis of a number of neurological disorders, including epilepsy, Alzheimer's and Parkinson's diseases. Therefore, the study of Glu receptors (GluRs) is of therapeutical importance. We showed here by fluorescence monitoring of transmembrane Ca2+ ion fluxes in response to (S)-alpha-amino-3-hidroxi-5-metil-4-izoxazol propionic acid ((S)-AMPA) on the time scale of 0.00004-10 s that Ca2+ ion influx proceeds through faster and slower desensitizing receptors. Pharmacological isolation of the slower and faster desensitizing AMPA receptor was possible by fluorescence monitoring of Ca2+ ion translocation in response to (S)-AMPA in the presence and absence of various 2-methyl-4-oxo-3H-quinazoline-3-alkyl-carboxilic acid derivatives (Qxs): the acetic acid Q1 inhibits the slower desensitizing receptor response specifically, while the acetyl-piperidine Q5 is a more potent inhibitor of the faster desensitizing receptor response. In addition, spontaneous interictal activity, as induced by high [K+] conditions in hippocampal slices, was reduced significantly by Q5, suggesting a possible anticonvulsant property of Q5. Substitutions of Qxs into the GluR2 S1S2 binding core were consistent with their effect by causing variable degree of S1S2 bridging interaction as one of the main determinants of AMPA receptor agonist activity. The exploitation of differences between similar receptors will be important in the development and use of drugs with high pharmacological specificity.

Differential Induction of Cytoplasmic Vacuolization and Methuosis by Novel 2-Indolyl-Substituted Pyridinylpropenones.

Because many cancers harbor mutations that confer resistance to apoptosis, there is a need for therapeutic agents that can trigger alternative forms of cell death. Methuosis is a novel form of non-apoptotic cell death characterized by accumulation of vacuoles derived from macropinosomes and endosomes. Previous studies identified an indole-based chalcone, 3-(5-methoxy-2-methylindol-3-yl)-1-(4-pyridinyl)-2-propen-1-one (MOMIPP), that induces methuosis in human cancer cells. Herein, we describe the synthesis of related 2-indolyl substituted pyridinylpropenones and their effects on U251 glioblastoma cells. Increasing the size of the 2-indolyl substituent substantially reduces growth inhibitory activity and cytotoxicity, but does not prevent cell vacuolization. Computational models suggest that the results are not due to steric-driven conformational effects. The unexpected uncoupling of vacuolization and cell death implies that the relationship between endosomal perturbations and methuotic cell death is more complex than previously realized. The new series of compounds will be useful in further defining the molecular and cellular mechanisms underlying methuosis.