Search

This study explored the inhibitory potential of gallium maltolate against severe acute respiratory syndrome coronavirus 2 and main and papain-like proteases. Computational methods, including density functional theory and molecular docking, were used to assess gallium maltolate reactivity and binding interactions. Density functional theory calculations revealed gallium maltolate’s high electron-capturing capacity, particularly around the gallium metal atom, which may contribute to their activity. Molecular docking demonstrated that gallium maltolate can form strong hydrogen bonds with key amino acid residues like glutamate-166 and cysteine-145, tightly binding to main and papain-like proteases. The binding energy and interactions of gallium maltolate were comparable to known SARS-CoV-2 inhibitors like N-[(5-methyl-1,2-oxazol-3-yl)carbonyl]-L-alanyl-L-valyl-N-{(2S,3E)-5-(benzyloxy)-5-oxo-1-[(3S)-2-oxopyrrolidin-3-yl]pent-3-en-2-yl}-L-leucinamide, indicating its potential as an antiviral agent. However, further experimental validation is required to confirm its effectiveness in inhibiting SARS-CoV-2 replication and treating COVID-19.

This study scrutinizes the complexities of designing and exploring the potential energy surfaces of systems containing more than twenty atoms with planar tetracoordinate carbons (ptCs). To tackle this issue, we utilized an established design rule to design a Naphtho [1,2-b:3,4-b′:5,6-b″:7,8-b′′′]tetrathiophene derivative computationally. This process began with substituting S atoms with CH− units, then replacing three sequential protons with two Si2+ units in the resultant polycyclic aromatic hydrocarbon polyanion. Despite not representing the global minimum, the newly designed Si8C22 system with four ptCs provided valuable insights into strategic design and potential energy surface exploration. Our results underscore the importance of employing adequate methodologies to confirm the stability of newly designed molecular structures containing planar hypercoordinate carbons.

The Michaelis–Menten mechanism, which describes the binding of a substrate to an enzyme, is a simplification of the process on a molecular scale. A more detailed model should include the binding of the substrate to precatalytic binding sites (PCBSs) prior to the transition to the catalytic site. Our work shows that the incorporation of PCBSs, in steady-state conditions, generates a Michaelis–Menten-type expression, in which the kinetic parameters KM and Vmax adopt more complex expressions than in the model without PCBSs. The equations governing reaction kinetics can be seen as generalized symmetries, relative to time translation actions over the state space of the underlying chemical system. The study of their structure and defining parameters can be interpreted as looking for invariants associated with these time evolution actions. The expression of KM decreases as the number of PCBSs increases, while Vmax reaches a minimum when the first PCBSs are incorporated into the model. To evaluate the trend of the dynamic behavior of the system, numerical simulations were performed based on schemes with different numbers of PCBSs and six conditions of kinetic constants. From these simulations, with equal kinetic constants for the formation of the Substrate/PCBS complex, it is observed that KM and Vmax are lower than those obtained with the Michaelis–Menten model. For the model with PCBSs, the Vmax reaches a minimum at one PCBS and that value is maintained for all of the systems evaluated. Since KM decreases with the number of PCBSs, the catalytic efficiency increases for enzymes fitting this model. All of these observations are consistent with the general equation obtained. This study allows us to explain, on the basis of the PCBS to KM and Vmax ratios, the effect on enzyme parameters due to mutations far from the catalytic site, at sites involved in the first enzyme/substrate interaction. In addition, it incorporates a new mechanism of enzyme activity regulation that could be fundamental to search for new activity-modulating sites or for the design of mutants with modified enzyme parameters.

The minimum energy structures of the Si3C5 and Si4C8 clusters are planar and contain planar tetracoordinate carbons (ptCs). These species have been classified, qualitatively, as global (π) and local (σ) aromatics according to the adaptive natural density partitioning (AdNDP) method, which is an orbital localization method. This work evaluates these species’ aromaticity, focusing on confirming and quantifying their global and local aromatic character. For this purpose, we use an orbital localization method based on the partitioning of the molecular space according to the topology of the electronic localization function (LOC-ELF). In addition, the magnetically induced current density is analyzed. The LOC-ELF-based analysis coincides with the AdNDP study (double aromaticity, global, and local). Moreover, the current density analysis detects global and local ring currents. The strength of the global and local current circuit is significant, involving 4n + 2 π- and σ-electrons, respectively. The latter implicates the Si-ptC-Si fragment, which would be related to the 3c-2e σ-bond detected by the orbital localization methods in this fragment.

