Autoantibodies against chemokines post-SARS-CoV-2 infection correlate with disease course.

Infection with severe acute respiratory syndrome coronavirus 2 associates with diverse symptoms, which can persist for months. While antiviral antibodies are protective, those targeting interferons and other immune factors are associated with adverse coronavirus disease 2019 (COVID-19) outcomes. Here we discovered that antibodies against specific chemokines were omnipresent post-COVID-19, were associated with favorable disease outcome and negatively correlated with the development of long COVID at 1 yr post-infection. Chemokine antibodies were also present in HIV-1 infection and autoimmune disorders, but they targeted different chemokines compared with COVID-19. Monoclonal antibodies derived from COVID-19 convalescents that bound to the chemokine N-loop impaired cell migration. Given the role of chemokines in orchestrating immune cell trafficking, naturally arising chemokine antibodies may modulate the inflammatory response and thus bear therapeutic potential.

A first-in-class Wiskott-Aldrich syndrome protein activator with anti-tumor activity in hematologic cancers.

Hematological cancers are among the most common cancers in adults and children. Despite significant improvements in therapies, many patients still succumb to the disease. Therefore, novel therapies are needed. The Wiskott-Aldrich syndrome protein (WASp) family regulates actin assembly in conjunction with the Arp2/3 complex, a ubiquitous nucleation factor. WASp is expressed exclusively in hematopoietic cells and exists in two allosteric conformations: autoinhibited or activated. Here, we describe the development of EG-011, a first-in-class small molecule activator of the WASp auto-inhibited form. EG-011 possesses in vitro and in vivo anti-tumor activity as a single agent in lymphoma, leukemia, and multiple myeloma, including models of secondary resistance to PI3K, BTK, and proteasome inhibitors. The in vitro activity was confirmed in a lymphoma xenograft. Actin polymerization and WASp binding was demonstrated using multiple techniques. Transcriptome analysis highlighted homology with drugs-inducing actin polymerization.

Computational Design of Novel Cyclic Peptides Endowed with Autophagy-Inhibiting Activity on Cancer Cell Lines.

(1) Autophagy plays a significant role in development and cell proliferation. This process is mainly accomplished by the LC3 protein, which, after maturation, builds the nascent autophagosomes. The inhibition of LC3 maturation results in the interference of autophagy activation. (2) In this study, starting from the structure of a known LC3B binder (LIR2-RavZ peptide), we identified new LC3B ligands by applying an in silico drug design strategy. The most promising peptides were synthesized, biophysically assayed, and biologically evaluated to ascertain their potential antiproliferative activity on five humans cell lines. (3) A cyclic peptide (named Pep6), endowed with high conformational stability (due to the presence of a disulfide bridge), displayed a Kd value on LC3B in the nanomolar range. Assays accomplished on PC3, MCF-7, and A549 cancer cell lines proved that Pep6 exhibited cytotoxic effects comparable to those of the peptide LIR2-RavZ, a reference LC3B ligand. Furthermore, it was ineffective on both normal prostatic epithelium PNT2 and autophagy-defective prostate cancer DU145 cells. (4) Pep6 can be considered a new autophagy inhibitor that can be employed as a pharmacological tool or even as a template for the rational design of new small molecules endowed with autophagy inhibitory activity.

ETS1 phosphorylation at threonine 38 is associated with the cell of origin of diffuse large B cell lymphoma and sustains the growth of tumour cells.

The transcriptional factor ETS1 is upregulated in 25% of diffuse large B cell lymphoma (DLBCL). Here, we studied the role of ETS1 phosphorylation at threonine 38, a marker for ETS1 activation, in DLBCL cellular models and clinical specimens. p-ETS1 was detected in activated B cell-like DLBCL (ABC), not in germinal centre B-cell-like DLBCL (GCB) cell lines and, accordingly, it was more common in ABC than GCB DLBCL diagnostic biopsies. MEK inhibition decreased both baseline and IgM stimulation-induced p-ETS1 levels. Genetic inhibition of phosphorylation of ETS1 at threonine 38 affected the growth and the BCR-mediated transcriptome program in DLBCL cell lines. Our data demonstrate that ETS1 phosphorylation at threonine 38 is important for the growth of DLBCL cells and its pharmacological inhibition could benefit lymphoma patients.

