Correction: de Oliveira et al. Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation. Toxins 2024, 16, 71.

In the published publication [...].

Interaction of major saffron constituent safranal with trypsin: An experimental and computational investigation.

Trypsin is a serine protease, an important digestive enzyme that digests the proteins in the small intestine. In the present study, we have investigated the interaction of safranal, a major saffron metabolite, with trypsin using spectroscopic and molecular docking analyses. Fluorescence emission spectra of trypsin were largely affected by the inner filter effect from safranal; that's why these were corrected using the standard procedure. The corrected fluorescence spectra have shown that the safranal quenched the intrinsic fluorescence of trypsin with a blue shift in the wavelength of emission maximum, which revealed that the microenvironment of the fluorophore became more hydrophobic. There was approximately 1: 1 fair binding between them, which increased with a rise in temperature. The interaction was favored, principally, by hydrophobic forces, and there was an efficient energy transfer from the fluorophore to the safranal. Synchronous fluorescence spectra suggested that the tryptophan residues were the major ones taking part in the fluorescence quenching of trypsin. Safranal also influenced the secondary structure of trypsin and caused partial unfolding. Molecular Docking and the Molecular Dynamics simulation of the free and complexed trypsin was also carried out. Safranal formed a stable, non-covalent complex within the S2'-S5' subsite. Moreover, two nearby tyrosine residues (Tyr39 and Tyr151) stabilized safranal through π-π interactions. Additionally, the presence of safranal led to changes in the protein flexibility and compactness, which could indicate changes in the surrounding of tryptophan residues, impacting their fluorescence. Furthermore, a loss in compactness is in line with the partial unfolding observed experimentally. Thus, both experimental and computational studies were in good agreement with each other.

Revisiting the reaction pathways for phospholipid hydrolysis catalyzed by phospholipase A2 with QM/MM methods.

Secreted phospholipase A2 (sPLA2) is a Ca2+-dependent, widely distributed enzyme superfamily in almost all mammalian tissues and bacteria. It is also a critical component of the venom of nearly all snakes, as well as many invertebrate species. In non-venomous contexts, sPLA2 assumes significance in cellular signaling pathways by binding cell membranes and catalyzing ester bond hydrolysis at the sn-2 position of phospholipids. Elevated levels of GIIA sPLA2 have been detected in the synovial fluid of arthritis patients, where it exhibits a pro-inflammatory function. Consequently, identifying sPLA2 inhibitors holds promise for creating highly effective pharmaceutical treatments. Beyond arthritis, the similarities among GIIA sPLA2s offer an opportunity for developing treatments against snakebite envenoming, the deadliest neglected tropical disease. Despite decades of study, the details of PLA2 membrane-binding, substrate-binding, and reaction mechanism remain elusive, demanding a comprehensive understanding of the sPLA2 catalytic mechanism. This study explores two reaction mechanism hypotheses, involving one or two water molecules, and distinct roles for the Ca2+ cofactor. Our research focuses on the human synovial sPLA2 enzyme bound to lipid bilayers of varying phospholipid compositions, and employing adiabatic QM/MM and QM/MM MD umbrella sampling methods to energetically and geometrically characterize the structures found along both reaction pathways. Our studies demonstrate the higher frequency of productive conformations within the single-water pathway, also revealing a lower free energy barrier for hydrolyzing POPC. Furthermore, we observe that the TS of this concerted one-step reaction closely resembles transition state geometries observed in X-ray crystallography complexes featuring high-affinity transition state analogue inhibitors.

DszA Catalyzes C-S Bond Cleavage through N5-Hydroperoxyl Formation.

Due to its detrimental impact on human health and the environment, regulations demand ultralow sulfur levels on fossil fuels, in particular in diesel. However, current desulfurization techniques are expensive and cannot efficiently remove heteroaromatic sulfur compounds, which are abundant in crude oil and concentrate in the diesel fraction after distillation. Biodesulfurization via the four enzymes of the metabolic 4S pathway of the bacterium Rhodococcus erythropolis (DszA-D) is a possible solution. However, the 4S pathway needs to operate at least 500 times faster for industrial applicability, a goal currently pursued through enzyme engineering. In this work, we unveil the catalytic mechanism of the flavin monooxygenase DszA. Surprisingly, we found that this enzyme follows a recently proposed atypical mechanism that passes through the formation of an N5OOH intermediate at the re side of the cofactor, aided by a well-defined, predominantly hydrophobic O2 pocket. Besides clarifying the unusual chemical mechanism of the complex DszA enzyme, with obvious implications for understanding the puzzling chemistry of flavin-mediated catalysis, the result is crucial for the rational engineering of DszA, contributing to making biodesulfurization attractive for the oil refining industry.

