Lupeol Acetate and α-Amyrin Terpenes Activity against Trypanosoma cruzi: Insights into Toxicity and Potential Mechanisms of Action.

BACKGROUND: Chagas disease is a potentially fatal disease caused by the parasite Trypanosoma cruzi. There is growing scientific interest in finding new and better therapeutic alternatives for this disease's treatment. METHODS: A total of 81 terpene compounds with potential trypanocidal activity were screened and found to have potential T. cruzi cysteine synthase (TcCS) inhibition using molecular docking, molecular dynamics, ADME and PAIN property analyses and in vitro susceptibility assays. RESULTS: Molecular docking analyses revealed energy ranges from -10.5 to -4.9 kcal/mol in the 81 tested compounds, where pentacyclic triterpenes were the best. Six compounds were selected to assess the stability of the TcCS-ligand complexes, of which lupeol acetate (ACLUPE) and α-amyrin (AMIR) exhibited the highest stability during 200 ns of molecular dynamics analysis. Such stability was primarily due to their hydrophobic interactions with the amino acids located in the enzyme's active site. In addition, ACLUPE and AMIR exhibited lipophilic characteristics, low intestinal absorption and no structural interferences or toxicity. Finally, selective index for ACLUPE was >5.94, with moderate potency in the trypomastigote stage (EC50 = 15.82 ± 3.7 µg/mL). AMIR's selective index was >9.36 and it was moderately potent in the amastigote stage (IC50 = 9.08 ± 23.85 µg/mL). CONCLUSIONS: The present study proposes a rational approach for exploring lupeol acetate and α-amyrin terpene compounds to design new drugs candidates for Chagas disease.

A Theoretical Analysis of Interaction Energies and Intermolecular Interactions between Amphotericin B and Potential Bioconjugates Used in the Modification of Nanocarriers for Drug Delivery.

Amphotericin B (AmB) is an antibiotic with a wide spectrum of action and low multidrug resistance, although it exhibits self-aggregation, low specificity, and solubility in aqueous media. An alternative for its oral administration is its encapsulation in polymers modified with bioconjugates. The aim of the present computational research is to determine the affinity between AmB and six bioconjugates to define which one could be more suitable. The CAM-B3LYP-D3/6-31+G(d,p) method was used for all computational calculations. The dimerization enthalpy of the most stable and abundant systems at pH = 7 allows obtaining this affinity order: AmB_1,2-distearoyl-sn-glycerol-3-phosphorylethanolamine (DSPE) > AmB_γ-cyclodextrin > AmB_DSPEc > AmB_retinol > AmB_cholesterol > AmB_dodecanol, where DSPEc is a DSPE analog. Quantum theory of atoms in molecules, the non-covalent interactions index, and natural bond orbital analysis revealed the highest abundance of noncovalent interactions for AmB-DSPE (51), about twice the number of interactions of the other dimers. Depending on the interactions' strength and abundance of the AmB-DSPE dimer, these are classified as strong: O-H---O (2), N-H---O (3) and weak: C-H---O (25), H---H (18), C-H---C (3). Although the C-H---O hydrogen bond is weak, the number of interactions involved in all dimers cannot be underestimated. Thus, non-covalent interactions drive the stabilization of copolymers, and from our analysis, the most promising candidates for encapsulating are DSPE and γ-cyclodextrin.

Experimental and theoretical study of the mechanism and rate constants of the sequential 5-exo-trig spirocyclization involving vinyl, aryl and N-alkoxyaminyl radicals.