Esmethadone-HCl (REL-1017): a promising rapid antidepressant.

This review article presents select recent studies that form the basis for the development of esmethadone into a potential new drug. Esmethadone is a promising member of the pharmacological class of uncompetitive N-methyl-D-aspartate receptor (NMDAR) antagonists that have shown efficacy for major depressive disorder (MDD) and other diseases and disorders, such as Alzheimer's dementia and pseudobulbar affect. The other drugs in the novel class of NMDAR antagonists with therapeutic uses that are discussed for comparative purposes in this review are esketamine, ketamine, dextromethorphan, and memantine. We present in silico, in vitro, in vivo, and clinical data for esmethadone and other uncompetitive NMDAR antagonists that may advance our understanding of the role of these receptors in neural plasticity in health and disease. The efficacy of NMDAR antagonists as rapid antidepressants may advance our understanding of the neurobiology of MDD and other neuropsychiatric diseases and disorders.

Focused Design of Novel Cyclic Peptides Endowed with GABARAP-Inhibiting Activity.

(1) Background: Disfunctions in autophagy machinery have been identified in various conditions, including neurodegenerative diseases, cancer, and inflammation. Among mammalian autophagy proteins, the Atg8 family member GABARAP has been shown to be greatly involved in the autophagy process of prostate cancer cells, supporting the idea that GABARAP inhibitors could be valuable tools to fight the progression of tumors. (2) Methods: In this paper, starting from the X-ray crystal structure of GABARAP in a complex with an AnkirinB-LIR domain, we identify two new peptides by applying in silico drug design techniques. The two ligands are synthesized, biophysically assayed, and biologically evaluated to ascertain their potential anticancer profile. (3) Results: Two cyclic peptides (WC8 and WC10) displayed promising biological activity, high conformational stability (due to the presence of disulfide bridges), and Kd values in the low micromolar range. The anticancer assays, performed on PC-3 cells, proved that both peptides exhibit antiproliferative effects comparable to those of peptide K1, a known GABARAP inhibitor. (4) Conclusions: WC8 and WC10 can be considered new GABARAP inhibitors to be employed as pharmacological tools or even templates for the rational design of new small molecules.

How Computational Chemistry and Drug Delivery Techniques Can Support the Development of New Anticancer Drugs.

The early and late development of new anticancer drugs, small molecules or peptides can be slowed down by some issues such as poor selectivity for the target or poor ADME properties. Computer-aided drug design (CADD) and target drug delivery (TDD) techniques, although apparently far from each other, are two research fields that can give a significant contribution to overcome these problems. Their combination may provide mechanistic understanding resulting in a synergy that makes possible the rational design of novel anticancer based therapies. Herein, we aim to discuss selected applications, some also from our research experience, in the fields of anticancer small organic drugs and peptides.

Structural Biology of STAT3 and Its Implications for Anticancer Therapies Development.

Transcription factors are proteins able to bind DNA and induce the transcription of specific genes. Consequently, they play a pivotal role in multiple cellular pathways and are frequently over-expressed or dysregulated in cancer. Here, we will focus on a specific "signal transducer and activator of transcription" (STAT3) factor that is involved in several pathologies, including cancer. For long time, the mechanism by which STAT3 exerts its cellular functions has been summarized by a three steps process: (1) Protein phosphorylation by specific kinases, (2) dimerization promoted by phosphorylation, (3) activation of gene expression by the phosphorylated dimer. Consequently, most of the inhibitors reported in literature aimed at blocking phosphorylation and dimerization. However, recent observations reopened the debate and the entire functional mechanism has been revisited stimulating the scientific community to pursue new inhibition strategies. In particular, the dimerization of the unphosphorylated species has been experimentally demonstrated and specific roles proposed also for these dimers. Despite difficulties in the expression and purification of the full length STAT3, structural biology investigations allowed the determination of atomistic structures of STAT3 dimers and several protein domains. Starting from this information, computational methods have been used both to improve the understanding of the STAT3 functional mechanism and to design new inhibitors to be used as anticancer drugs. In this review, we will focus on the contribution of structural biology to understand the roles of STAT3, to design new inhibitors and to suggest new strategies of pharmacological intervention.

