Application of rational enzyme engineering in a new route to etonogestrel and levonorgestrel: carbonyl reductase bioreduction of ethyl secodione.

Women in developing countries still face enormous challenges when accessing reproductive health care. Access to voluntary family planning empowers women allowing them to complete their education and join the paid workforce. This effectively helps to end poverty, hunger and promotes good health for all. According to the United Nations (UN) organization, in 2022, an estimated 257 million women still lacked access to safe and effective family planning methods globally. One of the main barriers is the associated cost of modern contraceptive methods. Funded by the Bill & Melinda Gates Foundation, Almac Group worked on the development of a novel biocatalytic route to etonogestrel and levonorgestrel, two modern contraceptive APIs, with the goal of substantially decreasing the cost of production and so enabling their use in developing nations. This present work combines the selection and engineering of a carbonyl reductase (CRED) enzyme from Almac's selectAZyme™ panel, with process development, to enable efficient and economically viable bioreduction of ethyl secodione to (13R,17S)-secol, the key chirality introducing intermediate en route to etonogestrel and levonorgestrel API. CRED library screening returned a good hit with an Almac CRED from Bacillus weidmannii, which allowed for highly stereoselective bioreduction at low enzyme loading of less than 1% w/w under screening assay conditions. However, the only co-solvent tolerated was DMSO up to ∼30% v/v, and it was impossible to achieve reaction completion with any enzyme loading at substrate titres of 20 g L-1 and above, due to the insolubility of the secodione. This triggered a rapid enzyme engineering program fully based on computational mutant selection. A small panel of 93 CRED mutants was rationally designed to increase the catalytic activity as well as thermal and solvent stability. The best mutant, Mutant-75, enabled a reaction at 45 °C to go to completion at 90 g L-1 substrate titre in a buffer/DMSO/heptane reaction medium fed over 6 h with substrate DMSO stock solution, with a low enzyme loading of 3.5% w/w wrt substrate. In screening assay conditions, Mutant-75 also showed a 2.2-fold activity increase. Our paper shows which computations and rational decisions enabled this outcome.

Designing Cell Delivery Peptides and SARS-CoV-2-Targeting Small Interfering RNAs: A Comprehensive Bioinformatics Study with Generative Adversarial Network-Based Peptide Design and In Vitro Assays.

Nucleic acid technologies with designed intracellular delivery systems are some of the most promising therapies of the future. Small interfering (si)RNAs inhibit gene expression and protein synthesis and may complement current vaccines with faster design and production. Although successful delivery remains an issue, delivery peptides may help to fill this gap. Here, we address this issue by applying bioinformatic approaches to design new putative cell delivery peptides and siRNAs for COVID-19 variants and other related viral diseases. Of the 29,880 RNA sequences analyzed, 62 were identified in silico as able to target the virus mRNA sequence, and from the 9,984 peptide sequences analyzed, 10 were selected as delivery peptides. From the latter, we further performed in vitro studies of the two best-ranked peptides and compared them with the broadly used TAT delivery peptide. One of them, seq5, displayed better internalization results with about double intensity signal compared to TAT after a 1 h incubation time in GFP-HeLa cells. This peptide has, thus, the features of a delivery peptide and could be used for cargo intracellular delivery.

An injectable, naproxen-conjugated, supramolecular hydrogel with ultra-low critical gelation concentration-prepared from a known folate receptor ligand.

