Substitution Models of Protein Evolution with Selection on Enzymatic Activity.

Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

Dual action of benzaldehydes: Inhibiting quorum sensing and enhancing antibiotic efficacy for controlling Pseudomonas aeruginosa biofilms.

Quorum sensing (QS) has a central role in biofilm lifestyle and antimicrobial resistance, and disrupting these signaling pathways is a promising strategy to control bacterial pathogenicity and virulence. In this study, the efficacy of three structurally related benzaldehydes (4-hydroxybenzaldehyde, 4-hydroxy-3-methoxybenzaldehyde (vanillin) and 4-hydroxy-3,5-dimethoxybenzaldehyde (syringaldehyde)) in disrupting the las and pqs systems of Pseudomonas aeruginosa was investigated using bioreporter strains and computational simulations. Additionally, these benzaldehydes were combined with tobramycin and ciprofloxacin antibiotics to evaluate their ability to increase antibiotic efficacy in preventing and eradicating P. aeruginosa biofilms. To this end, the total biomass, metabolic activity and culturability of the biofilm cells were determined. In vitro assays results indicated that the aromatic aldehydes have potential to inhibit the las and pqs systems by > 80 %. Molecular docking studies supported these findings, revealing the aldehydes binding in the same pocket as the natural ligands or receptor proteins (LasR, PQSA, PQSE, PQSR). Benzaldehydes were shown to act as virulence factor attenuators, with vanillin achieving a 48 % reduction in pyocyanin production. The benzaldehyde-tobramycin combination led not only to a 60 % reduction in biomass production but also to a 90 % reduction in the metabolic activity of established biofilms. A similar result was observed when benzaldehydes were combined with ciprofloxacin. 4-Hydroxybenzaldehyde demonstrated relevant action in increasing biofilm susceptibility to ciprofloxacin, resulting in a 65 % reduction in biomass. This study discloses, for the first time, that the benzaldehydes studied are potent QS inhibitors and also enhancers of antibiotics antibiofilm activity against P. aeruginosa.

Montelukast and cefoperazone act as antiquorum sensing and antibiofilm agents against Pseudomonas aeruginosa.

AIMS: Drug repurposing is an attractive strategy to control biofilm-related infectious diseases. In this study, two drugs (montelukast and cefoperazone) with well-established therapeutic applications were tested on Pseudomonas aeruginosa quorum sensing (QS) inhibition and biofilm control. METHODS AND RESULTS: The activity of montelukast and cefoperazone was evaluated for Pqs signal inhibition, pyocyanin synthesis, and prevention and eradication of Ps. aeruginosa biofilms. Cefoperazone inhibited the Pqs system by hindering the production of the autoinducer molecules 2-heptyl-4-hydroxyquinoline (HHQ) and 2-heptyl-3-hydroxy-4(1H)-quinolone (the Pseudomonas quinolone signal or PQS), corroborating in silico results. Pseudomonas aeruginosa pyocyanin production was reduced by 50%. The combination of the antibiotics cefoperazone and ciprofloxacin was synergistic for Ps. aeruginosa biofilm control. On the other hand, montelukast had no relevant effects on the inhibition of the Pqs system and against Ps. aeruginosa biofilm. CONCLUSION: This study provides for the first time strong evidence that cefoperazone interacts with the Pqs system, hindering the formation of the autoinducer molecules HHQ and PQS, reducing Ps. aeruginosa pathogenicity and virulence. Cefoperazone demonstrated a potential to be used in combination with less effective antibiotics (e.g. ciprofloxacin) to potentiate the biofilm control action.

Heterologous expression and structure prediction of a xylanase identified from a compost metagenomic library.

Xylanases are key biocatalysts in the degradation of the ß-1,4-glycosidic linkages in the xylan backbone of hemicellulose. These enzymes are potentially applied in a wide range of bioprocessing industries under harsh conditions. Metagenomics has emerged as powerful tools for the bioprospection and discovery of interesting bioactive molecules from extreme ecosystems with unique features, such as high temperatures. In this study, an innovative combination of function-driven screening of a compost metagenomic library and automatic extraction of halo areas with in-house MATLAB functions resulted in the identification of a promising clone with xylanase activity (LP4). The LP4 clone proved to be an effective xylanase producer under submerged fermentation conditions. Sequence and phylogenetic analyses revealed that the xylanase, Xyl4, corresponded to an endo-1,4-ß-xylanase belonging to glycosyl hydrolase family 10 (GH10). When xyl4 was expressed in Escherichia coli BL21(DE3), the enzyme activity increased about 2-fold compared to the LP4 clone. To get insight on the interaction of the enzyme with the substrate and establish possible strategies to improve its activity, the structure of Xyl4 was predicted, refined, and docked with xylohexaose. Our data unveiled, for the first time, the relevance of the amino acids Glu133 and Glu238 for catalysis, and a close inspection of the catalytic site suggested that the replacement of Phe316 by a bulkier Trp may improve Xyl4 activity. Our current findings contribute to enhancing the catalytic performance of Xyl4 towards industrial applications. KEY POINTS: • A GH10 endo-1,4-ß-xylanase (Xyl4) was isolated from a compost metagenomic library • MATLAB's in-house functions were developed to identify the xylanase-producing clones • Computational analysis showed that Glu133 and Glu238 are crucial residues for catalysis.

