Exploring Proprotein Convertase Subtilisin/Kexin 9 (PCSK9) Autoproteolysis Process by Molecular Simulations: Hints for Drug Design

Proprotein convertase subtilisin/kexin 9 (PCSK9) is a notable target for the treatment of hypercholesterolemia because it regulates the population of the low-density lipoprotein receptor (LDLR) on liver cells. The PCSK9 zymogen is a serine protease that spontaneously undergoes a double self-cleavage step. The available X-ray structures depict the PCSK9 mature state, but the atomic details of the zymogen state of the enzyme are still unknown. Additionally, why the protease activity of PCSK9 is blocked after the second autoprocessing step remains unclear, as this deviates from other members of the PCSK family. By performing constant-pH molecular dynamics (MD) simulations, we investigated the protonation state of the catalytic triad of PCSK9 and found that it strongly influences the catalytic properties of the enzyme. Moreover, we determined the final step of the maturation process by classical and steered MD simulations. This study could facilitate the identification of ligands capable of interfering with the PCSK9 maturation process

Biological Characterization of Computationally Designed Analogs of peptide TVFTSWEEYLDWV (Pep2-8) with Increased PCSK9 Antagonistic Activity

The inhibition of the PCSK9/LDLR protein-protein interaction (PPI) is a promising strategy for developing new hypocholesterolemic agents. Recently, new antibodies have been approved for therapy, but the high cost and low patients\u2019 compliance stimulate the development of alternatives. Starting from the structural information available for the complex between PCSK9 and TVFTSWEEYLDWV (Pep2-8) peptide inhibitor and using computational methods, in this work we identified two Pep2-8 analogs as potential inhibitors of the PCSK9/LDLR PPI. Their biological characterization confirmed the theoretical outcomes. Remarkably, the treatment of HepG2 cells with these peptides increased the LDLR protein level on the cellular membrane, with activities that were 100 and 50 times better than the one of Pep2-8 tested at a 50 \u3bcM concentration. Moreover, they were 50 and 5 times more active than Pep2-8 in improving the functional ability of HepG2 cells to uptake extracellular LDL

A Computational Assay of Estrogen Receptor alpha Antagonists Reveals the Key Common Structural Traits of Drugs Effectively Fighting Refractory Breast Cancers

Somatic mutations of the Estrogen Receptor alpha (ER alpha) occur with an up to 40% incidence in ER sensitive breast cancer (BC) patients undergoing prolonged endocrine treatments. These polymorphisms are implicated in acquired resistance, disease relapse, and increased mortality rates, hence representing a current major clinical challenge. Here, multi-microseconds (12.5 mu s) molecular dynamics simulations revealed that recurrent ER alpha. polymorphisms (i.e. L536Q, Y5375, Y537N, D538G) (mER alpha) are constitutively active in their apo form and that they prompt the selection of an agonist (active)-like conformation even upon antagonists binding. Interestingly, our simulations rationalize, for thefirst time, the efficacy profile of (pre)clinically used Selective Estrogen Receptor Modulators/Downregulators (SERMs/SERDs) against these variants, enlightening, at atomistic level of detail, the key common structural traits needed by drugs able to effectively fight refractory BC types. This knowledge represents a key advancement for mechanism-based therapeutics targeting resistant ER alpha isoforms, potentially allowing the community to move a step closer to 'precision medicine' calibrated on patients' genetic profiles and disease progression

A candidate ion-retaining state in the inward-facing conformation of sodium/galactose symporter: Clues from atomistic simulations

The recent Vibrio parahaemolyticus sodium/galactose (vSGLT) symporter crystal structure captures the protein in an inward-facing substrate-bound conformation, with the sodium ion placed, by structural alignment, in a site equivalent to the Na2 site of the leucine transporter (LeuT). A recent study, based on molecular dynamics simulations, showed that the sodium ion spontaneously leaves its initial position diffusing outside vSGLT, toward the intracellular space. This suggested that the crystal structure corresponds to an ion-releasing state of the transporter. Here, using metadynamics, we identified a more stable Na+ binding site corresponding to a putative ion-retaining state of the transporter. In addition, our simulations, consistently with mutagenesis studies, highlight the importance of D189 that, without being one of the NA(+)-coordinating residues, regulates its binding/release