Short peptides capped on the N-terminus with aromatic groups are often able to form supramolecular hydrogels-self-assembled networks of fibrils able to trap water molecules. Typically, these hydrogelators can form stiff gels at concentrations of 0.1 to 1.0 wt%-i.e. they consist of mainly water. The properties of these soft materials mimic those of the extracellular matrix (ECM) of biological tissue and therefore they have found many biomedical uses in tissue engineering, wound healing, drug delivery, biosensing and bioprinting applications. In drug delivery strategies related to cancer therapy, injectable hydrogels can serve as a depot formulation, where a sustained release of the chemotherapeutic from near the tumour site allows reduced doses and, therefore, decreased side effects. To further target cancer cells, folic acid-conjugated hydrogels and nanostructures are often sought, to exploit the overexpression of folate receptors on cancer cells-an approach which can allow the selective cellular uptake of an encapsulated drug. In this present study, two known dipeptide folate receptor ligands (1 and 2) recently identified from a screen of a DNA-encoded compound library, were synthesised and investigated for their hydrogelation ability and cytotoxicity. Compound 1, containing a naproxen capping group, rapidly forms hydrogels at concentrations as low as 0.03 wt%-one of the lowest critical gelation concentrations (CGCs) known for a supramolecular hydrogelator. In contrast, compound 2, which contains a 3-indolepropionic acid capping group, was unable to form hydrogels under a range of conditions and concentrations, instead forming nanospheres with diameters of 0.5 µm. Hydrogels of 1 were characterised by STEM microscopy, rheology, fluorescence spectroscopy and circular dichroism. Both compounds 1 and 2 had no impact on the proliferation of kerotinocytes (HaCaT cells) at concentrations up to 100 µM. Compound 1, containing the NSAID, was tested for anti-inflammatory activity in a human cyclooxygenase-1/2 model. The rate of the release of model drug compounds from within hydrogels of 1 was also investigated.

Correction to: In vitro mechanistic studies on α-amanitin and its putative antidotes.

In the original publication of the article.

In vitro mechanistic studies on α-amanitin and its putative antidotes.

α-Amanitin plays a key role in Amanita phalloides intoxications. The liver is a major target of α-amanitin toxicity, and while RNA polymerase II (RNA Pol II) transcription inhibition is a well-acknowledged mechanism of α-amanitin toxicity, other possible toxicological pathways remain to be elucidated. This study aimed to assess the mechanisms of α-amanitin hepatotoxicity in HepG2 cells. The putative protective effects of postulated antidotes were also tested in this cell model and in permeabilized HeLa cells. α-Amanitin (0.1-20 µM) displayed time- and concentration-dependent cytotoxicity, when evaluated through the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) reduction and neutral red uptake assays. Additionally, α-amanitin decreased nascent RNA synthesis in a concentration- and time-dependent manner. While α-amanitin did not induce changes in mitochondrial membrane potential, it caused a significant increase in intracellular ATP levels, which was not prevented by incubation with oligomycin, an ATP synthetase inhibitor. Concerning the cell redox status, α-amanitin did not increase reactive species production, but caused a significant increase in total and reduced glutathione, which was abolished by pre-incubation with the inhibitor of gamma-glutamylcysteine synthase, buthionine sulfoximine. None of the tested antidotes [N-acetyl cysteine, silibinin, benzylpenicillin, and polymyxin B (PolB)] conferred any protection against α-amanitin-induced cytotoxicity in HepG2 cells or reversed the inhibition of nascent RNA caused by the toxin in permeabilized HeLa cells. Still, PolB interfered with RNA Pol II activity at high concentrations, though not impacting on α-amanitin observed cytotoxicity. New hepatotoxic mechanisms of α-amanitin were described herein, but the lack of protection observed in clinically used antidotes may reflect the lack of knowledge on their true protection mechanisms and may explain their relatively low clinical efficacy.

SLMP53-1 Inhibits Tumor Cell Growth through Regulation of Glucose Metabolism and Angiogenesis in a P53-Dependent Manner.

