1
|
Piskorz T, Perez-Chirinos L, Qiao B, Sasselli IR. Tips and Tricks in the Modeling of Supramolecular Peptide Assemblies. ACS OMEGA 2024; 9:31254-31273. [PMID: 39072142 PMCID: PMC11270692 DOI: 10.1021/acsomega.4c02628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Revised: 06/17/2024] [Accepted: 06/19/2024] [Indexed: 07/30/2024]
Abstract
Supramolecular peptide assemblies (SPAs) hold promise as materials for nanotechnology and biomedicine. Although their investigation often entails adapting experimental techniques from their protein counterparts, SPAs are fundamentally distinct from proteins, posing unique challenges for their study. Computational methods have emerged as indispensable tools for gaining deeper insights into SPA structures at the molecular level, surpassing the limitations of experimental techniques, and as screening tools to reduce the experimental search space. However, computational studies have grappled with issues stemming from the absence of standardized procedures and relevant crystal structures. Fundamental disparities between SPAs and protein simulations, such as the absence of experimentally validated initial structures and the importance of the simulation size, number of molecules, and concentration, have compounded these challenges. Understanding the roles of various parameters and the capabilities of different models and simulation setups remains an ongoing endeavor. In this review, we aim to provide readers with guidance on the parameters to consider when conducting SPA simulations, elucidating their potential impact on outcomes and validity.
Collapse
Affiliation(s)
| | - Laura Perez-Chirinos
- Center
for Cooperative Research in Biomaterials (CIC biomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramón 182, 20014 Donostia-San Sebastián, Spain
| | - Baofu Qiao
- Department
of Natural Sciences, Baruch College, City
University of New York, New York, New York 10010, United States
| | - Ivan R. Sasselli
- Centro
de Física de Materiales (CFM), CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 San Sebastián, Spain
| |
Collapse
|
2
|
Kienlein M, Zacharias M. How arginine inhibits substrate-binding domain 2 elucidated using molecular dynamics simulations. Protein Sci 2024; 33:e5077. [PMID: 38888275 PMCID: PMC11184577 DOI: 10.1002/pro.5077] [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: 02/08/2024] [Revised: 04/19/2024] [Accepted: 05/28/2024] [Indexed: 06/20/2024]
Abstract
The substrate-binding domain 2 (SBD2) is an important part of the bacterial glutamine (GLN) transporter and mediates binding and delivery of GLN to the transporter translocation subunit. The SBD2 consists of two domains, D1 and D2, that bind GLN in the space between domains in a closed structure. In the absence of ligand, the SBD2 adopts an open conformation with larger space between domains. The GLN binding and closing are essential for the subsequent transport into the cell. Arginine (ARG) can also bind to SBD2 but does not induce closing and inhibits GLN transport. We use atomistic molecular dynamics (MD) simulations in explicit solvent to study ARG binding in the presence of the open SBD2 structure and observed reversible binding to the native GLN binding site with similar contacts but no transition to a closed SBD2 state. Absolute binding free energy simulations predict a considerable binding affinity of ARG and GLN to the binding site on the D1 domain. Free energy simulations to induce subsequent closing revealed a strong free energy penalty in case of ARG binding in contrast to GLN binding that favors the closed SBD2 state but still retains a free energy barrier for closing. The simulations allowed the identification of the molecular origin of the closing penalty in case of bound ARG and suggested a mutation of lysine at position 373 to alanine that strongly reduced the penalty and allowed closing even in the presence of bound ARG. The study offers an explanation of the molecular mechanism of how ARG competitively inhibits GLN from binding to SBD2 and from triggering the transition to a closed conformation. The proposed Lys373Ala mutation shows promise as a potential tool to validate whether a conformational mismatch between open SBD2 and the translocator is responsible for preventing ARG uptake to the cell.
Collapse
Affiliation(s)
- Maximilian Kienlein
- Center for Functional Protein Assemblies (CPA)Technical University of MunichGarchingGermany
| | - Martin Zacharias
- Center for Functional Protein Assemblies (CPA)Technical University of MunichGarchingGermany
| |
Collapse
|
3
|
Andrews B, Schweitzer-Stenner R, Urbanc B. Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data. J Phys Chem B 2024; 128:6217-6231. [PMID: 38877893 PMCID: PMC11215781 DOI: 10.1021/acs.jpcb.4c02278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/29/2024] [Accepted: 05/30/2024] [Indexed: 06/28/2024]
Abstract
Molecular dynamics (MD) is a great tool for elucidating conformational dynamics of proteins and peptides in water at the atomistic level that often surpasses the level of detail available experimentally. Structure predictions, however, are limited by the accuracy of the underlying MD force field. This limitation is particularly stark in the case of intrinsically disordered peptides and proteins, which are characterized by solvent-accessible and disordered peptide regions and domains. Recent studies show that most additive MD force fields, including CHARMM36m, do not reproduce the intrinsic conformational distributions of guest amino acid residues x in cationic GxG peptides in water in line with experimental data. Positing that a lack of polarizability in additive MD force fields may be the culprit for the reported discrepancies, we here examine the conformational dynamics of guest glycine and alanine residues in cationic GxG peptides in water using two polarizable MD force fields, CHARMM Drude and AMOEBA. Our results indicate that while AMOEBA captures the experimental data better than CHARMM Drude, neither of the two polarizable force fields offers an improvement of the Ramachandran distributions of glycine and alanine residues in cationic GGG and GAG peptides, respectively, over CHARMM36m.
Collapse
Affiliation(s)
- Brian Andrews
- Department
of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | | | - Brigita Urbanc
- Department
of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
4
|
Emperador A. PACSAB Server: A Web-Based Tool for the Study of Aggregation and the Conformational Ensemble of Disordered and Folded Proteins. Int J Mol Sci 2024; 25:6021. [PMID: 38892222 PMCID: PMC11172606 DOI: 10.3390/ijms25116021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Revised: 05/20/2024] [Accepted: 05/24/2024] [Indexed: 06/21/2024] Open
Abstract
We present in this article the PACSAB server, which is designed to provide information about the structural ensemble and interactions of both stable and disordered proteins to researchers in the field of molecular biology. The use of this tool does not require any computational skills as the user just needs to upload the structure of the protein to be studied; the server runs a simulation with the PACSAB model, a highly accurate coarse-grained model that is much more efficient than standard molecular dynamics for the exploration of the conformational space of multiprotein systems. The trajectories generated by the simulations based on this model reveal the propensity of the protein under study for aggregation, identify the residues playing a central role in the aggregation process, and reproduce the whole conformational space of disordered proteins. All of this information is shown and can be downloaded from the web page.
Collapse
Affiliation(s)
- Agustí Emperador
- Department of Physics, Universitat Politècnica de Catalunya, B4-B5 Campus Nord, Jordi Girona 1-3, 08034 Barcelona, Spain
| |
Collapse
|
5
|
Grassmann G, Miotto M, Desantis F, Di Rienzo L, Tartaglia GG, Pastore A, Ruocco G, Monti M, Milanetti E. Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments. Chem Rev 2024; 124:3932-3977. [PMID: 38535831 PMCID: PMC11009965 DOI: 10.1021/acs.chemrev.3c00550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 02/20/2024] [Accepted: 02/21/2024] [Indexed: 04/11/2024]
Abstract
Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.
Collapse
Affiliation(s)
- Greta Grassmann
- Department
of Biochemical Sciences “Alessandro Rossi Fanelli”, Sapienza University of Rome, Rome 00185, Italy
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Mattia Miotto
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Fausta Desantis
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- The
Open University Affiliated Research Centre at Istituto Italiano di
Tecnologia, Genoa 16163, Italy
| | - Lorenzo Di Rienzo
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Gian Gaetano Tartaglia
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
- Center
for Human Technologies, Genoa 16152, Italy
| | - Annalisa Pastore
- Experiment
Division, European Synchrotron Radiation
Facility, Grenoble 38043, France
| | - Giancarlo Ruocco
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| | - Michele Monti
- RNA
System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
| | - Edoardo Milanetti
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| |
Collapse
|
6
|
Sha H, Zhu F. Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model. J Phys Chem B 2024; 128:960-972. [PMID: 38251836 DOI: 10.1021/acs.jpcb.3c06881] [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: 01/23/2024]
Abstract
HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.
Collapse
Affiliation(s)
- Hao Sha
- Department of Physics, Indiana University─Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Fangqiang Zhu
- Department of Physics, Indiana University─Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
7
|
Pitaloka DAE, Arfan A, Ramadhan DSF, Chaidir L. Insights from the molecular mechanism of pyrazinamide to mutated pyrazinamidase linked to the pncA gene in clinical isolates of Mycobacterium tuberculosis. J Biomol Struct Dyn 2024; 42:759-765. [PMID: 37096659 DOI: 10.1080/07391102.2023.2195002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 03/18/2023] [Indexed: 04/26/2023]
Abstract
This study aims to conduct a comprehensive molecular dynamics strategy to evaluate whether mutations found in pyrazinamide monoresistant (PZAMR) strains of Mycobacterium tuberculosis (MTB) can potentially reduce the effectiveness of pyrazinamide (PZA) for tuberculosis (TB) treatment. Five single point mutations of pyrazinamidase (PZAse), an enzyme which is responsible for the activation of prodrug PZA into pyrazinoic acid, found in MTB clinical isolates, namely His82Arg, Thr87Met, Ser66Pro, Ala171Val, and Pro62Leu, were analyzed by the dynamics simulations both in the apo state (unbound state) and in the PZA bound state. The results showed that the mutation of His82 to Arg, Thr87 to Met, and Ser66 to Pro in PZAse affects the coordination state of the Fe2+ ion, which is a cofactor required for enzyme activity. These mutations change the flexibility, stability, and fluctuation of His51, His57, and ASP49 amino acid residues around the Fe2+ ion, culminating in an unstable complex and dissociation of PZA from the PZAse binding site. However, mutations of Ala171 to Val and Pro62 to Leu were found to have no effect on the complex's stability. Based on the results, PZAse mutations of His82Arg, Thr87Met, and Ser66Pro culminated in weak binding affinity for PZA and caused significant structural deformations that led to PZA resistance. Future structural and functional studies, as well as investigations into other aspects of drug resistance in PZAse, will require experimental clarification.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Dian Ayu Eka Pitaloka
- Department of Pharmacology and Clinical Pharmacy, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
- Center for Translational Biomarker Research, Universitas Padjadjaran, Bandung, Indonesia
| | - Arfan Arfan
- Department of Medicinal Chemistry, Faculty of Pharmacy, Universitas Halu Oleo, Kendari, Indonesia
| | - Dwi Syah Fitra Ramadhan
- Department of Pharmaceutical Chemistry, Sekolah Tinggi Ilmu Kesehatan Mandala Waluya, Kendari, Indonesia
| | - Lidya Chaidir
- Center for Translational Biomarker Research, Universitas Padjadjaran, Bandung, Indonesia
- Department of Biomedical Sciences, Faculty of Medicine, Universitas Padjadjaran, Sumedang, Indonesia
| |
Collapse
|
8
|
Linse JB, Hub JS. Scrutinizing the protein hydration shell from molecular dynamics simulations against consensus small-angle scattering data. Commun Chem 2023; 6:272. [PMID: 38086909 PMCID: PMC10716392 DOI: 10.1038/s42004-023-01067-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 11/20/2023] [Indexed: 06/09/2024] Open
Abstract
Biological macromolecules in solution are surrounded by a hydration shell, whose structure differs from the structure of bulk solvent. While the importance of the hydration shell for numerous biological functions is widely acknowledged, it remains unknown how the hydration shell is regulated by macromolecular shape and surface composition, mainly because a quantitative probe of the hydration shell structure has been missing. We show that small-angle scattering in solution using X-rays (SAXS) or neutrons (SANS) provide a protein-specific probe of the protein hydration shell that enables quantitative comparison with molecular simulations. Using explicit-solvent SAXS/SANS predictions, we derived the effect of the hydration shell on the radii of gyration Rg of five proteins using 18 combinations of protein force field and water model. By comparing computed Rg values from SAXS relative to SANS in D2O with consensus SAXS/SANS data from a recent worldwide community effort, we found that several but not all force fields yield a hydration shell contrast in remarkable agreement with experiments. The hydration shell contrast captured by Rg values depends strongly on protein charge and geometric shape, thus providing a protein-specific footprint of protein-water interactions and a novel observable for scrutinizing atomistic hydration shell models against experimental data.
Collapse
Affiliation(s)
- Johanna-Barbara Linse
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, 66123, Germany
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, 66123, Germany.
| |
Collapse
|
9
|
Ojha AA, Votapka LW, Amaro RE. QMrebind: incorporating quantum mechanical force field reparameterization at the ligand binding site for improved drug-target kinetics through milestoning simulations. Chem Sci 2023; 14:13159-13175. [PMID: 38023523 PMCID: PMC10664576 DOI: 10.1039/d3sc04195f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 10/22/2023] [Indexed: 12/01/2023] Open
Abstract
Understanding the interaction of ligands with biomolecules is an integral component of drug discovery and development. Challenges for computing thermodynamic and kinetic quantities for pharmaceutically relevant receptor-ligand complexes include the size and flexibility of the ligands, large-scale conformational rearrangements of the receptor, accurate force field parameters, simulation efficiency, and sufficient sampling associated with rare events. Our recently developed multiscale milestoning simulation approach, SEEKR2 (Simulation Enabled Estimation of Kinetic Rates v.2), has demonstrated success in predicting unbinding (koff) kinetics by employing molecular dynamics (MD) simulations in regions closer to the binding site. The MD region is further subdivided into smaller Voronoi tessellations to improve the simulation efficiency and parallelization. To date, all MD simulations are run using general molecular mechanics (MM) force fields. The accuracy of calculations can be further improved by incorporating quantum mechanical (QM) methods into generating system-specific force fields through reparameterizing ligand partial charges in the bound state. The force field reparameterization process modifies the potential energy landscape of the bimolecular complex, enabling a more accurate representation of the intermolecular interactions and polarization effects at the bound state. We present QMrebind (Quantum Mechanical force field reparameterization at the receptor-ligand binding site), an ORCA-based software that facilitates reparameterizing the potential energy function within the phase space representing the bound state in a receptor-ligand complex. With SEEKR2 koff estimates and experimentally determined kinetic rates, we compare and interpret the receptor-ligand unbinding kinetics obtained using the newly reparameterized force fields for model host-guest systems and HSP90-inhibitor complexes. This method provides an opportunity to achieve higher accuracy in predicting receptor-ligand koff rate constants.
Collapse
Affiliation(s)
- Anupam Anand Ojha
- Department of Chemistry and Biochemistry, University of California San Diego La Jolla California 92093 USA
| | - Lane William Votapka
- Department of Chemistry and Biochemistry, University of California San Diego La Jolla California 92093 USA
| | - Rommie Elizabeth Amaro
- Department of Molecular Biology, University of California San Diego La Jolla California 92093 USA
| |
Collapse
|
10
|
Prass T, Garidel P, Blech M, Schäfer LV. Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations. J Chem Inf Model 2023; 63:6129-6140. [PMID: 37757589 PMCID: PMC10565822 DOI: 10.1021/acs.jcim.3c00947] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Indexed: 09/29/2023]
Abstract
The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.
Collapse
Affiliation(s)
- Tobias
M. Prass
- Center
for Theoretical Chemistry, Ruhr University
Bochum, D-44780 Bochum, Germany
| | - Patrick Garidel
- Boehringer
Ingelheim Pharma GmbH & Co. KG, Innovation Unit, PDB, D-88397 Biberach
an der Riss, Germany
| | - Michaela Blech
- Boehringer
Ingelheim Pharma GmbH & Co. KG, Innovation Unit, PDB, D-88397 Biberach
an der Riss, Germany
| | - Lars V. Schäfer
- Center
for Theoretical Chemistry, Ruhr University
Bochum, D-44780 Bochum, Germany
| |
Collapse
|
11
|
Kasmawati H, Ruslin R, Arfan A, Sida NA, Saputra DI, Halimah E, Mustarichie R. Antibacterial Potency of an Active Compound from Sansevieria trifasciata Prain: An Integrated In Vitro and In Silico Study. Molecules 2023; 28:6096. [PMID: 37630348 PMCID: PMC10457997 DOI: 10.3390/molecules28166096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 08/02/2023] [Accepted: 08/04/2023] [Indexed: 08/27/2023] Open
Abstract
Sansevieria trifasciata Prain holds great potential as a valuable asset in pharmaceutical development. In this study, our focus is to explore and assess the antibacterial activity of various components derived from this plant, including extracts, fractions, subfractions, and isolates, explicitly targeting two common bacteria: Escherichia coli and Streptococcus aureus. The isolated compound, identified as a derivative pyridone alkaloid (5-methyl-11-(2-oxopyridin-1(2H)-yl)undecaneperoxoicacid), demonstrates notable antibacterial effects. The extracts, fractions, subfractions, and isolates reveal significant bacterial growth reductions (p < 0.05). The minimum inhibitory concentration (MIC) values for Escherichia coli were 1.95 ppm, 3.9 ppm, 15.62 ppm, and 7.81 ppm, respectively, while the MIC values for Streptococcus aureus were 1.95 ppm, 1.95 ppm, 15.62 ppm, and 7.81 ppm, respectively. Computational analysis showed the isolates' interaction with key residues on the active site of β-ketoacyl-ACP synthase from Escherichia coli and TyrRS from Streptococcus aureus. The findings indicate that the isolates exhibit a strong affinity for specific residues, including His333, Cys163, and Phe392 in β-ketoacyl-ACP synthase, as well as Arg88, His117, Glu160, and Gln213 in TyrRS. Comparative energy calculations using MMPBSA demonstrate the isolates' favorable binding energy (-104,101 kJ/mol for β-ketoacyl-ACP synthase and -81,060 kJ/mol for TyrRS) compared to ciprofloxacin. The elucidated antibacterial activity and molecular interactions of the isolates present valuable knowledge for future in vitro studies, facilitating the development of novel antibacterial agents targeting diverse bacterial strains.
Collapse
Affiliation(s)
- Henny Kasmawati
- Department of Pharmacy, Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.); (D.I.S.)
| | - Ruslin Ruslin
- Department of Pharmacy, Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.); (D.I.S.)
| | - Arfan Arfan
- Department of Pharmacy, Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.); (D.I.S.)
| | - Nurramadhani A. Sida
- Department of Pharmacy, Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.); (D.I.S.)
| | - Dimas Isnu Saputra
- Department of Pharmacy, Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.); (D.I.S.)
| | - Eli Halimah
- Department of Pharmacology and Clinical Pharmacy, Faculty of Pharmacy, Universitas Padjadjaran, Bandung 45363, Indonesia;
| | - Resmi Mustarichie
- Department of Analytical Pharmacy and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Bandung 45363, Indonesia
| |
Collapse
|
12
|
Saurabh S, Li Z, Hollowell P, Waigh T, Li P, Webster J, Seddon JM, Kalonia C, Lu JR, Bresme F. Structure and interaction of therapeutic proteins in solution: a combined simulation and experimental study. Mol Phys 2023; 121:e2236248. [PMID: 38107421 PMCID: PMC10721229 DOI: 10.1080/00268976.2023.2236248] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Accepted: 06/30/2023] [Indexed: 12/19/2023]
Abstract
The aggregation of therapeutic proteins in solution has attracted significant interest, driving efforts to understand the relationship between microscopic structural changes and protein-protein interactions determining aggregation processes in solution. Additionally, there is substantial interest in being able to predict aggregation based on protein structure as part of molecular developability assessments. Molecular Dynamics provides theoretical tools to complement experimental studies and to interrogate and identify the microscopic mechanisms determining aggregation. Here we perform all-atom MD simulations to study the structure and inter-protein interaction of the Fab and Fc fragments of the monoclonal antibody (mAb) COE3. We unravel the role of ion-protein interactions in building the ionic double layer and determining effective inter-protein interaction. Further, we demonstrate, using various state-of-the-art force fields (charmm, gromos, amber, opls/aa), that the protein solvation, ionic structure and protein-protein interaction depend significantly on the force field parameters. We perform SANS and Static Light Scattering experiments to assess the accuracy of the different forcefields. Comparison of the simulated and experimental results reveal significant differences in the forcefields' performance, particularly in their ability to predict the protein size in solution and inter-protein interactions quantified through the second virial coefficients. In addition, the performance of the forcefields is correlated with the protein hydration structure.
Collapse
Affiliation(s)
- Suman Saurabh
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London, United Kingdom
| | - Zongyi Li
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Peter Hollowell
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Thomas Waigh
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
- Photon Science Institute, The University of Manchester, Manchester, UK
| | - Peixun Li
- STFC ISIS Facility, Rutherford Appleton Laboratory, Didcot, UK
| | - John Webster
- STFC ISIS Facility, Rutherford Appleton Laboratory, Didcot, UK
| | - John M. Seddon
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London, United Kingdom
| | - Cavan Kalonia
- Dosage Form Design and Development, BioPharmaceutical Development, BioPharmaceuticals R&D, AstraZeneca, Gaithersburg, MD, USA
| | - Jian R. Lu
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London, United Kingdom
| |
Collapse
|
13
|
Panigrahy S, Sahu R, Reddy SK, Nayar D. Structure, energetics and dynamics in crowded amino acid solutions: a molecular dynamics study. Phys Chem Chem Phys 2023; 25:5430-5442. [PMID: 36744506 DOI: 10.1039/d2cp04238j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
A comprehensive understanding of crowding effects on biomolecular processes necessitates investigating the bulk thermodynamic and kinetic properties of the solutions with an accurate molecular representation of the crowded milieu. Recent studies have reparameterized the non-bonded dispersion interaction of solutes to precisely model intermolecular interactions, which would circumvent artificial aggregation as shown by the original force-fields. However, the performance of this reparameterization is yet to be assessed for concentrated crowded solutions in terms of investigating the hydration shell structure, energetics and dynamics. In this study, we perform molecular dynamics simulations of crowded aqueous solutions of five zwitterionic neutral amino acids (Gly, Ala, Thr, Pro, and Ser), mimicking the molecular crowding environment, using a modified AMBER ff99SB-ILDN force-field. We systematically examine and show that the reproducibility of the osmotic coefficients, density, viscosity and self-diffusivity of amino acids improves using the modified force-field in crowded concentrations. The modified force-field also improves the structuring of the solute solvation shells, solute interaction energy and convergence of tails of radial distribution functions, indicating reduction in the artificial aggregation. Our results also indicate that the hydrogen bonding network of water weakens and water molecules anomalously diffuse at small time scales in the crowded solutions. These results underscore the significance of examining the solution properties and anomalous hydration behaviour of water in crowded solutions, which have implications in shaping the structure and dynamics of biomolecules. The findings also illustrate the improvement in predicting bulk solution properties using the modified force-field, thereby providing an approach towards accurate modeling of crowded molecular solutions.
Collapse
Affiliation(s)
- Sibasankar Panigrahy
- Department of Materials Science and Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India.
| | - Rahul Sahu
- Center for Computational and Data Sciences, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Sandeep K Reddy
- Center for Computational and Data Sciences, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Divya Nayar
- Department of Materials Science and Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India.
| |
Collapse
|
14
|
Hajihassan Z, Afsharian NP, Ansari-Pour N. In silico engineering a CD80 variant with increased affinity to CTLA-4 and decreased affinity to CD28 for optimized cancer immunotherapy. J Immunol Methods 2023; 513:113425. [PMID: 36638881 DOI: 10.1016/j.jim.2023.113425] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Revised: 11/20/2022] [Accepted: 01/08/2023] [Indexed: 01/11/2023]
Abstract
CD80 or cluster of differentiation 80, also known as B7-1, is a member of the immunoglobulin super family, which binds to CTLA-4 and CD28 T cell receptors and induces inhibitory and inductive signals respectively. Although CTLA-4 and CD28 receptors belong to the same protein family, slight differences in their structures leads to CD80 having a higher binding affinity to CTLA-4 (-14.55 kcal/mol) compared with CD28(-12.51 kcal/mol). In this study, we constructed a variant of CD80 protein with increased binding affinity to CTLA-4 and decreased binding affinity to CD28. This variant has no signaling capability, and can act as a cap for these receptors to protect them from natural CD80 proteins existing in the body. The first step was the evolutionary and alanine scanning analysis of CD80 protein to determine conserved regions in this protein. Next, complex alanine scanning technique was employed to determine CD80 protein hotspots in CD80-CTLA-4 and CD80-CD28 protein complexes. This information was fed into a computational model developed in R for in silico mutagenesis and CD80 variant library construction. The 3D structures of variants were modeled using the Swiss model webserver. After modeling the 3D structures, HADDOCK server was employed to build all protein-protein complexes, which contain CTLA-4-CD80 variant complexes, Wild type CD80-CD28 complexes and CD28-CD80 variant complexes. Protein-protein binding free energy was determined using FoldX and the variant number 316 with mutations at 29, 31, 33 positions showed increased binding affinity to CTLA-4 (-21.43 kcal/mol) and decreased binding affinity to CD28 (- 9.54 kcal/mol). Finally, molecular dynamics (MD) simulations confirmed the stability of variant 316. In conclusion, we designed a new CD80 protein variant with potential immunotherapeutic applications.
Collapse
Affiliation(s)
- Zahra Hajihassan
- Department of Life Science Engineering, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran.
| | - Nessa Pesaran Afsharian
- Department of Life Science Engineering, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran
| | - Naser Ansari-Pour
- Department of Life Science Engineering, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran; MRC Molecular Haematology Unit, Weatherall Institute of Molecular Medicine, University of Oxford, Oxford, UK.
| |
Collapse
|
15
|
Pedersen KB, Flores-Canales JC, Schiøtt B. Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins. Proteins 2023; 91:47-61. [PMID: 35950933 PMCID: PMC10087257 DOI: 10.1002/prot.26409] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Revised: 06/17/2022] [Accepted: 07/12/2022] [Indexed: 12/29/2022]
Abstract
Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC, ff19SB/TIP4P-D, ff03CMAP/TIP4P-D, and a99SB-disp/TIP4P-disp, with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and 3 J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.
Collapse
Affiliation(s)
| | | | - Birgit Schiøtt
- Department of Chemistry, Aarhus University, Aarhus C, Denmark.,Interdisciplinary Nanoscience Center, Aarhus University, Aarhus C, Denmark
| |
Collapse
|
16
|
Chang L, Mondal A, Perez A. Towards rational computational peptide design. FRONTIERS IN BIOINFORMATICS 2022; 2:1046493. [PMID: 36338806 PMCID: PMC9634169 DOI: 10.3389/fbinf.2022.1046493] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2022] [Accepted: 10/11/2022] [Indexed: 11/16/2022] Open
Abstract
Peptides are prevalent in biology, mediating as many as 40% of protein-protein interactions, and involved in other cellular functions such as transport and signaling. Their ability to bind with high specificity make them promising therapeutical agents with intermediate properties between small molecules and large biologics. Beyond their biological role, peptides can be programmed to self-assembly, and they are already being used for functions as diverse as oligonuclotide delivery, tissue regeneration or as drugs. However, the transient nature of their interactions has limited the number of structures and knowledge of binding affinities available-and their flexible nature has limited the success of computational pipelines that predict the structures and affinities of these molecules. Fortunately, recent advances in experimental and computational pipelines are creating new opportunities for this field. We are starting to see promising predictions of complex structures, thermodynamic and kinetic properties. We believe in the following years this will lead to robust rational peptide design pipelines with success similar to those applied for small molecule drug discovery.
Collapse
Affiliation(s)
- Liwei Chang
- Department of Chemistry, University of Florida, Gainesville, FL, United States,Quantum Theory Project, University of Florida, Gainesville, FL, United States
| | - Arup Mondal
- Department of Chemistry, University of Florida, Gainesville, FL, United States,Quantum Theory Project, University of Florida, Gainesville, FL, United States
| | - Alberto Perez
- Department of Chemistry, University of Florida, Gainesville, FL, United States,Quantum Theory Project, University of Florida, Gainesville, FL, United States,*Correspondence: Alberto Perez,
| |
Collapse
|
17
|
McIvor JAP, Larsen DS, Mercadante D. Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM. J Phys Chem B 2022; 126:7833-7846. [PMID: 36125334 DOI: 10.1021/acs.jpcb.2c03837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the αR region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.
Collapse
Affiliation(s)
- Jordan A P McIvor
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Danaé S Larsen
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Davide Mercadante
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| |
Collapse
|
18
|
Sinha S, Tam B, Wang SM. Applications of Molecular Dynamics Simulation in Protein Study. MEMBRANES 2022; 12:membranes12090844. [PMID: 36135863 PMCID: PMC9505860 DOI: 10.3390/membranes12090844] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2022] [Revised: 08/24/2022] [Accepted: 08/25/2022] [Indexed: 05/29/2023]
Abstract
Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.
Collapse
|
19
|
Kasmawati H, Mustarichie R, Halimah E, Ruslin R, Arfan A, Sida NA. Unrevealing the Potential of Sansevieria trifasciata Prain Fraction for the Treatment of Androgenetic Alopecia by Inhibiting Androgen Receptors Based on LC-MS/MS Analysis, and In-Silico Studies. Molecules 2022; 27:molecules27144358. [PMID: 35889232 PMCID: PMC9318048 DOI: 10.3390/molecules27144358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Revised: 06/24/2022] [Accepted: 07/05/2022] [Indexed: 11/17/2022] Open
Abstract
Androgenetic Alopecia (AGA) occurs due to over-response to androgens causing severe hair loss on the scalp, and requires the development of new and efficient drugs to treat this condition. This study explores and identifies secondary metabolites from Sansevieriatrifasciata Prain using the LC-MS/MS and in-silico method. The inhibitory activity of bioactive compounds from S. trifasciata Prain against androgen receptors (PDB ID: 4K7A) was evaluated molecularly using docking and dynamics studies by comparing their binding energies, interactions, and stability with minoxidil. The results of the LC-MS/MS analysis identified Methyl pyrophaeophorbide A (1), Oliveramine (2), (2S)-3′, 4′-Methylenedioxy-5, 7-dimethoxyflavane (3), 1-Acetyl-β-carboline (4), Digiprolactone (5), Trichosanic acid (6) and Methyl gallate (7) from the leaves subfraction of this plant. Three alkaloid compounds (compounds 1, 3, and 4), and one flavonoid (compound 2), had lower docking scores of −7.0, −5.8, −5.2, and −6.3 kcal/mol, respectively. The prediction of binding energy using the MM-PBSA approach ensured that the potency of the four compounds was better than minoxidil, with energies of −66.13, −59.36, −40.39, and −40.25 kJ/mol for compounds 1, 3, 2, and 4, respectively. The dynamics simulation shows the stability of compound 1 based on the trajectory analysis for the 100 ns simulation. This research succeeded in identifying the compound and assessing the anti-alopecia activity of Sansevieria trifasciata Prain. Seven compounds were identified as new compounds never reported in Sansevieria trifasciata Prain. Four compounds were predicted to have better anti-alopecia activity than minoxidil in inhibiting androgen receptors through an in silico approach.
Collapse
Affiliation(s)
- Henny Kasmawati
- Doctoral Program in Pharmacy, Faculty of Pharmacy, Universitas Padjadjaran, Bandung 45363, Indonesia
- Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.)
- Correspondence: (H.K.); (R.M.)
| | - Resmi Mustarichie
- Department of Analytical Pharmacy and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Bandung 45363, Indonesia
- Correspondence: (H.K.); (R.M.)
| | - Eli Halimah
- Department of Pharmacology and Clinical Pharmacy, Faculty of Pharmacy, Universitas Padjadjaran, Bandung 45363, Indonesia;
| | - Ruslin Ruslin
- Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.)
| | - Arfan Arfan
- Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.)
| | - Nurramadhani A. Sida
- Faculty of Pharmacy, Universitas Halu Oleo, Kendari 93232, Indonesia; (R.R.); (A.A.); (N.A.S.)
| |
Collapse
|
20
|
Liu D, Qiu Y, Li Q, Zhang H. Atomistic Simulation of Lysozyme in Solutions Crowded by Tetraethylene Glycol: Force Field Dependence. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27072110. [PMID: 35408509 PMCID: PMC9000840 DOI: 10.3390/molecules27072110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/23/2022] [Revised: 03/21/2022] [Accepted: 03/22/2022] [Indexed: 11/16/2022]
Abstract
The behavior of biomolecules in crowded environments remains largely unknown due to the accuracy of simulation models and the limited experimental data for comparison. Here we chose a small crowder of tetraethylene glycol (PEG-4) to investigate the self-crowding of PEG-4 solutions and molecular crowding effects on the structure and diffusion of lysozyme at varied concentrations from dilute water to pure PEG-4 liquid. Two Amber-like force fields of Amber14SB and a99SB-disp were examined with TIP3P (fast diffusivity and low viscosity) and a99SB-disp (slow diffusivity and high viscosity) water models, respectively. Compared to the Amber14SB protein simulations, the a99SB-disp model yields more coordinated water and less PEG-4 molecules, less intramolecular hydrogen bonds (HBs), more protein-water HBs, and less protein-PEG HBs as well as stronger interactions and more hydrophilic and less hydrophobic contacts with solvent molecules. The a99SB-disp model offers comparable protein-solvent interactions in concentrated PEG-4 solutions to that in pure water. The PEG-4 crowding leads to a slow-down in the diffusivity of water, PEG-4, and protein, and the decline in the diffusion from atomistic simulations is close to or faster than the hard sphere model that neglects attractive interactions. Despite these differences, the overall structure of lysozyme appears to be maintained well at different PEG-4 concentrations for both force fields, except a slightly large deviation at 370 K at low concentrations with the a99SB-disp model. This is mainly attributed to the strong intramolecular interactions of the protein in the Amber14SB force field and to the large viscosity of the a99SB-disp water model. The results indicate that the protein force fields and the viscosity of crowder solutions affect the simulation of biomolecules under crowding conditions.
Collapse
|
21
|
Blanco MA. Computational models for studying physical instabilities in high concentration biotherapeutic formulations. MAbs 2022; 14:2044744. [PMID: 35282775 PMCID: PMC8928847 DOI: 10.1080/19420862.2022.2044744] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
Abstract
Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.
Collapse
Affiliation(s)
- Marco A. Blanco
- Materials and Biophysical Characterization, Analytical R & D, Merck & Co., Inc, Kenilworth, NJ USA
| |
Collapse
|
22
|
Heo L, Sugita Y, Feig M. Protein assembly and crowding simulations. Curr Opin Struct Biol 2022; 73:102340. [PMID: 35219215 PMCID: PMC8957576 DOI: 10.1016/j.sbi.2022.102340] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 01/07/2022] [Accepted: 01/18/2022] [Indexed: 11/17/2022]
Abstract
Proteins encounter frequent molecular interactions in biological environments. Computer simulations have become an increasingly important tool in providing mechanistic insights into how such interactions in vivo relate to their biological function. The review here focuses on simulations describing protein assembly and molecular crowding effects as two important aspects that are distinguished mainly by how specific and long-lived protein contacts are. On the topic of crowding, recent simulations have provided new insights into how crowding affects protein folding and stability, modulates enzyme activity, and affects diffusive properties. Recent studies of assembly processes focus on assembly pathways, especially for virus capsids, amyloid aggregation pathways, and the role of multivalent interactions leading to phase separation. Also, discussed are technical challenges in achieving increasingly realistic simulations of interactions in cellular environments.
Collapse
Affiliation(s)
- Lim Heo
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA. https://twitter.com/huhlim
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Computational Biophysics Research Team, RIKEN Center for Computational Science, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
23
|
Andrews B, Guerra J, Schweitzer-Stenner R, Urbanc B. Do molecular dynamics force fields accurately model Ramachandran distributions of amino acid residues in water? Phys Chem Chem Phys 2022; 24:3259-3279. [PMID: 35048087 DOI: 10.1039/d1cp05069a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Molecular dynamics (MD) is a powerful tool for studying intrinsically disordered proteins, however, its reliability depends on the accuracy of the force field. We assess Amber ff19SB, Amber ff14SB, OPLS-AA/M, and CHARMM36m with respect to their capacity to capture intrinsic conformational dynamics of 14 guest residues x (=G, A, L, V, I, F, Y, DP, EP, R, C, N, S, T) in GxG peptides in water. The MD-derived Ramachandran distribution of each guest residue is used to calculate 5 J-coupling constants and amide I' band profiles to facilitate a comparison to spectroscopic data through reduced χ2 functions. We show that the Gaussian model, optimized to best fit the experimental data, outperforms all MD force fields by an order of magnitude. The weaknesses of the MD force fields are: (i) insufficient variability of the polyproline II (pPII) population among the guest residues; (ii) oversampling of antiparallel at the expense of transitional β-strand region; (iii) inadequate sampling of turn-forming conformations for ionizable and polar residues; and (iv) insufficient guest residue-specificity of the Ramachandran distributions. Whereas Amber ff19SB performs worse than the other three force fields with respect to χ2 values, it accounts for residue-specific pPII content better than the other three force fields. Additional testing of residue-specific RSFF1 and Amber ff14SB combined with TIP4P/2005 on six guest residues x (=A, I, F, DP, R, S) reveals that residue specificity derived from protein coil libraries or an improved water model alone do not result in significantly lower χ2 values.
Collapse
Affiliation(s)
- Brian Andrews
- Department of Physics, Drexel University, Philadelphia, Pennsylvania, USA.
| | - Jose Guerra
- Department of Chemistry, Stony Brook University, Stony Brook, New York, USA
| | | | - Brigita Urbanc
- Department of Physics, Drexel University, Philadelphia, Pennsylvania, USA.
| |
Collapse
|
24
|
Chu WT, Yan Z, Chu X, Zheng X, Liu Z, Xu L, Zhang K, Wang J. Physics of biomolecular recognition and conformational dynamics. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2021; 84:126601. [PMID: 34753115 DOI: 10.1088/1361-6633/ac3800] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 11/09/2021] [Indexed: 06/13/2023]
Abstract
Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.
Collapse
Affiliation(s)
- Wen-Ting Chu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Zhiqiang Yan
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Xiakun Chu
- Department of Chemistry & Physics, State University of New York at Stony Brook, Stony Brook, NY 11794, United States of America
| | - Xiliang Zheng
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Zuojia Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Li Xu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Kun Zhang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China
| | - Jin Wang
- Department of Chemistry & Physics, State University of New York at Stony Brook, Stony Brook, NY 11794, United States of America
| |
Collapse
|
25
|
Ostrowska N, Feig M, Trylska J. Crowding affects structural dynamics and contributes to membrane association of the NS3/4A complex. Biophys J 2021; 120:3795-3806. [PMID: 34270995 PMCID: PMC8456185 DOI: 10.1016/j.bpj.2021.07.008] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2021] [Revised: 06/14/2021] [Accepted: 07/07/2021] [Indexed: 01/01/2023] Open
Abstract
Using molecular dynamics simulations, we describe how crowded environments affect the internal dynamics and diffusion of the hepatitis C virus proteases NS3/4A. This protease plays a key role in viral replication and is successfully used as a target for antiviral treatment. The NS3 enzyme requires a peptide cofactor, called NS4A, with its central part interacting with the NS3 β-sheet, and flexible, protruding terminal tails that are unstructured in water solution. The simulations describe the enzyme and water molecules at atomistic resolution, whereas crowders are modeled via either all-atom or coarse-grained models to emphasize different aspects of crowding. Crowders reflect the polyethylene glycol (PEG) molecules used in the experiments to mimic the crowded surrounding. A bead-shell model of folded coarse-grained PEG molecules considers mainly the excluded volume effect, whereas all-atom PEG models afford more protein-like crowder interactions. Circular dichroism spectroscopy experiments of the NS4A N-terminal tail show that a helical structure is formed in the presence of PEG crowders. The simulations suggest that crowding may assist in the formation of an NS4A helical fragment, positioned exactly where a transmembrane helix would fold upon the NS4A contact with the membrane. In addition, partially interactive PEGs help the NS4A N-tail to detach from the protease surface, thus enabling the process of helix insertion and potentially helping the virus establish a replication machinery needed to produce new viruses. Results point to an active role of crowding in assisting structural changes in disordered protein fragments that are necessary for their biological function.
Collapse
Affiliation(s)
- Natalia Ostrowska
- Centre of New Technologies, University of Warsaw, Warsaw, Poland,College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences, University of Warsaw, Warsaw, Poland
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, Warsaw, Poland,Corresponding author
| |
Collapse
|
26
|
Andrews B, Long K, Urbanc B. Soluble State of Villin Headpiece Protein as a Tool in the Assessment of MD Force Fields. J Phys Chem B 2021; 125:6897-6911. [PMID: 34143637 DOI: 10.1021/acs.jpcb.1c04589] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Protein self-assembly plays an important role in cellular processes. Whereas molecular dynamics (MD) represents a powerful tool in studying assembly mechanisms, its predictions depend on the accuracy of underlying force fields, which are known to overly promote protein assembly. We here examine villin headpiece domain, HP36, which remains soluble at concentrations amenable to MD studies. The experimental characterization of soluble HP36 at concentrations of 0.05 to 1 mM reveals concentration-independent 90% monomeric and 10% dimeric populations. Extensive all-atom MD simulations at two protein concentrations, 0.9 and 8.5 mM, probe the HP36 dimer population, stability, and kinetics of dimer formation within two MD force fields, Amber ff14SB and CHARMM36m. MD results demonstrate that whereas CHARMM36m captures experimental HP36 monomer populations at the lower concentration, both force fields overly promote HP36 association at the higher concentration. Moreover, contacts stabilizing HP36 dimers are force-field-dependent. CHARMM36m produces consistently higher HP36 monomer populations, lower association rates, and weaker dependence of these quantities on the protein concentration than Amber ff14SB. Nonetheless, the highest monomer populations and dissociation constants are observed when the TIP3P water model in Amber ff14SB is replaced by TIP4P/2005, showcasing the critical role of the water model in addressing the protein solubility problem in MD.
Collapse
Affiliation(s)
- Brian Andrews
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Kaho Long
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Brigita Urbanc
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
27
|
Qiu Y, Shan W, Zhang H. Force Field Benchmark of Amino Acids. 3. Hydration with Scaled Lennard-Jones Interactions. J Chem Inf Model 2021; 61:3571-3582. [PMID: 34185520 DOI: 10.1021/acs.jcim.1c00339] [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/13/2022]
Abstract
Classical protein force fields were reported with too weak protein-water interactions relative to protein-protein interactions, leading to more compact structures and artificial protein aggregation. Here we investigated the impacts of scaled Lennard-Jones (LJ) interactions on the hydration of amino acids and the simulation of folded and intrinsically disordered proteins (IDPs). The obtained optimal scaling parameters reproduce accurately hydration free energies of neutral amino acid side chain analogues and do not affect the compactness and structural stability of folded proteins significantly. The scaling leads to less compact IDPs and varies from case to case. Strengthening the interactions between protein and water oxygen or hydrogen atoms by increasing the interacting LJ well depth (ε) appears more effective than weakening protein-protein interactions by reducing the interacting dispersion coefficients (C6). We demonstrate that weakening water-water interactions is a solution as well to obtaining more favorable protein-water interactions in an indirect way, although modern force fields like Amber ff19SB and a99SB-disp tend to use water models with strong water-water interactions. This is likely a compromise between strong protein-protein interactions and strong water-water interactions. Independent optimization of protein force fields and water models is therefore needed to make both interactions more close to reality, leading to good accuracy without bias or scaling.
Collapse
Affiliation(s)
- Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Wenjie Shan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
28
|
Fakih TM, Kurniawan F, Yusuf M, Mudasir M, Tjahjono DH. Molecular Dynamics of Cobalt Protoporphyrin Antagonism of the Cancer Suppressor REV-ERBβ. Molecules 2021; 26:3251. [PMID: 34071361 PMCID: PMC8198987 DOI: 10.3390/molecules26113251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Revised: 05/15/2021] [Accepted: 05/26/2021] [Indexed: 11/17/2022] Open
Abstract
Nuclear receptor REV-ERBβ is an overexpressed oncoprotein that has been used as a target for cancer treatment. The metal-complex nature of its ligand, iron protoporphyrin IX (Heme), enables the REV-ERBβ to be used for multiple therapeutic modalities as a photonuclease, a photosensitizer, or a fluorescence imaging agent. The replacement of iron with cobalt as the metal center of protoporphyrin IX changes the ligand from an agonist to an antagonist of REV-ERBβ. The mechanism behind that phenomenon is still unclear, despite the availability of crystal structures of REV-ERBβ in complex with Heme and cobalt protoporphyrin IX (CoPP). This study used molecular dynamic simulations to compare the effects of REV-ERBβ binding to Heme and CoPP, respectively. The initial poses of Heme and CoPP in complex with agonist and antagonist forms of REV-ERBβ were predicted using molecular docking. The binding energies of each ligand were calculated using the MM/PBSA method. The computed binding affinity of Heme to REV-ERBβ was stronger than that of CoPP, in agreement with experimental results. CoPP altered the conformation of the ligand-binding site of REV-ERBβ, disrupting the binding site for nuclear receptor corepressor, which is required for REV-ERBβ to regulate the transcription of downstream target genes. Those results suggest that a subtle change in the metal center of porphyrin can change the behavior of porphyrin in cancer cell signaling. Therefore, modification of porphyrin-based agents for cancer therapy should be conducted carefully to avoid triggering unfavorable effects.
Collapse
Affiliation(s)
- Taufik Muhammad Fakih
- School of Pharmacy, Bandung Institute of Technology, Jalan Ganesha 10, Bandung 40135, Indonesia; (T.M.F.); (F.K.)
- Department of Pharmacy, Faculty of Mathematics and Natural Sciences, Universitas Islam Bandung, Jalan Rangga Gading 8, Bandung 40116, Indonesia
| | - Fransiska Kurniawan
- School of Pharmacy, Bandung Institute of Technology, Jalan Ganesha 10, Bandung 40135, Indonesia; (T.M.F.); (F.K.)
| | - Muhammad Yusuf
- Department of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Jalan Raya Bandung Sumedang Km 21, Sumedang 45363, Indonesia;
| | - Mudasir Mudasir
- Department of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Gadjah Mada, Sekip Utara BLS 21, Yogyakarta 55281, Indonesia;
| | - Daryono Hadi Tjahjono
- School of Pharmacy, Bandung Institute of Technology, Jalan Ganesha 10, Bandung 40135, Indonesia; (T.M.F.); (F.K.)
| |
Collapse
|
29
|
Docking-Based Virtual Screening and Molecular Dynamics Simulations of Quercetin Analogs as Enoyl-Acyl Carrier Protein Reductase (InhA) Inhibitors of Mycobacterium tuberculosis. Sci Pharm 2021. [DOI: 10.3390/scipharm89020020] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
The emergence of multidrug-resistant Mycobacterium tuberculosis (MTB) has become a major problem in treating tuberculosis (TB) and shows the need to develop new and efficient drugs for better TB control. This study aimed to use in silico techniques to discover potential inhibitors to the Enoyl-[acyl-carrier-protein] reductase (InhA), which controls mycobacterial cell wall construction. Initially, 391 quercetin analogs present in the KNApSAck_3D database were selected, filters were sequentially applied by docking-based virtual screening. After recategorizing the variables (bond energy prediction and molecular interaction, including hydrogen bond and hydrophobic bond), compounds C00013874, C00006532, and C00013887 were selected as hit ligands. These compounds showed great hydrophobic contributions, and for each hit ligand, 100 ns of molecular dynamic simulations were performed, and the binding free energy was calculated. C00013874 demonstrated the greatest capacity for the InhA enzyme inhibition with ΔGbind = −148.651 kcal/mol compare to NAD (native ligand) presented a ΔGbind = −87.570 kcal/mol. These data are preliminary studies and might be a suitable candidate for further experimental analysis.
Collapse
|
30
|
Assessment of transferable forcefields for protein simulations attests improved description of disordered states and secondary structure propensities, and hints at multi-protein systems as the next challenge for optimization. Comput Struct Biotechnol J 2021; 19:2626-2636. [PMID: 34025949 PMCID: PMC8120800 DOI: 10.1016/j.csbj.2021.04.050] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Revised: 04/20/2021] [Accepted: 04/22/2021] [Indexed: 01/02/2023] Open
Abstract
Continuous assessment of transferable forcefields for molecular simulations is essential to identify their weaknesses and direct improvement efforts. The latest efforts focused on better describing disordered proteins while retaining proper description of folded domains, important because forcefields of the previous generations produce overly compact disordered states. Such improvements should additionally alleviate the related problem of over-stabilized protein–protein interactions, which has been largely overlooked. Here we evaluated three state-of-the-art forcefields, current flagships of their respective developers, optimized for ordered and disordered proteins: CHARMM36m with its recommended corrected TIP3P* water, ff19SB with the recommended OPC water, and the 2019 a99SBdisp forcefield by D. E. Shaw Research with its modified TIP4P water; plus ff14SB with TIP3P as an example of the former generation of forcefields. Our evaluation entailed simulations of (i) multiple copies of a protein that is highly soluble yet undergoes weak dimerization, (ii) a disordered peptide with low, well-characterized alpha helical propensity, and (iii) a peptide known to form insoluble β-aggregates. Our results recapitulate ff14SB-TIP3P over-stabilizing aggregates and secondary structures and place a99SBdisp-TIP4PD at the other end i.e. predicting overly weak intermolecular interactions despite reasonably predicting secondary structure propensities. In-between, CHARMM36m-TIP3P* still over-stabilizes aggregates but predicts residue-wise alpha helical propensities in solution slightly better than ff19SB-OPC, while ff19SB-OPC poses the best prediction of weak dimerization of the soluble protein still predicting aggregation of the β-peptides. This independent assessment shows that the claimed forcefield improvements are real, but also that a right balance between noncovalent attraction and repulsion has not yet been reached. We thus propose developers to consider systems like those tested here in their forcefield tuning protocols. Last, the good performance of CHARMM36m-TIP3P* further shows that tuning 3-point water models might still be an alternative to the more costly 4-point models like OPC and TIP4PD.
Collapse
|
31
|
Souza PCT, Alessandri R, Barnoud J, Thallmair S, Faustino I, Grünewald F, Patmanidis I, Abdizadeh H, Bruininks BMH, Wassenaar TA, Kroon PC, Melcr J, Nieto V, Corradi V, Khan HM, Domański J, Javanainen M, Martinez-Seara H, Reuter N, Best RB, Vattulainen I, Monticelli L, Periole X, Tieleman DP, de Vries AH, Marrink SJ. Martini 3: a general purpose force field for coarse-grained molecular dynamics. Nat Methods 2021; 18:382-388. [PMID: 33782607 DOI: 10.1038/s41592-021-01098-3] [Citation(s) in RCA: 435] [Impact Index Per Article: 145.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 02/22/2021] [Indexed: 01/31/2023]
Abstract
The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.
Collapse
Affiliation(s)
- Paulo C T Souza
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands. .,Molecular Microbiology and Structural Biochemistry, UMR 5086 CNRS and University of Lyon, Lyon, France.
| | - Riccardo Alessandri
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Jonathan Barnoud
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands.,Intangible Realities Laboratory, University of Bristol, School of Chemistry, Bristol, UK
| | - Sebastian Thallmair
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands.,Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
| | - Ignacio Faustino
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Fabian Grünewald
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Ilias Patmanidis
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Haleh Abdizadeh
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Bart M H Bruininks
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Peter C Kroon
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Josef Melcr
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Vincent Nieto
- Molecular Microbiology and Structural Biochemistry, UMR 5086 CNRS and University of Lyon, Lyon, France
| | - Valentina Corradi
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada
| | - Hanif M Khan
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada.,Department of Chemistry and Computational Biology Unit, University of Bergen, Bergen, Norway
| | - Jan Domański
- Department of Biochemistry, University of Oxford, Oxford, UK.,Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, USA
| | - Matti Javanainen
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic.,Computational Physics Laboratory, Tampere University, Tampere, Finland
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Nathalie Reuter
- Department of Chemistry and Computational Biology Unit, University of Bergen, Bergen, Norway
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, USA
| | - Ilpo Vattulainen
- Computational Physics Laboratory, Tampere University, Tampere, Finland.,Department of Physics, University of Helsinki, Helsinki, Finland
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry, UMR 5086 CNRS and University of Lyon, Lyon, France
| | - Xavier Periole
- Department of Chemistry, Aarhus University, Aarhus C, Denmark
| | - D Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada
| | - Alex H de Vries
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Material, University of Groningen, Groningen, the Netherlands.
| |
Collapse
|
32
|
Zou J, Xiao S, Simmerling C, Raleigh DP. Quantitative Analysis of Protein Unfolded State Energetics: Experimental and Computational Studies Demonstrate That Non-Native Side-Chain Interactions Stabilize Local Native Backbone Structure. J Phys Chem B 2021; 125:3269-3277. [PMID: 33779182 DOI: 10.1021/acs.jpcb.0c08922] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Proteins fold on relatively smooth free energy landscapes which are biased toward the native state, but even simple topologies which fold rapidly can experience roughness on their free energy landscape. The details of these interactions are difficult to elucidate experimentally. Closely related to the problem of deciphering the details of the free energy landscape is the problem of defining the interactions in the denatured state ensemble (DSE) which is populated under native conditions, that is, under conditions where the native state is stable. The DSE of many proteins deviates from random coil models, but quantifying and defining the energetics of the transiently populated interactions in this ensemble is extremely challenging. Characterization of the DSE of proteins which fold to compact structures is also relevant to studies of intrinsically disordered proteins (IDPs) since interactions in the dynamic ensemble populated by IDPs can modulate their behavior. Here we show how experimental thermodynamic and pKa measurements can be combined with computational thermodynamic integration to quantify interactions in the DSE. We show that non-native side chain interactions can stabilize native backbone structure in the DSE and demonstrate that that even rapidly folding proteins can form energetically significant non-native interactions in their DSE. As an example, we characterize a non-native salt bridge that stabilizes local native backbone structure in the DSE of a widely studied fast-folding protein, the villin headpiece helical domain. The combined computational experimental approach is applicable to other protein unfolded states and provides insight that is impossible to obtain with either method alone.
Collapse
Affiliation(s)
- Junjie Zou
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794-3400, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-3400, United States
| | - Shifeng Xiao
- Shenzhen Key Laboratory of Marine Biotechnology and Ecology, College of Life Sciences and Oceanography, Shenzhen University, Shenzhen 518060, China
| | - Carlos Simmerling
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794-3400, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-3400, United States
| | - Daniel P Raleigh
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794-3400, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-3400, United States
| |
Collapse
|
33
|
Gotsmy M, Escalona Y, Oostenbrink C, Petrov D. Exploring the structure and dynamics of proteins in soil organic matter. Proteins 2021; 89:925-936. [PMID: 33675059 PMCID: PMC8360018 DOI: 10.1002/prot.26070] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Revised: 01/27/2021] [Accepted: 02/23/2021] [Indexed: 01/17/2023]
Abstract
Alongside inorganic materials, water, and air, soil organic matter (SOM) is one of the major components of soil and has tremendous influence on the environment given its vital role in the carbon cycle. Many soil dwelling organisms like plants, fungi and bacteria excrete proteins, whose interaction with SOM is poorly understood on an atomistic level. In this study, molecular dynamics simulations were used to investigate selected proteins in soil models of different complexity from simple co-solvent molecules to Leonardite humic acids (LHA). We analyzed the proteins in terms of their structural stability, the nature and strength of the interactions with their surroundings, as well as their aggregation behavior. Upon insertion of proteins in complex SOM models, their structural stability decreased, although no unfolding or disruption of secondary structure was observed. The interactions of proteins and SOM were primarily governed by electrostatic forces, often in form of hydrogen bonds. However, also weaker van der Waals forces made a significant contribution to the total interaction energies. Moreover, we showed that even though the molecular structure and size of SOM molecules varied, the functional groups of SOM ordered around the protein in a similar pattern. Finally, the number of aggregates formed by proteins and SOM molecules was shown to be primarily proportional to the size of the latter. Strikingly, for varying protein net charges no changes in the formation of aggregates with the strongly negatively charged LHA were observed.
Collapse
Affiliation(s)
- Mathias Gotsmy
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | - Yerko Escalona
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | - Drazen Petrov
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| |
Collapse
|
34
|
Pedebos C, Smith IPS, Boags A, Khalid S. The hitchhiker's guide to the periplasm: Unexpected molecular interactions of polymyxin B1 in E. coli. Structure 2021; 29:444-456.e2. [PMID: 33577754 DOI: 10.1016/j.str.2021.01.009] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2020] [Revised: 12/11/2020] [Accepted: 01/21/2021] [Indexed: 12/19/2022]
Abstract
The periplasm of Gram-negative bacteria is a complex, highly crowded molecular environment. Little is known about how antibiotics move across the periplasm and the interactions they experience. Here, atomistic molecular dynamics simulations are used to study the antibiotic polymyxin B1 within models of the periplasm, which are crowded to different extents. We show that PMB1 is likely to be able to "hitchhike" within the periplasm by binding to lipoprotein carriers-a previously unreported passive transport route. The simulations reveal that PMB1 forms both transient and long-lived interactions with proteins, osmolytes, lipids of the outer membrane, and the cell wall, and is rarely uncomplexed when in the periplasm. Furthermore, it can interfere in the conformational dynamics of native proteins. These are important considerations for interpreting its mechanism of action and are likely to also hold for other antibiotics that rely on diffusion to cross the periplasm.
Collapse
Affiliation(s)
- Conrado Pedebos
- School of Chemistry, University of Southampton, Highfield Campus, Southampton SO17 1BJ, UK
| | - Iain Peter Shand Smith
- School of Chemistry, University of Southampton, Highfield Campus, Southampton SO17 1BJ, UK
| | - Alister Boags
- School of Chemistry, University of Southampton, Highfield Campus, Southampton SO17 1BJ, UK
| | - Syma Khalid
- School of Chemistry, University of Southampton, Highfield Campus, Southampton SO17 1BJ, UK.
| |
Collapse
|
35
|
Effect of the Water Model in Simulations of Protein-Protein Recognition and Association. Polymers (Basel) 2021; 13:polym13020176. [PMID: 33419008 PMCID: PMC7825341 DOI: 10.3390/polym13020176] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Revised: 12/24/2020] [Accepted: 12/29/2020] [Indexed: 12/24/2022] Open
Abstract
We study self-association of ubiquitin and the disordered protein ACTR using the most commonly used water models. We find that dissociation events are found only with TIP4P-EW and TIP4P/2005, while the widely used TIP3P water model produces straightforward aggregation of the molecules due to the absence of dissociation events. We also find that TIP4P/2005 is the only water model that reproduces the fast association/dissociation dynamics of ubiquitin and best identifies its binding surface. Our results show the critical role of the water model in the description of protein–protein interactions and binding.
Collapse
|
36
|
Gopan G, Gruebele M, Rickard M. In-cell protein landscapes: making the match between theory, simulation and experiment. Curr Opin Struct Biol 2020; 66:163-169. [PMID: 33254078 DOI: 10.1016/j.sbi.2020.10.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Accepted: 10/10/2020] [Indexed: 11/26/2022]
Abstract
Theory, computation and experiment have matched up for the folding of small proteins in vitro, a difficult feat because folding energy landscapes are fairly smooth and free energy differences between states are small. Smoothness means that protein structure and folding are susceptible to the local environment inside living cells. Theory, computation and experiment are now exploring cellular modulation of energy landscapes. Interesting concepts have emerged, such as co-evolution of protein surfaces with their cellular environment to reduce detrimental interactions. Here we look at very recent work beginning to bring together theory, simulations and experiments in the area of protein landscape modulation, to see what problems might be solved in the near future by combining these approaches.
Collapse
Affiliation(s)
- Gopika Gopan
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Martin Gruebele
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA; Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
| | - Meredith Rickard
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
37
|
Popielec A, Ostrowska N, Wojciechowska M, Feig M, Trylska J. Crowded environment affects the activity and inhibition of the NS3/4A protease. Biochimie 2020; 176:169-180. [DOI: 10.1016/j.biochi.2020.07.009] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2020] [Revised: 07/17/2020] [Accepted: 07/17/2020] [Indexed: 12/18/2022]
|
38
|
Banerjee P, Lipowsky R, Santer M. Coarse-Grained Molecular Model for the Glycosylphosphatidylinositol Anchor with and without Protein. J Chem Theory Comput 2020; 16:3889-3903. [PMID: 32392421 PMCID: PMC7303967 DOI: 10.1021/acs.jctc.0c00056] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2020] [Indexed: 01/17/2023]
Abstract
Glycosylphosphatidylinositol (GPI) anchors are a unique class of complex glycolipids that anchor a great variety of proteins to the extracellular leaflet of plasma membranes of eukaryotic cells. These anchors can exist either with or without an attached protein called GPI-anchored protein (GPI-AP) both in vitro and in vivo. Although GPIs are known to participate in a broad range of cellular functions, it is to a large extent unknown how these are related to GPI structure and composition. Their conformational flexibility and microheterogeneity make it difficult to study them experimentally. Simplified atomistic models are amenable to all-atom computer simulations in small lipid bilayer patches but not suitable for studying their partitioning and trafficking in complex and heterogeneous membranes. Here, we present a coarse-grained model of the GPI anchor constructed with a modified version of the MARTINI force field that is suited for modeling carbohydrates, proteins, and lipids in an aqueous environment using MARTINI's polarizable water. The nonbonded interactions for sugars were reparametrized by calculating their partitioning free energies between polar and apolar phases. In addition, sugar-sugar interactions were optimized by adjusting the second virial coefficients of osmotic pressures for solutions of glucose, sucrose, and trehalose to match with experimental data. With respect to the conformational dynamics of GPI-anchored green fluorescent protein, the accessible time scales are now at least an order of magnitude larger than for the all-atom system. This is particularly important for fine-tuning the mutual interactions of lipids, carbohydrates, and amino acids when comparing to experimental results. We discuss the prospective use of the coarse-grained GPI model for studying protein-sorting and trafficking in membrane models.
Collapse
Affiliation(s)
- Pallavi Banerjee
- Max
Planck Institute of Colloids and Interfaces, Potsdam 14476, Germany
- Institute
of Biochemistry and Biology, University
of Potsdam, Potsdam 14469, Germany
| | - Reinhard Lipowsky
- Max
Planck Institute of Colloids and Interfaces, Potsdam 14476, Germany
- Institute
of Biochemistry and Biology, University
of Potsdam, Potsdam 14469, Germany
| | - Mark Santer
- Max
Planck Institute of Colloids and Interfaces, Potsdam 14476, Germany
| |
Collapse
|
39
|
Timr S, Madern D, Sterpone F. Protein thermal stability. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:239-272. [PMID: 32145947 DOI: 10.1016/bs.pmbts.2019.12.007] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.
Collapse
Affiliation(s)
- Stepan Timr
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | | | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
40
|
Honegger P, Steinhauser O. The nuclear Overhauser Effect (NOE) as a tool to study macromolecular confinement: Elucidation and disentangling of crowding and encapsulation effects. J Chem Phys 2020; 152:024120. [PMID: 31941328 DOI: 10.1063/1.5135816] [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/14/2022] Open
Abstract
We propose a methodology to capture short-lived but biophysically important contacts of biomacromolecules using the biomolecule-water nuclear Overhauser effect as an indirect microscope. Thus, instead of probing the direct correlation with the foreign biomolecule, we detect its presence by the disturbance it causes in the surrounding water. In addition, this information obtained is spatially resolved and can thus be attributed to specific sites. We extend this approach to the influence of more than one change in chemical environment and show a methodological way of resolution. This is achieved by taking double differences of corresponding σNOE/σROE ratios of the systems studied and separating specific, unspecific, and intermediate influence. While applied to crowding and encapsulation in this study, this method is generally suitable for any combination of changes in chemical environment.
Collapse
Affiliation(s)
- Philipp Honegger
- University of Vienna, Faculty of Chemistry, Department of Computational Biological Chemistry, Währingerstraße 17, A-1090 Vienna, Austria
| | - Othmar Steinhauser
- University of Vienna, Faculty of Chemistry, Department of Computational Biological Chemistry, Währingerstraße 17, A-1090 Vienna, Austria
| |
Collapse
|
41
|
Honegger P, Steinhauser O. The protein-water nuclear Overhauser effect (NOE) as an indirect microscope for molecular surface mapping of interaction patterns. Phys Chem Chem Phys 2019; 22:212-222. [PMID: 31799520 DOI: 10.1039/c9cp04752b] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In this computational study, the intermolecular solute-solvent Nuclear Overhauser Effect (NOE) of the model protein ubiquitin in different chemical environments (free, bound to a partner protein and encapsulated) is investigated. Short-ranged NOE observables such as the NOE/ROE ratio reveal hydration phenomena on absolute timescales such as fast hydration sites and slow water clefts. We demonstrate the ability of solute-solvent NOE differences measured of the same protein in different chemical environments to reveal hydration changes on the relative timescale. The resulting NOE/ROE-surface maps are shown to be a central key for analyzing biologically relevant chemical influences such as complexation and confinement: the presence of a complexing macromolecule or a confining surface wall modulates the water mobility in the vicinity of the probe protein, hence revealing which residues of said protein are proximate to the foreign interface and which are chemically unaffected. This way, hydration phenomena can serve to indirectly map the precise influence (position) of other molecules or interfaces onto the protein surface. This proposed one-protein many-solvents approach may offer experimental benefits over classical one-protein other-protein pseudo-intermolecular transient NOEs. Furthermore, combined influences such as complexation and confinement may exert non-additive influences on the protein compared to a reference state. We offer a mathematical method to disentangle the influence of these two different chemical environments.
Collapse
Affiliation(s)
- Philipp Honegger
- University of Vienna, Faculty of Chemistry, Department of Computational Biological Chemistry, Währingerstr. 17, A-1090 Vienna, Austria.
| | | |
Collapse
|
42
|
Emani PS, Yimer YY, Davidowski SK, Gebhart RN, Ferreira HE, Kuprov I, Pfaendtner J, Drobny GP. Combining Molecular and Spin Dynamics Simulations with Solid-State NMR: A Case Study of Amphiphilic Lysine-Leucine Repeat Peptide Aggregates. J Phys Chem B 2019; 123:10915-10929. [PMID: 31769684 DOI: 10.1021/acs.jpcb.9b09245] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Interpreting dynamics in solid-state molecular systems requires characterization of the potentially heterogeneous environmental contexts of molecules. In particular, the analysis of solid-state nuclear magnetic resonance (ssNMR) data to elucidate molecular dynamics (MD) involves modeling the restriction to overall tumbling by neighbors, as well as the concentrations of water and buffer. In this exploration of the factors that influence motion, we utilize atomistic MD trajectories of peptide aggregates with varying hydration to mimic an amorphous solid-state environment and predict ssNMR relaxation rates. We also account for spin diffusion in multiply spin-labeled (up to 19 nuclei) residues, with several models of dipolar-coupling networks. The framework serves as a general approach to determine essential spin couplings affecting relaxation, benchmark MD force fields, and reveal the hydration dependence of dynamics in a crowded environment. We demonstrate the methodology on a previously characterized amphiphilic 14-residue lysine-leucine repeat peptide, LKα14 (Ac-LKKLLKLLKKLLKL-c), which has an α-helical secondary structure and putatively forms leucine-burying tetramers in the solid state. We measure the R1 relaxation rates of uniformly 13C-labeled and site-specific 2H-labeled leucines in the hydrophobic core of LKα14 at multiple hydration levels. Studies of 9 and 18 tetramer bundles reveal the following: (a) for the incoherent component of 13C relaxation, the nearest-neighbor spin interactions dominate, while the 1H-1H interactions have minimal impact; (b) the AMBER ff14SB dihedral barriers for the leucine Cγ-Cδ bond ("methyl rotation barriers") must be lowered by a factor of 0.7 to better match the 2H data; (c) proton-driven spin diffusion explains some of the discrepancy between experimental and simulated rates for the Cβ and Cα nuclei; and (d) 13C relaxation rates are mostly underestimated in the MD simulations at all hydrations, and the discrepancies identify likely motions missing in the 50 ns MD trajectories.
Collapse
Affiliation(s)
- Prashant S Emani
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Yeneneh Y Yimer
- Department of Chemical Engineering , University of Washington , 105 Benson Hall, Box 351750 , Seattle , Washington 98195-1750 , United States
| | - Stephen K Davidowski
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Rachel N Gebhart
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Helen E Ferreira
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Ilya Kuprov
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , U.K
| | - Jim Pfaendtner
- Department of Chemical Engineering , University of Washington , 105 Benson Hall, Box 351750 , Seattle , Washington 98195-1750 , United States
| | - Gary P Drobny
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| |
Collapse
|
43
|
Bashardanesh Z, Elf J, Zhang H, van der Spoel D. Rotational and Translational Diffusion of Proteins as a Function of Concentration. ACS OMEGA 2019; 4:20654-20664. [PMID: 31858051 PMCID: PMC6906769 DOI: 10.1021/acsomega.9b02835] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/02/2019] [Accepted: 11/14/2019] [Indexed: 05/21/2023]
Abstract
Atomistic simulations of three different proteins at different concentrations are performed to obtain insight into protein mobility as a function of protein concentration. We report on simulations of proteins from diluted to the physiological water concentration (about 70% of the mass). First, the viscosity was computed and found to increase by a factor of 7-9 going from pure water to the highest protein concentration, in excellent agreement with in vivo nuclear magnetic resonance results. At a physiological concentration of proteins, the translational diffusion is found to be slowed down to about 30% of the in vitro values. The slow-down of diffusion found here using atomistic models is slightly more than that of a hard sphere model that neglects the electrostatic interactions. Interestingly, rotational diffusion of proteins is slowed down somewhat more (by about 80-95% compared to in vitro values) than translational diffusion, in line with experimental findings and consistent with the increased viscosity. The finding that rotation is retarded more than translation is attributed to solvent-separated clustering. No direct interactions between the proteins are found, and the clustering can likely be attributed to dispersion interactions that are stronger between proteins than between protein and water. Based on these simulations, we can also conclude that the internal dynamics of the proteins in our study are affected only marginally under crowding conditions, and the proteins become somewhat more stable at higher concentrations. Simulations were performed using a force field that was tuned for dealing with crowding conditions by strengthening the protein-water interactions. This force field seems to lead to a reproducible partial unfolding of an α-helix in one of the proteins, an effect that was not observed in the unmodified force field.
Collapse
Affiliation(s)
- Zahedeh Bashardanesh
- Uppsala
Center for Computational Chemistry, Science for Life Laboratory, Department
of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| | - Johan Elf
- Uppsala
Center for Computational Chemistry, Science for Life Laboratory, Department
of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| | - Haiyang Zhang
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology
Beijing, 100083 Beijing, China
| | - David van der Spoel
- Uppsala
Center for Computational Chemistry, Science for Life Laboratory, Department
of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
- E-mail: . Phone: +46 18 4714205
| |
Collapse
|
44
|
Oliveira Bortot L, Bashardanesh Z, van der Spoel D. Making Soup: Preparing and Validating Models of the Bacterial Cytoplasm for Molecular Simulation. J Chem Inf Model 2019; 60:322-331. [DOI: 10.1021/acs.jcim.9b00971] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Leandro Oliveira Bortot
- Laboratory of Biological Physics, School of Pharmaceutical Sciences of Ribeirão Preto, University of São Paulo, Avenida do Café s/n, 14040-903 Ribeirão Preto-SP, Brazil
| | - Zahedeh Bashardanesh
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
45
|
Tian C, Kasavajhala K, Belfon KAA, Raguette L, Huang H, Migues AN, Bickel J, Wang Y, Pincay J, Wu Q, Simmerling C. ff19SB: Amino-Acid-Specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution. J Chem Theory Comput 2019; 16:528-552. [PMID: 31714766 DOI: 10.1021/acs.jctc.9b00591] [Citation(s) in RCA: 834] [Impact Index Per Article: 166.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Molecular dynamics (MD) simulations have become increasingly popular in studying the motions and functions of biomolecules. The accuracy of the simulation, however, is highly determined by the molecular mechanics (MM) force field (FF), a set of functions with adjustable parameters to compute the potential energies from atomic positions. However, the overall quality of the FF, such as our previously published ff99SB and ff14SB, can be limited by assumptions that were made years ago. In the updated model presented here (ff19SB), we have significantly improved the backbone profiles for all 20 amino acids. We fit coupled φ/ψ parameters using 2D φ/ψ conformational scans for multiple amino acids, using as reference data the entire 2D quantum mechanics (QM) energy surface. We address the polarization inconsistency during dihedral parameter fitting by using both QM and MM in aqueous solution. Finally, we examine possible dependency of the backbone fitting on side chain rotamer. To extensively validate ff19SB parameters, and to compare to results using other Amber models, we have performed a total of ∼5 ms MD simulations in explicit solvent. Our results show that after amino-acid-specific training against QM data with solvent polarization, ff19SB not only reproduces the differences in amino-acid-specific Protein Data Bank (PDB) Ramachandran maps better but also shows significantly improved capability to differentiate amino-acid-dependent properties such as helical propensities. We also conclude that an inherent underestimation of helicity is present in ff14SB, which is (inexactly) compensated for by an increase in helical content driven by the TIP3P bias toward overly compact structures. In summary, ff19SB, when combined with a more accurate water model such as OPC, should have better predictive power for modeling sequence-specific behavior, protein mutations, and also rational protein design. Of the explicit water models tested here, we recommend use of OPC with ff19SB.
Collapse
Affiliation(s)
- Chuan Tian
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Koushik Kasavajhala
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Kellon A A Belfon
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Lauren Raguette
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - He Huang
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Angela N Migues
- Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - John Bickel
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Yuzhang Wang
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Jorge Pincay
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Qin Wu
- Center for Functional Nanomaterials , Brookhaven National Laboratory , Upton , New York 11973 , United States
| | - Carlos Simmerling
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| |
Collapse
|
46
|
Tolmachev DA, Boyko OS, Lukasheva NV, Martinez-Seara H, Karttunen M. Overbinding and Qualitative and Quantitative Changes Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Force Fields with and without NBFIX and ECC Corrections. J Chem Theory Comput 2019; 16:677-687. [DOI: 10.1021/acs.jctc.9b00813] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- D. A. Tolmachev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - O. S. Boyko
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - N. V. Lukasheva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - H. Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Mikko Karttunen
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| |
Collapse
|
47
|
Perez CP, Elmore DE, Radhakrishnan ML. Computationally Modeling Electrostatic Binding Energetics in a Crowded, Dynamic Environment: Physical Insights from a Peptide–DNA System. J Phys Chem B 2019; 123:10718-10734. [DOI: 10.1021/acs.jpcb.9b09478] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
48
|
Alessandri R, Souza PCT, Thallmair S, Melo MN, de Vries AH, Marrink SJ. Pitfalls of the Martini Model. J Chem Theory Comput 2019; 15:5448-5460. [PMID: 31498621 PMCID: PMC6785803 DOI: 10.1021/acs.jctc.9b00473] [Citation(s) in RCA: 123] [Impact Index Per Article: 24.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The computational and conceptual simplifications realized by coarse-grain (CG) models make them a ubiquitous tool in the current computational modeling landscape. Building block based CG models, such as the Martini model, possess the key advantage of allowing for a broad range of applications without the need to reparametrize the force field each time. However, there are certain inherent limitations to this approach, which we investigate in detail in this work. We first study the consequences of the absence of specific cross Lennard-Jones parameters between different particle sizes. We show that this lack may lead to artificially high free energy barriers in dimerization profiles. We then look at the effect of deviating too far from the standard bonded parameters, both in terms of solute partitioning behavior and solvent properties. Moreover, we show that too weak bonded force constants entail the risk of artificially inducing clustering, which has to be taken into account when designing elastic network models for proteins. These results have implications for the current use of the Martini CG model and provide clear directions for the reparametrization of the Martini model. Moreover, our findings are generally relevant for the parametrization of any other building block based force field.
Collapse
Affiliation(s)
- Riccardo Alessandri
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| | - Paulo C T Souza
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| | - Sebastian Thallmair
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| | - Manuel N Melo
- Instituto de Tecnologia Qúımica Bioĺogica Ant́onio Xavier , Universidade Nova de Lisboa , Av. da Reṕublica , 2780-157 Oeiras , Portugal
| | - Alex H de Vries
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| |
Collapse
|
49
|
Rickard MM, Zhang Y, Gruebele M, Pogorelov TV. In-Cell Protein-Protein Contacts: Transient Interactions in the Crowd. J Phys Chem Lett 2019; 10:5667-5673. [PMID: 31483661 DOI: 10.1021/acs.jpclett.9b01556] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Proteins in vivo are immersed in a crowded environment of water, ions, metabolites, and macromolecules. In-cell experiments highlight how transient weak protein-protein interactions promote (via functional "quinary structure") or hinder (via competitive binding or "sticking") complex formation. Computational models of the cytoplasm are expensive. We tackle this challenge with an all-atom model of a small volume of the E. coli cytoplasm to simulate protein-protein contacts up to the 5 μs time scale on the special-purpose supercomputer Anton 2. We use three CHARMM-derived force fields: C22*, C36m, and C36mCU (with CUFIX corrections). We find that both C36m and C36mCU form smaller contact surfaces than C22*. Although CUFIX was developed to reduce protein-protein sticking, larger contacts are observed with C36mCU than C36m. We show that the lifespan Δt of protein-protein contacts obeys a power law distribution between 0.03 and 3 μs, with ∼90% of all contacts lasting <1 μs (similar to the time scale for downhill folding).
Collapse
Affiliation(s)
- Meredith M Rickard
- Department of Chemistry , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
| | - Yi Zhang
- Center for Biophysics and Computational Biology , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
| | - Martin Gruebele
- Department of Chemistry , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
- Center for Biophysics and Computational Biology , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
- Department of Physics , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
| | - Taras V Pogorelov
- Department of Chemistry , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
- Center for Biophysics and Computational Biology , University of Illinois, Urbana-Champaign , Urbana , Illinois 61801 , United States
| |
Collapse
|
50
|
Ostrowska N, Feig M, Trylska J. Modeling Crowded Environment in Molecular Simulations. Front Mol Biosci 2019; 6:86. [PMID: 31572730 PMCID: PMC6749006 DOI: 10.3389/fmolb.2019.00086] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Accepted: 08/27/2019] [Indexed: 01/09/2023] Open
Abstract
Biomolecules perform their various functions in living cells, namely in an environment that is crowded by many macromolecules. Thus, simulating the dynamics and interactions of biomolecules should take into account not only water and ions but also other binding partners, metabolites, lipids and macromolecules found in cells. In the last decade, research on how to model macromolecular crowders around proteins in order to simulate their dynamics in models of cellular environments has gained a lot of attention. In this mini-review we focus on the models of crowding agents that have been used in computer modeling studies of proteins and peptides, especially via molecular dynamics simulations.
Collapse
Affiliation(s)
- Natalia Ostrowska
- Centre of New Technologies, University of Warsaw, Warsaw, Poland.,College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences, University of Warsaw, Warsaw, Poland
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, United States
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| |
Collapse
|