Evolution of the Quinoline Scaffold for the Treatment of Leishmaniasis: A Structural Perspective.

Since the beginning of the XXI century, Leishmaniasis has been integrated into the World Health Organization's list of the 20 neglected tropical diseases, being considered a public health issue in more than 88 countries, especially in the tropics, subtropics, and the Mediterranean area. Statistically, this disease presents a world prevalence of 12 million cases worldwide, with this number being expected to increase shortly due to the 350 million people considered at risk and the 2-2.5 million new cases appearing every year. The lack of an appropriate and effective treatment against this disease has intensified the interest of many research groups to pursue the discovery and development of novel treatments in close collaboration with the WHO, which hopes to eradicate it shortly. This paper intends to highlight the quinoline scaffold's potential for developing novel antileishmanial agents and provide a set of structural guidelines to help the research groups in the medicinal chemistry field perform more direct drug discovery and development programs. Thus, this review paper presents a thorough compilation of the most recent advances in the development of new quinoline-based antileishmanial agents, with a particular focus on structure-activity relationship studies that should be considerably useful for the future of the field.

Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation.

Viper venom phospholipase A2 enzymes (vvPLA2s) and phospholipase A2-like (PLA2-like) proteins are two of the principal toxins in viper venom that are responsible for the severe myotoxic and neurotoxic effects caused by snakebite envenoming, among other pathologies. As snakebite envenoming is the deadliest neglected tropical disease, a complete understanding of these proteins' properties and their mechanisms of action is urgently needed. Therefore, we created a database comprising information on the holo-form, cofactor-bound 3D structure of 217 vvPLA2 and PLA2-like proteins in their physiologic environment, as well as 79 membrane-bound viper species from 24 genera, which we have made available to the scientific community to accelerate the development of new anti-snakebite drugs. In addition, the analysis of the sequenced, 3D structure of the database proteins reveals essential aspects of the anatomy of the proteins, their toxicity mechanisms, and the conserved binding site areas that may anchor universal interspecific inhibitors. Moreover, it pinpoints hypotheses for the molecular origin of the myotoxicity of the PLA2-like proteins. Altogether, this study provides an understanding of the diversity of these toxins and how they are conserved, and it indicates how to develop broad, interspecies, efficient small-molecule inhibitors to target the toxin's many mechanisms of action.

Structural and dynamical changes of the Streptococcus gordonii metalloregulatory ScaR protein induced by Mn2+ ion binding.

Divalent metal ions are essential micronutrients for many intercellular reactions. Maintaining their homeostasis is necessary for the survival of bacteria. In Streptococcus gordonii, one of the primary colonizers of the tooth surface, the cellular concentration of manganese ions (Mn2+) is regulated by the manganese-sensing transcriptional factor ScaR which controls the expression of proteins involved in manganese homeostasis. To resolve the molecular mechanism through which the binding of Mn2+ ions increases the binding affinity of ScaR to DNA, a variety of computational (QM and MD) and experimental (ITC, DSC, EMSA, EPR, and CD) methods were applied. The computational results showed that Mn2+ binding induces a conformational change in ScaR that primarily affects the position of the DNA binding domains and, consequently, the DNA binding affinity of the protein. In addition, experimental results revealed a 1:4 binding stoichiometry between ScaR dimer and Mn2+ ions, while the computational results showed that the binding of Mn2+ ions in the primary binding sites is sufficient to induce the observed conformational change of ScaR.

Deciphering the Catalytic Mechanism of Virginiamycin B Lyase with Multiscale Methods and Molecular Dynamics Simulations.

