Gaussian basis set of triple zeta valence quality for the atoms from K to Kr: Application in DFT and CCSD(T) calculations of molecular properties

Exploring the Therapeutic Potential of Petiveria alliacea L. Phytochemicals: A Computational Study on Inhibiting SARS-CoV-2's Main Protease (Mpro)

The outbreak of SARS-CoV-2, also known as the COVID-19 pandemic, is still a critical risk factor for both human life and the global economy. Although, several promising therapies have been introduced in the literature to inhibit SARS-CoV-2, most of them are synthetic drugs that may have some adverse effects on the human body. Therefore, the main objective of this study was to carry out an in-silico investigation into the medicinal properties of Petiveria alliacea L. (P. alliacea L.)-mediated phytocompounds for the treatment of SARS-CoV-2 infections since phytochemicals have fewer adverse effects compared to synthetic drugs. To explore potential phytocompounds from P. alliacea L. as candidate drug molecules, we selected the infection-causing main protease (Mpro) of SARS-CoV-2 as the receptor protein. The molecular docking analysis of these receptor proteins with the different phytocompounds of P. alliacea L. was performed using AutoDock Vina. Then, we selected the three top-ranked phytocompounds (myricitrin, engeletin, and astilbin) as the candidate drug molecules based on their highest binding affinity scores of -8.9, -8.7 and -8.3 (Kcal/mol), respectively. Then, a 100 ns molecular dynamics (MD) simulation study was performed for their complexes with Mpro using YASARA software, computed RMSD, RMSF, PCA, DCCM, MM/PBSA, and free energy landscape (FEL), and found their almost stable binding performance. In addition, biological activity, ADME/T, DFT, and drug-likeness analyses exhibited the suitable pharmacokinetics properties of the selected phytocompounds. Therefore, the results of this study might be a useful resource for formulating a safe treatment plan for SARS-CoV-2 infections after experimental validation in wet-lab and clinical trials.

Correlation consistent basis sets designed for density functional theory: Third-row atoms (Ga-Br)

The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga-Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn-Sham limit is improved relative to other basis sets where convergence to the Kohn-Sham limit is typically not observed.

Interactions between carbon nanotubes and external structures of SARS-CoV-2 using molecular docking and molecular dynamics

Molecular modeling techniques are used to describe the process of interaction between nanotubes and the main structures of the Covid-19 virus: the envelope protein, the main protease, and the Spike glycoprotein. Molecular docking studies show that the ligands have interaction characteristics capable of adsorbing the structures. Molecular dynamics simulations provide information on the mean squared deviation of atomic positions ​​between 0.5 and 3.0 Å. The Gibbs free energy model and solvent accessible surface area approaches are used. Through the results obtained through molecular dynamics simulations, it is noted that the zig-zag nanotube prefers to interact with E-pro, M-pro, and S-gly, respectively. Molecular couplings and free energy showed that the S-gly active site residues strongly interact with zigzag, chiral, and armchair nanotubes, in this order. The interactions demonstrated in this manuscript may predict some promising candidates for virus antagonists, which may be confirmed through experimental approaches.

Can We Merge the Weak and Strong Tetrel Bonds? Electronic Features of Tetrahedral Molecules Interacted with Halide Anions

Using the orbital-free quantum crystallography approach, we have disclosed the quantitative trends in electronic features for bonds of different strengths formed by tetrel (Tt) atoms in stable molecular complexes consisting of electrically neutral tetrahedral molecules and halide anions. We have revealed the role of the electrostatic and exchange-correlation components of the total one-electron static potential that are determined by the equilibrium atomic structure and by kinetic Pauli potential, which reflects the spin-dependent electron motion features of the weak and strong bonds. The gap between the extreme positions in the electrostatic and total static potentials along the line linking the Tt atom and halide anion is wide for weak bonds and narrow for strong ones. It is in very good agreement with the number of minima in the Pauli potential between the bounded atoms. This gap exponentially correlates with the exchange-correlation potential in various series with a fixed nucleophilic fragment. A criterion for categorizing the noncovalent tetrel bonds (TtB) based on the potential features is suggested.

The Distance between Minima of Electron Density and Electrostatic Potential as a Measure of Halogen Bond Strength

In this study, we present results of a detailed topological analysis of electron density (ED) of 145 halogen-bonded complexes formed by various fluorine-, chlorine-, bromine-, and iodine-containing compounds with trimethylphosphine oxide, Me3PO. To characterize the halogen bond (XB) strength, we used the complexation enthalpy, the interatomic distance between oxygen and halogen, as well as the typical set of electron density properties at the bond critical points calculated at B3LYP/jorge-ATZP level of theory. We show for the first time that it is possible to predict the XB strength based on the distance between the minima of ED and molecular electrostatic potential (ESP) along the XB path. The gap between ED and ESP minima exponentially depends on local electronic kinetic energy density at the bond critical point and tends to be a common limiting value for the strongest halogen bond.

