1
|
De Sciscio M, Nardi AN, Centola F, Rossi M, Guarnera E, D’Abramo M. Molecular Modeling of the Deamidation Reaction in Solution: A Theoretical-Computational Study. J Phys Chem B 2023; 127:9550-9559. [PMID: 37903302 PMCID: PMC10641835 DOI: 10.1021/acs.jpcb.3c04662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Revised: 09/15/2023] [Accepted: 10/06/2023] [Indexed: 11/01/2023]
Abstract
In this work, a theoretical-computational method is applied to study the deamidation reaction, a critical post-translational modification in proteins, using a simple model molecule in solution. The method allows one to comprehensively address the environmental effect, thereby enabling one to accurately derive the kinetic rate constants for the three main steps of the deamidation process. The results presented, in rather good agreement with the available experimental data, underline the necessity for a rigorous treatment of environmental factors and a precise kinetic model to correctly assess the overall kinetics of the deamidation reaction.
Collapse
Affiliation(s)
- Maria
Laura De Sciscio
- Department
of Chemistry, University of Rome, Sapienza, P.le A. Moro 5, 00185 Rome, Italy
| | | | - Fabio Centola
- Global
Analytical Development, Merck Serono S.p.A., 00012 Guidonia Montecelio, Italy
| | - Mara Rossi
- Global
Analytical Development, Merck Serono S.p.A., 00012 Guidonia Montecelio, Italy
| | - Enrico Guarnera
- Global
Analytical Development, Merck Serono S.p.A., 00012 Guidonia Montecelio, Italy
- Antibody
Discovery and Protein Engineering, Merck
Healthcare KGaA, 64293 Darmstadt, Germany
| | - Marco D’Abramo
- Department
of Chemistry, University of Rome, Sapienza, P.le A. Moro 5, 00185 Rome, Italy
| |
Collapse
|
2
|
Fernández-Quintero ML, Ljungars A, Waibl F, Greiff V, Andersen JT, Gjølberg TT, Jenkins TP, Voldborg BG, Grav LM, Kumar S, Georges G, Kettenberger H, Liedl KR, Tessier PM, McCafferty J, Laustsen AH. Assessing developability early in the discovery process for novel biologics. MAbs 2023; 15:2171248. [PMID: 36823021 PMCID: PMC9980699 DOI: 10.1080/19420862.2023.2171248] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 01/18/2023] [Indexed: 02/25/2023] Open
Abstract
Beyond potency, a good developability profile is a key attribute of a biological drug. Selecting and screening for such attributes early in the drug development process can save resources and avoid costly late-stage failures. Here, we review some of the most important developability properties that can be assessed early on for biologics. These include the influence of the source of the biologic, its biophysical and pharmacokinetic properties, and how well it can be expressed recombinantly. We furthermore present in silico, in vitro, and in vivo methods and techniques that can be exploited at different stages of the discovery process to identify molecules with liabilities and thereby facilitate the selection of the most optimal drug leads. Finally, we reflect on the most relevant developability parameters for injectable versus orally delivered biologics and provide an outlook toward what general trends are expected to rise in the development of biologics.
Collapse
Affiliation(s)
- Monica L. Fernández-Quintero
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Anne Ljungars
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Franz Waibl
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Victor Greiff
- Department of Immunology, University of Oslo, Oslo, Norway
| | - Jan Terje Andersen
- Department of Immunology, University of Oslo, Oslo University Hospital Rikshospitalet, Oslo, Norway
- Institute of Clinical Medicine and Department of Pharmacology, University of Oslo, Oslo, Norway
| | | | - Timothy P. Jenkins
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Bjørn Gunnar Voldborg
- National Biologics Facility, Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Lise Marie Grav
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Sandeep Kumar
- Biotherapeutics Discovery, Boehringer Ingelheim Pharmaceuticals Inc, Ridgefield, CT, USA
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Hubert Kettenberger
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Klaus R. Liedl
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Peter M. Tessier
- Department of Chemical Engineering, Pharmaceutical Sciences and Biomedical Engineering, Biointerfaces Institute, University of Michigan, Ann Arbor, Michigan, USA
| | - John McCafferty
- Department of Medicine, Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge, Cambridge, UK
- Maxion Therapeutics, Babraham Research Campus, Cambridge, UK
| | - Andreas H. Laustsen
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| |
Collapse
|
3
|
Abstract
Monoclonal antibodies are susceptible to chemical and enzymatic modifications during manufacturing, storage, and shipping. Deamidation, isomerization, and oxidation can compromise the potency, efficacy, and safety of therapeutic antibodies. Recently, in silico tools have been used to identify liable residues and engineer antibodies with better chemical stability. Computational approaches for predicting deamidation, isomerization, oxidation, glycation, carbonylation, sulfation, and hydroxylation are reviewed here. Although liable motifs have been used to improve the chemical stability of antibodies, the accuracy of in silico predictions can be improved using machine learning and molecular dynamic simulations. In addition, there are opportunities to improve predictions for specific stress conditions, develop in silico prediction of novel modifications in antibodies, and predict the impact of modifications on physical stability and antigen-binding.
Collapse
Affiliation(s)
- Shabdita Vatsa
- Development Services, Lonza Biologics, Singapore, Singapore
| |
Collapse
|
4
|
Lawson KE, Dekle JK, Evans MN, Adamczyk AJ. Deamidation reaction network mapping of pharmacologic and related proteins: impact of solvation dielectric on the degradation energetics of asparagine dipeptides. REACT CHEM ENG 2022. [DOI: 10.1039/d2re00110a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Asn-X deamidation pathways in the FV region of the monoclonal antibody (mAb).
Collapse
Affiliation(s)
| | - Joseph K. Dekle
- Department of Chemical Engineering, Auburn University, Auburn, AL, USA
| | - Megan N. Evans
- Department of Chemical Engineering, Auburn University, Auburn, AL, USA
| | | |
Collapse
|
5
|
Exploration of tautomerizations of succinimide and maleimide assisted by ammonia and methanol: a theoretical perspective. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2446-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
6
|
Shalaev E, Soper A, Zeitler JA, Ohtake S, Roberts CJ, Pikal MJ, Wu K, Boldyreva E. Freezing of Aqueous Solutions and Chemical Stability of Amorphous Pharmaceuticals: Water Clusters Hypothesis. J Pharm Sci 2018; 108:36-49. [PMID: 30055227 DOI: 10.1016/j.xphs.2018.07.018] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2018] [Revised: 07/13/2018] [Accepted: 07/17/2018] [Indexed: 11/24/2022]
Abstract
Molecular mobility has been traditionally invoked to explain physical and chemical stability of diverse pharmaceutical systems. Although the molecular mobility concept has been credited with creating a scientific basis for stabilization of amorphous pharmaceuticals and biopharmaceuticals, it has become increasingly clear that this approach represents only a partial description of the underlying fundamental principles. An additional mechanism is proposed herein to address 2 key questions: (1) the existence of unfrozen water (i.e., partial or complete freezing inhibition) in aqueous solutions at subzero temperatures and (2) the role of water in the chemical stability of amorphous pharmaceuticals. These apparently distant phenomena are linked via the concept of water clusters. In particular, freezing inhibition is associated with the confinement of water clusters in a solidified matrix of an amorphous solute, with nanoscaled water clusters being observed in aqueous glasses using wide-angle neutron scattering. The chemical instability is suggested to be directly related to the catalysis of proton transfer by water clusters, considering that proton transfer is the key elementary reaction in many chemical processes, including such common reactions as hydrolysis and deamidation.
Collapse
Affiliation(s)
- Evgenyi Shalaev
- Pharmaceutical Development, Allergan plc., Irvine, California 92612.
| | - Alan Soper
- ISIS Facility, UKRI-STFC Rutherford Appleton Laboratory, Harwell Campus, Didcot, Oxon OX11 OQX, UK
| | - J Axel Zeitler
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, UK
| | - Satoshi Ohtake
- Pfizer BioTherapeutics Pharmaceutical Sciences, Chesterfield, Missouri 63198
| | | | - Michael J Pikal
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Connecticut, Storrs, Connecticut 06269
| | - Ke Wu
- Pharmaceutical Development, Allergan plc., Irvine, California 92612
| | - Elena Boldyreva
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090, Russian Federation; Boreskov Institute of Catalysis SB RAS, Novosibirsk 630090, Russian Federation
| |
Collapse
|
7
|
Herr JE, Yao K, McIntyre R, Toth DW, Parkhill J. Metadynamics for training neural network model chemistries: A competitive assessment. J Chem Phys 2018; 148:241710. [PMID: 29960377 DOI: 10.1063/1.5020067] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Neural network model chemistries (NNMCs) promise to facilitate the accurate exploration of chemical space and simulation of large reactive systems. One important path to improving these models is to add layers of physical detail, especially long-range forces. At short range, however, these models are data driven and data limited. Little is systematically known about how data should be sampled, and "test data" chosen randomly from some sampling techniques can provide poor information about generality. If the sampling method is narrow, "test error" can appear encouragingly tiny while the model fails catastrophically elsewhere. In this manuscript, we competitively evaluate two common sampling methods: molecular dynamics (MD), normal-mode sampling, and one uncommon alternative, Metadynamics (MetaMD), for preparing training geometries. We show that MD is an inefficient sampling method in the sense that additional samples do not improve generality. We also show that MetaMD is easily implemented in any NNMC software package with cost that scales linearly with the number of atoms in a sample molecule. MetaMD is a black-box way to ensure samples always reach out to new regions of chemical space, while remaining relevant to chemistry near kbT. It is a cheap tool to address the issue of generalization.
Collapse
Affiliation(s)
- John E Herr
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - Ryker McIntyre
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - David W Toth
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - John Parkhill
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| |
Collapse
|
8
|
Theoretical analysis of tautomerization of succinimide and analogous compounds: insights from DFT approach. Struct Chem 2018. [DOI: 10.1007/s11224-018-1075-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
9
|
Ohtake S, Feng S, Shalaev E. Effect of Water on the Chemical Stability of Amorphous Pharmaceuticals: 2. Deamidation of Peptides and Proteins. J Pharm Sci 2018; 107:42-56. [DOI: 10.1016/j.xphs.2017.09.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2017] [Accepted: 09/11/2017] [Indexed: 01/14/2023]
|
10
|
Plotnikov NV, Singh SK, Rouse JC, Kumar S. Quantifying the Risks of Asparagine Deamidation and Aspartate Isomerization in Biopharmaceuticals by Computing Reaction Free-Energy Surfaces. J Phys Chem B 2017; 121:719-730. [DOI: 10.1021/acs.jpcb.6b11614] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Nikolay V. Plotnikov
- Pharmaceutical
Research and Development, Biotherapeutics Pharmaceutical Sciences, Pfizer Inc., 700 Chesterfield Pkwy West, Chesterfield, Missouri 63017, United States
| | - Satish Kumar Singh
- Pharmaceutical
Research and Development, Biotherapeutics Pharmaceutical Sciences, Pfizer Inc., 700 Chesterfield Pkwy West, Chesterfield, Missouri 63017, United States
| | - Jason C. Rouse
- Analytical
Research and Development, Biotherapeutics Pharmaceutical Sciences, Pfizer Inc., 1 Burtt Road, Andover, Massachusetts 01810, United States
| | - Sandeep Kumar
- Pharmaceutical
Research and Development, Biotherapeutics Pharmaceutical Sciences, Pfizer Inc., 700 Chesterfield Pkwy West, Chesterfield, Missouri 63017, United States
| |
Collapse
|
11
|
Ince HH, Konuklar FAS, Ugur I, Ozcan ÖA, Sayadi M, Feig M, Aviyente V. Role of then+1 amino acid residue on the deamidation of asparagine in pentapeptides. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1068394] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
12
|
|
13
|
Zheng S, Pfaendtner J. Enhanced sampling of chemical and biochemical reactions with metadynamics. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.923574] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
14
|
Halim MA, Almatarneh MH, Poirier RA. Mechanistic Study of the Deamidation Reaction of Glutamine: A Computational Approach. J Phys Chem B 2014; 118:2316-30. [DOI: 10.1021/jp4107266] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Mohammad A. Halim
- Department
of Chemistry, Memorial University, St. John’s, NL, Canada A1B 3X7
| | | | - Raymond A. Poirier
- Department
of Chemistry, Memorial University, St. John’s, NL, Canada A1B 3X7
| |
Collapse
|
15
|
Selvaraj ARK, Murugan NA, Ågren H. Solvent Polarity-Induced Conformational Unlocking of Asparagine. J Phys Chem A 2012; 116:11702-8. [DOI: 10.1021/jp307715n] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Ananda Rama Krishnan Selvaraj
- Division of Theoretical Chemistry
and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - N. Arul Murugan
- Division of Theoretical Chemistry
and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Hans Ågren
- Division of Theoretical Chemistry
and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| |
Collapse
|
16
|
Mujika JI, Lopez X, Mulholland AJ. Mechanism of C-terminal intein cleavage in protein splicing from QM/MM molecular dynamics simulations. Org Biomol Chem 2011; 10:1207-18. [PMID: 22179261 DOI: 10.1039/c1ob06444d] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Protein splicing is a post-translational process in which a biologically inactive protein is activated by the release of a segment denoted as an intein. The process involves four steps. In the third, the scission of the intein takes place after the cyclization of the last amino acid of the segment, an asparagine. Little is known about the chemical reaction necessary for this cyclization. Experiments demonstrate that two histidines (the penultimate amino acid of the intein, and a histidine located 10 amino acids upstream) are relevant in the cyclization of the asparagine. We have investigated the mechanism and determinants of reaction in the GyrA intein focusing on the requirements for asparagine activation for its cyclization. First, the influence that the protonation states of these two histidines have on the orientation of the asparagine side chain is investigated by means of molecular dynamics simulation. Molecular dynamics simulations using the CHARMM27 force field were carried out on the three possible protonation states for each of these two histidines. The results indicate that the only protonation state in which the conformation of the system is suitable for cyclization is when the penultimate histidine is fully protonated (positively charged), and the upstream histidine is in the His(ε) neutral tautomeric form. The free energy profile for the reaction in which the asparagine is activated by a proton transfer to the upstream histidine is presented, computed by hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics at the SCCDFTB/CHARMM27 level of theory. The calculated free energy barrier for the reaction is 19.0 kcal mol(-1). B3LYP/6-31+G(d) QM/MM single-point calculations give a qualitatively a similar energy profile, although with somewhat higher energy barriers, in good agreement with the value derived from experiment of 25 kcal mol(-1) at 60 °C. QM/MM molecular dynamics simulations of the reactant, activated reactant and intermediate states highlight the importance of the Arg181-Val182-Asp183 segment in catalysing the reaction. Overall, the results indicate that nucleophilic activation of the asparagine for its cyclization by the upstream histidine acting as the base is a plausible mechanism for the C-terminal cleavage in protein splicing.
Collapse
Affiliation(s)
- Jon I Mujika
- Kimika Fakultatea, Euskal Herriko Unibertsitatea and Donostia International Physics Center (DIPC, P. K. 1072, 20080, Donostia, Euskadi, Spain.
| | | | | |
Collapse
|
17
|
Pan K, Deem MW. Predicting Fixation Tendencies of the H3N2 Influenza Virus by Free Energy Calculation. J Chem Theory Comput 2011; 7:1259-1272. [PMID: 21691431 PMCID: PMC3117619 DOI: 10.1021/ct100540p] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Influenza virus evolves to escape from immune system antibodies that bind to it. We used free energy calculations with Einstein crystals as reference states to calculate the difference of antibody binding free energy (ΔΔG) induced by amino acid substitution at each position in epitope B of the H3N2 influenza hemagglutinin, the key target for antibody. A substitution with positive ΔΔG value decreases the antibody binding constant. On average an uncharged to charged amino acid substitution generates the highest ΔΔG values. Also on average, substitutions between small amino acids generate ΔΔG values near to zero. The 21 sites in epitope B have varying expected free energy differences for a random substitution. Historical amino acid substitutions in epitope B for the A/Aichi/2/1968 strain of influenza A show that most fixed and temporarily circulating substitutions generate positive ΔΔG values. We propose that the observed pattern of H3N2 virus evolution is affected by the free energy landscape, the mapping from the free energy landscape to virus fitness landscape, and random genetic drift of the virus. Monte Carlo simulations of virus evolution are presented to support this view.
Collapse
Affiliation(s)
- Keyao Pan
- Department of Bioengineering, Rice University, Houston, Texas 77005
| | - Michael W. Deem
- Department of Bioengineering, Rice University, Houston, Texas 77005
- Department of Physics & Astronomy, Rice University, Houston, Texas 77005
| |
Collapse
|
18
|
Barducci A, Bonomi M, Parrinello M. Metadynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.31] [Citation(s) in RCA: 712] [Impact Index Per Article: 54.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Alessandro Barducci
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Lugano, Switzerland
| | - Massimiliano Bonomi
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Lugano, Switzerland
| | - Michele Parrinello
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Lugano, Switzerland
| |
Collapse
|