In this research the sequential generation and cyclization of N-alkoxyaminyl radicals to produce 1-azaspiro[4.4]nonane, a prominent scaffold in organic and medicinal chemistry, was studied. Competition experiments in benzene at 80 °C with brominated oxime ethers using Bu3SnH as chain transfer and AIBN as the initiator generated vinyl or aryl radicals which were captured by oxime ethers, allowing approximate 5-exo-trig cyclization constants at 4.6 × 108 s-1 and 9.9 × 108 s-1 respectively to be established. Similar results were obtained by kinetic studies using the transition state theory (TST) from ab initio calculations with density functional theory (DFT) using the M06-2X, B3LYP, mPW1PW91 and TPSSh functionals in combination with the 6-311+G(d, p) basis set. Additionally, it was found that the 5-exo-trig cyclization of the N-alkoxyaminyl radical onto CC double bonds is a reversible process whose constants were determined to be in the range of 6.2 × 100 s-1 and 3.5 × 106 s-1 at 80 °C, depending on the nature of the substituents. The calculation of the radical stabilization energy (RSE) shows that the N-alkoxyaminyl radical is a very stable species and its reactivity in the addition on alkenes is governed by its nucleophilic character and the stability of the carbon-centered radical formed after cyclization. The reduction constant of the N-alkoxyaminyl radical with Bu3SnH in the gas phase at 80 °C was also estimated to be 3.4 × 100 M-1 s-1 through computational calculations. This information facilitates the rational planning of cascades and other methodologies applied to the construction of carbocyclic and aza-heterocyclic compounds.

Estudio teórico de la oxidación de CO con O2 usando catalizadores de Au-Pd y Au-Pt / Theoretical study of the CO Oxidation with O2 using Au-Pd and Au-Pt catalysts / Estudo teórico da oxidação de CO com O2 usando catalisadores de Au-Pd e Au-Pt

Resumen En el presente estudio se realizaron cálculos con base en la Teoría del Funcional de la Densidad Electrónica (DFT) con la aproximación B3PW91/LANL2DZ para optimizar los sistemas monometálicos y bimetálicos Au9, Au8Pd, Au8Pt, AuPd8, AuPt8, Pd9 y Pt9. Los materiales fueron teóricamente evaluados como catalizadores para la oxidación de monóxido de carbono (CO) y se determinó el sistema más favorable para la adsorción de esta molécula. La sustitución de átomos de Pt y Pd por átomos de Au en los nonámeros generó un cambio en la estructura tridimensional del sistema. El análisis de reactividad global mostró que el clúster más reactivo es PÍ9, seguido por AuPt s . Los índices de Fukui identificaron los sitios más susceptibles para un ataque nucleofílico de ambos clústeres. La adsorción de CO generó una cascada de oxidación que liberó ~4,5 eV, indicando que la reacción es altamente exotérmica y exergónica. Los clústeres AuPt s y Pt 9 mostraron los valores más bajos de energía de activación de la etapa determinante del mecanismo. En general, la sustitución de un átomo de platino (o paladio) por un átomo de oro no afecta la reactividad de los nonámeros y, por tanto, se infiere que el clúster AuPt s podría ser un catalizador promisorio en la oxidación de CO.

Abstract In the current study were development calculations based on Density Functional Theory (DFT) with the B3PW91/LANL2DZ approach for optimizing both monometallic and bimetallic systems: Au9, AusPd, Au8Pt, AuPds, AuPts, Pd9 y Pt9. Such materials were theoretically tested as catalyst for the oxidation of carbon monoxide (CO) and the most favorable system for its further adsorption was determined. The substitution of Pt and Pd by Au atoms in the nonamers generated a change in the tridimensional structure of the system. The global reactivity analysis showed that the most reactive cluster is Pt9 followed by AuPts. On the other hand, the Fukui indexes identified the most susceptible sites for a nucleophilic attack of both clusters. The CO adsorption generated an oxidation cascade which liberated ∼ 4.5 eV, indicating that the reaction is highly exothermic and exergonic. Both AuPt8 and Pt9 showed the lowest values of activation energy in the determining step of the mechanism. In general, the substitution of a Pt (Pd) atom by an Au atom does not affect the reactivity of the nonamers and then it is inferred that the AuPds cluster could be a promissory catalyst in the CO oxidation.

