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.

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.

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.

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.

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+/+.

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.

Role of the variable active site residues in the function of thioredoxin family oxidoreductases.

The enzymes of the thioredoxin family fulfill a wide range of physiological functions. Although they possess a similar CXYC active site motif, with identical environment and stereochemical properties, the redox potential and pK(a) of the cysteine pair varies widely across the family. As a consequence, each family member promotes oxidation or reduction reactions, or even isomerization reactions. The analysis of the three-dimensional structures gives no clues to identify the molecular source for the different active site properties. Therefore, we carried out a set of quantum mechanical calculations in active site models to gain more understanding on the elusive molecular-level origin of the differentiation of the properties across the family. The obtained results, together with earlier quantum mechanical calculations performed in our laboratories, gave rise to a consistent line of evidence, which points to the fact that both active site cysteines play an important role in the differentiation. In contrary to what was assumed, differentiation is not achieved through a different stabilization of the solvent exposed cysteine but, instead, through a fine tuning of the nucleophilicity of both active site cysteines. Reductant enzymes have both cysteine thiolates poorly stabilized, oxidant proteins have both cysteine thiolates highly stabilized, and isomerases have one thiolate (solvent exposed) poorly stabilized and the other (buried) thiolate highly stabilized. The feasibility of shifting the chemical equilibrium toward oxidation, reduction, or isomerization only through subtle electrostatic effects is quite unusual, and it relies on the inherent thermoneutrality of the catalytic steps carried out by a set of chemically equivalent entities all of which are cysteine thiolates. Such pattern of stabilization/destabilization, detected in our calculations is fully consistent with the observed physiological roles of this family of enzymes.

Mechanism of thioredoxin-catalyzed disulfide reduction. Activation of the buried thiol and role of the variable active-site residues.

Thioredoxins (Trx) are enzymes with a characteristic CXYC active-site motif that catalyze the reduction of disulfide bonds in other proteins. We have theoretically explored this reaction mechanism, both in the gas phase and in water, using density functional theory. The mechanism of disulfide reduction involves two consecutive thiol-disulfide exchange reactions, that is, nucleophilic substitutions at sulfur (S(N)2@S): first, by one Trx cysteine-thiolate group (Cys-32) at a sulfur atom of the disulfide substrate and, second, by the other Trx cysteine-thiolate group (the buried thiol of Cys-35) at the sulfur atom of the first Trx cysteine. We have investigated the intrinsic nature of such S(N)2@S substitution using the simple CH3S(-) + CH3SSCH3 model and how it is affected by solvation in aqueous solution. Next, we have examined how the behavior of the elementary S(N)2@S steps changes in the more realistic enzyme-substrate model CGPC + CH3SSCH3, which contains the active-site of Trx. In all model reactions, solvation turns the hypervalent trisulfide anion (i.e., the S(N)2@S transition species) from a stable complex into a transition state. Importantly, our analyses suggest that the deprotonation of the buried thiol (which is required before the latter can enter into the second S(N)2@S step) is done by the leaving group evolving from the first S(N)2@S step. Finally, molecular dynamics (MD) simulations, in the gas phase and in water, of CGPC, CGGC, and the corresponding wild-type Trx and P34G Trx show that the activity of the thioredoxin active-site motif (CXYC) is determined not only by the structural rigidity associated with the particular variable residues (XY) but also by the number of amide N-H groups. The latter are involved in the stabilization of the Cys-32 thiolate and thus affect the acidity and nucleophilicity of this residue.

Insights on resistance to reverse transcriptase: the different patterns of interaction of the nucleoside reverse transcriptase inhibitors in the deoxyribonucleotide triphosphate binding site relative to the normal substrate.

