Structural Factors That Determine the Activity of the Xenobiotic Reductase B Enzyme from Pseudomonas putida on Nitroaromatic Compounds.

Xenobiotic reductase B (XenB) catalyzes the reduction of the aromatic ring or nitro groups of nitroaromatic compounds with methyl, amino or hydroxyl radicals. This reaction is of biotechnological interest for bioremediation, the reuse of industrial waste or the activation of prodrugs. However, the structural factors that explain the binding of XenB to different substrates are unknown. Molecular dynamics simulations and quantum mechanical calculations were performed to identify the residues involved in the formation and stabilization of the enzyme/substrate complex and to explain the use of different substrates by this enzyme. Our results show that Tyr65 and Tyr335 residues stabilize the ligands through hydrophobic interactions mediated by the aromatic rings of these aminoacids. The higher XenB activity determined with the substrates 1,3,5-trinitrobenzene and 2,4,6-trinitrotoluene is consistent with the lower energy of the highest occupied molecular orbital (LUMO) orbitals and a lower energy of the homo orbital (LUMO), which favors electrophile and nucleophilic activity, respectively. The electrostatic potential maps of these compounds suggest that the bonding requires a large hydrophobic region in the aromatic ring, which is promoted by substituents in ortho and para positions. These results are consistent with experimental data and could be used to propose point mutations that allow this enzyme to process new molecules of biotechnological interest.

Modulation of glucan-enzyme interactions by domain V in GTF-SI from Streptococcus mutans.

Glucansucrase GTF-SI from Streptococcus mutans is a multidomain enzyme that catalyzes the synthesis of glucan polymers. Domain V locates 100 Å from the catalytic site and is required for an optimal activity. Nevertheless, the mechanism governing its functional role remains elusive. In this work, homology modeling and molecular dynamics simulations were employed to examine the effect of domain V in the structure and glucan-binding ability of GTF-SI in full and truncated enzyme models. Our results showed that domain V increases the flexibility of the α4'-loop-α4″ motif near the catalytic site resulting in a higher surface for glucan association, and modulates the orientation of a growing oligosaccharide (N=8-23) in glucan-enzyme complexes towards engaging in favorable contacts throughout the protein, whereas in the truncated model the glucan protrudes randomly from domain B towards the solvent. These results are valuable to increase understanding about the functional role of domain V in GH70 glucansucrases.

Odd Loop Regions of XenA and XenB Enzymes Modulate Their Interaction with Nitro-explosives Compounds and Provide Structural Support for Their Regioselectivity.

The nitro-explosive compounds 2,4,6-trinitrotoluene, 2,4,6-trinitrophenol, and 1,2,3-trinitroglycerol are persistent environmental contaminants. The presence of different functional groups in these molecules represents a great challenge to enzymatic catalysis. The chemical variety of these three substrates is such that they do not bind and interact with catalytic residues within an enzyme with the same affinity. In this context, two Xenobiotic Reductase enzymes produced by the bacteria Pseudomonas putida can catalyze the reduction of these compounds with different affinities and regioselectivity. The structural bases that support this substrate promiscuity and catalytic preferences are unknown. Therefore, through molecular dynamics simulations and free energy calculations, we explored the structural properties driving the specific interactions of these enzymes with their substrates and cofactor. Models of Xenobiotic Reductase A and B enzymes in complex with 2,4,6-trinitrotoluene, 2,4,6-trinitrophenol, or 1,2,3-trinitroglycerol were built, and the ligand enzyme interaction was simulated by molecular dynamics. The structural analysis of the molecular dynamics simulations shows that loops 3, 5, 7, 9, 11, and 13 of Xenobiotic Reductase B, and loops 4, 5, 7, 11, 13, and 15 Xenobiotic Reductase A, are in contact with the ligands during the first stages of the molecular recognition. These loops are the most flexible regions for both enzymes; however, Xenobiotic Reductase B presents a greater range of movement and a higher number of residues interacting with the ligands. Finally, the distance between the cofactor and the different reactive groups in the substrate reflects the regioselectivity of the enzymes, and the free energy calculations are consistent with the substrate specificity of both enzymes studied. The simulation shows a stable interaction between the aromatic ring of the substrates and Xenobiotic Reductase B. In contrast, a less stable interaction with the different nitro groups of the aromatic ligands was observed. In the case of 1,2,3-trinitroglycerol, Xenobiotic Reductase B interacts more closely with the nitro groups of carbon 1 or 3, while Xenobiotic Reductase A is more selective by nitro groups of carbon 2. The obtained results suggest that the flexibility of the loops in Xenobiotic Reductase B and the presence of polar and aromatic residues present in loops 5 and 7 are fundamental to determine the affinity of the enzyme with the different substrates, and they also contribute to the proper orientation of the ligands that directs the catalytic reaction.