Resumo No presente estudo, cálculos baseados na Teoria do Funcional da Densidade Eletrônica (DFT) com a abordagem B3PW91/LANL2DZ foram realizados para otimizar sistemas monometálicos e bimetálicos Au9, Au8Pd, Au8Pt, AuPd8, AuPt8, Pd9 y Pt9. Tais materiais foram teoricamente avaliados como catalisadores para a oxidação do monóxido de carbono (CO) e foi determinado o sistema mais favorável para a adsorção desta molécula. A substituição dos átomos de Pt e Pd por átomos de Au nós não-nomes gerou uma mudança na estrutura tridimensional do sistema. A análise de reatividade global mostrou que o cluster mais reativo é Pt9, seguido por AuPt8. Os índices de Fukui identificaram os sítios mais suscetíveis ao ataque nucleofílico de ambos os clusters. A adsorção de CO gerou uma cascata de oxidação que liberou ~4,5 eV, indicando que a reação é altamente exotérmica e exergônica. Os aglomerados AuPt 8 y Pt 9 apresentaram os menores valores de energia de ativação do estágio determinante do mecanismo. Em geral, a substituição de um átomo de platina (ou paládio) por um átomo de ouro não afeta a reatividade dos não-nomes e, portanto, infere-se que o aglomerado AuPt 8 pode ser um catalisador promissor na oxidação do CO.

Thermodynamics of the Isomerization of Monoterpene Epoxides.

In this contribution, the thermodynamic analysis of α- and ß-pinene epoxide isomerization over Fe and Cu supported on MCM-41 is presented using computational chemistry and group contribution methods (GCMs). Some physical-chemical data (T c, P c, v c, Z c, ω, T b, T fus) and thermodynamic (S°298.15, C p,298.15 °, C v,298.15 °, ΔH f,298.15 °, ΔG f,298.15 °, ΔH vb °, ΔH fus, C pL) properties obtained by different GCMs are reported for several monoterpenes and monoterpenoids, which significantly contribute to the knowledge of the properties of these compounds. Density functional theory (DFT), PBE-D3/6-311G(d,p), was employed for determining the Gibbs free energy and the heat of reaction associated with the transformation of monoterpene epoxides into aldehydes, ketones, and related oxygenated compounds in the presence of different solvents and at several temperatures. The calculations were compared with available data reported and the experimental results of the catalytic reactions. The transformation of α- and ß-pinene epoxides into aldehydes appears to be more spontaneous and favorable than their transformations into alcohols in a wide range of temperatures. These results are in agreement with the experiments over Fe/MCM-41 and Cu/MCM-41, where α-pinene epoxide isomerization yields campholenic aldehyde (50-80% selectivity) as the main product. The 1.7Fe/MCM-41 material was more active in all solvents than 1.3Cu/MCM-41 for both α- and ß-pinene epoxide isomerization. However, perillyl alcohol (20-70% selectivity) was the most favored for the isomerization reaction, except when ethyl acetate was the solvent. Enthalpy and Gibbs free energy of the studied reactions estimated by both GCMs and DFT calculations did not show large differences for most of the reactions at evaluated temperatures.

Looking for the Azeotrope: A Computational Study of (Ethanol)6-Water, (Methanol)6-Water, (Ethanol)7, and (Methanol)7 Heptamers.

Considering that a molecular-level understanding of the azeotropic ethanol-water system can contribute to the search of new methodologies and/or modifications of industrial separation methods, this study tries to provide some clues to understand why azeotropes should be expected for ethanol, but not for methanol. Our exploration of the potential energy surface of (ethanol)6-water heteroheptamers, carried out at the B3LYP-D3/6-311++G(d,p) level, shows these heteroclusters to exhibit a cyclic structure where the cooperativity effects between the OH···O HBs is a fundamental ingredient. An analysis of this cooperativity clearly indicates that ethanol-water systems will exhibit a similarly high stability as the heterocluster size approaches the azeotrope. However, a similar behavior should not be expected for the methanol-containing analogues. A comparison between (ethanol)7, (ethanol)6-water, (methanol)7, and (methanol)6-water shows the ethanol-containing systems to be significantly more stable than the methanol-containing analogues. This result is probably due to the fact that the OH···O HBs are weaker than those found between ethanol molecules. However, our atoms in molecule (AIM) and noncovalent interaction (NCI) analyses unambiguously show that important contributors to the enhanced stability of the ethanol-containing clusters are the secondary van der Waals interactions between ethyl groups, which are not observed between methyl groups. Hence, while the formation of stable azeotropes is expected for the case of ethanol, for the methanol-containing analogues, the relative stability of the clusters is significantly smaller, and its formation is accompanied by an increase of the free energy.

