1
|
Jin J, Reichman DR. Hierarchical Framework for Predicting Entropies in Bottom-Up Coarse-Grained Models. J Phys Chem B 2024; 128:3182-3199. [PMID: 38507575 DOI: 10.1021/acs.jpcb.3c07624] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
The thermodynamic entropy of coarse-grained (CG) models stands as one of the most important properties for quantifying the missing information during the CG process and for establishing transferable (or extendible) CG interactions. However, performing additional CG simulations on top of model construction often leads to significant additional computational overhead. In this work, we propose a simple hierarchical framework for predicting the thermodynamic entropies of various molecular CG systems. Our approach employs a decomposition of the CG interactions, enabling the estimation of the CG partition function and thermodynamic properties a priori. Starting from the ideal gas description, we leverage classical perturbation theory to systematically incorporate simple yet essential interactions, ranging from the hard sphere model to the generalized van der Waals model. Additionally, we propose an alternative approach based on multiparticle correlation functions, allowing for systematic improvements through higher-order correlations. Numerical applications to molecular liquids validate the high fidelity of our approach, and our computational protocols demonstrate that a reduced model with simple energetics can reasonably estimate the thermodynamic entropy of CG models without performing any CG simulations. Overall, our findings present a systematic framework for estimating not only the entropy but also other thermodynamic properties of CG models, relying solely on information from the reference system.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
2
|
Nam K, Shao Y, Major DT, Wolf-Watz M. Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development. ACS OMEGA 2024; 9:7393-7412. [PMID: 38405524 PMCID: PMC10883025 DOI: 10.1021/acsomega.3c09084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 01/15/2024] [Accepted: 01/19/2024] [Indexed: 02/27/2024]
Abstract
Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.
Collapse
Affiliation(s)
- Kwangho Nam
- Department
of Chemistry and Biochemistry, University
of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yihan Shao
- Department
of Chemistry and Biochemistry, University
of Oklahoma, Norman, Oklahoma 73019-5251, United States
| | - Dan T. Major
- Department
of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | | |
Collapse
|
3
|
Shakibaie MR, Modaresi F, Azizi O, Tadjrobehkar O, Ghaemi MM. Amphiphilic peptide Mastoparan-B induces conformational changes within the AdeB efflux pump, down-regulates adeB gene expression, and restores antibiotic susceptibility in an MDR strain of Acinetobacter baumannii. Proteins 2023; 91:1205-1221. [PMID: 37455426 DOI: 10.1002/prot.26539] [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: 12/16/2022] [Revised: 05/11/2023] [Accepted: 05/30/2023] [Indexed: 07/18/2023]
Abstract
Mastoparan B (MP-B) is an amphiphilic peptide with a potent antimicrobial activity against most Gram-negative bacteria. However, there is little information available on the inhibition of the Acinetobacter baumannii resistance-nodulation-cell-division (RND) efflux pump using this antimicrobial peptide. Here, we carried out a series of in-silico experiments to find the mechanisms underlying the anti-efflux activity of MP-B using a multi-drug resistant (MDR) strain of A. baumannii (AB). According to our findings, MP-B demonstrated a potent antibacterial activity against an MDR-AB (minimum inhibitory concentration [MIC] = 1 μg/mL) followed by a 20-fold reduction in the adeB gene expression in the presence of sub-MIC of this peptide. Using Groningen Machine for Chemicals Simulation (GROMACS) via PyMOL Graphical User Interface (GUI), (we observed that, the AdeB transporter had conserved helix-turn-helix regions and a tight pore rich in Phe and Ala residues. To understand how inhibition of the AdeB is achieved, we generated 20 apo-MP-B poses using the InterPep and SiteMap tools. The high-quality model was created by homology modeling and used for docking via AutoDock/Vina to identify the MP-B binding sites. We established that the most apo-MP-B formed H-bonds to the backbone of five amino acids in the Helix-5. As a result, the dihedral angles of the involved amino acids shift by 9.0-9.6 Ǻ, causing a change in the conformation of the AdeB protein. This led to helix conformation stereoisomerization and block the AdeB activity. MP-B presumably has dual mechanisms. (1) It blocks the AdeB transporter by changing its conformation. (2) MP-B influences the adeB gene expression by binding to G-protein which laterally controls efflux regulators like MarA, RamA, SoxS, and Rob proteins.
Collapse
Affiliation(s)
- Mohammad Reza Shakibaie
- Department of Microbiology and Virology, Kerman University of Medical Sciences, Kerman, Iran
- Gastroenterology Hepatology Research Center, Institute of Basic and Clinical Physiology, Kerman University of Medical Sciences, Kerman, Iran
| | - Farzan Modaresi
- Department of Microbiology, School of Medicine, Jahrom University of Medical Sciences, Jahrom, Iran
| | - Omid Azizi
- Department of Laboratory Sciences, and Health Sciences Research center, Torbat Heydariyeh University of Medical Sciences, Torbate Heydarieh, Iran
| | - Omid Tadjrobehkar
- Department of Microbiology and Virology, Kerman University of Medical Sciences, Kerman, Iran
| | - Mohammad Mehdi Ghaemi
- Medical Informatics Research Center, Institute for Futures Studies in Health, Kerman University of Medical Sciences, Kerman, Iran
| |
Collapse
|
4
|
Zhang Y, He L, Li S. Temperature dependence of DNA elasticity: An all-atom molecular dynamics simulation study. J Chem Phys 2023; 158:094902. [PMID: 36889965 DOI: 10.1063/5.0138940] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023] Open
Abstract
We used all-atom molecular dynamics simulation to investigate the elastic properties of double-stranded DNA (dsDNA). We focused on the influences of temperature on the stretch, bend, and twist elasticities, as well as the twist-stretch coupling, of the dsDNA over a wide range of temperature. The results showed that the bending and twist persistence lengths, together with the stretch and twist moduli, decrease linearly with temperature. However, the twist-stretch coupling behaves in a positive correction and enhances as the temperature increases. The potential mechanisms of how temperature affects dsDNA elasticity and coupling were investigated by using the trajectories from atomistic simulation, in which thermal fluctuations in structural parameters were analyzed in detail. We analyzed the simulation results by comparing them with previous simulation and experimental data, which are in good agreement. The prediction about the temperature dependence of dsDNA elastic properties provides a deeper understanding of DNA elasticities in biological environments and potentially helps in the further development of DNA nanotechnology.
Collapse
Affiliation(s)
- Yahong Zhang
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| | - Linli He
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| | - Shiben Li
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| |
Collapse
|
5
|
Jin J, Voth GA. Statistical Mechanical Design Principles for Coarse-Grained Interactions across Different Conformational Free Energy Surfaces. J Phys Chem Lett 2023; 14:1354-1362. [PMID: 36728761 PMCID: PMC9940719 DOI: 10.1021/acs.jpclett.2c03844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Accepted: 01/27/2023] [Indexed: 06/18/2023]
Abstract
Systematic bottom-up coarse-graining (CG) of molecular systems provides a means to explore different coupled length and time scales while treating the molecular-scale physics at a reduced level. However, the configuration dependence of CG interactions often results in CG models with limited applicability for exploring the parametrized configurations. We propose a statistical mechanical theory to design CG interactions across different configurations and conditions. In order to span wide ranges of conformational space, distinct classical CG free energy surfaces for characteristic configurations are identified using molecular collective variables. The coupling interaction between different CG free energy surfaces can then be systematically determined by analogy to quantum mechanical approaches describing coupled states. The present theory can accurately capture the underlying many-body potentials of mean force in the CG variables for various order parameters applied to liquids, interfaces, and in principle proteins, uncovering the complex nature underlying the coupling interaction and imparting a new protocol for the design of predictive multiscale models.
Collapse
Affiliation(s)
| | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
6
|
Gorai B, Vashisth H. Progress in Simulation Studies of Insulin Structure and Function. Front Endocrinol (Lausanne) 2022; 13:908724. [PMID: 35795141 PMCID: PMC9252437 DOI: 10.3389/fendo.2022.908724] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 04/28/2022] [Indexed: 01/02/2023] Open
Abstract
Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.
Collapse
Affiliation(s)
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, United States
| |
Collapse
|
7
|
Rizzuti B. Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2022; 1870:140757. [PMID: 35051666 DOI: 10.1016/j.bbapap.2022.140757] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Revised: 01/08/2022] [Accepted: 01/10/2022] [Indexed: 12/22/2022]
Abstract
Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.
Collapse
Affiliation(s)
- Bruno Rizzuti
- CNR-NANOTEC, SS Rende (CS), Department of Physics, University of Calabria, 87036 Rende, Italy; Institute for Biocomputation and Physics of Complex Systems (BIFI), Joint Unit GBsC-CSIC-BIFI, University of Zaragoza, 50018 Zaragoza, Spain.
| |
Collapse
|
8
|
Najafi H, Jafari M, Farahavar G, Abolmaali SS, Azarpira N, Borandeh S, Ravanfar R. Recent advances in design and applications of biomimetic self-assembled peptide hydrogels for hard tissue regeneration. Biodes Manuf 2021; 4:735-756. [PMID: 34306798 PMCID: PMC8294290 DOI: 10.1007/s42242-021-00149-0] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Accepted: 06/12/2021] [Indexed: 12/22/2022]
Abstract
Abstract The development of natural biomaterials applied for hard tissue repair and regeneration is of great importance, especially in societies with a large elderly population. Self-assembled peptide hydrogels are a new generation of biomaterials that provide excellent biocompatibility, tunable mechanical stability, injectability, trigger capability, lack of immunogenic reactions, and the ability to load cells and active pharmaceutical agents for tissue regeneration. Peptide-based hydrogels are ideal templates for the deposition of hydroxyapatite crystals, which can mimic the extracellular matrix. Thus, peptide-based hydrogels enhance hard tissue repair and regeneration compared to conventional methods. This review presents three major self-assembled peptide hydrogels with potential application for bone and dental tissue regeneration, including ionic self-complementary peptides, amphiphilic (surfactant-like) peptides, and triple-helix (collagen-like) peptides. Special attention is given to the main bioactive peptides, the role and importance of self-assembled peptide hydrogels, and a brief overview on molecular simulation of self-assembled peptide hydrogels applied for bone and dental tissue engineering and regeneration. Graphic abstract
Collapse
Affiliation(s)
- Haniyeh Najafi
- Pharmaceutical Nanotechnology Department, Shiraz University of Medical Sciences, 71345-1583 Shiraz, Iran
| | - Mahboobeh Jafari
- Pharmaceutical Nanotechnology Department, Shiraz University of Medical Sciences, 71345-1583 Shiraz, Iran
| | - Ghazal Farahavar
- Pharmaceutical Nanotechnology Department, Shiraz University of Medical Sciences, 71345-1583 Shiraz, Iran
| | - Samira Sadat Abolmaali
- Pharmaceutical Nanotechnology Department, Shiraz University of Medical Sciences, 71345-1583 Shiraz, Iran
- Center for Nanotechnology in Drug Delivery, Shiraz University of Medical Sciences, 71345-1583 Shiraz, Iran
| | - Negar Azarpira
- Transplant Research Center, Shiraz University of Medical Sciences, Mohammad Rasoul-Allah Research Tower, 7193711351 Shiraz, Iran
| | - Sedigheh Borandeh
- Center for Nanotechnology in Drug Delivery, Shiraz University of Medical Sciences, 71345-1583 Shiraz, Iran
- Polymer Technology Research Group, Department of Chemical and Metallurgical Engineering, Aalto University, 02152 Espoo, Finland
| | - Raheleh Ravanfar
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125 USA
| |
Collapse
|
9
|
Schlick T, Portillo-Ledesma S, Myers CG, Beljak L, Chen J, Dakhel S, Darling D, Ghosh S, Hall J, Jan M, Liang E, Saju S, Vohr M, Wu C, Xu Y, Xue E. Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field. Annu Rev Biophys 2021; 50:267-301. [PMID: 33606945 PMCID: PMC8105287 DOI: 10.1146/annurev-biophys-091720-102019] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, New York 10003, USA;
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200122, China
| | | | - Christopher G Myers
- Department of Chemistry, New York University, New York, New York 10003, USA;
| | - Lauren Beljak
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Justin Chen
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sami Dakhel
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Daniel Darling
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sayak Ghosh
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Joseph Hall
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mikaeel Jan
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Emily Liang
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sera Saju
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mackenzie Vohr
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Chris Wu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Yifan Xu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Eva Xue
- College of Arts and Science, New York University, New York, New York 10003, USA
| |
Collapse
|
10
|
Schlick T, Portillo-Ledesma S. Biomolecular modeling thrives in the age of technology. NATURE COMPUTATIONAL SCIENCE 2021; 1:321-331. [PMID: 34423314 PMCID: PMC8378674 DOI: 10.1038/s43588-021-00060-9] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 03/22/2021] [Indexed: 12/12/2022]
Abstract
The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, NY, USA
- Courant Institute of Mathematical Sciences, New York University, New York, NY, USA
- New York University–East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai, China
| | | |
Collapse
|
11
|
Jin J, Pak AJ, Han Y, Voth GA. A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). II. Temperature transferability and structural properties at low temperature. J Chem Phys 2021; 154:044105. [PMID: 33514078 PMCID: PMC7826166 DOI: 10.1063/5.0026652] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Accepted: 12/14/2020] [Indexed: 11/14/2022] Open
Abstract
A number of studies have constructed coarse-grained (CG) models of water to understand its anomalous properties. Most of these properties emerge at low temperatures, and an accurate CG model needs to be applicable to these low-temperature ranges. However, direct use of CG models parameterized from other temperatures, e.g., room temperature, encounters a problem known as transferability, as the CG potential essentially follows the form of the many-body CG free energy function. Therefore, temperature-dependent changes to CG interactions must be accounted for. The collective behavior of water at low temperature is generally a many-body process, which often motivates the use of expensive many-body terms in the CG interactions. To surmount the aforementioned problems, we apply the Bottom-Up Many-Body Projected Water (BUMPer) CG model constructed from Paper I to study the low-temperature behavior of water. We report for the first time that the embedded three-body interaction enables BUMPer, despite its pairwise form, to capture the growth of ice at the ice/water interface with corroborating many-body correlations during the crystal growth. Furthermore, we propose temperature transferable BUMPer models that are indirectly constructed from the free energy decomposition scheme. Changes in CG interactions and corresponding structures are faithfully recapitulated by this framework. We further extend BUMPer to examine its ability to predict the structure, density, and diffusion anomalies by employing an alternative analysis based on structural correlations and pairwise potential forms to predict such anomalies. The presented analysis highlights the existence of these anomalies in the low-temperature regime and overcomes potential transferability problems.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Alexander J. Pak
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Yining Han
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
12
|
Badu S, Melnik R, Singh S. Mathematical and computational models of RNA nanoclusters and their applications in data-driven environments. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1804564] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Shyam Badu
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, Ontario, Canada
| | - Roderick Melnik
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, Ontario, Canada
- BCAM-Basque Center for Applied Mathematics, Bilbao, Spain
| | - Sundeep Singh
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, Ontario, Canada
| |
Collapse
|
13
|
Seshan G, Kanagasabai S, Ananthasri S, Kannappan B, Suvitha A, Jaimohan SM, Kanagaraj S, Kothandan G. Insights of structure-based pharmacophore studies and inhibitor design against Gal3 receptor through molecular dynamics simulations. J Biomol Struct Dyn 2020; 39:6987-6999. [PMID: 32772816 DOI: 10.1080/07391102.2020.1804452] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Our present work studies the structure-based pharmacophore modeling and designing inhibitor against Gal3 receptor through molecular dynamics (MD) simulations extensively. Pharmacophore models play a key role in computer-aided drug discovery like in the case of virtual screening of chemical databases, de novo drug design and lead optimization. Structure-based methods for developing pharmacophore models are important, and there have been a number of studies combining such methods with the use of MD simulations to model protein's flexibility. The two potential antagonists SNAP 37889 and SNAP 398299 were docked and simulated for 250 ns and the results are analyzed and carried for the structure-based pharmacophore studies. This helped in identification of the subtype selectivity of the binding sites of the Gal3 receptor. Our work mainly focuses on identifying these binding site residues and to design more potent inhibitors compared to the previously available inhibitors through pharmacophore models. The study provides crucial insight into the binding site residues Ala2, Asp3, Ala4, Gln5, Phe24, Gln79, Ala80, Ile82, Tyr83, Trp88, His99, Ile102, Tyr103, Met106, Tyr157, Tyr161, Pro174, Trp176, Arg181, Ala183, Leu184, Asp185, Thr188, Trp248, His251, His252, Ile255, Leu256, Phe258, Trp259, Tyr270, Arg273, Leu274 and His277, which plays a significant role in the conformational changes of the receptor and helps to understand the inhibition mechanism. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Gunalan Seshan
- Biopolymer Modelling Laboratory, Centre of Advanced Study in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai, Tamil Nadu, India
| | - Somarathinam Kanagasabai
- Biopolymer Modelling Laboratory, Centre of Advanced Study in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai, Tamil Nadu, India
| | - Sailapathi Ananthasri
- Biopolymer Modelling Laboratory, Centre of Advanced Study in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai, Tamil Nadu, India
| | - Balaji Kannappan
- National Research Center for Dementia, Department of Life Science, Chosun University, Gwangju, South Korea
| | - A Suvitha
- Department of Physics, CMR Institute of Technology, Bangalore, Karnataka, India
| | - S M Jaimohan
- Advanced Materials Laboratory, CSIR-Central Leather Research Institute, Chennai, Tamil Nadu, India
| | - Sekar Kanagaraj
- Laboratory for Structural Biology and Biocomputing, Department of Computational and Data Sciences, Indian Institute of Science, Bangalore, Karnataka, India
| | - Gugan Kothandan
- Biopolymer Modelling Laboratory, Centre of Advanced Study in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai, Tamil Nadu, India
| |
Collapse
|
14
|
Hinkle KR, Phelan FR. Solvation Free Energy of Self-Assembled Complexes: Using Molecular Dynamics to Understand the Separation of ssDNA-Wrapped Single-Walled Carbon Nanotubes. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2020; 124:10.1021/acs.jpcc.0c00983. [PMID: 34136061 PMCID: PMC8204635 DOI: 10.1021/acs.jpcc.0c00983] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Molecular dynamics simulations were used to characterize the self-assembly of single-stranded DNA (ssDNA) on a (6,5) single-walled carbon nanotube (SWCNT) in aqueous solution for the purpose of gaining an improved theoretical understanding of separation strategies for SWCNTs using ssDNA as a dispersant. Four separate ssDNA sequences, ((TAT)4, TTA(TAT)2ATT, C12, (GTC)2GT), at various levels of loading, were chosen for study based on published experimental work showing selective extraction of particular SWCNT species based on the ssDNA dispersant sequence. We develop a unique workflow based on free energy perturbation (FEP) and use this to determine the relative solubility of these complexes due to the adsorption of the ssDNA on the SWCNT surface, and hence, rank the favorability of separations observed during experiments. Results qualitatively agree with experiments and indicate that the nucleobase sequence of the adsorbed ssDNA greatly affects the free energy of complex solvation which ultimately drives SWCNT separation. Further, to elucidate the underlying physics governing the ssDNA-SWCNT solubility rankings, we also present calculations for four structural characteristics of ssDNA adsorption. We demonstrate that a unique type of intra-strand hydrogen bonding is the most important factor contributing to the stability of the ssDNA-SWCNT complexes and show how these adsorption characteristics are coupled with the FEP results.
Collapse
Affiliation(s)
- Kevin R. Hinkle
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland
- Department of Chemical and Materials Engineering, University of Dayton, Dayton, OH
| | - Frederick R. Phelan
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland
| |
Collapse
|
15
|
Abstract
Low resolution coarse-grained (CG) models are widely adopted for investigating phenomena that cannot be effectively simulated with all-atom (AA) models. Since the development of the many-body dissipative particle dynamics method, CG models have increasingly supplemented conventional pair potentials with one-body potentials of the local density (LD) around each site. These LD potentials appear to significantly extend the transferability of CG models, while also enabling more accurate descriptions of thermodynamic properties, interfacial phenomena, and many-body correlations. In this work, we systematically examine the properties of LD potentials. We first derive and numerically demonstrate a nontrivial transformation of pair and LD potentials that leaves the total forces and equilibrium distribution invariant. Consequently, the pair and LD potentials determined via bottom-up methods are not unique. We then investigate the sensitivity of CG models for glycerol to the weighting function employed for defining the local density. We employ the multiscale coarse-graining (MS-CG) method to simultaneously parameterize both pair and LD potentials. When employing a short-ranged Lucy function that defines the local density from the first solvation shell, the MS-CG model accurately reproduces the pair structure, pressure-density equation of state, and liquid-vapor interfacial profile of the AA model. The accuracy of the model generally decreases as the range of the Lucy function increases further. The MS-CG model provides similar accuracy when a smoothed Heaviside function is employed to define the local density from the first solvation shell. However, the model performs less well when this function acts on either longer or shorter length scales.
Collapse
Affiliation(s)
- Michael R DeLyser
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W G Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
16
|
Lebold KM, Noid WG. Dual-potential approach for coarse-grained implicit solvent models with accurate, internally consistent energetics and predictive transferability. J Chem Phys 2019; 151:164113. [PMID: 31675902 DOI: 10.1063/1.5125246] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The dual-potential approach promises coarse-grained (CG) models that accurately reproduce both structural and energetic properties, while simultaneously providing predictive estimates for the temperature-dependence of the effective CG potentials. In this work, we examine the dual-potential approach for implicit solvent CG models that reflect large entropic effects from the eliminated solvent. Specifically, we construct implicit solvent models at various resolutions, R, by retaining a fraction 0.10 ≤ R ≤ 0.95 of the molecules from a simple fluid of Lennard-Jones spheres. We consider the dual-potential approach in both the constant volume and constant pressure ensembles across a relatively wide range of temperatures. We approximate the many-body potential of mean force for the remaining solutes with pair and volume potentials, which we determine via multiscale coarse-graining and self-consistent pressure-matching, respectively. Interestingly, with increasing temperature, the pair potentials appear increasingly attractive, while the volume potentials become increasingly repulsive. The dual-potential approach not only reproduces the atomic energetics but also quite accurately predicts this temperature-dependence. We also derive an exact relationship between the thermodynamic specific heat of an atomic model and the energetic fluctuations that are observable at the CG resolution. With this generalized fluctuation relationship, the approximate CG models quite accurately reproduce the thermodynamic specific heat of the underlying atomic model.
Collapse
Affiliation(s)
- Kathryn M Lebold
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W G Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
17
|
Jin J, Pak AJ, Voth GA. Understanding Missing Entropy in Coarse-Grained Systems: Addressing Issues of Representability and Transferability. J Phys Chem Lett 2019; 10:4549-4557. [PMID: 31319036 PMCID: PMC6782054 DOI: 10.1021/acs.jpclett.9b01228] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Coarse-grained (CG) models facilitate efficient simulation of complex systems by integrating out the atomic, or fine-grained (FG), degrees of freedom. Systematically derived CG models from FG simulations often attempt to approximate the CG potential of mean force (PMF), an inherently multidimensional and many-body quantity, using additive pairwise contributions. However, they currently lack fundamental principles that enable their extensible use across different thermodynamic state points, i.e., transferability. In this work, we investigate the explicit energy-entropy decomposition of the CG PMF as a means to construct transferable CG models. In particular, despite its high-dimensional nature, we find for liquid systems that the entropic component to the CG PMF can similarly be represented using additive pairwise contributions, which we show is highly coupled to the CG configurational entropy. This approach formally connects the missing entropy that is lost due to the CG representation, i.e., translational, rotational, and vibrational modes associated with the missing degrees of freedom, to the CG entropy. By design, the present framework imparts transferable CG interactions across different temperatures due to the explicit definition of an additive entropic contribution. Furthermore, we demonstrate that transferability across composition state points, such as between bulk liquids and their mixtures, is also achieved by designing combining rules to approximate cross-interactions from bulk CG PMFs. Using the predicted CG model for liquid mixtures, structural correlations of the fitted CG model were found to corroborate a high-fidelity combining rule. Our findings elucidate the physical nature and compact representation of CG entropy and suggest a new approach for overcoming the transferability problem. We expect that this approach will further extend the current view of CG modeling into predictive multiscale modeling.
Collapse
|
18
|
Marrink SJ, Corradi V, Souza PC, Ingólfsson HI, Tieleman DP, Sansom MS. Computational Modeling of Realistic Cell Membranes. Chem Rev 2019; 119:6184-6226. [PMID: 30623647 PMCID: PMC6509646 DOI: 10.1021/acs.chemrev.8b00460] [Citation(s) in RCA: 435] [Impact Index Per Article: 87.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Indexed: 12/15/2022]
Abstract
Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.
Collapse
Affiliation(s)
- Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Paulo C.T. Souza
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Mark S.P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
19
|
Atomistic molecular dynamics simulations of the LCST conformational transition in poly(N-vinylcaprolactam) in water. J Mol Graph Model 2019; 90:51-58. [PMID: 31009934 DOI: 10.1016/j.jmgm.2019.04.004] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2019] [Revised: 04/06/2019] [Accepted: 04/09/2019] [Indexed: 11/22/2022]
Abstract
Thermoresponsive poly(N-vinylcaprolactam) (PVCL) has received growing interest due to a temperature-induced phase transition, which switches its solubility in aqueous solutions. However, the lower critical solution temperature (LCST) of PVCL is greatly influenced by the molecular weight, morphology and the environment. Therefore, despite of numerous experimental studies of the thermal response of PVCL, a driving force and a molecular origin of conformation transitions in solution remain far less studied. To get a better understanding of the coil-to-globule conformation transition of PVCL in aqueous solution, we examined the structure and conformation dynamics of a single-chain PVCL30 in a temperature range of 280-360 K by using atomistic molecular dynamics (MD) simulations. The united-atom GROMOS G53a6 force field was re-parameterized and fine-tuned by DFT calculations to reproduce the experimental LCST transition of PVCL. Our MD model reproduces the LCST transition of PVCL30 to occur within a temperature range of 34.6-38.5°. MD simulation results suggest a significant difference between the hydration state of the carbonyl group of PVCL below and above the LCST threshold. The analysis of the number of hydrogen bonds of PVCL with water molecules demonstrates that dehydration of the polymer plays an important role and drives the temperature-induced polymer collapse. Finally, the developed MD model and FF parameters were successfully tested for large-scale systems, such as mixture PVCL30 oligomer and single-chain PVCL816 polymer, respectively.
Collapse
|
20
|
Bansal R, Care A, Lord MS, Walsh TR, Sunna A. Experimental and theoretical tools to elucidate the binding mechanisms of solid-binding peptides. N Biotechnol 2019; 52:9-18. [PMID: 30954671 DOI: 10.1016/j.nbt.2019.04.001] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Revised: 03/31/2019] [Accepted: 04/01/2019] [Indexed: 12/19/2022]
Abstract
The interactions between biomolecules and solid surfaces play an important role in designing new materials and applications which mimic nature. Recently, solid-binding peptides (SBPs) have emerged as potential molecular building blocks in nanobiotechnology. SBPs exhibit high selectivity and binding affinity towards a wide range of inorganic and organic materials. Although these peptides have been widely used in various applications, there is a need to understand the interaction mechanism between the peptide and its material substrate, which is challenging both experimentally and theoretically. This review describes the main characterisation techniques currently available to study SBP-surface interactions and their contribution to gain a better insight for designing new peptides for tailored binding.
Collapse
Affiliation(s)
- Rachit Bansal
- Department of Molecular Sciences, Macquarie University, Sydney, NSW 2109, Australia; ARC Centre of Excellence for Nanoscale Biophotonics, Macquarie University, Sydney, NSW 2109, Australia
| | - Andrew Care
- Department of Molecular Sciences, Macquarie University, Sydney, NSW 2109, Australia; ARC Centre of Excellence for Nanoscale Biophotonics, Macquarie University, Sydney, NSW 2109, Australia
| | - Megan S Lord
- Graduate School of Biomedical Engineering, University of New South Wales, Sydney, NSW 2052, Australia
| | - Tiffany R Walsh
- Institute for Frontier Materials, Deakin University, Geelong, Victoria 3216, Australia
| | - Anwar Sunna
- Department of Molecular Sciences, Macquarie University, Sydney, NSW 2109, Australia; ARC Centre of Excellence for Nanoscale Biophotonics, Macquarie University, Sydney, NSW 2109, Australia; Biomolecular Discovery and Design Research Centre, Macquarie University, Sydney, NSW 2109, Australia.
| |
Collapse
|
21
|
Sharp ME, Vázquez FX, Wagner JW, Dannenhoffer-Lafage T, Voth GA. Multiconfigurational Coarse-Grained Molecular Dynamics. J Chem Theory Comput 2019; 15:3306-3315. [PMID: 30897328 PMCID: PMC6660024 DOI: 10.1021/acs.jctc.8b01133] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
![]()
Standard low resolution
coarse-grained modeling techniques have difficulty capturing multiple
configurations of protein systems. Here, we present a method for creating
accurate coarse-grained (CG) models with multiple configurations using
a linear combination of functions or “states”. Individual
CG models are created to capture the individual states, and the approximate
coupling between the two states is determined from an all-atom potential
of mean force. We show that the resulting multiconfiguration coarse-graining
(MCCG) method accurately captures the transition state as well as
the free energy between the two states. We have tested this method
on the folding of dodecaalanine, as well as the amphipathic helix
of endophilin.
Collapse
Affiliation(s)
- Morris E Sharp
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Francisco X Vázquez
- Department of Chemistry , St. John's University , Queens , New York 11439 , United States
| | - Jacob W Wagner
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Thomas Dannenhoffer-Lafage
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
22
|
Pahal V. Quercetin, a secondary metabolite present in methanolic extract of Calendula officinalis, is a potent inhibitor of peptide deformylase, undecaprenyl pyrophosphate synthase and DNA primase enzymes of Staphylococcus aureus: an in vitro and in silico result analysis. ACTA ACUST UNITED AC 2018. [DOI: 10.15406/mojddt.2018.02.00050] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
23
|
Yuan X, Xu Y. Recent Trends and Applications of Molecular Modeling in GPCR⁻Ligand Recognition and Structure-Based Drug Design. Int J Mol Sci 2018; 19:ijms19072105. [PMID: 30036949 PMCID: PMC6073596 DOI: 10.3390/ijms19072105] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Revised: 07/12/2018] [Accepted: 07/12/2018] [Indexed: 01/14/2023] Open
Abstract
G protein-coupled receptors represent the largest family of human membrane proteins and are modulated by a variety of drugs and endogenous ligands. Molecular modeling techniques, especially enhanced sampling methods, have provided significant insight into the mechanism of GPCR–ligand recognition. Notably, the crucial role of the membrane in the ligand-receptor association process has earned much attention. Additionally, docking, together with more accurate free energy calculation methods, is playing an important role in the design of novel compounds targeting GPCRs. Here, we summarize the recent progress in the computational studies focusing on the above issues. In the future, with continuous improvement in both computational hardware and algorithms, molecular modeling would serve as an indispensable tool in a wider scope of the research concerning GPCR–ligand recognition as well as drug design targeting GPCRs.
Collapse
Affiliation(s)
- Xiaojing Yuan
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences (CAS), Shanghai 201203, China.
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing 100049, China.
| | - Yechun Xu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences (CAS), Shanghai 201203, China.
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing 100049, China.
| |
Collapse
|
24
|
Gupta A, Chaudhary N, Aparoy P. MM-PBSA and per-residue decomposition energy studies on 7-Phenyl-imidazoquinolin-4(5H)-one derivatives: Identification of crucial site points at microsomal prostaglandin E synthase-1 (mPGES-1) active site. Int J Biol Macromol 2018; 119:352-359. [PMID: 30031079 DOI: 10.1016/j.ijbiomac.2018.07.050] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Revised: 07/10/2018] [Accepted: 07/11/2018] [Indexed: 12/29/2022]
Abstract
The huge therapeutic potential and the market share of painkillers are well-known. Due to the side effects associated with traditional NSAIDs and selective cyclooxygenase (COX-2) inhibitors, a new generation of painkillers is the need of the hour. In this regard, microsomal prostaglandin E synthase-1 (mPGES-1) offers great potential as an alternative drug target against inflammatory disorders. The present study is aimed at identifying the amino acids crucial in effective inhibitor binding at the mPGES-1 active site by performing molecular dynamics based studies on a series of 7-Phenyl-imidazoquinolin-4(5H)-one derivatives. Molecular dynamics (MD) simulations, MM-PBSA, per-residue energy decomposition and Dimensionality Reduction through Covariance matrix Transformation for Identification of Differences in dynamics (DIRECT-ID) analysis were performed to get insights into the structural details that can aid in novel drug design against mPGES-1. The high correlations of electrostatic and polar energy terms with biological activity highlight their importance and applicability in in silico screening studies. Further, per-residue energy decomposition studies revealed that Lys42, Arg52, Arg122, Pro124, Ser127, Val128 and Thr131 were contributing more towards inhibitor binding energy. The results clearly show that MM-PBSA can act as a filter in virtual screening experiments and can play major role in facilitating various mPGES-1 drug discovery studies.
Collapse
Affiliation(s)
- Ashish Gupta
- Centre for Computational Biology and Bioinformatics, School of Life Sciences, Central University of Himachal Pradesh, Dharamshala, Himachal Pradesh 176215, India
| | - Neha Chaudhary
- Centre for Computational Biology and Bioinformatics, School of Life Sciences, Central University of Himachal Pradesh, Dharamshala, Himachal Pradesh 176215, India
| | - Polamarasetty Aparoy
- Centre for Computational Biology and Bioinformatics, School of Life Sciences, Central University of Himachal Pradesh, Dharamshala, Himachal Pradesh 176215, India.
| |
Collapse
|
25
|
Abstract
The structure of RNA has been a natural subject for mathematical modeling, inviting many innovative computational frameworks. This single-stranded polynucleotide chain can fold upon itself in numerous ways to form hydrogen-bonded segments, imperfect with single-stranded loops. Illustrating these paired and non-paired interaction networks, known as RNA's secondary (2D) structure, using mathematical graph objects has been illuminating for RNA structure analysis. Building upon such seminal work from the 1970s and 1980s, graph models are now used to study not only RNA structure but also describe RNA's recurring modular units, sample the conformational space accessible to RNAs, predict RNA's three-dimensional folds, and apply the combined aspects to novel RNA design. In this article, we outline the development of the RNA-As-Graphs (or RAG) approach and highlight current applications to RNA structure prediction and design.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, 100 Washington Square East, Silver Building, New York University, New York, NY 10003, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer St., New York, NY 10012, USA; New York University ECNU - Center for Computational Chemistry at NYU Shanghai, 3663 North Zhongshan Road, Shanghai, 200062, China.
| |
Collapse
|
26
|
Bacci M, Langini C, Vymětal J, Caflisch A, Vitalis A. Focused conformational sampling in proteins. J Chem Phys 2018; 147:195102. [PMID: 29166086 DOI: 10.1063/1.4996879] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.
Collapse
Affiliation(s)
- Marco Bacci
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Cassiano Langini
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Jiří Vymětal
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Andreas Vitalis
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| |
Collapse
|
27
|
Thurston BA, Ferguson AL. Machine learning and molecular design of self-assembling -conjugated oligopeptides. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1469754] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Bryce A. Thurston
- Department of Physics, University of Illinois at Urbana-Champaign , Urbana, IL, USA
| | - Andrew L. Ferguson
- Department of Physics, University of Illinois at Urbana-Champaign , Urbana, IL, USA
- Department of Materials Science and Engineering, University of Illinois at Urbana-Champaign , Urbana, IL, USA
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign , Urbana, IL, USA
| |
Collapse
|
28
|
Rungsung I, Ramaswamy A. Effects of Peutz-Jeghers syndrome (PJS) causing missense mutations L67P, L182P, G242V and R297S on the structural dynamics of LKB1 (Liver kinase B1) protein. J Biomol Struct Dyn 2018; 37:796-810. [PMID: 29447078 DOI: 10.1080/07391102.2018.1441070] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The liver kinase B1 (LKB1) is encoded by LKB1 gene. Several pathogenic mutations of LKB1 causing Peutz-Jeghers syndrome and also cancers in breast, gastric, pancreas, and colon have been reported. The present study is focused to analyze the effects on the structural dynamics of LKB1 caused by the 4 pathogenic missense mutations (L67P, L182P, G242V, and R297S), which are reported to reduce the catalytic activity. In this study, the structural changes of LKB1 in apo- and in heterotrimeric complex (LKB1-STRADα-MO25α) form with wild and mutated LKB1 are investigated using all atomistic molecular dynamic simulation. The present study reveals that these four mutations initiate local structural distortions and the solvent accessibility of the surrounding regions of ATP-binding pocket such as glycine-rich loop, αB and αC loop, activation and catalytic loops. The mutations of L67P, L182P, and G242 V induce distortions of the secondary structure of β1-β3 sheets, π - π interaction (observed between Phe204 of LKB1 and Phe243 of MO25α), and increase the helical properties (both helical twist and length) of the adjacent αH-helix, respectively. The active kinase features like the conformation of catalytic and activation loops, salt bridge and, finally, the formation of stable R- and C-hydrophobic spines are also found to be perturbed by these mutations. Hence, the observed mutation-induced structural distortions fail to coordinate the essential binding nature of LKB1 with STRADα and MO25α, which eventually affects the native function of LKB1. These observations are in line with the experimentally reported reduced kinase activity of LKB1.
Collapse
Affiliation(s)
- Ikrormi Rungsung
- a Centre for Bioinformatics, School of Life Sciences , Pondicherry University , Puducherry 605014 , India
| | - Amutha Ramaswamy
- a Centre for Bioinformatics, School of Life Sciences , Pondicherry University , Puducherry 605014 , India
| |
Collapse
|
29
|
|
30
|
|
31
|
Onufriev AV, Izadi S. Water models for biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1347] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Alexey V. Onufriev
- Department of Physics; Virginia Tech; Blacksburg VA USA
- Department of Computer Science; Virginia Tech; Blacksburg VA USA
- Center for Soft Matter and Biological Physics; Virginia Tech; Blacksburg VA USA
| | - Saeed Izadi
- Early Stage Pharmaceutical Development; Genentech Inc.; South San Francisco, CA USA
| |
Collapse
|
32
|
Ramezanghorbani F, Dalgicdir C, Sayar M. A multi-state coarse grained modeling approach for an intrinsically disordered peptide. J Chem Phys 2017; 147:094103. [DOI: 10.1063/1.5001087] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
33
|
Bope CD, Tong D, Li X, Lu L. Fluctuation matching approach for elastic network model and structure-based model of biomacromolecules. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2017; 128:100-112. [DOI: 10.1016/j.pbiomolbio.2016.12.006] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2016] [Revised: 11/25/2016] [Accepted: 12/19/2016] [Indexed: 10/24/2022]
|
34
|
HÖft TA, Alpert BK. FAST UPDATING MULTIPOLE COULOMBIC POTENTIAL CALCULATION. SIAM JOURNAL ON SCIENTIFIC COMPUTING : A PUBLICATION OF THE SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS 2017; 39:https://doi.org/10.1137/16M1096189. [PMID: 33088167 PMCID: PMC7574401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present a numerical method to efficiently and accurately re-compute the Coulomb potential of a large ensemble of charged particles after a subset of the particles undergoes a change of position. Errors are bounded even after a large number of such shifts, making it practical for use in Monte Carlo Markov chain methods in molecular dynamics, computational astrophysics, computational chemistry, and other applications. The method uses truncated multipole expansions of the potential energy functional and a tree decomposition of the computational domain to reduce the computational complexity. Computational costs scale logarithmically in the size of the problem. Scaling, accuracy, and efficiency are confirmed with numerical experiments. The new method outperforms a direct calculation for moderate problem sizes.
Collapse
Affiliation(s)
- Thomas A HÖft
- Department of Mathematics, University of St. Thomas, Saint Paul, MN 55105
| | - Bradley K Alpert
- National Institute of Standards and Technology, Boulder, Colorado 80305
| |
Collapse
|
35
|
Rudzinski JF, Lu K, Milner ST, Maranas JK, Noid WG. Extended Ensemble Approach to Transferable Potentials for Low-Resolution Coarse-Grained Models of Ionomers. J Chem Theory Comput 2017; 13:2185-2201. [DOI: 10.1021/acs.jctc.6b01160] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joseph F. Rudzinski
- Department
of Chemistry and ‡Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Keran Lu
- Department
of Chemistry and ‡Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Scott T. Milner
- Department
of Chemistry and ‡Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Janna K. Maranas
- Department
of Chemistry and ‡Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - William G. Noid
- Department
of Chemistry and ‡Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
36
|
Effect of solvent model when probing protein dynamics with molecular dynamics. J Mol Graph Model 2017; 71:80-87. [DOI: 10.1016/j.jmgm.2016.11.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2016] [Revised: 11/02/2016] [Accepted: 11/03/2016] [Indexed: 12/17/2022]
|
37
|
Tang Z, Roberts CC, Chang CEA. Understanding ligand-receptor non-covalent binding kinetics using molecular modeling. FRONT BIOSCI-LANDMRK 2017; 22:960-981. [PMID: 27814657 PMCID: PMC5470370 DOI: 10.2741/4527] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Kinetic properties may serve as critical differentiators and predictors of drug efficacy and safety, in addition to the traditionally focused binding affinity. However the quantitative structure-kinetics relationship (QSKR) for modeling and ligand design is still poorly understood. This review provides an introduction to the kinetics of drug binding from a fundamental chemistry perspective. We focus on recent developments of computational tools and their applications to non-covalent binding kinetics.
Collapse
Affiliation(s)
- Zhiye Tang
- Department of Chemistry, University of California, Riverside, CA 92521
| | | | - Chia-En A Chang
- Department of Chemistry, University of California, Riverside, CA 92521,
| |
Collapse
|
38
|
Rychkov GN, Ilatovskiy AV, Nazarov IB, Shvetsov AV, Lebedev DV, Konev AY, Isaev-Ivanov VV, Onufriev AV. Partially Assembled Nucleosome Structures at Atomic Detail. Biophys J 2016; 112:460-472. [PMID: 28038734 DOI: 10.1016/j.bpj.2016.10.041] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2016] [Revised: 10/06/2016] [Accepted: 10/28/2016] [Indexed: 11/29/2022] Open
Abstract
The evidence is now overwhelming that partially assembled nucleosome states (PANS) are as important as the canonical nucleosome structure for the understanding of how accessibility to genomic DNA is regulated in cells. We use a combination of molecular dynamics simulation and atomic force microscopy to deliver, in atomic detail, structural models of three key PANS: the hexasome (H2A·H2B)·(H3·H4)2, the tetrasome (H3·H4)2, and the disome (H3·H4). Despite fluctuations of the conformation of the free DNA in these structures, regions of protected DNA in close contact with the histone core remain stable, thus establishing the basis for the understanding of the role of PANS in DNA accessibility regulation. On average, the length of protected DNA in each structure is roughly 18 basepairs per histone protein. Atomistically detailed PANS are used to explain experimental observations; specifically, we discuss interpretation of atomic force microscopy, Förster resonance energy transfer, and small-angle x-ray scattering data obtained under conditions when PANS are expected to exist. Further, we suggest an alternative interpretation of a recent genome-wide study of DNA protection in active chromatin of fruit fly, leading to a conclusion that the three PANS are present in actively transcribing regions in a substantial amount. The presence of PANS may not only be a consequence, but also a prerequisite for fast transcription in vivo.
Collapse
Affiliation(s)
- Georgy N Rychkov
- Division of Molecular and Radiation Biophysics, B.P. Konstantinov Petersburg Nuclear Physics Institute, National Research Center "Kurchatov Institute", Orlova Roscha, Gatchina, Leningrad District, Russia; Institute of Physics, Nanotechnology and Telecommunications, NRU Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia
| | - Andrey V Ilatovskiy
- Division of Molecular and Radiation Biophysics, B.P. Konstantinov Petersburg Nuclear Physics Institute, National Research Center "Kurchatov Institute", Orlova Roscha, Gatchina, Leningrad District, Russia; Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, California
| | - Igor B Nazarov
- Institute of Cytology, Russian Academy of Sciences, St. Petersburg, Russia
| | - Alexey V Shvetsov
- Division of Molecular and Radiation Biophysics, B.P. Konstantinov Petersburg Nuclear Physics Institute, National Research Center "Kurchatov Institute", Orlova Roscha, Gatchina, Leningrad District, Russia; Institute of Applied Mathematics and Mechanics, NRU Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia
| | - Dmitry V Lebedev
- Division of Molecular and Radiation Biophysics, B.P. Konstantinov Petersburg Nuclear Physics Institute, National Research Center "Kurchatov Institute", Orlova Roscha, Gatchina, Leningrad District, Russia
| | - Alexander Y Konev
- Division of Molecular and Radiation Biophysics, B.P. Konstantinov Petersburg Nuclear Physics Institute, National Research Center "Kurchatov Institute", Orlova Roscha, Gatchina, Leningrad District, Russia
| | - Vladimir V Isaev-Ivanov
- Division of Molecular and Radiation Biophysics, B.P. Konstantinov Petersburg Nuclear Physics Institute, National Research Center "Kurchatov Institute", Orlova Roscha, Gatchina, Leningrad District, Russia
| | - Alexey V Onufriev
- Departments of Computer Science and Physics, Virginia Tech, Blacksburg, Virginia.
| |
Collapse
|
39
|
Yamamoto S, Nakazawa S, Sugisaki K, Maekawa K, Sato K, Toyota K, Shiomi D, Takui T. Structural Determination of a DNA Oligomer for a Molecular Spin Qubit Lloyd Model of Quantum Computers. Z PHYS CHEM 2016. [DOI: 10.1515/zpch-2016-0799] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Abstract
The global molecular and local spin-site structures of a DNA duplex 22-oligomer with site-directed four spin-labeling were simulated by molecular mechanics (MM) calculations combined with Q-band pulsed electron-electron double resonance (PELDOR) spectroscopy. This molecular-spin bearing DNA oligomer is designed to give a complex testing ground for the structural determination of molecular spins incorporated in the DNA duplex, which serves as a platform for 1D periodic arrays of two or three non-equivalent electron spin qubit systems, (AB)n or (ABC)n, respectively, enabling to execute quantum computing or quantum information processing (Lloyd model of electron spin versions): A, B and C designate non-equivalent addressable spin qubits for quantum operations. The non-equivalence originates in difference in the electronic g-tensor. It is not feasible to determine the optimal structures for such DNA oligomers having molecular flexibility only by the MM calculations because there are many local minima in energy for their possible molecular structures. The spin-distance information derived from the PELDOR spectroscopy helps determine the optimal structures out of the possible ones acquired by the MM calculations. Based on the MM searched structures, we suggest the optimal structures for semi-macromolecules having site-directed multi-spin qubits. We emphasize that for our four molecular spins embedded in the DNA oligomer the Fajer’s error analysis in PELDOR-based distance measurements was of essential importance.
Collapse
Affiliation(s)
- Satoru Yamamoto
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
| | - Shigeaki Nakazawa
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
- FIRST Project on “Quantum Information Processing”, The Cabinet Office, JSPS, Tokyo 101-8430, Japan
| | - Kenji Sugisaki
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
- FIRST Project on “Quantum Information Processing”, The Cabinet Office, JSPS, Tokyo 101-8430, Japan
| | - Kensuke Maekawa
- Department of Regulatory Bioorganic Chemistry, The Institute of Scientific Industrial Research (ISIR), Osaka University, Ibaraki 567-0047, Japan
| | - Kazunobu Sato
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
- FIRST Project on “Quantum Information Processing”, The Cabinet Office, JSPS, Tokyo 101-8430, Japan , Phone: +81-6605-2605, Fax: +81-6605-2522
| | - Kazuo Toyota
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
- FIRST Project on “Quantum Information Processing”, The Cabinet Office, JSPS, Tokyo 101-8430, Japan
| | - Daisuke Shiomi
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
- FIRST Project on “Quantum Information Processing”, The Cabinet Office, JSPS, Tokyo 101-8430, Japan
| | - Takeji Takui
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138, Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
- FIRST Project on “Quantum Information Processing”, The Cabinet Office, JSPS, Tokyo 101-8430, Japan , Phone: +81-6605-2605, Fax: +81-6605-2522
| |
Collapse
|
40
|
Izadi S, Onufriev AV. Accuracy limit of rigid 3-point water models. J Chem Phys 2016; 145:074501. [PMID: 27544113 PMCID: PMC4991989 DOI: 10.1063/1.4960175] [Citation(s) in RCA: 169] [Impact Index Per Article: 21.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2016] [Accepted: 07/19/2016] [Indexed: 11/14/2022] Open
Abstract
Classical 3-point rigid water models are most widely used due to their computational efficiency. Recently, we introduced a new approach to constructing classical rigid water models [S. Izadi et al., J. Phys. Chem. Lett. 5, 3863 (2014)], which permits a virtually exhaustive search for globally optimal model parameters in the sub-space that is most relevant to the electrostatic properties of the water molecule in liquid phase. Here we apply the approach to develop a 3-point Optimal Point Charge (OPC3) water model. OPC3 is significantly more accurate than the commonly used water models of same class (TIP3P and SPCE) in reproducing a comprehensive set of liquid bulk properties, over a wide range of temperatures. Beyond bulk properties, we show that OPC3 predicts the intrinsic charge hydration asymmetry (CHA) of water - a characteristic dependence of hydration free energy on the sign of the solute charge - in very close agreement with experiment. Two other recent 3-point rigid water models, TIP3PFB and H2ODC, each developed by its own, completely different optimization method, approach the global accuracy optimum represented by OPC3 in both the parameter space and accuracy of bulk properties. Thus, we argue that an accuracy limit of practical 3-point rigid non-polarizable models has effectively been reached; remaining accuracy issues are discussed.
Collapse
Affiliation(s)
- Saeed Izadi
- Department of Biomedical Engineering and Mechanics, Virginia Tech, Blacksburg, Virginia 24060, USA
| | - Alexey V Onufriev
- Departments of Computer Science and Physics, Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24060, USA
| |
Collapse
|
41
|
All-atom/coarse-grained hybrid predictions of distribution coefficients in SAMPL5. J Comput Aided Mol Des 2016; 30:969-976. [PMID: 27460060 PMCID: PMC5206257 DOI: 10.1007/s10822-016-9926-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2016] [Accepted: 07/16/2016] [Indexed: 12/14/2022]
Abstract
We present blind predictions submitted to the SAMPL5 challenge on calculating distribution coefficients. The predictions were based on estimating the solvation free energies in water and cyclohexane of the 53 compounds in the challenge. These free energies were computed using alchemical free energy simulations based on a hybrid all-atom/coarse-grained model. The compounds were treated with the general Amber force field, whereas the solvent molecules were treated with the Elba coarse-grained model. Considering the simplicity of the solvent model and that we approximate the distribution coefficient with the partition coefficient of the neutral species, the predictions are of good accuracy. The correlation coefficient, R is 0.64, 82 % of the predictions have the correct sign and the mean absolute deviation is 1.8 log units. This is on a par with or better than the other simulation-based predictions in the challenge. We present an analysis of the deviations to experiments and compare the predictions to another submission that used all-atom solvent.
Collapse
|
42
|
|
43
|
Kim T, Freudenthal BD, Beard WA, Wilson SH, Schlick T. Insertion of oxidized nucleotide triggers rapid DNA polymerase opening. Nucleic Acids Res 2016; 44:4409-24. [PMID: 27034465 PMCID: PMC4872097 DOI: 10.1093/nar/gkw174] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2016] [Accepted: 03/04/2016] [Indexed: 11/21/2022] Open
Abstract
A novel mechanism is unveiled to explain why a pro-mutagenic nucleotide lesion (oxidized guanine, 8-oxoG) causes the mammalian DNA repair polymerase-β (pol-β) to rapidly transition to an inactive open conformation. The mechanism involves unexpected features revealed recently in time-lapse crystallography. Specifically, a delicate water network associated with a lesion-stabilizing auxilliary product ion Mg(p) triggers a cascade of events that leads to poor active site geometry and the rupture of crucial molecular interactions between key residues in both the anti(8-oxoG:C) and syn(8-oxoG:A) systems. Once the base pairs in these lesioned systems are broken, dislocation of both Asp192 (a metal coordinating ligand) and the oxoG phosphate group (PO4) interfere with the hydrogen bonding between Asp192 and Arg258, whose rotation toward Asp192 is crucial to the closed-to-open enzyme transition. Energetically, the lesioned open states are similar in energy to those of the corresponding closed complexes after chemistry, in marked contrast to the unlesioned pol-β anti(G:C) system, whose open state is energetically higher than the closed state. The delicate surveillance system offers a fundamental protective mechanism in the cell that triggers DNA repair events which help deter insertion of oxidized lesions.
Collapse
Affiliation(s)
- Taejin Kim
- Department of Chemistry, New York University, 10th Floor Silver Center, 100 Washington Square East, New York, NY 10003, USA
| | - Bret D Freudenthal
- Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709, USA
| | - William A Beard
- Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709, USA
| | - Samuel H Wilson
- Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709, USA
| | - Tamar Schlick
- Department of Chemistry, New York University, 10th Floor Silver Center, 100 Washington Square East, New York, NY 10003, USA Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA
| |
Collapse
|
44
|
Rodrigues JPGLM, Melquiond ASJ, Bonvin AMJJ. Molecular dynamics characterization of the conformational landscape of small peptides: A series of hands-on collaborative practical sessions for undergraduate students. BIOCHEMISTRY AND MOLECULAR BIOLOGY EDUCATION : A BIMONTHLY PUBLICATION OF THE INTERNATIONAL UNION OF BIOCHEMISTRY AND MOLECULAR BIOLOGY 2016; 44:160-167. [PMID: 26751257 DOI: 10.1002/bmb.20941] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2015] [Accepted: 11/15/2015] [Indexed: 06/05/2023]
Abstract
Molecular modelling and simulations are nowadays an integral part of research in areas ranging from physics to chemistry to structural biology, as well as pharmaceutical drug design. This popularity is due to the development of high-performance hardware and of accurate and efficient molecular mechanics algorithms by the scientific community. These improvements are also benefitting scientific education. Molecular simulations, their underlying theory, and their applications are particularly difficult to grasp for undergraduate students. Having hands-on experience with the methods contributes to a better understanding and solidification of the concepts taught during the lectures. To this end, we have created a computer practical class, which has been running for the past five years, composed of several sessions where students characterize the conformational landscape of small peptides using molecular dynamics simulations in order to gain insights on their binding to protein receptors. In this report, we detail the ingredients and recipe necessary to establish and carry out this practical, as well as some of the questions posed to the students and their expected results. Further, we cite some examples of the students' written reports, provide statistics, and share their feedbacks on the structure and execution of the sessions. These sessions were implemented alongside a theoretical molecular modelling course but have also been used successfully as a standalone tutorial during specialized workshops. The availability of the material on our web page also facilitates this integration and dissemination and lends strength to the thesis of open-source science and education.
Collapse
Affiliation(s)
- João P G L M Rodrigues
- Computational Structural Biology Group, Bijvoet Center for Biomolecular Research, Faculty of Science-Chemistry, Utrecht University, Utrecht, The Netherlands
| | - Adrien S J Melquiond
- Biomedical Genomics, Hubrecht Insitute, KNAW & University Medical Center Utrecht, Utrecht, The Netherlands
| | - Alexandre M J J Bonvin
- Computational Structural Biology Group, Bijvoet Center for Biomolecular Research, Faculty of Science-Chemistry, Utrecht University, Utrecht, The Netherlands
| |
Collapse
|
45
|
Comparing pharmacophore models derived from crystal structures and from molecular dynamics simulations. MONATSHEFTE FUR CHEMIE 2016; 147:553-563. [PMID: 27069282 PMCID: PMC4785218 DOI: 10.1007/s00706-016-1674-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/15/2015] [Accepted: 01/14/2016] [Indexed: 01/23/2023]
Abstract
ABSTRACT Pharmacophore modeling is a widely used technique in computer-aided drug discovery. Structure-based pharmacophore models of a ligand in complex with a protein have proven to be useful for supporting in silico hit discovery, hit to lead expansion, and lead optimization. As a structure-based approach it depends on the correct interpretation of ligand-protein interactions. There are legitimate concerns about the fidelity of the bound ligand and about non-physiological contacts with parts of the crystal and the solvent effects that influence the protein structure. A possible way to refine the structure of a protein-ligand system is to use the final structure of a given MD simulation. In this study we compare pharmacophore models built using the initial protein-ligand structure obtained from the protein data bank (PDB) with pharmacophore models built with the final structure of a molecular dynamics simulation. We show that the pharmacophore models differ in feature number and feature type and that the pharmacophore models built from the last structure of a MD simulation shows in some cases better ability to distinguish between active and decoy ligand structures. GRAPHICAL ABSTRACT
Collapse
|
46
|
Anandakrishnan R, Drozdetski A, Walker RC, Onufriev AV. Speed of conformational change: comparing explicit and implicit solvent molecular dynamics simulations. Biophys J 2016; 108:1153-64. [PMID: 25762327 DOI: 10.1016/j.bpj.2014.12.047] [Citation(s) in RCA: 130] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Revised: 12/23/2014] [Accepted: 12/29/2014] [Indexed: 11/24/2022] Open
Abstract
Adequate sampling of conformation space remains challenging in atomistic simulations, especially if the solvent is treated explicitly. Implicit-solvent simulations can speed up conformational sampling significantly. We compare the speed of conformational sampling between two commonly used methods of each class: the explicit-solvent particle mesh Ewald (PME) with TIP3P water model and a popular generalized Born (GB) implicit-solvent model, as implemented in the AMBER package. We systematically investigate small (dihedral angle flips in a protein), large (nucleosome tail collapse and DNA unwrapping), and mixed (folding of a miniprotein) conformational changes, with nominal simulation times ranging from nanoseconds to microseconds depending on system size. The speedups in conformational sampling for GB relative to PME simulations, are highly system- and problem-dependent. Where the simulation temperatures for PME and GB are the same, the corresponding speedups are approximately onefold (small conformational changes), between ∼1- and ∼100-fold (large changes), and approximately sevenfold (mixed case). The effects of temperature on speedup and free-energy landscapes, which may differ substantially between the solvent models, are discussed in detail for the case of miniprotein folding. In addition to speeding up conformational sampling, due to algorithmic differences, the implicit solvent model can be computationally faster for small systems or slower for large systems, depending on the number of solute and solvent atoms. For the conformational changes considered here, the combined speedups are approximately twofold, ∼1- to 60-fold, and ∼50-fold, respectively, in the low solvent viscosity regime afforded by the implicit solvent. For all the systems studied, 1) conformational sampling speedup increases as Langevin collision frequency (effective viscosity) decreases; and 2) conformational sampling speedup is mainly due to reduction in solvent viscosity rather than possible differences in free-energy landscapes between the solvent models.
Collapse
Affiliation(s)
| | | | - Ross C Walker
- San Diego Supercomputer Center and Department of Chemistry and Biochemistry, University of California, San Diego, San Diego, California
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia; Department of Physics, Virginia Tech, Blacksburg, Virginia.
| |
Collapse
|
47
|
Dunn NJH, Noid WG. Bottom-up coarse-grained models that accurately describe the structure, pressure, and compressibility of molecular liquids. J Chem Phys 2015; 143:243148. [DOI: 10.1063/1.4937383] [Citation(s) in RCA: 76] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Nicholas J. H. Dunn
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - W. G. Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
48
|
Boniecki MJ, Lach G, Dawson WK, Tomala K, Lukasz P, Soltysinski T, Rother KM, Bujnicki JM. SimRNA: a coarse-grained method for RNA folding simulations and 3D structure prediction. Nucleic Acids Res 2015; 44:e63. [PMID: 26687716 PMCID: PMC4838351 DOI: 10.1093/nar/gkv1479] [Citation(s) in RCA: 236] [Impact Index Per Article: 26.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2014] [Accepted: 12/05/2015] [Indexed: 01/08/2023] Open
Abstract
RNA molecules play fundamental roles in cellular processes. Their function and interactions with other biomolecules are dependent on the ability to form complex three-dimensional (3D) structures. However, experimental determination of RNA 3D structures is laborious and challenging, and therefore, the majority of known RNAs remain structurally uncharacterized. Here, we present SimRNA: a new method for computational RNA 3D structure prediction, which uses a coarse-grained representation, relies on the Monte Carlo method for sampling the conformational space, and employs a statistical potential to approximate the energy and identify conformations that correspond to biologically relevant structures. SimRNA can fold RNA molecules using only sequence information, and, on established test sequences, it recapitulates secondary structure with high accuracy, including correct prediction of pseudoknots. For modeling of complex 3D structures, it can use additional restraints, derived from experimental or computational analyses, including information about secondary structure and/or long-range contacts. SimRNA also can be used to analyze conformational landscapes and identify potential alternative structures.
Collapse
Affiliation(s)
- Michal J Boniecki
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Grzegorz Lach
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Wayne K Dawson
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Konrad Tomala
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Pawel Lukasz
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Tomasz Soltysinski
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Kristian M Rother
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Janusz M Bujnicki
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| |
Collapse
|
49
|
Ozer G, Luque A, Schlick T. The chromatin fiber: multiscale problems and approaches. Curr Opin Struct Biol 2015; 31:124-39. [PMID: 26057099 DOI: 10.1016/j.sbi.2015.04.002] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Revised: 04/08/2015] [Accepted: 04/09/2015] [Indexed: 12/20/2022]
Abstract
The structure of chromatin, affected by many factors from DNA linker lengths to posttranslational modifications, is crucial to the regulation of eukaryotic cells. Combined experimental and computational methods have led to new insights into its structural and dynamical features, from interactions due to the flexible core histone tails or linker histones to the physical mechanism driving the formation of chromosomal domains. Here we present a perspective of recent advances in chromatin modeling techniques at the atomic, mesoscopic, and chromosomal scales with a view toward developing multiscale computational strategies to integrate such findings. Innovative modeling methods that connect molecular to chromosomal scales are crucial for interpreting experiments and eventually deciphering the complex dynamic organization and function of chromatin in the cell.
Collapse
Affiliation(s)
- Gungor Ozer
- Department of Chemistry, 100 Washington Square East, New York University, New York, NY 10003, USA
| | - Antoni Luque
- Department of Chemistry, 100 Washington Square East, New York University, New York, NY 10003, USA; Current address: Department of Mathematics & Statistics and Viral Information Institute, San Diego State University, 5500 Campanile Drive, San Diego, CA 92182-7720, USA
| | - Tamar Schlick
- Department of Chemistry, 100 Washington Square East, New York University, New York, NY 10003, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA; NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China.
| |
Collapse
|
50
|
Chandramouli B, Zazza C, Mancini G, Brancato G. Boundary condition effects on the dynamic and electric properties of hydration layers. J Phys Chem A 2015; 119:5465-75. [PMID: 25752804 DOI: 10.1021/jp511824t] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Water solvation has a central role in several biochemical processes ranging from protein folding to biomolecular recognition and enzyme catalysis. Because of its importance, the structure and dynamics of hydration layers around biological macromolecules have been the targets of a great number of experimental and computational studies. In the present contribution, we have investigated the effects of periodic boundary conditions (PBCs), as used in conjunction with molecular dynamics (MD) simulations, on the dynamic and electric properties of water layers. In particular, we have systematically performed MD simulations of neat water and biomolecules in aqueous solutions by imposing a different external dielectric constant, a generally overlooked parameter in PBC simulations. The effect of the system size has also been addressed. Overall, our results consistently indicate that the dipole moment properties of water layers, and specifically the dipole moment fluctuations and the reorientational correlation functions, can be sensitive to the choice of the external boundary conditions, whereas other molecular properties, such as the self-diffusion coefficient and the reorientational relaxation times, are not affected. We think that our investigation may help to assess appropriate simulation conditions for modeling the aqueous environment of relevant biochemical systems and processes.
Collapse
Affiliation(s)
| | - Costantino Zazza
- †Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Giordano Mancini
- †Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,‡Istituto Nazionale di Fisica Nucleare (INFN) Sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Giuseppe Brancato
- †Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,‡Istituto Nazionale di Fisica Nucleare (INFN) Sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|