Pharmacological Comparative Characterization of REL-1017 (Esmethadone-HCl) and Other NMDAR Channel Blockers in Human Heterodimeric N-Methyl-D-Aspartate Receptors

Excessive Ca2+ currents via N-methyl-D-aspartate receptors (NMDARs) have been implicated in many disorders. Uncompetitive NMDAR channel blockers are an emerging class of drugs in clinical use for major depressive disorder (MDD) and other neuropsychiatric diseases. The pharmacological characterization of uncompetitive NMDAR blockers in clinical use may improve our understanding of NMDAR function in physiology and pathology. REL-1017 (esmethadone-HCl), a novel uncompetitive NMDAR channel blocker in Phase 3 trials for the treatment of MDD, was characterized together with dextromethorphan, memantine, (±)-ketamine, and MK-801 in cell lines over-expressing NMDAR subtypes using fluorometric imaging plate reader (FLIPR), automated patch-clamp, and manual patch-clamp electrophysiology. In the absence of Mg2+, NMDAR subtypes NR1-2D were most sensitive to low, sub-μM glutamate concentrations in FLIPR experiments. FLIPR Ca2+ determination demonstrated low μM affinity of REL-1017 at NMDARs with minimal subtype preference. In automated and manual patch-clamp electrophysiological experiments, REL-1017 exhibited preference for the NR1-2D NMDAR subtype in the presence of 1 mM Mg2+ and 1 μM L-glutamate. Tau off and trapping characteristics were similar for (±)-ketamine and REL-1017. Results of radioligand binding assays in rat cortical neurons correlated with the estimated affinities obtained in FLIPR assays and in automated and manual patch-clamp assays. In silico studies of NMDARs in closed and open conformation indicate that REL-1017 has a higher preference for docking and undocking the open-channel conformation compared to ketamine. In conclusion, the pharmacological characteristics of REL-1017 at NMDARs, including relatively low affinity at the NMDAR, NR1-2D subtype preference in the presence of 1 mM Mg2+, tau off and degree of trapping similar to (±)-ketamine, and preferential docking and undocking of the open NMDAR, could all be important variables for understanding the rapid-onset antidepressant effects of REL-1017 without psychotomimetic side effects

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

Compartmentalized activities of the pyruvate dehydrogenase complex sustain lipogenesis in prostate cancer.

The mechanisms by which mitochondrial metabolism supports cancer anabolism remain unclear. Here, we found that genetic and pharmacological inactivation of pyruvate dehydrogenase A1 (PDHA1), a subunit of the pyruvate dehydrogenase complex (PDC), inhibits prostate cancer development in mouse and human xenograft tumor models by affecting lipid biosynthesis. Mechanistically, we show that in prostate cancer, PDC localizes in both the mitochondria and the nucleus. Whereas nuclear PDC controls the expression of sterol regulatory element-binding transcription factor (SREBF)-target genes by mediating histone acetylation, mitochondrial PDC provides cytosolic citrate for lipid synthesis in a coordinated manner, thereby sustaining anabolism. Additionally, we found that PDHA1 and the PDC activator pyruvate dehydrogenase phosphatase 1 (PDP1) are frequently amplified and overexpressed at both the gene and protein levels in prostate tumors. Together, these findings demonstrate that both mitochondrial and nuclear PDC sustain prostate tumorigenesis by controlling lipid biosynthesis, thus suggesting this complex as a potential target for cancer therapy

QM/MM MD simulations on the enzymatic pathway of the human flap endonuclease (hFEN1) elucidate common cleavage pathways to RNase H enzymes