Structural and functional computational analysis of nicotine analogs as potential neuroprotective compounds in Parkinson disease.

As the mechanism of interaction between nicotinic receptors with nicotine analogs is not yet fully understood, information at molecular level obtained from computational calculations is needed. In this sense, this work is a computational study of eight nicotine analogs, all with pyrrolidine ring modifications over a nicotine-based backbone optimized with B3LYP-D3/aug-cc-pVDZ. A molecular characterization was performed focusing on geometrical parameters such as pseudo-rotation angles, atomic charges, HOMO and LUMO orbitals, reactivity indexes and intermolecular interactions. Three analogs, A2 (3-(1,3-dimethyl-4,5-dihydro-1h-pirazole-5-yl) pyridine), A3 (3-(3-methyl-4,5-dihydro-1H-pyrazol-5-yl)-pyridine) and A8 (5-methyl-3-(pyridine-3-yl)-4,5-dihydroisoxazole), were filtered suggesting putative neuroprotective activity taking into account different reactivity values, such as their lowest hardness: 2.37 eV (A8), 2.43 eV (A2) and 2.56 eV (A3), compared to the highest hardness value found: 2.71 eV for A5 (3-((2S,4R)-4-(fluoromethyl)-1-methylpyrrolidine-2-il) pyridine), similar to the value of nicotine (2.70 eV). Additionally, molecular docking of all 8 nicotine analogs with the α 7 nicotinic acetylcholine receptor (α 7 nAChR) was performed. High values of interaction between the receptor and the three nicotine analogs were obtained: A3 (-7.1 kcal/mol), A2 (-6.9 kcal/mol) and A8 (-6.8 kcal/mol); whereas the affinity energy of nicotine was -6.4 kcal/mol. Leu116 and Trp145 are key residues in the binding site of α 7 nAChR interacting with nicotine analogs. Therefore, based upon these results, possible application of these nicotine analogs as neuroprotective compounds and potential implication at the design of novel Parkinson's treatments is evidenced.

Reaction Mechanism of the Isomerization of Monoterpene Epoxides with Fe3+ as Active Catalytic Specie: A Computational Approach.

The reaction mechanism of the isomerization of α and ß-pinene epoxides with Fe species as catalysts was studied with density functional theory (DFT) calculations and an experimental methodology. ß-pinene epoxide can be isomerized into myrtanal and myrtenol in four steps, while in the case of perillyl alcohol, two additional steps are necessary. On the other hand, high selectivity to myrtanal obtained experimentally can be explained by the number of steps and the kind of the hydrogen transference that is easier for this compound in comparison with myrtenol and perillyl alcohol. A thermodynamic analysis showed that transformation into myrtenol, myrtanal, and perillyl alcohol is spontaneous but transformation into myrtanal is the most favorable. In the case of α-pinene epoxide rearrangement, a mechanistic study was carried out toward the optimization of the possible intermediates. Synthesis of campholenic aldehyde and carveol from α-pinene epoxide was carried out through three steps after the coordination of oxygen to iron, showing that in contrast to carveol formation, campholenic aldehyde synthesis is spontaneous. Analysis of ∇2ρ, the total energy density (H = V + G), and the |V|/G ratio evaluated at the bond critical point of the Fe-O bond showed for both epoxides that such interaction is closed shell instead of covalent. Apparently, α-pinene epoxide can be isomerized faster that ß-pinene epoxide with Fe3+, which is in agreement with previous experimental results. This is the first report where a reaction mechanism of isomerization of monoterpenes epoxides is studied based on very detailed experimental and computational methodologies.

Effect of Charge and Phosphine Ligands on the Electronic Structure of the Au8 Cluster.