Due to the emergence of antibiotic resistance, the need to explore novel antibiotics and/or novel strategies to counter antibiotic resistance is of utmost importance. In this work, we explored the molecular and mechanistic details of the degradation of a streptogramin B antibiotic by virginiamycin B (Vgb) lyase of Staphylococcus aureus using classical molecular dynamics simulations and multiscale quantum mechanics/molecular mechanics methods. Our results were in line with available experimental kinetic information. Although we were able to identify a stepwise mechanism, in the wild-type enzyme, the intermediate is short-lived, showing a small barrier to decay to the product state. The impact of point mutations on the reaction was also assessed, showing not only the importance of active site residues to the reaction catalyzed by Vgb lyase but also of near positive and negative residues surrounding the active site. Using molecular dynamics simulations, we also predicted the most likely protonation state of the 3-hydroxypicolinic moiety of the antibiotic and the impact of mutants on antibiotic binding. All this information will expand our understanding of linearization reactions of cyclic antibiotics, which are crucial for the development of novel strategies that aim to tackle antibiotic resistance.

Rational Engineering of (S)-Norcoclaurine Synthase for Efficient Benzylisoquinoline Alkaloids Biosynthesis.

(S)-Norcoclaurine is synthesized in vivo through a metabolic pathway that ends with (S)-norcoclaurine synthase (NCS). The former constitutes the scaffold for the biosynthesis of all benzylisoquinoline alkaloids (BIAs), including many drugs such as the opiates morphine and codeine and the semi-synthetic opioids oxycodone, hydrocodone, and hydromorphone. Unfortunately, the only source of complex BIAs is the opium poppy, leaving the drug supply dependent on poppy crops. Therefore, the bioproduction of (S)-norcoclaurine in heterologous hosts, such as bacteria or yeast, is an intense area of research nowadays. The efficiency of (S)-norcoclaurine biosynthesis is strongly dependent on the catalytic efficiency of NCS. Therefore, we identified vital NCS rate-enhancing mutations through the rational transition-state macrodipole stabilization method at the Quantum Mechanics/Molecular Mechanics (QM/MM) level. The results are a step forward for obtaining NCS variants able to biosynthesize (S)-norcoclaurine on a large scale.

The antioxidant activity of polysaccharides: A structure-function relationship overview.

Over the last years, polysaccharides have been linked to antioxidant effects using both in vitro chemical and biological models. The reported structures, claimed to act as antioxidants, comprise chitosan, pectic polysaccharides, glucans, mannoproteins, alginates, fucoidans, and many others of all type of biological sources. The structural features linked to the antioxidant action include the polysaccharide charge, molecular weight, and the occurrence of non-carbohydrate substituents. The establishment of structure/function relationships can be, however, biased by secondary phenomena that tailor polysaccharides behavior in antioxidant systems. In this sense, this review confronts some basic concepts of polysaccharides chemistry with the current claim of carbohydrates as antioxidants. It critically discusses how the fine structure and properties of polysaccharides can define polysaccharides as antioxidants. Polysaccharides antioxidant action is highly dependent on their solubility, sugar ring structure, molecular weight, occurrence of positive or negatively charged groups, protein moieties and covalently linked phenolic compounds. However, the occurrence of phenolic compounds and protein as contaminants leads to misleading results in methodologies often used for screening and characterization purposes, as well as in vivo models. Despite falling in the concept of antioxidants, the role of polysaccharides must be well defined according with the matrices where they are involved.

Beyond the TPP+ "gold standard": a new generation mitochondrial delivery vector based on extended PN frameworks.

Mitochondrial targeting represents an attractive strategy for treating metabolic, degenerative and hyperproliferative diseases, since this organelle plays key roles in essential cellular functions. Triphenylphosphonium (TPP+) moieties - the current "gold standard" - have been widely used as mitochondrial targeting vectors for a wide range of molecular cargo. Recently, further optimisation of the TPP+ platform drew considerable interest as a way to enhance mitochondrial therapies. However, although the modification of this system appears promising, the core structure of the TPP+ moiety remains largely unchanged. Thus, this study explored the use of aminophosphonium (PN+) and phosphazenylphosphonium (PPN+) main group frameworks as novel mitochondrial delivery vectors. The PPN+ moiety was found to be a highly promising platform for this purpose, owing to its unique electronic properties and high lipophilicity. This has been demonstrated by the high mitochondrial accumulation of a PPN+-conjugated fluorophore relative to its TPP+-conjugated counterpart, and has been further supported by density functional theory and molecular dynamics calculations, highlighting the PPN+ moiety's unusual electronic properties. These results demonstrate the potential of novel phosphorus-nitrogen based frameworks as highly effective mitochondrial delivery vectors over traditional TPP+ vectors.