Here, we present evidence that the D2h M2C50/2+ (M = Li-K, Be-Ca, Al-In, and Zn) species comprises planar hexacoordinate carbon (phC) structures that exhibit four covalent and two electrostatic interactions. These findings have been made possible using evolutionary methods for exploring the potential energy surface (AUTOMATON program) and the Interacting Quantum Atoms (IQA) methodology, which support the observed bonding interactions. It is worth noting, however, that these structures are not the global minimum. Nonetheless, incorporating two cyclopentadienyl anion ligands (Cp) into the CaC52+ system has enhanced the relative stability of the phC isomer. Moreover, cycloparaphenylene ([8]CPP) provides system protection and kinetic stability. These results indicate that using appropriate ligands presents a promising approach for expanding the chemistry of phC species.

Here, it is shown that the M3B12 (M = Cu-Au) clusters’ global minima consist of an elongated planar B12 fragment connected by an in-plane linear M3 fragment. This result is striking since this B12 planar structure is not favored in the bare cluster, nor when one or two metals are added. The minimum energy structures were revealed by screening the potential energy surface using genetic algorithms and density functional theory calculations. Chemical bonding analysis shows that the strong electrostatic interactions with the metal compensate for the high energy spent in the M3 and B12 fragment distortion. Furthermore, metals participate in the delocalized π-bonds, which infers an aromatic character to these species.

We computationally explore an alternative to stabilize one-dimensional (1D) silicon-lithium nanowires (NWs). The Li12Si9 Zintl phase exhibits the NW [Li6Si5]∞1, combined with Y-shaped Si4 structures. Interestingly, this NW could be assembled from the stacking of the Li6Si5 aromatic cluster. The [Li6Si5]∞1@CNT nanocomposite has been investigated with density functional theory (DFT), including molecular dynamics simulations and electronic structure calculations. We found that van der Waals interaction between Li’s and CNT’s walls is relevant for stabilizing this hybrid nanocomposite. This work suggests that nanostructured confinement (within CNTs) may be an alternative to stabilize this free NW, cleaning its properties regarding Li12Si9 solid phase, i.e., metallic character, concerning the perturbation provided by their environment in the Li12Si7 compound.

Cocoa beans contain antioxidant molecules with the potential to inhibit type 2 coronavirus (SARS-CoV-2), which causes a severe acute respiratory syndrome (COVID-19). In particular, protease. Therefore, using in silico tests, 30 molecules obtained from cocoa were evaluated. Using molecular docking and quantum mechanics calculations, the chemical properties and binding efficiency of each ligand was evaluated, which allowed the selection of 5 compounds of this series. The ability of amentoflavone, isorhoifolin, nicotiflorin, naringin and rutin to bind to the main viral protease was studied by means of free energy calculations and structural analysis performed from molecular dynamics simulations of the enzyme/inhibitor complex. Isorhoifolin and rutin stand out, presenting a more negative binding ΔG than the reference inhibitor N-[(5-methylisoxazol-3-yl)carbonyl]alanyl-l-valyl-N~1~-((1R,2Z)−4-(benzyloxy)−4-oxo-1-{[(3R)−2-oxopyrrolidin-3-yl]methyl}but-2-enyl)-L-leucinamide (N3). These results are consistent with high affinities of these molecules for the major SARS-CoV-2. The results presented in this paper are a solid starting point for future in vitro and in vivo experiments aiming to validate these molecules and /or test similar substances as inhibitors of SARS-CoV-2 protease.

The pandemic that started in Wuhan (China) in 2019 has caused a large number of deaths, and infected people around the world due to the absence of effective therapy against coronavirus 2 of the severe acute respiratory syndrome (SARS-CoV-2). Viral maturation requires the activity of the main viral protease (Mpro), so its inhibition stops the progress of the disease. To evaluate possible inhibitors, a computational model of the SARS-CoV-2 enzyme Mpro was constructed in complex with 26 synthetic ligands derived from coumarins and quinolines. Analysis of simulations of molecular dynamics and molecular docking of the models show a high affinity for the enzyme (∆Ebinding between −5.1 and 7.1 kcal mol−1). The six compounds with the highest affinity show Kd between 6.26 × 10–6 and 17.2 × 10–6, with binding affinity between −20 and −25 kcal mol−1, with ligand efficiency less than 0.3 associated with possible inhibitory candidates. In addition to the high affinity of these compounds for SARS-CoV-2 Mpro, low toxicity is expected considering the Lipinski, Veber and Pfizer rules. Therefore, this novel study provides candidate inhibitors that would allow experimental studies which can lead to the development of new treatments for SARS-CoV-2.

Recently, we reported a series of global minima whose structures consist of carbon rings decorated with heavier group 14 elements. Interestingly, these structures feature planar tetracoordinate carbons (ptCs) and result from the replacement of five or six protons (H+) from the cyclopentadienyl anion (C5H5−) or the pentalene dianion (C8H62−) by three or four E2+ dications (E = Si–Pb), respectively. The silicon derivatives of these series are the Si3C5 and Si4C8 clusters. Here we show that ptC persists in some clusters with an equivalent number of C and Si atoms, i.e., Si5C5, Si8C8, and Si9C9. In all these species, the ptC is embedded in a pentagonal C5 ring and participates in a three-center, two-electron (3c-2e) Si-ptC-Si σ-bond. Furthermore, these clusters are π-aromatic species according to chemical bonding analysis and magnetic criteria.