Crystal structures of Streptomyces tsukubaensis sigma factor SigG1 and anti-sigma RsfG.

Transcription of specific genes in bacteria under environmental stress is frequently initiated by extracytoplasmic function (ECF) σ factors. ECFs σ factors harbour two conserved domains, σ2 and σ4, for transcription initiation by recognition of the promoter region and recruitment of RNA polymerase (RNAP). The crystal structure of Streptomyces tsukubaensis SigG1, an ECF56-family σ factor, was determined revealing σ2, σ4 and the additional carboxi-terminal domain SnoaL_2 tightly packed in a compact conformation. The structure of anti-sigma RsfG was also determined by X-ray crystallography and shows a rare ß-barrel fold. Analysis of the metal binding motifs inside the protein barrel are consistent with Fe(III) binding, which is in agreement with previous findings that the Streptomyces tsukubaensis ECF56 SigG1-RsfG system is involved in metal-ion homeostasis.

TargIDe: a machine-learning workflow for target identification of molecules with antibiofilm activity against Pseudomonas aeruginosa.

Bacterial biofilms are a source of infectious human diseases and are heavily linked to antibiotic resistance. Pseudomonas aeruginosa is a multidrug-resistant bacterium widely present and implicated in several hospital-acquired infections. Over the last years, the development of new drugs able to inhibit Pseudomonas aeruginosa by interfering with its ability to form biofilms has become a promising strategy in drug discovery. Identifying molecules able to interfere with biofilm formation is difficult, but further developing these molecules by rationally improving their activity is particularly challenging, as it requires knowledge of the specific protein target that is inhibited. This work describes the development of a machine learning multitechnique consensus workflow to predict the protein targets of molecules with confirmed inhibitory activity against biofilm formation by Pseudomonas aeruginosa. It uses a specialized database containing all the known targets implicated in biofilm formation by Pseudomonas aeruginosa. The experimentally confirmed inhibitors available on ChEMBL, together with chemical descriptors, were used as the input features for a combination of nine different classification models, yielding a consensus method to predict the most likely target of a ligand. The implemented algorithm is freely available at https://github.com/BioSIM-Research-Group/TargIDe under licence GNU General Public Licence (GPL) version 3 and can easily be improved as more data become available.

The Buzz on Insecticides: A Review of Uses, Molecular Structures, Targets, Adverse Effects, and Alternatives.

Insecticides play a critical role in controlling the spread of insect-borne diseases and preserving crop health. These chemical substances are specifically formulated to kill or manage insect populations. Over the years, various types of insecticides have been developed, including organophosphates, carbamates, pyrethroids, and neonicotinoids, each with unique modes of action, physiological targets, and efficacy. Despite the advantages that insecticides offer, it is imperative to recognize the potential consequences on non-target species, the environment, and human health. It is therefore crucial to follow recommended label instructions and employ integrated pest management practices for the judicious use of insecticides. This review article provides an in-depth examination of the various types of insecticides, including their modes of action, physiological targets, environmental and human health impacts, and alternatives. The aim is to furnish a comprehensive overview of insecticides and to emphasize the significance of responsible and sustainable utilization.

