Tunable optical properties of isoreticular UiO-67 MOFs for photocatalysis: a theoretical study.

A theoretical study of the reported photocatalytic systems based on Zr-based MOF (UiO-67) with biphenyl-4,4'-dicarboxylic acid (bpdc) and 2,2'-bipyridine-5,5'-dicarboxylic acid (bpydc) as linkers was performed. Quantum chemical calculations were carried out to understand the optical properties of the materials and to facilitate the rational design of new UiO-67 derivatives with potentially improved features as photocatalysts under ambient conditions. Hence, the effect of the structural modifications on the optical properties was studied considering different designs based on the nature of the linkers: in 1 only the bpdc linker was considered, or the mixture 1 : 1 between bpdc and bpydc linkers (labeled as 1A). Also, substituents R, -NH2, and -SH, were included in the 1A MOF only over the bpdc linker (labeled as 1A-bpdc-R) and on both bpdc and bpydc linkers (labeled as 1A-R). Thus a family of six isoreticular UiO-67 derivatives was theoretically characterized using Density Functional Theory (DFT) calculations on the ground singlet (S0) and first excited states (singlet and triplet) using Time-Dependent Density Functional Theory (TD-DFT), multiconfigurational post-Hartree-Fock method via Complete Active Space Self-Consistent Field (CASSCF). In addition, the use of periodic DFT calculations suggest that the energy transfer (ET) channel between bpdc and bpydc linkers might generate more luminescence quenching of 1A when compare to 1. Besides, the results suggest that the 1A-R (R: -SH and NH2) can be used under ambient conditions; however, the ET exhibited by 1A, cannot take place in the same magnitude in these systems. These ET can favor the photocatalytic reduction of a potential metal ion, that can coordinate with the bpydc ligand, via LMCT transition. Consequently, the MOF might be photocatalytically active against molecules of interest (such as H2, N2, CO2, among others) with photo-reduced metal ions. These theoretical results serve as a useful tool to guide experimental efforts in the design of new photocatalytic MOF-based systems.

Co-Formulation of Recombinant Porcine IL-18 Enhances the Onset of Immune Response in a New Lawsonia intracellularis Vaccine.

Pig is one of the most consumed meats worldwide. One of the main conditions for pig production is Porcine Enteropathy caused by Lawsonia intracellularis. Among the effects of this disease is chronic mild diarrhea, which affects the weight gain of pigs, generating economic losses. Vaccines available to prevent this condition do not have the desired effect, but this limitation can be overcome using adjuvants. Pro-inflammatory cytokines, such as interleukin 18 (IL-18), can improve an immune response, reducing the immune window of protection. In this study, recombinant porcine IL-18 was produced and expressed in Escherichia coli and Pichia pastoris. The protein's biological activity was assessed in vitro and in vivo, and we determined that the P. pastoris protein had better immunostimulatory activity. A vaccine candidate against L. intracellularis, formulated with and without IL-18, was used to determine the pigs' cellular and humoral immune responses. Animals injected with the candidate vaccine co-formulated with IL-18 showed a significant increase of Th1 immune response markers and an earlier increase of antibodies than those vaccinated without the cytokine. This suggests that IL-18 acts as an immunostimulant and vaccine adjuvant to boost the immune response against the antigens, reducing the therapeutic window of recombinant protein-based vaccines.

Expanding the Knowledge of the Selective-Sensing Mechanism of Nitro Compounds by Luminescent Terbium Metal-Organic Frameworks through Multiconfigurational ab Initio Calculations.

