1
|
Jia Y, Fernandez A, Sampath J. PEGylation of Insulin and Lysozyme To Stabilize against Thermal Denaturation: A Molecular Dynamics Simulation Study. J Phys Chem B 2023; 127:6856-6866. [PMID: 37498538 DOI: 10.1021/acs.jpcb.3c01289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Biologic drugs or "biologics" (proteins derived from living organisms) are one of the fastest-growing classes of FDA-approved therapeutics. These compounds are often fragile and require conjugation to polymers for stabilization, with many proteins too ephemeral for therapeutic use. During storage or administration, proteins tend to unravel and lose their secondary structure due to changes in solution temperature, pH, and other external stressors. To enhance their lifetime, protein drugs currently in the market are conjugated with polyethylene glycol (PEG), owing to its ability to increase the stability, solubility, and pharmacokinetics of protein drugs. Here, we perform all-atom molecular dynamics simulations to study the unfolding process of egg-white lysozyme and insulin at elevated temperatures. We test the validity of two force fields─CHARMM36 and Amber ff99SB-ILDN─in the unfolding process. By calculating global and local properties, we capture residues that deteriorate first─these are the "weak links" in the proteins. Next, we conjugate both proteins with PEG and find that PEG preserves the native structure of the proteins at elevated temperatures by blocking water molecules from entering the hydrophobic core, thereby causing the secondary structure to stabilize.
Collapse
Affiliation(s)
- Yinhao Jia
- Department of Chemical Engineering, University of Florida, Gainesville, Florida 32611, United States
| | - Adam Fernandez
- Department of Chemical Engineering, University of Florida, Gainesville, Florida 32611, United States
| | - Janani Sampath
- Department of Chemical Engineering, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
2
|
Wang L, O'Mara ML. Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein. J Chem Theory Comput 2021; 17:6491-6508. [PMID: 34506133 DOI: 10.1021/acs.jctc.1c00414] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.
Collapse
Affiliation(s)
- Lily Wang
- Research School of Chemistry, College of Science, Australian National University, Canberra, ACT 2601, Australia
| | - Megan L O'Mara
- Research School of Chemistry, College of Science, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
3
|
Man VH, He X, Derreumaux P, Ji B, Xie XQ, Nguyen PH, Wang J. Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ 16-22 Dimer. J Chem Theory Comput 2019; 15:1440-1452. [PMID: 30633867 PMCID: PMC6745714 DOI: 10.1021/acs.jctc.8b01107] [Citation(s) in RCA: 87] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ16-22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16-22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique UPR 9080, CNRS, Université Denis Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Beihong Ji
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xiang-Qun Xie
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Phuong H. Nguyen
- Laboratoire de Biochimie Théorique UPR 9080, CNRS, Université Denis Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Junmei Wang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
- Corresponding Author:
| |
Collapse
|
4
|
Thompson HN, Thompson CE, Andrade Caceres R, Dardenne LE, Netz PA, Stassen H. Prion protein conversion triggered by acidic condition: a molecular dynamics study through different force fields. J Comput Chem 2018; 39:2000-2011. [DOI: 10.1002/jcc.25380] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2018] [Revised: 05/15/2018] [Accepted: 05/26/2018] [Indexed: 11/06/2022]
Affiliation(s)
- Helen Nathalia Thompson
- Departamento de Físico-Química, Instituto de Química; Universidade Federal do Rio Grande do Sul; 91501-970 Porto Alegre Rio Grande do Sul Brazil
| | - Claudia Elizabeth Thompson
- Departamento de Farmacociências; Universidade Federal de Ciências da Saúde de Porto Alegre; 90050-170 Porto Alegre Rio Grande do Sul Brazil
| | - Rafael Andrade Caceres
- Departamento de Farmacociências; Universidade Federal de Ciências da Saúde de Porto Alegre; 90050-170 Porto Alegre Rio Grande do Sul Brazil
| | | | - Paulo Augusto Netz
- Departamento de Físico-Química, Instituto de Química; Universidade Federal do Rio Grande do Sul; 91501-970 Porto Alegre Rio Grande do Sul Brazil
| | - Hubert Stassen
- Departamento de Físico-Química, Instituto de Química; Universidade Federal do Rio Grande do Sul; 91501-970 Porto Alegre Rio Grande do Sul Brazil
| |
Collapse
|
5
|
Smith MD, Rao JS, Segelken E, Cruz L. Force-Field Induced Bias in the Structure of Aβ21-30: A Comparison of OPLS, AMBER, CHARMM, and GROMOS Force Fields. J Chem Inf Model 2015; 55:2587-95. [PMID: 26629886 DOI: 10.1021/acs.jcim.5b00308] [Citation(s) in RCA: 73] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In this work we examine the dynamics of an intrinsically disordered protein fragment of the amyloid β, the Aβ21-30, under seven commonly used molecular dynamics force fields (OPLS-AA, CHARMM27-CMAP, AMBER99, AMBER99SB, AMBER99SB-ILDN, AMBER03, and GROMOS53A6), and three water models (TIP3P, TIP4P, and SPC/E). We find that the tested force fields and water models have little effect on the measures of radii of gyration and solvent accessible surface area (SASA); however, secondary structure measures and intrapeptide hydrogen-bonding are significantly modified, with AMBER (99, 99SB, 99SB-ILDN, and 03) and CHARMM22/27 force-fields readily increasing helical content and the variety of intrapeptide hydrogen bonds. On the basis of a comparison between the population of helical and β structures found in experiments, our data suggest that force fields that suppress the formation of helical structure might be a better choice to model the Aβ21-30 peptide.
Collapse
Affiliation(s)
- Micholas Dean Smith
- Department of Physics, Drexel University , 3141 Chestnut Street, Philadelphia, Pennsylvania 19104, United States
| | - J Srinivasa Rao
- Department of Physics, Drexel University , 3141 Chestnut Street, Philadelphia, Pennsylvania 19104, United States.,Department of Physics, New Jersey Institute of Technology , University Heights, Newark, New Jersey 07102-1982, United States
| | - Elizabeth Segelken
- Department of Physics, Drexel University , 3141 Chestnut Street, Philadelphia, Pennsylvania 19104, United States
| | - Luis Cruz
- Department of Physics, Drexel University , 3141 Chestnut Street, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
6
|
Hu JP, He HQ, Jiao X, Chang S. Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.773431] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
7
|
Nguyen PH, Li MS, Derreumaux P. Effects of all-atom force fields on amyloid oligomerization: replica exchange molecular dynamics simulations of the Aβ16–22 dimer and trimer. Phys Chem Chem Phys 2011; 13:9778-88. [DOI: 10.1039/c1cp20323a] [Citation(s) in RCA: 148] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|