Unraveling the Reaction Mechanism of Russell's Viper Venom Factor X Activator: A Paradigm for the Reactivity of Zinc Metalloproteinases?

Snake venom metalloproteinases (SVMPs) are important drug targets against snakebite envenoming, the neglected tropical disease with the highest mortality worldwide. Here, we focus on Russell's viper (Daboia russelii), one of the "big four" snakes of the Indian subcontinent that, together, are responsible for ca. 50,000 fatalities annually. The "Russell's viper venom factor X activator" (RVV-X), a highly toxic metalloproteinase, activates the blood coagulation factor X (FX), leading to the prey's abnormal blood clotting and death. Given its tremendous public health impact, the WHO recognized an urgent need to develop efficient, heat-stable, and affordable-for-all small-molecule inhibitors, for which a deep understanding of the mechanisms of action of snake's principal toxins is fundamental. In this study, we determine the catalytic mechanism of RVV-X by using a density functional theory/molecular mechanics (DFT:MM) methodology to calculate its free energy profile. The results showed that the catalytic process takes place via two steps. The first step involves a nucleophilic attack by an in situ generated hydroxide ion on the substrate carbonyl, yielding an activation barrier of 17.7 kcal·mol-1, while the second step corresponds to protonation of the peptide nitrogen and peptide bond cleavage with an energy barrier of 23.1 kcal·mol-1. Our study shows a unique role played by Zn2+ in catalysis by lowering the pKa of the Zn2+-bound water molecule, enough to permit the swift formation of the hydroxide nucleophile through barrierless deprotonation by the formally much less basic Glu140. Without the Zn2+ cofactor, this step would be rate-limiting.

Brewer's Spent Yeast Cell Wall Polysaccharides as Vegan and Clean Label Additives for Mayonnaise Formulation.

Brewer's spent yeast (BSY) mannoproteins have been reported to possess thickening and emulsifying properties. The commercial interest in yeast mannoproteins might be boosted considering the consolidation of their properties supported by structure/function relationships. This work aimed to attest the use of extracted BSY mannoproteins as a clean label and vegan source of ingredients for the replacement of food additives and protein from animal sources. To achieve this, structure/function relationships were performed by isolating polysaccharides with distinct structural features from BSY, either by using alkaline extraction (mild treatment) or subcritical water extraction (SWE) using microwave technology (hard treatment), and assessment of their emulsifying properties. Alkaline extractions solubilized mostly highly branched mannoproteins (N-linked type; 75%) and glycogen (25%), while SWE solubilized mannoproteins with short mannan chains (O-linked type; 55%) and (1→4)- and (ß1→3)-linked glucans, 33 and 12%, respectively. Extracts with high protein content yielded the most stable emulsions obtained by hand shaking, while the extracts composed of short chain mannans and ß-glucans yielded the best emulsions by using ultraturrax stirring. ß-Glucans and O-linked mannoproteins were found to contribute to emulsion stability by preventing Ostwald ripening. When applied in mayonnaise model emulsions, BSY extracts presented higher stability and yet similar texture properties as the reference emulsifiers. When used in a mayonnaise formulation, the BSY extracts were also able to replace egg yolk and modified starch (E1422) at 1/3 of their concentration. This shows that BSY alkali soluble mannoproteins and subcritical water extracted ß-glucans can be used as replacers of animal protein and additives in sauces.

Therapeutic applications of snake venoms: An invaluable potential of new drug candidates.