The current research shows that the excited-state dynamics of the antenna ligand, both in the interacting system sensor/analyte and in the sensor without analyte, is a safe tool for elucidating the detection principle of the luminescent lanthanide-based metal-organic framework sensors. In this report the detection principle of the luminescence quenching mechanism in two Tb-based MOFs sensors is elucidated. The first system is a luminescent Tb-MOF [Tb(BTTA)1.5(H2O)4.5]n (H2BTTA = 2,5-bis(1H-1,2,4-triazol-1-yl) terephthalic acid) selective to nitrobenzene (NB), labeled as Tb-1. The second system is {[Tb(DPYT)(BPDC)1/2(NO3)]·H2O}n (DPYT = 2,5-di(pyridin-4-yl) terephthalic acid, BPDC = biphenyl-4,4'-dicarboxylic acid), reported as a selective chemical sensor to nitromethane (NM) in situ, labeled as Tb-2. The luminescence quenching of the MOFs is promoted by intermolecular interactions with the analytes that induce destabilization of the T1 electronic state of the linker "antenna", altering thus the sensitization pathways of the Tb atoms. This study demonstrates the value of host-guest interaction simulations and the rate constants of the radiative and nonradiative processes in understanding and elucidating the sensing mechanism in Ln-MOF sensors.

Machine learning approach to discovery of small molecules with potential inhibitory action against vasoactive metalloproteases.

With the advancement of combinatorial chemistry and big data, drug repositioning has boomed. In this sense, machine learning and artificial intelligence techniques offer a priori information to identify the most promising candidates. In this study, we combine QSAR and docking methodologies to identify compounds with potential inhibitory activity of vasoactive metalloproteases for the treatment of cardiovascular diseases. To develop this study, we used a database of 191 thermolysin inhibitor compounds, which is the largest as far as we know. First, we use Dragon's molecular descriptors (0-3D) to develop classification models using Bayesian networks (Naive Bayes) and artificial neural networks (Multilayer Perceptron). The obtained models are used for virtual screening of small molecules in the international DrugBank database. Second, docking experiments are carried out for all three enzymes using the Autodock Vina program, to identify possible interactions with the active site of human metalloproteases. As a result, high-performance artificial intelligence QSAR models are obtained for training and prediction sets. These allowed the identification of 18 compounds with potential inhibitory activity and an adequate oral bioavailability profile, which were evaluated using docking. Four of them showed high binding energies for the three enzymes, and we propose them as potential dual ACE/NEP inhibitors for the control of blood pressure. In summary, the in silico strategies used here constitute an important tool for the early identification of new antihypertensive drug candidates, with substantial savings in time and money.

Antiviral Activity of Metabolites from Peruvian Plants against SARS-CoV-2: An In Silico Approach.

(1) Background: The COVID-19 pandemic lacks treatments; for this reason, the search for potential compounds against therapeutic targets is still necessary. Bioinformatics tools have allowed the rapid in silico screening of possible new metabolite candidates from natural resources or repurposing known ones. Thus, in this work, we aimed to select phytochemical candidates from Peruvian plants with antiviral potential against three therapeutical targets of SARS-CoV-2. (2) Methods: We applied in silico technics, such as virtual screening, molecular docking, molecular dynamics simulation, and MM/GBSA estimation. (3) Results: Rutin, a compound present in Peruvian native plants, showed affinity against three targets of SARS-CoV-2. The molecular dynamics simulation demonstrated the high stability of receptor-ligand systems during the time of the simulation. Our results showed that the Mpro-Rutin system exhibited higher binding free energy than PLpro-Rutin and N-Rutin systems through MM/GBSA analysis. (4) Conclusions: Our study provides insight on natural metabolites from Peruvian plants with therapeutical potential. We found Rutin as a potential candidate with multiple pharmacological properties against SARS-CoV-2.

Theoretical Evaluation of Novel Thermolysin Inhibitors from Bacillus thermoproteolyticus. Possible Antibacterial Agents.

The search for new antibacterial agents that could decrease bacterial resistance is a subject in continuous development. Gram-negative and Gram-positive bacteria possess a group of metalloproteins belonging to the MEROPS peptidase (M4) family, which is the main virulence factor of these bacteria. In this work, we used the previous results of a computational biochemistry protocol of a series of ligands designed in silico using thermolysin as a model for the search of antihypertensive agents. Here, thermolysin from Bacillus thermoproteolyticus, a metalloprotein of the M4 family, was used to determine the most promising candidate as an antibacterial agent. Our results from docking, molecular dynamics simulation, molecular mechanics Poisson-Boltzmann (MM-PBSA) method, ligand efficiency, and ADME-Tox properties (Absorption, Distribution, Metabolism, Excretion, and Toxicity) indicate that the designed ligands were adequately oriented in the thermolysin active site. The Lig783, Lig2177, and Lig3444 compounds showed the best dynamic behavior; however, from the ADME-Tox calculated properties, Lig783 was selected as the unique antibacterial agent candidate amongst the designed ligands.

Kinetic study of removal heavy metal from aqueous solution using the synthetic aluminum silicate.

One of the problems that most affect humanity today is the wastewater discharge into different water bodies. It was estimated that more than 7 million tons of wastewater are generated worldwide and are discharged into rivers, lakes, and reservoirs. Among the most dangerous wastewaters are those from inorganic chemistry research laboratories, mainly due to heavy metals. These problems have become a highly relevant topic, and numerous researchers have tried to design wastewater treatment systems that will deal more efficiently with heavy metals elimination. In this work, the synthesis, characterization, and evaluation of hydrated aluminium silicate were performed as alternative wastewater treatment from chemistry research and teaching laboratories. The compound obtained was [Formula: see text], which was characterized by the determination of its physicochemical properties. These revealed a low density, very porous material, with low crystallinity, strong chemical resistance, a large surface area, and a high apparent ionic exchange capacity. Absorption kinetics studies of heavy metals in aqueous solutions, through more widespread models, have demonstrated that [Formula: see text] has excellent properties as absorbents of this material. The amorphous hydrated aluminium silicate achieves a decrease in the concentration of all the metal ions studied, reducing them to discharge levels permissible.

Evaluation of new antihypertensive drugs designed in silico using Thermolysin as a target.

The search for new therapies for the treatment of Arterial hypertension is a major concern in the scientific community. Here, we employ a computational biochemistry protocol to evaluate the performance of six compounds (Lig783, Lig1022, Lig1392, Lig2177, Lig3444 and Lig6199) to act as antihypertensive agents. This protocol consists of Docking experiments, efficiency calculations of ligands, molecular dynamics simulations, free energy, pharmacological and toxicological properties predictions (ADME-Tox) of the six ligands against Thermolysin. Our results show that the docked structures had an adequate orientation in the pocket of the Thermolysin enzymes, reproducing the X-ray crystal structure of Inhibitor-Thermolysin complexes in an acceptable way. The most promising candidates to act as antihypertensive agents among the series are Lig2177 and Lig3444. These compounds form the most stable ligand-Thermolysin complexes according to their binding free energy values obtained in the docking experiments as well as MM-GBSA decomposition analysis calculations. They present the lowest values of Ki, indicating that these ligands bind strongly to Thermolysin. Lig2177 was oriented in the pocket of Thermolysin in such a way that both OH of the dihydroxyl-amino groups to establish hydrogen bond interactions with Glu146 and Glu166. In the same way, Lig3444 interacts with Asp150, Glu143 and Tyr157. Additionally, Lig2177 and Lig3444 fulfill all the requirements established by Lipinski Veber and Pfizer 3/75 rules, indicating that these compounds could be safe compounds to be used as antihypertensive agents. We are confident that our computational biochemistry protocol can be used to evaluate and predict the behavior of a broad range of compounds designed in silicoagainst a protein target.

Interactions of 2-phenyl-benzotriazole xenobiotic compounds with human Cytochrome P450-CYP1A1 by means of docking, molecular dynamics simulations and MM-GBSA calculations.

2-phenyl-benzotriazole xenobiotic compounds (PBTA-4, PBTA-6, PBTA-7 and PBTA-8) that were previously isolated and identified in waters of the Yodo river, in Japan (Nukaya et al., 2001; Ohe et al., 2004; Watanabe et al., 2001) were characterized as powerful pro-mutagens. In order to predict the activation mechanism of these pro-mutagens, we designed a computational biochemistry protocol, which includes, docking experiments, molecular dynamics simulations and free energy decomposition calculations to obtain information about the interaction of 2-phenyl-benzotriazole molecules into the active center of cytochrome P450-CYP1A1 (CYP1A1). Molecular docking calculations using AutoDock Vina software shows that PBTAs are proportionally oriented in the pocket of CYP1A1, establishing π-π stacking attractive interactions between the triazole group and the Phe224, as well as, the hydrogen bonds of the terminal NH2 over the benzotriazole units with the Asn255 and Ser116 amino acids. Molecular dynamics simulations using NAMD package showed that these interactions are stable along 100.0 ns of trajectories. Into this context, free binding energy calculations employing the MM-GBSA approach, shows that some differences exists among the interaction of PBTAs with CYP1A1, regarding the solvation, electrostatic and van der Waals interaction energy components. These results suggest that PBTA molecules might be activated by CYP1A1. Thus, enhancing their mutagenicity when compared with the pro-mutagen parent species.

Study of the affinity between the protein kinase PKA and homoarginine-containing peptides derived from kemptide: Free energy perturbation (FEP) calculations.

Protein kinases (PKs) discriminate between closely related sequences that contain serine, threonine, and/or tyrosine residues. Such specificity is defined by the amino acid sequence surrounding the phosphorylatable residue, so that it is possible to identify an optimal recognition motif (ORM) for each PK. The ORM for the protein kinase A (PKA), a well-known member of the PK family, is the sequence RRX(S/T)X, where arginines at the -3 and -2 positions play a key role with respect to the primed phosphorylation site. In this work, differential affinities of PKA for the peptide substrate Kemptide (LRRASLG) and mutants that substitute the arginine residues by the unnatural peptide homoarginine were evaluated through molecular dynamics (MD) and free energy perturbation (FEP) calculations. The FEP study for the homoarginine mutants required previous elaboration of a CHARMM "arginine to homoarginine" (R2B) hybrid topology file which is available in this manuscript as Supporting Information. Mutants substituting the arginine residues by alanine, lysine, and histidine were also considered in the comparison by using the same protocol. FEP calculations allowed estimating the free energy changes from the free PKA to PKA-substrate complex (ΔΔGE→ES ) when Kemptide structure was mutated. Both ΔΔGS→ES values for homoarginine mutants were predicted with a difference below 1 kcal/mol. In addition, FEP correctly predicted that all the studied mutations decrease the catalytic efficiency of Kemptide for PKA. © 2018 Wiley Periodicals, Inc.

Calcium binding and voltage gating in Cx46 hemichannels.

The opening of connexin (Cx) hemichannels in the membrane is tightly regulated by calcium (Ca2+) and membrane voltage. Electrophysiological and atomic force microscopy experiments indicate that Ca2+ stabilizes the hemichannel closed state. However, structural data show that Ca2+ binding induces an electrostatic seal preventing ion transport without significant structural rearrangements. In agreement with the closed-state stabilization hypothesis, we found that the apparent Ca2+ sensitivity is increased as the voltage is made more negative. Moreover, the voltage and Ca2+ dependence of the channel kinetics indicate that the voltage sensor movement and Ca2+ binding are allosterically coupled. An allosteric kinetic model in which the Ca2+ decreases the energy necessary to deactivate the voltage sensor reproduces the effects of Ca2+ and voltage in Cx46 hemichannels. In agreement with the model and suggesting a conformational change that narrows the pore, Ca2+ inhibits the water flux through Cx hemichannels. We conclude that Ca2+ and voltage act allosterically to stabilize the closed conformation of Cx46 hemichannels.

Pharmacological Modulation of Proton Channel Hv1 in Cancer Therapy: Future Perspectives.

The pharmacological modulation of the immunosuppressive tumor microenvironment has emerged as a relevant component for cancer therapy. Several approaches aiming to deplete innate and adaptive suppressive populations, to circumvent the impairment in antigen presentation, and to ultimately increase the frequency of activated tumor-specific T cells are currently being explored. In this review, we address the potentiality of targeting the voltage-gated proton channel, Hv1, as a novel strategy to modulate the tumor microenvironment. The function of Hv1 in immune cells such as macrophages, neutrophils, dendritic cells, and T cells has been associated with the maintenance of NADPH oxidase activity and the generation of reactive oxygen species, which are required for the host defense against pathogens. We discuss evidence suggesting that the Hv1 proton channel could also be important for the function of these cells within the tumor microenvironment. Furthermore, as summarized here, tumor cells express Hv1 as a primary mechanism to extrude the increased amount of protons generated metabolically, thus maintaining physiologic values for the intracellular pH. Therefore, because this channel might be relevant for both tumor cells and immune cells supporting tumor growth, the pharmacological inhibition of Hv1 could be an innovative approach for cancer therapy. With that focus, we analyzed the available compounds that inhibit Hv1, highlighted the need to develop better drugs suitable for patients, and commented on the future perspectives of targeting Hv1 in the context of cancer therapy.

Study of the Differential Activity of Thrombin Inhibitors Using Docking, QSAR, Molecular Dynamics, and MM-GBSA.

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.

Decentralized peri-urban wastewater treatment technologies assessment integrating sustainability indicators.

Selection of treatment technologies without considering the environmental, economic and social factors associated with each geographical context risks the occurrence of negative impacts that were not properly foreseen, working against the sustainable performance of the technology. The principal aim of this study was to evaluate 12 technologies for decentralized treatment of domestic wastewater applicable to peri-urban communities using sustainability approaches and, at the same time, continuing a discussion about how to address a more integrated assessment of overall sustainability. For this, a set of 13 indicators that embody the environmental, economic and social approach for the overall sustainability assessment were used by means of a target plot diagram as a tool for integrating indicators that represent a holistic analysis of the technologies. The obtained results put forward different degrees of sustainability, which led to the selection of: septic tank+land infiltration; up-flow anaerobic reactor+high rate trickling filter and septic tank+anaerobic filter as the most sustainable and attractive technologies to be applied in peri-urban communities, according to the employed indicators.

Study of the affinity between the protein kinase PKA and peptide substrates derived from kemptide using molecular dynamics simulations and MM/GBSA.

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.

Insights into the interactions between maleimide derivates and GSK3ß combining molecular docking and QSAR.

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.