Exploring the Non-Covalent Bonding in Water Clusters.

QTAIM and source function analysis were used to explore the non-covalent bonding in twelve different water clusters (H2O)n obtained by considering n = 2-7 and various geometrical arrangements. A total of seventy-seven O-H⋯O hydrogen bonds (HBs) were identified in the systems under consideration, and the examination of the electron density at the bond critical point (BCP) of these HBs revealed the existence of a great diversity of O-H⋯O interactions. Furthermore, the analysis of quantities, such as |V(r)|/G(r) and H(r), allowed a further description of the nature of analogous O-H⋯O interactions within each cluster. In the case of 2-D cyclic clusters, the HBs are nearly equivalent between them. However, significant differences among the O-H⋯O interactions were observed in 3-D clusters. The assessment of the source function (SF) confirmed these findings. Finally, the ability of SF to decompose the electron density (ρ) into atomic contributions allowed the evaluation of the localized or delocalized character of these contributions to ρ at the BCP associated to the different HBs, revealing that weak O-H⋯O interactions have a significant spread of the atomic contributions, whereas strong interactions have more localized atomic contributions. These observations suggest that the nature of the O-H⋯O hydrogen bond in water clusters is determined by the inductive effects originated by the different spatial arrangements of the water molecules in the studied clusters.

Dehydrochlorination of PCDDs on SWCN-Supported Ni10 and Ni13 Clusters, a DFT Study.

Polychlorinated dibenzo-p-dioxins (PCDDs) are known to be a group of compounds of high toxicity for animals and, particularly, for humans. Given that the most common method to destroy these compounds is by high-temperature combustion, finding other routes to render them less toxic is of paramount importance. Taking advantage of the physisorption properties of nanotubes, we studied the reactions of atomic hydrogen on physisorbed PCDDs using DFT; likewise, we investigated the reaction of molecular hydrogen on PCDDs aided by Ni10 and Ni13 clusters adsorbed on single-wall carbon nanotubes. Because dihydrogen is an easily accessible reactant, we found these reactions to be quite relevant as dehydrohalogenation methods to address PCDD toxicity.

A Theoretical Study of the C-X Bond Cleavage Mediated by Cob(II)Aloxime.

The C-X bond cleavage in different methyl halides (CH3X; X = Cl, Br, I) mediated by 5,6-dimethylbenzimidazole-bis(dimethylglyoximate)cobalt(II) (CoIICbx) was theoretically investigated in the present work. An SN2-like mechanism was considered to simulate the chemical process where the cobalt atom acts as the nucleophile and the halogen as the leaving group. The reaction path was computed by means of the intrinsic reaction coordinate method and analyzed in detail through the reaction force formalism, the quantum theory of atoms in molecules (QTAIM), and the calculation of one-electron density derived quantities, such as the source function (SF) and the spin density. A thorough comparison of the results with those obtained in the same reaction occurring in presence of 5,6-dimethylbenzimidazole-bis(dimethylglyoximate)cobalt(I) (CoICbx) was conducted to reveal the main differences between the two cases. The reactions mediated by CoIICbx were observed to be endothermic and possess higher activation energies in contrast to the reactions where the CoICbx complex is present. The latter was supported by the reaction force results, which suggest a relationship between the activation energy and the ionization potentials of the different nucleophiles present in the cleavage reaction. Moreover, the SF results indicates that the lower axial ligand (i.e., 5,6-dimethylbenzimidazole) exclusively participates on the first stage of the reaction mediated by the CoIICbx complex, while for the CoICbx case, it appears to have an important role along the whole process. Finally, the QTAIM charge analysis indicates that oxidation of the cobalt atom occurs in both cases; at the same time, it suggests the formation of an uncommon two-center one-electron bond in the CoIICbx case. The latter was confirmed by means of electron localization calculations, which resulted in a larger electron count at the Co-C interatomic region for the CoICbx case upon comparison with its CoIICbx counterpart.

When global and local molecular descriptors are more than the sum of its parts: Simple, But Not Simpler?

