1
|
Bajpai S, Petkov BK, Tong M, Abreu CRA, Nair NN, Tuckerman ME. An interoperable implementation of collective-variable based enhanced sampling methods in extended phase space within the OpenMM package. J Comput Chem 2023; 44:2166-2183. [PMID: 37464902 DOI: 10.1002/jcc.27182] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 05/30/2023] [Accepted: 06/06/2023] [Indexed: 07/20/2023]
Abstract
Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. In this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.
Collapse
Affiliation(s)
- Shitanshu Bajpai
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India
| | - Brian K Petkov
- Department of Chemistry, New York University (NYU), New York, New York, USA
| | - Muchen Tong
- Department of Chemistry, New York University (NYU), New York, New York, USA
| | - Charlles R A Abreu
- Chemical Engineering Department, Escola de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York, USA
- Courant Institute of Mathematical Sciences, New York University (NYU), New York, New York, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- Simons Center for Computational Physical Chemistry, New York University, New York, New York, USA
| |
Collapse
|
2
|
Spanò A, Fanton L, Pizzolato D, Moi J, Vinci F, Pesce A, Dongmo Foumthuim CJ, Giacometti A, Simeoni M. Rinmaker: a fast, versatile and reliable tool to determine residue interaction networks in proteins. BMC Bioinformatics 2023; 24:336. [PMID: 37697267 PMCID: PMC10496328 DOI: 10.1186/s12859-023-05466-y] [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/02/2023] [Accepted: 09/04/2023] [Indexed: 09/13/2023] Open
Abstract
BACKGROUND Residue Interaction Networks (RINs) map the crystallographic description of a protein into a graph, where amino acids are represented as nodes and non-covalent bonds as edges. Determination and visualization of a protein as a RIN provides insights on the topological properties (and hence their related biological functions) of large proteins without dealing with the full complexity of the three-dimensional description, and hence it represents an invaluable tool of modern bioinformatics. RESULTS We present RINmaker, a fast, flexible, and powerful tool for determining and visualizing RINs that include all standard non-covalent interactions. RINmaker is offered as a cross-platform and open source software that can be used either as a command-line tool or through a web application or a web API service. We benchmark its efficiency against the main alternatives and provide explicit tests to show its performance and its correctness. CONCLUSIONS RINmaker is designed to be fully customizable, from a simple and handy support for experimental research to a sophisticated computational tool that can be embedded into a large computational pipeline. Hence, it paves the way to bridge the gap between data-driven/machine learning approaches and numerical simulations of simple, physically motivated, models.
Collapse
Affiliation(s)
- Alvise Spanò
- Department of Environmental Science, Computer Science and Statistics, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Lorenzo Fanton
- Department of Environmental Science, Computer Science and Statistics, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Davide Pizzolato
- Department of Environmental Science, Computer Science and Statistics, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Jacopo Moi
- Department of Molecular Science and Nanosystems, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Francesco Vinci
- Department of Environmental Science, Computer Science and Statistics, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Alberto Pesce
- Department of Environmental Science, Computer Science and Statistics, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Cedrix J Dongmo Foumthuim
- Department of Molecular Science and Nanosystems, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
| | - Achille Giacometti
- Department of Molecular Science and Nanosystems, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy
- European Centre for Living Technology (ECLT), Dorsoduro 3246, 30123, Venice, Italy
| | - Marta Simeoni
- Department of Environmental Science, Computer Science and Statistics, University Ca' Foscari of Venice, Via Torino 155, 30172, Venice, Italy.
- European Centre for Living Technology (ECLT), Dorsoduro 3246, 30123, Venice, Italy.
| |
Collapse
|
3
|
Kapakayala AB, Nair NN. Boosting the conformational sampling by combining replica exchange with solute tempering and well-sliced metadynamics. J Comput Chem 2021; 42:2233-2240. [PMID: 34585768 DOI: 10.1002/jcc.26752] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2021] [Revised: 08/30/2021] [Accepted: 09/12/2021] [Indexed: 01/22/2023]
Abstract
Methods that combine collective variable (CV) based enhanced sampling and global tempering approaches are used in speeding-up the conformational sampling and free energy calculation of large and soft systems with a plethora of energy minima. In this paper, a new method of this kind is proposed in which the well-sliced metadynamics approach (WSMTD) is united with replica exchange with solute tempering (REST2) method. WSMTD employs a divide-and-conquer strategy wherein high-dimensional slices of a free energy surface are independently sampled and combined. The method enables one to accomplish a controlled exploration of the CV-space with a restraining bias as in umbrella sampling, and enhance-sampling of one or more orthogonal CVs using a metadynamics like bias. The new hybrid method proposed here enables boosting the sampling of more slow degrees of freedom in WSMTD simulations, without the need to specify associated CVs, through a replica exchange scheme within the framework of REST2. The high-dimensional slices of the probability distributions of CVs computed from the united WSMTD and REST2 simulations are subsequently combined using the weighted histogram analysis method to obtain the free energy surface. We show that the new method proposed here is accurate, improves the conformational sampling, and achieves quick convergence in free energy estimates. We demonstrate this by computing the conformational free energy landscapes of solvated alanine tripeptide and Trp-cage mini protein in explicit water.
Collapse
Affiliation(s)
- Anji Babu Kapakayala
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India.,School of Pharmacy and Biomedical Sciences, Curtin University, Perth, Australia
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| |
Collapse
|
4
|
Meuzelaar H, Panman MR, van Dijk CN, Woutersen S. Folding of a Zinc-Finger ββα-Motif Investigated Using Two-Dimensional and Time-Resolved Vibrational Spectroscopy. J Phys Chem B 2016; 120:11151-11158. [PMID: 27723346 DOI: 10.1021/acs.jpcb.6b08883] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Small proteins provide good model systems for studying the fundamental forces that control protein folding. Here, we investigate the folding dynamics of the 28-residue zinc-finger mutant FSD-1, which is designed to form a metal-independent folded ββα-motif, and which provides a testing ground for proteins containing a mixed α/β fold. Although the folding of FSD-1 has been actively studied, the folding mechanism remains largely unclear. In particular, it is unclear in what stage of folding the α-helix is formed. To address this issue we investigate the folding mechanism of FSD-1 using a combination of temperature-dependent UV circular dichroism (UV-CD), Fourier transform infrared (FTIR) spectroscopy, two-dimensional infrared (2D-IR) spectroscopy, and temperature-jump (T-jump) transient-IR spectroscopy. Our UV-CD and FTIR data show different thermal melting transitions, indicating multistate folding behavior. Temperature-dependent 2D-IR spectra indicate that the α-helix is the most stable structural element of FSD-1. To investigate the folding/unfolding re-equilibration dynamics of FSD-1, the conformational changes induced by a nanosecond T-jump are probed with transient-IR and transient dispersed-pump-probe (DPP) IR spectroscopy. We observe biexponential T-jump relaxation kinetics (with time constants of 80 ± 13 ns and 1300 ± 100 ns at 322 K), confirming that the folding involves an intermediate state. The IR and dispersed-pump-probe IR spectra associated with the two kinetic components suggest that the folding of FSD-1 involves early formation of the α-helix, followed by the formation of the β-hairpin and hydrophobic contacts.
Collapse
Affiliation(s)
- Heleen Meuzelaar
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Matthijs R Panman
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Chris N van Dijk
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Sander Woutersen
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| |
Collapse
|
5
|
Gupta M, Nayar D, Chakravarty C, Bandyopadhyay S. Comparison of hydration behavior and conformational preferences of the Trp-cage mini-protein in different rigid-body water models. Phys Chem Chem Phys 2016; 18:32796-32813. [DOI: 10.1039/c6cp04634g] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Trp-cage unfolds at different temperatures in different water models revealing the sensitivity of conformational order metrics to the choice of water models.
Collapse
Affiliation(s)
- Madhulika Gupta
- Department of Chemistry
- Indian Institute of Technology-Delhi
- New Delhi 110016
- India
| | - Divya Nayar
- Department of Chemistry
- Indian Institute of Technology-Delhi
- New Delhi 110016
- India
| | | | - Sanjoy Bandyopadhyay
- Molecular Modeling Laboratory
- Department of Chemistry
- Indian Institute of Technology
- Kharagpur 721302
- India
| |
Collapse
|
6
|
Vymětal J, Bathula SR, Černý J, Chaloupková R, Žídek L, Sklenář V, Vondrášek J. Retro operation on the Trp-cage miniprotein sequence produces an unstructured molecule capable of folding similar to the original only upon 2,2,2-trifluoroethanol addition. Protein Eng Des Sel 2014; 27:463-72. [DOI: 10.1093/protein/gzu046] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
7
|
Qiu W, Skafidas E. Detection of protein conformational changes with multilayer graphene nanopore sensors. ACS APPLIED MATERIALS & INTERFACES 2014; 6:16777-81. [PMID: 25185959 DOI: 10.1021/am5040279] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Detecting conformational change in protein or peptide is imperative in understanding their dynamic function and diagnosing diseases. Existing techniques either rely on ensemble average that lacks the necessary sensitivity or require florescence labeling. Here we propose to discriminate between different protein conformations with multiple layers of graphene nanopore sensors by measuring the effect of protein-produced electrostatic potential (EP) on electric transport. Using conformations of the octapeptide Angiotensin II obtained through molecular dynamics simulations, we show that the EP critically depends on the geometries of constituent atoms and each conformation carries a unique EP signature. We then, using quantum transport simulations, reveal that these characteristic EP profiles cause distinctive modulation to electric charge densities of the graphene nanopores, leading to distinguishable changes in conductivity. Our results open the potential of label-free, single-molecule, and real-time detection of protein conformational changes.
Collapse
Affiliation(s)
- Wanzhi Qiu
- Centre for Neural Engineering, The University of Melbourne , 203 Bouverie Street, Carlton, Victoria 3053, Australia
| | | |
Collapse
|
8
|
Singh P, Sarkar SK, Bandyopadhyay P. Wang-Landau density of states based study of the folding-unfolding transition in the mini-protein Trp-cage (TC5b). J Chem Phys 2014; 141:015103. [DOI: 10.1063/1.4885726] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Priya Singh
- School of Physical Sciences, Jawaharlal Nehru University, New Delhi - 110 067, India
| | - Subir K. Sarkar
- School of Physical Sciences, Jawaharlal Nehru University, New Delhi - 110 067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi - 110 067, India
| |
Collapse
|
9
|
Wu X, Yang G, Zu Y, Zhou L. Molecular dynamics studies of β-hairpin folding with the presence of the sodium ion. Comput Biol Chem 2012; 38:1-9. [PMID: 22487489 DOI: 10.1016/j.compbiolchem.2012.02.003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2011] [Revised: 01/20/2012] [Accepted: 02/14/2012] [Indexed: 11/19/2022]
Abstract
Metal ions are ubiquitous in protein systems and play a significant role during their folding processes. Nineteen independent structures were determined for the Na(+)/β-hairpin interacting systems, and their folding pathways are different. (i) For Na(S47), the turn is rapidly shaped with the help of Na(+) and acts as the folding nucleus for the rest regions. Two intermediate states are observed and the resulted structure is the most folded. (ii) For Na(B41), Na(B52), Na(B54), Na(S55) and Na(B56), the inclusive Na(+) ions are anchored by β-strands. The local structures around the Na(+) ions and the turn regions fold simultaneously and serve as two independent folding nuclei. (iii) The other systems have no folding nuclei and correspond to low-folded structures. Long-range electrostatic interactions contribute a lot to the folding, especially from the four negatively charged residues (Glu42, Asp46, Asp47 and Glu56). The initial positions of the Na(+) ions are largely responsible for the different folding behaviors. The interactions with sidechain- rather than backbone-O atoms generally lead to more compact structures. Another factor affecting the folding is whether the O atoms are associated with native H-bonds, and those involved show decreased affinities to metal ions. The addition of water solvent does not induce obvious folding and conformational transitions to the Na(+)/β-hairpin interacting systems.
Collapse
Affiliation(s)
- Xiaomin Wu
- Key Laboratory of Forest Plant Ecology, Ministry of Education, Northeast Forestry University, Harbin 150040, PR China
| | | | | | | |
Collapse
|