The Warburg effect is an emerging hallmark of cancer, which has the tumor suppressor p53 as its major regulator. Herein, we unveiled that p53 activation by (S)-tryptophanol-derived oxazoloisoindolinone (SLMP53-1) mediated the reprograming of glucose metabolism in cancer cells and xenograft human tumor tissue, interfering with angiogenesis and migration. Particularly, we showed that SLMP53-1 regulated glycolysis by downregulating glucose transporter 1 (GLUT1), hexokinase-2 (HK2), and phosphofructokinase-2 isoform 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase-3 (PFKFB3) (key glycolytic enzymes), while upregulating the mitochondrial markers synthesis of cytochrome c oxidase 2 (SCO2), cytochrome c oxidase subunit 4 (COX4), and OXPHOS mitochondrial complexes. SLMP53-1 also downregulated the monocarboxylate transporter 4 (MCT4), causing the subsequent reduction of lactate export by cancer cells. Besides the acidification of the extracellular environment, SLMP53-1 further increased E-cadherin and reduced metalloproteinase-9 (MMP-9) expression levels in both cancer cells and xenograft human tumor tissue, which suggested the interference of SLMP53-1 in extracellular matrix remodeling and epithelial-to-mesenchymal transition. Consistently, SLMP53-1 depleted angiogenesis, decreasing endothelial cell tube formation and vascular endothelial growth factor (VEGF) expression levels. SLMP53-1 also exhibited synergistic growth inhibitory activity in combination with the metabolic modulator dichloroacetic acid. These data reinforce the promising application of the p53-activating agent SLMP53-1 in cancer therapy, by targeting p53-mediated pathways of growth and dissemination.

Why the Flavin Adenine Dinucleotide (FAD) Cofactor Needs To Be Covalently Linked to Complex†II of the Electron-Transport Chain for the Conversion of FADH2 into FAD.

A covalently bound flavin cofactor is predominant in the succinate-ubiquinone oxidoreductase (SQR; Complex II), an essential component of aerobic electron transport, and in the menaquinol-fumarate oxidoreductase (QFR), the anaerobic counterpart, although it is only present in approximately 10 % of the known flavoenzymes. This work investigates the role of this 8α-N3-histidyl linkage between the flavin adenine dinucleotide (FAD) cofactor and the respiratory Complex II. After parameterization with DFT calculations, classical molecular-dynamics simulations and quantum-mechanics calculations for Complex II:FAD and Complex II:FADH2 , with and without the covalent bond, were performed. It was observed that the covalent bond is essential for the active-center arrangement of the FADH2 /FAD cofactor. Removal of this bond causes a displacement of the isoalloxazine group, which influences interactions with the protein, flavin solvation, and possible proton-transfer pathways. Specifically, for the noncovalently bound FADH2 cofactor, the N1 atom moves away from the His-A365 and His-A254 residues and the N5 atom moves away from the glutamine-62A residue. Both of the histidine and glutamine residues interact with a chain of water molecules that cross the enzyme, which is most likely involved in proton transfer. Breaking this chain of water molecules could thereby compromise proton transfer across the two active sites of Complex II.

Spatial requirement for PAMO for transformation of non-native linear substrates.

Phenylacetone monooxygenase is the most stable and thermo-tolerant member of the Baeyer-Villiger monooxygenases family, and therefore it is an ideal candidate for the synthesis of industrially relevant ester or lactone compounds. However, its limited substrate scope has largely limited its industrial applications. Linear substrates are interesting from an industrial point of view, it is thus necessary to identify the essential spatial requirement for achieving high conversions for non-native linear substrates. Here using molecular dynamics simulations, we compared the conversion of a non-native linear substrate 2-octanone and the native substrate phenylacetone, catalyzed by the WT enzyme and a quadruple variant P253F/G254A/R258M/L443F that exhibits significantly improved activity towards 2-octanone. We uncovered that a remarkable movement of L289 is crucial for a reshaping of the active site of the quadruple variant so as to prevent the aliphatic substrate from moving away from the C4a-peroxyflavin, thus enabling it to keep a catalytically relevant pose during the oxygenation process. By performing steady-state kinetic analysis of two single-mutation variants at position 258, we further validated that the L289 reposition is attributed to the combined effect of quadruple mutations. In order to further explore the substrate scope of PAMO we also studied the binding of cyclopentanone and 2-phenylcyclohexanone, which are the typical substrates of CPMO in group I and CHMO in group III, respectively. Our study provides fundamental atomic-level insights in rational engineering of PAMO for wide applications in industrial biocatalysis, in particular, in the biotransformation of long-chain aliphatic oils into potential biodiesels.

