1
|
Truong A, Barton M, Tran U, Mellody M, Berger D, Madory D, Hitch E, Jibrael B, Nikolaidis N, Luchko T, Keppetipola N. Unstructured linker regions play a role in the differential splicing activities of paralogous RNA binding proteins PTBP1 and PTBP2. J Biol Chem 2024; 300:105733. [PMID: 38336291 PMCID: PMC10914480 DOI: 10.1016/j.jbc.2024.105733] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Revised: 01/23/2024] [Accepted: 02/02/2024] [Indexed: 02/12/2024] Open
Abstract
RNA Binding Proteins regulate, in part, alternative pre-mRNA splicing and, in turn, gene expression patterns. Polypyrimidine tract binding proteins PTBP1 and PTBP2 are paralogous RNA binding proteins sharing 74% amino acid sequence identity. Both proteins contain four structured RNA-recognition motifs (RRMs) connected by linker regions and an N-terminal region. Despite their similarities, the paralogs have distinct tissue-specific expression patterns and can regulate discrete sets of target exons. How two highly structurally similar proteins can exert different splicing outcomes is not well understood. Previous studies revealed that PTBP2 is post-translationally phosphorylated in the unstructured N-terminal, Linker 1, and Linker 2 regions that share less sequence identity with PTBP1 signifying a role for these regions in dictating the paralog's distinct splicing activities. To this end, we conducted bioinformatics analysis to determine the evolutionary conservation of RRMs versus linker regions in PTBP1 and PTBP2 across species. To determine the role of PTBP2 unstructured regions in splicing activity, we created hybrid PTBP1-PTBP2 constructs that had counterpart PTBP1 regions swapped to an otherwise PTBP2 protein and assayed on differentially regulated exons. We also conducted molecular dynamics studies to investigate how negative charges introduced by phosphorylation in PTBP2 unstructured regions can alter their physical properties. Collectively, results from our studies reveal an important role for PTBP2 unstructured regions and suggest a role for phosphorylation in the differential splicing activities of the paralogs on certain regulated exons.
Collapse
Affiliation(s)
- Anthony Truong
- Department of Chemistry and Biochemistry, California State University Fullerton, Fullerton, California, USA
| | - Michael Barton
- Department of Physics and Astronomy, California State University, Northridge, Northridge, California, USA
| | - Uyenphuong Tran
- Department of Chemistry and Biochemistry, California State University Fullerton, Fullerton, California, USA
| | - Montana Mellody
- Department of Chemistry and Biochemistry, California State University Fullerton, Fullerton, California, USA
| | - Devon Berger
- Department of Biological Sciences, California State University Fullerton, Fullerton, California, USA
| | - Dean Madory
- Department of Biological Science, Santa Ana College, Santa Ana, California, USA
| | - Elizabeth Hitch
- Department of Biological Sciences, California State University Fullerton, Fullerton, California, USA
| | - Basma Jibrael
- Department of Chemistry and Biochemistry, California State University Fullerton, Fullerton, California, USA
| | - Nikolas Nikolaidis
- Department of Biological Sciences, California State University Fullerton, Fullerton, California, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Northridge, Northridge, California, USA.
| | - Niroshika Keppetipola
- Department of Chemistry and Biochemistry, California State University Fullerton, Fullerton, California, USA.
| |
Collapse
|
2
|
Case D, Aktulga HM, Belfon K, Cerutti DS, Cisneros GA, Cruzeiro VD, Forouzesh N, Giese TJ, Götz AW, Gohlke H, Izadi S, Kasavajhala K, Kaymak MC, King E, Kurtzman T, Lee TS, Li P, Liu J, Luchko T, Luo R, Manathunga M, Machado MR, Nguyen HM, O’Hearn KA, Onufriev AV, Pan F, Pantano S, Qi R, Rahnamoun A, Risheh A, Schott-Verdugo S, Shajan A, Swails J, Wang J, Wei H, Wu X, Wu Y, Zhang S, Zhao S, Zhu Q, Cheatham TE, Roe DR, Roitberg A, Simmerling C, York DM, Nagan MC, Merz KM. AmberTools. J Chem Inf Model 2023; 63:6183-6191. [PMID: 37805934 PMCID: PMC10598796 DOI: 10.1021/acs.jcim.3c01153] [Citation(s) in RCA: 42] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Indexed: 10/10/2023]
Abstract
AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.
Collapse
Affiliation(s)
- David
A. Case
- Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway 08854, New Jersey, United States
| | - Hasan Metin Aktulga
- Department
of Computer Science and Engineering, Michigan
State University, East Lansing 48824-1322, Michigan, United States
| | - Kellon Belfon
- FOG
Pharmaceuticals Inc., Cambridge 02140, Massachusetts, United States
| | - David S. Cerutti
- Psivant, 451 D Street, Suite 205, Boston 02210, Massachusetts, United States
| | - G. Andrés Cisneros
- Department
of Physics, Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson 75801, Texas, United States
| | - Vinícius
Wilian D. Cruzeiro
- Department
of Chemistry and The PULSE Institute, Stanford
University, Stanford 94305, California, United States
| | - Negin Forouzesh
- Department
of Computer Science, California State University, Los Angeles 90032, California, United States
| | - Timothy J. Giese
- Laboratory
for Biomolecular Simulation Research, Institute for Quantitative Biomedicine
and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway 08854, New Jersey, United States
| | - Andreas W. Götz
- San
Diego Supercomputer Center, University of
California San Diego, La Jolla 92093-0505, California, United States
| | - Holger Gohlke
- Institute
for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf 40225, Germany
- Institute
of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Jülich 52425, Germany
| | - Saeed Izadi
- Pharmaceutical
Development, Genentech, Inc., South San Francisco 94080, California, United
States
| | - Koushik Kasavajhala
- Laufer
Center for Physical and Quantitative Biology, Department of Chemistry, Stony Brook University, Stony Brook 11794, New York, United States
| | - Mehmet C. Kaymak
- Department
of Computer Science and Engineering, Michigan
State University, East Lansing 48824-1322, Michigan, United States
| | - Edward King
- Departments
of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering,
Materials Science and Engineering, and Biomedical Engineering, Graduate
Program in Chemical and Materials Physics, University of California, Irvine 92697, California, United States
| | - Tom Kurtzman
- Ph.D.
Programs in Chemistry, Biochemistry, and Biology, The Graduate Center of the City University of New York, 365 Fifth Avenue, New York 10016, New York, United States
- Department
of Chemistry, Lehman College, 250 Bedford Park Blvd West, Bronx 10468, New York, United States
| | - Tai-Sung Lee
- Laboratory
for Biomolecular Simulation Research, Institute for Quantitative Biomedicine
and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway 08854, New Jersey, United States
| | - Pengfei Li
- Department
of Chemistry and Biochemistry, Loyola University
Chicago, Chicago 60660, Illinois, United States
| | - Jian Liu
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Tyler Luchko
- Department
of Physics and Astronomy, California State
University, Northridge, Northridge 91330, California, United States
| | - Ray Luo
- Departments
of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering,
Materials Science and Engineering, and Biomedical Engineering, Graduate
Program in Chemical and Materials Physics, University of California, Irvine 92697, California, United States
| | - Madushanka Manathunga
- Department
of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing 48824-1322, Michigan, United States
| | | | - Hai Minh Nguyen
- Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway 08854, New Jersey, United States
| | - Kurt A. O’Hearn
- Department
of Computer Science and Engineering, Michigan
State University, East Lansing 48824-1322, Michigan, United States
| | - Alexey V. Onufriev
- Departments
of Computer Science and Physics, Virginia
Tech, Blacksburg 24061, Virginia, United
States
| | - Feng Pan
- Department
of Statistics, Florida State University, Tallahassee 32304, Florida, United States
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Ruxi Qi
- Cryo-EM
Center, Southern University of Science and
Technology, Shenzhen 518055, China
| | - Ali Rahnamoun
- Department
of Computer Science and Engineering, Michigan
State University, East Lansing 48824-1322, Michigan, United States
| | - Ali Risheh
- Department
of Computer Science, California State University, Los Angeles 90032, California, United States
| | - Stephan Schott-Verdugo
- Institute
of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, Jülich 52425, Germany
| | - Akhil Shajan
- Department
of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing 48824-1322, Michigan, United States
| | - Jason Swails
- Entos, 4470 W Sunset
Blvd, Suite 107, Los Angeles 90027, California, United States
| | - Junmei Wang
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh 15261, Pennsylvania, United States
| | - Haixin Wei
- Departments
of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering,
Materials Science and Engineering, and Biomedical Engineering, Graduate
Program in Chemical and Materials Physics, University of California, Irvine 92697, California, United States
| | - Xiongwu Wu
- Laboratory
of Computational Biology, NHLBI, NIH, Bethesda 20892, Maryland, United States
| | - Yongxian Wu
- Departments
of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering,
Materials Science and Engineering, and Biomedical Engineering, Graduate
Program in Chemical and Materials Physics, University of California, Irvine 92697, California, United States
| | - Shi Zhang
- Laboratory
for Biomolecular Simulation Research, Institute for Quantitative Biomedicine
and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway 08854, New Jersey, United States
| | - Shiji Zhao
- Departments
of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering,
Materials Science and Engineering, and Biomedical Engineering, Graduate
Program in Chemical and Materials Physics, University of California, Irvine 92697, California, United States
- Nurix Therapeutics, Inc., San Francisco 94158, California, United States
| | - Qiang Zhu
- Departments
of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering,
Materials Science and Engineering, and Biomedical Engineering, Graduate
Program in Chemical and Materials Physics, University of California, Irvine 92697, California, United States
| | - Thomas E. Cheatham
- Department
of Medicinal Chemistry, The University of
Utah, 30 South 2000 East, Salt Lake City 84112, Utah, United
States
| | - Daniel R. Roe
- Laboratory
of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda 20892, Maryland, United States
| | - Adrian Roitberg
- Department
of Chemistry, The University of Florida, 440 Leigh Hall, Gainesville 32611-7200, Florida, United States
| | - Carlos Simmerling
- Laufer
Center for Physical and Quantitative Biology, Department of Chemistry, Stony Brook University, Stony Brook 11794, New York, United States
| | - Darrin M. York
- Laboratory
for Biomolecular Simulation Research, Institute for Quantitative Biomedicine
and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway 08854, New Jersey, United States
| | - Maria C. Nagan
- Department
of Chemistry, Stony Brook University, Stony Brook 11794, New York, United States
| | - Kenneth M. Merz
- Department
of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing 48824-1322, Michigan, United States
| |
Collapse
|
3
|
Casillas L, Grigorian VM, Luchko T. Identifying Systematic Force Field Errors Using a 3D-RISM Element Counting Correction. Molecules 2023; 28:molecules28030925. [PMID: 36770599 PMCID: PMC9921782 DOI: 10.3390/molecules28030925] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Revised: 01/09/2023] [Accepted: 01/11/2023] [Indexed: 01/19/2023] Open
Abstract
Hydration free energies of small molecules are commonly used as benchmarks for solvation models. However, errors in predicting hydration free energies are partially due to the force fields used and not just the solvation model. To address this, we have used the 3D reference interaction site model (3D-RISM) of molecular solvation and existing benchmark explicit solvent calculations with a simple element count correction (ECC) to identify problems with the non-bond parameters in the general AMBER force field (GAFF). 3D-RISM was used to calculate hydration free energies of all 642 molecules in the FreeSolv database, and a partial molar volume correction (PMVC), ECC, and their combination (PMVECC) were applied to the results. The PMVECC produced a mean unsigned error of 1.01±0.04kcal/mol and root mean squared error of 1.44±0.07kcal/mol, better than the benchmark explicit solvent calculations from FreeSolv, and required less than 15 s of computing time per molecule on a single CPU core. Importantly, parameters for PMVECC showed systematic errors for molecules containing Cl, Br, I, and P. Applying ECC to the explicit solvent hydration free energies found the same systematic errors. The results strongly suggest that some small adjustments to the Lennard-Jones parameters for GAFF will lead to improved hydration free energy calculations for all solvent models.
Collapse
|
4
|
Greene D, Luchko T, Shiferaw Y. The role of subunit cooperativity on ryanodine receptor 2 calcium signaling. Biophys J 2023; 122:215-229. [PMID: 36348625 PMCID: PMC9822801 DOI: 10.1016/j.bpj.2022.11.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2022] [Revised: 10/09/2022] [Accepted: 11/03/2022] [Indexed: 11/09/2022] Open
Abstract
The ryanodine receptor type 2 (RyR2) is composed of four subunits that control calcium (Ca) release in cardiac cells. RyR2 serves primarily as a Ca sensor and can respond to rapid sub-millisecond pulses of Ca while remaining shut at resting concentrations. However, it is not known how the four subunits interact for the RyR2 to function as an effective Ca sensor. To address this question, and to understand the role of subunit cooperativity in Ca-mediated signal transduction, we have developed a computational model of the RyR2 composed of four interacting subunits. We first analyze the statistical properties of a single RyR2 tetramer, where each subunit can exist in a closed or open conformation. Our findings indicate that the number of subunits in the open state is a crucial parameter that dictates RyR2 kinetics. We find that three or four open subunits are required for the RyR2 to harness cooperative interactions to respond to sub-millisecond changes in Ca, while at the same time remaining shut at the resting Ca levels in the cardiac cell. If the required number of open subunits is lowered to one or two, the RyR2 cannot serve as a robust Ca sensor, as the large cooperativity required to stabilize the closed state prevents channel activation. Using this four-subunit model, we analyze the kinetics of Ca release from a RyR2 cluster. We show that the closure of a cluster of RyR2 channels is highly sensitive to the balance of cooperative interactions between closed and open subunits. Based on this result, we analyze how specific interactions between RyR2 subunits can induce persistent Ca leak from the sarcoplasmic reticulum (SR), which is believed to be arrhythmogenic. Thus, these results provide a framework to analyze how a pharmacologic or genetic modification of RyR2 subunit cooperativity can induce abnormal Ca cycling that can potentially lead to life-threatening arrhythmias.
Collapse
Affiliation(s)
- D'Artagnan Greene
- Department of Physics & Astronomy, California State University, Northridge
| | - Tyler Luchko
- Department of Physics & Astronomy, California State University, Northridge
| | - Yohannes Shiferaw
- Department of Physics & Astronomy, California State University, Northridge.
| |
Collapse
|
5
|
Greene D, Barton M, Luchko T, Shiferaw Y. Molecular Dynamics Simulations of the Cardiac Ryanodine Receptor Type 2 (RyR2) Gating Mechanism. J Phys Chem B 2022; 126:9790-9809. [PMID: 36384028 PMCID: PMC9720719 DOI: 10.1021/acs.jpcb.2c03031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to fatal cardiac arrhythmias such as catecholaminergic polymorphic ventricular tachycardia (CPVT). While many CPVT mutations are associated with an increase in Ca2+ leak from the sarcoplasmic reticulum, the mechanistic details of RyR2 channel gating are not well understood, and this poses a barrier in the development of new pharmacological treatments. To address this, we explore the gating mechanism of the RyR2 using molecular dynamics (MD) simulations. We test the effect of changing the conformation of certain structural elements by constructing chimera RyR2 structures that are derived from the currently available closed and open cryo-electron microscopy (cryo-EM) structures, and we then use MD simulations to relax the system. Our key finding is that the position of the S4-S5 linker (S4S5L) on a single subunit can determine whether the channel as a whole is open or closed. Our analysis reveals that the position of the S4S5L is regulated by interactions with the U-motif on the same subunit and with the S6 helix on an adjacent subunit. We find that, in general, channel gating is crucially dependent on high percent occupancy interactions between adjacent subunits. We compare our interaction analysis to 49 CPVT1 mutations in the literature and find that 73% appear near a high percent occupancy interaction between adjacent subunits. This suggests that disruption of cooperative, high percent occupancy interactions between adjacent subunits is a primary cause of channel leak and CPVT in mutant RyR2 channels.
Collapse
|
6
|
Wilson L, Krasny R, Luchko T. Accelerating the 3D reference interaction site model theory of molecular solvation with treecode summation and cut-offs. J Comput Chem 2022; 43:1251-1270. [PMID: 35567580 DOI: 10.1002/jcc.26889] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 02/25/2022] [Accepted: 04/15/2022] [Indexed: 11/07/2022]
Abstract
The 3D reference interaction site model (3D-RISM) of molecular solvation is a powerful tool for computing the equilibrium thermodynamics and density distributions of solvents, such as water and co-ions, around solute molecules. However, 3D-RISM solutions can be expensive to calculate, especially for proteins and other large molecules where calculating the potential energy between solute and solvent requires more than half the computation time. To address this problem, we have developed and implemented treecode summation for long-range interactions and analytically corrected cut-offs for short-range interactions to accelerate the potential energy and long-range asymptotics calculations in non-periodic 3D-RISM in the AmberTools molecular modeling suite. For the largest single protein considered in this work, tubulin, the total computation time was reduced by a factor of 4. In addition, parallel calculations with these new methods scale almost linearly and the iterative solver remains the largest impediment to parallel scaling. To demonstrate the utility of our approach for large systems, we used 3D-RISM to calculate the solvation thermodynamics and density distribution of 7-ring microtubule, consisting of 910 tubulin dimers, over 1.2 million atoms.
Collapse
Affiliation(s)
- Leighton Wilson
- Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Robert Krasny
- Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Los Angeles, California, USA
| |
Collapse
|
7
|
Gray JG, Giambaşu GM, Case DA, Luchko T. Integral equation models for solvent in macromolecular crystals. J Chem Phys 2022; 156:014801. [PMID: 34998331 PMCID: PMC8889494 DOI: 10.1063/5.0070869] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The solvent can occupy up to ∼70% of macromolecular crystals, and hence, having models that predict solvent distributions in periodic systems could improve the interpretation of crystallographic data. Yet, there are few implicit solvent models applicable to periodic solutes, and crystallographic structures are commonly solved assuming a flat solvent model. Here, we present a newly developed periodic version of the 3D-reference interaction site model (RISM) integral equation method that is able to solve efficiently and describe accurately water and ion distributions in periodic systems; the code can compute accurate gradients that can be used in minimizations or molecular dynamics simulations. The new method includes an extension of the Ornstein–Zernike equation needed to yield charge neutrality for charged solutes, which requires an additional contribution to the excess chemical potential that has not been previously identified; this is an important consideration for nucleic acids or any other charged system where most or all the counter- and co-ions are part of the “disordered” solvent. We present several calculations of proteins, RNAs, and small molecule crystals to show that x-ray scattering intensities and the solvent structure predicted by the periodic 3D-RISM solvent model are in closer agreement with the experiment than are intensities computed using the default flat solvent model in the refmac5 or phenix refinement programs, with the greatest improvement in the 2 to 4 Å range. Prospects for incorporating integral equation models into crystallographic refinement are discussed.
Collapse
Affiliation(s)
- Jonathon G Gray
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - George M Giambaşu
- Institute for Quantitative Biomedicine, Rutgers University, Piscataway, New Jersey 08854, USA
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Northridge, California 91330, USA
| |
Collapse
|
8
|
Greene D, Barton M, Luchko T, Shiferaw Y. Computational Analysis of Binding Interactions between the Ryanodine Receptor Type 2 and Calmodulin. J Phys Chem B 2021; 125:10720-10735. [PMID: 34533024 DOI: 10.1021/acs.jpcb.1c03896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to a variety of cardiac arrhythmias, such as catecholaminergic polymorphic ventricular tachycardia (CPVT). RyR2 is regulated by calmodulin (CaM), and mutations that disrupt their interaction can cause aberrant calcium release, leading to an arrhythmia. It was recently shown that increasing the RyR2-CaM binding affinity could rescue a defective CPVT-related RyR2 channel to near wild-type behavior. However, the interactions that determine the binding affinity at the RyR2-CaM binding interface are not well understood. In this study, we identify the key domains and interactions, including several new interactions, involved in the binding of CaM to RyR2. Also, our comparison between the wild-type and V3599K mutant suggests how the RyR2-CaM binding affinity can be increased via a change in the central and N-terminal lobe binding contacts for CaM. This computational approach provides new insights into the effect of a mutation at the RyR2-CaM binding interface, and it may find utility in drug design for the future treatment of cardiac arrhythmias.
Collapse
Affiliation(s)
- D'Artagnan Greene
- Department of Physics, California State University, Northridge, California 91330, United States
| | - Michael Barton
- Department of Physics, California State University, Northridge, California 91330, United States
| | - Tyler Luchko
- Department of Physics, California State University, Northridge, California 91330, United States
| | - Yohannes Shiferaw
- Department of Physics, California State University, Northridge, California 91330, United States
| |
Collapse
|
9
|
Olson B, Cruz A, Chen L, Ghattas M, Ji Y, Huang K, Ayoub S, Luchko T, McKay DJ, Kurtzman T. An online repository of solvation thermodynamic and structural maps of SARS-CoV-2 targets. J Comput Aided Mol Des 2020; 34:1219-1228. [PMID: 32918236 PMCID: PMC7486166 DOI: 10.1007/s10822-020-00341-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Accepted: 08/29/2020] [Indexed: 12/01/2022]
Abstract
SARS-CoV-2 recently jumped species and rapidly spread via human-to-human transmission to cause a global outbreak of COVID-19. The lack of effective vaccine combined with the severity of the disease necessitates attempts to develop small molecule drugs to combat the virus. COVID19_GIST_HSA is a freely available online repository to provide solvation thermodynamic maps of COVID-19-related protein small molecule drug targets. Grid inhomogeneous solvation theory maps were generated using AmberTools cpptraj-GIST, 3D reference interaction site model maps were created with AmberTools rism3d.snglpnt and hydration site analysis maps were created using SSTMap code. The resultant data can be applied to drug design efforts: scoring solvent displacement for docking, rational lead modification, prioritization of ligand- and protein- based pharmacophore elements, and creation of water-based pharmacophores. Herein, we demonstrate the use of the solvation thermodynamic mapping data. It is hoped that this freely provided data will aid in small molecule drug discovery efforts to defeat SARS-CoV-2.
Collapse
Affiliation(s)
- Brian Olson
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA
- Department of Biology and Chemistry, County College of Morris, 214 Center Grove Rd, Randolph, NJ, 07869, USA
| | - Anthony Cruz
- Lehman College Department of Chemistry, 205 W Bedford Park Blvd, Bronx, NY, 10468, USA
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA
| | - Lieyang Chen
- Lehman College Department of Chemistry, 205 W Bedford Park Blvd, Bronx, NY, 10468, USA
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA
| | - Mossa Ghattas
- Lehman College Department of Chemistry, 205 W Bedford Park Blvd, Bronx, NY, 10468, USA
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA
| | - Yeonji Ji
- Lehman College Department of Chemistry, 205 W Bedford Park Blvd, Bronx, NY, 10468, USA
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA
| | - Kunhui Huang
- Lehman College Department of Chemistry, 205 W Bedford Park Blvd, Bronx, NY, 10468, USA
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA
| | - Steven Ayoub
- Department of Chemistry and Biochemistry, California State University, Northridge, 18111 Nordhoff Street, Northridge, CA, 91330, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, 18111 Nordhoff Street, Northridge, CA, 91330, USA
| | - Daniel J McKay
- Ventus Therapeutics, Frederick-Banting, Montreal, QC, H9S 2A1, Canada
| | - Tom Kurtzman
- Lehman College Department of Chemistry, 205 W Bedford Park Blvd, Bronx, NY, 10468, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA.
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, 365 5th Avenue, New York, NY, 10016, USA.
| |
Collapse
|
10
|
McMillin PJ, Alegrete M, Peric M, Luchko T. Electron Paramagnetic Resonance Measurements of Four Nitroxide Probes in Supercooled Water Explained by Molecular Dynamics Simulations. J Phys Chem B 2020; 124:3962-3972. [PMID: 32301326 DOI: 10.1021/acs.jpcb.0c00684] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Electron paramagnetic resonance (EPR) measurements of the rotational diffusion of small nitroxide probes have been demonstrated to be a powerful technique for experimentally investigating the properties of supercooled liquids, such as water. However, since only the rotational diffusion of the probe molecules is measured and EPR measurements are indirect, it is not clear what the relationship between the behavior of water and the probe molecule is. To address this, we have performed molecular dynamics simulations of four nitroxide probes in TIP4P-Ew and OPC water models to directly compare with EPR experiments and to determine the behavior of the water and the underlying microscopic coupling between the water and the probes. In all, 200 ns simulations were run for 23 temperatures between 253 and 283 K for all four probes with each water model for an aggregate of 36.8 μs of simulation time. Simulations for both water models systematically underestimated the rotational diffusion coefficients for both water and probes, though OPC simulations were generally in better agreement with the experiments than TIP4P-Ew simulations. Despite this, when the temperature dependence of the data was fit to a power law, fit parameters for TIP4P-Ew were generally in better agreement with the experiments than OPC. For probe molecules, the singular temperature was found to be T0 = 226.5 ± 0.4 K from experiments, T0 = 208 ± 2 K for OPC water, and T0 = 215 ± 2 K for TIP4P-Ew water. While for water molecules, the singular temperature was found to be T0 = 220.3 ± 0.2 K from experiments, T0 = 208 ± 2 K for OPC water, and T0 = 220 ± 1 K for TIP4P-Ew water. Systematic underestimation of the rotational diffusion coefficients was most pronounced at lower temperatures and was clearly observed in changes to the Arrhenius activation energy. Above the maximum density temperature of Tρmax = 277 K, an activation energy of EA ≈ 16.7 kJ/mol was observed for the probes from experiments, while OPC had EA ≈ 15.2 kJ/mol and TIP4P-Ew had EA ≈ 14.6 kJ/mol. Below the maximum density temperature, the activation energy jumped to EA ≈ 32.5 kJ/mol for experiments but only EA ≈ 23 kJ/mol for OPC and EA ≈ 22 kJ/mol for TIP4P-Ew. In all cases, we saw good agreement between the behavior of the probe molecules and water. To understand why, we calculated the average number of hydrogen bonds between the probe molecules and water. From this, we were able to explain the rotational diffusion times for all of the probes. These results show that current molecular models are sufficient to capture physical phenomena observed with EPR and to help elucidate why the probes provide accurate insights into the behavior of supercooled water.
Collapse
Affiliation(s)
- Patrick J McMillin
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| | - Matthew Alegrete
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| | - Miroslav Peric
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| | - Tyler Luchko
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| |
Collapse
|
11
|
Nguyen C, Yamazaki T, Kovalenko A, Case DA, Gilson MK, Kurtzman T, Luchko T. A molecular reconstruction approach to site-based 3D-RISM and comparison to GIST hydration thermodynamic maps in an enzyme active site. PLoS One 2019; 14:e0219473. [PMID: 31291328 PMCID: PMC6619770 DOI: 10.1371/journal.pone.0219473] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2019] [Accepted: 06/24/2019] [Indexed: 11/25/2022] Open
Abstract
Computed, high-resolution, spatial distributions of solvation energy and entropy can provide detailed information about the role of water in molecular recognition. While grid inhomogeneous solvation theory (GIST) provides rigorous, detailed thermodynamic information from explicit solvent molecular dynamics simulations, recent developments in the 3D reference interaction site model (3D-RISM) theory allow many of the same quantities to be calculated in a fraction of the time. However, 3D-RISM produces atomic-site, rather than molecular, density distributions, which are difficult to extract physical meaning from. To overcome this difficulty, we introduce a method to reconstruct molecular density distributions from atomic-site density distributions. Furthermore, we assess the quality of the resulting solvation thermodynamics density distributions by analyzing the binding site of coagulation Factor Xa with both GIST and 3D-RISM. We find good qualitative agreement between the methods for oxygen and hydrogen densities as well as direct solute-solvent energetic interactions. However, 3D-RISM predicts lower energetic and entropic penalties for moving water from the bulk to the binding site.
Collapse
Affiliation(s)
- Crystal Nguyen
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California, United States of America
| | | | - Andriy Kovalenko
- National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta, Canada
- Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - David A. Case
- Department of Chemistry, Lehman College, The City University of New York, Bronx, New York, United States of America
| | - Michael K. Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California, United States of America
| | - Tom Kurtzman
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, United States of America
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Northridge, California, United States of America
| |
Collapse
|
12
|
Abstract
The Ornstein-Zernike (OZ) integral equation theory is a powerful approach to simple liquids due to its low computational cost and the fact that, when combined with an appropriate closure equation, the theory is thermodynamically complete. However, approximate closures proposed to date exhibit pressure or free energy inconsistencies that produce inaccurate or ambiguous results, limiting the usefulness of the Ornstein-Zernike approach. To address this problem, we combine methods to enforce both pressure and free energy consistency to create a new closure approximation and test it for a single-component Lennard-Jones fluid. The closure is a simple power series in the direct and total correlation functions for which we have derived analytical formulas for the excess Helmholtz free energy and chemical potential. These expressions contain a partial molar volumelike term, similar to excess chemical potential correction terms recently developed. Using our bridge approximation, we have calculated the pressure, Helmholtz free energy, and chemical potential for the Lennard-Jones fluid using the Kirkwood charging, thermodynamic integration techniques, and analytic expressions. These results are compared with those from the hypernetted chain equation and the Verlet-modified closure against Monte Carlo and equations-of-state data for reduced densities of ρ^{*}<1 and temperatures of T^{*}=1.5, 2.74, and 5. Our closure shows consistency among all thermodynamic paths, except for one expression of the Gibbs-Duhem relation, whereas the hypernetted chain equation and the Verlet-modified closure exhibit consistency between only a few relations. Accuracy of the closure is comparable to the Verlet-modified closure and a significant improvement to results obtained from the hypernetted chain equation.
Collapse
Affiliation(s)
- Tsogbayar Tsednee
- Department of Physics and Astronomy, California State University Northridge, 18111 Nordhoff Street, Northridge, California 91330, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University Northridge, 18111 Nordhoff Street, Northridge, California 91330, USA
| |
Collapse
|
13
|
Abstract
Implicit solvent models offer an attractive way to estimate the effects of a solvent environment on the properties of small or large solutes without the complications of explicit simulations. One common test of accuracy is to compute the free energy of transfer from gas to liquid for a variety of small molecules, since many of these values have been measured. Studies of the temperature dependence of these values (i.e. solvation enthalpies and entropies) can provide additional insights into the performance of implicit solvent models. Here, we show how to compute temperature derivatives of hydration free energies for the 3D-RISM integral equation approach. We have computed hydration free energies of 1123 small drug-like molecules (both neutral and charged). Temperature derivatives were also used to calculate hydration energies and entropies of 74 of these molecules (both neutral and charged) for which experimental data is available. While direct results have rather poor agreement with experiment, we have found that several previously proposed linear hydration free energy correction schemes give good agreement with experiment. These corrections also provide good agreement for hydration energies and entropies though simple extensions are required in some cases.
Collapse
Affiliation(s)
- J Johnson
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854
| | - D A Case
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854
| | - T Yamazaki
- Vancouver Prostate Centre, 2660 Oak Street, Vancouver, BC, V6H 3Z6, Canada
| | - S Gusarov
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Dr., Edmonton, AB, T6G 2M9, Canada
| | - A Kovalenko
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Dr., Edmonton, AB, T6G 2M9, Canada
- Department of Mechanical Engineering, University of Alberta, 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Str., Edmonton, AB, T6G 1H9, Canada
| | - T Luchko
- Department of Physics and Astronomy, California State University, Northridge, CA 91330
| |
Collapse
|
14
|
Giambaşu GM, Gebala MK, Panteva MT, Luchko T, Case DA, York DM. Competitive interaction of monovalent cations with DNA from 3D-RISM. Nucleic Acids Res 2015; 43:8405-15. [PMID: 26304542 PMCID: PMC4787805 DOI: 10.1093/nar/gkv830] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2015] [Accepted: 08/06/2015] [Indexed: 12/15/2022] Open
Abstract
The composition of the ion atmosphere surrounding nucleic acids affects their folding, condensation and binding to other molecules. It is thus of fundamental importance to gain predictive insight into the formation of the ion atmosphere and thermodynamic consequences when varying ionic conditions. An early step toward this goal is to benchmark computational models against quantitative experimental measurements. Herein, we test the ability of the three dimensional reference interaction site model (3D-RISM) to reproduce preferential interaction parameters determined from ion counting (IC) experiments for mixed alkali chlorides and dsDNA. Calculations agree well with experiment with slight deviations for salt concentrations >200 mM and capture the observed trend where the extent of cation accumulation around the DNA varies inversely with its ionic size. Ion distributions indicate that the smaller, more competitive cations accumulate to a greater extent near the phosphoryl groups, penetrating deeper into the grooves. In accord with experiment, calculated IC profiles do not vary with sequence, although the predicted ion distributions in the grooves are sequence and ion size dependent. Calculations on other nucleic acid conformations predict that the variation in linear charge density has a minor effect on the extent of cation competition.
Collapse
Affiliation(s)
- George M Giambaşu
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Magdalena K Gebala
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Maria T Panteva
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Tyler Luchko
- Department of Physics & Astronomy, California State University, Northridge, CA 91330, USA
| | - David A Case
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Darrin M York
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| |
Collapse
|
15
|
Giambaşu GM, Luchko T, Herschlag D, York DM, Case DA. Ion counting from explicit-solvent simulations and 3D-RISM. Biophys J 2014; 106:883-94. [PMID: 24559991 PMCID: PMC3944826 DOI: 10.1016/j.bpj.2014.01.021] [Citation(s) in RCA: 83] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Revised: 12/19/2013] [Accepted: 01/10/2014] [Indexed: 12/30/2022] Open
Abstract
The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20-25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na(+) binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.
Collapse
Affiliation(s)
- George M Giambaşu
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey
| | - Tyler Luchko
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California
| | - Darrin M York
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey.
| | - David A Case
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey.
| |
Collapse
|
16
|
Joung IS, Luchko T, Case DA. Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions. J Chem Phys 2013; 138:044103. [PMID: 23387564 DOI: 10.1063/1.4775743] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using the dielectrically consistent reference interaction site model (DRISM) of molecular solvation, we have calculated structural and thermodynamic information of alkali-halide salts in aqueous solution, as a function of salt concentration. The impact of varying the closure relation used with DRISM is investigated using the partial series expansion of order-n (PSE-n) family of closures, which includes the commonly used hypernetted-chain equation (HNC) and Kovalenko-Hirata closures. Results are compared to explicit molecular dynamics (MD) simulations, using the same force fields, and to experiment. The mean activity coefficients of ions predicted by DRISM agree well with experimental values at concentrations below 0.5 m, especially when using the HNC closure. As individual ion activities (and the corresponding solvation free energies) are not known from experiment, only DRISM and MD results are directly compared and found to have reasonably good agreement. The activity of water directly estimated from DRISM is nearly consistent with values derived from the DRISM ion activities and the Gibbs-Duhem equation, but the changes in the computed pressure as a function of salt concentration dominate these comparisons. Good agreement with experiment is obtained if these pressure changes are ignored. Radial distribution functions of NaCl solution at three concentrations were compared between DRISM and MD simulations. DRISM shows comparable water distribution around the cation, but water structures around the anion deviate from the MD results; this may also be related to the high pressure of the system. Despite some problems, DRISM-PSE-n is an effective tool for investigating thermodynamic properties of simple electrolytes.
Collapse
Affiliation(s)
- In Suk Joung
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey 08854, USA
| | | | | |
Collapse
|
17
|
Abstract
The so-called three-dimensional version (3D-RISM) can be used to describe the interactions of solvent components (here we treat water and ions) with a chemical or biomolecular solute of arbitrary size and shape. Here we give an overview of the current status of such models, describing some aspects of “pure” electrolytes (water plus simple ions) and of ionophores, proteins and nucleic acids in the presence of water and salts. Here we focus primarily on interactions with water and dissolved salts; as a practical matter, the discussion is mostly limited to monovalent ions, since studies of divalent ions present many difficult problems that have not yet been addressed. This is not a comprehensive review, but covers a few recent examples that illustrate current issues.
Collapse
Affiliation(s)
- Tyler Luchko
- Department of Chemistry and Chemical Biology and BioMaPS Institute Rutgers University Piscataway NJ 08854, USA
| | - In Suk Joung
- Department of Chemistry and Chemical Biology and BioMaPS Institute Rutgers University Piscataway NJ 08854, USA
| | - David A. Case
- Department of Chemistry and Chemical Biology and BioMaPS Institute Rutgers University Piscataway NJ 08854, USA
| |
Collapse
|
18
|
Freedman H, Luchko T, Luduena RF, Tuszynski JA. Molecular dynamics modeling of tubulin C-terminal tail interactions with the microtubule surface. Proteins 2011; 79:2968-82. [PMID: 21905119 DOI: 10.1002/prot.23155] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2011] [Revised: 05/27/2011] [Accepted: 07/17/2011] [Indexed: 11/10/2022]
Abstract
Tubulin, an α/β heterodimer, has had most of its 3D structure analyzed; however, the carboxy (C)-termini remain elusive. Importantly, the C-termini play critical roles in regulating microtubule structure and function. They are sites of most of the post-translational modifications of tubulin and interaction sites with molecular motors and microtubule-associated proteins. Simulated annealing was used in our molecular dynamics modeling to predict the interactions of the C-terminal tails with the tubulin dimer. We examined differences in their flexibility, interactions with the body of tubulin, and the existence of structural motifs. We found that the α-tubulin tail interacts with the H11 helix of β-tubulin, and the β-tubulin tail interacts with the H11 helix of α-tubulin. Tail domains and H10/B9 loops interact with each other and compete for interactions with positively-charged residues of the H11 helix on the neighboring monomer. In a simulation in which α-tubulin's H10/B9 loop switches on sub-nanosecond intervals between interactions with the C-terminal tail of α-tubulin and the H11 helix of β-tubulin, the intermediate domain of α-tubulin showed more fluctuations compared to those in the other simulations, indicating that tail domains may cause shifts in the position of this domain. This suggests that C-termini may affect the conformation of the tubulin dimer which may explain their essential function in microtubule formation and effects on ligand binding to microtubules. Our modeling also provides evidence for a disordered-helical/helical double-state system of the T3/H3 region of the microtubule, which could be linked to depolymerization following GTP hydrolysis.
Collapse
Affiliation(s)
- Holly Freedman
- CCMAR, FCT, University of Algarve, Campus de Gambelas, Faro, Portugal
| | | | | | | |
Collapse
|
19
|
Abstract
We have modified the popular MM/PBSA or MM/GBSA approaches (molecular mechanics for a biomolecule, combined with a Poisson-Boltzmann or generalized Born electrostatic and surface area nonelectrostatic solvation energy) by employing instead the statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D reference interaction site model, or 3D-RISM-KH) coupled with molecular mechanics or molecular dynamics ( Blinov , N. ; et al. Biophys. J. 2010 ; Luchko , T. ; et al. J. Chem. Theory Comput. 2010 ). Unlike the PBSA or GBSA semiempirical approaches, the 3D-RISM-KH theory yields a full molecular picture of the solvation structure and thermodynamics from the first principles, with proper account of chemical specificities of both solvent and biomolecules, such as hydrogen bonding, hydrophobic interactions, salt bridges, etc. We test the method on the binding of seven biotin analogues to avidin in aqueous solution and show it to work well in predicting the ligand-binding affinities. We have compared the results of 3D-RISM-KH with four different generalized Born and two Poisson-Boltzmann methods. They give absolute binding energies that differ by up to 208 kJ/mol and mean absolute deviations in the relative affinities of 10-43 kJ/mol.
Collapse
Affiliation(s)
- Samuel Genheden
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | | | | | | | | |
Collapse
|
20
|
Luchko T, Gusarov S, Roe DR, Simmerling C, Case DA, Tuszynski J, Kovalenko A. Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber. J Chem Theory Comput 2010; 6:607-624. [PMID: 20440377 PMCID: PMC2861832 DOI: 10.1021/ct900460m] [Citation(s) in RCA: 201] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present the three-dimensional molecular theory of solvation (also known as 3D-RISM) coupled with molecular dynamics (MD) simulation by contracting solvent degrees of freedom, accelerated by extrapolating solvent-induced forces and applying them in large multi-time steps (up to 20 fs) to enable simulation of large biomolecules. The method has been implemented in the Amber molecular modeling package, and is illustrated here on alanine dipeptide and protein G.
Collapse
Affiliation(s)
- Tyler Luchko
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada
- Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada
| | - Sergey Gusarov
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada
| | - Daniel R. Roe
- National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, MD 20899-8443
| | - Carlos Simmerling
- Department of Chemistry, Graduate Program in Biochemistry and Structural Biology, and Center for Structural Biology, Stony Brook University, Stony Brook, New York 11794-3400
- Computational Science Center, Brookhaven National Laboratory, Upton, New York 11973
| | - David A. Case
- BioMaPS Institute, Rutgers University, Piscataway, NJ
- Department of Chemistry & Chemical Biology, Rutgers University, Piscataway, NJ
| | - Jack Tuszynski
- Department of Oncology, University of Alberta, Edmonton, Alberta, Canada
| | - Andriy Kovalenko
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada
- Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada
| |
Collapse
|
21
|
Bennett M, Chik JK, Slysz GW, Luchko T, Tuszynski J, Sackett DL, Schriemer DC. Correction to Structural Mass Spectrometry of the αβ-Tubulin Dimer Supports a Revised Model of Microtubule Assembly. Biochemistry 2009. [DOI: 10.1021/bi901540r] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
22
|
Bennett MJ, Chik JK, Slysz GW, Luchko T, Tuszynski J, Sackett DL, Schriemer DC. Structural mass spectrometry of the alpha beta-tubulin dimer supports a revised model of microtubule assembly. Biochemistry 2009; 48:4858-70. [PMID: 19388626 DOI: 10.1021/bi900200q] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The molecular basis of microtubule lattice instability derives from the hydrolysis of GTP to GDP in the lattice-bound state of alphabeta-tubulin. While this has been appreciated for many years, there is ongoing debate over the molecular basis of this instability and the possible role of altered nucleotide occupancy in the induction of a conformational change in tubulin. The debate has organized around seemingly contradictory models. The allosteric model invokes nucleotide-dependent states of curvature in the free tubulin dimer, such that hydrolysis leads to pronounced bending and thus disruption of the lattice. The more recent lattice model describes a predominant role for the lattice in straightening free dimers that are curved regardless of their nucleotide state. In this model, lattice-bound GTP-tubulin provides the necessary force to straighten an incoming dimer. Interestingly, there is evidence for both models. The enduring nature of this debate stems from a lack of high-resolution data on the free dimer. In this study, we have prepared alphabeta-tubulin samples at high dilution and characterized the nature of nucleotide-induced conformational stability using bottom-up hydrogen/deuterium exchange mass spectrometry (H/DX-MS) coupled with isothermal urea denaturation experiments. These experiments were accompanied by molecular dynamics simulations of the free dimer. We demonstrate an intermediate state unique to GDP-tubulin, suggestive of the curved colchicine-stabilized structure at the intradimer interface but show that intradimer flexibility is an important property of the free dimer regardless of nucleotide occupancy. Our results indicate that the assembly properties of the free dimer may be better described on the basis of this flexibility. A blended model of assembly emerges in which free-dimer allosteric effects retain importance, in an assembly process dominated by lattice-induced effects.
Collapse
Affiliation(s)
- Melissa J Bennett
- Department of Biochemistry and Molecular Biology, University of Calgary, 3330 Hospital Drive NW, Calgary, Alberta, Canada T2N 4N1
| | | | | | | | | | | | | |
Collapse
|
23
|
Freedman H, Huzil JT, Luchko T, Ludueña RF, Tuszynski JA. Identification and characterization of an intermediate taxol binding site within microtubule nanopores and a mechanism for tubulin isotype binding selectivity. J Chem Inf Model 2009; 49:424-36. [PMID: 19434843 DOI: 10.1021/ci8003336] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Tubulin, the primary subunit of microtubules, is remarkable for the variety of small molecules to which it binds. Many of these are very useful or promising agents in cancer chemotherapy. One of the most useful of these is paclitaxel. The tubulin molecule is itself an alpha/beta heterodimer, both alpha- and beta-tubulin monomers existing as multiple isotypes. Despite the success of paclitaxel as an anticancer drug, resistance often occurs in cancer cells and has been associated with variations in tubulin isotype expression, most notably with the increased expression of betaIII-tubulin. Paclitaxel is thought to reach its binding site on beta-tubulin by diffusion through nanopores in the microtubule wall. It has been suggested that a transitional step in this process may be the binding of paclitaxel to an intermediate site within a nanopore, from which it moves directly to its binding site in the microtubule interior facing the lumen. To test this hypothesis, we have computationally docked paclitaxel within a microtubule nanopore and simulated its passage to the intermediate binding site. Targeted molecular dynamics was then used to test the hypothesis that paclitaxel utilizes the H6/H7 loop as a hinge to move directly from this intermediate binding site to its final position in the luminal binding site. We observed that this motion appears to be stabilized by the formation of a hydrogen bond involving serine 275 in beta-tubulin isotypes I, IIa, IIb, IVa, IVb, V, VII, and VIII. Interestingly, this residue is replaced by alanine in the betaIII and VI isotypes. This observation raises the possibility that the observed isotype difference in paclitaxel binding may be a kinetic effect arising from the isotype difference at this residue. We are now able to suggest derivatives of paclitaxel that may reverse the isotype-specificity or lead to an alternate stabilizing hydrogen-bond interaction with tubulin, thus increasing the rate of passage to the luminal binding site and hopefully offering a therapeutic advantage in paclitaxel resistant cases.
Collapse
Affiliation(s)
- Holly Freedman
- Department of Oncology, University of Alberta, Cross Cancer Institute, Edmonton, Alberta, Canada
| | | | | | | | | |
Collapse
|
24
|
Barakat KH, Torin Huzil J, Luchko T, Jordheim L, Dumontet C, Tuszynski J. Characterization of an inhibitory dynamic pharmacophore for the ERCC1-XPA interaction using a combined molecular dynamics and virtual screening approach. J Mol Graph Model 2009; 28:113-30. [PMID: 19473860 DOI: 10.1016/j.jmgm.2009.04.009] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2008] [Revised: 04/17/2009] [Accepted: 04/22/2009] [Indexed: 10/20/2022]
Abstract
Combination chemotherapy involving Cisplatin is a standard treatment for many cancers. However, following an initial positive response, patients will often relapse, presenting with Cisplatin-resistant disease. One possible mechanism for the acquired resistance to Cisplatin is an increase in DNA repair through the up-regulation of ERCC1, an essential component of the nucleotide excision repair complex. Recruitment of ERCC1 to the site of DNA damage is coordinated through its interaction with a protein known as XPA. As there are currently no effective inhibitors of this interaction, inhibition of the ERCC1/XPA interaction may provide an effective strategy for overcoming the development of Cisplatin-resistant cancers. To discover small molecule inhibitors of this interaction, we have screened both the NCI diversity set of ligands and DrugBank-small molecules against the XPA binding site in ERCC1. These compounds were screened using two different techniques in AUTODOCK to account for receptor flexibility. First, using a set of flexible residues, as determined from MD simulations of the XPA/ERCC1 complex and second, using the relaxed complex scheme implemented by performing independent docking experiments against an ensemble of target conformations that were generated from MD simulations. Lowest energy poses from the two different methods were then used to construct a pharmacophore model, which was then validated by comparison to UCN-01, a weak inhibitor of ERCC1 mediated nucleotide excision.
Collapse
Affiliation(s)
- Khaled H Barakat
- Department of Physics, University of Alberta, Edmonton, AB, Canada
| | | | | | | | | | | |
Collapse
|
25
|
Tuszynski JA, Carpenter EJ, Huzil JT, Malinski W, Luchko T, Luduena RF. The evolution of the structure of tubulin and its potential consequences for the role and function of microtubules in cells and embryos. Int J Dev Biol 2009; 50:341-58. [PMID: 16479502 DOI: 10.1387/ijdb.052063jt] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
This paper discusses the results of homology modeling and resulting calculation of key structural and physical properties for close to 300 tubulin sequences, including alpha, beta, gamma, delta and epsilon -tubulins. The basis for our calculations was the structure of the tubulin dimer published several years ago by Nogales et al. (1998), later refined to 3.5 resolution by Lowe et al. (2001). While, it appears that the alpha, beta and gamma-tubulins segregate into distinct structural families, we have found several differences in the physical properties within each group. Each of the alpha, beta and gamma- tubulin groups exhibit major differences in their net electric charge, dipole moments and dipole vector orientations. These properties could influence functional characteristics such as microtubule stability and assembly kinetics, due to their effects on the strength of protein-protein interactions. In addition to the general structural trends between tubulin isoforms, we have observed that the carboxy-termini of alpha and beta-tubulin exists in at least two stable configurations, either projecting away from the tubulin (or microtubule) surface, or collapsed onto the surface. In the latter case, the carboxy-termini form a lattice distinctly different from that of the well-known A and B lattices formed by the tubulin subunits. However, this C-terminal lattice is indistinguishable from the lattice formed when the microtubule-associated protein tau binds to the microtubule surface. Finally, we have discussed how tubulin sequence diversity arose in evolution giving rise to its particular phylogeny and how it may be used in cell- and tissue-specific expression including embryonal development.
Collapse
Affiliation(s)
- Jack A Tuszynski
- Department of Physics, University of Alberta, Edmonton, Alberta, Canada.
| | | | | | | | | | | |
Collapse
|
26
|
Tuszyński JA, Luchko T, Portet S, Dixon JM. Anisotropic elastic properties of microtubules. Eur Phys J E Soft Matter 2005; 17:29-35. [PMID: 15864724 DOI: 10.1140/epje/i2004-10102-5] [Citation(s) in RCA: 70] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2004] [Accepted: 02/10/2005] [Indexed: 05/02/2023]
Abstract
We review and model the experimental parameters which characterize elastic properties of microtubules. Three macroscopic estimates are made of the anisotropic elastic moduli, accounting for the molecular forces between tubulin dimers: for a longitudinal compression of a microtubule, for a lateral force and for a shearing force. These estimates reflect the anisotropies in these parameters observed in several recent experiments.
Collapse
Affiliation(s)
- J A Tuszyński
- Department of Physics, University of Alberta, T6G 2J1, Edmonton, Alberta, Canada,
| | | | | | | |
Collapse
|
27
|
Tuszynski J, Luchko T, Carpenter E, Crawford E. Results of Molecular Dynamics Computations of the Structural and Electrostatic Properties of Tubulin and Their Consequences for Microtubules. ACTA ACUST UNITED AC 2004. [DOI: 10.1166/jctn.2004.042] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|