A QM/MM approach on the structural and stereoelectronic factors governing glycosylation by GTF-SI from Streptococcus mutans.

In this work, QM/MM calculations were employed to examine the catalytic mechanism of the retaining glucosyltransferase GTF-SI enzyme, which participates in the process of caries formation. Our goal was to characterize, with atomistic details, the mechanism of sucrose hydrolysis and the catalytic factors that modulate this reaction. Our results suggest a concerted mechanism for sucrose hydrolysis in which the first event corresponds to the glycosidic bond breakage assisted by Glu515, followed by the nucleophilic attack of Asp477, leading to the formation of the Covalent Glycosyl Enzyme (CGE) intermediate. A novel conformational itinerary of the glucosyl moiety along the reaction mechanism was identified: 2H3 → 2H3-E3 → 4C1, and the calculated energy barrier is 16.4 kcal mol-1, which is in good agreement with experimental evidence showing a major contribution coming from the glycosidic bond breakage. Our calculations also revealed that Arg475 and Asp588 play a critical role as TS-stabilizers by electrostatic and charge transfer mechanisms, respectively. This is the first report dealing with the specific features of the mechanism and catalytic residues involved in GTF-SI hydrolysis of sucrose, which is a matter of relevance in enzyme catalysis and could be valuable to aid the design of novel and specific inhibitors targeting GTF-SI.

Antioxidant activity and enzymatic of lichen substances: A study based on cyclic voltammetry and theoretical.

The antioxidant activity of nine lichen substances, including methylatrarate (1), methyl haematommate (2), lobaric acid (3), fumarprotocetraric acid (4), sphaerophorin (5), subsphaeric acid (6), diffractaic acid (7), barbatolic acid (8) and salazinic acid (9) has been determined through cyclic voltammetry. The compounds 1-4 presented slopes close to the Nernst constant of 0.059 V, indicating a 2H+/2e- relation between protons and electrons, as long as the compounds 5, 6, 7, 8, and 9 present slopes between 0.037 V and 0.032 V, indicating a 1H+/2e- relation between protons and electrons. These results show a high free radical scavenging activity by means of the release of H+, suggesting an important antioxidant capacity of these molecules. Theoretical calculations of hydrogen bond dissociation enthalpies (BDE), proton affinities (PA), and Proton Transfer (PT) mechanisms, at M06-2x/6-311+G(d,p) level complement the experimental results. Computations support that the best antioxidant activity is obtained for the molecules (3, 4, 5, 6, 7 and 8), that have a carboxylic acid group close to a phenolic hydroxyl group, through hydrogen atomic transfer (HAT) and sequential proton loss electron transfer (SPLET) mechanisms. Additional computations were performed for modelling binding affinity of the lichen substances with CYPs enzymes, mainly CYP1A2, CYP51, and CYP2C9*2 isoforms, showing strong affinity for all the compounds described in this study.

Search for Novel Potent Inhibitors of the SARS-CoV-2 Papain-like Enzyme: A Computational Biochemistry Approach.

The rapid emergence and spread of new variants of coronavirus type 2, as well as the emergence of zoonotic viruses, highlights the need for methodologies that contribute to the search for new pharmacological treatments. In the present work, we searched for new SARS-CoV-2 papain-like protease inhibitors in the PubChem database, which has more than 100 million compounds. Based on the ligand efficacy index obtained by molecular docking, 500 compounds with higher affinity than another experimentally tested inhibitor were selected. Finally, the seven compounds with ADME parameters within the acceptable range for such a drug were selected. Next, molecular dynamics simulation studies at 200 ns, ΔG calculations using molecular mechanics with generalized Born and surface solvation, and quantum mechanical calculations were performed with the selected compounds. Using this in silico protocol, seven papain-like protease inhibitors are proposed: three compounds with similar free energy (D28, D04, and D59) and three compounds with higher binding free energy (D60, D99, and D06) than the experimentally tested inhibitor, plus one compound (D24) that could bind to the ubiquitin-binding region and reduce the effect on the host immune system. The proposed compounds could be used in in vitro assays, and the described protocol could be used for smart drug design.