Catalytic mechanism of phenylacetone monooxygenases for non-native linear substrates.

Phenylacetone monooxygenase (PAMO) is the most stable and thermo-tolerant member of the Baeyer-Villiger monooxygenase family, and therefore it is an ideal candidate for the synthesis of industrially relevant compounds. However, its limited substrate scope has largely limited its industrial applications. In the present work, we provide, for the first time, the catalytic mechanism of PAMO for the native substrate phenylacetone as well as for a linear non-native substrate 2-octanone, using molecular dynamics simulations, quantum mechanics and quantum mechanics/molecular mechanics calculations. We provide a theoretical basis for the preference of the enzyme for the native aromatic substrate over non-native linear substrates. Our study provides fundamental atomic-level insights that can be employed in the rational engineering of PAMO for wide applications in industrial biocatalysis, in particular, in the biotransformation of long-chain aliphatic oils into potential biodiesels.

Cooperative Electrostatic Interactions Drive Functional Evolution in the Alkaline Phosphatase Superfamily.

It is becoming widely accepted that catalytic promiscuity, i.e., the ability of a single enzyme to catalyze the turnover of multiple, chemically distinct substrates, plays a key role in the evolution of new enzyme functions. In this context, the members of the alkaline phosphatase superfamily have been extensively studied as model systems in order to understand the phenomenon of enzyme multifunctionality. In the present work, we model the selectivity of two multiply promiscuous members of this superfamily, namely the phosphonate monoester hydrolases from Burkholderia caryophylli and Rhizobium leguminosarum. We have performed extensive simulations of the enzymatic reaction of both wild-type enzymes and several experimentally characterized mutants. Our computational models are in agreement with key experimental observables, such as the observed activities of the wild-type enzymes, qualitative interpretations of experimental pH-rate profiles, and activity trends among several active site mutants. In all cases the substrates of interest bind to the enzyme in similar conformations, with largely unperturbed transition states from their corresponding analogues in aqueous solution. Examination of transition-state geometries and the contribution of individual residues to the calculated activation barriers suggest that the broad promiscuity of these enzymes arises from cooperative electrostatic interactions in the active site, allowing each enzyme to adapt to the electrostatic needs of different substrates. By comparing the structural and electrostatic features of several alkaline phosphatases, we suggest that this phenomenon is a generalized feature driving selectivity and promiscuity within this superfamily and can be in turn used for artificial enzyme design.

Modeling the mechanisms of biological GTP hydrolysis.

Enzymes that hydrolyze GTP are currently in the spotlight, due to their molecular switch mechanism that controls many cellular processes. One of the best-known classes of these enzymes are small GTPases such as members of the Ras superfamily, which catalyze the hydrolysis of the γ-phosphate bond in GTP. In addition, the availability of an increasing number of crystal structures of translational GTPases such as EF-Tu and EF-G have made it possible to probe the molecular details of GTP hydrolysis on the ribosome. However, despite a wealth of biochemical, structural and computational data, the way in which GTP hydrolysis is activated and regulated is still a controversial topic and well-designed simulations can play an important role in resolving and rationalizing the experimental data. In this review, we discuss the contributions of computational biology to our understanding of GTP hydrolysis on the ribosome and in small GTPases.

Understanding thio-effects in simple phosphoryl systems: role of solvent effects and nucleophile charge.