It is presently known that the long-term failure in the treatment of AIDS with the currently available nucleotide reverse transcriptase inhibitors (NRTIs) is related to the development of resistance by reverse transcriptase (RT) at the binding or incorporation level or both, or subsequent to the nucleotide incorporation (excision). To achieve greater insight on the differential interactions of two NRTIs that are mainly discriminated by different mechanisms, 2',3'-didehydro-2',3'-dideoxythymidine-5'-triphosphate (d4TTP, that is, phosphorylated stavudine) and 2',3'-dideoxycytidine-5'-triphosphate (ddCTP, that is, phosphorylated zalcitabine), with the primer/template (p/t) and with the N binding site of reverse transcriptase (RT) in relation to the normal substrate (dNTP), we have conducted a series of molecular dynamics (MD) simulations. We propose that the different resistance profiles arise from the different conformations adopted by the inhibitors at the N site. d4TTP adopts an ideal conformation for catalysis because it forms an ion-dipole intramolecular interaction with the beta-phosphate oxygen of the triphosphate, as does the normal substrate. In ddCTP, the lack of this essential interaction results in a different, noncatalytic conformation.

Determination of the DeltapKa between the active site cysteines of thioredoxin and DsbA.

Thioredoxin superfamily members share a considerable degree of structural similarity, with a conserved CX(i)X(j)C motif at the active site, where C stand for two cysteines that alternate between a reduced thiol and oxidized disulfide states, and X(i)and X(j) are two amino acids different in each family member. Despite these similarities, they display very different redox potentials and pKas for the active site dithiol, and fulfill different physiological roles. Thioredoxin, for example, promotes the reduction of disulfide bonds, while DsbA promotes their oxidation in prokaryotic cells. The factors that promote these differences are still not fully understood. However, it is generally accepted that the different stabilities of the redox active disulfide bond depends on the degree of stabilization, in the reduced state, of the thiolate of one of the active site cysteines (nucleophilic cysteine). In this work we have used QM/MM methods to compare and characterize the active site dithiols of both enzymes, and to shed some light on the structural features responsible for the large differences in pKa and redox potential between two homologous enzymes, thioredoxin and DsbA. We have also analyzed the main factors pointed out in the literature as responsible for their different properties. We obtained the value of 4.5 for pKa difference (DeltapKa) between the nucleophilic cysteines of both enzymes, which is in excellent agreement with most of the experimental values. Additionally, we found that the principal differentiating factor responsible for this observed DeltapKa are the alpha2-alpha helices, which greatly contribute to the mentioned value, by stabilizing the DsbA thiolate in a much greater extend than the thioredoxin thiolate. A double mutation of the conserved residues Asp26 and Lys57, in thioredoxin, and Glu24 Lys58, in DsbA, by alanines did not change the DeltapKa value; this supports the hypothesis that these residues are not involved in the differentiation of the properties of the active centre dithiol. However, we found out that these residues are important for the stabilization of the nucleophilic thiolate. The X(i) and X(j) residues also do not seem to promote the stabilization of the thiolates. In fact, the corresponding double alanine mutants are more stable than the wild-type enzymes. However, these residues are involved in the differentiation between thioredoxin and DsbA, stabilizing the DsbA thiolate by a larger extent than the thioredoxin thiolate.

Theoretical study of the unusual protonation properties of the active site cysteines in thioredoxin.

The thiol/disulfide oxidoreductases of the thioredoxin family have, in the active site, two cysteines that can be in a reduced or an oxidized form. One of the cysteines in the reduced state is deprotonated, and it is called nucleophilic cysteine. The pK(a) of this cysteine is different from that of a normal cysteine and varies widely among the different enzymes of this family. However, the factors responsible for the different degrees of stabilization of nucleophilic cysteine thiolate are not fully understood. Here, we have studied the well-known hypothesis of proton sharing between the active site thiols by performing a linear transit scan for the transfer of the proton between the active site cysteines. We used a two-layered (DFT/MM) ONIOM formalism, with the active site region treated at the B3LYP/6-31+G(d) level and the remains of the protein treated with the Amber Parm94 force field. The solvation free energy was accounted for with a continuum solvent model, by solving the Poisson-Boltzmann equation using the program Delphi. We have obtained excellent agreement with the experimental data available in the literature. Besides refuting the proton sharing hypothesis, our results include a value of 14.0 for the pK(a) of the buried cysteine, a quantity that has not been possible to obtain experimentally but which has been proven to be higher than 11. Additionally, this study also provides detailed information on the very interesting and so far unknown fact that the contribution of the enzymatic structure (8.3 kcal/mol) prevails in relation to that of the solvent (0.60 kcal/mol) concerning the differential stabilization of the active site thiolates.