Animal venoms and their chemical compounds have aroused both empirical and scientific attention for ages. However, there has been a significant increase in scientific investigations in recent decades, allowing the production of various formulations that are helping in the development of many important tools for biotechnological, diagnostic, or therapeutic use, both in human and animal health, as well as in plants. Venoms are composed of biomolecules and inorganic compounds that may have physiological and pharmacological activities that are not related to their principal actions (prey immobilization, digestion, and defense). Snake venom toxins, mainly enzymatic and non-enzymatic proteins, and peptides have been identified as potential prototypes for new drugs and/or models for the development of pharmacologically active structural domains for the treatment of cancer, cardiovascular diseases, neurodegenerative and autoimmune diseases, pain, and infectious-parasitic diseases. This minireview aims to provide an overview of the biotechnological potential of animal venoms, with a focus on snakes, and to introduce the reader to the fascinating world of Applied Toxinology, where animal biodiversity can be used to develop therapeutic and diagnostic applications for humans.

Catalytically Active Snake Venom PLA2 Enzymes: An Overview of Its Elusive Mechanisms of Reaction.

Snake venom-secreted phospholipase A2 (svPLA2) enzymes, both catalytically active and inactive, are a central component in envenoming. These are responsible for disrupting the cell membrane's integrity, inducing a wide range of pharmacological effects, such as the necrosis of the bitten limb, cardiorespiratory arrest, edema, and anticoagulation. Although extensively characterized, the reaction mechanisms of enzymatic svPLA2 are still to be thoroughly understood. This review presents and analyses the most plausible reaction mechanisms for svPLA2, such as the "single-water mechanism" or the "assisted-water mechanism" initially proposed for the homologous human PLA2. All of the mechanistic possibilities are characterized by a highly conserved Asp/His/water triad and a Ca2+ cofactor. The extraordinary increase in activity induced by binding to a lipid-water interface, known as "interfacial activation," critical for the PLA2s activity, is also discussed. Finally, a potential catalytic mechanism for the postulated noncatalytic PLA2-like proteins is anticipated.

Structure-Activity Relationship Studies of 9-Alkylamino-1,2,3,4-tetrahydroacridines against Leishmania (Leishmania) infantum Promastigotes.

Leishmaniasis is one of the most neglected diseases in modern times, mainly affecting people from developing countries of the tropics, subtropics and the Mediterranean basin, with approximately 350 million people considered at risk of developing this disease. The incidence of human leishmaniasis has increased over the past decades due to failing prevention and therapeutic measures-there are no vaccines and chemotherapy, which is problematic. Acridine derivatives constitute an interesting group of nitrogen-containing heterocyclic compounds associated with numerous bioactivities, with emphasis to their antileishmanial potential. The present work builds on computational studies focusing on a specific enzyme of the parasite, S-adenosylmethionine decarboxylase (AdoMet DC), with several 1,2,3,4-tetrahydro-acridines emerging as potential inhibitors, evidencing this scaffold as a promising building block for novel antileishmanial pharmaceuticals. Thus, several 1,2,3,4-tetrahydroacridine derivatives have been synthesized, their activity against Leishmania (Leishmania) infantum promastigotes evaluated and a structure-activity relationship (SAR) study was developed based on the results obtained. Even though the majority of the 1,2,3,4-tetrahydroacridines evaluated presented high levels of toxicity, the structural information gathered in this work allowed its application with another scaffold (quinoline), leading to the obtention of N1,N12-bis(7-chloroquinolin-4-yl)dodecane-1,12-diamine (12) as a promising novel antileishmanial agent (IC50 = 0.60 ± 0.11 µM, EC50 = 11.69 ± 3.96 µM and TI = 19.48).

Engineering DszC Mutants from Transition State Macrodipole Considerations and Evolutionary Sequence Analysis.

We describe an approach to identify enzyme mutants with increased turnover using the enzyme DszC as a case study. Our approach is based on recalculating the barriers of alanine mutants through single-point energy calculations at the hybrid QM/MM level in the wild-type reactant and transition state geometries. We analyze the difference in the electron density between the reactant and transition state to identify sites/residues where electrostatic interactions stabilize the transition state over the reactants. We also assess the insertion of a unit probe charge to identify positions in which the introduction of charged residues lowers the barrier.

Development of Nanoscale Graphene Oxide Models for the Adsorption of Biological Molecules.