Interactions of Co, Cu, and non-metal phthalocyanines with external structures of SARS-CoV-2 using docking and molecular dynamics

The new coronavirus, SARS-CoV-2, caused the COVID-19 pandemic, characterized by its high rate of contamination, propagation capacity, and lethality rate. In this work, we approach the use of phthalocyanines as an inhibitor of SARS-CoV-2, as they present several interactive properties of the phthalocyanines (Pc) of Cobalt (CoPc), Copper (CuPc) and without a metal group (NoPc) can interact with SARS-CoV-2, showing potential be used as filtering by adsorption on paints on walls, masks, clothes, and air conditioning filters. Molecular modeling techniques through Molecular Docking and Molecular Dynamics were used, where the target was the external structures of the virus, but specifically the envelope protein, main protease, and Spike glycoprotein proteases. Using the g_MM-GBSA module and with it, the molecular docking studies show that the ligands have interaction characteristics capable of adsorbing the structures. Molecular dynamics provided information on the root-mean-square deviation of the atomic positions provided values between 1 and 2.5. The generalized Born implicit solvation model, Gibbs free energy, and solvent accessible surface area approach were used. Among the results obtained through molecular dynamics, it was noticed that interactions occur since Pc could bind to residues of the active site of macromolecules, demonstrating good interactions; in particular with CoPc. Molecular couplings and free energy showed that S-gly active site residues interacted strongly with phthalocyanines with values ​​of - 182.443 kJ/mol (CoPc), 158.954 kJ/mol (CuPc), and - 129.963 kJ/mol (NoPc). The interactions of Pc's with SARS-CoV-2 may predict some promising candidates for antagonists to the virus, which if confirmed through experimental approaches, may contribute to resolving the global crisis of the COVID-19 pandemic.

Coupling of fluoroborylene ligands in manganese carbonyl chemistry to give a difluorodiborene ligand

A bridging difluorodiborene μ-B₂F₂ ligand has been observed in two of the three lowest energy Mn₂(BF)₂(CO)₇ structures as well as in the lowest energy Mn₂(BF)₂(CO)₆ structure.

Analyzing ZnO clusters through the density-functional theory

The potential energy surface of Zn _n O _n clusters (n = 2, 4, 6, 8) has been explored by using a simulated annealing method. For n = 2, 4, and 6, the CCSD(T)/TZP method was used as the reference, and from here it is shown that the M06-2X/TZP method gives the lowest deviations over PBE, PBE0, B3LYP, M06, and MP2 methods. Thus, with the M06-2X method we predict isomers of Zn _n O _n clusters, which coincide with some isomers reported previously. By using the atoms in molecules analysis, possible contacts between Zn and O atoms were found for all structures studied in this article. The bond paths involved in several clusters suggest that Zn _n O _n clusters can be obtained from the zincite (ZnO crystal), such an observation was confirmed for clusters with n = 2 - 9,18 and 20. The structure with n = 23 was obtained by the procedure presented here, from crystal information, which could be important to confirm experimental data delivered for n = 18 and 23.

Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential

Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K⁺ and Na⁺, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li⁺, Na⁺, K⁺, or Rb⁺) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K⁺ or two Rb⁺ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.

Manganese-centered tubular boron cluster - MnB16 (-): A new class of transition-metal molecules

We report the observation of a manganese-centered tubular boron cluster (MnB16 (-)), which is characterized by photoelectron spectroscopy and ab initio calculations. The relatively simple pattern of the photoelectron spectrum indicates the cluster to be highly symmetric. Ab initio calculations show that MnB16 (-) has a Mn-centered tubular structure with C4v symmetry due to first-order Jahn-Teller effect, while neutral MnB16 reduces to C2v symmetry due to second-order Jahn-Teller effect. In MnB16 (-), two unpaired electrons are observed, one on the Mn 3dz(2) orbital and another on the B16 tube, making it an unusual biradical. Strong covalent bonding is found between the Mn 3d orbitals and the B16 tube, which helps to stabilize the tubular structure. The current result suggests that there may exist a whole class of metal-stabilized tubular boron clusters. These metal-doped boron clusters provide a new bonding modality for transition metals, as well as a new avenue to design boron-based nanomaterials.

Characterization of B-H agostic compounds involved in the dehydrogenation of amine-boranes by group 4 metallocenes

