1
|
Kharche S, Yadav M, Hande V, Prakash S, Sengupta D. Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model. J Chem Inf Model 2024; 64:837-850. [PMID: 38291973 DOI: 10.1021/acs.jcim.3c00802] [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: 02/01/2024]
Abstract
The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.
Collapse
Affiliation(s)
- Shalmali Kharche
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Manjul Yadav
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Vrushali Hande
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Shikha Prakash
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Durba Sengupta
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| |
Collapse
|
2
|
Singh AP, Kumar R, Gupta D. Structural insights into the mechanism of human methyltransferase hPRMT4. J Biomol Struct Dyn 2022; 40:10821-10834. [PMID: 34308797 DOI: 10.1080/07391102.2021.1950567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
Human PRMT4, also known as CARM1, is a type I arginine methyltransferase protein that catalyse the formation of asymmetrical dimethyl arginine product. Structural studies done to date on PRMT4 have shown that the N-terminal region, Rossmann fold and dimerization arm play an important role in PRMT4 activity. Elucidating the functions of these regions in catalysis remains to be explored. Studies have shown the existence of communication pathways in PRMT4, which need further elucidation. The molecular dynamics (MD) simulations performed in this study show differences in different monomeric and dimeric forms of hPRMT4, revealing the role of the N-terminal region, Rossmann fold and dimerization arm. The study shows the conformational changes that occur during dimerization and SAM binding. Our cross-correlation analysis showed a correlation between these regions. Further, we performed PSN and network analysis to establish the existence of communication networks and an allosteric pathway. This study shows the use of MD simulations and network analysis to explore the aspects of PRMT4 dimerization, SAM binding and demonstrates the existence of an allosteric network. These findings shed novel insights into the conformational changes associated with hPRMT4, the mechanism of its dimerization, SAM binding and clues for better inhibitor designs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Amar Pratap Singh
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, New Delhi, India
| | - Rakesh Kumar
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, New Delhi, India
| | - Dinesh Gupta
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, New Delhi, India
| |
Collapse
|
3
|
Sharma P, Desikan R, Ayappa KG. Evaluating Coarse-Grained MARTINI Force-Fields for Capturing the Ripple Phase of Lipid Membranes. J Phys Chem B 2021; 125:6587-6599. [PMID: 34081861 DOI: 10.1021/acs.jpcb.1c03277] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (Pβ') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolutions. Here, with an aggregate ∼2.5 μs all-atom and ∼122 μs CGMD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the Pβ' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (>2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated Pβ' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the Pβ' phase, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically trapped structures caused by finite size effects rather than being representative of true Pβ' phases. We showed that a MARTINI FF variant that could capture the tilted Lβ' gel phase, a prerequisite for stabilizing the Pβ' phase, was unable to capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific reparametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.
Collapse
Affiliation(s)
- Pradyumn Sharma
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru 560012, India
| | - Rajat Desikan
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru 560012, India
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru 560012, India.,Centre for BioSystems Science and Engineering, Indian Institute of Science, Bengaluru 560012, India
| |
Collapse
|
4
|
Desikan R, Behera A, Maiti PK, Ayappa KG. Using multiscale molecular dynamics simulations to obtain insights into pore forming toxin mechanisms. Methods Enzymol 2021; 649:461-502. [PMID: 33712196 DOI: 10.1016/bs.mie.2021.01.021] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Pore forming toxins (PFTs) are virulent proteins released by several species, including many strains of bacteria, to attack and kill host cells. In this article, we focus on the utility of molecular dynamics (MD) simulations and the molecular insights gleaned from these techniques on the pore forming pathways of PFTs. In addition to all-atom simulations which are widely used, coarse-grained MARTINI models and structure-based models have also been used to study PFTs. Here, the emphasis is on methods and techniques involved while setting up, monitoring, and evaluating properties from MD simulations of PFTs in a membrane environment. We draw from several case studies to illustrate how MD simulations have provided molecular insights into protein-protein and protein-lipid interactions, lipid dynamics, conformational transitions and structures of both the oligomeric intermediates and assembled pore structures.
Collapse
Affiliation(s)
- Rajat Desikan
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India
| | - Amit Behera
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India
| | - Prabal K Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bengaluru, India
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India; Centre for Biosystems Science and Engineering, Indian Institute of Science, Bengaluru, India.
| |
Collapse
|
5
|
Pore-forming proteins: From defense factors to endogenous executors of cell death. Chem Phys Lipids 2020; 234:105026. [PMID: 33309552 DOI: 10.1016/j.chemphyslip.2020.105026] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Revised: 12/07/2020] [Accepted: 12/07/2020] [Indexed: 02/06/2023]
Abstract
Pore-forming proteins (PFPs) and small antimicrobial peptides (AMPs) represent a large family of molecules with the common ability to punch holes in cell membranes to alter their permeability. They play a fundamental role as infectious bacteria's defensive tools against host's immune system and as executors of endogenous machineries of regulated cell death in eukaryotic cells. Despite being highly divergent in primary sequence and 3D structure, specific folds of pore-forming domains have been conserved. In fact, pore formation is considered an ancient mechanism that takes place through a general multistep process involving: membrane partitioning and insertion, oligomerization and pore formation. However, different PFPs and AMPs assemble and form pores following different mechanisms that could end up either in the formation of protein-lined or protein-lipid pores. In this review, we analyze the current findings in the mechanism of action of different PFPs and AMPs that support a wide role of membrane pore formation in nature. We also provide the newest insights into the development of state-of-art techniques that have facilitated the characterization of membrane pores. To understand the physiological role of these peptides/proteins or develop clinical applications, it is essential to uncover the molecular mechanism of how they perforate membranes.
Collapse
|
6
|
|
7
|
Desikan R, Maiti PK, Ayappa KG. Predicting interfacial hot-spot residues that stabilize protein-protein interfaces in oligomeric membrane-toxin pores through hydrogen bonds and salt bridges. J Biomol Struct Dyn 2020; 39:20-34. [DOI: 10.1080/07391102.2020.1711806] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Rajat Desikan
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, India
| | - Prabal K. Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore, India
| | - K. Ganapathy Ayappa
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, India
- Centre for Biosystems Science and Engineering, Indian Institute of Science, Bangalore, India
| |
Collapse
|
8
|
Ionic transport through a protein nanopore: a Coarse-Grained Molecular Dynamics Study. Sci Rep 2019; 9:15740. [PMID: 31673049 PMCID: PMC6823379 DOI: 10.1038/s41598-019-51942-y] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2018] [Accepted: 07/17/2019] [Indexed: 11/09/2022] Open
Abstract
The MARTINI coarse-grained (CG) force field is used to test the ability of CG models to simulate ionic transport through protein nanopores. The ionic conductivity of CG ions in solution was computed and compared with experimental results. Next, we studied the electrostatic behavior of a solvated CG lipid bilayer in salt solution under an external electric field. We showed this approach correctly describes the experimental conditions under a potential bias. Finally, we performed CG molecular dynamics simulations of the ionic transport through a protein nanopore (α-hemolysin) inserted in a lipid bilayer, under different electric fields, for 2-3 microseconds. The resulting I - V curve is qualitatively consistent with experiments, although the computed current is one order of magnitude smaller. Current saturation was observed for potential biases over ±350 mV. We also discuss the time to reach a stationary regime and the role of the protein flexibility in our CG simulations.
Collapse
|
9
|
Barrera EE, Machado MR, Pantano S. Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics. J Chem Theory Comput 2019; 15:5674-5688. [PMID: 31433946 DOI: 10.1021/acs.jctc.9b00435] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.
Collapse
Affiliation(s)
- Exequiel E Barrera
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Matías R Machado
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Sergio Pantano
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay.,Shanghai Institute for Advanced Immunochemical Studies , ShanghaiTech University , Shanghai 201210 , China
| |
Collapse
|
10
|
Abstract
Pneumolysin (PLY), a major virulence factor of Streptococcus pneumoniae, perforates cholesterol-rich lipid membranes. PLY protomers oligomerize as rings on the membrane and then undergo a structural transition that triggers the formation of membrane pores. Structures of PLY rings in prepore and pore conformations define the beginning and end of this transition, but the detailed mechanism of pore formation remains unclear. With atomistic and coarse-grained molecular dynamics simulations, we resolve key steps during PLY pore formation. Our simulations confirm critical PLY membrane-binding sites identified previously by mutagenesis. The transmembrane β-hairpins of the PLY pore conformation are stable only for oligomers, forming a curtain-like membrane-spanning β-sheet. Its hydrophilic inner face draws water into the protein-lipid interface, forcing lipids to recede. For PLY rings, this zone of lipid clearance expands into a cylindrical membrane pore. The lipid plug caught inside the PLY ring can escape by lipid efflux via the lower leaflet. If this path is too slow or blocked, the pore opens by membrane buckling, driven by the line tension acting on the detached rim of the lipid plug. Interestingly, PLY rings are just wide enough for the plug to buckle spontaneously in mammalian membranes. In a survey of electron cryo-microscopy (cryo-EM) and atomic force microscopy images, we identify key intermediates along both the efflux and buckling pathways to pore formation, as seen in the simulations.
Collapse
|
11
|
Ray S, Holden S, Martin LL, Panwar AS. Mechanistic insight into the early stages of amyloid formation using an anuran peptide. Pept Sci (Hoboken) 2019. [DOI: 10.1002/pep2.24120] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Sourav Ray
- IITB‐Monash Research AcademyIndian Institute of Technology Bombay Powai Mumbai India
- School of ChemistryMonash University Clayton Victoria Australia
- Department of Metallurgical Engineering and Materials ScienceIndian Institute of Technology Bombay Powai Mumbai India
| | | | | | - Ajay Singh Panwar
- Department of Metallurgical Engineering and Materials ScienceIndian Institute of Technology Bombay Powai Mumbai India
| |
Collapse
|
12
|
Desikan R, Maiti PK, Ayappa KG. Assessing the Structure and Stability of Transmembrane Oligomeric Intermediates of an α-Helical Toxin. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11496-11510. [PMID: 28930630 DOI: 10.1021/acs.langmuir.7b02277] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Protein membrane interactions play an important role in our understanding of diverse phenomena ranging from membrane-assisted protein aggregation to oligomerization and folding. Pore-forming toxins (PFTs) are the primary vehicle for infection by several strains of bacteria. These proteins which are expressed in a water-soluble form (monomers) bind to the target membrane and conformationally transform (protomers) and self-assemble to form a multimer transmembrane pore complex through a process of oligomerization. On the basis of the structure of the transmembrane domains, PFTs are broadly classified into β or α toxins. In contrast to β-PFTs, the paucity of available crystal structures coupled with the amphipathic nature of the transmembrane domains has hindered our understanding of α-PFT pore formation. In this article, we use molecular dynamics (MD) simulations to examine the process of pore formation of the bacterial α-PFT, cytolysin A from Escherichia coli (ClyA) in lipid bilayer membranes. Using atomistic MD simulations ranging from 50 to 500 ns, we show that transmembrane oligomeric intermediates or "arcs" form stable proteolipidic complexes consisting of protein arcs with toroidal lipids lining the free edges. By creating initial conditions where the lipids are contained within the arcs, we study the dynamics of spontaneous lipid evacuation and toroidal edge formation. This process occurs on the time scale of tens of nanoseconds, suggesting that once protomers oligomerize, transmembrane arcs are rapidly stabilized to form functional water channels capable of leakage. Using umbrella sampling with a coarse-grained molecular model, we obtain the free energy of insertion of a single protomer into the membrane. A single inserted protomer has a stabilization free energy of -52.9 ± 1.2 kJ/mol and forms a stable transmembrane water channel capable of leakage. Our simulations reveal that arcs are stable and viable intermediates that can occur during the pore-formation pathway for ClyA.
Collapse
Affiliation(s)
- Rajat Desikan
- Department of Chemical Engineering, ‡Centre for Condensed Matter Theory, Department of Physics, and §Centre for Biosystems Science and Engineering, Indian Institute of Science , Bengaluru, India 560012
| | - Prabal K Maiti
- Department of Chemical Engineering, ‡Centre for Condensed Matter Theory, Department of Physics, and §Centre for Biosystems Science and Engineering, Indian Institute of Science , Bengaluru, India 560012
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, ‡Centre for Condensed Matter Theory, Department of Physics, and §Centre for Biosystems Science and Engineering, Indian Institute of Science , Bengaluru, India 560012
| |
Collapse
|