1
|
Substrates and Cyclic Peptide Inhibitors of the Oligonucleotide-Activated Sirtuin 7. Angew Chem Int Ed Engl 2023; 62:e202314597. [PMID: 37873919 DOI: 10.1002/anie.202314597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 10/19/2023] [Accepted: 10/23/2023] [Indexed: 10/25/2023]
Abstract
The sirtuins are NAD+ -dependent lysine deacylases, comprising seven isoforms (SIRT1-7) in humans, which are involved in the regulation of a plethora of biological processes, including gene expression and metabolism. The sirtuins share a common hydrolytic mechanism but display preferences for different ϵ-N-acyllysine substrates. SIRT7 deacetylates targets in nuclei and nucleoli but remains one of the lesser studied of the seven isoforms, in part due to a lack of chemical tools to specifically probe SIRT7 activity. Here we expressed SIRT7 and, using small-angle X-ray scattering, reveal SIRT7 to be a monomeric enzyme with a low degree of globular flexibility in solution. We developed a fluorogenic assay for investigation of the substrate preferences of SIRT7 and to evaluate compounds that modulate its activity. We report several mechanism-based SIRT7 inhibitors as well as de novo cyclic peptide inhibitors selected from mRNA-display library screening that exhibit selectivity for SIRT7 over other sirtuin isoforms, stabilize SIRT7 in cells, and cause an increase in the acetylation of H3 K18.
Collapse
|
2
|
Combining Molecular Dynamics Simulations and Biophysical Characterization to Investigate Protein-Specific Excipient Effects on Reteplase during Freeze Drying. Pharmaceutics 2023; 15:1854. [PMID: 37514040 PMCID: PMC10384596 DOI: 10.3390/pharmaceutics15071854] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Revised: 06/24/2023] [Accepted: 06/26/2023] [Indexed: 07/30/2023] Open
Abstract
We performed molecular dynamics simulations of Reteplase in the presence of different excipients to study the stabilizing mechanisms and to identify the role of excipients during freeze drying. To simulate the freeze-drying process, we divided the process into five distinct steps: (i) protein-excipient formulations at room temperature, (ii) the ice-growth process, (iii)-(iv) the partially solvated and fully dried formulations, and (v) the reconstitution. Furthermore, coarse-grained (CG) simulations were employed to explore the protein-aggregation process in the presence of arginine. By using a coarse-grained representation, we could observe the collective behavior and interactions between protein molecules during the aggregation process. The CG simulations revealed that the presence of arginine prevented intermolecular interactions of the catalytic domain of Reteplase, thus reducing the aggregation propensity. This suggests that arginine played a stabilizing role by interacting with protein-specific regions. From the freeze-drying simulations, we could identify several protein-specific events: (i) collapse of the domain structure, (ii) recovery of the drying-induced damages during reconstitution, and (iii) stabilization of the local aggregation-prone region via direct interactions with excipients. Complementary to the simulations, we employed nanoDSF, size-exclusion chromatography, and CD spectroscopy to investigate the effect of the freeze-drying process on the protein structure and stability.
Collapse
|
3
|
Chemometrics in Protein Formulation: Stability Governed by Repulsion and Protein Unfolding. Mol Pharm 2023. [PMID: 37146162 DOI: 10.1021/acs.molpharmaceut.3c00013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Therapeutic proteins can be challenging to develop due to their complexity and the requirement of an acceptable formulation to ensure patient safety and efficacy. To date, there is no universal formulation development strategy that can identify optimal formulation conditions for all types of proteins in a fast and reliable manner. In this work, high-throughput characterization, employing a toolbox of five techniques, was performed on 14 structurally different proteins formulated in 6 different buffer conditions and in the presence of 4 different excipients. Multivariate data analysis and chemometrics were used to analyze the data in an unbiased way. First, observed changes in stability were primarily determined by the individual protein. Second, pH and ionic strength are the two most important factors determining the physical stability of proteins, where there exists a significant statistical interaction between protein and pH/ionic strength. Additionally, we developed prediction methods by partial least-squares regression. Colloidal stability indicators are important for prediction of real-time stability, while conformational stability indicators are important for prediction of stability under accelerated stress conditions at 40 °C. In order to predict real-time storage stability, protein-protein repulsion and the initial monomer fraction are the most important properties to monitor.
Collapse
|
4
|
Structural characterization of human tryptophan hydroxylase 2 reveals that L-Phe is superior to L-Trp as the regulatory domain ligand. Structure 2023:S0969-2126(23)00127-2. [PMID: 37119821 DOI: 10.1016/j.str.2023.04.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2022] [Revised: 03/03/2023] [Accepted: 04/04/2023] [Indexed: 05/01/2023]
Abstract
Tryptophan hydroxylase 2 (TPH2) catalyzes the rate-limiting step in serotonin biosynthesis in the brain. Consequently, regulation of TPH2 is relevant for serotonin-related diseases, yet the regulatory mechanism of TPH2 is poorly understood and structural and dynamical insights are missing. We use NMR spectroscopy to determine the structure of a 47 N-terminally truncated variant of the regulatory domain (RD) dimer of human TPH2 in complex with L-Phe, and show that L-Phe is the superior RD ligand compared with the natural substrate, L-Trp. Using cryo-EM, we obtain a low-resolution structure of a similarly truncated variant of the complete tetrameric enzyme with dimerized RDs. The cryo-EM two-dimensional (2D) class averages additionally indicate that the RDs are dynamic in the tetramer and likely exist in a monomer-dimer equilibrium. Our results provide structural information on the RD as an isolated domain and in the TPH2 tetramer, which will facilitate future elucidation of TPH2's regulatory mechanism.
Collapse
|
5
|
Abstract
While α-, β-, and γ-cyclodextrin (CD) are ubiquitous hosts employed by supramolecular chemists, δ-CD (formed from nine α-1,4-linked glucopyranose units) has received very little attention. α-, β-, and γ-CD are the major products of the enzymatic breakdown of starch by cyclodextrin glucanotransferase (CGTase), but δ-CD forms only transiently in this reaction, as a minor component of a complex mixture of linear and cyclic glucans. In this work, we show how δ-CD can be synthesized in unprecedented yields by employing a bolaamphiphile template in an enzyme-mediated dynamic combinatorial library of cyclodextrins. NMR spectroscopy studies revealed that δ-CD can thread up to three bolaamphiphiles forming [2]-, [3]-, or [4]-pseudorotaxanes, depending on the size of the hydrophilic headgroup and the length of the alkyl chain axle. Threading of the first bolaamphiphile occurs in fast exchange on the NMR chemical shift time scale, while subsequent threading occurs in slow exchange. To extract quantitative information for 1:2 and 1:3 binding events occurring in mixed exchange regimes, we derived equations for nonlinear curve fitting that take into consideration both the chemical shift changes for species in fast exchange and the integrals for species in slow exchange to determine Ka1, Ka2, and Ka3. Template T1 could be used to direct the enzymatic synthesis of δ-CD due to the cooperative formation of a 1:2 complex─the [3]-pseudorotaxane δ-CD·T12. Importantly, T1 is recyclable. It can be readily recovered from the enzymatic reaction by precipitation and reused in subsequent syntheses enabling preparative-scale synthesis of δ-CD.
Collapse
|
6
|
Molecular insights into alginate β-lactoglobulin A multivalencies-The foundation for their amorphous aggregates and coacervation. Protein Sci 2023; 32:e4556. [PMID: 36571497 PMCID: PMC9847093 DOI: 10.1002/pro.4556] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Revised: 12/06/2022] [Accepted: 12/22/2022] [Indexed: 12/27/2022]
Abstract
For improved control of biomaterial property design, a better understanding of complex coacervation involving anionic polysaccharides and proteins is needed. Here, we address the initial steps in condensate formation of β-lactoglobulin A (β-LgA) with nine defined alginate oligosaccharides (AOSs) and describe their multivalent interactions in structural detail. Binding of AOSs containing four, five, or six uronic acid residues (UARs), either all mannuronate (M), all guluronate (G), or alternating M and G embodying the block structural components of alginates, was characterized by isothermal titration calorimetry, nuclear magnetic resonance spectroscopy (NMR), and molecular docking. β-LgA was highly multivalent exhibiting binding stoichiometries decreasing from five to two AOSs with increasing degree of polymerization (DP) and similar affinities in the mid micromolar range. The different AOS binding sites on β-LgA were identified by NMR chemical shift perturbation analyses and showed diverse compositions of charged, polar and hydrophobic residues. Distinct sites for the shorter AOSs merged to accommodate longer AOSs. The AOSs bound dynamically to β-LgA, as concluded from saturation transfer difference and 1 H-ligand-targeted NMR analyses. Molecular docking using Glide within the Schrödinger suite 2016-1 revealed the orientation of AOSs to only vary slightly at the preferred β-LgA binding site resulting in similar XP glide scores. The multivalency coupled with highly dynamic AOS binding with lack of confined conformations in the β-LgA complexes may help explain the first steps toward disordered β-LgA alginate coacervate structures.
Collapse
|
7
|
Dynamics of Human Serum Transferrin in Varying Physicochemical Conditions Explored by Using Molecular Dynamics Simulations. Mol Pharm 2022; 19:2795-2806. [PMID: 35776490 DOI: 10.1021/acs.molpharmaceut.2c00158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Conformational stability of human serum transferrin (Tf) at varying pH values and salt and excipient concentrations were investigated using molecular dynamics (MD) simulations, and the results are compared with previously published small-angle X-ray scattering (SAXS) experiments. SAXS study showed that at pH 5, Tf is predominantly present in a partially open (PO) form, and the factions of PO differ based on the physicochemical condition and drift toward the closed form (HO) as the pH increases. Tf is a bilobal glycoprotein that is composed of homologous halves termed the N- and C-lobes. The current study shows that the protonation of Y188 and K206 at pH 5 is the primary conformational drive into PO, which shifts toward the closed (HO) conformer as the pH increases. Furthermore, at pH 6.5, PO is unfavorable due to negative charge-charge repulsion at the N/C-lobe interface linker region causing increased hinge distance when compared to HO, which has favorable attractive electrostatic interactions in this region. Subsequently, the effect of salt concentration was studied at 70 and 140 mM NaCl. At 70 mM NaCl and pH 5, chloride ions bind strongly in the N-lobe iron-binding site, whereas these interactions are weak at pH 6.5. With increasing salt concentration at pH 5, the regions surrounding the N-lobe iron-binding site are saturated, and as a consequence, sodium and chloride ions accumulate into the bulk. Additionally, protein-excipient interactions were investigated. At pH 5, the excipients interact in similar loop regions, E89-T93, and D416-D420, located in the N- and C-lobes of the HO conformer, respectively. It is anticipated that interactions of additives in these two loop regions cause conformational changes that lead to iron-coordinating residues in the N-lobe to drift away from iron and thus drive HO to PO conversion. Furthermore, at pH 6.5 and 140 mM histidine, these interactions are negligible leading to the stabilization of HO.
Collapse
|
8
|
pH- and concentration-dependent supramolecular assembly of a fungal defensin plectasin variant into helical non-amyloid fibrils. Nat Commun 2022; 13:3162. [PMID: 35672293 PMCID: PMC9174238 DOI: 10.1038/s41467-022-30462-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Accepted: 04/28/2022] [Indexed: 02/04/2023] Open
Abstract
Self-assembly and fibril formation play important roles in protein behaviour. Amyloid fibril formation is well-studied due to its role in neurodegenerative diseases and characterized by refolding of the protein into predominantly β-sheet form. However, much less is known about the assembly of proteins into other types of supramolecular structures. Using cryo-electron microscopy at a resolution of 1.97 Å, we show that a triple-mutant of the anti-microbial peptide plectasin, PPI42, assembles into helical non-amyloid fibrils. The in vitro anti-microbial activity was determined and shown to be enhanced compared to the wildtype. Plectasin contains a cysteine-stabilised α-helix-β-sheet structure, which remains intact upon fibril formation. Two protofilaments form a right-handed protein fibril. The fibril formation is reversible and follows sigmoidal kinetics with a pH- and concentration dependent equilibrium between soluble monomer and protein fibril. This high-resolution structure reveals that α/β proteins can natively assemble into fibrils. Here the authors report the cryo-EM structure of a triple-mutant of the anti-microbial peptide plectasin, PPI42, assembling in a pH- and concentration dependent manner into helical non-amyloid fibrils. The fibrils formation is reversible, and follows a sigmoidal kinetics. The fibrils adopt a right-handed helical superstructure composed by two protofilaments, stabilized by an outer hydrophobic ring and an inner hydrophobic centre. These findings reveal that α/β proteins can natively assemble into fibrils.
Collapse
|
9
|
In-silico study of the interactions between acylated glucagon like-peptide-1 analogues and the native receptor. J Biomol Struct Dyn 2022:1-15. [PMID: 35612899 DOI: 10.1080/07391102.2022.2078409] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
We have performed a series of multiple molecular dynamics (MD) simulations of glucagon-like peptide-1 (GLP-1) and acylated GLP-1 analogues in complex with the endogenous receptor (GLP-1R) to obtain a molecular understanding of how fatty acid (FA) chain structure, acylation position on the peptide, and presence of a linker affect the binding. MD simulations were analysed to extract heatmaps of receptor-peptide interaction patterns and to determine the free energy of binding using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach. The extracted free energies from MM-PBSA calculations are in qualitative agreement with experimentally determined potencies. Furthermore, the interaction patterns seen in the receptor-GLP-1 complex simulations resemble previously reported binding interactions validating the simulations. Analysing the receptor-GLP-1 analogue complex simulations, we found that the major differences between the systems stem from FA interactions and positioning of acylation in the peptide. Hydrophobic interactions between the FA chain and a hydrophobic patch on the extracellular domain contribute significantly to the binding affinity. Acylation on Lys26 resulted in noticeably more interactions between the FA chain and the extracellular domain hydrophobic patch than found for acylation on Lys34 and Lys38, respectively. The presence of a charged linker between the peptide and FA chain can potentially stabilise the complex by forming hydrogen bonds to arginine residues in the linker region between the extracellular domain and the transmembrane domain. A molecular understanding of the fatty acid structure and its effect on binding provides important insights into designing acylated agonists for GLP-1R.Communicated by Ramaswamy H. Sarma.
Collapse
|
10
|
Self-Interactions of Two Monoclonal Antibodies: Small-Angle X-ray Scattering, Light Scattering, and Coarse-Grained Modeling. Mol Pharm 2021; 19:508-519. [PMID: 34939811 DOI: 10.1021/acs.molpharmaceut.1c00627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using light scattering (LS), small-angle X-ray scattering (SAXS), and coarse-grained Monte Carlo (MC) simulations, we studied the self-interactions of two monoclonal antibodies (mAbs), PPI03 and PPI13. With LS measurements, we obtained the osmotic second virial coefficient, B22, and the molecular weight, Mw, of the two mAbs, while with SAXS measurements, we studied the mAbs' self-interaction behavior in the high protein concentration regime up to 125 g/L. Through SAXS-derived coarse-grained representations of the mAbs, we performed MC simulations with either a one-protein or a two-protein model to predict B22. By comparing simulation and experimental results, we validated our models and obtained insights into the mAbs' self-interaction properties, highlighting the role of both ion binding and charged patches on the mAb surfaces. Our models provide useful information about mAbs' self-interaction properties and can assist the screening of conditions driving to colloidal stability.
Collapse
|
11
|
Abstract
![]()
Aggregation is a
common phenomenon in the field of protein therapeutics
and can lead to function loss or immunogenic patient responses. Two
strategies are currently used to reduce aggregation: (1) finding a
suitable formulation, which is labor-intensive and requires large
protein quantities, or (2) engineering the protein, which requires
extensive knowledge about the protein aggregation pathway. We present
a biophysical characterization of the oligomerization and aggregation
processes by Interferon alpha-2a (IFNα-2a), a protein drug with
antiviral and immunomodulatory properties. This study combines experimental
high throughput screening with detailed investigations by small-angle
X-ray scattering and analytical ultracentrifugation. Metropolis Monte
Carlo simulations are used to gain insight into the underlying intermolecular
interactions. IFNα-2a forms soluble oligomers that are controlled
by a fast pH and concentration-dependent equilibrium. Close to the
isoelectric point of 6, IFNα-2a forms insoluble aggregates which
can be prevented by adding salt. We show that monomer attraction is driven mainly by molecular anisotropic dipole–dipole interactions
that increase with increasing pH. Repulsion is due
to monopole–monopole interactions and depends on the charge
of IFNα-2a. The study highlights how combining multiple methods
helps to systematically dissect the molecular mechanisms driving oligomer
formation and to design ultimately efficient strategies for preventing
detrimental protein aggregation.
Collapse
|
12
|
Combination of high throughput and structural screening to assess protein stability - a screening perspective. Eur J Pharm Biopharm 2021; 171:1-10. [PMID: 34826593 DOI: 10.1016/j.ejpb.2021.08.018] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 08/23/2021] [Accepted: 08/27/2021] [Indexed: 11/04/2022]
Abstract
High throughput screening to measure the stability of industrially relevant proteins and their variants is necessary for quality assessment in the development process. Advances in automation, measurement time and sample consumption for many techniques allow rapid measurements with minimal amount of protein. However, many methods include automated data analysis, potentially neglecting important aspects of the proteińs behavior in certain conditions. In this study we implement small angle X-ray scattering (SAXS), typically not used to assess protein behavior in industrial screening, in a high throughput screening workflow to address problems of contradicting results and reproducibility among different high throughput methods. As a case study we use the lipases of Thermomyces lanuginosus and Rhizomucor miehei, widely used industrial biocatalysts. We show that even the initial analysis of the SAXS data without performing any time-consuming modelling provide valuable information on interparticle interactions. We conclude that recent advances in automation and data processing, have enabled SAXS to be used more widely as a tool to gain in-depth knowledge highly useful for protein formulation development. This is especially relevant in light of increasing accessibility to SAXS due to the commercial availability of benchtop instruments.
Collapse
|
13
|
Abstract
Enzyme reactions, both in Nature and technical applications, commonly occur at the interface of immiscible phases. Nevertheless, stringent descriptions of interfacial enzyme catalysis remain sparse, and this is partly due to a shortage of coherent experimental data to guide and assess such work. In this work, we produced and kinetically characterized 83 cellulases, which revealed a conspicuous linear free energy relationship (LFER) between the substrate binding strength and the activation barrier. The scaling occurred despite the investigated enzymes being structurally and mechanistically diverse. We suggest that the scaling reflects basic physical restrictions of the hydrolytic process and that evolutionary selection has condensed cellulase phenotypes near the line. One consequence of the LFER is that the activity of a cellulase can be estimated from its substrate binding strength, irrespectively of structural and mechanistic details, and this appears promising for in silico selection and design within this industrially important group of enzymes.
Collapse
|
14
|
Computing Cellulase Kinetics with a Two-Domain Linear Interaction Energy Approach. ACS OMEGA 2021; 6:1547-1555. [PMID: 33490814 PMCID: PMC7818601 DOI: 10.1021/acsomega.0c05361] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2020] [Accepted: 12/24/2020] [Indexed: 05/21/2023]
Abstract
While heterogeneous enzyme reactions play an essential role in both nature and green industries, computational predictions of their catalytic properties remain scarce. Recent experimental work demonstrated the applicability of the Sabatier principle for heterogeneous biocatalysis. This provides a simple relationship between binding strength and the catalytic rate and potentially opens a new way for inexpensive computational determination of kinetic parameters. However, broader implementation of this approach will require fast and reliable prediction of binding free energies of complex two-phase systems, and computational procedures for this are still elusive. Here, we propose a new framework for the assessment of the binding strengths of multidomain proteins, in general, and interfacial enzymes, in particular, based on an extended linear interaction energy (LIE) method. This two-domain LIE (2D-LIE) approach was successfully applied to predict binding and activation free energies of a diverse set of cellulases and resulted in robust models with high accuracy. Overall, our method provides a fast computational screening tool for cellulases that have not been experimentally characterized, and we posit that it may also be applicable to other heterogeneously acting biocatalysts.
Collapse
|
15
|
Small angle X-ray scattering and molecular dynamic simulations provide molecular insight for stability of recombinant human transferrin. JOURNAL OF STRUCTURAL BIOLOGY-X 2020; 4:100017. [PMID: 32647821 PMCID: PMC7337065 DOI: 10.1016/j.yjsbx.2019.100017] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 11/15/2019] [Accepted: 11/25/2019] [Indexed: 12/20/2022]
Abstract
Combination of SAXS and MD simulations can explain the molecular background for the stability studies. Destabilizing excipients in general induce opening of the N-lobe by binding to specific residues in the protein. The fully open conformation of transferrin is not seen without aggregation.
Transferrin is an attractive candidate for drug delivery due to its ability to cross the blood brain barrier. However, in order to be able to use it for therapeutic purposes, it is important to investigate how its stability depends on different formulation conditions. Combining high-throughput thermal and chemical denaturation studies with small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations, it was possible to connect the stability of transferrin with its conformational changes. Lowering pH induces opening of the transferrin N-lobe, which results in a negative effect on the stability. Presence of NaCl or arginine at low pH enhances the opening and has a negative impact on the overall protein stability. Statement of Significance Protein-based therapeutics have become an essential part of medical treatment. They are highly specific, have high affinity and fewer off-target effects. However, stabilization of proteins is critical, time-consuming, and expensive, and it is not yet possible to predict the behavior of proteins under different conditions. The current work is focused on a molecular understanding of the stability of human serum transferrin; a protein which is abundant in blood serum, may pass the blood brain barrier and therefore with high potential in drug delivery. Combination of high throughput unfolding techniques and structural studies, using small angle X-ray scattering and molecular dynamic simulations, allows us to understand the behavior of transferrin on a molecular level.
Collapse
|
16
|
The Effect of Point Mutations on the Biophysical Properties of an Antimicrobial Peptide: Development of a Screening Protocol for Peptide Stability Screening. Mol Pharm 2020; 17:3298-3313. [DOI: 10.1021/acs.molpharmaceut.0c00406] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
17
|
Albumin-neprilysin fusion protein: understanding stability using small angle X-ray scattering and molecular dynamic simulations. Sci Rep 2020; 10:10089. [PMID: 32572086 PMCID: PMC7308280 DOI: 10.1038/s41598-020-67002-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Accepted: 05/27/2020] [Indexed: 12/16/2022] Open
Abstract
Fusion technology is widely used in protein-drug development to increase activity, stability, and bioavailability of protein therapeutics. Fusion proteins, like any other type of biopharmaceuticals, need to remain stable during production and storage. Due to the high complexity and additional intramolecular interactions, it is not possible to predict the behavior of fusion proteins based on the behavior the individual proteins. Therefore, understanding the stability of fusion proteins on the molecular level is crucial for the development of biopharmaceuticals. The current study on the albumin-neprilysin (HSA-NEP) fusion protein uses a combination of thermal and chemical unfolding with small angle X-ray scattering and molecular dynamics simulations to show a correlation between decreasing stability and increasing repulsive interactions, which is unusual for most biopharmaceuticals. It is also seen that HSA-NEP is not fully flexible: it is present in both compact and extended conformations. Additionally, the volume fraction of each conformation changes with pH. Finally, the presence of NaCl and arginine increases stability at pH 6.5, but decreases stability at pH 5.0.
Collapse
|
18
|
Water-Intake and Water-Molecule Paths to the Active Site of Secretory Phospholipase A 2 Studied Using MD Simulations and the Tracking Tool AQUA-DUCT. J Phys Chem B 2020; 124:1881-1891. [PMID: 32064878 DOI: 10.1021/acs.jpcb.9b10837] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Secretory phospholipases A2 (sPLA2s) are a subclass of enzymes that catalyze the hydrolysis at the sn-2 position of glycerophospholipids, producing free fatty acids and lysophospholipids. In this study, different phospholipids with structural modifications close to the scissile sn-2 ester bond were studied to determine the effect of the structural changes on the formation of the Michaelis-Menten complex and the water entry/exit pathways using molecular dynamics simulations and the computational tracking tool AQUA-DUCT. Structural modifications include methylation, dehydrogenation, and polarization close to the sn-2 scissile bond. We found that all water molecules reaching the active site of sPLA2-IIA pass by the aromatic residues Phe5 and Tyr51 and enter the active site through an active-site cleft. The relative amount of water available for the enzymatic reaction of the different phospholipid-sPLA2 complexes was determined together with the distance between key atoms in the catalytic machinery. The results showed that (Z)-unsaturated phospholipid is a good substrate for sPLA2-IIA. The computational results are in good agreement with previously reported experimental data on the ability of sPLA2-IIA to hydrolyze liposomes made from the different phospholipids, and the results provide new insights into the necessary active-site solvation of the Michaelis-Menten complex and can pave the road for rational design in engineering applications.
Collapse
|
19
|
Water Distribution and Clustering on the Lyophilized IgG1 Surface: Insight from Molecular Dynamics Simulations. Mol Pharm 2020; 17:900-908. [PMID: 31990562 DOI: 10.1021/acs.molpharmaceut.9b01150] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Water has a critical role in the stability of the higher-order structure of proteins. In addition, it is considered to be a major destabilization factor for the physical and chemical stability of freeze-dried proteins and peptides. Physical and chemical aspects of protein/water relationships are commonly studied with the use of water vapor sorption isotherms for amorphous lyophilized proteins, which, in turn, are commonly analyzed using the Brunauer-Emmett-Teller (BET) equation to obtain the parameters, Wm and CB. The parameter Wm is generally referred to as the "monolayer limit of adsorption" and has a narrow range of 6-8% for most proteins. In this study, the water distribution on an IgG1 surface is investigated by molecular dynamics (MD) simulations at different water contents. The monolayer of water molecules was found to have limited coverage of the protein surface, and the true monolayer coverage of the protein globule actually occurs at a hydration level above 30%. The distribution of water molecules on the IgG1 surface is also highly heterogeneous, and the heterogeneity is not considered in the BET theory. In this study, a mechanistic model has been developed to describe the water vapor sorption isotherm. This model is based on the analysis of the hydrogen bonding network extracted from the MD simulations. The model is consistent with the experimental Type-II isotherm, which is usually observed for proteins. The physical meaning of the BET monolayer was redefined as the onset of water cluster formation. A simple model to calculate the onset water level, Wm, is proposed based on the hydration of different amino acids, as determined from the MD simulations.
Collapse
|
20
|
The Catalytic Acid–Base in GH109 Resides in a Conserved GGHGG Loop and Allows for Comparable α-Retaining and β-Inverting Activity in an N-Acetylgalactosaminidase from Akkermansia muciniphila. ACS Catal 2020. [DOI: 10.1021/acscatal.9b04474] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
21
|
Advancing Therapeutic Protein Discovery and Development through Comprehensive Computational and Biophysical Characterization. Mol Pharm 2020; 17:426-440. [PMID: 31790599 DOI: 10.1021/acs.molpharmaceut.9b00852] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.
Collapse
|
22
|
Removal of N-linked glycans in cellobiohydrolase Cel7A from Trichoderma reesei reveals higher activity and binding affinity on crystalline cellulose. BIOTECHNOLOGY FOR BIOFUELS 2020; 13:136. [PMID: 32782472 PMCID: PMC7412794 DOI: 10.1186/s13068-020-01779-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2020] [Accepted: 07/29/2020] [Indexed: 05/15/2023]
Abstract
BACKGROUND Cellobiohydrolase from glycoside hydrolase family 7 is a major component of commercial enzymatic mixtures for lignocellulosic biomass degradation. For many years, Trichoderma reesei Cel7A (TrCel7A) has served as a model to understand structure-function relationships of processive cellobiohydrolases. The architecture of TrCel7A includes an N-glycosylated catalytic domain, which is connected to a carbohydrate-binding module through a flexible, O-glycosylated linker. Depending on the fungal expression host, glycosylation can vary not only in glycoforms, but also in site occupancy, leading to a complex pattern of glycans, which can affect the enzyme's stability and kinetics. RESULTS Two expression hosts, Aspergillus oryzae and Trichoderma reesei, were utilized to successfully express wild-types TrCel7A (WT Ao and WT Tr ) and the triple N-glycosylation site deficient mutants TrCel7A N45Q, N270Q, N384Q (ΔN-glyc Ao and ΔN-glyc Tr ). Also, we expressed single N-glycosylation site deficient mutants TrCel7A (N45Q Ao , N270Q Ao , N384Q Ao ). The TrCel7A enzymes were studied by steady-state kinetics under both substrate- and enzyme-saturating conditions using different cellulosic substrates. The Michaelis constant (K M ) was consistently found to be lowered for the variants with reduced N-glycosylation content, and for the triple deficient mutants, it was less than half of the WTs' value on some substrates. The ability of the enzyme to combine productively with sites on the cellulose surface followed a similar pattern on all tested substrates. Thus, site density (number of sites per gram cellulose) was 30-60% higher for the single deficient variants compared to the WT, and about twofold larger for the triple deficient enzyme. Molecular dynamic simulation of the N-glycan mutants TrCel7A revealed higher number of contacts between CD and cellulose crystal upon removal of glycans at position N45 and N384. CONCLUSIONS The kinetic changes of TrCel7A imposed by removal of N-linked glycans reflected modifications of substrate accessibility. The presence of N-glycans with extended structures increased K M and decreased attack site density of TrCel7A likely due to steric hindrance effect and distance between the enzyme and the cellulose surface, preventing the enzyme from achieving optimal conformation. This knowledge could be applied to modify enzyme glycosylation to engineer enzyme with higher activity on the insoluble substrates.
Collapse
|
23
|
Investigations of Albumin–Insulin Detemir Complexes Using Molecular Dynamics Simulations and Free Energy Calculations. Mol Pharm 2019; 17:132-144. [DOI: 10.1021/acs.molpharmaceut.9b00839] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
24
|
Structural and functional aspects of mannuronic acid-specific PL6 alginate lyase from the human gut microbe Bacteroides cellulosilyticus. J Biol Chem 2019; 294:17915-17930. [PMID: 31530640 PMCID: PMC6879350 DOI: 10.1074/jbc.ra119.010206] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Revised: 09/16/2019] [Indexed: 01/28/2023] Open
Abstract
Alginate is a linear polysaccharide from brown algae consisting of 1,4-linked β-d-mannuronic acid (M) and α-l-guluronic acid (G) arranged in M, G, and mixed MG blocks. Alginate was assumed to be indigestible in humans, but bacteria isolated from fecal samples can utilize alginate. Moreover, genomes of some human gut microbiome-associated bacteria encode putative alginate-degrading enzymes. Here, we genome-mined a polysaccharide lyase family 6 alginate lyase from the gut bacterium Bacteroides cellulosilyticus (BcelPL6). The structure of recombinant BcelPL6 was solved by X-ray crystallography to 1.3 Å resolution, revealing a single-domain, monomeric parallel β-helix containing a 10-step asparagine ladder characteristic of alginate-converting parallel β-helix enzymes. Substitutions of the conserved catalytic site residues Lys-249, Arg-270, and His-271 resulted in activity loss. However, imidazole restored the activity of BcelPL6-H271N to 2.5% that of the native enzyme. Molecular docking oriented tetra-mannuronic acid for syn attack correlated with M specificity. Using biochemical analyses, we found that BcelPL6 initially releases unsaturated oligosaccharides of a degree of polymerization of 2-7 from alginate and polyM, which were further degraded to di- and trisaccharides. Unlike other PL6 members, BcelPL6 had low activity on polyMG and none on polyG. Surprisingly, polyG increased BcelPL6 activity on alginate 7-fold. LC-electrospray ionization-MS quantification of products and lack of activity on NaBH4-reduced octa-mannuronic acid indicated that BcelPL6 is an endolyase that further degrades the oligosaccharide products with an intact reducing end. We anticipate that our results advance predictions of the specificity and mode of action of PL6 enzymes.
Collapse
|
25
|
Solution structures of long-acting insulin analogues and their complexes with albumin. ACTA CRYSTALLOGRAPHICA SECTION D-STRUCTURAL BIOLOGY 2019; 75:272-282. [PMID: 30950398 DOI: 10.1107/s2059798318017552] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/08/2018] [Accepted: 12/11/2018] [Indexed: 12/15/2022]
Abstract
The lipidation of peptide drugs is one strategy to obtain extended half-lives, enabling once-daily or even less frequent injections for patients. The half-life extension results from a combination of self-association and association with human serum albumin (albumin). The self-association and association with albumin of two insulin analogues, insulin detemir and insulin degludec, were investigated by small-angle X-ray scattering (SAXS) and dynamic light scattering (DLS) in phenolic buffers. Detemir shows concentration-dependent self-association, with an equilibrium between hexamer, dihexamer, trihexamer and larger species, while degludec appears as a dihexamer independent of concentration. The solution structure of the detemir trihexamer has a bent shape. The stoichiometry of the association with albumin was studied using DLS. For albumin-detemir the molar stoichiometry was determined to be 1:6 (albumin:detemir ratio) and for albumin-degludec it was between 1:6 and 1:12 (albumin:degludec ratio). Batch SAXS measurements of a 1:6 albumin:detemir concentration series revealed a concentration dependence of complex formation. The data allowed the modelling of a complex between albumin and a detemir hexamer and a complex consisting of two albumins binding to opposite ends of a detemir dihexamer. Measurements of size-exclusion chromatography coupled to SAXS revealed a complex between a degludec dihexamer and albumin. Based on the results, equilibria for the albumin-detemir and albumin-degludec mixtures are proposed.
Collapse
|
26
|
Self-Interaction of Human Serum Albumin: A Formulation Perspective. ACS OMEGA 2018; 3:16105-16117. [PMID: 30556026 PMCID: PMC6288999 DOI: 10.1021/acsomega.8b02245] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2018] [Accepted: 11/13/2018] [Indexed: 05/12/2023]
Abstract
In the present study, small-angle X-ray scattering (SAXS) and static light scattering (SLS) have been used to study the solution properties and self-interaction of recombinant human serum albumin (rHSA) molecules in three pharmaceutically relevant buffer systems. Measurements are carried out up to high protein concentrations and as a function of ionic strength by adding sodium chloride to probe the role of electrostatic interactions. The effective structure factors (S eff) as a function of the scattering vector magnitude q have been extracted from the scattering profiles and fit to the solution of the Ornstein-Zernike equation using a screened Yukawa potential to describe the double-layer force. Although only a limited q range is used, accurate fits required including an electrostatic repulsion element in the model at low ionic strength, while only a hard sphere model with a tunable diameter is necessary for fitting to high-ionic-strength data. The fit values of net charge agree with available data from potentiometric titrations. Osmotic compressibility data obtained by extrapolating the SAXS profiles or directly from SLS measurements has been fit to a 10-term virial expansion for hard spheres and an equation of state for hard biaxial ellipsoids. We show that modeling rHSA as an ellipsoid, rather than a sphere, provides a much more accurate fit for the thermodynamic data over the entire concentration range. Osmotic virial coefficient data, derived at low protein concentration, can be used to parameterize the model for predicting the behavior up to concentrations as high as 450 g/L. The findings are especially important for the biopharmaceutical sector, which require approaches for predicting concentrated protein solution behavior using minimal sample consumption.
Collapse
|
27
|
Dual Nicotinic Acetylcholine Receptor α4β2 Antagonists/α7 Agonists: Synthesis, Docking Studies, and Pharmacological Evaluation of Tetrahydroisoquinolines and Tetrahydroisoquinolinium Salts. J Med Chem 2018; 61:1719-1729. [DOI: 10.1021/acs.jmedchem.7b01895] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
28
|
Interaction between structurally different heteroexopolysaccharides and β-lactoglobulin studied by solution scattering and analytical ultracentrifugation. Int J Biol Macromol 2018; 111:746-754. [PMID: 29329814 DOI: 10.1016/j.ijbiomac.2018.01.050] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2017] [Revised: 12/24/2017] [Accepted: 01/08/2018] [Indexed: 12/28/2022]
Abstract
Despite a very large number of bacterial exopolysaccharides have been reported, detailed knowledge on their molecular structures and associative interactions with proteins is lacking. Small-angle X-ray scattering, dynamic light scattering and analytical ultracentrifugation (AUC) were used to characterize the interactions of six lactic acid bacterial heteroexopolysaccharides (HePS-1-HePS-6) with β-lactoglobulin (BLG). Compared to free HePSs, a large increase in the X-ray radius of gyration RG, maximum length L and hydrodynamic diameter dH of HePS-1-HePS-4 mixed with BLG revealed strong aggregation, the extent of which depended on the compact conformation and degree of branching of these HePSs. No significant effects were observed with HePS-5 and HePS-6. Turbidity and AUC analyses showed that both soluble and insoluble BLG-HePS complexes were formed. The findings provide new insights into the role of molecular structures in associative interactions between HePSs and BLG which has relevance for various industrial applications.
Collapse
|
29
|
The influence of different linker modifications on the catalytic activity and cellulose affinity of cellobiohydrolase Cel7A from Hypocrea jecorina. Protein Eng Des Sel 2017; 30:495-501. [PMID: 28873985 DOI: 10.1093/protein/gzx036] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2017] [Accepted: 07/04/2017] [Indexed: 12/13/2022] Open
Abstract
Various cellulases consist of a catalytic domain connected to a carbohydrate-binding module (CBM) by a flexible linker peptide. The linker if often strongly O-glycosylated and typically has a length of 20-50 amino acid residues. Functional roles, other than connecting the two folded domains, of the linker and its glycans, have been widely discussed, but experimental evidence remains sparse. One of the most studied cellulose degrading enzymes is the multi-domain cellobiohydrolase Cel7A from Hypocrea jecorina. Here, we designed variants of Cel7A with mutations in the linker region to elucidate the role of the linker. We found that moderate modification of the linker could result in significant changes in substrate affinity and catalytic efficacy. These changes were quite different for different linker variants. Thus, deletion of six residues near the catalytic domain had essentially no effects on enzyme function. Conversely, a substitution of four glycosylation sites near the middle of the linker reduced substrate affinity and increased maximal turnover. The observation of weaker binding provides some support of recent suggestions that linker glycans may be directly involved in substrate interactions. However, a variant with several inserted glycosylation sites near the CBM also showed lower affinity for the substrate compared to the wild-type, and we suggest that substrate interactions of the glycans depend on their exact location as well as other factors such as changes in structure and dynamics of the linker peptide.
Collapse
|
30
|
Small-Angle X-ray Scattering Data in Combination with RosettaDock Improves the Docking Energy Landscape. J Chem Inf Model 2017; 57:2463-2475. [DOI: 10.1021/acs.jcim.6b00789] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
31
|
Glucagon-like Peptide 1 Conjugated to Recombinant Human Serum Albumin Variants with Modified Neonatal Fc Receptor Binding Properties. Impact on Molecular Structure and Half-Life. Biochemistry 2017; 56:4860-4870. [PMID: 28799326 DOI: 10.1021/acs.biochem.7b00492] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Glucagon-like peptide 1 (GLP-1) is a small incretin hormone stimulated by food intake, resulting in an amplification of the insulin response. Though GLP-1 is interesting as a drug candidate for the treatment of type 2 diabetes mellitus, its short plasma half-life of <3 min limits its clinical use. A strategy for extending the half-life of GLP-1 utilizes the long half-life of human serum albumin (HSA) by combining the two via chemical conjugation or genetic fusion. HSA has a plasma half-life of around 21 days because of its interaction with the neonatal Fc receptor (FcRn) expressed in endothelial cells of blood vessels, which rescues circulating HSA from lysosomal degradation. We have conjugated GLP-1 to C34 of native sequence recombinant HSA (rHSA) and two rHSA variants, one with increased and one with decreased binding affinity for human FcRn. We have investigated the impact of conjugation on FcRn binding affinities, GLP-1 potency, and pharmacokinetics, combined with the solution structure of the rHSA variants and GLP-1-albumin conjugates. The solution structures, determined by small-angle X-ray scattering, show the GLP-1 pointing away from the surface of rHSA. Combining the solution structures with the available structural information about the FcRn and GLP-1 receptor obtained from X-ray crystallography, we can explain the observed in vitro and in vivo behavior. We conclude that the conjugation of GLP-1 to rHSA does not affect the interaction between rHSA and FcRn, while the observed decrease in the potency of GLP-1 can be explained by a steric hindrance of binding of GLP-1 to its receptor.
Collapse
|
32
|
Revealing the Compact Structure of Lactic Acid Bacterial Heteroexopolysaccharides by SAXS and DLS. Biomacromolecules 2017; 18:747-756. [DOI: 10.1021/acs.biomac.6b01597] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|