In this work, we use density functional theory calculations with a hybrid exchange-correlation functional and effective core pseudopotentials to determine the geometry of bare and phosphine-protected Au8 nanoclusters and characterize their electronic structure. Au8 clusters were bonded to four and eight PH3 ligands in order to evaluate the effect of ligand concentration on the electronic structure, while different positional configurations were also tried for four ligands attached to the cluster. We show that the neutral clusters become more nucleophilic as the ligands bind to the clusters at stable sites. The ground-state planar configuration of Au8 is maintained depending on the concentration and position of ligands. The effect of ionizing to the +2 charge state results in disruption of planar geometry in some cases because of inoccupation of a molecular orbital with the Au-Au bonding character. Natural bond order charge analyses showed that Au atoms oxidize upon ionization, instead of phosphine. The net positive charge makes the clusters more electrophilic with a capacity to absorb electrons from nucleophiles depending on the concentration and position of ligands and on the concentration of low-coordinated gold atoms. Besides, ionization energies and electron affinities were calculated through different mechanisms, finding that both variables are much higher for charged systems and change inversely with the concentration of ligands.

The role of OH O and CH O hydrogen bonds and H H interactions in ethanol/methanol-water heterohexamers.

Bioethanol is one of the world's most extensively produced biofuels. However, it is difficult to purify due to the formation of the ethanol-water azeotrope. Knowledge of the azeotrope structure at the molecular level can help to improve existing purification methods. In order to achieve a better understanding of this azeotrope structure, the characterization of (ethanol)5-water heterohexamers was carried out by analyzing the results of electronic structure calculations performed at the B3LYP/6-31+G(d) level. Hexamerization energies were found to range between -36.8 and -25.8 kcal/mol. Topological analysis of the electron density confirmed the existence of primary (OH…O) hydrogen bonds (HBs), secondary (CH…O) HBs, and H…H interactions in these clusters. Comparison with three different solvated alcohol systems featuring the same types of atom-atom interactions permitted the following order of stability to be determined: (methanol)5-water > (methanol)6 > (ethanol)5-water > (ethanol)6. These findings, together with accompanying geometric and spectroscopic analyses, show that similar cooperative effects exist among the primary HBs for structures with the same arrangement of primary HBs, regardless of the nature of the molecules involved. This result provides an indication that the molecular ratio can be considered to determine the unusual behavior of the ethanol-water system. The investigation also highlights the presence of several types of weak interaction in addition to primary HBs. Graphical Abstract Water-ethanol clusters exhibit a variety of interaction types between their atoms, such as primary OH...O (blue), secondary CH...O (green) and H...H (yellow) interactions as revealed by Quantum Chemical Topology.

Understanding the azeotropic diethyl carbonate-water mixture by structural and energetic characterization of DEC(H2O)(n) heteroclusters.

Diethyl carbonate (DEC) is an oxygenated fuel additive. During its synthesis through a promising green process, a DEC-water azeotrope is formed, which decreases DEC production efficiency in the gas phase. Molecular information about this system is scarce but could be of benefit in understanding (and potentially improving) the synthetic process. Therefore, we report a detailed computational study of the conformers of DEC, and their microsolvation with up to four water molecules, with the goal of understanding the observed 1:3 DEC:H2O molar ratio. The most stable DEC conformers (with mutual energy differences < 1.5 kcal mol(-1)) contribute to the energetic and structural properties of the complexes. An exhaustive stochastic exploration of each potential energy surface of DEC-(H2O)n, (where n = 1, 2, 3, 4) heteroclusters discovered 3, 8, 7, and 4 heterodimers, heterotrimers, heterotetramers, and heteropentamers, respectively, at the MP2/6-311++G(d,p) level of theory. DEC conformers and energies of the most stable structures at each heterocluster size were refined using CCSD(T)/6-311++G(d,p). Energy decomposition, electron density topology, and cooperative effects analyses were carried out to determine the relationship between the geometrical features of the heteroclusters and the non-covalent interaction types responsible for their stabilization. Our findings show that electrostatic and exchange energies are responsible for heterocluster stabilization, and also suggest a mutual weakening among hydrogen bonds when more than three water molecules are present. All described results are complementary and suggest a structural and energetic explanation at the molecular level for the experimental molar ratio of 1:3 (DEC:H2O) for the DEC-water azeotrope.