Selection of Bis-Indolyl Pyridines and Triphenylamines as New Inhibitors of SARS-CoV-2 Cellular Entry by Modulating the Spike Protein/ACE2 Interfaces.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the infectious agent that has caused the current coronavirus disease (COVID) pandemic. Viral infection relies on the viral S (spike) protein/cellular receptor ACE2 interaction. Disrupting this interaction would lead to early blockage of viral replication. To identify chemical tools to further study these functional interfaces, 139,146 compounds from different chemical libraries were screened through an S/ACE2 in silico virtual molecular model. The best compounds were selected for further characterization using both cellular and biochemical approaches, reiterating SARS-CoV-2 entry and the S/ACE2 interaction. We report here two selected hits, bis-indolyl pyridine AB-00011778 and triphenylamine AB-00047476. Both of these compounds can block the infectivity of lentiviral vectors pseudotyped with the SARS-CoV-2 S protein as well as wild-type and circulating variant SARS-CoV-2 strains in various human cell lines, including pulmonary cells naturally susceptible to infection. AlphaLISA and biolayer interferometry confirmed a direct inhibitory effect of these drugs on the S/ACE2 association. A specific study of the AB-00011778 inhibitory properties showed that this drug inhibits viral replication with a 50% effective concentration (EC50) between 0.1 and 0.5 µM depending on the cell lines. Molecular docking calculations of the interaction parameters of the molecules within the S/ACE2 complex from both wild-type and circulating variants of the virus showed that the molecules may target multiple sites within the S/ACE2 interface. Our work indicates that AB-00011778 constitutes a good tool for modulating this interface and a strong lead compound for further therapeutic purposes.

The Catalytic Mechanism of Pdx2 Glutaminase Driven by a Cys-His-Glu Triad: A Computational Study.

The catalytic mechanism of Pdx2 was studied with atomic detail employing the computational ONIOM hybrid QM/MM methodology. Pdx2 employs a Cys-His-Glu catalytic triad to deaminate glutamine to glutamate and ammonia - the source of the nitrogen of pyridoxal 5'-phosphate (PLP). This enzyme is, therefore, a rate-limiting step in the PLP biosynthetic pathway of Malaria and Tuberculosis pathogens that rely on this mechanism to obtain PLP. For this reason, Pdx2 is considered a novel and promising drug target to treat these diseases. The results obtained show that the catalytic mechanism of Pdx2 occurs in six steps that can be divided into four stages: (i) activation of Cys87 , (ii) deamination of glutamine with the formation of the glutamyl-thioester intermediate, (iii) hydrolysis of the formed intermediate, and (iv) enzymatic turnover. The kinetic data available in the literature (19.1-19.5 kcal mol-1 ) agree very well with the calculated free energy barrier of the hydrolytic step (18.2 kcal.mol-11 ), which is the rate-limiting step of the catalytic process when substrate is readily available in the active site. This catalytic mechanism differs from other known amidases in three main points: i) it requires the activation of the nucleophile Cys87 to a thiolate; ii) the hydrolysis occurs in a single step and therefore does not require the formation of a second tetrahedral reaction intermediate, as it is proposed, and iii) Glu198 does not have a direct role in the catalytic process. Together, these results can be used for the synthesis of new transition state analogue inhibitors capable of inhibiting Pdx2 and impair diseases like Malaria and Tuberculosis.

New insights into the catalytic mechanism of the SARS-CoV-2 main protease: an ONIOM QM/MM approach.

SARS-CoV-2 Mpro, also known as the main protease or 3C-like protease, is a key enzyme involved in the replication process of the virus that is causing the COVID-19 pandemic. It is also the most promising antiviral drug target targeting SARS-CoV-2 virus. In this work, the catalytic mechanism of Mpro was studied using the full model of the enzyme and a computational QM/MM methodology with a 69/72-atoms QM region treated at DLPNO-CCSD(T)/CBS//B3LYP/6-31G(d,p):AMBER level and including the catalytic important oxyanion-hole residues. The transition state of each step was fully characterized and described together with the related reactants and products. The rate-limiting step of the catalytic process is the hydrolysis of the thioester-enzyme adduct, and the calculated barrier closely agrees with the available kinetic data. The calculated Gibbs free energy profile, together with the full atomistic detail of the structures involved in catalysis, can now serve as valuable models for the rational drug design of transition state analogs as new inhibitors targeting the SARS-CoV-2 virus.

Special Issue on "Enzymes as Biocatalysts: Current Research Trends and Applications".

Enzymes are able to catalyze a wide diversity of chemical reactions in nature, and they do it at an amazing level [...].

A Molecular Mechanics Energy Partitioning Software for Biomolecular Systems.

The partitioning of the molecular mechanics (MM) energy in calculations involving biomolecular systems is important to identify the source of major stabilizing interactions, e.g., in ligand-protein interactions, or to identify residues with considerable contributions in hybrid multiscale calculations, i.e., quantum mechanics/molecular mechanics (QM/MM). Here, we describe Energy Split, a software program to calculate MM energy partitioning considering the AMBER Hamiltonian and parameters. Energy Split includes a graphical interface plugin for VMD to facilitate the selection of atoms and molecules belonging to each part of the system. Energy Split is freely available at or can be easily installed through the VMD Store.

Development of Neuropeptide Y and Cell-Penetrating Peptide MAP Adsorbed onto Lipid Nanoparticle Surface.

