In Silico Comparative Analysis of Ivermectin and Nirmatrelvir Inhibitors Interacting with the SARS-CoV-2 Main Protease.

Exploring therapeutic options is crucial in the ongoing COVID-19 pandemic caused by SARS-CoV-2. Nirmatrelvir, which is a potent inhibitor that targets the SARS-CoV-2 Mpro, shows promise as an antiviral treatment. Additionally, Ivermectin, which is a broad-spectrum antiparasitic drug, has demonstrated effectiveness against the virus in laboratory settings. However, its clinical implications are still debated. Using computational methods, such as molecular docking and 100 ns molecular dynamics simulations, we investigated how Nirmatrelvir and Ivermectin interacted with SARS-CoV-2 Mpro(A). Calculations using density functional theory were instrumental in elucidating the behavior of isolated molecules, primarily by analyzing the frontier molecular orbitals. Our analysis revealed distinct binding patterns: Nirmatrelvir formed strong interactions with amino acids, like MET49, MET165, HIS41, HIS163, HIS164, PHE140, CYS145, GLU166, and ASN142, showing stable binding, with a root-mean-square deviation (RMSD) of around 2.0 Å. On the other hand, Ivermectin interacted with THR237, THR239, LEU271, LEU272, and LEU287, displaying an RMSD of 1.87 Å, indicating enduring interactions. Both ligands stabilized Mpro(A), with Ivermectin showing stability and persistent interactions despite forming fewer hydrogen bonds. These findings offer detailed insights into how Nirmatrelvir and Ivermectin bind to the SARS-CoV-2 main protease, providing valuable information for potential therapeutic strategies against COVID-19.

Comparison of elasticity changes in the paraspinal muscles of adolescent patients with scoliosis treated with surgery and bracing.

Scoliosis is a three-dimensional spinal deformity, and paraspinal muscles play an important role as stabilizers of the spinal curve. In this prospective study, we compared elasticity changes in the paraspinal muscles of adolescent patients with scoliosis after surgery or bracing. Elasticity was measured on the concave and convex sides of the paraspinal muscles at the apex of the curve at the beginning of treatment and 6 and 12 months after treatment. Twenty-six patients with correction surgery (n = 15) or bracing (n = 11) were included. At initial evaluation, the Cobb angle was larger in the surgery group (72.3 ± 20.2° in surgery vs. 30.6 ± 5.1° in brace, p < 0.001). The estimated mean elasticity value of the paraspinal muscles was lower in the surgery group at baseline on the convex side (15.8 vs. 22.8 kPa, p = 0.037) and 6 months on both the concave (12.1 vs. 22.7 kPa, p = 0.004) and convex (13.4 vs. 23.8 kPa, p = 0.005) sides. There was a significant stiffness decrease from baseline to 6 months on the concave side in the surgery group (5.9 kPa, p = 0.025). However, the elasticity change recovered at 12 months without significant differences between the two groups.

Structural, theoretical and biological activity of mono and binuclear nickel(II) complexes with symmetrical and asymmetrical 4,6-diacetylresorcinol-dithiocarbazate ligands.

The present work reports the synthesis and a structural study of two novel dithiocarbazate, the 4,6-diacetylresorcinol-S-benzyldithiocarbazate (H3L1) and the 4,6-diacetylresorcinol-bis(S-benzyldithiocarbazate) (H4L2), and their Ni(II) complexes, [Ni(HL1)(Py)] (1) and [Ni2(L2)(PPh3)2] (2). Single crystal X-ray analyzes reveal mono and binuclear complexes and the metal centers with distorted square planar geometry. The analyses of the Hirshfeld surface and fingerprints plots revealed intermolecular contacts attributed to the H···H and C···H/H···C bonds. The Density Functional Theory (DFT), with the B3LYP functional and 6-311-G(d,p)/LanL2DZ basis sets, was employed to optimize the geometries of synthesized compounds. From the resulting geometries, the highest occupied and lowest unoccupied molecular orbital maps (HOMO-LUMO), orbital energy gap, electron localization function (ELF), electron density, natural bond orbital (NBO) analysis, and complexation of the ligands with Ni(II) were calculated supporting the experimental data. The ESI (+)-MS/MS data indicated the presence in solution of the characteristic fragmentation with the [H3L1]+ and [H4L2]+ molecular ions for the ligands. The pharmacological potential of the dithiocarbazate ligands and their Ni(II) complexes were evaluated in vitro against MDA-MB-231 human breast cancer cells. A remarkable cytotoxic activity was observed, more evident for free ligands than complexes at low concentrations; however, this latter showed a better dose-response pattern, being more attractive in terms of pharmacokinetics and therapeutic window.