In this report, we introduce a set of aggregation operators (AOs) to calculate global and local (group and atom type) molecular descriptors (MDs) as a generalization of the classical approach of molecular encoding using the sum of the atomic (or fragment) contributions. These AOs are implemented in a new and free software denominated MD-LOVIs ( http://tomocomd.com/md-lovis ), which allows for the calculation of MDs from atomic weights vector and LOVIs (local vertex invariants). This software was developed in Java programming language and employed the Chemical Development Kit (CDK) library for handling chemical structures and the calculation of atomic weights. An analysis of the complexities of the algorithms presented herein demonstrates that these aspects were efficiently implemented. The calculation speed experiments show that the MD-LOVIs software has satisfactory behavior when compared to software such as Padel, CDKDescriptor, DRAGON and Bluecal software. Shannon's entropy (SE)-based variability studies demonstrate that MD-LOVIs yields indices with greater information content when compared to those of popular academic and commercial software. A principal component analysis reveals that our approach captures chemical information orthogonal to that codified by the DRAGON, Padel and Mold2 software, as a result of the several generalizations in MD-LOVIs not used in other programs. Lastly, three QSARs were built using multiple linear regression with genetic algorithms, and the statistical parameters of these models demonstrate that the MD-LOVIs indices obtained with AOs yield better performance than those obtained when the summation operator is used exclusively. Moreover, it is also revealed that the MD-LOVIs indices yield models with comparable to superior performance when compared to other QSAR methodologies reported in the literature, despite their simplicity. The studies performed herein collectively demonstrated that MD-LOVIs software generates indices as simple as possible, but not simpler and that use of AOs enhances the diversity of the chemical information codified, which consequently improves the performance of traditional MDs.

Stability of "No-Pair Ferromagnetic" Lithium Clusters.

High-spin lithium clusters, n+1Lin (n = 2-21), have been systematically studied by using density functional theory. Although these high-spin clusters have no bonding electron pairs, they are stable with respect to isolated atoms. A set of 42 density functional theory functionals were benchmarked against CCSD(T)/cc-pVQZ results for clusters from the dimer to the hexamer. For these clusters, the strong non-additivity on the binding energy is analyzed employing a many-body energy decomposition scheme, concluding that most of the binding energy is due to a balance between the three- and four-body contributions. After a quality parameter had been defined, the LC-BP86 functional was identified as the most promising one for the description of high-spin lithium clusters. We employ the dependence of the second energy difference on cluster size to predict the formation of a higher-stability cluster.

Physico-Chemical and Structural Interpretation of Discrete Derivative Indices on N-Tuples Atoms.

This report examines the interpretation of the Graph Derivative Indices (GDIs) from three different perspectives (i.e., in structural, steric and electronic terms). It is found that the individual vertex frequencies may be expressed in terms of the geometrical and electronic reactivity of the atoms and bonds, respectively. On the other hand, it is demonstrated that the GDIs are sensitive to progressive structural modifications in terms of: size, ramifications, electronic richness, conjugation effects and molecular symmetry. Moreover, it is observed that the GDIs quantify the interaction capacity among molecules and codify information on the activation entropy. A structure property relationship study reveals that there exists a direct correspondence between the individual frequencies of atoms and Hückel's Free Valence, as well as between the atomic GDIs and the chemical shift in NMR, which collectively validates the theory that these indices codify steric and electronic information of the atoms in a molecule. Taking in consideration the regularity and coherence found in experiments performed with the GDIs, it is possible to say that GDIs possess plausible interpretation in structural and physicochemical terms.

Complex Networks Analyses of Antibiofilm Peptides: An Emerging Tool for Next-Generation Antimicrobials' Discovery.

Microbial biofilms cause several environmental and industrial issues, even affecting human health. Although they have long represented a threat due to their resistance to antibiotics, there are currently no approved antibiofilm agents for clinical treatments. The multi-functionality of antimicrobial peptides (AMPs), including their antibiofilm activity and their potential to target multiple microbes, has motivated the synthesis of AMPs and their relatives for developing antibiofilm agents for clinical purposes. Antibiofilm peptides (ABFPs) have been organized in databases that have allowed the building of prediction tools which have assisted in the discovery/design of new antibiofilm agents. However, the complex network approach has not yet been explored as an assistant tool for this aim. Herein, a kind of similarity network called the half-space proximal network (HSPN) is applied to represent/analyze the chemical space of ABFPs, aiming to identify privileged scaffolds for the development of next-generation antimicrobials that are able to target both planktonic and biofilm microbial forms. Such analyses also considered the metadata associated with the ABFPs, such as origin, other activities, targets, etc., in which the relationships were projected by multilayer networks called metadata networks (METNs). From the complex networks' mining, a reduced but informative set of 66 ABFPs was extracted, representing the original antibiofilm space. This subset contained the most central to atypical ABFPs, some of them having the desired properties for developing next-generation antimicrobials. Therefore, this subset is advisable for assisting the search for/design of both new antibiofilms and antimicrobial agents. The provided ABFP motifs list, discovered within the HSPN communities, is also useful for the same purpose.

The closed-edge structure of graphite and the effect of electrostatic charging.

The properties of graphite, and of few-layer graphene, can be strongly influenced by the edge structure of the graphene planes, but there is still much that we do not understand about the geometry and stability of these edges. We present an experimental and theoretical study of the closed edges of graphite crystals, and of the effect of an electric field on their structure. High-resolution transmission electron microscopy is used to image the edge structure of fresh graphite and of graphite that has been exposed to an electric field, which experiences a separation of the graphene layers. Computer simulations based on density functional theory are used to rationalise and quantify the preference for the formation of multiple concentric loops at the edges. A model is also presented to explain how the application of an electric field leads to the separation of the folded edges.

Higher-Order and Mixed Discrete Derivatives such as a Novel Graph- Theoretical Invariant for Generating New Molecular Descriptors.

BACKGROUND: Recently, some authors have defined new molecular descriptors (MDs) based on the use of the Graph Discrete Derivative, known as Graph Derivative Indices (GDI). This new approach about discrete derivatives over various elements from a graph takes as outset the formation of subgraphs. Previously, these definitions were extended into the chemical context (N-tuples) and interpreted in structural/physicalchemical terms as well as applied into the description of several endpoints, with good results. OBJECTIVE: A generalization of GDIs using the definitions of Higher Order and Mixed Derivative for molecular graphs is proposed as a generalization of the previous works, allowing the generation of a new family of MDs. METHODS: An extension of the previously defined GDIs is presented, and for this purpose, the concept of Higher Order Derivatives and Mixed Derivatives is introduced. These novel approaches to obtaining MDs based on the concepts of discrete derivatives (finite difference) of the molecular graphs use the elements of the hypermatrices conceived from 12 different ways (12 events) of fragmenting the molecular structures. The result of applying the higher order and mixed GDIs over any molecular structure allows finding Local Vertex Invariants (LOVIs) for atom-pairs, for atoms-pairs-pairs and so on. All new families of GDIs are implemented in a computational software denominated DIVATI (acronym for Discrete DeriVAtive Type Indices), a module of KeysFinder Framework in TOMOCOMD-CARDD system. RESULTS: QSAR modeling of the biological activity (Log 1/K) of 31 steroids reveals that the GDIs obtained using the higher order and mixed GDIs approaches yield slightly higher performance compared to previously reported approaches based on the duplex, triplex and quadruplex matrix. In fact, the statistical parameters for models obtained with the higher-order and mixed GDI method are superior to those reported in the literature by using other 0-3D QSAR methods. CONCLUSION: It can be suggested that the higher-order and mixed GDIs, appear as a promissory tool in QSAR/QSPRs, similarity/dissimilarity analysis and virtual screening studies.

Theoretical investigation of the mechanism for the reductive dehalogenation of methyl halides mediated by the CoI-based compounds cobalamin and cobaloxime.

Theoretical calculations focusing on the cleavage of the C-X bond in methyl halides (CH3X; X = Cl, Br, I) as mediated by CoI-based systems have been carried out using the hybrid functional ωB97-XD together with the basis set 6-311++G(2d,2p). A total of seven CoI-based compounds were evaluated: cob[I]alamin (CoICbl) in its base-on form and cobaloxime (CoICbx) with either no ligand or different ligands (either pyridine (PYR), tributylphosphine (TBP), dimethyl sulfide (DMS), cyclohexylisocyanide (CI), or 5,6-dimethylbenzimidazole (DMB)) at the lower axial position. For the large CoICbl system, an ONIOM scheme was employed, where the high layer was described at the DFT level and the low layer was computed using the semi-empirical method PM6. A full DFT model was employed for the CoICbx cases. An SN2-like mechanism was evaluated in all cases. The intrinsic reaction coordinate profiles suggested early transition states with activation energies of ≈ 12 kcal/mol, ≈ 10 kcal/mol, and ≈ 5 kcal/mol for C-Cl, C-Br, and C-I cleavage, respectively, which is consistent with the leaving group abilities of these halides. The evolutions of the atomic charges in and the bond orders of Co-C and C-X were computed, and the results confirmed the existence of early transition states (δBav≈ 40%), where the polarization Cδ+-Xδ- (%Ev ≈ 43%) is the determining factor in the reaction process. Finally, a comparison of all the determined parameters showed that the reaction in the DMB-CoICbx system resembles the process that occurs in the larger CoICbl, suggesting that the former system could be a reliable model for the study of reductive dehalogenation mediated by vitamin B12, which is key to the anaerobic microbiological treatment of halocarbon contaminants.

2,2,3,3'-Tetra-phenyl-7,7'-biquinoxaline.

In the crystal structure of the title compound, C(40)H(26)N(4), mol-ecules reside on crystallographic centers of inversion and are linked via C-H⋯N inter-actions about inversion centers into one-dimensional chains: longer C-H⋯π(arene) inter-actions complete the inter-molecular inter-actions.

Generalized Molecular Descriptors Derived From Event-Based Discrete Derivative.

In the present study, a generalized approach for molecular structure characterization is introduced, based on the relation frequency matrix (F) representation of the molecular graph and the subsequent calculation of the corresponding discrete derivative (finite difference) over a pair of elements (atoms). In earlier publications (22- 24), an unique event, named connected subgraphs, (based on the Kier-Hall's subgraphs) was systematically employed for the computation of the matrix F. The present report is a generalization of this notion, in which eleven additional events are introduced, classified in three categories, namely, topological (terminal paths, vertex path incidence, quantum subgraphs, walks of length k, Sach's subgraphs), fingerprints (MACCs, E-state and substructure fingerprints) and atomic contributions (Ghose and Crippen atom-types for hydrophobicity and refractivity) for F generation. The events are intended to capture diverse information by the generation or search of different kinds of substructures from the graph representation of a molecule. The discrete derivative over duplex atom relations are calculated for each event, and the resulting derivatives, local vertex invariants (LOVIs) are finally obtained. These LOVIs are subsequently employed as the basis for the calculation of global and local indices over groups of atoms (heteroatoms, halogens, methyl carbons, etc.), by using norms, means, statistics and classical algorithms as aggregator (fusion) operators. These indices were implemented in our house software DIVATI (Derivative Type Indices, a new module of TOMOCOMDCARDD system). DIVATI provides a friendly and cross-platform graphical user interface, developed in the Java programming language and is freely available at: http: //www.tomocomd.com. Factor analysis shows that the presented events are rather orthogonal and collect diverse information about the chemical structure. Finally, QSPR models were built to describe the logP and logK of 34 furylethylenes derivatives using the eleven events. Generally, the equations obtained according to these events showed high correlations, with the Sach's sub-graphs and Multiplicity events showing the best behavior in the description of logK (Q2 LOO value of 99.06%) and logP (Q2 LOO value of 98.1 %), respectively. These results show that these new eventbased indices constitute a powerful approach for chemoinformatics studies.