1
|
Koca Fındık B, Jafari M, Song LF, Li Z, Aviyente V, Merz KM. Binding of Phosphate Species to Ca 2+ and Mg 2+ in Aqueous Solution. J Chem Theory Comput 2024; 20:4298-4307. [PMID: 38718258 PMCID: PMC11137831 DOI: 10.1021/acs.jctc.4c00218] [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/21/2024] [Revised: 04/27/2024] [Accepted: 05/01/2024] [Indexed: 05/29/2024]
Abstract
Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H2PO4-, HPO42-, PO43- coordinating with Ca2+ and Mg2+. First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca2+ and Mg2+ generally show similar trends at longer metal-phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation-anion interactions.
Collapse
Affiliation(s)
- Basak Koca Fındık
- Department
of Chemistry, Bogazici University, Bebek, Istanbul 34342, Turkey
| | - Majid Jafari
- Department
of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
| | - Lin Frank Song
- Biochemical
and Biophysical Systems Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Zhen Li
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Viktorya Aviyente
- Department
of Chemistry, Bogazici University, Bebek, Istanbul 34342, Turkey
| | - Kenneth M. Merz
- Department
of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
2
|
Forma A, Grunwald A, Zembala P, Januszewski J, Brachet A, Zembala R, Świątek K, Baj J. Micronutrient Status and Breast Cancer: A Narrative Review. Int J Mol Sci 2024; 25:4968. [PMID: 38732186 PMCID: PMC11084730 DOI: 10.3390/ijms25094968] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2024] [Revised: 04/28/2024] [Accepted: 04/30/2024] [Indexed: 05/13/2024] Open
Abstract
Breast cancer is one of the most common cancers worldwide, at the same time being one of the most prevalent causes of women's death. Many factors such as alcohol, weight fluctuations, or hormonal replacement therapy can potentially contribute to breast cancer development and progression. Another important factor in breast cancer onset includes micronutrient status. In this narrative review, we analyzed 23 micronutrients and their possible influence on breast cancer onset and progression. Further, the aim of this study was to investigate the impact of micronutrient status on the prevention of breast cancer and its possible influence on various therapeutic pathways. We researched meta-analyses, systemic and narrative reviews, retrospective studies, as well as original studies on human and animal models. The results of these studies indicate a possible correlation between the different levels of micronutrients and a decreased risk of breast cancer as well as a better survival rate. However, further studies are necessary to establish adequate doses of supplementation of the chosen micronutrients and the exact mechanisms of micronutrient impact on breast cancer therapy.
Collapse
Affiliation(s)
- Alicja Forma
- Chair and Department of Forensic Medicine, Medical University of Lublin, Jaczewskiego 8b, 20-090 Lublin, Poland; (A.G.); (A.B.)
| | - Arkadiusz Grunwald
- Chair and Department of Forensic Medicine, Medical University of Lublin, Jaczewskiego 8b, 20-090 Lublin, Poland; (A.G.); (A.B.)
| | - Patryk Zembala
- Faculty of Medicine, Medical University of Warsaw, Banacha 1A, 02-097 Warsaw, Poland;
| | - Jacek Januszewski
- Department of Correct, Clinical and Imaging Anatomy, Chair of Fundamental Sciences, Medical University of Lublin, Jaczewskiego 4, 20-090 Lublin, Poland; (J.J.); (K.Ś.); (J.B.)
| | - Adam Brachet
- Chair and Department of Forensic Medicine, Medical University of Lublin, Jaczewskiego 8b, 20-090 Lublin, Poland; (A.G.); (A.B.)
| | - Roksana Zembala
- Faculty of Medicine, Cardinal Stefan Wyszynski University in Warsaw, Wóycickiego 1/3, 01-938 Warsaw, Poland;
| | - Kamila Świątek
- Department of Correct, Clinical and Imaging Anatomy, Chair of Fundamental Sciences, Medical University of Lublin, Jaczewskiego 4, 20-090 Lublin, Poland; (J.J.); (K.Ś.); (J.B.)
| | - Jacek Baj
- Department of Correct, Clinical and Imaging Anatomy, Chair of Fundamental Sciences, Medical University of Lublin, Jaczewskiego 4, 20-090 Lublin, Poland; (J.J.); (K.Ś.); (J.B.)
| |
Collapse
|
3
|
Mondal B, Chakraborty D, Hori N, Nguyen HT, Thirumalai D. Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences. J Chem Theory Comput 2024; 20:2934-2946. [PMID: 38498914 DOI: 10.1021/acs.jctc.3c01193] [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: 03/20/2024]
Abstract
Interplay between divalent cations (Mg2+ and Ca2+) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT30 and dA30, with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg2+ and Ca2+, we find that Mg2+ binds to the minor grooves and phosphate groups, whereas Ca2+ binds specifically to the minor groove. Both Mg2+ and Ca2+ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA30 conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT30 populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.
Collapse
Affiliation(s)
- Balaka Mondal
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - Debayan Chakraborty
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - Naoto Hori
- School of Pharmacy, University of Nottingham, Nottingham NG7 2RD, U.K
| | - Hung T Nguyen
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - D Thirumalai
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
- Department of Physics, The University of Texas, Austin, Texas 78712, United States
| |
Collapse
|
4
|
Pietrek LM, Stelzl LS, Hummer G. Hierarchical Assembly of Single-Stranded RNA. J Chem Theory Comput 2024; 20:2246-2260. [PMID: 38361440 PMCID: PMC10938505 DOI: 10.1021/acs.jctc.3c01049] [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/22/2023] [Revised: 12/09/2023] [Accepted: 01/25/2024] [Indexed: 02/17/2024]
Abstract
Single-stranded RNA (ssRNA) plays a major role in the flow of genetic information-most notably, in the form of messenger RNA (mRNA)-and in the regulation of biological processes. The highly dynamic nature of chains of unpaired nucleobases challenges structural characterizations of ssRNA by experiments or molecular dynamics (MD) simulations alike. Here, we use hierarchical chain growth (HCG) to construct ensembles of ssRNA chains. HCG assembles the structures of protein and nucleic acid chains from fragment libraries created by MD simulations. Applied to homo- and heteropolymeric ssRNAs of different lengths, we find that HCG produces structural ensembles that overall are in good agreement with diverse experiments, including nuclear magnetic resonance (NMR), small-angle X-ray scattering (SAXS), and single-molecule Förster resonance energy transfer (FRET). The agreement can be further improved by ensemble refinement using Bayesian inference of ensembles (BioEn). HCG can also be used to assemble RNA structures that combine base-paired and base-unpaired regions, as illustrated for the 5' untranslated region (UTR) of SARS-CoV-2 RNA.
Collapse
Affiliation(s)
- Lisa M. Pietrek
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Lukas S. Stelzl
- Faculty
of Biology, Johannes Gutenberg University
Mainz, Gresemundweg 2, 55128 Mainz, Germany
- KOMET
1, Institute of Physics, Johannes Gutenberg
University Mainz, 55099 Mainz, Germany
- Institute
of Molecular Biology (IMB), 55128 Mainz, Germany
| | - Gerhard Hummer
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
- Institute
for Biophysics, Goethe University, Max-von-Laue-Straße 9, 60438 Frankfurt am Main, Germany
| |
Collapse
|
5
|
Rinaldi S, Moroni E, Rozza R, Magistrato A. Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation. J Chem Theory Comput 2024; 20:993-1018. [PMID: 38287883 DOI: 10.1021/acs.jctc.3c01239] [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: 01/31/2024]
Abstract
Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.
Collapse
Affiliation(s)
- Silvia Rinaldi
- National Research Council of Italy (CNR) - Institute of Chemistry of OrganoMetallic Compounds (ICCOM), c/o Area di Ricerca CNR di Firenze Via Madonna del Piano 10, 50019 Sesto Fiorentino, Florence, Italy
| | - Elisabetta Moroni
- National Research Council of Italy (CNR) - Institute of Chemical Sciences and Technologies (SCITEC), via Mario Bianco 9, 20131 Milano, Italy
| | - Riccardo Rozza
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| | - Alessandra Magistrato
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| |
Collapse
|
6
|
Ibrahim M, Wenzel C, Lallemang M, Balzer BN, Schwierz N. Adsorbing DNA to Mica by Cations: Influence of Valency and Ion Type. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:15553-15562. [PMID: 37877163 DOI: 10.1021/acs.langmuir.3c01835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2023]
Abstract
Ion-mediated attraction between DNA and mica plays a crucial role in biotechnological applications and molecular imaging. Here, we combine molecular dynamics simulations and single-molecule atomic force microscopy experiments to characterize the detachment forces of single-stranded DNA at mica surfaces mediated by the metal cations Li+, Na+, K+, Cs+, Mg2+, and Ca2+. Ion-specific adsorption at the mica/water interface compensates (Li+ and Na+) or overcompensates (K+, Cs+, Mg2+, and Ca2+) the bare negative surface charge of mica. In addition, direct and water-mediated contacts are formed between the ions, the phosphate oxygens of DNA, and mica. The different contact types give rise to low- and high-force pathways and a broad distribution of detachment forces. Weakly hydrated ions, such as Cs+ and water-mediated contacts, lead to low detachment forces and high mobility of the DNA on the surface. Direct ion-DNA or ion-surface contacts lead to significantly higher forces. The comprehensive view gained from our combined approach allows us to highlight the most promising cations for imaging in physiological conditions: K+, which overcompensates the negative mica charge and induces long-ranged attractions. Mg2+ and Ca2+, which form a few specific and long-lived contacts to bind DNA with high affinity.
Collapse
Affiliation(s)
- Mohd Ibrahim
- Institute of Physics, University of Augsburg, Universitätsstraße 1, 86159 Augsburg, Germany
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Christiane Wenzel
- Institute of Physical Chemistry, University of Freiburg, Albertstraße 21, 79104 Freiburg, Germany
- Cluster of Excellence livMatS @ FIT-Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
| | - Max Lallemang
- Institute of Physical Chemistry, University of Freiburg, Albertstraße 21, 79104 Freiburg, Germany
- Cluster of Excellence livMatS @ FIT-Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
| | - Bizan N Balzer
- Institute of Physical Chemistry, University of Freiburg, Albertstraße 21, 79104 Freiburg, Germany
- Cluster of Excellence livMatS @ FIT-Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
- Freiburg Materials Research Center (FMF), University of Freiburg, Stefan-Meier-Str. 21, 79104 Freiburg, Germany
| | - Nadine Schwierz
- Institute of Physics, University of Augsburg, Universitätsstraße 1, 86159 Augsburg, Germany
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
7
|
Schofield P, Taylor AI, Rihon J, Peña Martinez CD, Zinn S, Mattelaer CA, Jackson J, Dhaliwal G, Schepers G, Herdewijn P, Lescrinier E, Christ D, Holliger P. Characterization of an HNA aptamer suggests a non-canonical G-quadruplex motif. Nucleic Acids Res 2023; 51:7736-7748. [PMID: 37439359 PMCID: PMC10450178 DOI: 10.1093/nar/gkad592] [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: 03/16/2021] [Revised: 06/09/2023] [Accepted: 07/05/2023] [Indexed: 07/14/2023] Open
Abstract
Nucleic acids not only form the basis of heredity, but are increasingly a source of novel nano-structures, -devices and drugs. This has spurred the development of chemically modified alternatives (xeno nucleic acids (XNAs)) comprising chemical configurations not found in nature to extend their chemical and functional scope. XNAs can be evolved into ligands (XNA aptamers) that bind their targets with high affinity and specificity. However, detailed investigations into structural and functional aspects of XNA aptamers have been limited. Here we describe a detailed structure-function analysis of LYS-S8-19, a 1',5'-anhydrohexitol nucleic acid (HNA) aptamer to hen egg-white lysozyme (HEL). Mapping of the aptamer interaction interface with its cognate HEL target antigen revealed interaction epitopes, affinities, kinetics and hot-spots of binding energy similar to protein ligands such as anti-HEL-nanobodies. Truncation analysis and molecular dynamics (MD) simulations suggest that the HNA aptamer core motif folds into a novel and not previously observed HNA tertiary structure, comprising non-canonical hT-hA-hT/hT-hT-hT triplet and hG4-quadruplex structures, consistent with its recognition by two different G4-specific antibodies.
Collapse
Affiliation(s)
- Peter Schofield
- Garvan Institute of Medical Research, Darlinghurst, Sydney, NSW 2010, Australia
- St Vincent's Clinical School, Faculty of Medicine, University of New South Wales, Kensington, Sydney, NSW 2010, Australia
| | - Alexander I Taylor
- MRC Laboratory of Molecular Biology, Cambridge CB2 2QH, UK
- Cambridge Institute of Therapeutic Immunology & Infectious Disease (CITIID), University of Cambridge, Cambridge CB2 0AW, UK
| | - Jérôme Rihon
- Rega Institute, Laboratory of Medicinal Chemistry, Katholieke Universiteit Leuven, Herestraat 49, B 3000, Leuven, Belgium
| | - Cristian D Peña Martinez
- Garvan Institute of Medical Research, Darlinghurst, Sydney, NSW 2010, Australia
- St Vincent's Clinical School, Faculty of Medicine, University of New South Wales, Kensington, Sydney, NSW 2010, Australia
| | - Sacha Zinn
- Garvan Institute of Medical Research, Darlinghurst, Sydney, NSW 2010, Australia
- St Vincent's Clinical School, Faculty of Medicine, University of New South Wales, Kensington, Sydney, NSW 2010, Australia
| | - Charles-Alexandre Mattelaer
- Rega Institute, Laboratory of Medicinal Chemistry, Katholieke Universiteit Leuven, Herestraat 49, B 3000, Leuven, Belgium
| | - Jennifer Jackson
- Garvan Institute of Medical Research, Darlinghurst, Sydney, NSW 2010, Australia
| | - Gurpreet Dhaliwal
- Cambridge Institute of Therapeutic Immunology & Infectious Disease (CITIID), University of Cambridge, Cambridge CB2 0AW, UK
| | - Guy Schepers
- Rega Institute, Laboratory of Medicinal Chemistry, Katholieke Universiteit Leuven, Herestraat 49, B 3000, Leuven, Belgium
| | - Piet Herdewijn
- Rega Institute, Laboratory of Medicinal Chemistry, Katholieke Universiteit Leuven, Herestraat 49, B 3000, Leuven, Belgium
| | - Eveline Lescrinier
- Rega Institute, Laboratory of Medicinal Chemistry, Katholieke Universiteit Leuven, Herestraat 49, B 3000, Leuven, Belgium
| | - Daniel Christ
- Garvan Institute of Medical Research, Darlinghurst, Sydney, NSW 2010, Australia
- St Vincent's Clinical School, Faculty of Medicine, University of New South Wales, Kensington, Sydney, NSW 2010, Australia
| | | |
Collapse
|
8
|
Mu ZC, Tan YL, Liu J, Zhang BG, Shi YZ. Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding. Molecules 2023; 28:4833. [PMID: 37375388 DOI: 10.3390/molecules28124833] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 06/11/2023] [Accepted: 06/14/2023] [Indexed: 06/29/2023] Open
Abstract
DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.
Collapse
Affiliation(s)
- Zi-Chun Mu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| |
Collapse
|
9
|
Sule K, Anikovskiy M, Prenner EJ. Lipid Structure Determines the Differential Impact of Single Metal Additions and Binary Mixtures of Manganese, Calcium and Magnesium on Membrane Fluidity and Liposome Size. Int J Mol Sci 2023; 24:ijms24021066. [PMID: 36674581 PMCID: PMC9860990 DOI: 10.3390/ijms24021066] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Revised: 12/23/2022] [Accepted: 01/02/2023] [Indexed: 01/08/2023] Open
Abstract
Unilamellar vesicles of the biologically relevant lipids phosphatidic acid (PA) and phosphatidylserine (PS) with fully saturated (DM-) or partly unsaturated (PO-) acyl side chains were exposed to Ca, Mn and Mg in single metal additions; in equimolar mixtures or by sequential additions of one metal at a time. Laurdan generalized polarization measured the membrane fluidity, while dynamic light scattering reported liposome size changes complemented by zeta potential. All metals induced membrane rigidity and increased liposome sizes across all systems. Mn had the strongest effect overall, but Mg was comparable for DMPS. Lipid side chain architecture was important as GP values for binary mixtures were higher than expected from the sum of values for single additions added to POPS but smaller for DMPS. Sequential additions were predominantly different for Ca:Mg mixtures. Mn induced the strongest increase of liposome size in saturated lipids whereas Ca effects dominated unsaturated matrices. Binary additions induced larger sizes than the sum of single additions for POPS, but much lower changes in DMPA. The order of addition was relevant for PS systems. Thus, lipid structure determines metal effects, but their impact is modulated by other ions. Thus, metal effects may differ with the local lipid architecture and metal concentrations within cells.
Collapse
Affiliation(s)
- Kevin Sule
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Max Anikovskiy
- Department of Chemistry, Nanoscience Program, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Elmar J. Prenner
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
- Correspondence:
| |
Collapse
|
10
|
Cruz-León S, Schwierz N. RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations. J Phys Chem B 2022; 126:8646-8654. [PMID: 36260822 PMCID: PMC9639116 DOI: 10.1021/acs.jpcb.2c04488] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 μs of simulation time. For the mono- and divalent cations Li+, Na+, K+, Cs+, Ca2+, Sr2+, and Ba2+, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Str. 3, 60438Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Str. 3, 60438Frankfurt am Main, Germany,Institute
of Physics, University of Augsburg, Universitätsstraße 1, 86159Augsburg, Germany,E-mail:
| |
Collapse
|
11
|
Potential effects of metal ion induced two-state allostery on the regulatory mechanism of add adenine riboswitch. Commun Biol 2022; 5:1120. [PMID: 36273041 PMCID: PMC9588036 DOI: 10.1038/s42003-022-04096-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Accepted: 10/11/2022] [Indexed: 11/09/2022] Open
Abstract
Riboswitches normally regulate gene expression through structural changes in response to the specific binding of cellular metabolites or metal ions. Taking add adenine riboswitch as an example, we explore the influences of metal ions (especially for K+ and Mg2+ ions) on the structure and dynamics of riboswitch aptamer (with and without ligand) by using molecular dynamic (MD) simulations. Our results show that a two-state transition marked by the structural deformation at the connection of J12 and P1 (CJ12-P1) is not only related to the binding of cognate ligands, but also strongly coupled with the change of metal ion environments. Moreover, the deformation of the structure at CJ12-P1 can be transmitted to P1 directly connected to the expression platform in multiple ways, which will affect the structure and stability of P1 to varying degrees, and finally change the regulation state of this riboswitch. Molecular dynamic simulations are employed to assess the influence of metal ions on riboswitch structure and dynamics, suggesting a conformational control of riboswitch aptamers by metal ions before ligand binding.
Collapse
|
12
|
Borišek J, Aupič J, Magistrato A. Establishing the catalytic and regulatory mechanism of
RNA
‐based machineries. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Jure Borišek
- Theory Department National Institute of Chemistry Ljubljana Slovenia
| | | | | |
Collapse
|
13
|
Transformation characteristics of A-DNA in salt solution revealed through molecular dynamics simulations. Biophys Chem 2022; 288:106845. [DOI: 10.1016/j.bpc.2022.106845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 05/25/2022] [Accepted: 06/06/2022] [Indexed: 11/17/2022]
|
14
|
Cruz-León S, Vanderlinden W, Müller P, Forster T, Staudt G, Lin YY, Lipfert J, Schwierz N. Twisting DNA by salt. Nucleic Acids Res 2022; 50:5726-5738. [PMID: 35640616 PMCID: PMC9177979 DOI: 10.1093/nar/gkac445] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2021] [Revised: 05/06/2022] [Accepted: 05/10/2022] [Indexed: 12/11/2022] Open
Abstract
The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Willem Vanderlinden
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Peter Müller
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Tobias Forster
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Georgina Staudt
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Yi-Yun Lin
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Jan Lipfert
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
15
|
Puyo-Fourtine J, Juillé M, Hénin J, Clavaguéra C, Duboué-Dijon E. Consistent Picture of Phosphate-Divalent Cation Binding from Models with Implicit and Explicit Electronic Polarization. J Phys Chem B 2022; 126:4022-4034. [PMID: 35608554 DOI: 10.1021/acs.jpcb.2c01158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The binding of divalent cations to the ubiquitous phosphate group is essential for a number of key biological processes, such as DNA compaction, RNA folding, or interactions of some proteins with membranes. Yet, probing their binding sites, modes, and associated binding free energy is a challenge for both experiments and simulations. In simulations, standard force fields strongly overestimate the interaction between phosphate groups and divalent cations. Here, we examine how different strategies to include electronic polarization effects in force fields─implicitly, through the use of scaled charges or pair-specific Lennard-Jones parameters, or explicitly, with the polarizable force fields Drude and AMOEBA─capture the interactions of a model phosphate compound, dimethyl phosphate, with calcium and magnesium divalent cations. We show that both implicit and explicit approaches, when carefully parameterized, are successful in capturing the overall binding free energy and that common trends emerge from the comparison of different simulation approaches. Overall, the binding is very moderate, slightly weaker for Ca2+ than Mg2+, and the solvent-shared ion pair is slightly more stable than the contact monodentate ion pair. The bidentate ion pair is higher in energy (or even fully unstable for Mg2+). Our results thus suggest practical ways to capture the divalent cations with biomolecular phosphate groups in complex biochemical systems. In particular, the computational efficiency of implicit models makes them ideally suited for large-scale simulations of biological assemblies, with improved accuracy compared to state-of-the-art fixed-charge force fields.
Collapse
Affiliation(s)
- Julie Puyo-Fourtine
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Marie Juillé
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Jérôme Hénin
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Elise Duboué-Dijon
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| |
Collapse
|
16
|
Grotz KK, Schwierz N. Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E. J Chem Phys 2022; 156:114501. [DOI: 10.1063/5.0087292] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.
Collapse
Affiliation(s)
- Kara K. Grotz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| | - Nadine Schwierz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| |
Collapse
|
17
|
Neumann J, Schwierz N. Artificial Intelligence Resolves Kinetic Pathways of Magnesium Binding to RNA. J Chem Theory Comput 2022; 18:1202-1212. [PMID: 35084846 PMCID: PMC8830046 DOI: 10.1021/acs.jctc.1c00752] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Magnesium is an indispensable cofactor in countless vital processes. In order to understand its functional role, the characterization of the binding pathways to biomolecules such as RNA is crucial. Despite the importance, a molecular description is still lacking since the transition from the water-mediated outer-sphere to the direct inner-sphere coordination is on the millisecond time scale and therefore out of reach for conventional simulation techniques. To fill this gap, we use transition path sampling to resolve the binding pathways and to elucidate the role of the solvent in the binding process. The results reveal that the molecular void provoked by the leaving phosphate oxygen of the RNA is immediately filled by an entering water molecule. In addition, water molecules from the first and second hydration shell couple to the concerted exchange. To capture the intimate solute-solvent coupling, we perform a committor analysis as the basis for a machine learning algorithm that derives the optimal deep learning model from thousands of scanned architectures using hyperparameter tuning. The results reveal that the properly optimized deep network architecture recognizes the important solvent structures, extracts the relevant information, and predicts the commitment probability with high accuracy. Our results provide detailed insights into the solute-solvent coupling which is ubiquitous for kosmotropic ions and governs a large variety of biochemical reactions in aqueous solutions.
Collapse
Affiliation(s)
- Jan Neumann
- Allianz Global Investors GmbH, Bockenheimer Landstrasse 42, 60323 Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, 60438 Frankfurt am Main, Germany
| |
Collapse
|
18
|
Grotz KK, Schwierz N. Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D. J Chem Theory Comput 2022; 18:526-537. [PMID: 34881568 PMCID: PMC8757469 DOI: 10.1021/acs.jctc.1c00791] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Indexed: 12/14/2022]
Abstract
Magnesium is essential in many vital processes. To correctly describe Mg2+ in physiological processes by molecular dynamics simulations, accurate force fields are fundamental. Despite the importance, force fields based on the commonly used 12-6 Lennard-Jones potential showed significant shortcomings. Recently progress was made by an optimization procedure that implicitly accounts for polarizability. The resulting microMg and nanoMg force fields (J. Chem. Theory Comput. 2021, 17, 2530-2540) accurately reproduce a broad range of experimental solution properties and the binding affinity to nucleic acids in TIP3P water. Since countless simulation studies rely on available water models and ion force fields, we here extend the optimization and provide Mg2+ parameters in combination with the SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D water models. For each water model, the Mg2+ force fields reproduce the solvation free energy, the distance to oxygens in the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, and the binding affinity and distance to the phosphate oxygens on nucleic acids. We present two parameter sets: MicroMg yields water exchange on the microsecond time scale and matches the experimental exchange rate. Depending on the water model, nanoMg yields accelerated water exchange in the range of 106 to 108 exchanges per second. The nanoMg parameters can be used to enhance the sampling of binding events, to obtain converged distributions of Mg2+, or to predict ion binding sites in biomolecular simulations. The parameter files are freely available at https://github.com/bio-phys/optimizedMgFFs.
Collapse
Affiliation(s)
- Kara K. Grotz
- Department of Theoretical
Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Nadine Schwierz
- Department of Theoretical
Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| |
Collapse
|