It was recently shown that, when tested in cellular systems, quercetin oxidized products (Qox) have significantly better antioxidant activity than quercetin (Q) itself. The main Qox identified in the experiments are either 2,5,7,3′,4′-pentahydroxy-3,4-flavandione (Fl) or its tautomer, 2-(3,4-dihydroxybenzoyl)-2,4,6-trihydroxy-3(2H)-benzofuranone (Bf). We have now performed a theoretical evaluation of different physicochemical properties using density functional theory (DFT) calculations on Q and its main Qox species. The most stable structures (for Q and Qox) were identified after a structural search on their potential energy surface. Since proton affinities (PAs) are much lower than the bond dissociation enthalpies (BDEs) of phenolic hydrogens, we consider that direct antioxidant activity in these species is mainly due to the sequential proton loss electron transfer (SPLET) mechanism. Moreover, our kinetic studies, according to transition state theory, show that Q is more favored by this mechanism. However, Qox have lower PAs than Q, suggesting that antioxidant activity by the SPLET mechanism should be a result of a balance between proclivity to transfer protons (which favors Qox) and the reaction kinetics of the conjugated base in the sequential electron transfer mechanism (which favors Q). Therefore, our results support the idea that Q is a better direct antioxidant than its oxidized derivatives due to its kinetically favored SPLET reactions. Moreover, our molecular docking calculations indicate a stabilizing interaction between either Q or Qox and the kelch-like ECH-associated protein-1 (Keap1), in the nuclear factor erythroid 2-related factor 2 (Nrf2)-binding site. This should favor the release of the Nrf2 factor, the master regulator of anti-oxidative responses, promoting the expression of the antioxidant responsive element (ARE)-dependent genes. Interestingly, the computed Keap1-metabolite interaction energy is most favored for the Bf compound, which in turn is the most stable oxidized tautomer, according to their computed energies. These results provide further support for the hypothesis that Qox species may be better indirect antioxidants than Q, reducing reactive oxygen species in animal cells by activating endogenous antioxidants.

Autistic Spectrum Disorder (ASD) is characterized by the impairment of socio-communicative skills and the presence of restricted and stereotyped behavior patterns. Recent researches have revealed the influence of mitochondrial physiology on the development of ASD. Several research groups have identified defects in respiratory complexes, coenzyme-Q10 deficiency, increased oxidative damage, decreased of superoxide dismutase (SOD2). A study on the influence of mitochondrial physiology on the development of ASD can provide new alternatives and challenges. That is why we set ourselves the general objective to initiate studies of mitochondrial physiology in Chilean children with ASD. A sample of oral mucosa was collected in a group of 12 children diagnosed with ASD and 12 children without ASD. In children with ASD, we found a significant increase in mitochondrial DNA levels. Likewise, in these children, an increase in the protein oxidation was observed. Finally, a downward trend in the expression of the HIGD2A and SOD2 genes was observed, while DRP1, FIS1, MFN1, MFN2, and OPA1 gene expression show an upward trend. The increment of mitochondrial DNA, high oxidative stress, and high expression of the MFN2 gene could help as a scanner of the mitochondrial function in children with ASD.

En este trabajo se realizó un estudio experimental y teórico sobre la N-alquilación de nitroimidazoles. Los N-alquil-2-metil-nitroimidazoles corresponden a moléculas con actividad biológica, que se obtienen por la reacción del 2-metil-5-nitroimidazol y diferentes halogenuros de alquilo. En esta reacción se observó la formación de una mezcla de productos isoméricos en diferente proporción, denominados como N-alquil-2-metil-4-nitroimidazol y N-alquil-2-metil-5-nitroimidazol, respectivamente. La reacción sugiere la formación de un equilibrio tautomérico, que genera dos sitios nucleofílicos susceptibles a un ataque electrofílico por parte del halogenuro de alquilo. La reactividad local nucleofílica del anillo nitroimidazólico se determinó utilizando índices de reactividad local, como la función de Fukui y el potencial electrostático, además de la función de localización electrónica (ELF). La función de Fukui se integra en los átomos utilizando los esquemas de partición basados en el análisis de cargas de Mulliken y orbitales naturales de enlace (NBO). Por último, se evaluaron los perfiles de reacción. Los resultados mostraron una leve diferencia en la reactividad local, sin embargo, se observan diferencias significativas en las barreras energéticas, que explican la formación de un producto isomérico sobre el otro, encontrándo así, concordancia entre los resultados experimentales y los teóricos.In this work we realized and experimental and theoretical study of the N-alkylation of nitroimidazoles. The N-alkyl-2-methyl-nitroimidazoles correspond to biologically active molecules, obtained by reaction of 2-methyl-5-nitroimidazole and different alkyl halides. This reaction showed the formation of a mixture of isomeric products in different proportions, denominated like N-alkyl-2-methyl-4-nitroimidazole and N-alkyl-2-methyl-5-nitroimidazole, respectively. The reaction suggestes the formation of a tautomeric equilibrium, which generates two nucleophilic sites susceptible to electrophilic attack by the alkyl halide. The local nucleophilic reactivity of the nitroimidazole ring is determined using local reactivity indices such as the Fukui function and the electrostatic potential, besides the electronic localization function (ELF). The Fukui function was integrated for each atom using partition schemes based on analysis of Mulliken charges and natural bond orbital (NBO). Finally the reaction profiles were assessed. The results show a minor difference in the local reactivity. Nevertheless a significant difference in energy barriers is observed explaining the formation of an isomeric product over another. These results agree quite well with the experimental data.

