1
|
Karami Y, Bignon E. Cysteine hyperoxidation rewires communication pathways in the nucleosome and destabilizes the dyad. Comput Struct Biotechnol J 2024; 23:1387-1396. [PMID: 38596314 PMCID: PMC11001638 DOI: 10.1016/j.csbj.2024.03.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 03/27/2024] [Accepted: 03/27/2024] [Indexed: 04/11/2024] Open
Abstract
Gene activity is tightly controlled by reversible chemical modifications called epigenetic marks, which are of various types and modulate gene accessibility without affecting the DNA sequence. Despite an increasing body of evidence demonstrating the role of oxidative-type modifications of histones in gene expression regulation, there remains a complete absence of structural data at the atomistic level to understand the molecular mechanisms behind their regulatory action. Owing to μs time-scale MD simulations and protein communication networks analysis, we describe the impact of histone H3 hyperoxidation (i.e., S-sulfonylation) on the nucleosome core particle dynamics. Our results reveal the atomic-scale details of the intrinsic structural networks within the canonical histone core and their perturbation by hyperoxidation of the histone H3 C110. We show that this modification involves local rearrangements of the communication networks and destabilizes the dyad, and that one modification is enough to induce a maximal structural signature. Our results suggest that cysteine hyperoxidation in the nucleosome core particle might favor its disassembly.
Collapse
Affiliation(s)
- Yasaman Karami
- Université de Lorraine, CNRS, Inria, LORIA, F-54000 Nancy, France
| | | |
Collapse
|
2
|
Marien J, Prévost C, Sacquin-Mora S. nP-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat. J Chem Inf Model 2024; 64:6570-6582. [PMID: 39092904 DOI: 10.1021/acs.jcim.4c00742] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Tau is an intrinsically disordered (IDP) microtubule-associated protein (MAP) that plays a key part in microtubule assembly and organization. The function of tau can be regulated by multiple phosphorylation sites. These post-translational modifications are known to decrease the binding affinity of tau for microtubules, and abnormal tau phosphorylation patterns are involved in Alzheimer's disease. Using all-atom molecular dynamics simulations, we compared the conformational landscapes explored by the tau R2 repeat domain (which comprises a strong tubulin binding site) in its native state and with multiple phosphorylations on the S285, S289, and S293 residues, with four different standard force field (FF)/water model combinations. We find that the different parameters used for the phosphate groups (which can be more or less flexible) in these FFs and the specific interactions between bulk cations and water lead to the formation of a specific type of counterion bridge, termed nP-collab (for nphosphate collaboration, with n being an integer), where counterions form stable structures binding with two or three phosphate groups simultaneously. The resulting effect of nP-collabs on the tau-R2 conformational space differs when using sodium or potassium cations and is likely to impact the peptide overall dynamics and how this MAP interacts with tubulins. We also investigated the effect of phosphoresidue spacing and ionic concentration by modeling polyalanine peptides containing two phosphoserines located one-six residues apart. Three new metrics specifically tailored for IDPs (proteic Menger curvature, local curvature, and local flexibility) were introduced, which allow us to fully characterize the impact of nP-collabs on the dynamics of disordered peptides at the residue level.
Collapse
Affiliation(s)
- Jules Marien
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Chantal Prévost
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Sophie Sacquin-Mora
- Laboratoire de Biochimie Théorique, Université Paris-Cité, CNRS, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
3
|
Hicks CW, Rahman S, Gloor SL, Fields JK, Husby NL, Vaidya A, Maier KE, Morgan M, Wolberger C. Ubiquitinated histone H2B as gatekeeper of the nucleosome acidic patch. Nucleic Acids Res 2024:gkae698. [PMID: 39149911 DOI: 10.1093/nar/gkae698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Revised: 07/15/2024] [Accepted: 08/01/2024] [Indexed: 08/17/2024] Open
Abstract
Monoubiquitination of histones H2B-K120 (H2BK120ub) and H2A-K119 (H2AK119ub) play opposing roles in regulating transcription and chromatin compaction. H2BK120ub is a hallmark of actively transcribed euchromatin, while H2AK119ub is highly enriched in transcriptionally repressed heterochromatin. Whereas H2BK120ub is known to stimulate the binding or activity of various chromatin-modifying enzymes, this post-translational modification (PTM) also interferes with the binding of several proteins to the nucleosome H2A/H2B acidic patch via an unknown mechanism. Here, we report cryoEM structures of an H2BK120ub nucleosome showing that ubiquitin adopts discrete positions that occlude the acidic patch. Molecular dynamics simulations show that ubiquitin remains stably positioned over this nucleosome region. By contrast, our cryoEM structures of H2AK119ub nucleosomes show ubiquitin adopting discrete positions that minimally occlude the acidic patch. Consistent with these observations, H2BK120ub, but not H2AK119ub, abrogates nucleosome interactions with acidic patch-binding proteins RCC1 and LANA, and single-domain antibodies specific to this region. Our results suggest a mechanism by which H2BK120ub serves as a gatekeeper to the acidic patch and point to distinct roles for histone H2AK119 and H2BK120 ubiquitination in regulating protein binding to nucleosomes.
Collapse
Affiliation(s)
- Chad W Hicks
- Department of Biophysics & Biophysical Chemistry, Johns Hopkins University School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
| | - Sanim Rahman
- Department of Biophysics & Biophysical Chemistry, Johns Hopkins University School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
| | - Susan L Gloor
- EpiCypher Inc., 6 Davis Drive, Suite 755, Durham, NC 27709, USA
| | - James K Fields
- Department of Biophysics & Biophysical Chemistry, Johns Hopkins University School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
| | | | - Anup Vaidya
- EpiCypher Inc., 6 Davis Drive, Suite 755, Durham, NC 27709, USA
| | - Keith E Maier
- EpiCypher Inc., 6 Davis Drive, Suite 755, Durham, NC 27709, USA
| | - Michael Morgan
- Department of Biophysics & Biophysical Chemistry, Johns Hopkins University School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
| | - Cynthia Wolberger
- Department of Biophysics & Biophysical Chemistry, Johns Hopkins University School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
| |
Collapse
|
4
|
Srivastava DK, Navratna V, Tosh DK, Chinn A, Sk MF, Tajkhorshid E, Jacobson KA, Gouaux E. Structure of the human dopamine transporter and mechanisms of inhibition. Nature 2024; 632:672-677. [PMID: 39112705 PMCID: PMC11324517 DOI: 10.1038/s41586-024-07739-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Accepted: 06/19/2024] [Indexed: 08/16/2024]
Abstract
The neurotransmitter dopamine has central roles in mood, appetite, arousal and movement1. Despite its importance in brain physiology and function, and as a target for illicit and therapeutic drugs, the human dopamine transporter (hDAT) and mechanisms by which it is inhibited by small molecules and Zn2+ are without a high-resolution structural context. Here we determine the structure of hDAT in a tripartite complex with the competitive inhibitor and cocaine analogue, (-)-2-β-carbomethoxy-3-β-(4-fluorophenyl)tropane2 (β-CFT), the non-competitive inhibitor MRS72923 and Zn2+ (ref. 4). We show how β-CFT occupies the central site, approximately halfway across the membrane, stabilizing the transporter in an outward-open conformation. MRS7292 binds to a structurally uncharacterized allosteric site, adjacent to the extracellular vestibule, sequestered underneath the extracellular loop 4 (EL4) and adjacent to transmembrane helix 1b (TM1b), acting as a wedge, precluding movement of TM1b and closure of the extracellular gate. A Zn2+ ion further stabilizes the outward-facing conformation by coupling EL4 to EL2, TM7 and TM8, thus providing specific insights into how Zn2+ restrains the movement of EL4 relative to EL2 and inhibits transport activity.
Collapse
Affiliation(s)
| | - Vikas Navratna
- Vollum Institute, Oregon Health and Science University, Portland, OR, USA
- Life Sciences Institute, University of Michigan, Ann Arbor, MI, USA
| | - Dilip K Tosh
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, USA
| | - Audrey Chinn
- Vollum Institute, Oregon Health and Science University, Portland, OR, USA
| | - Md Fulbabu Sk
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Biochemistry University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Biochemistry University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Kenneth A Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, USA.
| | - Eric Gouaux
- Vollum Institute, Oregon Health and Science University, Portland, OR, USA.
- Howard Hughes Medical Institute, Oregon Health and Science University, Portland, OR, USA.
| |
Collapse
|
5
|
Zlobin A, Smirnov I, Golovin A. Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases. Protein Sci 2024; 33:e5100. [PMID: 39022909 PMCID: PMC11255601 DOI: 10.1002/pro.5100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Revised: 05/28/2024] [Accepted: 06/19/2024] [Indexed: 07/20/2024]
Abstract
Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.
Collapse
Affiliation(s)
- Alexander Zlobin
- Institute for Drug DiscoveryLeipzig University Medical SchoolLeipzigGermany
- Faculty of Bioengineering and BioinformaticsLomonosov Moscow State UniversityMoscowRussia
| | - Ivan Smirnov
- Shemyakin and Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of SciencesMoscowRussia
| | - Andrey Golovin
- Faculty of Bioengineering and BioinformaticsLomonosov Moscow State UniversityMoscowRussia
- Shemyakin and Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of SciencesMoscowRussia
- Belozersky Institute of Physico‐Chemical BiologyLomonosov Moscow State UniversityMoscowRussia
| |
Collapse
|
6
|
Sankaranarayanan SA, Yadav DN, Yadav S, Srivastava A, Pramatha SR, Kotagiri VR, Joshi H, Rengan AK. Tailoring Phage Nanosomes for Enhanced Theranostic Properties of Near Infrared Dyes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024. [PMID: 39074245 DOI: 10.1021/acs.langmuir.4c01010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
Near-infrared (NIR) phototherapies offer noninvasive, cost-effective solutions for treating tumors and microbial infections. However, organic NIR dyes commonly used suffer from solubility and stability issues requiring frequent dosing. We address this challenge by exploring the bacteriophage-mediated enhancement of NIR dye properties. Upon encapsulation within phage nanosomes, IR780 and Indocyanine green (ICG), with similar optical properties but distinct water solubility and exhibit enhanced UV-vis absorbance and photothermal transduction efficacy compared to liposomes. Experimental characterization corroborated with all-atom molecular dynamics simulations imprints the nanoscale structure, solubility, dynamics, and binding of these NIR dye molecules to the membrane and protein molecules present in Phage capsid. These NIR dye-loaded phage nanosomes, coencapsulated with mitoxantrone, demonstrate enhanced anticancer activity, and when combined with amphotericin B, these dye molecules exhibit superior photothermal effects against fungal infections. Our findings present a simple and efficient approach for tuning the photothermal performance of existing NIR dyes through a rational design for enhanced therapeutic outcomes.
Collapse
Affiliation(s)
| | | | - Saanya Yadav
- Department of Biotechnology, IIT Hyderabad, Kandi, Telangana 502284, India
| | - Aditya Srivastava
- Department of Biomedical Engineering, IIT Hyderabad, Kandi, Telangana 502284, India
| | | | | | - Himanshu Joshi
- Department of Biotechnology, IIT Hyderabad, Kandi, Telangana 502284, India
| | - Aravind Kumar Rengan
- Department of Biomedical Engineering, IIT Hyderabad, Kandi, Telangana 502284, India
| |
Collapse
|
7
|
Biriukov D, Vácha R. Pathways to a Shiny Future: Building the Foundation for Computational Physical Chemistry and Biophysics in 2050. ACS PHYSICAL CHEMISTRY AU 2024; 4:302-313. [PMID: 39069976 PMCID: PMC11274290 DOI: 10.1021/acsphyschemau.4c00003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Revised: 03/15/2024] [Accepted: 03/18/2024] [Indexed: 07/30/2024]
Abstract
In the last quarter-century, the field of molecular dynamics (MD) has undergone a remarkable transformation, propelled by substantial enhancements in software, hardware, and underlying methodologies. In this Perspective, we contemplate the future trajectory of MD simulations and their possible look at the year 2050. We spotlight the pivotal role of artificial intelligence (AI) in shaping the future of MD and the broader field of computational physical chemistry. We outline critical strategies and initiatives that are essential for the seamless integration of such technologies. Our discussion delves into topics like multiscale modeling, adept management of ever-increasing data deluge, the establishment of centralized simulation databases, and the autonomous refinement, cross-validation, and self-expansion of these repositories. The successful implementation of these advancements requires scientific transparency, a cautiously optimistic approach to interpreting AI-driven simulations and their analysis, and a mindset that prioritizes knowledge-motivated research alongside AI-enhanced big data exploration. While history reminds us that the trajectory of technological progress can be unpredictable, this Perspective offers guidance on preparedness and proactive measures, aiming to steer future advancements in the most beneficial and successful direction.
Collapse
Affiliation(s)
- Denys Biriukov
- CEITEC
− Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- National
Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
| | - Robert Vácha
- CEITEC
− Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- National
Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- Department
of Condensed Matter Physics, Faculty of Science, Masaryk University, Kotlářská 267/2, 611 37 Brno, Czech
Republic
| |
Collapse
|
8
|
Wallerstein J, Han X, Levkovets M, Lesovoy D, Malmodin D, Mirabello C, Wallner B, Sun R, Sandalova T, Agback P, Karlsson G, Achour A, Agback T, Orekhov V. Insights into mechanisms of MALT1 allostery from NMR and AlphaFold dynamic analyses. Commun Biol 2024; 7:868. [PMID: 39014105 PMCID: PMC11252132 DOI: 10.1038/s42003-024-06558-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 07/05/2024] [Indexed: 07/18/2024] Open
Abstract
Mucosa-associated lymphoid tissue lymphoma-translocation protein 1 (MALT1) is an attractive target for the development of modulatory compounds in the treatment of lymphoma and other cancers. While the three-dimensional structure of MALT1 has been previously determined through X-ray analysis, its dynamic behaviour in solution has remained unexplored. We present here dynamic analyses of the apo MALT1 form along with the E549A mutation. This investigation used NMR 15N relaxation and NOE measurements between side-chain methyl groups. Our findings confirm that MALT1 exists as a monomer in solution, and demonstrate that the domains display semi-independent movements in relation to each other. Our dynamic study, covering multiple time scales, along with the assessment of conformational populations by Molecular Dynamic simulations, Alpha Fold modelling and PCA analysis, put the side chain of residue W580 in an inward position, shedding light at potential mechanisms underlying the allosteric regulation of this enzyme.
Collapse
Affiliation(s)
- Johan Wallerstein
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden
| | - Xiao Han
- Science for Life Laboratory, Department of Medicine, Solna, Karolinska Institute, SE-17165, Solna, Sweden
- Division of Infectious Diseases, Karolinska University Hospital, SE‑171 76, Stockholm, Sweden
| | - Maria Levkovets
- Swedish NMR Centre, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden
| | - Dmitry Lesovoy
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS, 117997, Moscow, Russia
| | - Daniel Malmodin
- Swedish NMR Centre, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden
| | - Claudio Mirabello
- Dept of Physics, Chemistry and Biology, Linköping University, 581 83, Linköping, Sweden
- National Bioinformatics Infrastructure Sweden, Science for Life Laboratory, Solna, Sweden
| | - Björn Wallner
- National Bioinformatics Infrastructure Sweden, Science for Life Laboratory, Solna, Sweden
| | - Renhua Sun
- Science for Life Laboratory, Department of Medicine, Solna, Karolinska Institute, SE-17165, Solna, Sweden
- Division of Infectious Diseases, Karolinska University Hospital, SE‑171 76, Stockholm, Sweden
| | - Tatyana Sandalova
- Science for Life Laboratory, Department of Medicine, Solna, Karolinska Institute, SE-17165, Solna, Sweden
- Division of Infectious Diseases, Karolinska University Hospital, SE‑171 76, Stockholm, Sweden
| | - Peter Agback
- Department of Molecular Sciences, Swedish University of Agricultural Sciences, PO Box 7015, SE-750 07, Uppsala, Sweden
| | - Göran Karlsson
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden
- Swedish NMR Centre, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden
| | - Adnane Achour
- Science for Life Laboratory, Department of Medicine, Solna, Karolinska Institute, SE-17165, Solna, Sweden
- Division of Infectious Diseases, Karolinska University Hospital, SE‑171 76, Stockholm, Sweden
| | - Tatiana Agback
- Department of Molecular Sciences, Swedish University of Agricultural Sciences, PO Box 7015, SE-750 07, Uppsala, Sweden.
| | - Vladislav Orekhov
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden.
- Swedish NMR Centre, University of Gothenburg, Box 465, SE-40530, Gothenburg, Sweden.
| |
Collapse
|
9
|
Chakraborty R, Crawford-Eng HT, Leburton JP. Asymmetric ion transport through "Janus" MoSSe sub-nanometer pores. NANOSCALE 2024; 16:13106-13120. [PMID: 38912547 DOI: 10.1039/d4nr00589a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
We conduct all-atom molecular dynamics simulations to systematically investigate the underlying mechanisms governing ion transport through a sub-nanometer pore decorated with negative charges in a "Janus" MoSSe membrane. The charge imbalance between S and Se atoms on each side of the membrane induces different types of ion adsorption processes depending on the pore inner charge configuration, and the polarity of external biases, which leads to asymmetry in ionic I-V characteristics. Statistical analysis of the total translocation times including adsorption-desorption processes, and ion dwell times indicates that potassium ions predominantly remain adsorbed during their interaction with the membrane before undertaking a quick translocation through the pore. High applied biases suppress cation adsorption, which results in fast translocation with the current flow boosted by negative inner charges around the pore. We also show that in a membrane consisting of several "Janus" layers, the applied bias necessary to overcome the sub-nm pore barrier increases with the number of layers, providing control over the ionic current.
Collapse
Affiliation(s)
- Rajat Chakraborty
- Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
- Nick Holonyak, Jr. Micro and Nanotechnology Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Henry T Crawford-Eng
- Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
| | - Jean-Pierre Leburton
- Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
- Nick Holonyak, Jr. Micro and Nanotechnology Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
10
|
Mehrafrooz B, Yu L, Pandey L, Siwy ZS, Wanunu M, Aksimentiev A. Electro-osmotic Flow Generation via a Sticky Ion Action. ACS NANO 2024; 18:17521-17533. [PMID: 38832758 PMCID: PMC11233251 DOI: 10.1021/acsnano.4c00829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
Selective transport of ions through nanometer-sized pores is fundamental to cell biology and central to many technological processes such as water desalination and electrical energy storage. Conventional methods for generating ion selectivity include placement of fixed electrical charges at the inner surface of a nanopore through either point mutations in a protein pore or chemical treatment of a solid-state nanopore surface, with each nanopore type requiring a custom approach. Here, we describe a general method for transforming a nanoscale pore into a highly selective, anion-conducting channel capable of generating a giant electro-osmotic effect. Our molecular dynamics simulations and reverse potential measurements show that exposure of a biological nanopore to high concentrations of guanidinium chloride renders the nanopore surface positively charged due to transient binding of guanidinium cations to the protein surface. A comparison of four biological nanopores reveals the relationship between ion selectivity, nanopore shape, composition of the nanopore surface, and electro-osmotic flow. Guanidinium ions are also found to produce anion selectivity and a giant electro-osmotic flow in solid-state nanopores via the same mechanism. Our sticky-ion approach to generate electro-osmotic flow can have numerous applications in controlling molecular transport at the nanoscale and for detection, identification, and sequencing of individual proteins.
Collapse
Affiliation(s)
- Behzad Mehrafrooz
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Luning Yu
- Department of Physics, Northeastern University, Boston, Massachusetts 02115, United States
| | - Laxmi Pandey
- Department of Physics, Northeastern University, Boston, Massachusetts 02115, United States
| | - Zuzanna S Siwy
- Department of Physics, University of California at Irvine, Irvine, California 92697, United States
| | - Meni Wanunu
- Department of Physics, Northeastern University, Boston, Massachusetts 02115, United States
- Department of Bioengineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Aleksei Aksimentiev
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
11
|
Kienlein M, Zacharias M. How arginine inhibits substrate-binding domain 2 elucidated using molecular dynamics simulations. Protein Sci 2024; 33:e5077. [PMID: 38888275 PMCID: PMC11184577 DOI: 10.1002/pro.5077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Revised: 04/19/2024] [Accepted: 05/28/2024] [Indexed: 06/20/2024]
Abstract
The substrate-binding domain 2 (SBD2) is an important part of the bacterial glutamine (GLN) transporter and mediates binding and delivery of GLN to the transporter translocation subunit. The SBD2 consists of two domains, D1 and D2, that bind GLN in the space between domains in a closed structure. In the absence of ligand, the SBD2 adopts an open conformation with larger space between domains. The GLN binding and closing are essential for the subsequent transport into the cell. Arginine (ARG) can also bind to SBD2 but does not induce closing and inhibits GLN transport. We use atomistic molecular dynamics (MD) simulations in explicit solvent to study ARG binding in the presence of the open SBD2 structure and observed reversible binding to the native GLN binding site with similar contacts but no transition to a closed SBD2 state. Absolute binding free energy simulations predict a considerable binding affinity of ARG and GLN to the binding site on the D1 domain. Free energy simulations to induce subsequent closing revealed a strong free energy penalty in case of ARG binding in contrast to GLN binding that favors the closed SBD2 state but still retains a free energy barrier for closing. The simulations allowed the identification of the molecular origin of the closing penalty in case of bound ARG and suggested a mutation of lysine at position 373 to alanine that strongly reduced the penalty and allowed closing even in the presence of bound ARG. The study offers an explanation of the molecular mechanism of how ARG competitively inhibits GLN from binding to SBD2 and from triggering the transition to a closed conformation. The proposed Lys373Ala mutation shows promise as a potential tool to validate whether a conformational mismatch between open SBD2 and the translocator is responsible for preventing ARG uptake to the cell.
Collapse
Affiliation(s)
- Maximilian Kienlein
- Center for Functional Protein Assemblies (CPA)Technical University of MunichGarchingGermany
| | - Martin Zacharias
- Center for Functional Protein Assemblies (CPA)Technical University of MunichGarchingGermany
| |
Collapse
|
12
|
Deng Y, Carnevale V, Ditchfield R, Pletneva EV. Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c. J Phys Chem B 2024; 128:5935-5949. [PMID: 38864552 DOI: 10.1021/acs.jpcb.4c00625] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.
Collapse
Affiliation(s)
- Yunling Deng
- Department of Chemistry, Dartmouth College, Hanover, New Hampshire 03755, United States
| | - Vincenzo Carnevale
- Institute for Genomics and Evolutionary Medicine, Institute for Computational Molecular Science, and Department of Biology, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Robert Ditchfield
- Department of Chemistry, Dartmouth College, Hanover, New Hampshire 03755, United States
| | - Ekaterina V Pletneva
- Department of Chemistry, Dartmouth College, Hanover, New Hampshire 03755, United States
| |
Collapse
|
13
|
Kariyawasam NL, Wereszczynski J. The Influence of Ionic Environment on Nucleosome-Mica Interactions Revealed via Molecular Dynamics Simulations. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.25.600666. [PMID: 38979319 PMCID: PMC11230366 DOI: 10.1101/2024.06.25.600666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Mica serves as a crucial substrate in Atomic Force Microscopy (AFM) studies for visualizing and characterizing nucleosomes. Nucleosomes interact with the negatively charged mica surface via adsorbed cations. However, the specific influences of monovalent and divalent cations on nucleosome adsorption to the mica surface remain unclear. In this study, we investigated the binding of nucleosomes to the mica surface in the presence of monovalent potassium ions and divalent magnesium ions using molecular dynamics simulations. We also explored the impact of pre-treated mica surfaces on nucleosome binding and structure. Our findings reveal that nucleosome-mica interactions vary depending on the cations present, resulting in distinct effects on nucleosome structure. Notably, nucleosomes bind effectively to a mica surface in the presence of potassium ions with minimal structural perturbations.
Collapse
Affiliation(s)
- Nilusha L Kariyawasam
- Department of Physics, Illinois Institute of Technology, Chicago, USA
- Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, USA
| | - Jeff Wereszczynski
- Departments of Physics and Biology, Illinois Institute of Technology, Chicago, USA
- Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, USA
| |
Collapse
|
14
|
Siachouli P, Karadima KS, Mavrantzas VG, Pandis SN. The effect of functional groups on the glass transition temperature of atmospheric organic compounds: a molecular dynamics study. SOFT MATTER 2024; 20:4783-4794. [PMID: 38847330 DOI: 10.1039/d4sm00405a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
Organic compounds constitute a substantial part of atmospheric particulate matter not only in terms of mass concentration but also in terms of distinct functional groups. The glass transition temperature provides an indirect way to investigate the phase state of the organic compounds, playing a crucial role in understanding their behavior and influence on aerosol processes. Molecular dynamics (MD) simulations were implemented here to predict the glass transition temperature (Tg) of atmospherically relevant organic compounds as well as the influence of their functional groups and length of their carbon chain. The cooling step used in the simulations was chosen to be neither too low (to supress crystallization) nor too high (to avoid Tg overprediction). According to the MD simulations, the predicted Tg is sensitive to the functional groups as follows: carboxylic acid (-COOH) > hydroxyl (-OH) and (-COOH) > carbonyls (-CO). Increasing the number of carbon atoms leads to higher Tg for the linearly structured compounds. Linear compounds with lower molecular weight were found to exhibit a lower Tg. No clear correlation between O : C and Tg was observed. The architecture of the carbon chain (linear, or branched, or ring) was also found to impact the glass transition temperature. Compounds containing a non-aromatic carbon ring are characterized by a higher Tg compared to linear and branched ones with the same number of carbon atoms.
Collapse
Affiliation(s)
- Panagiota Siachouli
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
| | - Katerina S Karadima
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
| | - Vlasis G Mavrantzas
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
- Particle Technology Laboratory, Department of Mechanical and Process Engineering, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Spyros N Pandis
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
| |
Collapse
|
15
|
Bošković F, Maffeo C, Patiño-Guillén G, Tivony R, Aksimentiev A, Keyser UF. Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix. ACS NANO 2024; 18:15013-15024. [PMID: 38822455 PMCID: PMC11171748 DOI: 10.1021/acsnano.4c01466] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Revised: 05/13/2024] [Accepted: 05/22/2024] [Indexed: 06/03/2024]
Abstract
Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.
Collapse
Affiliation(s)
- Filip Bošković
- Cavendish
Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| | - Christopher Maffeo
- Department
of Physics, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
- Beckman
Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | | | - Ran Tivony
- Cavendish
Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| | - Aleksei Aksimentiev
- Department
of Physics, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
- Beckman
Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
- Department
of Bioengineering, University of Illinois
at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Ulrich F. Keyser
- Cavendish
Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| |
Collapse
|
16
|
Sun T, Korolev N, Minhas V, Mirzoev A, Lyubartsev AP, Nordenskiöld L. Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles. Biophys J 2024; 123:1414-1434. [PMID: 37915169 PMCID: PMC11163297 DOI: 10.1016/j.bpj.2023.10.030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Revised: 10/05/2023] [Accepted: 10/27/2023] [Indexed: 11/03/2023] Open
Abstract
Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.
Collapse
Affiliation(s)
- Tiedong Sun
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Nikolay Korolev
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Vishal Minhas
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Alexander Mirzoev
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden.
| | - Lars Nordenskiöld
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore.
| |
Collapse
|
17
|
Bošković F, Maffeo C, Patiño-Guillén G, Tivony R, Aksimentiev A, Keyser UF. Nanopore translocation reveals electrophoretic force on non-canonical RNA:DNA double helix. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.09.12.557357. [PMID: 37745457 PMCID: PMC10515835 DOI: 10.1101/2023.09.12.557357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/26/2023]
Abstract
Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of non-canonical RNA:DNA hybrids in electrophoretic transport with the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a non-canonical helix with distinct transport properties from B-form DD molecules. We find RD and DD molecules with the same contour length move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same within experimental errors. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and molecular understanding of electrophoretic transport.
Collapse
|
18
|
Nan Y, MacKerell AD. Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field. J Chem Theory Comput 2024; 20:3242-3257. [PMID: 38588064 PMCID: PMC11039353 DOI: 10.1021/acs.jctc.3c01380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a protocol to enhance the quality of group 1 monatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting preexisting individual parameters for ions, model compounds, and water, we employ atom-pair specific Lennard-Jones (LJ) (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methylacetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.
Collapse
Affiliation(s)
- Yiling Nan
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201 MD
| | - Alexander D. MacKerell
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201 MD
| |
Collapse
|
19
|
Grassmann G, Miotto M, Desantis F, Di Rienzo L, Tartaglia GG, Pastore A, Ruocco G, Monti M, Milanetti E. Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments. Chem Rev 2024; 124:3932-3977. [PMID: 38535831 PMCID: PMC11009965 DOI: 10.1021/acs.chemrev.3c00550] [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: 07/31/2023] [Revised: 02/20/2024] [Accepted: 02/21/2024] [Indexed: 04/11/2024]
Abstract
Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.
Collapse
Affiliation(s)
- Greta Grassmann
- Department
of Biochemical Sciences “Alessandro Rossi Fanelli”, Sapienza University of Rome, Rome 00185, Italy
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Mattia Miotto
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Fausta Desantis
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- The
Open University Affiliated Research Centre at Istituto Italiano di
Tecnologia, Genoa 16163, Italy
| | - Lorenzo Di Rienzo
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Gian Gaetano Tartaglia
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
- Center
for Human Technologies, Genoa 16152, Italy
| | - Annalisa Pastore
- Experiment
Division, European Synchrotron Radiation
Facility, Grenoble 38043, France
| | - Giancarlo Ruocco
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| | - Michele Monti
- RNA
System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
| | - Edoardo Milanetti
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| |
Collapse
|
20
|
Coshic K, Maffeo C, Winogradoff D, Aksimentiev A. The structure and physical properties of a packaged bacteriophage particle. Nature 2024; 627:905-914. [PMID: 38448589 PMCID: PMC11196859 DOI: 10.1038/s41586-024-07150-4] [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: 09/01/2023] [Accepted: 02/01/2024] [Indexed: 03/08/2024]
Abstract
A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.
Collapse
Affiliation(s)
- Kush Coshic
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Christopher Maffeo
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - David Winogradoff
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Aleksei Aksimentiev
- Center for Biophysics and Quantitative Biology, 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.
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
21
|
Smirnova E, Bignon E, Schultz P, Papai G, Ben Shem A. Binding to nucleosome poises human SIRT6 for histone H3 deacetylation. eLife 2024; 12:RP87989. [PMID: 38415718 PMCID: PMC10942634 DOI: 10.7554/elife.87989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024] Open
Abstract
Sirtuin 6 (SIRT6) is an NAD+-dependent histone H3 deacetylase that is prominently found associated with chromatin, attenuates transcriptionally active promoters and regulates DNA repair, metabolic homeostasis and lifespan. Unlike other sirtuins, it has low affinity to free histone tails but demonstrates strong binding to nucleosomes. It is poorly understood how SIRT6 docking on nucleosomes stimulates its histone deacetylation activity. Here, we present the structure of human SIRT6 bound to a nucleosome determined by cryogenic electron microscopy. The zinc finger domain of SIRT6 associates tightly with the acidic patch of the nucleosome through multiple arginine anchors. The Rossmann fold domain binds to the terminus of the looser DNA half of the nucleosome, detaching two turns of the DNA from the histone octamer and placing the NAD+ binding pocket close to the DNA exit site. This domain shows flexibility with respect to the fixed zinc finger and moves with, but also relative to, the unwrapped DNA terminus. We apply molecular dynamics simulations of the histone tails in the nucleosome to show that in this mode of interaction, the active site of SIRT6 is perfectly poised to catalyze deacetylation of the H3 histone tail and that the partial unwrapping of the DNA allows even lysines close to the H3 core to reach the enzyme.
Collapse
Affiliation(s)
- Ekaterina Smirnova
- Department of Integrated Structural Biology, IGBMC, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC)IllkirchFrance
- Université de Strasbourg, IGBMC UMR 7104-UMR-S 1258IllkirchFrance
- CNRS, UMR 7104IllkirchFrance
- Inserm, UMR-S 1258IllkirchFrance
| | | | - Patrick Schultz
- Department of Integrated Structural Biology, IGBMC, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC)IllkirchFrance
- Université de Strasbourg, IGBMC UMR 7104-UMR-S 1258IllkirchFrance
- CNRS, UMR 7104IllkirchFrance
- Inserm, UMR-S 1258IllkirchFrance
| | - Gabor Papai
- Department of Integrated Structural Biology, IGBMC, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC)IllkirchFrance
- Université de Strasbourg, IGBMC UMR 7104-UMR-S 1258IllkirchFrance
- CNRS, UMR 7104IllkirchFrance
- Inserm, UMR-S 1258IllkirchFrance
| | - Adam Ben Shem
- Department of Integrated Structural Biology, IGBMC, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC)IllkirchFrance
- Université de Strasbourg, IGBMC UMR 7104-UMR-S 1258IllkirchFrance
- CNRS, UMR 7104IllkirchFrance
- Inserm, UMR-S 1258IllkirchFrance
| |
Collapse
|
22
|
Song Y, Kim JS. Structure and dynamics of double-stranded DNA rotaxanes. NANOSCALE 2024; 16:4317-4324. [PMID: 38353661 DOI: 10.1039/d3nr05846h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
A DNA rotaxane, with its unique mechanically interlocked architecture consisting of a circular DNA molecule threaded onto a linear DNA axle, holds promise as a fundamental component for nanoscale functional devices. Nevertheless, its structural and dynamic behaviors, essential for advancing molecular machinery, remain largely unexplored. Using extensive all-atom molecular dynamics simulations, we investigated the behaviors of double-stranded DNA (dsDNA) rotaxanes, concentrating on the effects of shape distortion induced by torsional stress in small circular dsDNA containing 70-90 base pairs. We analyzed structural characteristics, including shape, intermolecular distances, and tilt angles, while also exploring dynamic properties such as translational diffusion and toroidal rotation. Our results indicate that shape distortion brings the circular and linear dsDNA components into closer proximity and causes a slight increase in translational diffusion yet a minor decrease in toroidal rotation. Nevertheless, there is no apparent evidence of coupling between translation and rotation. Overall, the insights from this study indicate that such shape distortion does not significantly alter their structure and dynamics. This finding provides flexibility for the design of DNA rotaxanes in nanoscale applications.
Collapse
Affiliation(s)
- Yeonho Song
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea.
| | - Jun Soo Kim
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea.
| |
Collapse
|
23
|
Liu J, Aksimentiev A. Molecular Determinants of Current Blockade Produced by Peptide Transport Through a Nanopore. ACS NANOSCIENCE AU 2024; 4:21-29. [PMID: 38406313 PMCID: PMC10885333 DOI: 10.1021/acsnanoscienceau.3c00046] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 10/28/2023] [Accepted: 11/03/2023] [Indexed: 02/27/2024]
Abstract
The nanopore sensing method holds the promise of delivering a single molecule technology for identification of biological proteins, direct detection of post-translational modifications, and perhaps de novo determination of a protein's amino acid sequence. The key quantity measured in such nanopore sensing experiments is the magnitude of the ionic current passing through a nanopore blocked by a polypeptide chain. Establishing a relationship between the amino acid sequence of a peptide fragment confined within a nanopore and the blockade current flowing through the nanopore remains a major challenge for realizing the nanopore protein sequencing. Using the results of all-atom molecular dynamics simulations, here we compare nanopore sequencing of DNA with nanopore sequencing of proteins. We then delineate the factors affecting the blockade current modulation by the peptide sequence, showing that the current can be determined by (i) the steric footprint of an amino acid, (ii) its interactions with the pore wall, (iii) the local stretching of a polypeptide chain, and (iv) the local enhancement of the ion concentration at the nanopore constriction. We conclude with a brief discussion of the prospects for purely computational prediction of the blockade currents.
Collapse
Affiliation(s)
- Jingqian Liu
- Center
for Biophysics and Quantitative Biology, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
- Beckman
Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Aleksei Aksimentiev
- Center
for Biophysics and Quantitative Biology, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
- Beckman
Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
- Department
of Physics, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
24
|
Chandrasekhar S, Swope TP, Fadaei F, Hollis DR, Bricker R, Houser D, Portman J, Schmidt TL. Bending Unwinds DNA. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.14.579968. [PMID: 38405957 PMCID: PMC10888926 DOI: 10.1101/2024.02.14.579968] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/27/2024]
Abstract
In biology, DNA is often tightly bent to small radii. Solely based on the groove asymmetry, a 30-year-old theoretical paper predicted that such bending should unwind DNA, but this effect has not been directly experimentally quantified so far. We developed a ligation-based assay with nicked DNA circles of variable length, thereby decoupling the twist-dependent ligation efficiency from the large bending strain which dominates conventional circularization assays. We demonstrate that tightly bent DNA indeed unwinds to over 11 base pairs/turn, exactly as predicted. Our discovery requires reassessing the molecular mechanisms and energetics of all processes where DNA is tightly bent or relaxed again, including DNA packaging, gene regulation and expression.
Collapse
Affiliation(s)
- Soumya Chandrasekhar
- Department of Physics, Kent State University, Kent, OH, 44242, USA
- Contributed equally
| | - Thomas P. Swope
- Department of Physics, Kent State University, Kent, OH, 44242, USA
- Contributed equally
| | - Fatemeh Fadaei
- Department of Physics, Kent State University, Kent, OH, 44242, USA
| | - Daniel R. Hollis
- Department of Physics, Kent State University, Kent, OH, 44242, USA
| | - Rachel Bricker
- Department of Physics, Kent State University, Kent, OH, 44242, USA
| | - Draven Houser
- Department of Physics, Kent State University, Kent, OH, 44242, USA
| | - John Portman
- Department of Physics, Kent State University, Kent, OH, 44242, USA
| | | |
Collapse
|
25
|
Riopedre-Fernandez M, Biriukov D, Dračínský M, Martinez-Seara H. Hyaluronan-arginine enhanced and dynamic interaction emerges from distinctive molecular signature due to electrostatics and side-chain specificity. Carbohydr Polym 2024; 325:121568. [PMID: 38008475 DOI: 10.1016/j.carbpol.2023.121568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/01/2023] [Accepted: 11/05/2023] [Indexed: 11/28/2023]
Abstract
Hyaluronan is a natural carbohydrate polymer with a negative charge that fosters gel-like conditions crucial for its cellular functions and industrial applications. As a recognized ligand for proteins, understanding their mutual interactions provides solid ground to tune hyaluronan's gel properties using biocompatible peptides. This work employs NMR and molecular dynamics simulations to identify molecular motifs relevant to hyaluronan-peptide interactions using arginine, lysine, and glycine oligopeptides. Arginine-rich peptides exhibit the strongest binding to hyaluronan according to chemical shift perturbation measurements, followed distantly by the similarly charged lysine. This difference highlights the significance of electrostatics and the peculiarities of the guanidinium side chain in arginine, capable of non-polar interactions that further stabilize the binding. Additional nuclear Overhauser effect measurements do not show stable interaction partners, precluding strong and well-defined complexes. Finally, molecular simulations support our findings and show an extended but significant interaction region, especially for arginine, responsible for the observed enhanced binding, which can also promote cross-linking of hyaluronan polymers. Our findings pave the way for optimizing biocompatible peptides to alter hyaluronan gels' properties efficiently and also explain why hyaluronan-protein interaction typically involves positively charged arginine-rich regions also capable of forming hydrogen bonds and non-polar interactions.
Collapse
Affiliation(s)
- Miguel Riopedre-Fernandez
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 16000, Czech Republic
| | - Denys Biriukov
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 16000, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 16000, Czech Republic
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 16000, Czech Republic.
| |
Collapse
|
26
|
Love O, Winkler L, Cheatham TE. van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD. J Chem Theory Comput 2024; 20:625-643. [PMID: 38157247 PMCID: PMC10809421 DOI: 10.1021/acs.jctc.3c01164] [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: 10/20/2023] [Revised: 12/11/2023] [Accepted: 12/14/2023] [Indexed: 01/03/2024]
Abstract
Molecular dynamics simulations can be used in combination with experimental techniques to uncover the intricacies of biomolecular structure, dynamics, and the resulting interactions. However, many noncanonical nucleic acid structures have proven to be challenging to replicate in accurate agreement with experimental data, often attributed to known force field deficiencies. A common force field criticism is the handling of van der Waals (vdW) parameters, which have not been updated since the regular use of Ewald's methods became routine. This work dives into the effects of minute vdW radii shifts on RNA tetranucleotide, B-DNA, and Z-DNA model systems described by commonly used Amber force fields. Using multidimensional replica exchange molecular dynamics (M-REMD), the GACC RNA tetranucleotide demonstrated changes in the structural distribution between the NMR minor and anomalous structure populations based on the O2' vdW radii scanning. However, no significant change in the NMR Major conformation population was observed. There were minimal changes in the B-DNA structure but there were more substantial improvements in Z-DNA structural descriptions, specifically with the Tumuc1 force field. This occurred with both LJbb vdW radii adjustments and incorporation of the CUFIX nonbonded parameter modifications. Though the limited vdW modifications tested did not provide a universal fix to the challenge of simulating the various known nucleic acid structures, they do provide direction and a greater understanding for future force field development efforts.
Collapse
Affiliation(s)
| | | | - Thomas E. Cheatham
- Department of Medicinal Chemistry,
College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| |
Collapse
|
27
|
Zhou H, Fu H, Shao X, Cai W. Binding Thermodynamics of Fourth-Generation EGFR Inhibitors Revealed by Absolute Binding Free Energy Calculations. J Chem Inf Model 2023; 63:7837-7846. [PMID: 38054791 DOI: 10.1021/acs.jcim.3c01636] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
The overexpression or mutation of the kinase domain of the epidermal growth factor receptor (EGFR) is strongly associated with non-small-cell lung cancer (NSCLC). EGFR tyrosine kinase inhibitors (TKIs) have proven to be effective in treating NSCLC patients. However, EGFR mutations can result in drug resistance. To elucidate the mechanisms underlying this resistance and inform future drug development, we examined the binding affinities of BLU-945, a recently reported fourth-generation TKI, to wild-type EGFR (EGFRWT) and its double-mutant (L858R/T790M; EGFRDM) and triple-mutant (L858R/T790M/C797S; EGFRTM) forms. We compared the binding affinities of BLU-945, BLU-945 analogues, CH7233163 (another fourth-generation TKI), and erlotinib (a first-generation TKI) using absolute binding free energy calculations. Our findings reveal that BLU-945 and CH7233163 exhibit binding affinities to both EGFRDM and EGFRTM stronger than those of erlotinib, corroborating experimental data. We identified K745 and T854 as the key residues in the binding of fourth-generation EGFR TKIs. Electrostatic forces were the predominant driving force for the binding of fourth-generation TKIs to EGFR mutants. Furthermore, we discovered that the incorporation of piperidinol and sulfone groups in BLU-945 substantially enhanced its binding capacity to EGFR mutants. Our study offers valuable theoretical insights for optimizing fourth-generation EGFR TKIs.
Collapse
Affiliation(s)
- Huaxin Zhou
- Research Center for Analytical Sciences, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| | - Haohao Fu
- Research Center for Analytical Sciences, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
- School of Materials Science and Engineering, Smart Sensing Interdisciplinary Science Center, Nankai University, Tianjin 300350, China
| | - Xueguang Shao
- Research Center for Analytical Sciences, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
- School of Materials Science and Engineering, Smart Sensing Interdisciplinary Science Center, Nankai University, Tianjin 300350, China
| | - Wensheng Cai
- Research Center for Analytical Sciences, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
- School of Materials Science and Engineering, Smart Sensing Interdisciplinary Science Center, Nankai University, Tianjin 300350, China
| |
Collapse
|
28
|
Yamada T, Miyazaki Y, Harada S, Kumar A, Vanni S, Shinoda W. Improved Protein Model in SPICA Force Field. J Chem Theory Comput 2023; 19:8967-8977. [PMID: 37989551 DOI: 10.1021/acs.jctc.3c01016] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.
Collapse
Affiliation(s)
- Teppei Yamada
- Graduate School of Natural Science and Technology, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Shogo Harada
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Ashutosh Kumar
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Stefano Vanni
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| |
Collapse
|
29
|
Linse JB, Hub JS. Scrutinizing the protein hydration shell from molecular dynamics simulations against consensus small-angle scattering data. Commun Chem 2023; 6:272. [PMID: 38086909 PMCID: PMC10716392 DOI: 10.1038/s42004-023-01067-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 11/20/2023] [Indexed: 06/09/2024] Open
Abstract
Biological macromolecules in solution are surrounded by a hydration shell, whose structure differs from the structure of bulk solvent. While the importance of the hydration shell for numerous biological functions is widely acknowledged, it remains unknown how the hydration shell is regulated by macromolecular shape and surface composition, mainly because a quantitative probe of the hydration shell structure has been missing. We show that small-angle scattering in solution using X-rays (SAXS) or neutrons (SANS) provide a protein-specific probe of the protein hydration shell that enables quantitative comparison with molecular simulations. Using explicit-solvent SAXS/SANS predictions, we derived the effect of the hydration shell on the radii of gyration Rg of five proteins using 18 combinations of protein force field and water model. By comparing computed Rg values from SAXS relative to SANS in D2O with consensus SAXS/SANS data from a recent worldwide community effort, we found that several but not all force fields yield a hydration shell contrast in remarkable agreement with experiments. The hydration shell contrast captured by Rg values depends strongly on protein charge and geometric shape, thus providing a protein-specific footprint of protein-water interactions and a novel observable for scrutinizing atomistic hydration shell models against experimental data.
Collapse
Affiliation(s)
- Johanna-Barbara Linse
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, 66123, Germany
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, 66123, Germany.
| |
Collapse
|
30
|
Mlýnský V, Kührová P, Stadlbauer P, Krepl M, Otyepka M, Banáš P, Šponer J. Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations. J Chem Theory Comput 2023; 19:8423-8433. [PMID: 37944118 PMCID: PMC10687871 DOI: 10.1021/acs.jctc.3c00990] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2023] [Revised: 10/18/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023]
Abstract
Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Petra Kührová
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Petr Stadlbauer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Miroslav Krepl
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Michal Otyepka
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB−Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB−Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| |
Collapse
|
31
|
Samuel Russell PP, Alaeen S, Pogorelov TV. In-Cell Dynamics: The Next Focus of All-Atom Simulations. J Phys Chem B 2023; 127:9863-9872. [PMID: 37793083 PMCID: PMC10874638 DOI: 10.1021/acs.jpcb.3c05166] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
The cell is a crowded space where large biomolecules and metabolites are in continuous motion. Great strides have been made in in vitro studies of protein dynamics, folding, and protein-protein interactions, and much new data are emerging of how they differ in the cell. In this Perspective, we highlight the current progress in atomistic modeling of in-cell environments, both bacteria and mammals, with emphasis on classical all-atom molecular dynamics simulations. These simulations have been recently used to capture and characterize functional and non-functional protein-protein interactions, protein folding dynamics of small proteins with varied topologies, and dynamics of metabolites. We further discuss the challenges and efforts for updating modern force fields critical to the progress of cellular environment simulations. We also briefly summarize developments in relevant state-of-the-art experimental techniques. As computational and experimental methodologies continue to progress and produce more directly comparable data, we are poised to capture the complex atomistic picture of the cell.
Collapse
Affiliation(s)
- Premila P Samuel Russell
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Sepehr Alaeen
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Taras V Pogorelov
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- School of Chemical Sciences, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- National Center for Supercomputing Applications, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
32
|
Wilson C, Karttunen M, de Groot BL, Gapsys V. Accurately Predicting Protein p Ka Values Using Nonequilibrium Alchemy. J Chem Theory Comput 2023; 19:7833-7845. [PMID: 37820376 PMCID: PMC10653114 DOI: 10.1021/acs.jctc.3c00721] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Indexed: 10/13/2023]
Abstract
The stability, solubility, and function of a protein depend on both its net charge and the protonation states of its individual residues. pKa is a measure of the tendency for a given residue to (de)protonate at a specific pH. Although pKa values can be resolved experimentally, theory and computation provide a compelling alternative. To this end, we assess the applicability of a nonequilibrium (NEQ) alchemical free energy method to the problem of pKa prediction. On a data set of 144 residues that span 13 proteins, we report an average unsigned error of 0.77 ± 0.09, 0.69 ± 0.09, and 0.52 ± 0.04 pK for aspartate, glutamate, and lysine, respectively. This is comparable to current state-of-the-art predictors and the accuracy recently reached using free energy perturbation methods (e.g., FEP+). Moreover, we demonstrate that our open-source, pmx-based approach can accurately resolve the pKa values of coupled residues and observe a substantial performance disparity associated with the lysine partial charges in Amber14SB/Amber99SB*-ILDN, for which an underused fix already exists.
Collapse
Affiliation(s)
- Carter
J. Wilson
- Department
of Mathematics, The University of Western
Ontario, N6A 5B7 London, Canada
- Centre
for Advanced Materials and Biomaterials Research (CAMBR), The University of Western Ontario, N6A 5B7 London, Canada
| | - Mikko Karttunen
- Centre
for Advanced Materials and Biomaterials Research (CAMBR), The University of Western Ontario, N6A 5B7 London, Canada
- Department
of Physics & Astronomy, The University
of Western Ontario, N6A
5B7 London, Canada
- Department
of Chemistry, The University of Western
Ontario, N6A 5B7 London, Canada
| | - Bert L. de Groot
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Multidisciplinary
Sciences, 37077 Göttingen, Germany
| | - Vytautas Gapsys
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Multidisciplinary
Sciences, 37077 Göttingen, Germany
- Computational
Chemistry, Janssen Research & Development, Janssen Pharmaceutica N. V., Turnhoutseweg 30, B-2340 Beerse, Belgium
| |
Collapse
|
33
|
Farafonov VS, Stich M, Nerukh DA. Complete Virion Simulated: All-Atom Model of an MS2 Bacteriophage with Native Genome. J Chem Theory Comput 2023; 19:7924-7933. [PMID: 37856311 DOI: 10.1021/acs.jctc.3c00846] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2023]
Abstract
For the first time, a complete all-atom molecular dynamics (MD) model of a virus, bacteriophage MS2, in its entirety, including a protein outer shell, native genomic RNA with necessary divalent ions, and surrounding explicit aqueous solution with ions at physiological concentration, was built. The model is based on an experimentally measured cryo-EM structure, which was substantially augmented by reconstructing missing or low-resolution parts of the measured density (where the atomistic structure cannot be fit unambiguously). The model was tested by a quarter of a microsecond MD run, and various biophysical characteristics are obtained and analyzed. The developed methodology of building the model can be used for reconstructing other large biomolecular structures when experimental data are fragmented and/or of varying resolution, while the model itself can be used for studying the biology of MS2, including the dynamics of its interaction with the host bacteria.
Collapse
Affiliation(s)
- Vladimir S Farafonov
- Department of Physical Chemistry, V. N. Karazin Kharkiv National University, Kharkiv 61022, Ukraine
- Department of Mathematics, Aston University, Birmingham B4 7ET, U.K
| | - Michael Stich
- Department of Mathematics, Aston University, Birmingham B4 7ET, U.K
- Área de Matemática Aplicada, Departamento de Matemática Aplicada, Ciencia e Ingeniería de Materiales y Tecnología Electrónica, Universidad Rey Juan Carlos, C. Tulipán s/n, Móstoles, Madrid 28933, Spain
| | - Dmitry A Nerukh
- Department of Mathematics, Aston University, Birmingham B4 7ET, U.K
| |
Collapse
|
34
|
Russell PPS, Rickard MM, Boob M, Gruebele M, Pogorelov TV. In silico protein dynamics in the human cytoplasm: Partial folding, misfolding, fold switching, and non-native interactions. Protein Sci 2023; 32:e4790. [PMID: 37774143 PMCID: PMC10578126 DOI: 10.1002/pro.4790] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 08/10/2023] [Accepted: 09/25/2023] [Indexed: 10/01/2023]
Abstract
We examine the influence of cellular interactions in all-atom models of a section of the Homo sapiens cytoplasm on the early folding events of the three-helix bundle protein B (PB). While genetically engineered PB is known to fold in dilute water box simulations in three microseconds, the three initially unfolded PB copies in our two cytoplasm models using a similar force field did not reach the native state during 30-microsecond simulations. We did however capture the formation of all three helices in a compact native-like topology. Folding in vivo is delayed because intramolecular contact formation within PB is in direct competition with intermolecular contacts between PB and surrounding macromolecules. In extreme cases, intermolecular beta-sheets are formed. Interactions with other macromolecules are also observed to promote structure formation, for example when a PB helix in our simulations is shielded from solvent by macromolecular crowding. Sticking and crowding in our models initiate sampling of helix/sheet structural plasticity of PB. Relatedly, in past in vitro experiments, similar GA domains were shown to switch between two different folds. Finally, we also observed that stickiness between PB and the cellular environment can be modulated in our simulations through the reduction in protein hydrophobicity when we reversed PB back to the wild-type sequence. This study demonstrates that even fast-folding proteins can get stuck in non-native states in the cell, making them useful models for protein-chaperone interactions and early stages of aggregate formation relevant to cellular disease.
Collapse
Affiliation(s)
| | - Meredith M. Rickard
- Department of ChemistryUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
| | - Mayank Boob
- Center for Biophysics and Quantitative BiologyUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
| | - Martin Gruebele
- Department of ChemistryUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
- Center for Biophysics and Quantitative BiologyUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
- Beckman Institute for Advanced Science and TechnologyUniversity of Illinois at Urbana‐ChampaignUrbanaIllinoisUSA
- Department of PhysicsUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
| | - Taras V. Pogorelov
- Department of ChemistryUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
- Center for Biophysics and Quantitative BiologyUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
- Beckman Institute for Advanced Science and TechnologyUniversity of Illinois at Urbana‐ChampaignUrbanaIllinoisUSA
- National Center for Supercomputing ApplicationsUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
- School of Chemical SciencesUniversity of Illinois Urbana‐ChampaignUrbanaIllinoisUSA
| |
Collapse
|
35
|
Oleinikov PD, Fedulova AS, Armeev GA, Motorin NA, Singh-Palchevskaia L, Sivkina AL, Feskin PG, Glukhov GS, Afonin DA, Komarova GA, Kirpichnikov MP, Studitsky VM, Feofanov AV, Shaytan AK. Interactions of Nucleosomes with Acidic Patch-Binding Peptides: A Combined Structural Bioinformatics, Molecular Modeling, Fluorescence Polarization, and Single-Molecule FRET Study. Int J Mol Sci 2023; 24:15194. [PMID: 37894874 PMCID: PMC10606924 DOI: 10.3390/ijms242015194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 10/08/2023] [Accepted: 10/11/2023] [Indexed: 10/29/2023] Open
Abstract
In eukaryotic organisms, genomic DNA associates with histone proteins to form nucleosomes. Nucleosomes provide a basis for genome compaction, epigenetic markup, and mediate interactions of nuclear proteins with their target DNA loci. A negatively charged (acidic) patch located on the H2A-H2B histone dimer is a characteristic feature of the nucleosomal surface. The acidic patch is a common site in the attachment of various chromatin proteins, including viral ones. Acidic patch-binding peptides present perspective compounds that can be used to modulate chromatin functioning by disrupting interactions of nucleosomes with natural proteins or alternatively targeting artificial moieties to the nucleosomes, which may be beneficial for the development of new therapeutics. In this work, we used several computational and experimental techniques to improve our understanding of how peptides may bind to the acidic patch and what are the consequences of their binding. Through extensive analysis of the PDB database, histone sequence analysis, and molecular dynamic simulations, we elucidated common binding patterns and key interactions that stabilize peptide-nucleosome complexes. Through MD simulations and FRET measurements, we characterized changes in nucleosome dynamics conferred by peptide binding. Using fluorescence polarization and gel electrophoresis, we evaluated the affinity and specificity of the LANA1-22 peptide to DNA and nucleosomes. Taken together, our study provides new insights into the different patterns of intermolecular interactions that can be employed by natural and designed peptides to bind to nucleosomes, and the effects of peptide binding on nucleosome dynamics and stability.
Collapse
Affiliation(s)
- Pavel D. Oleinikov
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| | | | - Grigoriy A. Armeev
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| | - Nikita A. Motorin
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| | | | - Anastasiia L. Sivkina
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
- Laboratory of Structural-Functional Organization of Chromosomes, Institute of Gene Biology, Russian Academy of Sciences, 119334 Moscow, Russia
| | - Pavel G. Feskin
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| | - Grigory S. Glukhov
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
- Faculty of Biology, MSU-BIT Shenzhen University, Shenzhen 518172, China
| | - Dmitry A. Afonin
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| | - Galina A. Komarova
- Department of Physics, Lomonosov Moscow State University, 119234 Moscow, Russia
| | - Mikhail P. Kirpichnikov
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, 117997 Moscow, Russia
| | - Vasily M. Studitsky
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
- Fox Chase Cancer Center, Philadelphia, PA 19111, USA
| | - Alexey V. Feofanov
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| | - Alexey K. Shaytan
- Department of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
| |
Collapse
|
36
|
Klein F, Soñora M, Helene Santos L, Nazareno Frigini E, Ballesteros-Casallas A, Rodrigo Machado M, Pantano S. The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels. J Struct Biol 2023; 215:107985. [PMID: 37331570 DOI: 10.1016/j.jsb.2023.107985] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 05/18/2023] [Accepted: 06/13/2023] [Indexed: 06/20/2023]
Abstract
The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.
Collapse
Affiliation(s)
- Florencia Klein
- Laboratoire de Biochimie Théorique, UPR9080, CNRS, Paris, France
| | - Martín Soñora
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay
| | | | - Ezequiel Nazareno Frigini
- Instituto Multidisciplinario de Investigaciones Biológicas de San Luis (IMIBIO-SL), Universidad Nacional de San Luis - CONICET, San Luis, Argentina
| | - Andrés Ballesteros-Casallas
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay; Area Bioinformática, DETEMA, Facultad de Química, Universidad de la República, General Flores 2124, Montevideo, 11600, Uruguay
| | | | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay; Area Bioinformática, DETEMA, Facultad de Química, Universidad de la República, General Flores 2124, Montevideo, 11600, Uruguay.
| |
Collapse
|
37
|
Yu L, Kang X, Li F, Mehrafrooz B, Makhamreh A, Fallahi A, Foster JC, Aksimentiev A, Chen M, Wanunu M. Unidirectional single-file transport of full-length proteins through a nanopore. Nat Biotechnol 2023; 41:1130-1139. [PMID: 36624148 PMCID: PMC10329728 DOI: 10.1038/s41587-022-01598-3] [Citation(s) in RCA: 32] [Impact Index Per Article: 32.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 11/02/2022] [Indexed: 01/10/2023]
Abstract
The electrical current blockade of a peptide or protein threading through a nanopore can be used as a fingerprint of the molecule in biosensor applications. However, threading of full-length proteins has only been achieved using enzymatic unfolding and translocation. Here we describe an enzyme-free approach for unidirectional, slow transport of full-length proteins through nanopores. We show that the combination of a chemically resistant biological nanopore, α-hemolysin (narrowest part is ~1.4 nm in diameter), and a high concentration guanidinium chloride buffer enables unidirectional, single-file protein transport propelled by an electroosmotic effect. We show that the mean protein translocation velocity depends linearly on the applied voltage and translocation times depend linearly on length, resembling the translocation dynamics of ssDNA. Using a supervised machine-learning classifier, we demonstrate that single-translocation events contain sufficient information to distinguish their threading orientation and identity with accuracies larger than 90%. Capture rates of protein are increased substantially when either a genetically encoded charged peptide tail or a DNA tag is added to a protein.
Collapse
Affiliation(s)
- Luning Yu
- Department of Physics, Northeastern University, Boston, MA, USA
| | - Xinqi Kang
- Department of Bioengineering, Northeastern University, Boston, MA, USA
| | - Fanjun Li
- Department of Chemistry, University of Massachusetts at Amherst, Amherst, MA, USA
| | - Behzad Mehrafrooz
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Amr Makhamreh
- Department of Bioengineering, Northeastern University, Boston, MA, USA
| | - Ali Fallahi
- Department of Bioengineering, Northeastern University, Boston, MA, USA
| | - Joshua C Foster
- Molecular and Cellular Biology Program, University of Massachusetts at Amherst, Amherst, MA, USA
| | - Aleksei Aksimentiev
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Physics, 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
| | - Min Chen
- Department of Chemistry, University of Massachusetts at Amherst, Amherst, MA, USA
- Molecular and Cellular Biology Program, University of Massachusetts at Amherst, Amherst, MA, USA
| | - Meni Wanunu
- Department of Physics, Northeastern University, Boston, MA, USA.
- Department of Bioengineering, Northeastern University, Boston, MA, USA.
- Chemistry and Chemical Biology, Northeastern University, Boston, MA, USA.
| |
Collapse
|
38
|
Liebl K, Zacharias M. The development of nucleic acids force fields: From an unchallenged past to a competitive future. Biophys J 2023; 122:2841-2851. [PMID: 36540025 PMCID: PMC10398263 DOI: 10.1016/j.bpj.2022.12.022] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 11/08/2022] [Accepted: 12/15/2022] [Indexed: 12/24/2022] Open
Abstract
Molecular dynamics simulations have strongly matured as a method to study biomolecular processes. Their validity, however, is determined by the accuracy of the underlying force fields that describe the forces between all atoms. In this article, we review the development of nucleic acids force fields. We describe the early attempts in the 1990s and emphasize their strong influence on recent force fields. State-of-the-art force fields still use the same Lennard-Jones parameters derived 25 years ago in spite of the fact that these parameters were in general not fitted for nucleic acids. In addition, electrostatic parameters also are deprecated, which may explain some of the current force field deficiencies. We compare different force fields for various systems and discuss new tests of the recently developed Tumuc1 force field. The OL-force fields and Tumuc1 are arguably the best force fields to describe the DNA double helix. However, no force field is flawless. In particular, the description of sugar-puckering remains a problem for nucleic acids force fields. Future refinements are required, so we review methods for force field refinement and give an outlook to the future of force fields.
Collapse
Affiliation(s)
- Korbinian Liebl
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois.
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Munich, Germany
| |
Collapse
|
39
|
Yang D, Zhao Z, Tajkhorshid E, Gouaux E. Structures and membrane interactions of native serotonin transporter in complexes with psychostimulants. Proc Natl Acad Sci U S A 2023; 120:e2304602120. [PMID: 37436958 PMCID: PMC10629533 DOI: 10.1073/pnas.2304602120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 05/03/2023] [Indexed: 07/14/2023] Open
Abstract
The serotonin transporter (SERT) is a member of the SLC6 neurotransmitter transporter family that mediates serotonin reuptake at presynaptic nerve terminals. SERT is the target of both therapeutic antidepressant drugs and psychostimulant substances such as cocaine and methamphetamines, which are small molecules that perturb normal serotonergic transmission by interfering with serotonin transport. Despite decades of studies, important functional aspects of SERT such as the oligomerization state of native SERT and its interactions with potential proteins remain unresolved. Here, we develop methods to isolate SERT from porcine brain (pSERT) using a mild, nonionic detergent, utilize fluorescence-detection size-exclusion chromatography to investigate its oligomerization state and interactions with other proteins, and employ single-particle cryo-electron microscopy to elucidate the structures of pSERT in complexes with methamphetamine or cocaine, providing structural insights into psychostimulant recognition and accompanying pSERT conformations. Methamphetamine and cocaine both bind to the central site, stabilizing the transporter in an outward open conformation. We also identify densities attributable to multiple cholesterol or cholesteryl hemisuccinate (CHS) molecules, as well as to a detergent molecule bound to the pSERT allosteric site. Under our conditions of isolation, we find that pSERT is best described as a monomeric entity, isolated without interacting proteins, and is ensconced by multiple cholesterol or CHS molecules.
Collapse
Affiliation(s)
- Dongxue Yang
- Vollum Institute, Oregon Health and Science University, Portland, OR97239
| | - Zhiyu Zhao
- Department of Biochemistry, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL61801
| | - Emad Tajkhorshid
- Department of Biochemistry, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL61801
| | - Eric Gouaux
- Vollum Institute, Oregon Health and Science University, Portland, OR97239
- HHMI, Oregon Health and Science University, Portland, OR97239
| |
Collapse
|
40
|
Sarthak K, Winogradoff D, Ge Y, Myong S, Aksimentiev A. Benchmarking Molecular Dynamics Force Fields for All-Atom Simulations of Biological Condensates. J Chem Theory Comput 2023; 19:3721-3740. [PMID: 37134270 PMCID: PMC11169342 DOI: 10.1021/acs.jctc.3c00148] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Proteins containing intrinsically disordered regions are integral parts of the cellular signaling pathways and common components of biological condensates. Point mutations in the protein sequence, genetic at birth or acquired through aging, can alter the properties of the condensates, marking the onset of neurodegenerative diseases such as ALS and dementia. While the all-atom molecular dynamics method can, in principle, elucidate the conformational changes that arise from point mutations, the applications of this method to protein condensate systems is conditioned upon the availability of molecular force fields that can accurately describe both structured and disordered regions of such proteins. Using the special-purpose Anton 2 supercomputer, we benchmarked the efficacy of nine presently available molecular force fields in describing the structure and dynamics of a Fused in sarcoma (FUS) protein. Five-microsecond simulations of the full-length FUS protein characterized the effect of the force field on the global conformation of the protein, self-interactions among its side chains, solvent accessible surface area, and the diffusion constant. Using the results of dynamic light scattering as a benchmark for the FUS radius of gyration, we identified several force fields that produced FUS conformations within the experimental range. Next, we used these force fields to perform ten-microsecond simulations of two structured RNA binding domains of FUS bound to their respective RNA targets, finding the choice of the force field to affect stability of the RNA-FUS complex. Taken together, our data suggest that a combination of protein and RNA force fields sharing a common four-point water model provides an optimal description of proteins containing both disordered and structured regions and RNA-protein interactions. To make simulations of such systems available beyond the Anton 2 machines, we describe and validate implementation of the best performing force fields in a publicly available molecular dynamics program NAMD. Our NAMD implementation enables simulations of large (tens of millions of atoms) biological condensate systems and makes such simulations accessible to a broader scientific community.
Collapse
Affiliation(s)
- Kumar Sarthak
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61820, United States
| | - David Winogradoff
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61820, United States
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, Illinois 61820, United States
| | - Yingda Ge
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Sua Myong
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, Illinois 61820, United States
- Program in Cell, Molecular, Developmental Biology, and Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Aleksei Aksimentiev
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, Illinois 61820, United States
| |
Collapse
|
41
|
Baltrukevich H, Bartos P. RNA-protein complexes and force field polarizability. Front Chem 2023; 11:1217506. [PMID: 37426330 PMCID: PMC10323139 DOI: 10.3389/fchem.2023.1217506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Accepted: 06/14/2023] [Indexed: 07/11/2023] Open
Abstract
Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.
Collapse
Affiliation(s)
| | - Piia Bartos
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, Kuopio, Finland
| |
Collapse
|
42
|
Shen K, Nguyen M, Sherck N, Yoo B, Köhler S, Speros J, Delaney KT, Shell MS, Fredrickson GH. Predicting surfactant phase behavior with a molecularly informed field theory. J Colloid Interface Sci 2023; 638:84-98. [PMID: 36736121 DOI: 10.1016/j.jcis.2023.01.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2022] [Revised: 12/24/2022] [Accepted: 01/03/2023] [Indexed: 01/09/2023]
Abstract
HYPOTHESIS The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations. SIMULATION WORK We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements. FINDINGS The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.
Collapse
Affiliation(s)
- Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States; Materials Research Laboratory, University of California, Santa Barbara, Santa Barbara 93106, CA, United States.
| | - My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States
| | - Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States
| | - Brian Yoo
- BASF Corporation, Tarrytown 10591, NY, United States
| | | | - Joshua Speros
- California Research Alliance (CARA) by BASF, Berkeley 94720, CA, United States
| | - Kris T Delaney
- Materials Research Laboratory, University of California, Santa Barbara, Santa Barbara 93106, CA, United States
| | - M Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States.
| | - Glenn H Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States; Materials Research Laboratory, University of California, Santa Barbara, Santa Barbara 93106, CA, United States; Department of Materials Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States.
| |
Collapse
|
43
|
Zhang Z, Šponer J, Bussi G, Mlýnský V, Šulc P, Simmons CR, Stephanopoulos N, Krepl M. Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations. J Chem Inf Model 2023; 63:2794-2809. [PMID: 37126365 PMCID: PMC10170514 DOI: 10.1021/acs.jcim.3c00358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.
Collapse
Affiliation(s)
- Zhengyue Zhang
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- CEITEC─Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
- National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Petr Šulc
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Chad R Simmons
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Nicholas Stephanopoulos
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 241/27, 783 71 Olomouc, Czech Republic
| |
Collapse
|
44
|
Jeong KB, Ryu M, Kim JS, Kim M, Yoo J, Chung M, Oh S, Jo G, Lee SG, Kim HM, Lee MK, Chi SW. Single-molecule fingerprinting of protein-drug interaction using a funneled biological nanopore. Nat Commun 2023; 14:1461. [PMID: 37015934 PMCID: PMC10073129 DOI: 10.1038/s41467-023-37098-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Accepted: 03/01/2023] [Indexed: 04/06/2023] Open
Abstract
In drug discovery, efficient screening of protein-drug interactions (PDIs) is hampered by the limitations of current biophysical approaches. Here, we develop a biological nanopore sensor for single-molecule detection of proteins and PDIs using the pore-forming toxin YaxAB. Using this YaxAB nanopore, we demonstrate label-free, single-molecule detection of interactions between the anticancer Bcl-xL protein and small-molecule drugs as well as the Bak-BH3 peptide. The long funnel-shaped structure and nanofluidic characteristics of the YaxAB nanopore enable the electro-osmotic trapping of diverse folded proteins and high-resolution monitoring of PDIs. Distinctive nanopore event distributions observed in the two-dimensional (ΔI/Io-versus-IN) plot illustrate the ability of the YaxAB nanopore to discriminate individual small-molecule drugs bound to Bcl-xL from non-binders. Taken together, our results present the YaxAB nanopore as a robust platform for label-free, ultrasensitive, single-molecule detection of PDIs, opening up a possibility for low-cost, highly efficient drug discovery against diverse drug targets.
Collapse
Affiliation(s)
- Ki-Baek Jeong
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Critical Diseases Diagnostics Convergence Research Center, KRIBB, Daejeon, 34141, Republic of Korea
| | - Minju Ryu
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Department of Proteome Structural Biology, KRIBB School of Bioscience, University of Science and Technology, Daejeon, 34113, Republic of Korea
| | - Jin-Sik Kim
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Critical Diseases Diagnostics Convergence Research Center, KRIBB, Daejeon, 34141, Republic of Korea
| | - Minsoo Kim
- Department of Physics, Sungkyunkwan University, Suwon, Gyeonggi, 16419, Republic of Korea
| | - Jejoong Yoo
- Department of Physics, Sungkyunkwan University, Suwon, Gyeonggi, 16419, Republic of Korea
| | - Minji Chung
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Sohee Oh
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Gyunghee Jo
- Center for Biomolecular and Cellular Structure, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea
| | - Seong-Gyu Lee
- Center for Biomolecular and Cellular Structure, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea
| | - Ho Min Kim
- Center for Biomolecular and Cellular Structure, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea
- Graduate School of Medical Science and Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Mi-Kyung Lee
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea.
- Critical Diseases Diagnostics Convergence Research Center, KRIBB, Daejeon, 34141, Republic of Korea.
- Department of Proteome Structural Biology, KRIBB School of Bioscience, University of Science and Technology, Daejeon, 34113, Republic of Korea.
| | - Seung-Wook Chi
- Disease Target Structure Research Center, Division of Biomedical Research, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea.
- Department of Proteome Structural Biology, KRIBB School of Bioscience, University of Science and Technology, Daejeon, 34113, Republic of Korea.
- School of Pharmacy, Sungkyunkwan University, Suwon, Gyeonggi, 16419, Republic of Korea.
| |
Collapse
|
45
|
Strelnikov IA, Kovaleva NA, Klinov AP, Zubova EA. C-B-A Test of DNA Force Fields. ACS OMEGA 2023; 8:10253-10265. [PMID: 36969447 PMCID: PMC10034787 DOI: 10.1021/acsomega.2c07781] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 02/24/2023] [Indexed: 06/18/2023]
Abstract
The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.
Collapse
|
46
|
Liebl K, Zacharias M. Toward Force Fields with Improved Base Stacking Descriptions. J Chem Theory Comput 2023; 19:1529-1536. [PMID: 36795949 DOI: 10.1021/acs.jctc.2c01121] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.
Collapse
Affiliation(s)
- Korbinian Liebl
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Garching 85748, Germany
| |
Collapse
|
47
|
Lee SN, Cho HJ, Jeong H, Ryu B, Lee HJ, Kim M, Yoo J, Woo JS, Lee HH. Cryo-EM structures of human Cx36/GJD2 neuronal gap junction channel. Nat Commun 2023; 14:1347. [PMID: 36906653 PMCID: PMC10008584 DOI: 10.1038/s41467-023-37040-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2022] [Accepted: 02/28/2023] [Indexed: 03/13/2023] Open
Abstract
Connexin 36 (Cx36) is responsible for signal transmission in electrical synapses by forming interneuronal gap junctions. Despite the critical role of Cx36 in normal brain function, the molecular architecture of the Cx36 gap junction channel (GJC) is unknown. Here, we determine cryo-electron microscopy structures of Cx36 GJC at 2.2-3.6 Å resolutions, revealing a dynamic equilibrium between its closed and open states. In the closed state, channel pores are obstructed by lipids, while N-terminal helices (NTHs) are excluded from the pore. In the open state with pore-lining NTHs, the pore is more acidic than those in Cx26 and Cx46/50 GJCs, explaining its strong cation selectivity. The conformational change during channel opening also includes the α-to-π-helix transition of the first transmembrane helix, which weakens the protomer-protomer interaction. Our structural analyses provide high resolution information on the conformational flexibility of Cx36 GJC and suggest a potential role of lipids in the channel gating.
Collapse
Affiliation(s)
- Seu-Na Lee
- Department of Life Sciences, Korea University, Seoul, 02841, Republic of Korea
| | - Hwa-Jin Cho
- Department of Chemistry, College of Natural Sciences, Seoul National University, Seoul, 08826, Korea
| | - Hyeongseop Jeong
- Center for Research Equipment, Korea Basic Science Institute, Chungcheongbuk-do, 28119, Korea
| | - Bumhan Ryu
- Research Solution Center, Institute for Basic Science, Daejeon, 34126, Republic of Korea
| | - Hyuk-Joon Lee
- Department of Life Sciences, Korea University, Seoul, 02841, Republic of Korea
| | - Minsoo Kim
- Department of Physics, Sungkyunkwan University, Suwon, 16419, South Korea
| | - Jejoong Yoo
- Department of Physics, Sungkyunkwan University, Suwon, 16419, South Korea
| | - Jae-Sung Woo
- Department of Life Sciences, Korea University, Seoul, 02841, Republic of Korea.
| | - Hyung Ho Lee
- Department of Chemistry, College of Natural Sciences, Seoul National University, Seoul, 08826, Korea.
| |
Collapse
|
48
|
Sarthak K, Winogradoff D, Ge Y, Myong S, Aksimentiev A. Benchmarking Molecular Dynamics Force Fields for All-Atom Simulations of Biological Condensates. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.09.527891. [PMID: 36798393 PMCID: PMC9934651 DOI: 10.1101/2023.02.09.527891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
Abstract
Proteins containing intrinsically disordered regions are integral components of the cellular signaling pathways and common components of biological condensates. Point mutations in the protein sequence, genetic at birth or acquired through aging, can alter the properties of the condensates, marking the onset of neurodegenerative diseases such as ALS and dementia. While all-atom molecular dynamics method can, in principle, elucidate the conformational changes responsible for the aging of the condensate, the applications of this method to protein condensate systems is conditioned by the availability of molecular force fields that can accurately describe both structured and disordered regions of such proteins. Using the special-purpose Anton 2 supercomputer, we benchmarked the efficacy of nine presently available molecular force fields in describing the structure and dynamics of a Fused in sarcoma (FUS) protein. Five-microsecond simulations of the full-length FUS protein characterized the effect of the force field on the global conformation of the protein, self-interactions among its side chains, solvent accessible surface area and the diffusion constant. Using the results of dynamic light scattering as a benchmark for the FUS radius of gyration, we identified several force field that produced FUS conformations within the experimental range. Next, we used these force fields to perform ten-microsecond simulations of two structured RNA binding domains of FUS bound to their respective RNA targets, finding the choice of the force field to affect stability of the RNA-FUS complex. Taken together, our data suggest that a combination of protein and RNA force fields sharing a common four-point water model provides an optimal description of proteins containing both disordered and structured regions and RNA-protein interactions. To make simulations of such systems available beyond the Anton 2 machines, we describe and validate implementation of the best performing force fields in a publicly available molecular dynamics program NAMD. Our NAMD implementation enables simulations of large (tens of millions of atoms) biological condensate systems and makes such simulations accessible to a broader scientific community. Graphical TOC Entry
Collapse
|
49
|
Conformational changes in the human Cx43/GJA1 gap junction channel visualized using cryo-EM. Nat Commun 2023; 14:931. [PMID: 36805660 PMCID: PMC9938869 DOI: 10.1038/s41467-023-36593-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 02/07/2023] [Indexed: 02/20/2023] Open
Abstract
Connexin family proteins assemble into hexameric hemichannels in the cell membrane. The hemichannels dock together between two adjacent membranes to form gap junction intercellular channels (GJIChs). We report the cryo-electron microscopy structures of Cx43 GJICh, revealing the dynamic equilibrium state of various channel conformations in detergents and lipid nanodiscs. We identify three different N-terminal helix conformations of Cx43-gate-covering (GCN), pore-lining (PLN), and flexible intermediate (FIN)-that are randomly distributed in purified GJICh particles. The conformational equilibrium shifts to GCN by cholesteryl hemisuccinates and to PLN by C-terminal truncations and at varying pH. While GJIChs that mainly comprise GCN protomers are occluded by lipids, those containing conformationally heterogeneous protomers show markedly different pore sizes. We observe an α-to-π-helix transition in the first transmembrane helix, which creates a side opening to the membrane in the FIN and PLN conformations. This study provides basic structural information to understand the mechanisms of action and regulation of Cx43 GJICh.
Collapse
|
50
|
Westberry BP, Mansel BW, Lundin L, Williams MAK. Molecular dynamics simulations and X-ray scattering show the κ-carrageenan disorder-to-order transition to be the formation of double-helices. Carbohydr Polym 2023; 302:120417. [PMID: 36604079 DOI: 10.1016/j.carbpol.2022.120417] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2022] [Revised: 10/20/2022] [Accepted: 11/23/2022] [Indexed: 12/12/2022]
Abstract
Recent molecular dynamics simulations, verified experimentally by solution-state x-ray scattering experiments, have found that κ-carrageenan chains contain helical secondary structure, akin to that found in the solid-state, even in aqueous solution. Furthermore, upon the addition of ions to single chains the simulations found no evidence that any conformational transitions take place. These findings challenge the long-held assumption that the so-called disorder-to-order transition in carrageenan systems involves a uni-molecular 'coil-to-helix transition'. Herein, the results of further molecular dynamics simulations undertaken using pairs of κ-carrageenan chains in 0.1 M NaI solutions are reported, and are validated experimentally using state-of-the-art solution-state WAXS experiments. From initially separated chains double-helices are shown to form, leading the authors to propose 'two single helices-to-stabilized double-helix' as a description of the molecular events taking place during the disorder-to-order transition.
Collapse
Affiliation(s)
- Benjamin P Westberry
- School of Natural Sciences, Massey University, Palmerston North, New Zealand; MacDiarmid Institute for Advanced Materials and Nanotechnology, Wellington, New Zealand.
| | - Bradley W Mansel
- National Synchrotron Radiation Research Center, 101 Hsin-Ann Road, Hsinchu Science Park, Hsinchu, 30076, Taiwan, ROC
| | - Leif Lundin
- CSIRO Agriculture and Food, 671 Sneydes Road, Werribee, Victoria 3030, Australia
| | - M A K Williams
- School of Natural Sciences, Massey University, Palmerston North, New Zealand; MacDiarmid Institute for Advanced Materials and Nanotechnology, Wellington, New Zealand
| |
Collapse
|