Phosphine-stabilised Au9 clusters interacting with titania and silica surfaces: the first evidence for the density of states signature of the support-immobilised cluster.

Chemically made, atomically precise phosphine-stabilized clusters Au9(PPh3)8(NO3)3 were deposited on titania and silica from solutions at various concentrations and the samples heated under vacuum to remove the ligands. Metastable induced electron spectroscopy was used to determine the density of states at the surface, and X-ray photoelectron spectroscopy for analysing the composition of the surface. It was found for the Au9 cluster deposited on titania that the ligands react with the titania substrate. Based on analysis using the singular value decomposition algorithm, the series of MIE spectra can be described as a linear combination of 3 base spectra that are assigned to the spectra of the substrate, the phosphine ligands on the substrate, and the Au clusters anchored to titania after removal of the ligands. On silica, the Au clusters show significant agglomeration after heat treatment and no interaction of the ligands with the substrate can be identified.

An orbital and electron density analysis of weak interactions in ethanol-water, methanol-water, ethanol and methanol small clusters.

A computational study of (ethanol)(n)-water, n = 1 to 5 heteroclusters was carried out employing the B3LYP∕6-31+G(d) approach. The molecular (MO) and atomic (AO) orbital analysis and the topological study of the electron density provided results that were successfully correlated. Results were compared with those obtained for (ethanol)(n), (methanol)(n), n = 1 to 6 clusters and (methanol)(n)-water, n = 1 to 5 heteroclusters. These systems showed the same trends observed in the (ethanol)(n)-water, n = 1 to 5 heteroclusters such as an O---O distance of 5 Å to which the O-H---O hydrogen bonds (HBs) can have significant influence on the constituent monomers. The HOMO of the hetero(clusters) is less stable than the HOMO of the isolated alcohol monomer as the hetero(cluster) size increases, that destabilization is higher for linear geometries than for cyclic geometries. Changes of the occupancy and energy of the AO are correlated with the strength of O-H---O and C-H---O HBs as well as with the proton donor and/or acceptor character of the involved molecules. In summary, the current MO and AO analysis provides alternative ways to characterize HBs. However, this analysis cannot be applied to the study of H---H interactions observed in the molecular graphs.

The dynamic behavior of a liquid ethanol-water mixture: a perspective from quantum chemical topology.

Quantum Chemical Topology (QCT) is used to reveal the dynamics of atom-atom interactions in a liquid. A molecular dynamics simulation was carried out on an ethanol-water liquid mixture at its azeotropic concentration (X(ethanol)=0.899), using high-rank multipolar electrostatics. A thousand (ethanol)(9)-water heterodecamers, respecting the water-ethanol ratio of the azeotropic mixture, were extracted from the simulation. Ab initio electron densities were computed at the B3LYP/6-31+G(d) level for these molecular clusters. A video shows the dynamical behavior of a pattern of bond critical points and atomic interaction lines, fluctuating over 1 ns. A bond critical point distribution revealed the fluctuating behavior of water and ethanol molecules in terms of O-H···O, C-H···O and H···H interactions. Interestingly, the water molecule formed one to six C-H···O and one to four O-H···O interactions as a proton acceptor. We found that the more localized a dynamical bond critical point distribution, the higher the average electron density at its bond critical points. The formation of multiple C-H···O interactions affected the shape of the oxygen basin of the water molecule, which is shown in three dimensions. The hydrogen atoms of water strongly preferred to form H···H interactions with ethanol's alkyl hydrogen atoms over its hydroxyl hydrogen.

Efectos cooperaiivos en heterotetrámeros (etanol)3-auua / Cooperatve effeccts in (ethanol)3-water heterotetramers / Efeitos cooperativos em heterotetrámeros (etanol)3-água