Functionalization of nanoparticles surfaces have been widely used to improve diagnostic and therapeutic biological outcome. Several methods can be applied to modify nanoparticle surface; however, in this article we focus toward a simple and less time-consuming method. We applied an adsorption method on already formulated nanostructured lipid carriers (NLC) to functionalize these nanoparticles with three distinct peptides sequences. We selected a cell-penetrating peptide (CPP), a lysine modified model amphipathic peptide (Lys(N3)-MAP), CPP/drug complex, and the neuropeptide Y. The aim of this work is to evaluate the effect of several parameters such as peptide concentration, different types of NLC, different types of peptides, and incubation medium on the physicochemical proprieties of NLC and determine if adsorption occurs. The preliminary results from zeta potential analysis indicate some evidence that this method was successful in adsorbing three types of peptides onto NLC. Several non-covalent interactions appear to be involved in peptide adsorption with the possibility of three adsorption peptide hypothesis that may occur with NLC in solution. Moreover, and for the first time, in silico docking analysis demonstrated strong interaction between CPP MAP and NPY Y1 receptor with high score values when compared to standard antagonist and NPY.

Compensatory epistasis explored by molecular dynamics simulations.

A non-negligible proportion of human pathogenic variants are known to be present as wild type in at least some non-human mammalian species. The standard explanation for this finding is that molecular mechanisms of compensatory epistasis can alleviate the mutations' otherwise pathogenic effects. Examples of compensated variants have been described in the literature but the interacting residue(s) postulated to play a compensatory role have rarely been ascertained. In this study, the examination of five human X-chromosomally encoded proteins (FIX, GLA, HPRT1, NDP and OTC) allowed us to identify several candidate compensated variants. Strong evidence for a compensated/compensatory pair of amino acids in the coagulation FIXa protein (involving residues 270 and 271) was found in a variety of mammalian species. Both amino acid residues are located within the 60-loop, spatially close to the 39-loop that performs a key role in coagulation serine proteases. To understand the nature of the underlying interactions, molecular dynamics simulations were performed. The predicted conformational change in the 39-loop consequent to the Glu270Lys substitution (associated with hemophilia B) appears to impair the protein's interaction with its substrate but, importantly, such steric hindrance is largely mitigated in those proteins that carry the compensatory residue (Pro271) at the neighboring amino acid position.

An Unsual Cys-Glu-Lys Catalytic Triad is Responsible for the Catalytic Mechanism of the Nitrilase Superfamily: A QM/MM Study on Nit2.

Nitrilase 2 (Nit2) is a representative member of the nitrilase superfamily that catalyzes the hydrolysis of α-ketosuccinamate into oxaloacetate. It has been associated with the metabolism of rapidly dividing cells like cancer cells. The catalytic mechanism of Nit2 employs a catalytic triad formed by Cys191, Glu81 and Lys150. The Cys191 and Glu81 play an active role during the catalytic process while the Lys150 is shown to play only a secondary role. The results demonstrate that the catalytic mechanism of Nit2 involves four steps. The nucleophilic attack of Cys191 to the α-ketosuccinamate, the formation of two tetrahedral enzyme adducts and the hydrolysis of a thioacyl-enzyme intermediate, from which results the formation of oxaloacetate and enzymatic turnover. The rate limiting step of the catalytic process is the formation of the first tetrahedral intermediate with a calculated activation free energy of 18.4 kcal/mol, which agrees very well with the experimental kcat (17.67 kcal/mol).

Combined in silico and in vitro studies to identify novel antidiabetic flavonoids targeting glycogen phosphorylase.

Novel pharmacological strategies for the treatment of diabetic patients are now focusing on inhibiting glycogenolysis steps. In this regard, glycogen phosphorylase (GP) is a validated target for the discovery of innovative antihyperglycemic molecules. Natural products, and in particular flavonoids, have been reported as potent inhibitors of GP at the cellular level. Herein, free-energy calculations and microscale thermophoresis approaches were performed to get an in-depth assessment of the binding affinities and elucidate intermolecular interactions of several flavonoids at the inhibitor site of GP. To our knowledge, this is the first study indicating genistein, 8-prenylgenistein, apigenin, 8-prenylapigenin, 8-prenylnaringenin, galangin and valoneic acid dilactone as natural molecules with high inhibitory potency toward GP. We identified: i) the residues Phe285, Tyr613, Glu382 and/or Arg770 as the most relevant for the binding of the best flavonoids to the inhibitor site of GP, and ii) the 5-OH, 7-OH, 8-prenyl substitutions in ring A and the 4'-OH insertion in ring B to favor flavonoid binding at this site. Our results are invaluable to plan further structural modifications through organic synthesis approaches and develop more effective pharmaceuticals for Type 2 Diabetes treatment, and serve as the starting point for the exploration of food products for therapeutic usage, as well as for the development of novel bio-functional food and dietary supplements/herbal medicines.