Exploring the structure-function relationship in a hypothetical membrane protein from a nucleotide sequence.

An educational activity is proposed that uses software for proteomic analysis and databases available for free on the Internet to estimate the structure and function of a hypothetical protein from its coding nucleotide sequence. This bioinformatics activity is recommended for integrated introductory courses that address the structure function relationship in proteins.

Biotransformation of 2,4,6-Trinitrotoluene by Pseudomonas sp. TNT3 isolated from Deception Island, Antarctica.

2,4,6-Trinitrotoluene (TNT) is a nitroaromatic explosive, highly toxic and mutagenic for organisms. In this study, we report for the first time the screening and isolation of TNT-degrading bacteria from Antarctic environmental samples with potential use as bioremediation agents. Ten TNT-degrading bacterial strains were isolated from Deception Island. Among them, Pseudomonas sp. TNT3 was selected as the best candidate since it showed the highest tolerance, growth, and TNT biotransformation capabilities. Our results showed that TNT biotransformation involves the reduction of the nitro groups. Additionally, Pseudomonas sp. TNT3 was capable of transforming 100 mg/L TNT within 48 h at 28 °C, showing higher biotransformation capability than Pseudomonas putida KT2440, a known TNT-degrading bacterium. Functional annotation of Pseudomonas sp. TNT3 genome revealed a versatile set of molecular functions involved in xenobiotic degradation pathways. Two putative xenobiotic reductases (XenA_TNT3 and XenB_TNT3) were identified by means of homology searches and phylogenetic relationships. These enzymes were also characterized at molecular level using homology modelling and molecular dynamics simulations. Both enzymes share different levels of sequence similarity with other previously described TNT-degrading enzymes and with their closest potential homologues in databases.

Local conditions for global stability in the space of codons of the genetic code.

The polar requirement is an attribute of amino acids that is a major determinant of the structure and function of the proteins, and it plays a role in the flexibility and robustness of the genetic code. The viability of an organism depends on flexibility, which allows the exploration of new functions. However, robustness is necessary to protect the organism from deleterious changes derived from misreading errors and single-point mutations. Compared with random codes, the standard genetic code is one of the most robust against such errors. Here, using analytical and numerical calculations and the set of amino acid-encoding codons, we have proposed some local conditions that are necessary for the optimal robustness of the genetic code, and we explored the association between the local conditions and the robustness. The localness of the proposed conditions and the underlying evolutionary mechanism, which begins with a random code and progresses toward more efficient codes (e.g., the standard code), might be biologically plausible.

Probable relationship between partitions of the set of codons and the origin of the genetic code.

Here we study the distribution of randomly generated partitions of the set of amino acid-coding codons. Some results are an application from a previous work, about the Stirling numbers of the second kind and triplet codes, both to the cases of triplet codes having four stop codons, as in mammalian mitochondrial genetic code, and hypothetical doublet codes. Extending previous results, in this work it is found that the most probable number of blocks of synonymous codons, in a genetic code, is similar to the number of amino acids when there are four stop codons, as well as it could be for a primigenious doublet code. Also it is studied the integer partitions associated to patterns of synonymous codons and it is shown, for the canonical code, that the standard deviation inside an integer partition is one of the most probable. We think that, in some early epoch, the genetic code might have had a maximum of the disorder or entropy, independent of the assignment between codons and amino acids, reaching a state similar to "code freeze" proposed by Francis Crick. In later stages, maybe deterministic rules have reassigned codons to amino acids, forming the natural codes, such as the canonical code, but keeping the numerical features describing the set partitions and the integer partitions, like a "fossil numbers"; both kinds of partitions about the set of amino acid-coding codons.

The most probable number of blocks for the partitions of the set of codons could have determined the number of standard amino acids.

Given a genetic code formed by 64 codons, we calculate the number of partitions of the set of encoding amino acid codons. When there are 0-3 stop codons, the results indicate that the most probable number of partitions is 19 and/or 20. Then, assuming that in the early evolution the genetic code could have had random variations, we suggest that the most probable number of partitions of the set of encoding amino acid codons determined the actual number 20 of standard amino acids.