La teoría de funcionales de la densidad (DFT: B3LYP/6-31+G(d)) fue empleada para la optimización de agregados estables sobre la superficie de energía potencial de los heterotetrámeros (etanol)3-agua. Las energías de tetramerización pueden llegar a valores hasta de -21,00 kcal/mol. Esta energía no se puede obtener considerando solo contribuciones de interacciones entre dos moléculas del agregado, lo cual sugiere la presencia de efectos cooperativos globales (positivos). Tales efectos son reflejados en distancias menores de los puentes de hidrógeno y distancias menores oxígeno-oxígeno, lo mismo que en elongaciones mayores del enlace O-Hen la molécula dadora de protón con un corrimiento hacia el rojo mayor en los heterotetrámeros, comparado con los heterodímeros de etanol-agua y el dímero de etanol. La mayor cooperatividad fue observada en los cuatro puentes de hidrógeno dispuestos en el patrón geométrico cíclico más grande posible, actuando todas las moléculas como aceptoras y dadoras de protón simultáneamente. Un análisis similar al de la caracterización de heterotetrámeros de (etanol)3-agua se llevó a cabo para los heterotetrámeros (metanol)3-agua y tetrámeros de etanol y metanol. La comparación de estos valores mostró que existe una gran similitud entre todos los parámetros analizados para agregados con el mismo patrón geométrico.

Density Functional Theory (DFT: B3LYP/6-31 + G(d)) was used for the optimization of clusters on the potential energy surface of (ethanol)3-water heterotetramers. The tetramerization energies can reach values up to -21.00 kcal/ mol. This energy can not be obtained by just considering the contributions from interactions between two cluster molecules, which suggests of the presence of global cooperative effects (positive). These effects are reflected in smaller hydrogen bond distances and smaller oxygen-oxygen distances, as well as in greater elongations of the O-H proton donor bond with a stronger "red-shift" in the heterotetramers compared to the ethanol-water heterodimers and the ethanol dimer. The largest cooperativity effect was observed in the four hydrogen bonds arranged in the largest possible cyclic geometric pattern, where all the molecules act as proton acceptor and donor simultaneously. A similar analysis to the characterization of (ethanol)3-water heterotetramers was carried out on (methanol)3-water heterotetramers, and ethanol and methanol tetramers, whose comparison showed a great similarity between all evaluated parameters for the clusters with equal geometric pattern.

A teoria de funcionais de densidade (DFT: B3LYP/6-31+G(d)) foi empregada para a otimização de agregados sobre a superfície de energia potencial dos heterotetrâmeros (etanol)3-água. As energias de tetramerização podem alcançar valores de até -21.00 kcal/mol. Esta energia não pode ser obtida por apenas considerando as contribuições das interações entre agregados de duas moléculas, o que sugere a presença global dos cooperativos efeitos (positivos). Tais efeitos são refletidos em menores comprimentos das pontes de hidrogênio e distâncias oxigênio-oxigênio, e também em maiores alongamentos da ligação O-H na molécula doadora de prótons, com um maior "red-shift" associado nos heterotetrâmeros do que nos heterodímeros de etanol-água e no dímero de etanol. A mais alta cooperatividade foi observada com as quatro pontes de hidrogênio dispostas no maior padrão geométrico cíclico possível, atuando simultaneamente todas as moléculas como aceptoras e doadoras de prótons. Uma análise similar ao da caracterização de heterotetrâmeros de (eta-nol)3-água se levou a cabo sobre heterotetrâmeros de (metanol)3-água e metanol, cuja comparação mostrou uma grande similaridade entre todos os parâmetros analisados para agregados com igual padrão geométrico.

Molecular interaction of (ethanol)2-water heterotrimers.

The potential energy surface of the (ethanol)2-water heterotrimers for the trans and gauche conformers of ethanol was studied using density functional theory. The same approximation was used for characterizing representative clusters of (ethanol)3, (methanol)3, and (methanol)2-water. Trimerization energies and enthalpies as well as the analysis of geometric parameters suggest that the structures with a cyclic pattern in the three hydrogen bonds of the type O-H---O (primary hydrogen bonds), where all molecules are proton donor-acceptor at the same time, are more stable than those with just two primary hydrogen bonds. Additionally, we propose the formation of "secondary hydrogen bonds" between hydrogen atoms of the methyl group of ethanol and the oxygen atom of water or other ethanol molecule (C-H---O), which were found to be weaker than the primary hydrogen bonds.