Flap endonucleases (FENs) are nucleic acid hydrolyzing enzymes in charge of excising S'-small DNA and RNA fragments (flaps) protruding from nucleic acid structures during the lagging strand DNA replication or the long-patch base excision repair (LP-BER) processes. In this work we report, for the first time, an atomistic and energetic rendering of the enzymatic catalysis promoted by the human FEN1. After reconstruction of a reactive hFEN/double strand (ds) DNA adduct we employed mixed quantum-classical (QM/MM) metadynamics and umbrella sampling free energy calculations, with the QM part treated with the AM1/d-PhoT Hamiltonian, to perform an extensive characterization of all possible reaction pathways underlying the enzymatic cycle. Our extensive investigation points to a most likely reaction pathway very similar to that recently proposed for ribonuclease H, in which the rate determining step is the nucleophilic attack of a water to the scissile phosphate, which occurs concomitantly with its activation by the pro-Rp oxygen of the nucleobase flanking the scissile phosphate. This step requires a free energy barrier in good agreement with experimental data (Delta G(exp)double dagger = 16.1 kcal/mol vs Delta F-calc double dagger = 16 +/- 2 kcal/mol). Due to the important role of FENs in maintaining nucleic acid fidelity and cell proliferation, a detailed understanding of its enzymatic mechanism has broad interest to elucidate a key enzymatic biological process for preserving genome integrity and has implications for medical and biotechnological applications

First-Principles Modeling of Biological Systems and Structure-Based Drug-Design

Molecular modeling techniques play a relevant role in drug design providing detailed information at atomistic level on the structural, dynamical, mechanistic and electronic properties of biological systems involved in diseases\u2019 onset, integrating and supporting commonly used experimental approaches. These information are often not accessible to the experimental techniques taken singularly, but are of crucial importance for drug design. Due to the enormous increase of the computer power in the last decades, quantum mechanical (QM) or first-principles-based methods have become often used to address biological issues of pharmaceutical relevance, providing relevant information for drug design. Due to their complexity and their size biological systems are often investigated by means of a mixed quantum-classical (QM/MM) approach, which treats at an accurate QM level a limited chemically relevant portion of the system and at the molecular mechanics (MM) level the remaining of the biomolecule and its environment. This method provides a good compromise between computational cost and accuracy, allowing to characterize the properties of the biological system and the (free) energy landscape of the process in study with the accuracy of a QM description. In this review, after a brief introduction of QM and QM/MM methods, we will discuss few representative examples, taken from our work, of the application of these methods in the study of metallo-enzymes of pharmaceutical interest, of metal-containing anticancer drugs targeting the DNA as well as of neurodegenerative diseases. The information obtained from these studies may provide the basis for a rationale structure-based drug design of new and more efficient inhibitors or drugs

Influence of the membrane lipophilic environment on the structure and on the substrate access/egress routes of the human aromatase enzyme. A computational study

Human aromatase (HA), an enzyme located on the membrane of the endoplasmatic reticulum, is of crucial biological importance in the biosynthesis of estrogens. High levels of estrogens are related with important pathologies, conferring to HA a key role as a pharmacological target. In this study we provide, for the first time, an atomistic model of HA embedded on a membrane model to understand the influence of the membrane lipophilic environment on the structural and dynamical properties of HA and on the access/egress pathways of the substrate (androstenedione, ASD) and of the oxygen molecule (involved in the enzymatic process) into/from the HA active site. To this end we used several computational techniques such as force field-based molecular dynamics (MD) simulations, Random Expulsion MD, Steered MD, and Implicit Ligand Sampling. Our results show that the membrane anchoring does not markedly affect the structural properties and the flexibility of the protein, but they clearly point out that the membrane has a marked effect on the access/egress routes of the reactants, stabilizing the formation of different channels for both ASD and O-2 with respect to those observed in pure water solution. Due to the importance of HA in medicine and since access/egress channels may influence its substrate selectivity, a detailed understanding of the role of the membrane in shaping these channels may be of valuable help in drug design