Fernando D. Suvire

Tetrahydroisoquinolines as dopaminergic ligands: 1-Butyl-7-chloro-6-hydroxy-tetrahydroisoquinoline, a new compound with antidepressant-like activity in mice.

Three series of 1-substituted-7-chloro-6-hydroxy-tetrahydroisoquinolines (1-butyl-, 1-phenyl- and 1-benzyl derivatives) were prepared to explore the influence of each of these groups at the 1-position on the affinity for dopamine receptors. All the compounds displayed affinity for D(1)-like and/or D(2)-like dopamine receptors in striatal membranes, and were unable to inhibit [(3)H]-dopamine uptake in striatal synaptosomes. Different structure requirements have been observed for adequate D(1) or D(2) affinities. This paper details the synthesis, structural elucidation, dopaminergic binding assays, structure-activity relationships (SAR) of these three series of isoquinolines. Moreover, 1-butyl-7-chloro-6-hydroxy-tetrahydroisoquinoline (1e) with the highest affinity towards D(2)-like receptors (K(i) value of 66nM) and the highest selectivity (49-fold D(2) vs D(1)) by in vitro binding experiments was then evaluated in behavioral assays (spontaneous activity and forced swimming test) in mice. Compound 1e increased locomotor activity in a large dose range (0.04-25mg/kg). Furthermore, this lead compound produced reduction in immobility time in the forced swimming test at a dose (0.01mg/kg) that did not modify locomotor activity. The haloperidol (0.03mg/kg), a D(2) receptor preferred antagonist, blocked the antidepressant-like effect of compound 1e.

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration.

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

Tetrahydroisoquinolines acting as dopaminergic ligands. A molecular modeling study using MD simulations and QM calculations

A molecular modeling study on 16 1-benzyl tetrahydroisoquinolines (BTHIQs) acting as dopaminergic ligands was carried out. By combining molecular dynamics simulations with ab initio and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of BTHIQs interacting with the human dopamine D2 receptor (D2 DR) is reported here, providing a clear picture of the binding interactions of BTHIQs from both structural and energetic viewpoints. Molecular aspects of the binding interactions between BTHIQs and the D2 DR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental pKi values was obtained, predicting the potential dopaminergic effect of non-synthesized BTHIQs.

Ab Initio and DFT Study of the Conformational Energy Hypersurface of Cyclic Gly-Gly-Gly

The multidimensional conformational potential energy hypersurface (PEHS) of cyclic Gly-Gly-Gly (1,4,7-triazonane-2,5,8-trione) was comprehensively investigated at the Hartree-Fock (RHF/6-31G(d)) level of theory. The equilibrium structures, their relative stability, and the transition state (TS) structures involved in the conformational interconversion pathways were analyzed. aug-cc-pVTZ//B3LYP/6-311++G** single point calculations predict a trans-cis-cis conformation as the energetically preferred form for this compound. However, all of the levels of theory employed here predicted that two forms, a trans-cis-cis and a cis-cis-cis (crown), of conformers contribute significantly to the equilibrium mixture at room temperature. The conformational interconversion between the global minimum and the symmetric cis-cis-cis crown form requires 12.49 kcal/mol at the RHF 6-31G(d) level of theory, whereas the conformational interconversion between the cis-cis-cis crown and cis-cis-cis boat form requires 18.70 kcal/mol. An exploratory topological analysis of the PEHS was also carried out. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of these cyclic tripeptides, describing the conformations as well as the conformational interconversion processes in these hypersurfaces.

A potentiometric and spectrophotometric study on acid-base equilibria in ethanol-aqueous solution of acetazolamide and related compounds

Acid-base equilibria in ethanol-aqueous solution of 5-acetamido-1,3,4-thiadiazole-2-sulfonamide (acetazolamide, H(2)acm), 5-tertbutyloxycarbonylamido-1,3,4-thiadiazole-2-sulfonamide (B-H(2)ats), 5-amino-1,3,4-thiadiazole-2-sulfonamide (Hats) and 5-amino-1,3,4-thiadiazole-2-thiol (Hatm) at 25 degrees C, 0.15 mol dm(-3) ionic strength (NaNO(3)), have been investigated by potentiometry and UV spectrophotometry. The ionization constants were calculated with SUPERQUAD program from potentiometric measurements and by a method according to Edsall et al. using the mole fractions determined by complementary tri-stimulus colorimetry (CTS). The constants obtained by potentiometry were: B-H(2)ats, pk(a(1))=7.33(3) and pk(a(2))=9.27(1); Hats, pk(a(1))=2.51(3) and pk(a(2))=8.49(1); Hatm, pk(a(1))=1.92(1) and pk(a(2))=6.81(1); whereas the constants determined by spectrophotometry were: H(2)acm, pk(a(1))=7.78(1) and pk(a(2))=9.57(2); B-H(2)ats, pk(a(1))=7.71(2) and pk(a(2))=9.61(2); Hats, pk(a(1))=2.19(3) and pk(a(2))=8.61(2); Hatm, pk(a(2))=6.90(2). Theoretical calculations using MO semiempirical and ab-initio RHF/6-31G* computations for the compounds were also performed. It was possible to clarify the preferred deprotonation mechanism of acetazolamide and B-H(2)ats in which the first deprotonation takes place at the carbonamido group.

Alkaloids from Hippeastrum argentinum and Their Cholinesterase-Inhibitory Activities: An in Vitro and in Silico Study