For over a decade, amine-borane has been considered as a potential chemical hydrogen vector in the context of a search for cleaner energy sources. When catalyzed by organometallic complexes, the reaction mechanisms currently considered involve the formation of β-BH agostic intermediates. A thorough understanding of these intermediates may constitute a crucial step toward the identification of ideal catalysts. Topological approaches such as QTAIM and ELF revealed to be particularly suitable for the description of β-agostic interactions. When studying model catalysts, accurate theoretical calculations may be carried out. However, for a comparison with experimental data, calculations should also be carried out on large organo-metallic species, often including transition metals belonging to the second or the third row. In such a case, DFT methods are particularly attractive. Unfortunately, triple-ζ all electrons basis sets are not easily available for heavy transition metal elements. Thus, a subtle balance should be reached between the affordable level of calculations and the required accuracy of the electronic description of the systems. Herein we propose the use of B3LYP functional in combination with the LanL2DZ pseudopotential for the metal atom and 6-311++G(2d,2p) basis set for the other atoms, followed by a single point using the DKH2 relativistic Hamiltonian in combination with the B3LYP/DZP-DKH level, as a "minimum level of theory" leading to a consistent topological description of the interaction within the ELF and QTAIM framework, in the context of isolated (gas-phase) group 4 metallocene catalysts.

The nature of inter- and intramolecular interactions in F2OXe(…)HX (X= F, Cl, Br, I) complexes

Electronic structure of the XeOF₂ molecule and its two complexes with HX (X= F, Cl, Br, I) molecules have been studied in the gas phase using quantum chemical topology methods: topological analysis of electron localization function (ELF), electron density, ρ(r), reduced gradient of electron density |RDG(r)| in real space, and symmetry adapted perturbation theory (SAPT) in the Hilbert space. The wave function has been approximated by the MP2 and DFT methods, using APF-D, B3LYP, M062X, and B2PLYP functionals, with the dispersion correction as proposed by Grimme (GD3). For the Xe-F and Xe=O bonds in the isolated XeOF₂ molecule, the bonding ELF-localization basins have not been observed. According to the ELF results, these interactions are not of covalent nature with shared electron density. There are two stable F₂OXe^…HF complexes. The first one is stabilized by the F-H^…F and Xe^…F interactions (type I) and the second by the F-H^…O hydrogen bond (type II). The SAPT analysis confirms the electrostatic term, E_elst⁽¹⁾ and the induction energy, E_ind⁽²⁾ to be the major contributors to stabilizing both types of complexes.

Unifying General and Segmented Contracted Basis Sets. Segmented Polarization Consistent Basis Sets

We propose a method, denoted P-orthogonalization, for converting a general contracted basis set to a computationally more efficient segmented contracted basis set, while inheriting the full accuracy of the general contracted basis set. The procedure can be used for any general contracted basis set to remove the redundancies between general contracted functions in terms of primitive functions. The P-orthogonalization procedure is used to construct a segmented contracted version of the polarization consistent basis sets, which are optimized for density functional theory calculations. Benchmark calculations show that the new pcs-n basis sets provide uniform error control of the basis set incompleteness for molecular systems composed of atoms from the first three rows in the periodic table (H-Kr) and for different exchange-correlation functionals. The basis set errors at a given zeta quality level are lower than other existing basis sets, and the pcs-n basis sets are furthermore shown to be among the computationally most efficient. The pcs-n basis sets are available in qualities ranging from (unpolarized) double-zeta to pentuple zeta quality and should therefore be well suited for both routine and benchmark calculations using density functional theory methods in general.

Synthesis and reactivity of 4,5-[1,2-dicarba-closo-dodecaborano(12)]-1,3-diselenacyclopentane: opening of the icosahedron to give a zwitterionic intermediate and conversion into 7,8-dicarba-nido-undecaborate(1-)

The reaction of the 1,2-diselenolato-1,2-dicarba-closo-dodecaborane(12) dianion [1,2-(1,2-C(2)B(10)H(10))Se(2)](2-) with dichloromethane (CH(2)Cl(2) or CD(2)Cl(2)) in the presence of donor solvents gave 4,5-[1,2-dicarba-closo-dodecaborano(12)]-1,3-diselenacyclopentane, the title compound, which was characterized by X-ray structural analysis and NMR spectroscopy ((1)H, (11)B, (13)C, and (77)Se). In the presence of pyridine, opening of the icosahedron took place, and a zwitterionic intermediate was isolated and fully characterized in the solid state by X-ray diffraction and in solution by multinuclear magnetic resonance techniques. Although such types of intermediates, prior to deboronation of the ortho-carborane cage, have been proposed several times, this is first example for which the structure has been confirmed unambiguously. This intermediate possesses a nido structure and contains a 7,8-dicarba-nido-undecaborate(1-) anion and a boronium cation, the latter with two pyridine rings linked to the boron atom, which has been extruded from the cage. It was shown that this process is reversible as long as the deboronation is not complete. The formation of the intermediate is accompanied by deboronation, which leads to the 7,8-dicarba-nido-undecaborate(1-) anion. The latter was prepared independently by conventional routes from the title compound, isolated as crystalline material as the tetrabutyl ammonium salt, and characterized by X-ray structural analysis and multinuclear magnetic resonance spectroscopy ((1)H, (11)B, (13)C, and (77)Se).