1
|
Iyer SS, Srivastava A. Membrane lateral organization from potential energy disconnectivity graph. Biophys Chem 2024; 313:107284. [PMID: 39002248 DOI: 10.1016/j.bpc.2024.107284] [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: 03/29/2024] [Revised: 06/17/2024] [Accepted: 06/19/2024] [Indexed: 07/15/2024]
Abstract
Understanding the thermodynamic and kinetic properties of biomolecules requires elucidation of their complex energy landscape. A disconnectivity graph analysis of the energy landscape provides a framework for mapping the multi-dimensional landscape onto a two-dimensional representation while preserving the key features of the energy landscape. Several studies show that the structure or shape of the disconnectity graph is directly associated with the function of protein and nucleic acid molecules. In this review, we discuss how disconnectivity analysis of the potential energy surface can be extended to lipid molecules to glean important information about membrane organization. The shape of the disconnectivity graphs can be used to predict the lateral organization of multi-component lipid bilayer. We hope that this review encourages the use of disconnectivity graphs routinely by membrane biophysicists to predict the lateral organization of lipids.
Collapse
Affiliation(s)
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, C. V. Raman Road, Bangalore, Karnataka 560012, India.
| |
Collapse
|
2
|
Jiang H, Li H, Wong WH, Fan X. Revealing Free Energy Landscape From MD Data via Conditional Angle Partition Tree. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2023; 20:1384-1394. [PMID: 35503836 DOI: 10.1109/tcbb.2022.3172352] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Deciphering the free energy landscape of biomolecular structure space is crucial for understanding many complex molecular processes, such as protein-protein interaction, RNA folding, and protein folding. A major source of current dynamic structure data is Molecular Dynamics (MD) simulations. Several methods have been proposed to investigate the free energy landscape from MD data, but all of them rely on the assumption that kinetic similarity is associated with global geometric similarity, which may lead to unsatisfactory results. In this paper, we proposed a new method called Conditional Angle Partition Tree to reveal the hierarchical free energy landscape by correlating local geometric similarity with kinetic similarity. Its application on the benchmark alanine dipeptide MD data showed a much better performance than existing methods in exploring and understanding the free energy landscape. We also applied it to the MD data of Villin HP35. Our results are more reasonable on various aspects than those from other methods and very informative on the hierarchical structure of its energy landscape.
Collapse
|
3
|
Mirth J, Zhai Y, Bush J, Alvarado EG, Jordan H, Heim M, Krishnamoorthy B, Pflaum M, Clark A, Z Y, Adams H. Representations of energy landscapes by sublevelset persistent homology: An example with n-alkanes. J Chem Phys 2021; 154:114114. [PMID: 33752361 DOI: 10.1063/5.0036747] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Encoding the complex features of an energy landscape is a challenging task, and often, chemists pursue the most salient features (minima and barriers) along a highly reduced space, i.e., two- or three-dimensions. Even though disconnectivity graphs or merge trees summarize the connectivity of the local minima of an energy landscape via the lowest-barrier pathways, there is much information to be gained by also considering the topology of each connected component at different energy thresholds (or sublevelsets). We propose sublevelset persistent homology as an appropriate tool for this purpose. Our computations on the configuration phase space of n-alkanes from butane to octane allow us to conjecture, and then prove, a complete characterization of the sublevelset persistent homology of the alkane CmH2m+2 Potential Energy Landscapes (PELs), for all m, in all homological dimensions. We further compare both the analytical configurational PELs and sampled data from molecular dynamics simulation using the united and all-atom descriptions of the intramolecular interactions. In turn, this supports the application of distance metrics to quantify sampling fidelity and lays the foundation for future work regarding new metrics that quantify differences between the topological features of high-dimensional energy landscapes.
Collapse
Affiliation(s)
- Joshua Mirth
- Department of Mathematics, Colorado State University, Fort Collins, Colorado 80524, USA
| | - Yanqin Zhai
- Department of Nuclear, Plasma, and Radiological Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Johnathan Bush
- Department of Mathematics, Colorado State University, Fort Collins, Colorado 80524, USA
| | - Enrique G Alvarado
- Department of Mathematics and Statistics, Washington State University, Pullman, Washington 99164, USA
| | - Howie Jordan
- Department of Mathematics, University of Colorado, Boulder, Colorado 80309, USA
| | - Mark Heim
- Department of Mathematics, Colorado State University, Fort Collins, Colorado 80524, USA
| | - Bala Krishnamoorthy
- Department of Mathematics and Statistics, Washington State University, Vancouver, Washington 98686, USA
| | - Markus Pflaum
- Department of Mathematics, University of Colorado, Boulder, Colorado 80309, USA
| | - Aurora Clark
- Department of Chemistry, Washington State University, Pullman, Washington 99164, USA
| | - Y Z
- Department of Nuclear, Plasma, and Radiological Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Henry Adams
- Department of Mathematics, Colorado State University, Fort Collins, Colorado 80524, USA
| |
Collapse
|
4
|
Srivastava R. Application of Optimization Algorithms in Clusters. Front Chem 2021; 9:637286. [PMID: 33777900 PMCID: PMC7994592 DOI: 10.3389/fchem.2021.637286] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Accepted: 01/21/2021] [Indexed: 12/23/2022] Open
Abstract
The structural characterization of clusters or nanoparticles is essential to rationalize their size and composition-dependent properties. As experiments alone could not provide complete picture of cluster structures, so independent theoretical investigations are needed to find out a detail description of the geometric arrangement and corresponding properties of the clusters. The potential energy surfaces (PES) are explored to find several minima with an ultimate goal of locating the global minima (GM) for the clusters. Optimization algorithms, such as genetic algorithm (GA), basin hopping method and its variants, self-consistent basin-to-deformed-basin mapping, heuristic algorithm combined with the surface and interior operators (HA-SIO), fast annealing evolutionary algorithm (FAEA), random tunneling algorithm (RTA), and dynamic lattice searching (DLS) have been developed to solve the geometrical isomers in pure elemental clusters. Various model or empirical potentials (EPs) as Lennard-Jones (LJ), Born-Mayer, Gupta, Sutton-Chen, and Murrell-Mottram potentials are used to describe the bonding in different type of clusters. Due to existence of a large number of homotops in nanoalloys, genetic algorithm, basin-hopping algorithm, modified adaptive immune optimization algorithm (AIOA), evolutionary algorithm (EA), kick method and Knowledge Led Master Code (KLMC) are also used. In this review the optimization algorithms, computational techniques and accuracy of results obtained by using these mechanisms for different types of clusters will be discussed.
Collapse
|
5
|
Farris ACK, Seaton DT, Landau DP. Effects of lattice constraints in coarse-grained protein models. J Chem Phys 2021; 154:084903. [PMID: 33639740 DOI: 10.1063/5.0038184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We compare and contrast folding behavior in several coarse-grained protein models, both on- and off-lattice, in an attempt to uncover the effect of lattice constraints in these kinds of models. Using modern, extended ensemble Monte Carlo methods-Wang-Landau sampling, multicanonical sampling, replica-exchange Wang-Landau sampling, and replica-exchange multicanonical sampling, we investigate the thermodynamic and structural behavior of the protein Crambin within the context of the hydrophobic-polar, hydrophobic-"neutral"-polar (H0P), and semi-flexible H0P model frameworks. We uncover the folding process in all cases; all models undergo, at least, the two major structural transitions observed in nature-the coil-globule collapse and the folding transition. As the complexity of the model increases, these two major transitions begin to split into multi-step processes, wherein the lattice coarse-graining has a significant impact on the details of these processes. The results show that the level of structural coarse-graining is coupled to the level of interaction coarse-graining.
Collapse
Affiliation(s)
- Alfred C K Farris
- Department of Physics and Astronomy, Oxford College of Emory University, Oxford, Georgia 30054, USA
| | - Daniel T Seaton
- Open Learning, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - David P Landau
- Center for Simulational Physics, Department of Physics and Astronomy, The University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
6
|
Neelamraju S, Wales DJ, Gosavi S. Protein energy landscape exploration with structure-based models. Curr Opin Struct Biol 2020; 64:145-151. [DOI: 10.1016/j.sbi.2020.07.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Revised: 06/30/2020] [Accepted: 07/15/2020] [Indexed: 12/11/2022]
|
7
|
Abstract
A vast array of natural phenomena can be understood through the long-established schema of chemical bonding. Conventional chemical bonds arise through local gradients resulting from the rearrangement of electrons; however, it is possible that the hallmark features of chemical bonding could arise through local gradients resulting from nonelectronic forms of mediation. If other forms of mediation give rise to "bonds" that act like conventional ones, recognizing them as bonds could open new forms of supramolecular descriptions of phenomena at the nano- and microscales. Here, we show via a minimal model that crowded hard-particle systems governed solely by entropy exhibit the hallmark features of bonding despite the absence of chemical interactions. We quantitatively characterize these features and compare them to those exhibited by chemical bonds to argue for the existence of entropic bonds. As an example of the utility of the entropic bond classification, we demonstrate the nearly equivalent tradeoff between chemical bonds and entropic bonds in the colloidal crystallization of hard hexagonal nanoplates.
Collapse
|
8
|
Neelamraju S, Wales DJ, Gosavi S. Go-Kit: A Tool To Enable Energy Landscape Exploration of Proteins. J Chem Inf Model 2019; 59:1703-1708. [PMID: 30977648 DOI: 10.1021/acs.jcim.9b00007] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Coarse-grained Go̅-like models, based on the principle of minimal frustration, provide valuable insight into fundamental questions in the field of protein folding and dynamics. In conjunction with commonly used molecular dynamics (MD) simulations, energy landscape exploration methods like discrete path sampling (DPS) with Go̅-like models can provide quantitative details of the thermodynamics and kinetics of proteins. Here we present Go-kit, a software that facilitates the setup of MD and DPS simulations of several flavors of Go̅-like models. Go-kit is designed for use with MD (GROMACS) and DPS (PATHSAMPLE) simulation engines that are open source. The Go-kit code is written in python2.7 and is also open source. A case study for the ribosomal protein S6 is discussed to illustrate the utility of the software, which is available at https://github.com/gokit1/gokit .
Collapse
Affiliation(s)
- Sridhar Neelamraju
- Simons Centre for the Study of Living Machines , National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bellary Road , Bangalore 560065 , India.,University Chemical Laboratories , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , U.K
| | - David J Wales
- University Chemical Laboratories , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , U.K
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines , National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bellary Road , Bangalore 560065 , India
| |
Collapse
|
9
|
Hey JC, Doyle EJ, Chen Y, Johnston RL. Isomers and energy landscapes of micro-hydrated sulfite and chlorate clusters. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2018; 376:rsta.2017.0154. [PMID: 29431682 PMCID: PMC5805918 DOI: 10.1098/rsta.2017.0154] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Accepted: 11/22/2017] [Indexed: 06/08/2023]
Abstract
We present putative global minima for the micro-hydrated sulfite SO32-(H2O) N and chlorate ClO3-(H2O) N systems in the range 3≤N≤15 found using basin-hopping global structure optimization with an empirical potential. We present a structural analysis of the hydration of a large number of minimized structures for hydrated sulfite and chlorate clusters in the range 3≤N≤50. We show that sulfite is a significantly stronger net acceptor of hydrogen bonding within water clusters than chlorate, completely suppressing the appearance of hydroxyl groups pointing out from the cluster surface (dangling OH bonds), in low-energy clusters. We also present a qualitative analysis of a highly explored energy landscape in the region of the global minimum of the eight water hydrated sulfite and chlorate systems.This article is part of the theme issue 'Modern theoretical chemistry'.
Collapse
Affiliation(s)
- John C Hey
- School of Chemistry, University of Birmingham, Birmingham B15 2TT, UK
| | - Emily J Doyle
- School of Chemistry, University of Birmingham, Birmingham B15 2TT, UK
| | - Yuting Chen
- School of Chemistry, University of Birmingham, Birmingham B15 2TT, UK
| | - Roy L Johnston
- School of Chemistry, University of Birmingham, Birmingham B15 2TT, UK
| |
Collapse
|
10
|
|
11
|
Koorehdavoudi H, Bogdan P. A Statistical Physics Characterization of the Complex Systems Dynamics: Quantifying Complexity from Spatio-Temporal Interactions. Sci Rep 2016; 6:27602. [PMID: 27297496 PMCID: PMC4906350 DOI: 10.1038/srep27602] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2016] [Accepted: 05/20/2016] [Indexed: 12/28/2022] Open
Abstract
Biological systems are frequently categorized as complex systems due to their capabilities of generating spatio-temporal structures from apparent random decisions. In spite of research on analyzing biological systems, we lack a quantifiable framework for measuring their complexity. To fill this gap, in this paper, we develop a new paradigm to study a collective group of N agents moving and interacting in a three-dimensional space. Our paradigm helps to identify the spatio-temporal states of the motion of the group and their associated transition probabilities. This framework enables the estimation of the free energy landscape corresponding to the identified states. Based on the energy landscape, we quantify missing information, emergence, self-organization and complexity for a collective motion. We show that the collective motion of the group of agents evolves to reach the most probable state with relatively lowest energy level and lowest missing information compared to other possible states. Our analysis demonstrates that the natural group of animals exhibit a higher degree of emergence, self-organization and complexity over time. Consequently, this algorithm can be integrated into new frameworks to engineer collective motions to achieve certain degrees of emergence, self-organization and complexity.
Collapse
Affiliation(s)
- Hana Koorehdavoudi
- Department of Aerospace and Mechanical Engineering, University of Southern California, Los Angeles, CA 90089-1453, USA
| | - Paul Bogdan
- Department of Electrical Engineering, University of Southern California, Los Angeles, CA 90089-2560, USA
| |
Collapse
|
12
|
Hey JC, Smeeton LC, Oakley MT, Johnston RL. Isomers and Energy Landscapes of Perchlorate-Water Clusters and a Comparison to Pure Water and Sulfate-Water Clusters. J Phys Chem A 2016; 120:4008-15. [PMID: 27223243 DOI: 10.1021/acs.jpca.6b01495] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Hydrated ions are crucially important in a wide array of environments, from biology to the atmosphere, and the presence and concentration of ions in a system can drastically alter its behavior. One way in which ions can affect systems is in their interactions with proteins. The Hofmeister series ranks ions by their ability to salt-out proteins, with kosmotropes, such as sulfate, increasing their stability and chaotropes, such as perchlorate, decreasing their stability. We study hydrated perchlorate clusters as they are strongly chaotropic and thus exhibit different properties than sulfate. In this study we simulate small hydrated perchlorate clusters using a basin-hopping geometry optimization search with empirical potentials. We compare topological features of these clusters to data from both computational and experimental studies of hydrated sulfate ions and draw some conclusions about ion effects in the Hofmeister series. We observe a patterning conferred to the water molecules within the cluster by the presence of the perchlorate ion and compare the magnitude of this effect to that observed in previous studies involving sulfate. We also investigate the influence of the overall ionic charge on the low-energy structures adopted by these clusters.
Collapse
Affiliation(s)
- John C Hey
- School of Chemistry, University of Birmingham , Edgbaston, Birmingham B15 2TT, U.K
| | - Lewis C Smeeton
- School of Chemistry, University of Birmingham , Edgbaston, Birmingham B15 2TT, U.K
| | - Mark T Oakley
- School of Chemistry, University of Birmingham , Edgbaston, Birmingham B15 2TT, U.K
| | - Roy L Johnston
- School of Chemistry, University of Birmingham , Edgbaston, Birmingham B15 2TT, U.K
| |
Collapse
|
13
|
Neelamraju S, Oakley MT, Johnston RL. Chiral effects on helicity studied via the energy landscape of short (d, l)-alanine peptides. J Chem Phys 2015; 143:165103. [DOI: 10.1063/1.4933428] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sridhar Neelamraju
- School of Chemistry, University of Birmingham, Edgbaston B15 2TT, United Kingdom
| | - Mark T. Oakley
- School of Chemistry, University of Birmingham, Edgbaston B15 2TT, United Kingdom
| | - Roy L. Johnston
- School of Chemistry, University of Birmingham, Edgbaston B15 2TT, United Kingdom
| |
Collapse
|
14
|
Smeeton LC, Farrell JD, Oakley MT, Wales DJ, Johnston RL. Structures and Energy Landscapes of Hydrated Sulfate Clusters. J Chem Theory Comput 2015; 11:2377-84. [DOI: 10.1021/acs.jctc.5b00151] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Lewis C. Smeeton
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| | - James D. Farrell
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Mark T. Oakley
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| | - David J. Wales
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Roy L. Johnston
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| |
Collapse
|
15
|
Heard CJ, Johnston RL, Schön JC. Energy Landscape Exploration of Sub-Nanometre Copper-Silver Clusters. Chemphyschem 2015; 16:1461-9. [DOI: 10.1002/cphc.201402887] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2014] [Indexed: 11/12/2022]
|
16
|
Smeeton LC, Oakley MT, Johnston RL. Visualizing energy landscapes with metric disconnectivity graphs. J Comput Chem 2014; 35:1481-90. [PMID: 24866379 PMCID: PMC4285870 DOI: 10.1002/jcc.23643] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Revised: 03/12/2014] [Accepted: 04/14/2014] [Indexed: 11/24/2022]
Abstract
The visualization of multidimensional energy landscapes is important, providing insight into the kinetics and thermodynamics of a system, as well the range of structures a system can adopt. It is, however, highly nontrivial, with the number of dimensions required for a faithful reproduction of the landscape far higher than can be represented in two or three dimensions. Metric disconnectivity graphs provide a possible solution, incorporating the landscape connectivity information present in disconnectivity graphs with structural information in the form of a metric. In this study, we present a new software package, PyConnect, which is capable of producing both disconnectivity graphs and metric disconnectivity graphs in two or three dimensions. We present as a test case the analysis of the 69-bead BLN coarse-grained model protein and show that, by choosing appropriate order parameters, metric disconnectivity graphs can resolve correlations between structural features on the energy landscape with the landscapes energetic and kinetic properties.
Collapse
Affiliation(s)
- Lewis C Smeeton
- School of Chemistry, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom
| | | | | |
Collapse
|