Accurate acid dissociation constant (pKa) calculation for the sulfachloropyridazine and similar molecules.

Accurate calculation of the acid dissociation constant (pKa) has fundamental importance for the description of molecular systems with pharmacological activities. The search for a more appropriate procedure for its determination is always welcome and has aroused increasing interest from the scientific community. In this sense, this work presents a computational study involving the combination of ten DFT functionals (M062X, M06L, B3LYP, BLYP, PBEPBE, BP86, LC-BLYP, SPBE, CAM-B3LYP, LC-PBEPBE) and HF method, eight basis set functions (6-311G, 6-311 + G, 6-311G(d,p), 6-311 + G(d,p), 6-311+ +G(d,p), 6-311(2d,2p), 6-311+ +G(2d,2p), and aug-cc-pVDZ), and three solvation models (SMD, PCM, and CPCM) for an accurate sulfachloropyridazine (SCR) pKa determination. It was found that the smallest deviation (0.02 unit of pKa) between the current study and experimental result was achieved with the BLYP/6-311 + G(d,p)/PCM combination. Therefore, this combination was extended to calculate the pKa of six SCR similar molecules selected through the eletroshape similarity method. For all these molecules, the difference between the obtained results and experimental data ranged between 0.14 and 0.69 units of pKa. This feature suggests that the obtained combination can determine pKa with experimental precision for complexes that are formed by sulfonamide functional group (SO2NHR). Graphical Abstract A computational study involving the combination of different levels of theory, basis sets and solvation models for an accurate sulfanamide pKa determination.

Relativistic four-component potential energy curves for the lowest 23 covalent states of molecular astatine (At2).

The potential energy curves (PECs) of all covalent states of Molecular Astatine (At2) have been investigated in this work within a four-component relativistic framework using the MOLFDIR program package. The ground state was determined using multireference configuration interaction with all single and double excitations including Davidson size-extensivity correction (MRCISD+Q) whereas the 22 excited states were treated by complete open shell configuration interaction (COSCI). Spectroscopic constants (Re,ωe,ωexe,ωeye, De,Be,αe,ße,Te) are presented for all states as well as vertical excitations obtained at COSCI, MRCISD and MRCISD+Q levels. In addition, it is also presented accurate extended Rydberg analytical form for the ground state X: (1)0g+.

Theoretical Investigation on H2O2-Ng (He, Ne, Ar, Kr, Xe, and Rn) Complexes Suitable for Stereodynamics: Interactions and Thermal Chiral Rate Consequences.

Although molecular collisions of noble gases (Ng) can be theoretically used to distinguish between the enantiomers of hydrogen peroxide - H2O2 (HP), little is known about the effects of HP-Ng interactions on the chiral rate. In this work, the chiral rate as a function of temperature (CRT) between enantiomeric conformations of HP and Ng (Ng=He, Ne, Ar, Kr, Xe, and Rn) are presented at MP2(full)/aug-cc-pVTZ level of theory through a fully basis set superposition error (BSSE) corrected potential energy surface. The results show that: (a) the CRT is highly affected even at a small decrease in the height of trans-barrier; (b) its smallest values occur with Ne for all temperatures between 100 and 4,000 K; (c) that the decrease of CRT shows an inverse correlation with respect to the average valence electron energy of the Ng and (d) Ne and He may be the noble gases more suitable for study the oriented collision dynamics of HP. In addition to binding energies, the electron density ρ and its Laplacian ∇2 ρ topological analyses were also performed within the atoms in molecules (AIM) theory in order to determine the nature of the HP-Ng interactions. The results of this work provide a more complete foundation on experiments to study HP's chirality using Ng in crossed molecular beams without a light source.