1
|
Mou X, Liu K, He L, Li S. Mechanical response of double-stranded DNA: Bend, twist, and overwind. J Chem Phys 2024; 161:085102. [PMID: 39177087 DOI: 10.1063/5.0216585] [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: 04/30/2024] [Accepted: 08/09/2024] [Indexed: 08/24/2024] Open
Abstract
We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.
Collapse
Affiliation(s)
- Xuankang Mou
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| | - Kai Liu
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| | - Linli He
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| | - Shiben Li
- Department of Physics, Wenzhou University, Wenzhou, Zhejiang 325035, China
| |
Collapse
|
2
|
Kar M, Vogel LT, Chauhan G, Felekyan S, Ausserwöger H, Welsh TJ, Dar F, Kamath AR, Knowles TPJ, Hyman AA, Seidel CAM, Pappu RV. Solutes unmask differences in clustering versus phase separation of FET proteins. Nat Commun 2024; 15:4408. [PMID: 38782886 PMCID: PMC11116469 DOI: 10.1038/s41467-024-48775-3] [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: 08/22/2023] [Accepted: 05/03/2024] [Indexed: 05/25/2024] Open
Abstract
Phase separation and percolation contribute to phase transitions of multivalent macromolecules. Contributions of percolation are evident through the viscoelasticity of condensates and through the formation of heterogeneous distributions of nano- and mesoscale pre-percolation clusters in sub-saturated solutions. Here, we show that clusters formed in sub-saturated solutions of FET (FUS-EWSR1-TAF15) proteins are affected differently by glutamate versus chloride. These differences on the nanoscale, gleaned using a suite of methods deployed across a wide range of protein concentrations, are prevalent and can be unmasked even though the driving forces for phase separation remain unchanged in glutamate versus chloride. Strikingly, differences in anion-mediated interactions that drive clustering saturate on the micron-scale. Beyond this length scale the system separates into coexisting phases. Overall, we find that sequence-encoded interactions, mediated by solution components, make synergistic and distinct contributions to the formation of pre-percolation clusters in sub-saturated solutions, and to the driving forces for phase separation.
Collapse
Affiliation(s)
- Mrityunjoy Kar
- Max Planck Institute of Cell Biology and Genetics, 01307, Dresden, Germany
| | - Laura T Vogel
- Department of Molecular Physical Chemistry, Heinrich Heine University, 40225, Düsseldorf, Germany
| | - Gaurav Chauhan
- Department of Biomedical Engineering and Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Suren Felekyan
- Department of Molecular Physical Chemistry, Heinrich Heine University, 40225, Düsseldorf, Germany
| | - Hannes Ausserwöger
- Centre for Misfolding Diseases, Yusuf Hamied Department of Chemistry, University of Cambridge, CB2 1EW, Cambridge, UK
| | - Timothy J Welsh
- Centre for Misfolding Diseases, Yusuf Hamied Department of Chemistry, University of Cambridge, CB2 1EW, Cambridge, UK
| | - Furqan Dar
- Department of Biomedical Engineering and Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Anjana R Kamath
- Max Planck Institute of Cell Biology and Genetics, 01307, Dresden, Germany
| | - Tuomas P J Knowles
- Centre for Misfolding Diseases, Yusuf Hamied Department of Chemistry, University of Cambridge, CB2 1EW, Cambridge, UK
| | - Anthony A Hyman
- Max Planck Institute of Cell Biology and Genetics, 01307, Dresden, Germany.
| | - Claus A M Seidel
- Department of Molecular Physical Chemistry, Heinrich Heine University, 40225, Düsseldorf, Germany.
| | - Rohit V Pappu
- Department of Biomedical Engineering and Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO, 63130, USA.
| |
Collapse
|
3
|
Minero GA, Møllebjerg A, Thiesen C, Johansen M, Jørgensen N, Birkedal V, Otzen DE, Meyer R. Extracellular G-quadruplexes and Z-DNA protect biofilms from DNase I, and G-quadruplexes form a DNAzyme with peroxidase activity. Nucleic Acids Res 2024; 52:1575-1590. [PMID: 38296834 PMCID: PMC10939358 DOI: 10.1093/nar/gkae034] [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: 06/30/2023] [Revised: 01/04/2024] [Accepted: 01/10/2024] [Indexed: 02/02/2024] Open
Abstract
Many bacteria form biofilms to protect themselves from predators or stressful environmental conditions. In the biofilm, bacteria are embedded in a protective extracellular matrix composed of polysaccharides, proteins and extracellular DNA (eDNA). eDNA most often is released from lysed bacteria or host mammalian cells, and it is the only matrix component most biofilms appear to have in common. However, little is known about the form DNA takes in the extracellular space, and how different non-canonical DNA structures such as Z-DNA or G-quadruplexes might contribute to its function in the biofilm. The aim of this study was to determine if non-canonical DNA structures form in eDNA-rich staphylococcal biofilms, and if these structures protect the biofilm from degradation by nucleases. We grew Staphylococcus epidermidis biofilms in laboratory media supplemented with hemin and NaCl to stabilize secondary DNA structures and visualized their location by immunolabelling and fluorescence microscopy. We furthermore visualized the macroscopic biofilm structure by optical coherence tomography. We developed assays to quantify degradation of Z-DNA and G-quadruplex DNA oligos by different nucleases, and subsequently investigated how these enzymes affected eDNA in the biofilms. Z-DNA and G-quadruplex DNA were abundant in the biofilm matrix, and were often present in a web-like structures. In vitro, the structures did not form in the absence of NaCl or mechanical shaking during biofilm growth, or in bacterial strains deficient in eDNA or exopolysaccharide production. We thus infer that eDNA and polysaccharides interact, leading to non-canonical DNA structures under mechanical stress when stabilized by salt. We also confirmed that G-quadruplex DNA and Z-DNA was present in biofilms from infected implants in a murine implant-associated osteomyelitis model. Mammalian DNase I lacked activity against Z-DNA and G-quadruplex DNA, while Micrococcal nuclease could degrade G-quadruplex DNA and S1 Aspergillus nuclease could degrade Z-DNA. Micrococcal nuclease, which originates from Staphylococcus aureus, may thus be key for dispersal of biofilm in staphylococci. In addition to its structural role, we show for the first time that the eDNA in biofilms forms a DNAzyme with peroxidase-like activity in the presence of hemin. While peroxidases are part of host defenses against pathogens, we now show that biofilms can possess intrinsic peroxidase activity in the extracellular matrix.
Collapse
Affiliation(s)
| | - Andreas Møllebjerg
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000 Aarhus, Denmark
| | - Celine Thiesen
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000 Aarhus, Denmark
| | - Mikkel Illemann Johansen
- Department Infectious Diseases, Aarhus University Hospital, Palle Juul-Jensens bvld 99, 8200 Aarhus N, Denmark
| | - Nis Pedersen Jørgensen
- Department Infectious Diseases, Aarhus University Hospital, Palle Juul-Jensens bvld 99, 8200 Aarhus N, Denmark
| | - Victoria Birkedal
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000 Aarhus, Denmark
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus, Denmark
| | - Daniel Erik Otzen
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000 Aarhus, Denmark
- Department of Molecular Biology and Genetics, Aarhus University, Universitetsbyen 81, 8000 Aarhus, Denmark
| | - Rikke Louise Meyer
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000 Aarhus, Denmark
- Department of Biology, Aarhus University, Ny Munkegade 114, 8000 Aarhus, Denmark
| |
Collapse
|
4
|
Nichols PJ, Krall JB, Henen MA, Welty R, MacFadden A, Vicens Q, Vögeli B. Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate. J Am Chem Soc 2024; 146:677-694. [PMID: 38131335 PMCID: PMC11155437 DOI: 10.1021/jacs.3c10406] [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: 12/23/2023]
Abstract
The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.
Collapse
Affiliation(s)
- Parker J. Nichols
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
| | - Jeffrey B. Krall
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
| | - Morkos A. Henen
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
- Faculty of Pharmacy, Mansoura University, Mansoura, 35516, Egypt
| | - Robb Welty
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
| | - Andrea MacFadden
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
| | - Quentin Vicens
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
- RNA Bioscience Initiative, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
- Present address: Department of Biology and Biochemistry, Center for Nuclear Receptors and Cellular Signaling, University of Houston, Houston, Texas 77204, USA
| | - Beat Vögeli
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado, 80045, USA
| |
Collapse
|
5
|
Bartas M, Slychko K, Červeň J, Pečinka P, Arndt-Jovin DJ, Jovin TM. Extensive Bioinformatics Analyses Reveal a Phylogenetically Conserved Winged Helix (WH) Domain (Zτ) of Topoisomerase IIα, Elucidating Its Very High Affinity for Left-Handed Z-DNA and Suggesting Novel Putative Functions. Int J Mol Sci 2023; 24:10740. [PMID: 37445918 DOI: 10.3390/ijms241310740] [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] [Received: 04/19/2023] [Revised: 06/13/2023] [Accepted: 06/22/2023] [Indexed: 07/15/2023] Open
Abstract
The dynamic processes operating on genomic DNA, such as gene expression and cellular division, lead inexorably to topological challenges in the form of entanglements, catenanes, knots, "bubbles", R-loops, and other outcomes of supercoiling and helical disruption. The resolution of toxic topological stress is the function attributed to DNA topoisomerases. A prominent example is the negative supercoiling (nsc) trailing processive enzymes such as DNA and RNA polymerases. The multiple equilibrium states that nscDNA can adopt by redistribution of helical twist and writhe include the left-handed double-helical conformation known as Z-DNA. Thirty years ago, one of our labs isolated a protein from Drosophila cells and embryos with a 100-fold greater affinity for Z-DNA than for B-DNA, and identified it as topoisomerase II (gene Top2, orthologous to the human UniProt proteins TOP2A and TOP2B). GTP increased the affinity and selectivity for Z-DNA even further and also led to inhibition of the isomerase enzymatic activity. An allosteric mechanism was proposed, in which topoII acts as a Z-DNA-binding protein (ZBP) to stabilize given states of topological (sub)domains and associated multiprotein complexes. We have now explored this possibility by comprehensive bioinformatic analyses of the available protein sequences of topoII representing organisms covering the whole tree of life. Multiple alignment of these sequences revealed an extremely high level of evolutionary conservation, including a winged-helix protein segment, here denoted as Zτ, constituting the putative structural homolog of Zα, the canonical Z-DNA/Z-RNA binding domain previously identified in the interferon-inducible RNA Adenosine-to-Inosine-editing deaminase, ADAR1p150. In contrast to Zα, which is separate from the protein segment responsible for catalysis, Zτ encompasses the active site tyrosine of topoII; a GTP-binding site and a GxxG sequence motif are in close proximity. Quantitative Zτ-Zα similarity comparisons and molecular docking with interaction scoring further supported the "B-Z-topoII hypothesis" and has led to an expanded mechanism for topoII function incorporating the recognition of Z-DNA segments ("Z-flipons") as an inherent and essential element. We further propose that the two Zτ domains of the topoII homodimer exhibit a single-turnover "conformase" activity on given G(ate) B-DNA segments ("Z-flipins"), inducing their transition to the left-handed Z-conformation. Inasmuch as the topoII-Z-DNA complexes are isomerase inactive, we infer that they fulfill important structural roles in key processes such as mitosis. Topoisomerases are preeminent targets of anti-cancer drug discovery, and we anticipate that detailed elucidation of their structural-functional interactions with Z-DNA and GTP will facilitate the design of novel, more potent and selective anti-cancer chemotherapeutic agents.
Collapse
Affiliation(s)
- Martin Bartas
- Department of Biology and Ecology, University of Ostrava, 710 00 Ostrava, Czech Republic
| | - Kristyna Slychko
- Department of Biology and Ecology, University of Ostrava, 710 00 Ostrava, Czech Republic
| | - Jiří Červeň
- Department of Biology and Ecology, University of Ostrava, 710 00 Ostrava, Czech Republic
| | - Petr Pečinka
- Department of Biology and Ecology, University of Ostrava, 710 00 Ostrava, Czech Republic
| | - Donna J Arndt-Jovin
- Emeritus Laboratory of Cellular Dynamics, Max Planck Institute for Multidisciplinary Sciences, 37077 Göttingen, Germany
| | - Thomas M Jovin
- Emeritus Laboratory of Cellular Dynamics, Max Planck Institute for Multidisciplinary Sciences, 37077 Göttingen, Germany
| |
Collapse
|
6
|
Duan M, Li Y, Zhang F, Huang Q. Assessing B-Z DNA Transitions in Solutions via Infrared Spectroscopy. Biomolecules 2023; 13:964. [PMID: 37371544 DOI: 10.3390/biom13060964] [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] [Received: 04/16/2023] [Revised: 05/25/2023] [Accepted: 06/06/2023] [Indexed: 06/29/2023] Open
Abstract
Z-DNA refers to the left-handed double-helix DNA that has attracted much attention because of its association with some specific biological functions. However, because of its low content and unstable conformation, Z-DNA is normally difficult to observe or identify. Up to now, there has been a lack of unified or standard analytical methods among diverse techniques for probing Z-DNA and its transformation conveniently. In this work, NaCl, MgCl2, and ethanol were utilized to induce d(GC)8 from B-DNA to Z-DNA in vitro, and Fourier transform infrared (FTIR) spectroscopy was employed to monitor the transformation of Z-DNA under different induction conditions. The structural changes during the transformation process were carefully examined, and the DNA chirality alterations were validated by the circular dichroism (CD) measurements. The Z-DNA characteristic signals in the 1450 cm-1-900 cm-1 region of the d(GC)8 infrared (IR) spectrum were observed, which include the peaks at 1320 cm-1, 1125 cm-1 and 925 cm-1, respectively. The intensity ratios of A1320/A970, A1125/A970, and A925/A970 increased with Z-DNA content in the transition process. Furthermore, compared with the CD spectra, the IR spectra showed higher sensitivity to Z-DNA, providing more information about the molecular structure change of DNA. Therefore, this study has established a more reliable FTIR analytical approach to assess BZ DNA conformational changes in solutions, which may help the understanding of the Z-DNA transition mechanism and promote the study of Z-DNA functions in biological systems.
Collapse
Affiliation(s)
- Mengmeng Duan
- Henan Key Laboratory of Ion-Beam Bioengineering, School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450052, China
| | - Yalin Li
- School of Food and Biological Engineering, Henan University of Animal Husbandry and Economy, Zhengzhou 450047, China
| | - Fengqiu Zhang
- Henan Key Laboratory of Ion-Beam Bioengineering, School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450052, China
| | - Qing Huang
- CAS Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Institute of Intelligent Machines, Hefei Institutes of Physical Sciences, Chinese Academy of Sciences, Hefei 230031, China
- Science Island Branch of Graduate School, University of Science and Technology, Hefei 230026, China
| |
Collapse
|
7
|
Luengo-Márquez J, Zalvide-Pombo J, Pérez R, Assenza S. Force-dependent elasticity of nucleic acids. NANOSCALE 2023; 15:6738-6744. [PMID: 36942727 DOI: 10.1039/d2nr06324g] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The functioning of double-stranded (ds) nucleic acids (NAs) in cellular processes is strongly mediated by their elastic response. These processes involve proteins that interact with dsDNA or dsRNA and distort their structures. The perturbation of the elasticity of NAs arising from these deformations is not properly considered by most theoretical frameworks. In this work, we introduce a novel method to assess the impact of mechanical stress on the elastic response of dsDNA and dsRNA through the analysis of the fluctuations of the double helix. Application of this approach to atomistic simulations reveals qualitative differences in the force dependence of the mechanical properties of dsDNA with respect to those of dsRNA, which we relate to structural features of these molecules by means of physically-sound minimalistic models.
Collapse
Affiliation(s)
- Juan Luengo-Márquez
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
- Instituto Nicolás Cabrera, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Juan Zalvide-Pombo
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
| | - Rubén Pérez
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
- Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain.
| | - Salvatore Assenza
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
- Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain.
- Instituto Nicolás Cabrera, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|
8
|
Nichols PJ, Krall JB, Henen MA, Vögeli B, Vicens Q. Z-RNA biology: a central role in the innate immune response? RNA (NEW YORK, N.Y.) 2023; 29:273-281. [PMID: 36596670 PMCID: PMC9945438 DOI: 10.1261/rna.079429.122] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Z-RNA is a higher-energy, left-handed conformation of RNA, whose function has remained elusive. A growing body of work alludes to regulatory roles for Z-RNA in the immune response. Here, we review how Z-RNA features present in cellular RNAs-especially containing retroelements-could be recognized by a family of winged helix proteins, with an impact on host defense. We also discuss how mutations to specific Z-contacting amino acids disrupt their ability to stabilize Z-RNA, resulting in functional losses. We end by highlighting knowledge gaps in the field, which, if addressed, would significantly advance this active area of research.
Collapse
Affiliation(s)
- Parker J Nichols
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Jeffrey B Krall
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Morkos A Henen
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
- Faculty of Pharmacy, Mansoura University, Mansoura, 35516, Egypt
| | - Beat Vögeli
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
- RNA Bioscience Initiative, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Quentin Vicens
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
- RNA Bioscience Initiative, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| |
Collapse
|