Databases for the study of biofilms: current status and potential applications.

Biofilms play an important role in health, being associated with >80% of all microbial infections in the body and in the development of antibiotic resistance. Research in this field has continuously produced large volumes of data. Being able to handle all this information will be paramount for progress in this field. However, this places a heavy burden on the development of strategies to gather, organize and make this information available in a way that can be readily and effectively used by those requiring it. Lately, efforts towards this goal have been reported, particularly with the development of Quorumpeps, BiofOmics, BaAMPs, QSPpred, dPABBs, aBiofilm and the Biofilms Structural Database. This work reviews these databases and highlights their applicability and potential, while stressing some of the challenges for the coming years in database development and usage brought about by the use of big data and machine learning.

Perspectives on the Role of Enzymatic Biocatalysis for the Degradation of Plastic PET.

Plastics are highly durable and widely used materials. Current methodologies of plastic degradation, elimination, and recycling are flawed. In recent years, biodegradation (the usage of microorganisms for material recycling) has grown as a valid alternative to previously used methods. The evolution of bioengineering techniques and the discovery of novel microorganisms and enzymes with degradation ability have been key. One of the most produced plastics is PET, a long chain polymer of terephthalic acid (TPA) and ethylene glycol (EG) repeating monomers. Many enzymes with PET degradation activity have been discovered, characterized, and engineered in the last few years. However, classification and integrated knowledge of these enzymes are not trivial. Therefore, in this work we present a summary of currently known PET degrading enzymes, focusing on their structural and activity characteristics, and summarizing engineering efforts to improve activity. Although several high potential enzymes have been discovered, further efforts to improve activity and thermal stability are necessary.

Identification of New Potential Inhibitors of Quorum Sensing through a Specialized Multi-Level Computational Approach.

Biofilms are aggregates of microorganisms anchored to a surface and embedded in a self-produced matrix of extracellular polymeric substances and have been associated with 80% of all bacterial infections in humans. Because bacteria in biofilms are less amenable to antibiotic treatment, biofilms have been associated with developing antibiotic resistance, a problem that urges developing new therapeutic options and approaches. Interfering with quorum-sensing (QS), an important process of cell-to-cell communication by bacteria in biofilms is a promising strategy to inhibit biofilm formation and development. Here we describe and apply an in silico computational protocol for identifying novel potential inhibitors of quorum-sensing, using CviR-the quorum-sensing receptor from Chromobacterium violaceum-as a model target. This in silico approach combines protein-ligand docking (with 7 different docking programs/scoring functions), receptor-based virtual screening, molecular dynamic simulations, and free energy calculations. Particular emphasis was dedicated to optimizing the discrimination ability between active/inactive molecules in virtual screening tests using a target-specific training set. Overall, the optimized protocol was used to evaluate 66,461 molecules, including those on the ZINC/FDA-Approved database and to the Mu.Ta.Lig Virtual Chemotheca. Multiple promising compounds were identified, yielding good prospects for future experimental validation and for drug repurposing towards QS inhibition.

In Silico Identification of Possible Inhibitors for Protein Kinase B (PknB) of Mycobacterium tuberculosis.

With tuberculosis still being one of leading causes of death in the world and the emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), researchers have been seeking to find further therapeutic strategies or more specific molecular targets. PknB is one of the 11 Ser/Thr protein kinases of Mtb and is responsible for phosphorylation-mediated signaling, mainly involved in cell wall synthesis, cell division and metabolism. With the amount of structural information available and the great interest in protein kinases, PknB has become an attractive target for drug development. This work describes the optimization and application of an in silico computational protocol to find new PknB inhibitors. This multi-level computational approach combines protein-ligand docking, structure-based virtual screening, molecular dynamics simulations and free energy calculations. The optimized protocol was applied to screen a large dataset containing 129,650 molecules, obtained from the ZINC/FDA-Approved database, Mu.Ta.Lig Virtual Chemotheca and Chimiothèque Nationale. It was observed that the most promising compounds selected occupy the adenine-binding pocket in PknB, and the main interacting residues are Leu17, Val26, Tyr94 and Met155. Only one of the compounds was able to move the active site residues into an open conformation. It was also observed that the P-loop and magnesium position loops change according to the characteristics of the ligand. This protocol led to the identification of six compounds for further experimental testing while also providing additional structural information for the design of more specific and more effective derivatives.