1
|
Mioduszewski Ł. Choosing the right density for a concentrated protein system like gluten in a coarse-grained model. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2023; 52:583-591. [PMID: 37378869 PMCID: PMC10618313 DOI: 10.1007/s00249-023-01667-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Revised: 05/05/2023] [Accepted: 05/24/2023] [Indexed: 06/29/2023]
Abstract
Large coarse-grained simulations are often conducted with an implicit solvent, which makes it hard to assess the water content of the sample and the effective concentration of the system. Here the number and the size of cavities and entanglements in the system, together with density profiles, are used to asses the homogeneity and interconnectedness of gluten. This is a continuation of an earlier article, "Viscoelastic properties of wheat gluten in a molecular dynamics study" (Mioduszewski and Cieplak 2021b). It turns out there is a wide range of densities (between 1 residue per cubic nanometer and 3 residues/nm[Formula: see text]) where the system is interconnected, but not homogeneous: there are still large empty spaces, surrounded by an entangled protein network. Those findings should be of importance to any coarse-grained simulation of large protein systems.
Collapse
Affiliation(s)
- Łukasz Mioduszewski
- Faculty of Mathematics and Natural Sciences, Cardinal Stefan Wyszyński University, Wóycickiego 1/3, 01-938, Warsaw, Poland.
| |
Collapse
|
2
|
Farnudi A, Ejtehadi MR, Everaers R. Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization. Phys Rev E 2023; 108:015301. [PMID: 37583159 DOI: 10.1103/physreve.108.015301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Accepted: 05/26/2023] [Indexed: 08/17/2023]
Abstract
Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.
Collapse
Affiliation(s)
- Ali Farnudi
- Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran
| | - Mohammad Reza Ejtehadi
- Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran
| | - Ralf Everaers
- Ecole Normale Supérieure (ENS) de Lyon, CNRS, Laboratoire de Physique and Centre Blaise Pascal de l'ENS de Lyon, F-69342 Lyon, France
| |
Collapse
|
3
|
Kumawat A, Namsani S, Pramanik D, Roy S, Singh JK. Integrated docking and enhanced sampling-based selection of repurposing drugs for SARS-CoV-2 by targeting host dependent factors. J Biomol Struct Dyn 2022; 40:9897-9908. [PMID: 34155961 PMCID: PMC8220434 DOI: 10.1080/07391102.2021.1937319] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Since the onset of global pandemic, the most focused research currently in progress is the development of potential drug candidates and clinical trials of existing FDA approved drugs for other relevant diseases, in order to repurpose them for the COVID-19. At the same time, several high throughput screenings of drugs have been reported to inhibit the viral components during the early course of infection but with little proven efficacies. Here, we investigate the drug repurposing strategies to counteract the coronavirus infection which involves several potential targetable host proteins involved in viral replication and disease progression. We report the high throughput analysis of literature-derived repurposing drug candidates that can be used to target the genetic regulators known to interact with viral proteins based on experimental and interactome studies. In this work we have performed integrated molecular docking followed by molecular dynamics (MD) simulations and free energy calculations through an expedite in silico process where the number of screened candidates reduces sequentially at every step based on physicochemical interactions. We elucidate that in addition to the pre-clinical and FDA approved drugs that targets specific regulatory proteins, a range of chemical compounds (Nafamostat, Chloramphenicol, Ponatinib) binds to the other gene transcription and translation regulatory proteins with higher affinity and may harbour potential for therapeutic uses. There is a rapid growing interest in the development of combination therapy for COVID-19 to target multiple enzymes/pathways. Our in silico approach would be useful in generating leads for experimental screening for rapid drug repurposing against SARS-CoV-2 interacting host proteins.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Amit Kumawat
- Prescience Insilico Private Limited, Bangalore, India,Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| | | | - Debabrata Pramanik
- Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| | - Sudip Roy
- Prescience Insilico Private Limited, Bangalore, India,CONTACT Sudip Roy ;
| | - Jayant K. Singh
- Prescience Insilico Private Limited, Bangalore, India,Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India,Jayant K. Singh
| |
Collapse
|
4
|
Abstract
Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.
Collapse
Affiliation(s)
- Serge Perez
- Centre de Recherche sur les Macromolecules Vegetales, University of Grenoble-Alpes, Centre National de la Recherche Scientifique, Grenoble F-38041, France
| | - Olga Makshakova
- FRC Kazan Scientific Center of Russian Academy of Sciences, Kazan Institute of Biochemistry and Biophysics, Kazan 420111, Russia
| |
Collapse
|
5
|
Martini 3 Model of Cellulose Microfibrils: On the Route to Capture Large Conformational Changes of Polysaccharides. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27030976. [PMID: 35164241 PMCID: PMC8838816 DOI: 10.3390/molecules27030976] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 01/25/2022] [Accepted: 01/28/2022] [Indexed: 12/18/2022]
Abstract
High resolution data from all-atom molecular simulations is used to parameterize a Martini 3 coarse-grained (CG) model of cellulose I allomorphs and cellulose type-II fibrils. In this case, elementary molecules are represented by four effective beads centred in the positions of O2, O3, C6, and O6 atoms in the D-glucose cellulose subunit. Non-bonded interactions between CG beads are tuned according to a low statistical criterion of structural deviation using the Martini 3 type of interactions and are capable of being indistinguishable for all studied cases. To maintain the crystalline structure of each single cellulose chain in the microfibrils, elastic potentials are employed to retain the ribbon-like structure in each chain. We find that our model is capable of describing different fibril-twist angles associated with each type of cellulose fibril in close agreement with atomistic simulation. Furthermore, our CG model poses a very small deviation from the native-like structure, making it appropriate to capture large conformational changes such as those that occur during the self-assembly process. We expect to provide a computational model suitable for several new applications such as cellulose self-assembly in different aqueous solutions and the thermal treatment of fibrils of great importance in bioindustrial applications.
Collapse
|
6
|
Cieplak M, Mioduszewski Ł, Chwastyk M. Contact-Based Analysis of Aggregation of Intrinsically Disordered Proteins. Methods Mol Biol 2022; 2340:105-120. [PMID: 35167072 DOI: 10.1007/978-1-0716-1546-1_6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We review the contact-based description of aggregation of intrinsically disordered proteins in coarse-grained and all-atom models. We consider polyglutamines and polyalanines at various concentrations of the peptides. We also study associations of two chains of α-synuclein and up to 20 chains of a 12-residue-long segment of protein tau. We demonstrate that the total number of two-chain association events (in an aggregate that comprises at least two chains) provides a useful measure of the propensity to aggregate. This measure is consistent, for instance, with the previously reported mass spectroscopy data. The distribution of the number of association events is given essentially by a power law as a function of the duration of these events. The corresponding exponent depends on the protein and the temperature but not on the concentration of the proteins.
Collapse
Affiliation(s)
- Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland.
| | | | - Mateusz Chwastyk
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| |
Collapse
|
7
|
Moreira RA, Baker JL, Guzman HV, Poma AB. Assessing the Stability of Biological Fibrils by Molecular-Scale Simulations. Methods Mol Biol 2022; 2340:357-378. [PMID: 35167082 DOI: 10.1007/978-1-0716-1546-1_16] [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] [Indexed: 06/14/2023]
Abstract
The nanomechanical characterization of several biological fibrils that are the result of protein aggregation via molecular dynamics simulation is nowadays feasible, and together with atomic force microscopy experiments has widened our understanding of the forces in the regime of pN-nN and system sizes of about hundreds of nanometers. Several methodologies have been developed to achieve this target, and they range from the atomistic representation via molecular force fields to coarse-grained strategies that provide comparable results with experiments in a systematic way. In this chapter, we discuss several methodologies for the calculation of mechanical parameters, such as the elastic constants of relevant biological systems. They are presented together with details about parameterization and current limitations. Then, we discuss some of the applications of such methodologies for the description of bacterial filament and β-amyloid systems. Finally, the latest lines of development are discussed.
Collapse
Affiliation(s)
- Rodrigo A Moreira
- Soft Matter and Biosystems, Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Joseph L Baker
- Department of Chemistry, The College of New Jersey, Ewing, NJ, USA
| | | | - Adolfo B Poma
- Soft Matter and Biosystems, Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland.
| |
Collapse
|
8
|
Poma AB, Thu TTM, Tri LTM, Nguyen HL, Li MS. Nanomechanical Stability of Aβ Tetramers and Fibril-like Structures: Molecular Dynamics Simulations. J Phys Chem B 2021; 125:7628-7637. [PMID: 34253022 PMCID: PMC8389904 DOI: 10.1021/acs.jpcb.1c02322] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Revised: 06/25/2021] [Indexed: 02/07/2023]
Abstract
Alzheimer's disease (AD) is a neurodegenerative disorder and one of the main causes of dementia. The disease is associated with amyloid beta (Aβ) peptide aggregation forming initial clusters and then fibril structure and plaques. Other neurodegenerative diseases such as type 2 diabetes, amyotrophic lateral sclerosis, and Parkinson's disease follow a similar mechanism. Therefore, inhibition of Aβ aggregation is considered an effective way to prevent AD. Recent experiments have provided evidence that oligomers are more toxic agents than mature fibrils, prompting researchers to investigate various factors that may influence their properties. One of these factors is nanomechanical stability, which plays an important role in the self-assembly of Aβ and possibly other proteins. This stability is also likely to be related to cell toxicity. In this work, we compare the mechanical stability of Aβ-tetramers and fibrillar structures using a structure-based coarse-grained (CG) approach and all-atom molecular dynamics simulation. Our results support the evidence for an increase in mechanical stability during the Aβ fibrillization process, which is consistent with in vitro AFM characterization of Aβ42 oligomers. Namely, using a CG model, we showed that the Young modulus of tetramers is lower than that of fibrils and, as follows from the experiment, is about 1 GPa. Hydrogen bonds are the dominant contribution to the detachment of one chain from the Aβ fibril fragment. They tend to be more organized along the pulling direction, whereas in the Aβ tetramers no preference is observed.
Collapse
Affiliation(s)
- Adolfo B. Poma
- Institute
of Fundamental Technological Research, Polish
Academy of Sciences, Pawińskiego 5B, 02-106 Warsaw, Poland
- International
Center for Research on Innovative Biobased Materials (ICRI-BioM)—International
Research Agenda, Lodz University of Technology, Żeromskiego 116, 90-924 Lodz, Poland
| | - Tran Thi Minh Thu
- Institute
for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh
Hiep Ward, District 12, Ho Chi Minh City, Vietnam
- Faculty
of Materials Science and Technology, Ho
Chi Minh City University of Science - VNUHCM, 227 Nguyen Van Cu Street, District 5, Ho Chi Minh City, Vietnam
- Vietnam
National University, Ho Chi Minh
City 700000, Vietnam
| | - Lam Tang Minh Tri
- Faculty
of Materials Science and Technology, Ho
Chi Minh City University of Science - VNUHCM, 227 Nguyen Van Cu Street, District 5, Ho Chi Minh City, Vietnam
- Vietnam
National University, Ho Chi Minh
City 700000, Vietnam
| | - Hoang Linh Nguyen
- Institute
for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh
Hiep Ward, District 12, Ho Chi Minh City, Vietnam
- Ho
Chi Minh City University of Technology (HCMUT), Ho Chi Minh City 700000, Vietnam
- Vietnam
National University, Ho Chi Minh
City 700000, Vietnam
| | - Mai Suan Li
- Institute
of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
9
|
Abstract
Cellulose is the most common biopolymer and widely used in our daily life. Due to its unique properties and biodegradability, it has been attracting increased attention in the recent years and various new applications of cellulose and its derivatives are constantly being found. The development of new materials with improved properties, however, is not always an easy task, and theoretical models and computer simulations can often help in this process. In this review, we give an overview of different coarse-grained models of cellulose and their applications to various systems. Various coarse-grained models with different mapping schemes are presented, which can efficiently simulate systems from the single cellulose fibril/crystal to the assembly of many fibrils/crystals. We also discuss relevant applications of these models with a focus on the mechanical properties, self-assembly, chiral nematic phases, conversion between cellulose allomorphs, composite materials and interactions with other molecules.
Collapse
|
10
|
Mioduszewski Ł, Cieplak M. Viscoelastic properties of wheat gluten in a molecular dynamics study. PLoS Comput Biol 2021; 17:e1008840. [PMID: 33760823 PMCID: PMC8021197 DOI: 10.1371/journal.pcbi.1008840] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2020] [Revised: 04/05/2021] [Accepted: 02/28/2021] [Indexed: 11/19/2022] Open
Abstract
Wheat (Triticum spp.) gluten consists mainly of intrinsincally disordered storage proteins (glutenins and gliadins) that can form megadalton-sized networks. These networks are responsible for the unique viscoelastic properties of wheat dough and affect the quality of bread. These properties have not yet been studied by molecular level simulations. Here, we use a newly developed α-C-based coarse-grained model to study ∼ 4000-residue systems. The corresponding time-dependent properties are studied through shear and axial deformations. We measure the response force to the deformation, the number of entanglements and cavities, the mobility of residues, the number of the inter-chain bonds, etc. Glutenins are shown to influence the mechanics of gluten much more than gliadins. Our simulations are consistent with the existing ideas about gluten elasticity and emphasize the role of entanglements and hydrogen bonding. We also demonstrate that the storage proteins in maize and rice lead to weaker elasticity which points to the unique properties of wheat gluten. During the breadmaking process, expanding gas bubbles cause the dough to increase volume. Gluten proteins act as an elastic scaffold in that process, allowing the wheat dough to grow more than other kinds of dough. Thus, explaining the unique viscoelastic properties of gluten at the molecular level may be of great interest to the baking industry. Assessing the structural properties of gluten is difficult because its proteins are disordered. We provide the first molecular dynamics model of gluten elasticity, that is able to distinguish gluten and proteins from different plants, like maize and rice. Our model shows the structural changes the gluten proteins undergo during their deformation, which mimics the mixing of dough during kneading. It also allows for a determination of the force required to extend gluten proteins, as during baking. The data confirms existing theories about gluten, but it also provides molecular-level information about the extraordinary elasticity of gluten.
Collapse
Affiliation(s)
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
- * E-mail:
| |
Collapse
|
11
|
Mahmood MI, Poma AB, Okazaki KI. Optimizing Gō-MARTINI Coarse-Grained Model for F-BAR Protein on Lipid Membrane. Front Mol Biosci 2021; 8:619381. [PMID: 33693028 PMCID: PMC7937874 DOI: 10.3389/fmolb.2021.619381] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 01/14/2021] [Indexed: 12/31/2022] Open
Abstract
Coarse-grained (CG) molecular dynamics (MD) simulations allow us to access much larger length and time scales than atomistic MD simulations, providing an attractive alternative to the conventional simulations. Based on the well-known MARTINI CG force field, the recently developed Gō-MARTINI model for proteins describes large-amplitude structural dynamics, which has not been possible with the commonly used elastic network model. Using the Gō-MARTINI model, we conduct MD simulations of the F-BAR Pacsin1 protein on lipid membrane. We observe that structural changes of the non-globular protein are largely dependent on the definition of the native contacts in the Gō model. To address this issue, we introduced a simple cutoff scheme and tuned the cutoff distance of the native contacts and the interaction strength of the Lennard-Jones potentials in the Gō-MARTINI model. With the optimized Gō-MARTINI model, we show that it reproduces structural fluctuations of the Pacsin1 dimer from atomistic simulations. We also show that two Pacsin1 dimers properly assemble through lateral interaction on the lipid membrane. Our work presents a first step towards describing membrane remodeling processes in the Gō-MARTINI CG framework by simulating a crucial step of protein assembly on the membrane.
Collapse
Affiliation(s)
- Md Iqbal Mahmood
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan
| | - Adolfo B Poma
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Kei-Ichi Okazaki
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan
| |
Collapse
|
12
|
Shivgan AT, Marzinek JK, Huber RG, Krah A, Henchman RH, Matsudaira P, Verma CS, Bond PJ. Extending the Martini Coarse-Grained Force Field to N-Glycans. J Chem Inf Model 2020; 60:3864-3883. [PMID: 32702979 DOI: 10.1021/acs.jcim.0c00495] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.
Collapse
Affiliation(s)
- Aishwary T Shivgan
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543.,Bioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Jan K Marzinek
- Bioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Roland G Huber
- Bioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Alexander Krah
- Bioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom.,Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| | - Paul Matsudaira
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543.,Centre for BioImaging Sciences, National University of Singapore, Singapore 117543
| | - Chandra S Verma
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543.,Bioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671.,School of Biological Sciences, Nanyang Technological University, 50 Nanyang Drive, Singapore 637551
| | - Peter J Bond
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543.,Bioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| |
Collapse
|
13
|
Mioduszewski Ł, Różycki B, Cieplak M. Pseudo-Improper-Dihedral Model for Intrinsically Disordered Proteins. J Chem Theory Comput 2020; 16:4726-4733. [PMID: 32436706 PMCID: PMC7588027 DOI: 10.1021/acs.jctc.0c00338] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
We present a new coarse-grained Cα-based protein model with a nonradial multibody pseudo-improper-dihedral potential that is transferable, time-independent, and suitable for molecular dynamics. It captures the nature of backbone and side-chain interactions between amino acid residues by adapting a simple improper dihedral term for a one-bead-per-residue model. It is parameterized for intrinsically disordered proteins and applicable to simulations of such proteins and their assemblies on millisecond time scales.
Collapse
Affiliation(s)
- Łukasz Mioduszewski
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Bartosz Różycki
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
14
|
Transient knots in intrinsically disordered proteins and neurodegeneration. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 174:79-103. [PMID: 32828471 DOI: 10.1016/bs.pmbts.2020.03.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
We provide a brief overview of the topological features found in structured proteins and of the dynamical processes that involve knots. We then discuss the knotted states that arise in the intrinsically disordered polyglutamine and α-synuclein. We argue that the existence of the knotted conformations stalls degradation by proteases and thus enhances aggregation. This mechanism works if the length of a peptide chain exceeds a threshold, as in the Huntington disease. We also study the cavities that form within the conformations of the disordered proteins. The volume of the cavities varies in time in a way that is different than that of the radius of gyration or the end-to-end distance. In addition, we study the traffic between the conformational basins and identify patterns associated with the deep and shallow knots. The results are obtained by molecular dynamics simulations that use coarse-grained and all-atom models (with and without the explicit solvent).
Collapse
|
15
|
Mioduszewski Ł, Cieplak M. Protein droplets in systems of disordered homopeptides and the amyloid glass phase. Phys Chem Chem Phys 2020; 22:15592-15599. [DOI: 10.1039/d0cp01635g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We study phase diagrams and liquid droplet formation for systems of polyglutamine and polyalanine chains in a coarse-grained model.
Collapse
Affiliation(s)
| | - Marek Cieplak
- Institute of Physics
- Polish Academy of Sciences
- 02-668 Warsaw
- Poland
| |
Collapse
|
16
|
Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications. Int J Mol Sci 2019; 20:ijms20153774. [PMID: 31375023 PMCID: PMC6696403 DOI: 10.3390/ijms20153774] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2019] [Revised: 07/28/2019] [Accepted: 07/30/2019] [Indexed: 12/23/2022] Open
Abstract
Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.
Collapse
|
17
|
Beltran-Villegas DJ, Intriago D, Kim KHC, Behabtu N, Londono JD, Jayaraman A. Coarse-grained molecular dynamics simulations of α-1,3-glucan. SOFT MATTER 2019; 15:4669-4681. [PMID: 31112203 DOI: 10.1039/c9sm00580c] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
In this paper we present a computational study of aggregation in aqueous solutions of α-1,3-glucan captured using a coarse-grained (CG) model that can be extended to other polysaccharides. This CG model captures atomistic geometry (i.e., relative placement of the hydrogen bonding donors and acceptors within the monomer) of the α-1,3-glucan monomer, the directional interactions due to the donor-acceptor hydrogen bonds, and their effect on aggregation of multiple α-1,3-glucan chains without the extensive computational resources needed for simulations with atomistic models. Using this CG model, we conduct molecular dynamics simulations to assess the effect of varying α-1,3-glucan chain length and hydrogen bond interaction strengths on the aggregation of multiple chains at finite concentrations in implicit solvent. We quantify the hydrogen bonding strength needed for multiple chains to aggregate, the distribution of inter- and intra-chain hydrogen bonds within the aggregate and in some cases, the shapes of the aggregate. We also explore the effect of substitution/silencing of some randomly selected or specific hydrogen bonding sites in the chain on the aggregation and aggregate structure. In the unmodified α-1,3-glucan solution, the inter-chain hydrogen bonds cause the chains to aggregate into sheets. Random silencing of hydrogen bonding donor sites only increases the hydrogen bond strength needed for aggregation but retains the same aggregate structure as the unmodified chains. Specific silencing of the hydrogen-bonding site on the C6 carbon leads to the chains aggregating into planar sheets that then fold over to form hollow cylinders at intermediate hydrogen bond strength - 4.7 to 5.3 kcal mol-1. These cylindrical aggregates assemble end-to-end to form larger aggregates at higher hydrogen bond strengths.
Collapse
|
18
|
Senapati S, Poma AB, Cieplak M, Filipek S, Park PSH. Differentiating between Inactive and Active States of Rhodopsin by Atomic Force Microscopy in Native Membranes. Anal Chem 2019; 91:7226-7235. [PMID: 31074606 DOI: 10.1021/acs.analchem.9b00546] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Membrane proteins, including G protein-coupled receptors (GPCRs), present a challenge in studying their structural properties under physiological conditions. Moreover, to better understand the activity of proteins requires examination of single molecule behaviors rather than ensemble averaged behaviors. Force-distance curve-based AFM (FD-AFM) was utilized to directly probe and localize the conformational states of a GPCR within the membrane at nanoscale resolution based on the mechanical properties of the receptor. FD-AFM was applied to rhodopsin, the light receptor and a prototypical GPCR, embedded in native rod outer segment disc membranes from photoreceptor cells of the retina in mice. Both FD-AFM and computational studies on coarse-grained models of rhodopsin revealed that the active state of the receptor has a higher Young's modulus compared to the inactive state of the receptor. Thus, the inactive and active states of rhodopsin could be differentiated based on the stiffness of the receptor. Differentiating the states based on the Young's modulus allowed for the mapping of the different states within the membrane. Quantifying the active states present in the membrane containing the constitutively active G90D rhodopsin mutant or apoprotein opsin revealed that most receptors adopt an active state. Traditionally, constitutive activity of GPCRs has been described in terms of two-state models where the receptor can achieve only a single active state. FD-AFM data are inconsistent with a two-state model but instead require models that incorporate multiple active states.
Collapse
Affiliation(s)
- Subhadip Senapati
- Department of Ophthalmology and Visual Sciences , Case Western Reserve University , Cleveland , Ohio 44106 , United States
| | - Adolfo B Poma
- Institute of Fundamental Technological Research , Polish Academy of Sciences , Pawińskiego 5B , 02-106 Warsaw , Poland.,Institute of Physics , Polish Academy of Sciences , Aleja Lotników 32/46 , 02-668 Warsaw , Poland
| | - Marek Cieplak
- Institute of Physics , Polish Academy of Sciences , Aleja Lotników 32/46 , 02-668 Warsaw , Poland
| | - Sławomir Filipek
- Faculty of Chemistry, Biological and Chemical Research Centre , University of Warsaw , 02-093 Warsaw , Poland
| | - Paul S H Park
- Department of Ophthalmology and Visual Sciences , Case Western Reserve University , Cleveland , Ohio 44106 , United States
| |
Collapse
|
19
|
Poma AB, Li MS, Theodorakis PE. Generalization of the elastic network model for the study of large conformational changes in biomolecules. Phys Chem Chem Phys 2019; 20:17020-17028. [PMID: 29904772 DOI: 10.1039/c8cp03086c] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Gō-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Gō-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for α-helix and β-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.
Collapse
Affiliation(s)
- Adolfo B Poma
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland.
| | | | | |
Collapse
|
20
|
Poma AB, Guzman HV, Li MS, Theodorakis PE. Mechanical and thermodynamic properties of Aβ 42, Aβ 40, and α-synuclein fibrils: a coarse-grained method to complement experimental studies. BEILSTEIN JOURNAL OF NANOTECHNOLOGY 2019; 10:500-513. [PMID: 30873322 PMCID: PMC6404408 DOI: 10.3762/bjnano.10.51] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Accepted: 02/08/2019] [Indexed: 05/07/2023]
Abstract
We perform molecular dynamics simulation on several relevant biological fibrils associated with neurodegenerative diseases such as Aβ40, Aβ42, and α-synuclein systems to obtain a molecular understanding and interpretation of nanomechanical characterization experiments. The computational method is versatile and addresses a new subarea within the mechanical characterization of heterogeneous soft materials. We investigate both the elastic and thermodynamic properties of the biological fibrils in order to substantiate experimental nanomechanical characterization techniques that are quickly developing and reaching dynamic imaging with video rate capabilities. The computational method qualitatively reproduces results of experiments with biological fibrils, validating its use in extrapolation to macroscopic material properties. Our computational techniques can be used for the co-design of new experiments aiming to unveil nanomechanical properties of biological fibrils from a point of view of molecular understanding. Our approach allows a comparison of diverse elastic properties based on different deformations , i.e., tensile (Y L), shear (S), and indentation (Y T) deformation. From our analysis, we find a significant elastic anisotropy between axial and transverse directions (i.e., Y T > Y L) for all systems. Interestingly, our results indicate a higher mechanostability of Aβ42 fibrils compared to Aβ40, suggesting a significant correlation between mechanical stability and aggregation propensity (rate) in amyloid systems. That is, the higher the mechanical stability the faster the fibril formation. Finally, we find that α-synuclein fibrils are thermally less stable than β-amyloid fibrils. We anticipate that our molecular-level analysis of the mechanical response under different deformation conditions for the range of fibrils considered here will provide significant insights for the experimental observations.
Collapse
Affiliation(s)
- Adolfo B Poma
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Pawińskiego 5B, 02-106 Warsaw, Poland
| | - Horacio V Guzman
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
- Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | | |
Collapse
|
21
|
Mioduszewski Ł, Cieplak M. Disordered peptide chains in an α-C-based coarse-grained model. Phys Chem Chem Phys 2018; 20:19057-19070. [PMID: 29972174 DOI: 10.1039/c8cp03309a] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We construct a one-bead-per-residue coarse-grained dynamical model to describe intrinsically disordered proteins at significantly longer timescales than in the all-atom models. In this model, inter-residue contacts form and disappear during the course of the time evolution. The contacts may arise between the sidechains, the backbones or the sidechains and backbones of the interacting residues. The model yields results that are consistent with many all-atom and experimental data on these systems. We demonstrate that the geometrical properties of various homopeptides differ substantially in this model. In particular, the average radius of gyration scales with the sequence length in a residue-dependent manner.
Collapse
Affiliation(s)
- Łukasz Mioduszewski
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland.
| | | |
Collapse
|
22
|
Zhao Y, Chwastyk M, Cieplak M. Structural entanglements in protein complexes. J Chem Phys 2018; 146:225102. [PMID: 29166058 DOI: 10.1063/1.4985221] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We consider multi-chain protein native structures and propose a criterion that determines whether two chains in the system are entangled or not. The criterion is based on the behavior observed by pulling at both termini of each chain simultaneously in the two chains. We have identified about 900 entangled systems in the Protein Data Bank and provided a more detailed analysis for several of them. We argue that entanglement enhances the thermodynamic stability of the system but it may have other functions: burying the hydrophobic residues at the interface and increasing the DNA or RNA binding area. We also study the folding and stretching properties of the knotted dimeric proteins MJ0366, YibK, and bacteriophytochrome. These proteins have been studied theoretically in their monomeric versions so far. The dimers are seen to separate on stretching through the tensile mechanism and the characteristic unraveling force depends on the pulling direction.
Collapse
Affiliation(s)
- Yani Zhao
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, PL-02668 Warsaw, Poland
| | - Mateusz Chwastyk
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, PL-02668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, PL-02668 Warsaw, Poland
| |
Collapse
|
23
|
Poma AB, Chwastyk M, Cieplak M. Elastic moduli of biological fibers in a coarse-grained model: crystalline cellulose and β-amyloids. Phys Chem Chem Phys 2018; 19:28195-28206. [PMID: 29022971 DOI: 10.1039/c7cp05269c] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
We study the mechanical response of cellulose and β-amyloid microfibrils to three types of deformation: tensile, indentational, and shear. The cellulose microfibrils correspond to the allomorphs Iα or Iβ whereas the β-amyloid microfibrils correspond to the polymorphs of either two- or three-fold symmetry. This response can be characterized by three elastic moduli, namely, YL, YT, and S. We use a structure-based coarse-grained model to analyze the deformations in a unified manner. We find that each of the moduli is almost the same for the two allomorphs of cellulose but YL is about 20 times larger than YT (140 GPa vs. 7 GPa), indicating the existence of significant anisotropy. For cellulose we note that the anisotropy results from the involvement of covalent bonds in stretching. For β-amyloid, the sense of anisotropy is opposite to that of cellulose. In the three-fold symmetry case, YL is about half of YT (3 vs. 7) whereas for two-fold symmetry the anisotropy is much larger (1.6 vs. 21 GPa). The S modulus is derived to be 1.2 GPa for three-fold symmetry and one half of it for the other symmetry and 3.0 GPa for cellulose. The values of the moduli reflect deformations in the hydrogen-bond network. Unlike in our theoretical approach, no experiment can measure all three elastic moduli with the same apparatus. However, our theoretical results are consistent with various measured values: typical YL for cellulose Iβ ranges from 133 to 155 GPa, YT from 2 to 25 GPa, and S from 1.8 to 3.8 GPa. For β-amyloid, the experimental values of S and YT are about 0.3 GPa and 3.3 GPa respectively, while the value of YL has not been reported.
Collapse
Affiliation(s)
- Adolfo B Poma
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, PL-02668 Warsaw, Poland.
| | | | | |
Collapse
|
24
|
Wołek K, Cieplak M. Self-assembly of model proteins into virus capsids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:474003. [PMID: 29027904 PMCID: PMC7104874 DOI: 10.1088/1361-648x/aa9351] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 09/29/2017] [Accepted: 10/13/2017] [Indexed: 06/07/2023]
Abstract
We consider self-assembly of proteins into a virus capsid by the methods of molecular dynamics. The capsid corresponds either to SPMV or CCMV and is studied with and without the RNA molecule inside. The proteins are flexible and described by the structure-based coarse-grained model augmented by electrostatic interactions. Previous studies of the capsid self-assembly involved solid objects of a supramolecular scale, e.g. corresponding to capsomeres, with engineered couplings and stochastic movements. In our approach, a single capsid is dissociated by an application of a high temperature for a variable period and then the system is cooled down to allow for self-assembly. The restoration of the capsid proceeds to various extent, depending on the nature of the dissociated state, but is rarely complete because some proteins depart too far unless the process takes place in a confined space.
Collapse
Affiliation(s)
- Karol Wołek
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
25
|
Lubecka EA, Liwo A. A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. II. Backbone-local potentials of coarse-grained O1→4-bonded polyglucose chains. J Chem Phys 2017; 147:115101. [DOI: 10.1063/1.4994130] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Emilia A. Lubecka
- Institute of Informatics, University of Gdańsk, Wita Stwosza 57, 80-308 Gdańsk, Poland
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
26
|
Różycki B, Cieplak M. Stiffness of the C-terminal disordered linker affects the geometry of the active site in endoglucanase Cel8A. MOLECULAR BIOSYSTEMS 2017; 12:3589-3599. [PMID: 27714009 DOI: 10.1039/c6mb00606j] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Cellulosomes are complex multi-enzyme machineries which efficiently degrade plant cell-wall polysaccharides. The multiple domains of the cellulosome proteins are often tethered together by intrinsically disordered regions. The properties and functions of these disordered linkers are not well understood. In this work, we study endoglucanase Cel8A, which is a relevant enzymatic component of the cellulosomes of Clostridium thermocellum. We use both all-atom and coarse-grained simulations to investigate how the conformations of the catalytic domain of Cel8A are affected by the disordered linker at its C terminus. We find that when the endoglucanase is bound to its substrate, the effective stiffness of the linker can influence the distances between groups of amino-acid residues throughout the entire enzymatic domain. In particular, variations in the linker stiffness can lead to small changes in the geometry of the active-site cleft. We suggest that such geometrical changes may have an effect on the catalytic activity of the enzyme.
Collapse
Affiliation(s)
- Bartosz Różycki
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland.
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland.
| |
Collapse
|
27
|
Zhao Y, Cieplak M. Structural Changes in Barley Protein LTP1 Isoforms at Air-Water Interfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:4769-4780. [PMID: 28457129 DOI: 10.1021/acs.langmuir.7b00791] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We use a coarse-grained model to study the conformational changes in two barley proteins, LTP1 and its ligand adduct isoform LTP1b, that result from their adsorption to the air-water interface. The model introduces the interface through hydropathy indices. We justify the model by all-atom simulations. The choice of the proteins is motivated by making attempts to understand formation and stability of foam in beer. We demonstrate that both proteins flatten out at the interface and can make a continuous stabilizing and denser film. We show that the degree of the flattening depends on the protein (the layers of LTP1b should be denser than those of LTP1) and on the presence of glycation. It also depends on the number (≤4) of the disulfide bonds in the proteins. The geometry of the proteins is sensitive to the specificity of the absent bonds. We provide estimates of the volume of cavities of the proteins when away from the interface.
Collapse
Affiliation(s)
- Yani Zhao
- Institute of Physics, Polish Academy of Sciences , Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences , Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
28
|
Poma AB, Cieplak M, Theodorakis PE. Combining the MARTINI and Structure-Based Coarse-Grained Approaches for the Molecular Dynamics Studies of Conformational Transitions in Proteins. J Chem Theory Comput 2017; 13:1366-1374. [PMID: 28195464 DOI: 10.1021/acs.jctc.6b00986] [Citation(s) in RCA: 106] [Impact Index Per Article: 15.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The application of coarse-grained (CG) models in biology is essential to access large length and time scales required for the description of many biological processes. The ELNEDIN protein model is based on the well-known MARTINI CG force-field and incorporates additionally harmonic bonds of a certain spring constant within a defined cutoff distance between pairs of residues, in order to preserve the native structure of the protein. In this case, the use of unbreakable harmonic bonds hinders the study of unfolding and folding processes. To overcome this barrier we have replaced the harmonic bonds with Lennard-Jones interactions based on the contact map of the native protein structure as is done in Go̅-like models. This model exhibits very good agreement with all-atom simulations and the ELNEDIN. Moreover, it can capture the structural motion linked to particular catalytic activity in the Man5B protein, in agreement with all-atom simulations. In addition, our model is based on the van der Waals radii, instead of a cutoff distance, which results in a smaller contact map. In conclusion, we anticipate that our model will provide further possibilities for studying biological systems based on the MARTINI CG force-field by using advanced-sampling methods, such as parallel tempering and metadynamics.
Collapse
Affiliation(s)
- Adolfo B Poma
- Institute of Physics, Polish Academy of Sciences , Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences , Al. Lotników 32/46, 02-668 Warsaw, Poland
| | | |
Collapse
|
29
|
Topological transformations in proteins: effects of heating and proximity of an interface. Sci Rep 2017; 7:39851. [PMID: 28051124 PMCID: PMC5209716 DOI: 10.1038/srep39851] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2016] [Accepted: 11/28/2016] [Indexed: 01/04/2023] Open
Abstract
Using a structure-based coarse-grained model of proteins, we study the mechanism of unfolding of knotted proteins through heating. We find that the dominant mechanisms of unfolding depend on the temperature applied and are generally distinct from those identified for folding at its optimal temperature. In particular, for shallowly knotted proteins, folding usually involves formation of two loops whereas unfolding through high-temperature heating is dominated by untying of single loops. Untying the knots is found to generally precede unfolding unless the protein is deeply knotted and the heating temperature exceeds a threshold value. We then use a phenomenological model of the air-water interface to show that such an interface can untie shallow knots, but it can also make knots in proteins that are natively unknotted.
Collapse
|
30
|
Wołek K, Cieplak M. Criteria for folding in structure-based models of proteins. J Chem Phys 2016; 144:185102. [DOI: 10.1063/1.4948783] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Karol Wołek
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|