Recent experimental work (J. Org. Chem., 2012, 77, 5829) demonstrated pronounced differences in measured thio-effects for the hydrolysis of (thio)phosphodichloridates by water and hydroxide nucleophiles. In the present work, we have performed detailed quantum chemical calculations of these reactions, with the aim of rationalizing the molecular bases for this discrimination. The calculations highlight the interplay between nucleophile charge and transition state solvation in SN2(P) mechanisms as the basis of these differences, rather than a change in mechanism.

A breakthrough on Amanita phalloides poisoning: an effective antidotal effect by polymyxin B.

Amanita phalloides is responsible for more than 90 % of mushroom-related fatalities, and no effective antidote is available. α-Amanitin, the main toxin of A. phalloides, inhibits RNA polymerase II (RNAP II), causing hepatic and kidney failure. In silico studies included docking and molecular dynamics simulation coupled to molecular mechanics with generalized Born and surface area method energy decomposition on RNAP II. They were performed with a clinical drug that shares chemical similarities to α-amanitin, polymyxin B. The results show that polymyxin B potentially binds to RNAP II in the same interface of α-amanitin, preventing the toxin from binding to RNAP II. In vivo, the inhibition of the mRNA transcripts elicited by α-amanitin was efficiently reverted by polymyxin B in the kidneys. Moreover, polymyxin B significantly decreased the hepatic and renal α-amanitin-induced injury as seen by the histology and hepatic aminotransferases plasma data. In the survival assay, all animals exposed to α-amanitin died within 5 days, whereas 50 % survived up to 30 days when polymyxin B was administered 4, 8, and 12 h post-α-amanitin. Moreover, a single dose of polymyxin B administered concomitantly with α-amanitin was able to guarantee 100 % survival. Polymyxin B protects RNAP II from inactivation leading to an effective prevention of organ damage and increasing survival in α-amanitin-treated animals. The present use of clinically relevant concentrations of an already human-use-approved drug prompts the use of polymyxin B as an antidote for A. phalloides poisoning in humans.

Electronic structure investigation and parametrization of biologically relevant iron-sulfur clusters.

The application of classical molecular dynamics simulations to the study of metalloenzymes has been hampered by the lack of suitable molecular mechanics force field parameters to treat the metal centers within standard biomolecular simulation packages. These parameters cannot be generalized, nor be easily automated, and hence should be obtained for each system separately. Here we present density functional theory calculations on [Fe2S2(SCH3)4]2+/+, [Fe3S4(SCH3)3]+/0 and [Fe4S4(SCH3)4]2+/+ and the derivation of parameters that are compatible with the AMBER force field. Molecular dynamics simulations performed using these parameters on respiratory Complex II of the electron transport chain showed that the reduced clusters are more stabilized by the protein environment, which leads to smaller changes in bond lengths and angles upon reduction. This effect is larger in the smaller iron-sulfur cluster, [Fe2S2(SCH3)4]2+/+.

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.

Parameters for molecular dynamics simulations of iron-sulfur proteins.

Iron-sulfur proteins involved in electron transfer reactions have finely tuned redox potentials, which allow them to be highly efficient and specific. Factors such as metal center solvent exposure, interaction with charged residues, or hydrogen bonds between the ligand residues and amide backbone groups have all been pointed out to cause such specific redox potentials. Here, we derived parameters compatible with the AMBER force field for the metal centers of iron-sulfur proteins and applied them in the molecular dynamics simulations of three iron-sulfur proteins. We used density-functional theory (DFT) calculations and Seminario's method for the parameterization. Parameter validation was obtained by matching structures and normal frequencies at the quantum mechanics and molecular mechanics levels of theory. Having guaranteed a correct representation of the protein coordination spheres, the amide H-bonds and the water exposure to the ligands were analyzed. Our results for the pattern of interactions with the metal centers are consistent to those obtained by nuclear magnetic resonance spectroscopy (NMR) experiments and DFT calculations, allowing the application of molecular dynamics to the study of those proteins.

Insights into the Degradation of Polymer-Drug Conjugates by an Overexpressed Enzyme in Cancer Cells.