Insight into the Mechanism of Hydrolysis of Meropenem by OXA-23 Serine-ß-lactamase Gained by Quantum Mechanics/Molecular Mechanics Calculations.

The fast and constant development of drug resistant bacteria represents a serious medical emergency. To overcome this problem, the development of drugs with new structures and modes of action is urgently needed. In this work, we investigated, at the atomistic level, the mechanisms of hydrolysis of Meropenem by OXA-23, a class D ß-lactamase, combining unbiased classical molecular dynamics and umbrella sampling simulations with classical force field-based and quantum mechanics/molecular mechanics potentials. Our calculations provide a detailed structural and dynamic picture of the molecular steps leading to the formation of the Meropenem-OXA-23 covalent adduct, the subsequent hydrolysis, and the final release of the inactive antibiotic. In this mechanistic framework, the predicted activation energy is in good agreement with experimental kinetic measurements, validating the expected reaction path.

Structure-based approach for identification of novel phenylboronic acids as serine-ß-lactamase inhibitors.

ß-Lactamases are bacterial enzymes conferring resistance to ß-lactam antibiotics in clinically-relevant pathogens, and represent relevant drug targets. Recently, the identification of new boronic acids (i.e. RPX7009) paved the way to the clinical application of these molecules as potential drugs. Here, we screened in silico a library of ~1400 boronic acids as potential AmpC ß-lactamase inhibitors. Six of the most promising candidates were evaluated in biochemical assays leading to the identification of potent inhibitors of clinically-relevant ß-lactamases like AmpC, KPC-2 and CTX-M-15. One of the selected compounds showed nanomolar K i value with the clinically-relevant KPC-2 carbapenemase, while another one exhibited broad spectrum inhibition, being also active on Enterobacter AmpC and the OXA-48 class D carbapenemase.

Molecular Determinants for Unphosphorylated STAT3 Dimerization Determined by Integrative Modeling.

Signal transducer and activator of transcription factors (STATs) are proteins that can translocate into the nucleus, bind DNA, and activate gene transcription. STAT proteins play a crucial role in cell proliferation, apoptosis, and differentiation. The prevalent view is that STAT proteins are able to form dimers and bind DNA only upon phosphorylation of specific tyrosine residues in the transactivation domain. However, this paradigm has been questioned recently by the observation of dimers of unphosphorylated STATs (USTATs) by X-ray, Förster resonance energy transfer, and site-directed mutagenesis. A more complex picture of the dimerization process and of the role of the dimers is, thus, emerging. Here we present an integrated modeling study of STAT3, a member of the STAT family of utmost importance in cancer development and therapy, in which we combine available experimental data with several computational methodologies such as homology modeling, protein-protein docking, and molecular dynamics to build reliable atomistic models of USTAT3 dimers. The models generated with the integrative approach presented here were then validated by performing computational alanine scanning for all the residues in the protein-protein interface. These results confirmed the experimental observation of the importance of some of these residues (in particular Leu78 and Asp19) in the USTAT3 dimerization process. Given the growing importance of USTAT3 dimers in several cellular pathways, our models provide an important tool for studying the effects of pathological mutations at the molecular and/or atomistic level, and in the rational design of new inhibitors of dimerization.

Molecular Dynamics of Biomolecules through Direct Analysis of Dipolar Couplings.

Residual dipolar couplings (RDCs) are important probes in structural biology, but their analysis is often complicated by the determination of an alignment tensor or its associated assumptions. We here apply the maximum entropy principle to derive a tensor-free formalism which allows for direct, dynamic analysis of RDCs and holds the classic tensor formalism as a special case. Specifically, the framework enables us to robustly analyze data regardless of whether a clear separation of internal and overall dynamics is possible. Such a separation is often difficult in the core subjects of current structural biology, which include multidomain and intrinsically disordered proteins as well as nucleic acids. We demonstrate the method is tractable and self-consistent and generalizes to data sets comprised of observations from multiple different alignment conditions.

