1
|
Tessmer MH, Stoll S. A novel approach to modeling side chain ensembles of the bifunctional spin label RX. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.05.24.542139. [PMID: 37292623 PMCID: PMC10245940 DOI: 10.1101/2023.05.24.542139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We introduce a novel approach to modeling side chain ensembles of bifunctional spin labels. This approach utilizes rotamer libraries to generate side chain conformational ensembles. Because the bifunctional label is constrained by two attachment sites, the label is split into two monofunctional rotamers which are first attached to their respective sites, then rejoined by a local optimization in dihedral space. We validate this method against a set of previously published experimental data using the bifunctional spin label, RX. This method is relatively fast and can readily be used for both experimental analysis and protein modeling, providing significant advantages over modeling bifunctional labels with molecular dynamics simulations. Use of bifunctional labels for site directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy dramatically reduces label mobility, which can significantly improve resolution of small changes in protein backbone structure and dynamics. Coupling the use of bifunctional labels with side chain modeling methods allows for improved quantitative application of experimental SDSL EPR data to protein modeling.
Collapse
Affiliation(s)
- Maxx H. Tessmer
- Department of Chemistry, University of Washington, Seattle, WA 98103, United States
| | - Stefan Stoll
- Department of Chemistry, University of Washington, Seattle, WA 98103, United States
| |
Collapse
|
2
|
Luo S, Wohl S, Zheng W, Yang S. Biophysical and Integrative Characterization of Protein Intrinsic Disorder as a Prime Target for Drug Discovery. Biomolecules 2023; 13:biom13030530. [PMID: 36979465 PMCID: PMC10046839 DOI: 10.3390/biom13030530] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2023] [Revised: 03/07/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023] Open
Abstract
Protein intrinsic disorder is increasingly recognized for its biological and disease-driven functions. However, it represents significant challenges for biophysical studies due to its high conformational flexibility. In addressing these challenges, we highlight the complementary and distinct capabilities of a range of experimental and computational methods and further describe integrative strategies available for combining these techniques. Integrative biophysics methods provide valuable insights into the sequence–structure–function relationship of disordered proteins, setting the stage for protein intrinsic disorder to become a promising target for drug discovery. Finally, we briefly summarize recent advances in the development of new small molecule inhibitors targeting the disordered N-terminal domains of three vital transcription factors.
Collapse
Affiliation(s)
- Shuqi Luo
- Center for Proteomics and Department of Nutrition, School of Medicine, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Samuel Wohl
- Department of Physics, Arizona State University, Tempe, AZ 85287, USA
| | - Wenwei Zheng
- College of Integrative Sciences and Arts, Arizona State University, Mesa, AZ 85212, USA
- Correspondence: (W.Z.); (S.Y.)
| | - Sichun Yang
- Center for Proteomics and Department of Nutrition, School of Medicine, Case Western Reserve University, Cleveland, OH 44106, USA
- Case Comprehensive Cancer Center, Case Western Reserve University, Cleveland, OH 44106, USA
- Correspondence: (W.Z.); (S.Y.)
| |
Collapse
|
3
|
Mittal S, Dutta S, Shukla D. Reconciling membrane protein simulations with experimental DEER spectroscopy data. Phys Chem Chem Phys 2023; 25:6253-6262. [PMID: 36757376 DOI: 10.1039/d2cp02890e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Spectroscopy experiments are crucial to study membrane proteins for which traditional structure determination methods still prove challenging. Double electron-electron resonance (DEER) spectroscopy experiments provide protein residue-pair distance distributions that are indicative of their conformational heterogeneity. Atomistic molecular dynamics (MD) simulations are another tool that have been proven to be vital to study the structural dynamics of membrane proteins such as to identify inward-open, occluded, and outward-open conformations of transporter membrane proteins, among other partially open or closed states of the protein. Yet, studies have reported that there is no direct consensus between the distributional data from DEER experiments and MD simulations, which has challenged validation of structures obtained from long-timescale simulations and using simulations to design experiments. Current coping strategies for comparisons rely on heuristics, such as mapping the nearest matching peaks between two ensembles or biased simulations. Here we examine the differences in residue-pair distance distributions arising due to the choice of membranes around the protein and covalent modification of a pair of residues to nitroxide spin labels in DEER experiments. Through comparing MD simulations of two proteins, PepTSo and LeuT-both of which have been characterized using DEER experiments previously-we show that the proteins' dynamics are similar despite the choice of the detergent micelle as a membrane mimetic in DEER experiments. On the other hand, covalently modified residues show slight local differences in their dynamics and a huge divergence when the oxygen atom pair distances between spin labeled residues are measured rather than protein backbone distances. Given the computational expense associated with pairwise MTSSL labeled MD simulations, we examine the use of biased simulations to explore the conformational dynamics of the spin labels only to reveal that such simulations alter the underlying protein dynamics. Our study identifies the main cause for the mismatch between DEER experiments and MD simulations and will accelerate the development of potential mitigation strategies to improve the match.
Collapse
Affiliation(s)
- Shriyaa Mittal
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| | - Soumajit Dutta
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Diwakar Shukla
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA. .,Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA.,National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL, USA.,Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.,NIH Center for Macromolecular Modeling and Bioinformatics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
4
|
Mechanism of voltage gating in the voltage-sensing phosphatase Ci-VSP. Proc Natl Acad Sci U S A 2022; 119:e2206649119. [PMID: 36279472 PMCID: PMC9636939 DOI: 10.1073/pnas.2206649119] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Conformational changes in voltage-sensing domains (VSDs) are driven by the transmembrane electric field acting on the protein charges. Yet, the overall energetics and detailed mechanism of this process are not fully understood. Here, we determined free energy and displacement charge landscapes as well as the major conformations visited during a complete functional gating cycle in the isolated VSD of the phosphatase Ci-VSP (Ci-VSD) comprising four transmembrane helices (segments S1 to S4). Molecular dynamics simulations highlight the extent of S4 movements. In addition to the crystallographically determined activated “Up” and resting “Down” states, the simulations predict two Ci-VSD conformations: a deeper resting state (“down-minus”) and an extended activated (“up-plus”) state. These additional conformations were experimentally probed via systematic cysteine mutagenesis with metal-ion bridges and the engineering of proton conducting mutants at hyperpolarizing voltages. The present results show that these four states are visited sequentially in a stepwise manner during voltage activation, each step translocating one arginine or the equivalent of ∼1
e
0
across the membrane electric field, yielding a transfer of ∼3
e
0
charges in total for the complete process.
Collapse
|
5
|
Vermaas JV, Mayne CG, Shinn E, Tajkhorshid E. Assembly and Analysis of Cell-Scale Membrane Envelopes. J Chem Inf Model 2022; 62:602-617. [PMID: 34910495 PMCID: PMC8903035 DOI: 10.1021/acs.jcim.1c01050] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.
Collapse
Affiliation(s)
- Josh V. Vermaas
- Biosciences Center, National Renewable Energy Laboratory, Golden, CO 80401,;
| | - Christopher G. Mayne
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| | - Eric Shinn
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801,;
| |
Collapse
|
6
|
Probing Structural Dynamics of Membrane Proteins Using Electron Paramagnetic Resonance Spectroscopic Techniques. BIOPHYSICA 2021. [DOI: 10.3390/biophysica1020009] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Membrane proteins are essential for the survival of living organisms. They are involved in important biological functions including transportation of ions and molecules across the cell membrane and triggering the signaling pathways. They are targets of more than half of the modern medical drugs. Despite their biological significance, information about the structural dynamics of membrane proteins is lagging when compared to that of globular proteins. The major challenges with these systems are low expression yields and lack of appropriate solubilizing medium required for biophysical techniques. Electron paramagnetic resonance (EPR) spectroscopy coupled with site directed spin labeling (SDSL) is a rapidly growing powerful biophysical technique that can be used to obtain pertinent structural and dynamic information on membrane proteins. In this brief review, we will focus on the overview of the widely used EPR approaches and their emerging applications to answer structural and conformational dynamics related questions on important membrane protein systems.
Collapse
|
7
|
Gamble Jarvi A, Sargun A, Bogetti X, Wang J, Achim C, Saxena S. Development of Cu 2+-Based Distance Methods and Force Field Parameters for the Determination of PNA Conformations and Dynamics by EPR and MD Simulations. J Phys Chem B 2020; 124:7544-7556. [PMID: 32790374 DOI: 10.1021/acs.jpcb.0c05509] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Peptide nucleic acids (PNAs) are a promising group of synthetic analogues of DNA and RNA that offer several distinct advantages over the naturally occurring nucleic acids for applications in biosensing, drug delivery, and nanoelectronics. Because of its structural differences from DNA/RNA, methods to analyze and assess the structure, conformations, and dynamics are needed. In this work, we develop synergistic techniques for the study of the PNA conformation. We use CuQ2, a Cu2+ complex with 8-hydroxyquinoline (HQ), as an alternative base pair and as a spin label in electron paramagnetic resonance (EPR) distance methods. We use molecular dynamics (MD) simulations with newly developed force field parameters for the spin labels to interpret the distance constraints determined by EPR. We complement these methods by UV-vis and circular dichroism measurements and assess the efficacy of the Cu2+ label on a PNA duplex whose backbone is based on aminoethylglycine and a duplex with a hydroxymethyl backbone modification. We show that the Cu2+ label functions efficiently within the standard PNA and the hydroxymethyl-modified PNA and that the MD parameters may be used to accurately reproduce our EPR findings. Through the combination of EPR and MD, we gain new insights into the PNA structure and conformations as well as into the mechanism of orientational selectivity in Cu2+ EPR at X-band. These results present for the first time a rigid Cu2+ spin label used for EPR distance measurements in PNA and the accompanying MD force fields for the spin label. Our studies also reveal that the spin labels have a low impact on the structure of the PNA duplexes. The combined MD and EPR approach represents an important new tool for the characterization of the PNA duplex structure and provides valuable information to aid in the rational application of PNA at large.
Collapse
Affiliation(s)
- Austin Gamble Jarvi
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Artur Sargun
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Xiaowei Bogetti
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15206, United States
| | - Catalina Achim
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Sunil Saxena
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
8
|
Phillips JC, Hardy DJ, Maia JDC, Stone JE, Ribeiro JV, Bernardi RC, Buch R, Fiorin G, Hénin J, Jiang W, McGreevy R, Melo MCR, Radak BK, Skeel RD, Singharoy A, Wang Y, Roux B, Aksimentiev A, Luthey-Schulten Z, Kalé LV, Schulten K, Chipot C, Tajkhorshid E. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J Chem Phys 2020; 153:044130. [PMID: 32752662 PMCID: PMC7395834 DOI: 10.1063/5.0014475] [Citation(s) in RCA: 1332] [Impact Index Per Article: 333.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Accepted: 07/01/2020] [Indexed: 02/06/2023] Open
Abstract
NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.
Collapse
Affiliation(s)
| | - David J. Hardy
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Julio D. C. Maia
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - John E. Stone
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - João V. Ribeiro
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Rafael C. Bernardi
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | | | - Giacomo Fiorin
- National Heart, Lung and Blood Institute, National
Institutes of Health, Bethesda, Maryland 20814,
USA
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique UPR 9080, CNRS
and Université de Paris, Paris, France
| | | | - Ryan McGreevy
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | | | - Brian K. Radak
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Robert D. Skeel
- School of Mathematical and Statistical Sciences,
Arizona State University, Tempe, Arizona 85281,
USA
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State
University, Tempe, Arizona 85281, USA
| | - Yi Wang
- Department of Physics, The Chinese University of
Hong Kong, Shatin, Hong Kong, China
| | - Benoît Roux
- Department of Biochemistry, University of
Chicago, Chicago, Illinois 60637, USA
| | | | | | | | | | - Christophe Chipot
- Authors to whom correspondence should be addressed:
and . URL: http://www.ks.uiuc.edu
| | - Emad Tajkhorshid
- Authors to whom correspondence should be addressed:
and . URL: http://www.ks.uiuc.edu
| |
Collapse
|
9
|
Sahu ID, Lorigan GA. Electron Paramagnetic Resonance as a Tool for Studying Membrane Proteins. Biomolecules 2020; 10:E763. [PMID: 32414134 PMCID: PMC7278021 DOI: 10.3390/biom10050763] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2020] [Revised: 04/24/2020] [Accepted: 04/24/2020] [Indexed: 12/13/2022] Open
Abstract
Membrane proteins possess a variety of functions essential to the survival of organisms. However, due to their inherent hydrophobic nature, it is extremely difficult to probe the structure and dynamic properties of membrane proteins using traditional biophysical techniques, particularly in their native environments. Electron paramagnetic resonance (EPR) spectroscopy in combination with site-directed spin labeling (SDSL) is a very powerful and rapidly growing biophysical technique to study pertinent structural and dynamic properties of membrane proteins with no size restrictions. In this review, we will briefly discuss the most commonly used EPR techniques and their recent applications for answering structure and conformational dynamics related questions of important membrane protein systems.
Collapse
Affiliation(s)
- Indra D. Sahu
- Natural Science Division, Campbellsville University, Campbellsville, KY 42718, USA
- Department of Chemistry and Biochemistry, Miami University, Oxford, OH 45056, USA
| | - Gary A. Lorigan
- Department of Chemistry and Biochemistry, Miami University, Oxford, OH 45056, USA
| |
Collapse
|
10
|
Gangopadhyay A, Chakraborty HJ, Datta A. Employing virtual screening and molecular dynamics simulations for identifying hits against the active cholera toxin. Toxicon 2019; 170:1-9. [PMID: 31494206 DOI: 10.1016/j.toxicon.2019.09.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Revised: 08/22/2019] [Accepted: 09/01/2019] [Indexed: 12/24/2022]
Abstract
Cholera is a major global threat, affecting millions each year. The ADP ribosyltransferase activity of the active cholera toxin catalyses the massive loss of water and electrolytes during cholera infections. The active toxin heterodimer comprises the A1 subunit from Vibrio cholerae and ARF (ADP Ribosylation Factor) from the human host. Although the active toxin is a potential target for drug discovery against cholera, it has been scarcely targeted to date. The A1-ARF interface contains a potential druggable site for small molecule inhibitors. By combining a sequential docking and scoring strategy with molecular dynamics (MD) simulations, this study identified hits against the protein-protein interface (PPI) of the active cholera toxin from an in-house library of 9,175 ADMET-screened alkaloids. The docking algorithms and scoring functions of Glide SP, Glide XP, and AutoDock were employed for initial library screening. Three alkaloids were initially selected by docking-based virtual screening. The stability of the hit-toxin complexes was validated by MD simulations. Two of the three hits, namely, A6225 (7-formyldehydrothalicsimidine) and A16503 (1,2,7,8-tetrahydroxy dibenz[cd,f]indol-4(5H)-one), formed stable complexes with the toxin. Analyses of the hydrogen bond occupancies revealed that the hits formed stable hydrogen bonds with the toxin PPI. The hits identified herein can serve as reference compounds for drug discovery against cholera in the future.
Collapse
Affiliation(s)
- Aditi Gangopadhyay
- Department of Chemical Technology, University of Calcutta, 92, APC Road, Kolkata 700009, West Bengal, India; DBT Centre for Bioinformatics, Presidency University, Kolkata 700073, West Bengal, India.
| | - Hirak Jyoti Chakraborty
- DBT Centre for Bioinformatics, Presidency University, Kolkata 700073, West Bengal, India; Central Inland Fisheries Research Institute, Barrackpore, Kolkata 700120, West Bengal, India
| | - Abhijit Datta
- DBT Centre for Bioinformatics, Presidency University, Kolkata 700073, West Bengal, India; Department of Botany, Jhargram Raj College, Jhargram 721507, Paschim Medinipur, India
| |
Collapse
|
11
|
Qi Y, Lee J, Cheng X, Shen R, Islam SM, Roux B, Im W. CHARMM-GUI DEER facilitator for spin-pair distance distribution calculations and preparation of restrained-ensemble molecular dynamics simulations. J Comput Chem 2019; 41:415-420. [PMID: 31329318 DOI: 10.1002/jcc.26032] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Revised: 06/11/2019] [Accepted: 06/29/2019] [Indexed: 12/21/2022]
Abstract
The double electron-electron resonance (DEER) is a powerful structural biology technique to obtain distance information in the range of 18 to 80 å by measuring the dipolar coupling between two unpaired electron spins. The distance distributions obtained from the experiment provide valuable structural information about the protein in its native environment that can be exploited using restrained ensemble molecular dynamics (reMD) simulations. We present a new tool DEER Facilitator in CHARMM-GUI that consists of two modules Spin-Pair Distributor and reMD Prepper to setup simulations that utilize information from DEER experiments. Spin-Pair Distributor provides a web-based interface to calculate the spin-pair distance distribution of labeled sites in a protein using MD simulations. The calculated distribution can be used to guide the selection of the labeling sites in experiments as well as validate different protein structure models. reMD Prepper facilities the setup of reMD simulations using different types of spin labels in four different environments including vacuum, solution, micelle, and bilayer. The applications of these two modules are demonstrated with several test cases. Spin-Pair Distributor and reMD Prepper are available at http://www.charmm-gui.org/input/deer and http://www.charmm-gui.org/input/deerre. DEER Facilitator is expected to facilitate advanced biomolecular modeling and simulation, thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems based on experimental DEER data. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Yifei Qi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Jumin Lee
- Departments of Biological Sciences and Bioengineering, Lehigh University, Bethlehem, Pennsylvania, 18015
| | - Xi Cheng
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Rong Shen
- Department of Biochemistry and Molecular Biology and Department of Chemistry, University of Chicago, Chicago, Illinois, 60637
| | - Shahidul M Islam
- Department of Chemistry, University of Illinois at Chicago, Chicago, Illinois, 60607
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology and Department of Chemistry, University of Chicago, Chicago, Illinois, 60637
| | - Wonpil Im
- Departments of Biological Sciences and Bioengineering, Lehigh University, Bethlehem, Pennsylvania, 18015
| |
Collapse
|
12
|
Poursoleiman A, Karimi-Jafari MH, Zolmajd-Haghighi Z, Bagheri M, Haertlé T, Behbehani GR, Ghasemi A, Stroylova YY, Muronetz VI, Saboury AA. Polymyxins interaction to the human serum albumin: A thermodynamic and computational study. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2019; 217:155-163. [PMID: 30933779 DOI: 10.1016/j.saa.2019.03.077] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2018] [Revised: 03/18/2019] [Accepted: 03/22/2019] [Indexed: 06/09/2023]
Abstract
Polymyxin B and E (colistin), are a group of cationic charged cyclic antibiotic lipopeptides that are frequently used in the clinics to treat infections caused by the multidrug-resistant gram-negative bacteria. Since the interactions with the blood plasma drug-transport proteins may play a critical role in determining their pharmacological and pharmacokinetic profiles, we studied the binding properties of polymyxins to the human serum albumin (HSA) under simulated physiological conditions by the combination of biophysical approaches, such as isothermal titration calorimetry (ITC), fluorescence anisotropy, circular dichroism (CD) buttressed by computational studies. The HSA binding to the polymyxins was relatively strong (Ka ≈ 1.0 × 107 M-1). Molecular docking indicated that polymyxins bind to the cleft of HSA between domains I and III via the electrostatic interactions. This evidence was further confirmed by the entropy-driven interaction for the polymyxins bound HSA. Far UV-CD experiments showed that the secondary structure of HSA doesn't alter and its stable structure is preserved. Collectively, these investigations revealed that the polymyxins bind preferentially to the partially unfolded intermediate forms of the protein structure; however, HSA molecule does not undergo any significant conformational changes upon binding. This is promising as it may limit the unfavorable side effects of the medicine. On the whole, the results provide quantitative and qualitative insight of the binding interaction between HSA and polymyxins, which is important in understanding their effect as therapeutic agents.
Collapse
Affiliation(s)
- A Poursoleiman
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - M H Karimi-Jafari
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Z Zolmajd-Haghighi
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - M Bagheri
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - T Haertlé
- Poznan University of Life Sciences, Department of Animal Nutrition, Poznan, Poland; Biopolymers, Interactions, Assemblies, UR 1268, Institute National de la Recherche Agronomique, Nantes, France
| | - G Rezaei Behbehani
- Chemistry Department, Imam Khomeini International University, Qazvin, Iran
| | - A Ghasemi
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Y Y Stroylova
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia
| | - V I Muronetz
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia
| | - A A Saboury
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
| |
Collapse
|
13
|
Garate JA, Bernardin A, Escalona Y, Yanez C, English NJ, Perez-Acle T. Orientational and Folding Thermodynamics via Electric Dipole Moment Restraining. J Phys Chem B 2019; 123:2599-2608. [PMID: 30831028 DOI: 10.1021/acs.jpcb.8b09374] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The projection of molecular processes onto a small set of relevant descriptors, the so-called reaction coordinates or collective variables (CVs), is a technique nowadays routinely employed by the biomolecular simulation community. In this work, we implemented two CVs to manipulate the orientation (i.e., angle) (μ⃗a) and magnitude (|μ⃗|) of the electric dipole moment. In doing so, we studied the thermodynamics of water orientation under the application of external voltages and the folding of two polypeptides at zero-field conditions. The projection of the free-energy [potential of mean force (PMF)] along water orientation defined an upper limit of around 0.3 V for irrelevant thermodynamic effects. On the other hand, sufficiently strong μ⃗a restraints applied on 12-alanine (Ala12) triggered structural effects because of the alignment of local dipoles; for lower restraints, a full-body rotation is achieved. The manipulation of |μ⃗| produced strong perturbations on the secondary structure of Ala12, promoting an enhanced sampling to its configurational space. Rigorous free-energy calculations in the form of 2-D PMFs for deca-alanine showed the utility of |μ⃗| as a reaction coordinate to study folding in small α helices. As a whole, we propose that the manipulation of both components of the dipole moment, μ⃗a and |μ⃗|, provides thermodynamics insights into the structural conformation and stability of biomolecules. These new CVs are implemented in the Colvars module, available for NAMD and LAMMPS.
Collapse
Affiliation(s)
- Jose Antonio Garate
- Centro Interdisciplinario de Neurociencia de Valparaiso , Universidad de Valparaiso , Pasaje Harrington 287 , Playa Ancha, Valparaiso 2381850 , Chile
| | - Alejandro Bernardin
- Centro Interdisciplinario de Neurociencia de Valparaiso , Universidad de Valparaiso , Pasaje Harrington 287 , Playa Ancha, Valparaiso 2381850 , Chile.,Computational Biology Lab , Fundacion Ciencia & Vida , Avenida Zanartu 1482, Nunoa , Santiago 7780272 , Chile
| | - Yerko Escalona
- Institute for Molecular Modeling and Simulation , Muthgasse 18 , Vienna 1190 , Austria
| | - Carlos Yanez
- Computational Biology Lab , Fundacion Ciencia & Vida , Avenida Zanartu 1482, Nunoa , Santiago 7780272 , Chile
| | - Niall J English
- School of Chemical and Bioprocess Engineering , University College Dublin , Belfield, Dublin 4 , Ireland
| | - Tomas Perez-Acle
- Centro Interdisciplinario de Neurociencia de Valparaiso , Universidad de Valparaiso , Pasaje Harrington 287 , Playa Ancha, Valparaiso 2381850 , Chile.,Computational Biology Lab , Fundacion Ciencia & Vida , Avenida Zanartu 1482, Nunoa , Santiago 7780272 , Chile
| |
Collapse
|
14
|
Peter EK, Černý J. A Hybrid Hamiltonian for the Accelerated Sampling along Experimental Restraints. Int J Mol Sci 2019; 20:E370. [PMID: 30654563 PMCID: PMC6359555 DOI: 10.3390/ijms20020370] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2018] [Revised: 01/07/2019] [Accepted: 01/10/2019] [Indexed: 12/25/2022] Open
Abstract
In this article, we present an enhanced sampling method based on a hybrid Hamiltonian which combines experimental distance restraints with a bias dependent from multiple path-dependent variables. This simulation method determines the bias-coordinates on the fly and does not require a priori knowledge about reaction coordinates. The hybrid Hamiltonian accelerates the sampling of proteins, and, combined with experimental distance information, the technique considers the restraints adaptively and in dependency of the system's intrinsic dynamics. We validate the methodology on the dipole relaxation of two water models and the conformational landscape of dialanine. Using experimental NMR-restraint data, we explore the folding landscape of the TrpCage mini-protein and in a second example apply distance restraints from chemical crosslinking/mass spectrometry experiments for the sampling of the conformation space of the Killer Cell Lectin-like Receptor Subfamily B Member 1A (NKR-P1A). The new methodology has the potential to adaptively introduce experimental restraints without affecting the conformational space of the system along an ergodic trajectory. Since only a limited number of input- and no-order parameters are required for the setup of the simulation, the method is broadly applicable and has the potential to be combined with coarse-graining methods.
Collapse
Affiliation(s)
- Emanuel K Peter
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, 252 50 Vestec, Prague West, Czech Republic.
| | - Jiří Černý
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, 252 50 Vestec, Prague West, Czech Republic.
| |
Collapse
|
15
|
Saranya V, Shankar R, Vijayakumar S. Structural exploration of viral matrix protein 40 interaction with the transition metal ions (Ag+ and Cu2+). J Biomol Struct Dyn 2018; 37:2875-2896. [DOI: 10.1080/07391102.2018.1498803] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- V. Saranya
- Department of Physics, Bharathiar University, Coimbatore, India
| | - R. Shankar
- Department of Physics, Bharathiar University, Coimbatore, India
| | - S. Vijayakumar
- Department of Medical Physics, Bharathiar University, Coimbatore, India
| |
Collapse
|
16
|
Samanta P, Wang Y, Fuladi S, Zou J, Li Y, Shen L, Weber C, Khalili-Araghi F. Molecular determination of claudin-15 organization and channel selectivity. J Gen Physiol 2018; 150:949-968. [PMID: 29915162 PMCID: PMC6028499 DOI: 10.1085/jgp.201711868] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 03/08/2018] [Accepted: 05/04/2018] [Indexed: 12/28/2022] Open
Abstract
Members of the claudin family form tight junctions between adjacent epithelial and endothelial cells. Samanta et al. build an atomic model of claudin-15 using molecular dynamics simulations and conclude that four claudin-15 molecules each contribute an aspartic acid residue to form a selectivity filter. Tight junctions are macromolecular structures that traverse the space between adjacent cells in epithelia and endothelia. Members of the claudin family are known to determine tight junction permeability in a charge- and size-selective manner. Here, we use molecular dynamics simulations to build and refine an atomic model of claudin-15 channels and study its transport properties. Our simulations indicate that claudin-15 forms well-defined channels for ions and molecules and otherwise “seals” the paracellular space through hydrophobic interactions. Ionic currents, calculated from simulation trajectories of wild-type as well as mutant channels, reflect in vitro measurements. The simulations suggest that the selectivity filter is formed by a cage of four aspartic acid residues (D55), contributed by four claudin-15 molecules, which creates a negative electrostatic potential to favor cation flux over anion flux. Charge reversal or charge ablation mutations of D55 significantly reduce cation permeability in silico and in vitro, whereas mutations of other negatively charged pore amino acid residues have a significantly smaller impact on channel permeability and selectivity. The simulations also indicate that water and small ions can pass through the channel, but larger cations, such as tetramethylammonium, do not traverse the pore. Thus, our model provides an atomic view of claudin channels, their transport function, and a potential three-dimensional organization of its selectivity filter.
Collapse
Affiliation(s)
| | - Yitang Wang
- Department of Pathology, The University of Chicago, Chicago, IL.,Department of Surgery, The University of Chicago, Chicago, IL
| | - Shadi Fuladi
- Department of Physics, University of Illinois, Chicago, IL
| | - Jinjing Zou
- Department of Pathology, The University of Chicago, Chicago, IL
| | - Ye Li
- Department of Pathology, The University of Chicago, Chicago, IL
| | - Le Shen
- Department of Pathology, The University of Chicago, Chicago, IL .,Department of Surgery, The University of Chicago, Chicago, IL
| | | | | |
Collapse
|
17
|
Cassidy CK, Himes BA, Luthey-Schulten Z, Zhang P. CryoEM-based hybrid modeling approaches for structure determination. Curr Opin Microbiol 2018; 43:14-23. [PMID: 29107896 PMCID: PMC5934336 DOI: 10.1016/j.mib.2017.10.002] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2017] [Revised: 10/04/2017] [Accepted: 10/09/2017] [Indexed: 12/21/2022]
Abstract
Recent advances in cryo-electron microscopy (cryoEM) have dramatically improved the resolutions at which vitrified biological specimens can be studied, revealing new structural and mechanistic insights over a broad range of spatial scales. Bolstered by these advances, much effort has been directed toward the development of hybrid modeling methodologies for the construction and refinement of high-fidelity atomistic models from cryoEM data. In this brief review, we will survey the key elements of cryoEM-based hybrid modeling, providing an overview of available computational tools and strategies as well as several recent applications.
Collapse
Affiliation(s)
- C Keith Cassidy
- Department of Physics, Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Benjamin A Himes
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Zaida Luthey-Schulten
- Department of Chemistry, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Peijun Zhang
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA; Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Roosevelt Drive, Oxford OX3 7BN, UK; Electron Bio-Imaging Centre, Diamond Light Sources, Harwell Science and Innovation Campus, Didcot OX11 0DE, UK.
| |
Collapse
|
18
|
Rydzewski J, Jakubowski R, Nicosia G, Nowak W. Conformational Sampling of a Biomolecular Rugged Energy Landscape. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 15:732-739. [PMID: 27913358 DOI: 10.1109/tcbb.2016.2634008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The protein structure refinement using conformational sampling is important in hitherto protein studies. In this paper, we examined the protein structure refinement by means of potential energy minimization using immune computing as a method of sampling conformations. The method was tested on the x-ray structure and 30 decoys of the mutant of [Leu]Enkephalin, a paradigmatic example of the biomolecular multiple-minima problem. In order to score the refined conformations, we used a standard potential energy function with the OPLSAA force field. The effectiveness of the search was assessed using a variety of methods. The robustness of sampling was checked by the energy yield function which measures quantitatively the number of the peptide decoys residing in an energetic funnel. Furthermore, the potential energy-dependent Pareto fronts were calculated to elucidate dissimilarities between peptide conformations and the native state as observed by x-ray crystallography. Our results showed that the probed potential energy landscape of [Leu]Enkephalin is self-similar on different metric scales and that the local potential energy minima of the peptide decoys are metastable, thus they can be refined to conformations whose potential energy is decreased by approximately 250 kJ/mol.
Collapse
|
19
|
Gianti E, Carnevale V. Computational Approaches to Studying Voltage-Gated Ion Channel Modulation by General Anesthetics. Methods Enzymol 2018; 602:25-59. [DOI: 10.1016/bs.mie.2018.01.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
20
|
Jo S, Cheng X, Lee J, Kim S, Park SJ, Patel DS, Beaven AH, Lee KI, Rui H, Park S, Lee HS, Roux B, MacKerell AD, Klauda JB, Qi Y, Im W. CHARMM-GUI 10 years for biomolecular modeling and simulation. J Comput Chem 2017; 38:1114-1124. [PMID: 27862047 PMCID: PMC5403596 DOI: 10.1002/jcc.24660] [Citation(s) in RCA: 179] [Impact Index Per Article: 25.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2016] [Revised: 10/04/2016] [Accepted: 10/18/2016] [Indexed: 12/16/2022]
Abstract
CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Sunhwan Jo
- Leadership Computing Facility, Argonne National Laboratory, 9700 Cass Ave, Argonne, Illinois
| | - Xi Cheng
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica Chinese Academy of Sciences, 555 Zuchongzhi Road, Pudong, Shanghai, 201203, China
| | - Jumin Lee
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Seonghoon Kim
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Sang-Jun Park
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Dhilon S Patel
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Andrew H Beaven
- Department of Chemistry, The University of Kansas, Lawrence, Kansas
| | - Kyu Il Lee
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Huan Rui
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois
| | - Soohyung Park
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Hui Sun Lee
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences School of Pharmacy, University of Maryland, Baltimore, Maryland
| | - Jeffrey B Klauda
- Department of Chemical and Biomolecular Engineering and the Biophysics Program, University of Maryland College Park, Maryland
| | - Yifei Qi
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| | - Wonpil Im
- Department of Biological Sciences and Bioengineering Program, Lehigh University, Pennsylvania
| |
Collapse
|
21
|
Metal binding mediated conformational change of XPA protein:a potential cytotoxic mechanism of nickel in the nucleotide excision repair. J Mol Model 2016; 22:156. [PMID: 27307058 DOI: 10.1007/s00894-016-3017-x] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2016] [Accepted: 05/26/2016] [Indexed: 12/13/2022]
Abstract
Nucleotide excision repair (NER) is a pivotal life process for repairing DNA nucleotide mismatch caused by chemicals, metal ions, radiation, and other factors. As the initiation step of NER, the xeroderma pigmentosum complementation group A protein (XPA) recognizes damaged DNA molecules, and recruits the replication protein A (RPA), another important player in the NER process. The stability of the Zn(2+)-chelated Zn-finger domain of XPA center core portion (i.e., XPA98-210) is the foundation of its biological functionality, while the displacement of the Zn(2+) by toxic metal ions (such as Ni(2+), a known human carcinogen and allergen) may impair the effectiveness of NER and hence elevate the chance of carcinogenesis. In this study, we first calculated the force field parameters for the bonded model in the metal center of the XPA98-210 system, showing that the calculated results, including charges, bonds, angles etc., are congruent with previously reported results measured by spectrometry experiments and quantum chemistry computation. Then, comparative molecular dynamics simulations using these parameters revealed the changes in the conformation and motion mode of XPA98-210 Zn-finger after the substitution of Zn(2+) by Ni(2+). The results showed that Ni(2+) dramatically disrupted the relative positions of the four Cys residues in the Zn-finger structure, forcing them to collapse from a tetrahedron into an almost planar structure. Finally, we acquired the binding mode of XPA98-210 with its ligands RPA70N and DNA based on molecular docking and structural alignment. We found that XPA98-210's Zn-finger domain primarily binds to a V-shaped cleft in RPA70N, while the cationic band in its C-terminal subdomain participates in the recognition of damaged DNA. In addition, this article sheds light on the multi-component interaction pattern among XPA, DNA, and other NER-related proteins (i.e., RPA70N, RPA70A, RPA70B, RPA70C, RPA32, and RPA14) based on previously reported structural biology information. Thus, we derived a putative cytotoxic mechanism associated with the nickel ion, where the Ni(2+) disrupts the conformation of the XPA Zn-finger, directly weakening its interaction with RPA70N, and thus lowering the effectiveness of the NER process. In sum, this work not only provides a theoretical insight into the multi-protein interactions involved in the NER process and potential cytotoxic mechanism associated with Ni(2+) binding in XPA, but may also facilitate rational anti-cancer drug design based on the NER mechanism.
Collapse
|