Intensive efforts have been made to provide better treatments to cancer patients. Currently, nanoparticle-based drug delivery systems have gained propulsion, as they can overcome the drawbacks of free drugs. However, drug stability inside the nanocapsule must be ensured to prevent burst release. To overcome this, drugs conjugated to amphiphilic copolymers, assembled into nanoparticles, can provide a sustained release if endogenously degraded. Thus, we have designed and assessed the drug release viability of polymer-drug conjugates by the human Carboxylesterase 2, for a targeted drug activation. We performed molecular dynamics simulations applying a quantum mechanics/molecular mechanics potential to study the degradation profiles of 30 designed conjugates, where six were predicted to be hydrolyzed by this enzyme. We further analyzed the enzyme-substrate environment to delve into what structural features may lead to successful hydrolysis. These findings contribute to the development of new medicines ensuring effective cancer treatments with fewer side effects.

Discovery and Rational Mutagenesis of Methionine Sulfoxide Reductase Biocatalysts To Expand the Substrate Scope of the Kinetic Resolution of Chiral Sulfoxides.

Methionine sulfoxide reductase A (MsrA) enzymes have recently found applications as nonoxidative biocatalysts in the enantioselective kinetic resolution of racemic sulfoxides. This work describes the identification of selective and robust MsrA biocatalysts able to catalyze the enantioselective reduction of a variety of aromatic and aliphatic chiral sulfoxides at 8-64 mM concentration with high yields and excellent ees (up to 99%). Moreover, with the aim to expand the substrate scope of MsrA biocatalysts, a library of mutant enzymes has been designed via rational mutagenesis utilizing in silico docking, molecular dynamics, and structural nuclear magnetic resonance (NMR) studies. The mutant enzyme MsrA33 was found to catalyze the kinetic resolution of bulky sulfoxide substrates bearing non-methyl substituents on the sulfur atom with ees up to 99%, overcoming a significant limitation of the currently available MsrA biocatalysts.

Introduction of a Glycine Linker Connecting the Heavy and Light Chains in Synthetic Cardosin B-Derived Rennet Changes the Specificity of Subpocket S3'.

The development of plant-based synthetic rennets is of high commercial interest, due to the current great consumer demand for animal product alternatives. A previously developed recombinant form of the aspartic protease cardosin B with a three-glycine linker showed great potential due to its good performance in milk coagulation. This enzyme was found to be more specific and less proteolytically active than the native form for milk clotting, but the underlying structural causes for these activity changes were not completely clear. Here, we have performed molecular dynamics simulations with the recombinant enzyme with and without the linker. Our results showed that the introduction of the linker changes the subpocket S3', which is located more than 4 nm away. These results showcase how small modifications in proteins can have significant effects in distant regions in the protein structure that affect their biotechnological applications.

Transcription factor allosteric regulation through substrate coordination to zinc.

The development of new synthetic biology circuits for biotechnology and medicine requires deeper mechanistic insight into allosteric transcription factors (aTFs). Here we studied the aTF UxuR, a homodimer of two domains connected by a highly flexible linker region. To explore how ligand binding to UxuR affects protein dynamics we performed molecular dynamics simulations in the free protein, the aTF bound to the inducer D-fructuronate or the structural isomer D-glucuronate. We then validated our results by constructing a sensor plasmid for D-fructuronate in Escherichia coli and performed site-directed mutagenesis. Our results show that zinc coordination is necessary for UxuR function since mutation to alanines prevents expression de-repression by D-fructuronate. Analyzing the different complexes, we found that the disordered linker regions allow the N-terminal domains to display fast and large movements. When the inducer is bound, UxuR can sample an open conformation with a more pronounced negative charge at the surface of the N-terminal DNA binding domains. In opposition, in the free and D-glucuronate bond forms the protein samples closed conformations, with a more positive character at the surface of the DNA binding regions. These molecular insights provide a new basis to harness these systems for biological systems engineering.