Non-peptidic thrombin inhibitors (TIs; 177 compounds) with diverse groups at motifs P1 (such as oxyguanidine, amidinohydrazone, amidine, amidinopiperidine), P2 (such as cyanofluorophenylacetamide, 2-(2-chloro-6-fluorophenyl)acetamide), and P3 (such as phenylethyl, arylsulfonate groups) were studied using molecular modeling to analyze their interactions with S1, S2, and S3 subsites of the thrombin binding site. Firstly, a protocol combining docking and three dimensional quantitative structure-activity relationship was performed. We described the orientations and preferred active conformations of the studied inhibitors, and derived a predictive CoMSIA model including steric, donor hydrogen bond, and acceptor hydrogen bond fields. Secondly, the dynamic behaviors of some selected TIs (compounds 26, 133, 147, 149, 162, and 177 in this manuscript) that contain different molecular features and different activities were analyzed by creating the solvated models and using molecular dynamics (MD) simulations. We used the conformational structures derived from MD to accomplish binding free energetic calculations using MM-GBSA. With this analysis, we theorized about the effect of van der Waals contacts, electrostatic interactions and solvation in the potency of TIs. In general, the contents reported in this article help to understand the physical and chemical characteristics of thrombin-inhibitor complexes.

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔG(E→ES)) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔG(E→TS)) determined by the catalytic efficiency (k(cat)/K(M)) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.

Many protein kinase (PK) inhibitors have been reported in recent years, but only a few have been approved for clinical use. The understanding of the available molecular information using computational tools is an alternative to contribute to this process. With this in mind, we studied the binding modes of 77 maleimide derivates inside the PK glycogen synthase kinase 3 beta (GSK3β) using docking experiments. We found that the orientations that these compounds adopt inside GSK3β binding site prioritize the formation of hydrogen bond (HB) interactions between the maleimide group and the residues at the hinge region (residues Val135 and Asp133), and adopt propeller-like conformations (where the maleimide is the propeller axis and the heterocyclic substituents are two slanted blades). In addition, quantitative structure-activity relationship (QSAR) models using CoMSIA methodology were constructed to explain the trend of the GSK3β inhibitory activities for the studied compounds. We found a model to explain the structure-activity relationship of non-cyclic maleimide (NCM) derivatives (54 compounds). The best CoMSIA model (training set included 44 compounds) included steric, hydrophobic, and HB donor fields and had a good Q(2) value of 0.539. It also predicted adequately the most active compounds contained in the test set. Furthermore, the analysis of the plots of the steric CoMSIA field describes the elements involved in the differential potency of the inhibitors that can be considered for the selection of suitable inhibitors.

In this work we realized and experimental and theoretical study of the N-alkylation of nitroimidazoles.The N-alkyl-2-methyl-nitroimidazoles correspond to biologically active molecules, obtained by reactionof 2-methyl-5-nitroimidazole and different alkyl halides. This reaction showed the formation of a mixtureof isomeric products in different proportions, denominated like N-alkyl-2-methyl-4-nitroimidazole andN-alkyl-2-methyl-5-nitroimidazole, respectively. The reaction suggestes the formation of a tautomericequilibrium, which generates two nucleophilic sites susceptible to electrophilic attack by the alkyl halide.The local nucleophilic reactivity of the nitroimidazole ring is determined using local reactivity indices suchas the Fukui function and the electrostatic potential, besides the electronic localization function (ELF). TheFukui function was integrated for each atom using partition schemes based on analysis of Mulliken chargesand natural bond orbital (NBO). Finally the reaction profiles were assessed. The results show a minordifference in the local reactivity. Nevertheless a significant difference in energy barriers is observed explainingthe formation of an isomeric product over another. These results agree quite well with the experimental data.