Transmembrane Protease Serine 2 Proteolytic Cleavage of the SARS-CoV-2 Spike Protein: A Mechanistic Quantum Mechanics/Molecular Mechanics Study to Inspire the Design of New Drugs To Fight the COVID-19 Pandemic.

Despite the development of vaccines against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, there is an urgent need for efficient drugs to treat infected patients. An attractive drug target is the human transmembrane protease serine 2 (TMPRSS2) because of its vital role in the viral infection mechanism of SARS-CoV-2 by activation of the virus spike protein (S protein). Having in mind that the information derived from quantum mechanics/molecular mechanics (QM/MM) studies could be an important tool in the design of transition-state (TS) analogue inhibitors, we resorted to adiabatic QM/MM calculations to determine the mechanism of the first step (acylation) of proteolytic cleavage of the S protein with atomistic details. Acylation occurred in two stages: (i) proton transfer from Ser441 to His296 concerted with the nucleophilic attack of Ser441 to the substrate's P1-Arg and (ii) proton transfer from His296 to the P1'-Ser residue concerted with the cleavage of the ArgP1-SerP1' peptide bond, with a Gibbs activation energy of 17.1 and 15.8 kcal mol-1, relative to the reactant. An oxyanion hole composed of two hydrogen bonds stabilized the rate-limiting TS by 8 kcal mol-1. An analysis of the TMPRSS2 interactions with the high-energy, short-lived tetrahedral intermediate highlighted the limitations of current clinical inhibitors and pointed out specific ways to develop higher-affinity TS analogue inhibitors. The results support the development of more efficient drugs against SARS-CoV-2 using a human target, free from resistance development.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase.

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol-1, and for the second stage, it was 7.2 kcal·mol-1, in excellent agreement with the experimental barrier (17.6 kcal·mol-1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol-1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

Structural and mechanistic aspects of S-S bonds in the thioredoxin-like family of proteins.

Disulfide bonds play a critical role in a variety of structural and mechanistic processes associated with proteins inside the cells and in the extracellular environment. The thioredoxin family of proteins like thioredoxin (Trx), glutaredoxin (Grx) and protein disulfide isomerase, are involved in the formation, transfer or isomerization of disulfide bonds through a characteristic thiol-disulfide exchange reaction. Here, we review the structural and mechanistic determinants behind the thiol-disulfide exchange reactions for the different enzyme types within this family, rationalizing the known experimental data in light of the results from computational studies. The analysis sheds new atomic-level insight into the structural and mechanistic variations that characterize the different enzymes in the family, helping to explain the associated functional diversity. Furthermore, we review here a pattern of stabilization/destabilization of the conserved active-site cysteine residues presented beforehand, which is fully consistent with the observed roles played by the thioredoxin family of enzymes.

The Catalytic Mechanism of Human Transketolase.

We have computationally determined the catalytic mechanism of human transketolase (hTK) using a cluster model approach and density functional theory calculations. We were able to determine all the relevant structures, bringing solid evidences to the proposed experimental mechanism, and to add important detail to the structure of the transition states and the energy profile associated with catalysis. Furthermore, we have established the existence of a crucial intermediate of the catalytic cycle, in agreement with experiments. The calculated data brought new insights to hTK's catalytic mechanism, providing free-energy values for the chemical reaction, as well as adding atomistic detail to the experimental mechanism.

Improving the Catalytic Power of the DszD Enzyme for the Biodesulfurization of Crude Oil and Derivatives.

The enhancement of the catalytic power of enzymes is a subject of enormous interest both for science and for industry. The latter, in particular, due to the vast applications enzymes can have in industrial processes, for instance in the desulfurization of crude oil, which is mandatory by law in many developed countries and is currently performed using costly chemical processes. In this work we sought to enhance the turnover rate of DszD from Rhodococcus erythropolis, a NADH-FMN oxidoreductase responsible for supplying FMNH2 to DszA and DszC in the biodesulfurization process of crude oil, the 4S pathway. For this purpose, we replaced the wild type spectator residue of the rate limiting step of the reduction of FMN to FMNH2 , a process catalysed by DszD and known to play an important role in the reaction energy profile. As replacements, we used all the naturally occurring amino acids, one at a time, using computational methodologies, and repeated the above-mentioned reaction with each mutant. To calculate the different free energy profiles, one for each mutated model, we applied quantum mechanics/molecular mechanics (QM/MM) methods within an ONIOM scheme. The free energy barriers obtained varied between 15.1 and 29.9 kcal mol-1 . Multiple factors contributed to the different ΔG values. The most relevant were electrostatic interactions and the induction of a favourable alignment between substrate and cofactor. These results confirm the great potential that chirurgic mutations have for increasing the catalytic power of DszD in relation to the wild type (wt) enzyme.

The Catalytic Mechanism of the Pyridoxal-5'-phosphate-Dependent Enzyme, Histidine Decarboxylase: A Computational Study.

The catalytic mechanism of histidine decarboxylase (HDC), a pyridoxal-5'-phosphate (PLP)-dependent enzyme, was studied by using a computational QM/MM approach following the scheme M06-2X/6-311++G(3df,2pd):Amber. The reaction involves two sequential steps: the decarboxylation of l-histidine and the protonation of the generated intermediate from which results histamine. The rate-limiting step is the first one (ΔG≠ =17.6 kcal mol-1 ; ΔGr =13.7 kcal mol-1 ) and agrees closely with the available experimental kcat (1.73 s-1 ), which corresponds to an activation barrier of 17.9 kcal mol-1 . In contrast, the second step is very fast (ΔG≠ =1.9 kcal mol-1 ) and exergonic (ΔGr =-33.2 kcal mol-1 ). Our results agree with the available experimental data and allow us to explain the role played by several active site residues that are considered relevant according to site-directed mutagenesis studies, namely Tyr334B, Asp273A, Lys305A, and Ser354B. These results can provide insights regarding the catalytic mechanism of other enzymes belonging to family II of PLP-dependent decarboxylases.

Periplasmic nitrate reductase and formate dehydrogenase: similar molecular architectures with very different enzymatic activities.

It is remarkable how nature has been able to construct enzymes that, despite sharing many similarities, have simple but key differences that tune them for completely different functions in living cells. Periplasmic nitrate reductase (Nap) and formate dehydrogenase (Fdh) from the DMSOr family are representative examples of this. Both enzymes share almost identical three-dimensional protein foldings and active sites, in terms of coordination number, geometry and nature of the ligands. The substrates of both enzymes (nitrate and formate) are polyatomic anions that also share similar charge and stereochemistry. In terms of the catalytic mechanism, both enzymes have a common activation mechanism (the sulfur-shift mechanism) that ensures a constant coordination number around the metal ion during the catalytic cycle. In spite of these similarities, they catalyze very different reactions: Nap abstracts an oxygen atom from nitrate releasing nitrite, whereas FdH catalyzes a hydrogen atom transfer from formate and releases carbon dioxide. In this Account, a critical analysis of structure, function, and catalytic mechanism of the molybdenum enzymes periplasmic nitrate reductase (Nap) and formate dehydrogenase (Fdh) is presented. We conclude that the main structural driving force that dictates the type of reaction, catalyzed by each enzyme, is a key difference on one active site residue that is located in the top region of the active sites of both enzymes. In both enzymes, the active site is centered on the metal ion of the cofactor (Mo in Nap and Mo or W in Fdh) that is coordinated by four sulfur atoms from two pyranopterin guanosine dinucleotide (PGD) molecules and by a sulfido. However, while in Nap there is a Cys directly coordinated to the Mo ion, in FdH there is a SeCys instead. In Fdh there is also an important His that interacts very closely with the SeCys, whereas in Nap the same position is occupied by a Met. The role of Cys in Nap and SeCys in FdH is similar in both enzymes; however, Met and His have different roles. His participates directly on catalysis, and it is therefore detrimental for the catalytic cycle of FdH. Met only participates in substrate binding. We concluded that this small but key difference dictates the type of reaction that is catalyzed by each enzyme. In addition, it allows explaining why formate can bind in the Nap active site in the same way as the natural substrate (nitrate), but the reaction becomes stalled afterward.

Calculation of distribution coefficients in the SAMPL5 challenge from atomic solvation parameters and surface areas.

In the context of SAMPL5, we submitted blind predictions of the cyclohexane/water distribution coefficient (D) for a series of 53 drug-like molecules. Our method is purely empirical and based on the additive contribution of each solute atom to the free energy of solvation in water and in cyclohexane. The contribution of each atom depends on the atom type and on the exposed surface area. Comparatively to similar methods in the literature, we used a very small set of atomic parameters: only 10 for solvation in water and 1 for solvation in cyclohexane. As a result, the method is protected from overfitting and the error in the blind predictions could be reasonably estimated. Moreover, this approach is fast: it takes only 0.5 s to predict the distribution coefficient for all 53 SAMPL5 compounds, allowing its application in virtual screening campaigns. The performance of our approach (submission 49) is modest but satisfactory in view of its efficiency: the root mean square error (RMSE) was 3.3 log D units for the 53 compounds, while the RMSE of the best performing method (using COSMO-RS) was 2.1 (submission 16). Our method is implemented as a Python script available at https://github.com/diogomart/SAMPL5-DC-surface-empirical .

Improving the Biodesulfurization of Crude Oil and Derivatives: A QM/MM Investigation of the Catalytic Mechanism of NADH-FMN Oxidoreductase (DszD).

The development of biocatalytic desulfurization strategies of petroleum and its derivatives could result in more economic alternatives than the widely used chemical desulfurization. The organism Rhodococcus erythropolis IGTS8 has been shown to metabolize organic sulfur compounds through a mechanism known as 4S pathway, which involves four enzymes (DszA, DszB, DszC, and DszD) and has been explored in biodesulfurization. Here we have applied QM/MM methods to study the catalytic mechanism of the enzyme DszD, a NADH-FMN oxidoreductase that occupies a central place on the 4S pathway by catalyzing the formation of the FMNH2 that is used by the two monooxynases in the cycle: DszA and DszC. In addition, to clarify the catalytic mechanism of this enzyme, this study analyzed in detail the role played by the active site Thr residue and of Asn and Ala enzyme mutants. The results help to explain previous experimental evidence and suggest new strategies for improving biodesulfurization through an increase in the activity of DszD.

Insights into the structural determinants of substrate specificity and activity in mouse aldehyde oxidases.

In this work, a combination of homology modeling and molecular dynamics (MD) simulations was used to investigate the factors that modulate substrate specificity and activity of the mouse AOX isoforms: mAOX1, mAOX2 (previously mAOX3l1), mAOX3 and mAOX4. The results indicate that the AOX isoform structures are highly preserved and even more conserved than the corresponding amino acid sequences. The only differences are at the protein surface and substrate-binding site region. The substrate-binding site of all isoforms consists of two regions: the active site, which is highly conserved among all isoforms, and a isoform-specific region located above. We predict that mAOX1 accepts a broader range of substrates of different shape, size and nature relative to the other isoforms. In contrast, mAOX4 appears to accept a more restricted range of substrates. Its narrow and hydrophobic binding site indicates that it only accepts small hydrophobic substrates. Although mAOX2 and mAOX3 are very similar to each other, we propose the following pairs of overlapping substrate specificities: mAOX2/mAOX4 and mAOX3/mAXO1. Based on these considerations, we propose that the catalytic activity between all isoforms should be similar but the differences observed in the binding site might influence the substrate specificity of each enzyme. These results also suggest that the presence of several AOX isoforms in mouse allows them to oxidize more efficiently a wider range of substrates. This contrasts with the same or other organisms that only express one isoform and are less efficient or incapable of oxidizing the same type of substrates.

PLP undergoes conformational changes during the course of an enzymatic reaction.

Numerous enzymes, such as the pyridoxal 5'-phosphate (PLP)-dependent enzymes, require cofactors for their activities. Using X-ray crystallography, structural snapshots of the L-serine dehydratase catalytic reaction of a bacterial PLP-dependent enzyme were determined. In the structures, the dihedral angle between the pyridine ring and the Schiff-base linkage of PLP varied from 18° to 52°. It is proposed that the organic cofactor PLP directly catalyzes reactions by active conformational changes, and the novel catalytic mechanism involving the PLP cofactor was confirmed by high-level quantum-mechanical calculations. The conformational change was essential for nucleophilic attack of the substrate on PLP, for concerted proton transfer from the substrate to the protein and for directing carbanion formation of the substrate. Over the whole catalytic cycle, the organic cofactor catalyzes a series of reactions, like the enzyme. The conformational change of the PLP cofactor in catalysis serves as a starting point for identifying the previously unknown catalytic roles of organic cofactors.

Divalent metal ion-based catalytic mechanism of the Nudix hydrolase Orf153 (YmfB) from Escherichia coli.

YmfB from Escherichia coli is the Nudix hydrolase involved in the metabolism of thiamine pyrophosphate, an important compound in primary metabolism and a cofactor of many enzymes. In addition, it hydrolyzes (d)NTPs to (d)NMPs and inorganic orthophosphates in a stepwise manner. The structures of YmfB alone and in complex with three sulfates and two manganese ions determined by X-ray crystallography, when compared with the structures of other Nudix hydrolases such as MutT, Ap4Aase and DR1025, provide insight into the unique hydrolysis mechanism of YmfB. Mass-spectrometric analysis confirmed that water attacks the terminal phosphates of GTP and GDP sequentially. Kinetic analysis of binding-site mutants showed that no individual residue is absolutely required for catalytic activity, suggesting that protein residues do not participate in the deprotonation of the attacking water. Thermodynamic integration calculations show that a hydroxyl ion bound to two divalent metal ions attacks the phosphate directly without the help of a nearby catalytic base.

AutoDock4(Zn): an improved AutoDock force field for small-molecule docking to zinc metalloproteins.

Zinc is present in a wide variety of proteins and is important in the metabolism of most organisms. Zinc metalloenzymes are therapeutically relevant targets in diseases such as cancer, heart disease, bacterial infection, and Alzheimer's disease. In most cases a drug molecule targeting such enzymes establishes an interaction that coordinates with the zinc ion. Thus, accurate prediction of the interaction of ligands with zinc is an important aspect of computational docking and virtual screening against zinc containing proteins. We have extended the AutoDock force field to include a specialized potential describing the interactions of zinc-coordinating ligands. This potential describes both the energetic and geometric components of the interaction. The new force field, named AutoDock4Zn, was calibrated on a data set of 292 crystal complexes containing zinc. Redocking experiments show that the force field provides significant improvement in performance in both free energy of binding estimation as well as in root-mean-square deviation from the crystal structure pose. The new force field has been implemented in AutoDock without modification to the source code.

Isomerization of Δ5-androstene-3,17-dione into Δ4-androstene-3,17-dione catalyzed by human glutathione transferase A3-3: a computational study identifies a dual role for glutathione.

Glutathione transferases (GSTs) are important enzymes in the metabolism of electrophilic xenobiotic and endobiotic toxic compounds. In addition, human GST A3-3 also catalyzes the double bond isomerization of Δ5-androstene-3,17-dione (Δ(5)-AD) and Δ(5)-pregnene-3,20-dione (Δ(5)-PD), which are the immediate precursors of testosterone and progesterone. In fact, GST A3-3 is the most efficient human enzyme known to exist in the catalysis of these reactions. In this work, we have used density functional theory (DFT) calculations to propose a refined mechanism for the isomerization of Δ(5)-AD catalyzed by GST A3-3. In this mechanism the glutathione (GSH) thiol and Tyr9 catalyze the proton transfer from the Δ(5)-AD C4 atom to the Δ(5)-AD C6 atom, with a rate limiting activation energy of 15.8 kcal · mol(-1). GSH has a dual function, because it is also responsible for stabilizing the negative charge that is formed in the O3 atom of the enolate intermediate. The catalytic role of Tyr9 depends on significant conformational rearrangements of its side chain. Neither of these contributions to catalysis has been observed before. Residues Phe10, Leu111, Ala 208, and Ala 216 complete the list of the important catalytic residues. The mechanism detailed here is based on the GST A3-3:GSH:Δ(4)-AD crystal structure and is consistent with all available experimental data.

Mining and engineering activity in catalytic amyloids.

This chapter describes how to test different amyloid preparations for catalytic properties. We describe how to express, purify, prepare and test two types of pathological amyloid (tau and α-synuclein) and two functional amyloid proteins, namely CsgA from Escherichia coli and FapC from Pseudomonas. We therefore preface the methods section with an introduction to these two examples of functional amyloid and their remarkable structural and kinetic properties and high physical stability, which renders them very attractive for a range of nanotechnological designs, both for structural, medical and catalytic purposes. The simplicity and high surface exposure of the CsgA amyloid is particularly useful for the introduction of new functional properties and we therefore provide a computational protocol to graft active sites from an enzyme of interest into the amyloid structure. We hope that the methods described will inspire other researchers to explore the remarkable opportunities provided by bacterial functional amyloid in biotechnology.

Bacillus pumilus cellulitis with bacteremia in a person who injects drugs, living with HIV-HCV co-infection: a case report.

Bacillus pumilus (B. pumilus) is a ubiquitous spore-forming bacteria that has rarely been implicated in extraintestinal infections, mostly in immunocompromised hosts. The authors report a case of B. pumilus cellulitis with bacteremia in a person who injects drugs living with human immunodeficiency virus-hepatitis C virus (HIV-HCV) co-infection. Although similar cases have been reported for some species of the genus, namely Bacillus anthracis (B. anthracis) and Bacillus cereus (B. cereus), this case reinforces the importance of considering other Bacillus spp. as potential pathogens in skin and soft tissue infections and bloodstream infections related to intravenous drug use.

What's new about the tumor microenvironment of urothelial carcinoma?

Urothelial carcinoma is a significant global health concern that accounts for a substantial part of cancer diagnoses and deaths worldwide. The tumor microenvironment is a complex ecosystem composed of stromal cells, soluble factors, and altered extracellular matrix, that mutually interact in a highly immunomodulated environment, with a prominent role in tumor development, progression, and treatment resistance. This article reviews the current state of knowledge of the different cell populations that compose the tumor microenvironment of urothelial carcinoma, its main functions, and distinct interactions with other cellular and non-cellular components, molecular alterations and aberrant signaling pathways already identified. It also focuses on the clinical implications of these findings, and its potential to translate into improved quality of life and overall survival. Determining new targets or defining prognostic signatures for urothelial carcinoma is an ongoing challenge that could be accelerated through a deeper understanding of the tumor microenvironment.

Beta-Adrenergic Blockade in Advanced Non-Small Cell Lung Cancer Patients Receiving Immunotherapy: A Multicentric Study.

Introduction The standard treatment of cancer has dramatically improved with immune checkpoint inhibitors (ICIs). Despite their proven advantage, many patients fail to exhibit a meaningful and lasting response. The beta-adrenergic signalling pathway may hold significant promise due to its role in promoting an immunosuppressive milieu within the tumour microenvironment. Inhibiting ß-adrenergic signalling could enhance ICI activity; however, blocking this pathway for this purpose has yielded conflicting results. The primary objective of this study was to evaluate the effect of beta-blocker use on overall survival and progression-free survival during ICI therapy. Methods A multicentric, retrospective, observational study was conducted in four Portuguese institutions. Patients with advanced non-small cell lung cancer treated with ICIs between January 2018 and December 2019 were included. Those using beta blockers for non-oncological reasons were compared with non-users. Results Among the 171 patients included, 36 concomitantly received beta blockers and ICIs. No significant increase was found in progression-free survival among patients who took ß-blockers (HR 0.74, 95% confidence interval (CI) 0.48-1.12, p = 0.151), and no statistically significant difference was found in overall survival. An apparent trend was observed towards better outcomes in the beta-blocker group, with a median overall survival of 9.93 months in the group not taking ß-blockers versus 14.90 months in the ß-blocker group (p = 0.291) and a median progression-free survival of 5.37 in the group not taking ß-blockers versus 10.87 months in the ß-blocker group (p = 0.151). Nine (25%) patients in the beta-blocker group and 16 (12%) in the non-beta-blocker group were progressive disease-free at the end of follow-up. This difference between the two groups is statistically significant (p = 0.047). Conclusion Our study found no statistically significant evidence that beta blockers enhance the effectiveness of immunotherapy. Using adrenergic blockade to modulate the immune system shows promise, warranting the need to develop prospective clinical studies.

Mechanism of glutathione transferase P1-1-catalyzed activation of the prodrug canfosfamide (TLK286, TELCYTA).

Canfosfamide (TLK286, TELCYTA) is a prodrug that upon activation by glutathione transferase P1-1 (GST P1-1) yields an anticancer alkylating agent and a glutathione derivative. The rationale underlying the use of TLK286 in chemotherapy is that tumor cells overexpressing GST P1-1 will be locally exposed to the released alkylating agent with limited collateral toxicity to the surrounding normal tissues. TLK286 has demonstrated clinical effects in phase II and III clinical trials for the treatment of malignancies, such as ovarian cancer, nonsmall cell lung cancer, and breast cancer, as a single agent and in combination with other chemotherapeutic agents. In spite of these promising results, the detailed mechanism of GST P1-1 activation of the prodrug has not been elucidated. Here, we propose a mechanism for the TLK286 activation by GST P1-1 on the basis of density functional theory (DFT) and on potential of mean force (PMF) calculations. A catalytic water molecule is instrumental to the activation by forming a network of intermolecular interactions between the active-site Tyr7 hydroxyl and the sulfone and COO(-) groups of TLK286. The results obtained are consistent with the available experimental kinetic data and provide an atomistic understanding of the TLK286 activation mechanism.

The first mammalian aldehyde oxidase crystal structure: insights into substrate specificity.

BACKGROUND: Aldehyde oxidases have pharmacological relevance, and AOX3 is the major drug-metabolizing enzyme in rodents. RESULTS: The crystal structure of mouse AOX3 with kinetics and molecular docking studies provides insights into its enzymatic characteristics. CONCLUSION: Differences in substrate and inhibitor specificities can be rationalized by comparing the AOX3 and xanthine oxidase structures. SIGNIFICANCE: The first aldehyde oxidase structure represents a major advance for drug design and mechanistic studies. Aldehyde oxidases (AOXs) are homodimeric proteins belonging to the xanthine oxidase family of molybdenum-containing enzymes. Each 150-kDa monomer contains a FAD redox cofactor, two spectroscopically distinct [2Fe-2S] clusters, and a molybdenum cofactor located within the protein active site. AOXs are characterized by broad range substrate specificity, oxidizing different aldehydes and aromatic N-heterocycles. Despite increasing recognition of its role in the metabolism of drugs and xenobiotics, the physiological function of the protein is still largely unknown. We have crystallized and solved the crystal structure of mouse liver aldehyde oxidase 3 to 2.9 Å. This is the first mammalian AOX whose structure has been solved. The structure provides important insights into the protein active center and further evidence on the catalytic differences characterizing AOX and xanthine oxidoreductase. The mouse liver aldehyde oxidase 3 three-dimensional structure combined with kinetic, mutagenesis data, molecular docking, and molecular dynamics studies make a decisive contribution to understand the molecular basis of its rather broad substrate specificity.