Two new alkaloids, 4-O-methylnangustine (1) and 7-hydroxyclivonine (2) (montanine and homolycorine types, respectively), and four known alkaloids were isolated from the bulbs of Hippeastrum argentinum, and their cholinesterase-inhibitory activities were evaluated. These compounds were identified using GC-MS, and their structures were defined by physical data analysis. Compound 2 showed weak butyrylcholinesterase (BuChE)-inhibitory activity, with a half-maximal inhibitory concentration (IC50) value of 67.3 ± 0.09 μM. To better understand the experimental results, a molecular modeling study was also performed. The combination of a docking study, molecular dynamics simulations, and quantum theory of atoms in molecules calculations provides new insight into the molecular interactions of compound 2 with BuChE, which were compared to those of galantamine.

Synthesis, dopaminergic profile, and molecular dynamics calculations of N-aralkyl substituted 2-aminoindans

Brain dopaminergic system has a crucial role in the etiology of several neuropsychiatric disorders, including Parkinsons disease, depression, and schizophrenia. Several dopaminergic drugs are used to treat these pathologies, but many problems are attributed to these therapies. Within this context, the search for new more efficient dopaminergic agents with less adverse effects represents a vast research field. The aim of the present study was to synthesize N-[2-(4,5-dihydroxyphenyl)-methyl-ethyl]-4,5-dihydroxy-2-aminoindan hydrobromide (3), planned to be a dopamine ligand, and to evaluate its dopaminergic action profile. This compound was assayed as a diastereoisomeric mixture in two experimental models: stereotyped behavior (gnaw) and renal urinary response, after central administration. The pharmacological results showed that compound 3 significantly blocked the apomorphine-induced stereotypy and dopamine-induced diuresis and natriuresis in rats. Thus, compound 3 demonstrated an inhibitory effect on dopaminergic-induced behavior and renal action. N-[2-(-Methyl-ethyl)]-4,5-dihydroxy-2-aminoindan hydrobromide (4) was previously reported as an inotropic agent, and in the present work it was also re-evaluated as a diastereoisomeric mixture for its possible central action on the behavior parameters such as stereotypy and dopamine-induced diuresis and natriuresis in rats. Our results indicate that compound 4 produces an agonistic response, possibly through dopaminergic mechanisms. To better understand the experimental results we performed molecular dynamics simulations of two complexes: compound 3/D(2)DAR (dopamine receptor) and compound 4/D(2)DAR. The differential binding mode obtained for these complexes could explain the antagonist and agonist activity obtained for compounds 3 and 4, respectively.

Dynamics of flexible cycloalkanes. Ab initio and DFT study of the conformational energy hypersurface of cyclononane.

The multidimensional Potential Energy Hypersurface (PEHS) for the cyclononane molecule was comprehensively investigated at the Hartree–Fock (HF), and Density Functional Theory (DFT) levels of theory. Second‐order Møller–Plesset perturbation theory (MP2) optimizations were also carried out to confirm the low‐energy conformations. The previously reported Geometrical Algorithm to Search Conformational Space (GASCOS) has been used to generate the starting geometries for the conformational analysis. The GASCOS algorithm combined with ab initio and DFT optimization permits searching of the potential energy hypersurface for all minimum‐energy conformations as well as transition structures connecting the low‐energy forms. The search located all previously reported structures together with 11 transition states, some of which were not found by earlier searching techniques. Altogether, 16 geometries (five low‐energy conformations and 11 transition states) were found to be important for a description of the conformational features of cyclononane. RB3LYP/aug‐cc‐pVTZ//RB3LYP/6‐31G(d) calculations suggest a conformational mixture between the twist boat–chair and twist chair–boat conformations as the preferred forms. In addition only the twist chair–chair conformation with 1.52 kcal/mol above the global minimum should contribute somewhat to the equilibrium mixture of conformations. Our results allow us to form a concise idea about the internal intricacies of the 9D vector space describing the conformation of cyclononane as well as the associated conformational potential energy hypersurface of nine independent variables.

Correlation Between Experimental and DFT/GIAO Computed 13C NMR Chemical Shifts of Organic Compounds. Effects of Asymmetry

The experimental 13 C chemical shifts of five different series of organic compounds are compared with pre- dicted 13 C NMR chemical shifts obtained via empirically scaled GIAO shieldings. Our results indicate that the inclusion of a scaling factor allow to obtain an excellent correlation between δcalc and δexp. Although the inclusion of asymmetry im- proves this correlation, such enhancement was not observed for all the 54 tested compounds. We found RMS, a parameter related with the structural feature of the whole molecule, which could indicate the benefits of including asymmetry in these calculations.

Ring inversion in 1,4,7 cyclononatriene and analogues: Ab initio and DFT calculations and topological analysis

The multidimensional conformational potential energy hypersurfaces (PEHSs) for cis‐cis‐cis 1,4,7 cyclononatriene (I), Tribenzocyclononatriene (TBCN) (II), and cis‐cis‐cis cyclic triglycine (III) were comprehensively investigated at the Hartree–Fock (HF/6‐31G(d)) and density functional theory (B3LYP/6‐31G(d,p)) levels of theory. The equilibrium structures, their relative stability, and the transition state (TS) structures involved in the conformational interconversion pathways were analyzed. Altogether, four geometries (two low‐energy conformations and two transition states) were found to be important for a description of the conformational features of compounds I–III. B3LYP/aug‐cc‐pvdz//B3lYP/6‐31G(d,p) and MP2/6‐31G(d,p)//B3LYP/6‐31G(d,p) single point calculations predict that the conformational interconversion between crown and twist forms requires 14.01, 26.71, and 17.79 kcal/mol for compounds I, II, and III, respectively, which is in agreement with the available experimental data. A topological study of the conformational PEHSs of compounds I–III was performed. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of compounds I–III, describing the conformations as well as the conformational interconversion process in these hypersurfaces.