Role of Water in the Puzzling Mechanism of the Final Aromatization Step Promoted by the Human Aromatase Enzyme. Insights from QM/MM MD Simulations.

The enzyme human aromatase (HA) catalyzes the conversion of androgens to estrogens via two hydroxylation reactions and a final unique aromatization step. Despite the great interest of HA as a drug target against breast cancer detailed structural and spectroscopic information on this enzyme became available only in the past few years. As such, the enigmatic mechanism of the final aromatization step is still a matter of debate. Here, we investigated the final step of the HA enzymatic cycle via hybrid quantum-classical (QM/MM) metadynamics and blue-moon ensemble simulations. Our results show that the rate-determining step of the aromatization process is the nucleophilic attack of the distal oxygen of a peroxo-ferric species on the formyl carbon of the enol-19-oxo-androstenedione, which occurs with a free energy barrier (ΔF(#)) of ∼ 16.7 ± 1.9 kcal/mol, in good agreement with experimental data. This reaction is followed by a water mediated 1ß-hydrogen abstraction (ΔF(#) = 7.9 ± 0.8 kcal/mol) and by the formation of a hydroxo-ferric moiety. This latter may be finally protonated by a hydrogen delivery channel involving Asp309 and Thr310, both residues pointed out as crucial for HA activity. In the absence of the catalytic water in the active site the substrate does not assume a position suitable to undergo the nucleophilic attack. Our data not only reveal a novel possible mechanism for the aromatization process consistent with some of the spectroscopic and kinetic data available in the literature, complementing current knowledge on the mechanism of this enzyme, but also point out a remarkable influence of the level of theory used on the calculated free energy barriers. The structural information obtained in this study may be used for the rational structure-based drug design of HA inhibitors to be employed in breast cancer therapy.

Allosteric Modulation of Alpha7 Nicotinic Receptors: Mechanistic Insight through Metadynamics and Essential Dynamics.

Increasing attention has recently been devoted to allosteric modulators, as they can provide inherent advantages over classic receptor agonists. In the field of nicotinic receptors (nAChRs), the main advantage is that allosteric modulators can trigger pharmacological responses, limiting receptor desensitization. Most of the known allosteric ligands are "positive allosteric modulators" (PAMs), which increase both sensitivity to receptor agonists and current amplitude. Intriguingly, some allosteric modulators are also able to activate the α7 receptor (α7-nAChR) even in the absence of orthosteric agonists. These compounds have been named "ago-allosteric modulators" and GAT107 has been studied in depth because of its unique mechanism of action. We here investigate by molecular dynamics simulations, metadynamics, and essential dynamics the activation mechanism of α7-nAChR, in the presence of different nicotinic modulators. We determine the free energy profiles associated with the closed-to-open motion of the loop C, and we highlight mechanistic differences observed in the presence of different modulators. In particular, we demonstrate that GAT107 triggers conformational motions and cross-talk similar to those observed when the α7-nACh receptor is in complex with both an agonist and an allosteric modulator.

Covalent docking of selected boron-based serine beta-lactamase inhibitors.

AmpC ß-lactamase is a hydrolytic enzyme conferring resistance to ß-lactam antibiotics in multiple Gram-negative bacteria. Therefore, identification of non-ß-lactam compounds able to inhibit the enzyme is crucial for the development of novel antibacterial therapies. In general, AmpC inhibitors have to engage the highly solvent-exposed catalytic site of the enzyme. Therefore, understanding the implications of ligand-protein induced-fit and water-mediated interactions behind the inhibitor-enzyme recognition process is fundamental for undertaking structure-based drug design process. Here, we focus on boronic acids, a promising class of beta-lactamase covalent inhibitors. First, we optimized a docking protocol able to reproduce the experimentally determined binding mode of AmpC inhibitors bearing a boronic group. This goal was pursued (1) performing rigid and flexible docking calculations aiming to establish the role of the side chain conformations; and (2) investigating the role of specific water molecules in shaping the enzyme active site and mediating ligand protein interactions. Our calculations showed that some water molecules, conserved in the majority of the considered X-ray structures, are needed to correctly predict the binding pose of known covalent AmpC inhibitors. On this basis, we formalized our findings in a docking and scoring protocol that could be useful for the structure-based design of new boronic acid AmpC inhibitors.

Inactivation of TEM-1 by avibactam (NXL-104): insights from quantum mechanics/molecular mechanics metadynamics simulations.

The fast and constant development of drug-resistant bacteria represents a serious medical emergence. To overcome this problem, the development of drugs with new structures and modes of action is urgently needed. In this context, avibactam represents a promising, innovative inhibitor of beta-lactamases with a novel molecular structure compared to previously developed inhibitors, showing a promising inhibitory activity toward a significant number of beta-lactamase enzymes. In this work, we studied, at the atomistic level, the mechanisms of formation of the covalent complex between avibactam and TEM-1, an experimentally well-characterized class A beta-lactamase, using classical and quantum mechanics/molecular mechanics (QM/MM) simulations combined with metadynamics. Our simulations provide a detailed structural and energetic picture of the molecular steps leading to the formation of the avibactam/TEM-1 covalent adduct. In particular, they support a mechanism in which the rate-determining step is the water-assisted Glu166 deprotonation by Ser70. In this mechanistic framework, the predicted activation energy is in good agreement with experimental kinetic measurements. Additionally, our simulations highlight the important role of Lys73 in assisting the Ser70 and Ser130 deprotonations. While based on the specific case of the avibactam/TEM-1, the simple protocol we present here can be immediately extended and applied to the study of covalent complex formation in different enzyme-inhibitor pairs.

Computational approaches elucidate the allosteric mechanism of human aromatase inhibition: a novel possible route to Small-molecule regulation of CYP450s activities?

Human aromatase (HA) is a P450 cytochrome (CYP) with an essential role in estrogen biosynthesis. Since more than 70% of breast cancers are positive for estrogenic receptor (ER), the reduction of estrogen physiological concentrations through HA inhibition is one of most important therapeutic strategies against this cancer type. Recently, experimental evidence showed that selected taxmoxifen metabolites, which are typically used as estrogen receptor modulators (SERMs), inhibit HA through an allosteric mechanism. In this work, we present a computational protocol to (i) characterize the structural framework and (ii) define the atomistic details of the determinants for the noncompetitive inhibition mechanism. Our calculations identify two putative binding sites able to efficiently bind all tamoxifen metabolites. Analysis of long-scale molecular dynamics simulations reveal that endoxifen, the most effective noncompetitive inhibitor, induces significant enzyme rigidity by binding in one of the possible peripheral sites. The consequence of this binding event is the suppression of one of the functional enzymatic collective motions associated with breathing of the substrate access channel. Moreover, an internal dynamics-based alignment of HA with six other human cytochromes shows that this collective motion is common to other members of the CYP450 protein family. On this basis, our findings may thus be of help for the development of new (pan)inhibitors for the therapeutic treatment of cancer, targeting and modulating the activity of HA and of estrogen receptor, and may also stimulate the development of new drug design strategies for chemoprevention and chemoprotection via allosteric inhibition of CYP450 proteins.

Targeting senescence induced by age or chemotherapy with a polyphenol-rich natural extract improves longevity and healthspan in mice.

Accumulating senescent cells within tissues contribute to the progression of aging and age-related diseases. Botanical extracts, rich in phytoconstituents, present a useful resource for discovering therapies that could target senescence and thus improve healthspan. Here, we show that daily oral administration of a standardized extract of Salvia haenkei (Haenkenium (HK)) extended lifespan and healthspan of naturally aged mice. HK treatment inhibited age-induced inflammation, fibrosis and senescence markers across several tissues, as well as increased muscle strength and fur thickness compared with age-matched controls. We also found that HK treatment reduced acutely induced senescence by the chemotherapeutic agent doxorubicin, using p16LUC reporter mice. We profiled the constituent components of HK by mass spectrometry, and identified luteolin-the most concentrated flavonoid in HK-as a senomorphic compound. Mechanistically, by performing surface plasmon resonance and in situ proximity ligation assay, we found that luteolin disrupted the p16-CDK6 interaction. This work demonstrates that administration of HK promotes longevity in mice, possibly by modulating cellular senescence and by disrupting the p16-CDK6 interaction.

Structural and dynamic properties of monoclonal antibodies immobilized on CNTs: a computational study.

Due to the widespread application of carbon nanotube (CNT)-based materials in nanomedicine, it is nowadays of paramount importance to unravel at the atomistic level of detail the structural properties of such bioconjugates in order to rationalize and predict the effect exerted by the graphitic framework on the bio-active counterpart. In this paper, we report for the first time all-atom explicit solvent molecular dynamics (MD) simulations investigating the structural and dynamic properties of a noncovalent bioconjugate in which the monoclonal Cetuximab antibody (Ctx) is adsorbed on a CNT surface. Upon selection of the three most representative adsorption modes as obtained by docking studies, force-field MD and DFT simulations unambiguously showed that hydrophobic interactions mainly govern the adsorption of the protein on the graphitic surface. Two main adsorption poses have been predicted: a pose-fab (p-fab) and pose-fc (p-fc) (fab = fragment antigen binding region; fc = fragment crystallizable region), the former being favored with small-diameter tubes (≤40 Å). In all the predicted poses, the secondary structure of Ctx is largely unaffected by the presence of the graphitic surface and, consistently with previous literature studies, our simulations reveal that positively charged amino acidic residues, such as Lys and Arg, predominantly contribute to the stabilization of the CNT⋅Ctx complex acting like surfactants. The predicted structural models are consistent with the experimental data, for which the immobilization of the antibody on CNTs does not disrupt the structural and recognition properties of the Ctx, consequently supporting the reliability of the used bioconjugation strategy for engineering stable and responsive hybrid nanomaterials for therapeutic applications. Moreover, a remarkable structural similarity of Ctx with antibodies of different isotypes suggests that in principle the CNT framework can interact in the same manner with all antibodies currently used in clinical applications.

Structural role of uracil DNA glycosylase for the recognition of uracil in DNA duplexes. Clues from atomistic simulations.

In the first stage of the base excision repair pathway the enzyme uracil DNA glycosylase (UNG) recognizes and excises uracil (U) from DNA filaments. U repair is believed to occur via a multistep base-flipping process, through which the damaged U base is initially detected and then engulfed into the enzyme active site, where it is cleaved. The subtle recognition mechanism by which UNG discriminates between U and the other similar pyrimidine nucleobases is still a matter of active debate. Detailed structural information on the different steps of the base-flipping pathway may provide insights on it. However, to date only two intermediates have been trapped crystallographically thanks to chemical modifications of the target and/or of its complementary base. Here, we performed force-field based molecular dynamics (MD) simulations to explore the structural and dynamical properties of distinct UNG/dsDNA adducts, containing A:U, A:T, G:U, or G:C base pairs, at different stages of the base-flipping pathway. Our simulations reveal that if U is present in the DNA sequence a short-lived extra-helical (EH) intermediate exists. This is stabilized by a water-mediated H-bond network, which connects U with His148, a residue pointed out by mutational studies to play a key role for U recognition and catalysis. Moreover, in this EH intermediate, UNG induces a remarkable overall axis bend to DNA. We believe this aspect may facilitate the flipping of U, with respect to other similar nucleobases, in the latter part of the base-extrusion process. In fact, a large DNA bend has been demonstrated to be associated with a lowering of the free energy barrier for base-flipping. A detailed comparison of our results with partially flipped intermediates identified crystallographically or computationally for other base-flipping enzymes allows us to validate our results and to formulate hypothesis on the recognition mechanism of UNG. Our study provides a first ground for a detailed understanding of the UNG repair pathway, which is necessary to devise new pharmaceutical strategies for targeting DNA-related pathologies.