Graphene oxide (GO), a nanomaterial with promising applications that range from water purification to enzyme immobilization, is actively present in scientific research since its discovery. GO studies with computational methodologies such as molecular dynamics are frequently reported in the literature; however, the models used often rely on approximations, such as randomly placing functional groups and the use of generalized force fields. Therefore, it is important to develop new MD models that provide a more accurate description of GO structures and their interaction with an aqueous solvent and other adsorbate molecules. In this paper, we derived new force field non-bonded parameters from linear-scaling density functional theory calculations of nanoscale GO sheets with more than 10,000 atoms through an atoms-in-molecules (AIM) partitioning scheme. The resulting GAFF2-AIM force field, derived from the bonded terms of GAFF2 parameterization, reproduces the solvent structure reported in ab initio MD simulations better than the force field nowadays widely used in the literature. Additionally, we analyzed the effect of the ionic strength of the medium and of the C/O ratio on the distribution of charges surrounding the GO sheets. Finally, we simulated the adsorption of natural amino acid molecules to a GO sheet and estimated their free energy of binding, which compared very favorably to their respective experimental values, validating the force field presented in this work.

QM/MM Study of the Reaction Mechanism of Thermophilic Glucuronoyl Esterase for Biomass Treatment.

Hydrolysis of lignocellulosic biomass, composed of a lignin-carbohydrate-complex (LCC) matrix, is critical for producing bioethanol from glucose. However, current methods for LCC processing require costly and polluting processes. The fungal Thermothelomyces thermophila glucuronoyl esterase (TtGE) is a promising thermophilic enzyme that hydrolyses LCC ester bonds. This study describes the TtGE catalytic mechanism using QM/MM methods. Two nearly-degenerate rate-determining transition states were found, with barriers of 16 and 17 kcal ⋅ mol-1 , both with a zwitterionic nature that results from a proton interplay from His346 to either the Ser213-hydroxyl or the lignin leaving group and the rehybridisation of the ester moiety of the substrate to an alkoxide. An oxyanion hole, characteristic of esterases, was provided by the conserved Arg214 through its backbone and sidechain. Our work further suggests that a mutation of Glu267 to a non-negative residue will decrease the energetic barrier in ca. -5 kcal ⋅ mol-1 , improving the catalytic rate of TtGE.

Role of Enzyme and Active Site Conformational Dynamics in the Catalysis by α-Amylase Explored with QM/MM Molecular Dynamics.

We assessed enzyme:substrate conformational dynamics and the rate-limiting glycosylation step of a human pancreatic α-amylase:maltopentose complex. Microsecond molecular dynamics simulations suggested that the distance of the catalytic Asp197 nucleophile to the anomeric carbon of the buried glucoside is responsible for most of the enzyme active site fluctuations and that both Asp197 and Asp300 interact the most with the buried glucoside unit. The buried glucoside binds either in a 4C1 chair or 2SO skew conformations, both of which can change to TS-like conformations characteristic of retaining glucosidases. Starting from four distinct enzyme:substrate complexes, umbrella sampling quantum mechanics/molecular mechanics simulations (converged within less than 1 kcal·mol-1 within a total simulation time of 1.6 ns) indicated that the reaction occurrs with a Gibbs barrier of 13.9 kcal·mol -1, in one asynchronous concerted step encompassing an acid-base reaction with Glu233 followed by a loose SN2-like nucleophilic substitution by the Asp197. The transition state is characterized by a 2H3 half-chair conformation of the buried glucoside that quickly changes to the E3 envelope conformation preceding the attack of the anomeric carbon by the Asp197 nucleophile. Thermodynamic analysis of the reaction supported that a water molecule tightly hydrogen bonded to the glycosidic oxygen of the substrate at the reactant state (∼1.6 Å) forms a short hydrogen bond with Glu233 at the transition state (∼1.7 Å) and lowers the Gibbs barrier in over 5 kcal·mol-1. The resulting Asp197-glycosyl was mostly found in the 4C1 conformation, although the more endergonic B3,O conformation was also observed. Altogether, the combination of short distances for the acid-base reaction with the Glu233 and for the nucleophilic attack by the Asp197 nucleophile and the availability of water within hydrogen bonding distance of the glycosidic oxygen provides a reliable criteria to identify reactive conformations of α-amylase complexes.