Next Article in Journal
Alpha-Synuclein and Lipids: The Elephant in the Room?
Next Article in Special Issue
Angiotensin II-Induced Long Non-Coding RNA Alivec Regulates Chondrogenesis in Vascular Smooth Muscle Cells
Previous Article in Journal
Dynamic Expression of Membrane Type 1-Matrix Metalloproteinase (Mt1-mmp/Mmp14) in the Mouse Embryo
Previous Article in Special Issue
Mechanisms of Fibroblast Activation and Myocardial Fibrosis: Lessons Learned from FB-Specific Conditional Mouse Models
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 10
Issue 9
10.3390/cells10092449
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Effects of Ractopamine on Cardiovascular, Respiratory, and Locomotory Physiology in Animal Model Zebrafish Larvae

by
Kumail Abbas
1,†,
Ferry Saputra
1,†,
Michael Edbert Suryanto
1,†,
Yu-Heng Lai
2,
Jong-Chin Huang
3,
Wen-Hao Yu
4,
Kelvin H.-C. Chen
3,*,
Ying-Ting Lin
4,5,* and
Chung-Der Hsiao
1,6,7,8,*
1
Department of Bioscience Technology, Chung Yuan Christian University, Chung-Li, Taoyuan City 3020314, Taiwan
2
Department of Chemistry, Chinese Culture University, Taipei 11114, Taiwan
3
Department of Applied Chemistry, National Pingtung University, Pingtung 900391, Taiwan
4
Department of Biotechnology, College of Life Science, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
5
Drug Development and Value Creation Research Center, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
6
Department of Chemistry, Chung Yuan Christian University, Chung-Li, Taoyuan City 320314, Taiwan
7
Center for Nanotechnology, Chung Yuan Christian University, Chung-Li, Taoyuan City 320314, Taiwan
8
Research Center for Aquatic Toxicology and Pharmacology, Chung Yuan Christian University, Chung-Li, Taoyuan City 320314, Taiwan
*
Authors to whom correspondence should be addressed.
†
Authors contributed equally to this work.
Cells 2021, 10(9), 2449; https://doi.org/10.3390/cells10092449
Submission received: 6 July 2021 / Revised: 8 September 2021 / Accepted: 10 September 2021 / Published: 17 September 2021
(This article belongs to the Collection Cardiovascular Disease: From Molecular and Cellular Mechanisms to Therapeutic Opportunities)
Browse Figures
Review Reports Versions Notes

Abstract

:
Ractopamine (RAC) is a beta-adrenoceptor agonist that is used to promote lean and increased food conversion efficiency in livestock. This compound has been considered to be causing behavioral and physiological alterations in livestock like pig. Few studies have addressed the potential non-target effect of RAC in aquatic animals. In this study, we aimed to explore the potential physiological response after acute RAC exposure in zebrafish by evaluating multiple endpoints like locomotor activity, oxygen consumption, and cardiovascular performance. Zebrafish larvae were subjected to waterborne RAC exposure at 0.1, 1, 2, 4, or 8 ppm for 24 h, and the corresponding cardiovascular, respiratory, and locomotion activities were monitored and quantified. In addition, we also performed in silico molecular docking for RAC with 10 zebrafish endogenous β-adrenergic receptors to elucidate the potential acting mechanism of RAC. Results show RAC administration can significantly boost locomotor activity, cardiac performance, oxygen consumption, and blood flow rate, but without affecting the cardiac rhythm regularity in zebrafish embryos. Based on structure-based flexible molecular docking, RAC display similar binding affinity to all ten subtypes of endogenous β-adrenergic receptors, from adra1aa to adra2db, which are equivalent to the human one. This result suggests RAC might act as high potency and broad spectrum β-adrenergic receptors agonist on boosting the locomotor activity, cardiac performance, and oxygen consumption in zebrafish. To validate our results, we co-incubated a well-known β-blocker of propranolol (PROP) with RAC. PROP exposure tends to minimize the locomotor hyperactivity, high oxygen consumption, and cardiac rate in zebrafish larvae. In silico structure-based molecular simulation and binding affinity tests show PROP has an overall lower binding affinity than RAC. Taken together, our studies provide solid in vivo evidence to support that RAC plays crucial roles on modulating cardiovascular, respiratory, and locomotory physiology in zebrafish for the first time. In addition, the versatile functions of RAC as β-agonist possibly mediated via receptor competition with PROP as β-antagonist.

Graphical Abstract

1. Introduction

The excessive usage of drugs in veterinary medicine affects non-targeted aquatic animals, which are also effecting human population through drinking water [1,2]. These drugs are excreted through manure, which is being used as fertilizers, thus polluting aquatic life as well [3]. One widely used drug is ractopamine (RAC), a synthetic β-adrenergic agonist, used as feed additive for increased effectiveness of feed, increased growth and muscle leanness, also termed as ‘lean meat agent’ [4,5]. β-adrenergic receptors are classified into β1, β2, and β3. β1- and β2-adrenergic receptors being most abundant and expressed in heart, kidney, lungs, and blood vessels, while β3-adrenergic receptors are found in the adipose tissue and a very few in the heart [6]. Adrenergic receptors help in modulating several metabolic functions in fish regarding oxygen uptake, cardiac rate, vascular resistance, and hemoglobin–oxygen binding affinity [7,8].
RAC modulates the metabolism and redirects the nutrients to muscles from the adipose tissue, increases lipolysis and protein synthesis [9,10]. The U.S. Food and Drug Administration (FDA) authorized RAC in 2000, and it has been used in selected food animals as feed additive in different countries worldwide. Although RAC is currently banned in 160 countries, including the European Union members [11] and Taiwan [12], use in livestock is still in practice in some countries, including USA, Canada, Japan, South Africa, and Mexico [13].
Over 88% of ingested RAC is excreted through animal feces [14], that is disposed of without treatment and becomes contaminated in subsurface, brooks, rivers, and other territorial waters [15]. Different concentrations of RAC were detected in groundwater near piggery (0.054 μg/L) [15], and agricultural catchment area (1.3 × 10−5 to 5.4 × 10−4 μg/L) [16]. RAC was also detected in sewerage ponds near pig farms in variable amounts ranging from 0.134 to 0.524 μg/L [15] and 0.138 μg/L to 30 μg/L [17]. In a previous investigative study, in hospital sewage and rivers in Taiwan, four β-agonists were detected, while 70% of the accumulated specimens were RAC [18].
However, recent studies demonstrated that the use of RAC may affect heart, thyroid, urethra, and prostate [19], and also causes cardiac toxicity [20]. Few studies have been performed on RAC exposure to non-target aquatic animals. It has been revealed that chronic exposure of 44 days, disrupts the endocrine system in Japanese Medaka (Oryzias latipes), exposed to different doses of RAC between 5 and 625 μg/L [21]. While chronic exposure of different concentrations of RAC in adult zebrafish (Danio rerio) at 0.1, 0.2, 0.85, 8.5, and 85 μg/L induced behavioral changes and oxidative stress [22], and increases the heart rate of zebrafish embryos and larvae [23].
In silico structure-based molecular simulation imparts a great and helpful method to demonstrate a possible ligand–receptor binding location. Due to tremendous advancements in computational scope and efficiency in recent years, molecular docking has become a valuable and accessible method to show the molecular interactivity between the ligand and its receptor at atomic resolution [24]. This method is exceptionally crucial for the academic community in drug discovery and toxicity forecasting. The consequences of the technique can come up with advice on searching out drugs and designing new concepts with more significant binding tendencies [25]. To see how RAC can bind to zebrafish endogenous β-adrenergic receptors, we employed a dry lab molecular simulation to anticipate further reliable receptor structure supporting their ligand–receptor interactions at an atomic scale. To construct a 3D receptor structure model, firstly, we fetched the required sequence from UniProt (https://www.uniprot.org/; accessed on 9 March 2021) and then spotted the most comparative corresponding structure as a template. By doing the sequence alignment between the receptor and template, homology modeling [26] builds up 3D protein structure models to see whether RAC can bind to the constructed receptor. We optimized the RAC structure, explored for probable binding cavities of the receptor, and performed molecular docking for the ligand–receptor binding situation.
Aquatic and environmental toxicity has been a great concern for the last decade. Due to excessive globalization and industrial growth, aquatic life is being affected indirectly. The use of zebrafish as an animal model has wide advantages, including transparency of embryo and larva, high throughput examination, short experimental duration, cost-effective, genetic resemblance with humans (around 85%), less quantity of compound is required, and it is also been approved by the U.S FDA and European Medical Agency (EMEA) for harmful and safety evaluation for Investigative New Drug (IND) compliance [27]. The objective of this study is to explore the potential physiological impacts of the feed additive RAC that is being excessively used in livestock and aquatic industry as a repartitioning agent, which in turn causes some physiological and behavioral alterations in non-target animals and ultimately humans. Although RAC is banned in different countries, its misuse is still a big concern, especially in Taiwan.

2. Materials and Methods

2.1. Housing and Breeding of Zebrafish

AB strain of zebrafish stock were acquired from the Taiwan Zebrafish Core Facility at Academia Sinica (TZCAS) and housed in laboratory for one month prior to the experiment under recirculating water supply at 26.5 °C and a 10/14 dark–light cycle. We fed artemia and commercial feed to fish twice. Female to male ratio in accordance to 1:2 were kept in a breeding tank for the duration of a night by separating them with a glassy separator following our previous published study [28,29]. The barrier was removed the next morning. After collection and sanitization, the embryos, were then put down at 28 °C in methylene blue water. Until the third day of post-fertilization, all the expired larvae were being taken out and methylene blue water medium was replaced daily. All of these zebrafish experimental protocols and principles were authorized by the Chung Yuan Christian University animal care and welfare committee (Number: CYCU109001, issue date: 20 January 2020).

2.2. Chemical Preparation

RAC was purchased from ANPEL Laboratory Technologies Inc., Shanghai, China and was prepared as 1000 ppm for stock solution using ddH2O. Then stock solution was further diluted into desired concentrations using ddH2O at the time of exposure. PROP was purchased from Sigma-Aldrich, Inc., St Louis, Missouri, USA with >99% purity and 1000 ppm stock solution was prepared using ddH2O, which then was further diluted into desired concentrations using ddH2O.

2.3. Automated Behavioral Analysis Using ZebraBox

In this study, we used a high-throughput monitoring enclosure named ‘ZebraBox’ developed by Viewpoint Company (ViewPoint 3.22.3.85, ViewPoint Life Sciences, Inc., Civrieux, France, 2014 (http://www.viewpointlifesciences.com; accessed on 9 June 2021) to analyze zebrafish behaviors. After 96 h post-fertilization (hpf), 48 zebrafish larvae were collected in a 9-cm petri dish for each control and RAC treated groups. The number of experimental animals was based on prior studies [29]. Zebrafish were exposed to RAC at 0.1, 1, 2, 4, and 8 ppm concentrations and were labelled as the treatment group and a control group kept under the same conditions for a 24-h incubation period. After the 24-h incubation period, each group of zebrafish larvae were transferred to 48-well plate for locomotion quantification. Acclimation time of almost 1 h was given before the 48-well plate was kept on the ZebraBox machine plate. After placing the 48-well plate in the ZebraBox, additionally a 10-min acclimation time before proceeding the procedure was given.
Video tracking analysis was performed in a 48-well plate individually for each group of both zebrafish under the infrared illuminating system of the ZebraBox, which has the ability of high-quality imaging in dark. Video was recorded for 80 min (10 min of light and dark period alternatively) for each group and the same video was then used for Rotation and Burst analysis. Overall locomotion activity was analyzed on the basis of total distance travelled, burst count, and rotation count. Movement categorization as cruising (normal speed), increased body activity (short, powerful, and intermittent activity), and no activity (freeze activity) was made in the previous study [30]. As per our previous protocol [31], we categorized the velocities in three ranges to measure the total distance chronology as: velocities exceeding 20 mm/s were named larger or increased body activity, 0.5 to 20 mm/s was called as cruising activity, and less than 0.5 mm/s as inactivity. Moreover, for total burst count, the simulation of video track parameters was set as burst 20 pixels/s and freeze (no movement). Furthermore, rotation chronology was determined by the clockwise and anti-clockwise movements in millimeters and contemplated a diameter over 2 mm as one rotation and rotation less than 2 mm was abandoned. Moreover, the back angle was set at 60° [29]. The tests were performed in duplicate.

2.4. Zebrafish Larvae Cardiac Physiology and Rhythm Assay

Zebrafish embryos were collected after 48 hpf and expose to RAC for 24 h at different concentration (0.1, 4, and 8 ppm). To keep the position of the larvae, 3% methylcellulose was used right before recording as the mounting medium. Recording was done using a high-speed digital charged coupling device (CCD) (AZ Instrument Corp., Taichung City, Taiwan) mounted on to an inverted microscope (ICX41, Sunny Optical Technology, Zhejiang, China) while a Hoffmann modulation contrast objective lens with 40× magnification and LPlans objective lens with 10× magnification was used to record videos with high resolution and contrast. HiBest Viewer software package (AZ Instrument, Taichung City, Taiwan) was deployed to record at a frame rate of 200 frames per second (fps) for 10 s. Ventricle movement, focused on the cardiac physiology parameters, was recorded using Hoffman modulation contrast objective lens, while LPlans objective lens was used to check the rhythm of the heart chamber. To check the variability of zebrafish heartbeat, Poincare plot plug-in from Origin Pro 2019 software (OriginLab Corporation, Northampton, MA, USA) was used to calculate the sd1 and sd2 of both heart chamber. All cardiac performance and rhythm assay was calculated according to our previously published protocol [32].

2.5. Zebrafish Oxygen Consumption Analysis

The oxygen consumption analysis was conducted using a 24-well plate kept in the Sensor Dish Reader (SDR) and software MicroResp® version 1.0.4 (Loligo Systems, Viborg, Denmark) was used to record the oxygen concentration in the well during experiment time. Approximately 46 zebrafish larvae at 96 hpf were exposed to RAC for 24 h. At experiment time, the larvae were put together with 80 μL of the compound. One zebrafish larva was put on each well, except one well was used as a blank and filled with ddH2O only. The oxygen consumption rate was recorded for 80 min for each group of zebrafish larvae and the experiment was done in triplicate.

2.6. Blood Flow Measurement

The blood flow measurement was done according to our previously published study reported by Santoso et al. [33]. Briefly, about 10 zebrafish larvae were placed under an inverted microscope (ICX41, Sunny Optical Technology, Zhejiang, China) in a 5-cm plastic petri dish and mounted with 3% methylcellulose to prevent animal movement. After that, Hoffmann modulation contrast objective Lens with 40× magnification was used to visualize the blood vessel in dorsal aorta (DA). High-speed CCD camera (AZ Instrument, Taichung City, Taiwan) capable to record at 200 fps was used to record the blood movement for 10 s. The video was further analyzed using “Trackmate” plugin at ImageJ Version 1.53e to get the position of blood cells in every frame [34,35]. The data then processed in Microsoft Excel (Excel version in Microsoft Office Professional Plus 2016 for Windows) in to get the velocity of detected blood cell. The experiment was done in duplication.

2.7. Molecular Docking Analysis

To simulate the binding between the inspected RAC and PROP with ten different subtypes of zebrafish endogenous β-adrenergic receptors, we performed homology modeling and molecular docking in Windows 7 Home Basic on an Asus computer (Intel® CoreTM i7 2.67 GHz processor), using Swiss Model web-server (https://swissmodel.expasy.org/; accessed on 14 March 2021) [36] and Modeller Software v9.20 [37] and flexible CDOCKER [38,39]. Homology modeling first constructs ten three-dimensional protein structures using the sequences of zebrafish endogenous β-adrenergic receptors from UniProt (https://www.uniprot.org; accessed on 9 March 2021) and the template homologous crystal templates from Protein Data Bank (PDB, https://www.rcsb.org; accessed on 9 March 2021). The related protein sequences for adra1aa, adra1ab, adra1ba, adra1bb, adra1d, adra2a, adra2b, adra2c, adra2da, and adra2db are used to locate four homologous crystal templates, 6kux (Crystal structures of the alpha2A adrenergic receptor in complex with an antagonist RSC), 5v54 (Crystal structure of 5-HT1B receptor in complex with methiothepin), 7e32 (Serotonin 1D (5-HT1D) receptor-Gi protein complex), and 7e2z (Aripiprazole-bound serotonin 1A (5-HT1A) receptor-Gi protein complex) by Swiss Model tools. Among four crystal binding sites, we overlapped RAC with the four crystal ligands with the identities of 85.3% (6kux), 79.9% (5v54), 73.7% (7e32), 79.8% (7e2z), and getting the best fit with the 6kux ligand. Thus, 6kux was chosen as the template structure. Best predicting structure models were identified by the two scoring functions, DOPE [40], and GA341 [41]. Flexible molecular docking then did the simulation of the binding between RAC or PROP with β-adrenergic receptors of zebrafish. Note that the 6kux (adra2a, Homo sapiens) binding site was also taken for the flexible docking as comparison. A cavity searching method, eraser algorithm [42], is used to locate the best binding pocket. The partial charges of all atoms within the receptor were determined, and all hydrogens were restored by CHARMm force field [43]. Molecular docking results showing the binding between RAC or PROP with zebrafish endogenous β-adrenergic receptors were reported in the Table 1 of results section.

2.8. Co-incubation of RAC and PROP Experiment

Based on the results, RAC exposure could induce the locomotion, cardiac physiology, respiratory performance, and blood flow velocity alterations in larval zebrafish. Additionally, based on molecular docking analysis, it showed that RAC can interact with β-adrenergic receptors with high affinity. To validate this hypothesis, we evaluated the physiological outcomes of zebrafish after co-incubating with RAC and β-adrenergic receptor antagonist of PROP. For all rescue tests (excluding oxygen consumption assay), the highest concentration of RAC, 8 ppm, was selected and co-incubated with 1 or 4 ppm of PROP. The same protocols were applied as before with 24-h co-incubation on 48 hpf for cardiac and blood flow assays and 96 hpf zebrafish larvae for locomotion and oxygen consumption analysis. The test was conducted in duplicate.

3. Results

3.1. Total Distance Traveled for Zebrafish Larvae after RAC Exposure

Initially, we performed photomotor response (PMR) assay to evaluate the potential effects of RAC on zebrafish locomotor activity responding to four alternative dark/light cycle switches. After PMR assay, we extracted data and compared the total distance traveled, burst, and rotation counts of zebrafish larvae between the control and RAC exposed groups. The total distance chronology was measured on the basis of the velocity range of our previous publication [29,31]. Basically, zebrafish embryos display higher mobility in the dark phase than light phase in the PMR assay, and this pattern was still well distributed even after RAC treatment (Appendix Figure A1). We observed RAC provoked significant increase (p < 0.05) of average total distance traveled at 0.1 ppm (7.367 ± 0.124 cm), 1 ppm (7.766 ± 0.161 cm), 2 ppm (7.370 ± 0.128 cm), and 8 ppm (8.477 ± 0.182 cm) concentration in the light cycle compared to the control (6.253 ± 0.096 cm) (Figure 1A and Video S1). Also, in the dark cycle, zebrafish exposed to RAC 0.1 ppm (11.090 ± 0.135 cm), 1 ppm (11.660 ± 0.131 cm), 2 ppm (12.460 ± 0.118 cm), and 8 ppm (13.120 ± 0.185 cm) displayed higher average distance traveled than the control (8.880 ± 0.099 cm) (Figure 1B). One data point at 4 ppm, however, showing significant lower total distance traveled (p < 0.05) in both light (5.468 ± 0.163 cm) and dark cycle (8.668 ± 0.177 cm). This result reveals that the administration of RAC to zebrafish larvae can trigger psychoactive effects at both dark and light phases for the PMR assay.

3.2. Total Rotation Movement for Zebrafish Larvae after RAC Exposure

The rotational chronology (clockwise or anti-clockwise) was influenced by many discrete agents such as feeding, brightness, or external pressure [44]. In response to RAC exposure, similar with previous locomotor results, zebrafish larvae displayed increment of rotational movements at both light and dark cycles. Results of RAC exposure showed significant increased (p < 0.05) on zebrafish rotational movement activity during both light and dark cycles in the range from 1 to 8 ppm concentration (Figure 1C,D). However, the RAC 0.1 ppm group showed hypoactivity in the dark cycle rotation test. In conclusion, the administration of RAC to zebrafish embryos can trigger high rotation movement at both dark and light phases in the PMR assay.

3.3. Total Burst Count for Zebrafish Larvae after RAC Exposure

Total burst count was determined by calculating the rapid change in the body activity of an experimental animal with higher than 20 pixels per second, that is an effective and sensitive index in evaluating the irregular movement and anxiety in zebrafish [29]. Statistical analysis showed that RAC induced higher average burst movement in zebrafish larvae in light cycle at all treated concentrations i.e., 0.1 ppm (2.822 ± 0.129 burst count/min), 1 ppm (6.960 ± 0.262 burst count/min), 2 ppm (8.398 ± 0.286 burst count/min), 4 ppm (7.659 ± 0.262 burst count/min), and 8 ppm (5.732 ± 0.163 burst count/min) compared to the control group (1.946 ± 0.130 burst count/min) with p < 0.0001 (Figure 1E). In comparison to light phase, zebrafish larvae displayed more burst movement in dark phase at 0.1 ppm (19.840 ± 0.274) 1 ppm (31.270 ± 0.424 burst count/min), 2 ppm (36.610 ± 0.451 burst count/min), 4 ppm concentrations (24.750 ± 0.333 burst count/min), and 8 ppm (21.460 ± 0.298) comparing to the control (15.400 ± 0.368 burst count/min) (Figure 1F). In conclusion, the administration of RAC to zebrafish embryos can trigger high burst movement at both dark and light phases for the PMR assay.

3.4. Respiratory Rate Analysis in Response to RAC Exposure in Zebrafish Larvae

As previous locomotor data showed that RAC exposure to zebrafish larvae cause hyperactivity. This intriguing observation leads us to ask whether the locomotor hyperactivity is associated with oxygen consumption alteration. In order to validate this hypothesis, we conducted oxygen consumption measurement on zebrafish embryos after RAC treatment by using Loligo Microplate Respirometry System (Loligo Systems, Viborg, Denmark) (Figure 2A). When the blue LED irradiated on the oxygen sensitive sensor embedded in the microplate, the blue light energy will be absorbed by O2 and emitted as different intensity of red-light signals. By monitoring the relative oxygen level in the microplate over time, we found the dissolved oxygen level is sharply declined in embryos treated with RAC with significant p < 0.001 (Figure 2B). By statistical quantification, we identified a significant increase of oxygen consumption in zebrafish embryos after exposed to different concentrations of RAC starting at the lowest 0.1 ppm (3.081 ± 0.158 ppm), 4 ppm (3.258 ± 0.121 ppm), and to the highest 8 ppm (3.411 ± 0.144 ppm). This result supported our hypothesis that acute RAC exposure in zebrafish caused hyperactivity and associated with high oxygen consumption (Figure 2C).

3.5. Cardiovascular Performance Assay

Next, we performed several important assays to elucidate the corresponding alteration of cardiovascular system after RAC exposure in zebrafish embryos at 72 hpf after 24-h exposure to RAC from 48 hpf onwards. To measure cardiac performance, several cardiac physiology endpoints were analyzed, such as heart rate in atrium, heart rate in ventricle, stroke volume, cardiac output, ejection fraction, and shortening fraction. We also checked the cardiac rhythm by analyzing the heart rate variability in both atrium and ventricle chambers and time interval between atrium to ventricle relaxation and vice versa. Furthermore, blood flow velocity in dorsal aorta was analyzed to measure the effect of RAC exposure to vascular system in zebrafish.
Heart rates show the overall condition of the heart. Lower heart rate usually showed that the person has a weaker health condition. However, in some special case like athlete, they will have lower resting heart rate which come from the increase power of the heart muscle [45]. In contrary to the lower heart rate, a higher heart rate is not preferable as it in many cases related to stress conditions and other cardiovascular related diseases [46]. It is intriguing to find that although locomotor hyperactivity and high oxygen consumption were detected, heart rate still maintains a regular level after RAC exposure at RAC doses of 0.1 and 4 ppm (Figure 3A,B). At a high RAC dose of 8 ppm, significant heart rate elevation was detected in both atrium and ventricle chambers (Figure 3A,B).
Stroke volume showed the volume of blood that leaves the heart chamber for each contraction. Stroke volume was calculated by subtracting the heart volume at the end-diastolic phase by the end-systolic phase that assumed that the heart chamber is an ellipsoid shape [47]. After incubation in RAC, we observed significant elevation of stroke volume after incubation in 4 ppm (93.260 ± 8.070 pL/beat) and 8 ppm (100.800 ± 6.529 pL/beat) (Figure 3C), which suggested that RAC could increase the heart contractility.
Cardiac output is defined as the total volume of blood pumped from the ventricle per minute, which is calculated by multiplying the stroke volume with the heart rate. Just as with stroke volume result, significant elevation also observed in zebrafish cardiac output after 4 ppm (14805 ± 1437 pL/min) and 8 ppm (16618 ± 942 pL/min) incubation of RAC (Figure 3D). Those effects suggest the cardiac output elevation after RAC exposure is largely contributed by stroke volume rather than heart rate elevation.
The ejection fraction is the measurement that how much blood is pumped out of the ventricle and is calculated by dividing the stroke volume with heart volume at the end diastolic phase. Ejection fraction refers the ability of heart to pump blood [48]. Therefore, the more blood volume pumped out, the increased ability of heart when treated with 8 ppm RAC (45.94 ± 1.64 EF%) was displayed (Figure 3E). Shortening fraction is measured by calculating the changing percentage in the dimension of ventricle during systolic phase [49]. Shortening fraction is correlated with ejection fraction [50]. To pump a large volume of blood (stroke volume) with high force (Ejection fraction), the lost dimension and percentage changes in ventricular dimension occurred when treated with 4 ppm (16.86 ± 1.41 SF%) and 8 ppm RAC (18.42 ± 0.96 SF%) (Figure 3F).
Finally, we also measured heart rate variability to evaluate whether acute RAC exposure can trigger arrhythmia in zebrafish embryos or not. Furthermore, we also checked the time interval between atrium to ventricle muscle relaxation and vice versa. By using Poincare plot, sd1 and sd2 value showed no significant difference in both the atrium and ventricle chamber (Figure 4A–D). No significant change was observed in the time interval between each chamber beat in response to RAC exposure, which suggested that RAC exposure did not change the heart rate variability and rhythm of zebrafish larvae (Figure 4E,F).

3.6. Blood Flow Measurement after RAC Exposure

In this endpoint, we observed the blood flow velocity in the dorsal aorta of zebrafish in correlation with cardiac analysis (The measurement position was highlighted in Figure 5A,B). The blood flow images were captured by high-speed CCD setting and later the maximal and average blood flow velocity was analyzed by using ImageJ-based method according to our previously published protocol described by Santoso et al. [33]. Heart rate and blood flow are tightly associated each other and could be influenced by chemical exposure like RAC. Results show RAC at 4 ppm can significantly enhance the average and maximal blood flow velocity with statistical means ± SEM value of 448.70 ± 24.18 µm/s and 963.40 ± 65.73 µm/s, respectively (Figure 5C–E and Video S2). Based on data collected from Figure 3 and Figure 4, we found RAC exposure affects neither heart rate nor cardiac rhythm, but indeed can influence the stroke volume, cardiac output and ejection, and shortening fractions in zebrafish embryos. This alteration triggered by RAC can enhance the pumping capacity and blood volume without effecting size of heart chamber and cardiac rhythm in zebrafish embryos. Relating our blood flow results with cardiac performance in zebrafish, we concluded that RAC can increase the blood flow velocity in zebrafish without increasing heart rate too much (comparing Figure 3A,B and Figure 5D,E).

3.7. Molecular Docking for RAC and Endogenous β-adrenergic Receptor

In zebrafish genome, 10 genes were annotated as β-adrenergic receptors. We tried to validate which receptor has the higher affinity with RAC and might play a role on mediating RAC physiological activity, also comparing with PROP (a β-adrenergic receptor blocker). First, we searched ZFIN database (https://zfin.org/; accessed on 14 June 2021) for gene annotated with adrenergic receptors and later download the amino acid sequences for structure fitting. Briefly, after homology modeling structures of zebrafish endogenous β-adrenergic receptor subtypes were made. The amino acid sequence (aa) of zebrafish (Danio rerio) ten subtypes were selected from the Uniprot database (https://www.uniprot.org/; accessed on 9 March 2021) and their corresponding IDs and chromosomal position were listed in Table 1. The template structures for modeling structures were picked by NCBI blastp with the four PDB crystal protein structures (with 29.66% to 69.15% protein sequence identities). Ultimately, the crystal structure, 6kux, is chosen because the 6kux ligand has the best overlap with RAC among the four ones. The MODELLER software v9.20 constructed all homology modeling structures for flexible molecular docking. CDOCKER scores were assigned by the molecular docking module, flexible CDOCKER, indicating the binding affinity between each target protein and RAC [36,37]. Results from molecular docking showed that RAC may be able to interact with ten subtypes of zebrafish adrenergic receptors similar with the human homolog with CDOCKER scores of either 47.62 (adra1aa), 44.77 (adra1ab), 42.85 (adra1ba), 44.83 (adra1bb), 42.50 (adra1d), 43.68 (adra2a), 46.91 (adra2b), 41.90 (adra2c), 44.27 (adra2da), 42.74 (adra2db), or 44.15 (adra2a, Homo sapiens). PROP interacts with ten subtypes of zebrafish adrenergic receptors with CDOCKER scores of either 25.37, 25.97 (adra1aa), 22.98, 22.40 (adra1ab), 22.51, 22.41 (adra1ba), 25.58, 27.33 (adra1bb), 23.17, 23.37 (adra1d), 21.53, 23.25 (adra2a), 24.93, 24.48 (adra2b), 24.20, 23.57 (adra2c), 24.42, 22.44 (adra2da), 25.17, 22.50 (adra2db), or 21.83, 20.62 (adra2a, Homo sapiens) in R form and S form, respectively. Based on in silico data, PROP overall has about 50% lower binding affinity than RAC. Take zebrafish adra1aa as an example, the endogenous β-adrenergic receptor binding pocket for RAC was identified in the middle of the hollow cylinder of the seven transmembrane domains by homology modeling (Figure 6A). Three-dimensional (3D) (Figure 6B,C) structures illustrate three hydrogen bond formations between RAC and zebrafish adra1aa at the positions Glu101, Ser418, and Asp120 (highlighted by the green and blue dotted line). In conclusion, based on molecular docking, RAC might be able to bind to the ten zebrafish β-adrenergic receptors subtypes, from adra1aa to adra2db, with strong binding affinity supported compared to their human homolog.

3.8. Test Physiological Effects of β-Blocker ‘Propranolol’

We also tested potential effects of PROP exposure on zebrafish larvae using different behavioral and physiological endpoints after 24-h incubation. PROP has been found to be effective in performance and somatic anxiety [51]. Our PROP data supports the previous study showing PROP failed to affect the anxiolytic activity in a specific manner and pattern [52]. At the highest tested dose of PROP (8 ppm), zebrafish larvae showed lowered total distance in both light/dark cycle and this pattern was the same in other endpoints like rotational and burst movement (data shown in Table A1). But other tested concentration exhibited an irregular pattern in locomotory behavior. This shows that PROP can be only effective at high dose to treat the acute anxiety disorders. Previous studies reported the PROP toxicity after chronic exposure of 14 days in Japanese Medaka resulted in reduced growth rate [53] and significant changes in plasma steroid level [54]. In vertebrates, β-adrenoreceptors are mainly located at cardiac and skeletal muscles [55] and are also found in fish myocardium [56].
In our study, PROP exposure can significantly lower the heart rate at 4 and 8 ppm concentrations (for both atrium and ventricle chambers, Table A2) support previous literature show PROP lowered heart rate in zebrafish larvae [53,57]. Interestingly, we did not observe any significant alteration in average and maximum blood flow velocity in zebrafish larvae exposed to different doses of PROP (Table A2) despite lowered atrium and ventricle heart rate was detected. In addition, oxygen utilization was significantly higher at PROP 8 ppm dose (Table A2) as that of RAC. As our molecular docking data demonstrates that PROP has less binding affinity than RAC (Figure 6 and Table 1). Previous literature states that PROP has been found to treat cardiovascular diseases and lowering the heart rate [8,58] and β-adrenoreceptors are primarily located in cardiac tissues [55]. Basically, our cardiac results from PROP exposed zebrafish larvae support above previous studies. However, PROP failed to exhibit and set any specific pattern for other behavioral and physiological parameters we studied. Further acute and chronic studies of PROP exposure could help to better understand the effect on behavioral and physiological parameters.

3.9. Rescue Physiological Alterations Induced by RAC with PROP

β-blockers are extensively used for the treatment of cardiovascular diseases in human like hypertension, heart failure, arrhythmia, angina, migraines, and other related conditions [8,58,59,60], which works by competitive inhibition of β-adrenergic receptors. One of these drugs is PROP, which is being extensively used since its discovery in 1960s [61]. Considering the previous literature [62], we tested two PROP doses (1 and 4 ppm) at our RAC highest tested dose (8 ppm) and analyzed different behavioral and physiological endpoints like locomotion, oxygen consumption, cardiac physiology, and blood flow.
As a result, co-incubation of RAC and PROP indeed can significantly reduce locomotion hyperactivity triggered by 8 ppm RAC. The total distance travelled in both light and dark cycles for RAC exposed fish larvae can be significantly reduced after receiving either 1 or 4 ppm PROP administration (Figure 7A,B). For rotation behavior, we found PROP exposure also can effectively rescue high rotation behavior triggered by RAC in both light and dark cycles (Figure 7C,D). For burst movement behavior, the relative high burst activity can effectively restore after PROP exposure in both light and dark cycles (Figure 7E,F). Taking together, all evidences collected from light–dark switching photomotor assay suggests RAC and PROP function by competing with β-adrenoreceptor in the opposite direction.
As per previous literature, β-blockers are primarily used for cardiovascular diseases on fighting against arrhythmia [8,58,60]. Co-incubation of PROP with RAC significantly reduced the tachycardia as per our hypothesis at highest tested concentration. Figure 8A,B shows tachycardia phenotype triggered by 8 ppm RAC exposure can be effectively rescued by either 1 or 4 ppm PROP administration. Other cardiac physiological endpoints like stroke volume, cardiac output, ejection and shortening fraction display no significant alterations after PROP co-incubation (Figure 8C–F). Like our previous results show in Figure 4, RAC exposure did not alter the cardiac rhythm endpoints i.e., A–V and V–A interval, SD 1 and SD 2 for either atrium or ventricle. In contrary, co-incubation of both RAC and PROP significantly altered the cardiac rhythm endpoints by showing more higher levels of SD1 and SD2 in atrium (Figure 8I,J). The SD1 and SD2 levels in ventricle, on the contrary, display no significant different on RAC–PROP exposure compared to control (Figure 8K,L).
Since cardiac function alterations are tightly associated with blood flow dynamics, we aimed to analyze possible blood flow alteration after PROP rescue in RAC exposed fish. We used the same protocol as mentioned in Section 3.6 for blood flow measurement. RAC significantly increased the maximum and average blood flow velocity at a 4 ppm dose as shown in Figure 5D,E. Co-incubating PROP with RAC, however, reduced the maximum and average blood flow velocity as regarding with only RAC exposure. Significant differences in maximum and average blood flow velocity after co-treated RAC with PROP 1 and 4 ppm were displayed in Figure 9A,B when compared to RAC only. As per our results from oxygen consumption after RAC exposure in zebrafish larvae, all of the tested concentrations consumed more oxygen than the control (Figure 2C). PROP exposure alone did not seem to be lowering the oxygen consumption rate in zebrafish larvae (Table A2). However, the RAC–PROP co-incubation groups displayed significant lower oxygen consumption level (p < 0.0001) compared to the RAC-treated group (Figure 9C). Further, co-incubation with PROP 1 ppm could lowered the oxygen uptake to similar level with control (p = 0.2658). Taking together, all these data support RAC and PROP act through β-adrenergic receptors on mediating blood flow control in the opposite direction.

4. Discussion

Previously, it was first reported that RAC was associated with arterial and cardiac damage in canine. Based on cardiac anatomy, apoptosis within myocyte was the key reason that showed heart failure due to RAC exposure [20]. However, there was lack of extensive investigation of mechanism when exposing RAC. Therefore, in this study, we demonstrated serial experiments that associated with locomotion, oxygen consumption, and cardiovascular physiology endpoints, which were dedicated to identify the possible physiological alterations by using zebrafish as a simple animal model. The most important finding for this paper is that we identified the RAC function as a highly potent chemical reagent to boost the cardiovascular, oxygen consumption and locomotor activities in zebrafish embryos. This finding basically fit with the typical function of RAC as metabolism modifier to promote the metabolic rate in livestock [9,63] and fish [23]. Further, we performed a rescue trial in response to RAC induced alterations, using PROP as β-blocker or β-adrenergic receptor antagonist, which is known as non-selective β-blocker agent and is chiral in nature, to validate our results. We also analyzed the effects of only PROP exposure on zebrafish larvae using locomotory, cardiac, oxygen consumption, and blood flow endpoints on different concentrations. By utilizing this combinational approach, we were able to validate the potential mechanism of RAC on mediating zebrafish physiological alterations.
Considering the behavioral effects of acute RAC exposure to zebrafish larvae, we found the hyperactivity (increased total distance traveled, total rotational movements, and total burst count) in zebrafish larvae, which is similar to some previous studies in different animal models including zebrafish larvae [64]. In pigs, RAC treated animals were found to be more active and alert just after the onset of RAC feeding [63], but after the third week there was no difference in pig behavior. RAC was also found to be interacting with stress in rats [65]. Different studies of zebrafish exposure to RAC have stated increased locomotor activity and exploratory behavior, showing hyperactivity in zebrafish larvae [23]. Another study demonstrates the unaltered behavior in adult zebrafish at highest tested dose [22]. Taken these studies together, it can be said that RAC affects locomotory behavior differently depending upon developmental stages of zebrafish. While co-incubation of RAC with PROP, significantly increased the anxiolytic behavior of different locomotory endpoints in zebrafish larvae at the highest tested dose, which supports the previous studies that PROP has a significant role in reducing anxiety and hypertension [66,67,68]. Compared to RAC 8 ppm, co-incubation of PROP with RAC significantly reduced different locomotory endpoints, especially for burst and rotation.
RAC exposure significantly affected the cardiac physiology endpoints, i.e., heart rate, stroke volume, cardiac output, ejection fraction, and shortening fraction without significantly altering the cardiac rhythm endpoints. Significant results of cardiac physiology endpoints relate the findings with increased blood flow in zebrafish larvae exposed to RAC. It is intriguing to notice the faster heart rate and blood flow in RAC exposed zebrafish embryos did not trigger any sign of heart rate irregularity. This result is consistent with previous findings reported in finishing steers [23,69]. While co-incubation of PROP with RAC, significantly reduced the heart rate in the atrium and ventricle, which confirms the previous studies and supports the working mechanism of PROP in reducing heart rate [53,57,62,70]. The co-incubation of RAC and PROP also can rescue the high blood velocity and high oxygen consumption rate triggered by RAC even PROP are given at low dose as 1 ppm.
The alteration effects showed in this study are mainly administered by RAC that served as β-adrenergic agonist. The β-adrenergic agonist has been previously reported to regulate multiple physiological and metabolic functions including muscle system, heart rate, oxygen uptake, and also hemoglobin–oxygen affinity in fish [31,53,71]. In this study, we provided evidence of RAC physiological effects using well-established animal model, zebrafish by using multiple parameters. In line with previous study, RAC is able to increase the distance traveled and heart rate in zebrafish embryo and larvae [23]. Further, the increased oxygen consumption in RAC treatment might be due to additional muscle mass accretion. The over-abundance of F-actin-capping protein, which involves in skeletal muscle organization steers the oxygen transportation acceleration [72]. It also has been reported increasing anxiety, aggressive behavior, blood pressure, and palpitation in rats and pigs fed RAC [63,65].
For mechanism study, we have done the molecular docking analysis to demonstrate RAC have generally strong binding affinity to ten zebrafish subtypes of endogenous β-adrenergic receptors comparable to a human subtype. Therefore, based on molecular docking, we proposed the high potency on elevating cardiovascular and locomotion physiology in zebrafish might be mediated by multiple β-adrenergic receptor activation. It has been reported that S-propranolol has more β-adrenergic blocking activity than R-propranolol in mammals [73]. However, our in silico molecular docking study shows that the zebrafish subtype, adra1aa, binding to PROP with a chiral center, gains the highest CDOCKER scores for S form, 27.33, while the human adra2a (6kux) gains 21.83 and 20.62, in R form and S form, respectively. PROP overall has lower binding affinity than RAC (Figure 6 and Table 1). According to zebrafish genetic database adra1ba, adra1bb, adra2a, adra2b, and adra2db are expressed in heart of zebrafish while adra2a, adra2da, and adra2c are expressed in CNS [74]. adra1d and adra1ab genes are responsible for blood circulation while adra2da is expressed in the musculature system [75]. Our molecular docking results show that adra1aa has the highest binding affinity (47.63) among all tested subtypes (Table 1), which in human is responsible for Alzheimer’s disease and hypertension [76].
Based on this data, we were intrigued to evaluate the competition effect of RAC and PROP. We believed despite PROP having a lower binding affinity than RAC, it might still be able to display an antagonistic effect to inhibit RAC. In accordance with a prior study, PROP as β-receptor antagonist responsible for the neutralization of RAC effects [77]. Co-incubation of PROP and RAC in the same doses prevented RAC-stimulated cAMP production in myotubes, which is highly related in skeletal muscle adaptation. Another study showed PROP co-administered with β-adrenergic agonist, clenbuterol, is able to block the rise in cardiac muscle growth response, but PROP itself did not inhibit the ability of clenbuterol to stimulate protein accretion [78]. As of yet, the exact explanation remains unclear. In this case, whether PROP is literally able to inhibit RAC directly or it acts by working through a different β-receptors may need further investigations. In the future, administration of β-adrenergic receptor subtype-specific antagonist or perform systematically gene knockout/knockdown experiment for zebrafish endogenous β-adrenergic receptors might be able to provide more direct evidence to elucidate this hypothesis. Literature also shows that β-adrenergic agonists are associated with increased cardiac rate in different animal species like dogs along with pigs and fish [20,79]. The persistent increased cardiac rate in response to β-adrenoceptor agonist can results in cardiomyocyte death due to inflammation and impaired calcium handling [80]. Further studies are required to explore the potential toxicity of chronic RAC exposure of zebrafish on cardiovascular physiology.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cells10092449/s1, Video S1: Photomotor assay for zebrafish larvae exposed to ractopamine at 8 ppm (deposited at YouTube: https://www.youtube.com/watch?v=irfQ4ykIOpE&t=2s; accessed on 29 August 2021), Video S2: Blood flow tracking for zebrafish larvae exposed to ractopamine at 4 ppm (deposited at YouTube: https://www.youtube.com/watch?v=TDNzbCEGyTA; accessed on 29 August 2021).

Author Contributions

Conceptualization, K.A. and C.-D.H.; Data curation, K.A., F.S., M.E.S., Y.-H.L., J.-C.H. and W.-H.Y.; Funding acquisition and Investigation, K.H.-C.C. and C.-D.H.; Project administration, C.-D.H.; Resources, Y.-T.L.; Supervision, C.-D.H.; Validation and Visualization, Y.-T.L. and F.S.; Writing—original draft, K.A., K.H.-C.C., Y.-T.L. and C.-D.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the grants sponsored by the Ministry of Science and Technology MOST 108-2313-B-033 -001-MY3 to C.-D.H. and grant number MOST 108-2113-M153-003 to K.H.-C.C.

Institutional Review Board Statement

All experiments involving zebrafish were performed following the guidelines approved by the Institutional Animal Care and Use Committees (IACUCs) of the Chung Yuan Christian University (Approval No. 109001, issue date 15 January 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available directed to the corresponding authors.

Acknowledgments

We appreciate Taiwan Zebrafish Stock Center at Academia Sinica (TZCAS) for providing zebrafish AB strain.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Cells 10 02449 g0a1 550
Figure A1. Comparison of locomotion behavioral pattern of zebrafish exposed to ractopamine at 0, 0.1, 1, 2, 4, and 8 ppm during both light and dark cycles. (A) Average total distance traveled, (B) average rotation, and (C) average burst movement count from 0 to 80 min were analyzed. The data are presented as mean with SEM (n = 96 for all groups).
Figure A1. Comparison of locomotion behavioral pattern of zebrafish exposed to ractopamine at 0, 0.1, 1, 2, 4, and 8 ppm during both light and dark cycles. (A) Average total distance traveled, (B) average rotation, and (C) average burst movement count from 0 to 80 min were analyzed. The data are presented as mean with SEM (n = 96 for all groups).
Cells 10 02449 g0a1
Table
Table A1. Comparison of zebrafish larvae locomotor activity endpoints after treated with 5 different doses of Propranolol.
Table A1. Comparison of zebrafish larvae locomotor activity endpoints after treated with 5 different doses of Propranolol.
Behavior EndpointsPropranolol Concentrations
0 ppm (control)0.1 ppm1 ppm2 ppm4 ppm8 ppm
Mean ± SEM Mean ± SEM Mean ± SEM Mean ± SEM Mean ± SEM Mean ± SEM
Average Total
Distance (cm)		Light cycle7.35 ± 0.126a6.35 ± 0.128b7.79 ± 0.171c7.42 ± 0.293e10.09 ± 0.324d6.59 ± 0.163d
Dark cycle10.06 ± 0.106a11.05 ± 0.107b11.68 ± 0.162b9.88 ± 0.286c12.43 ± 0.336d7.49 ± 0.185d
Average Rotation (count)Light cycle0.032 ± 0.002a0.072 ± 0.005b0.041 ± 0.003c0.029 ± 0.003d0.023 ± 0.002e0.009 ± 0.002e
Dark cycle0.214 ± 0.006a0.422 ± 0.013b0.253 ± 0.008c0.147 ± 0.007d0.164 ± 0.008e0.034 ± 0.003e
Average Burst Movement (count)Light cycle0.654 ± 0.037a1.598 ± 0.091b1.485 ± 0.078c1.049 ± 0.064b1.124 ± 0.062a0.469 ± 0.025a
Dark cycle7.25 ± 0.125a29.74 ± 0.441b21.41 ± 0.320c11.77 ± 0.275a9.847 ± 0.257e1.797 ± 0.069e
Significant differences (p < 0.05) of each locomotor activity between different concentrations of propranolol were tested by Kruskal–Wallis test followed by Dunn’s test and the significances were indicated by small letters a, b, c, d, e (in rows). Small letters with different label indicate they are statistically different.
Table
Table A2. Comparison of cardiac performance, blood flow velocity, and oxygen consumption in larval zebrafish after treated with 5 different doses of Propranolol.
Table A2. Comparison of cardiac performance, blood flow velocity, and oxygen consumption in larval zebrafish after treated with 5 different doses of Propranolol.
Cardiac Related ParametersPropranolol Concentrations
0 ppm (control) 0.1 ppm 4 ppm 8 ppm
Mean ± SEM Mean ± SEM Mean ± SEM Mean ± SEM
Heart Rate Ventricle (bpm)145.3 ± 2.956a140.6 ± 3.316ab134.7 ± 2.701bc130.8 ± 1.975c
Heart Rate Atrium (bpm)145.4 ± 2.888a140.4 ± 3.427ab134.9 ± 2.664bc130.7 ± 2.039c
Stroke Volume (pL/beat)80.94 ± 7.274a78.68 ± 7.275a71.37 ± 7.992a75.75 ± 4.338a
Cardiac Output (pL/min)11663 ± 101a10978 ± 983.0a9424 ± 993.9a9926 ± 568.2a
Ejection Fraction (%)37.36 ± 2.101a35.85 ± 2.148a32.37 ± 2.826a34.57 ± 1.508a
Shortening Fraction (%)14.68 ± 1.280a13.84 ± 1.078a12.36 ± 1.582a12.56 ± 0.728a
AV–Interval (s)0.1864 ± 0.022a0.161 ± 0.0206a0.149 ± 0.015a0.1785 ± 0.021a
VA–Interval (s)0.228 ± 0.0196a0.272 ± 0.023a0.299 ± 0.020b0.2821 ± 0.021a
SD1 Atrium (s)0.036 ± 0.006a0.034 ± 0.007a0.042 ± 0.006a0.041 ± 0.007a
SD2 Atrium (s)0.0212 ± 0.003a0.0214 ± 0.004ac0.027 ± 0.004a0.034 ± 0.005c
SD1 Ventricle (s)0.0219 ± 0.002a0.0319 ± 0.005ab0.033 ± 0.003b0.0274 ± 0.003ab
SD2 Ventricle (s)0.0136 ± 0.002a0.0204 ± 0.003b0.021 ± 0.002b0.017 ± 0.002ab
Max. Blood Flow Velocity (µm/s)856.2 ± 44.20a763.0 ± 71.85a769.1 ± 57.77a876.2 ± 51.92a
Avg. Blood Flow Velocity (µm/s)405.3 ± 15.55a348.5 ± 27.02a360.6 ± 24.48a400.6 ± 20.90a
Oxygen Consumption Rate (ppm)2.081 ± 0.182a2.485 ± 0.209a2.306 ± 0.171a2.685 ± 0.176b
Significant differences (p < 0.05) of each locomotor activity between different concentrations of propranolol were indicated by small letters a, b, and c (in rows) and analyzed by one-way ANOVA with Fisher’s LDS test or with Kruskal–Wallis followed by Dunn’s test.

References

  1. Mennigen, J.A.; Stroud, P.; Zamora, J.M.; Moon, T.W.; Trudeau, V.L. Pharmaceuticals as neuroendocrine disruptors: Lessons learned from fish on Prozac. J. Toxicol. Environ. Health Part B 2011, 14, 387–412. [Google Scholar] [CrossRef]
  2. Houtman, C.J. Emerging contaminants in surface waters and their relevance for the production of drinking water in Europe. J. Integr. Environ. Sci. 2010, 7, 271–295. [Google Scholar] [CrossRef]
  3. Brumovský, M.; Bečanová, J.; Kohoutek, J.; Borghini, M.; Nizzetto, L. Contaminants of emerging concern in the open sea waters of the Western Mediterranean. Environ. Pollut. 2017, 229, 976–983. [Google Scholar] [CrossRef] [PubMed]
  4. Hakk, H.; Shelver, W.L.; Casey, F.X. Fate and transport of the β-adrenergic agonist ractopamine hydrochloride in soil–water systems. J. Environ. Sci. 2016, 45, 40–48. [Google Scholar] [CrossRef] [PubMed]
  5. Kriewald, R.D. Effects of Ractopamine HCL on Physical and Reproductive Parameters in the Horse. PhD diss., A&M University, Houston, TX, USA, 2010. [Google Scholar]
  6. Ali, D.C.; Naveed, M.; Gordon, A.; Majeed, F.; Saeed, M.; Ogbuke, M.I.; Atif, M.; Zubair, H.M.; Changxing, L. β-Adrenergic receptor, an essential target in cardiovascular diseases. Hear. Fail. Rev. 2019, 25, 343–354. [Google Scholar] [CrossRef]
  7. Fent, K.; Weston, A.; Caminada, D. Ecotoxicology of human pharmaceuticals. Aquat. Toxicol. 2006, 76, 122–159. [Google Scholar] [CrossRef] [PubMed]
  8. Massarsky, A.; Trudeau, V.L.; Moon, T.W. β-blockers as endocrine disruptors: The potential effects of human β-blockers on aquatic organisms. J. Exp. Zool. A Ecol. Genet. Physiol. 2011, 315, 251–265. [Google Scholar] [CrossRef]
  9. de Almeida, V.V.; Nuñez, A.J.C.; Miyada, V.S. Ractopamine as a metabolic modifier feed additive for finishing pigs: A review. Braz. Arch. Biol. Technol. 2012, 55, 445–456. [Google Scholar] [CrossRef]
  10. Liu, X.; Grandy, D.K.; Janowsky, A. Ractopamine, a livestock feed additive, is a full agonist at trace amine–associated receptor 1. J. Pharmacol. Exp. Ther. 2014, 350, 124–129. [Google Scholar] [CrossRef] [PubMed]
  11. Centner, T.J.; Alvey, J.C.; Stelzleni, A.M. Beta agonists in livestock feed: Status, health concerns, and international trade. J. Anim. Sci. 2014, 92, 4234–4240. [Google Scholar] [CrossRef]
  12. Lee, H.-C.; Chen, C.-M.; Wei, J.-T.; Chiu, H.-Y. Analysis of veterinary drug residue monitoring results for commercial livestock products in Taiwan between 2011 and 2015. J. Food Drug Anal. 2018, 26, 565–571. [Google Scholar] [CrossRef]
  13. Anderson, D.B.; Moody, D.E.; Hancock, D.L. Beta adrenergic agonists. In Encycl. Anim. Sci.; Marcel-Dekker Inc.: New York, NY, USA, 2004; pp. 104–108. [Google Scholar]
  14. Dalidowicz, J.; Thomson, T.; Babbitt, G. Ractopamine hydrochloride, a phenethanolamine repartitioning agent: Metabolism and tissue residues. In Xenobiotics and Food-Producing Animals; ACS Publications: Washington, DC, USA, 1992; pp. 234–243. [Google Scholar]
  15. Bartelt-Hunt, S.; Snow, D.; Damon-Powell, T.; Miesbach, D. Occurrence of steroid hormones and antibiotics in shallow groundwater impacted by livestock waste control facilities. J. Contam. Hydrol. 2011, 123, 94–103. [Google Scholar] [CrossRef] [PubMed]
  16. Jaimes-Correa, J.C.; Snow, D.D.; Bartelt-Hunt, S.L. Seasonal occurrence of antibiotics and a beta agonist in an agriculturally-intensive watershed. Environ. Pollut. 2015, 205, 87–96. [Google Scholar] [CrossRef] [PubMed]
  17. Sakai, N.; Sakai, M.; Haron, D.E.M.; Yoneda, M.; Mohd, M.A. Beta-agonist residues in cattle, chicken and swine livers at the wet market and the environmental impacts of wastewater from livestock farms in Selangor State, Malaysia. Chemosphere 2016, 165, 183–190. [Google Scholar] [CrossRef] [PubMed]
  18. Yu, T.-H.; Lin, A.; Wang, X.-H.; Lin, C.-F. Occurrence of β-blockers and β-agonists in hospital effluents and their receiving rivers in southern Taiwan. DESALINATION Water Treat. 2011, 32, 49–56. [Google Scholar] [CrossRef]
  19. Catalano, D.; Odore, R.; Amedeo, S.; Bellino, C.; Biasibetti, E.; Miniscalco, B.; Perona, G.; Pollicino, P.; Savarino, P.; Tomassone, L.; et al. Physiopathological changes related to the use of ractopamine in swine: Clinical and pathological investigations. Livest. Sci. 2012, 144, 74–81. [Google Scholar] [CrossRef]
  20. Yaeger, M.; Mullin, K.; Ensley, S.; Ware, W.; Slavin, R. Myocardial toxicity in a group of greyhounds administered ractopamine. Vet. Pathol. 2012, 49, 569–573. [Google Scholar] [CrossRef]
  21. Sun, L.; Wang, S.; Lin, X.; Tan, H.; Fu, Z. Early Life Exposure to Ractopamine Causes Endocrine-Disrupting Effects in Japanese Medaka (Oryzias latipes). Bull. Environ. Contam. Toxicol. 2015, 96, 150–155. [Google Scholar] [CrossRef]
  22. Sachett, A.; Bevilaqua, F.; Chitolina, R.; Garbinato, C.; Gasparetto, H.; Dal Magro, J.; Conterato, G.M.; Siebel, A.M. Ractopamine hydrochloride induces behavioral alterations and oxidative status imbalance in zebrafish. J. Toxicol. Environ. Health Part A 2018, 81, 194–201. [Google Scholar] [CrossRef]
  23. Garbinato, C.; Schneider, S.E.; Sachett, A.; Decui, L.; Conterato, G.M.; Müller, L.G.; Siebel, A.M. Exposure to ractopamine hydrochloride induces changes in heart rate and behavior in zebrafish embryos and larvae. Environ. Sci. Pollut. Res. 2020, 27, 21468–21475. [Google Scholar] [CrossRef]
  24. Ferreira, L.G.; Dos Santos, R.N.; Oliva, G.; Andricopulo, A.D. Molecular Docking and Structure-Based Drug Design Strategies. Molecules 2015, 20, 13384–13421. [Google Scholar] [CrossRef]
  25. Li, J.; Fu, A.; Zhang, L. An Overview of Scoring Functions Used for Protein–Ligand Interactions in Molecular Docking. Interdiscip. Sci. Comput. Life Sci. 2019, 11, 320–328. [Google Scholar] [CrossRef] [PubMed]
  26. Lohning, A.E.; Levonis, S.M.; Williams-Noonan, B.; Schweiker, S. A Practical Guide to Molecular Docking and Homology Modelling for Medicinal Chemists. Curr. Top. Med. Chem. 2017, 17, 2023–2040. [Google Scholar] [CrossRef]
  27. He, J.-H.; Gao, J.-M.; Huang, C.-J.; Li, C.-Q. Zebrafish models for assessing developmental and reproductive toxicity. Neurotoxicology Teratol. 2014, 42, 35–42. [Google Scholar] [CrossRef] [PubMed]
  28. Avdesh, A.; Chen, M.; Martin-Iverson, M.; Mondal, A.; Ong, D.; Rainey-Smith, S.; Taddei, K.; Lardelli, M.; Groth, D.M.; Verdile, G.; et al. Regular Care and Maintenance of a Zebrafish (Danio rerio) Laboratory: An Introduction. J. Vis. Exp. 2012, e4196. [Google Scholar] [CrossRef]
  29. Hussain, A.; Audira, G.; Malhotra, N.; Uapipatanakul, B.; Chen, J.-R.; Lai, Y.-H.; Huang, J.-C.; Chen, K.H.-C.; Lai, H.-T.; Hsiao, C.-D. Multiple Screening of Pesticides Toxicity in Zebrafish and Daphnia Based on Locomotor Activity Alterations. Biomolecules 2020, 10, 1224. [Google Scholar] [CrossRef]
  30. Kristofco, L.A.; Cruz, L.C.; Haddad, S.P.; Behra, M.L.; Chambliss, C.K.; Brooks, B.W. Age matters: Developmental stage of Danio rerio larvae influences photomotor response thresholds to diazinion or diphenhydramine. Aquat. Toxicol. 2015, 170, 344–354. [Google Scholar] [CrossRef]
  31. Suryanto, M.E.; Audira, G.; Uapipatanakul, B.; Hussain, A.; Saputra, F.; Siregar, P.; Chen, K.H.-C.; Hsiao, C.-D. Antidepressant screening demonstrated non-monotonic responses to amitriptyline, amoxapine and sertraline in locomotor activity assay in larval zebrafish. Cells 2021, 10, 738. [Google Scholar] [CrossRef]
  32. Sampurna, B.P.; Audira, G.; Juniardi, S.; Lai, Y.-H.; Hsiao, C.-D. A simple image-based method to measure cardiac rhythm in zebrafish embryos. Inventions 2018, 3, 21. [Google Scholar] [CrossRef]
  33. Santoso, F.; Sampurna, B.P.; Lai, Y.-H.; Liang, S.-T.; Hao, E.; Chen, J.-R.; Hsiao, C.-D. Development of a simple image-based method for dynamic blood flow tracking in zebrafish embryos and its application in drug toxicity evaluation. Inventions 2019, 4, 65. [Google Scholar] [CrossRef]
  34. Tinevez, J.-Y.; Perry, N.; Schindelin, J.; Hoopes, G.; Reynolds, G.; Laplantine, E.; Bednarek, S.Y.; Shorte, S.; Eliceiri, K. TrackMate: An open and extensible platform for single-particle tracking. Methods 2017, 115, 80–90. [Google Scholar] [CrossRef]
  35. Rueden, C.; Schindelin, J.E.; Hiner, M.C.; Dezonia, B.E.; Walter, A.E.; Arena, E.T.; Eliceiri, K.W. ImageJ2: ImageJ for the next generation of scientific image data. BMC Bioinform. 2017, 18, 1–26. [Google Scholar] [CrossRef]
  36. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [PubMed]
  37. Webb, B.; Sali, A. Comparative protein structure modeling using modeller. Curr. Protoc. Bioinform. 2016, 54, 561–567. [Google Scholar] [CrossRef]
  38. Koska, J.; Spassov, V.Z.; Maynard, A.J.; Yan, L.; Austin, N.; Flook, P.; Venkatachalam, C.M. Fully Automated Molecular Mechanics Based Induced Fit Protein−Ligand Docking Method. J. Chem. Inf. Model. 2008, 48, 1965–1973. [Google Scholar] [CrossRef] [PubMed]
  39. Gagnon, J.K.; Law, S.M.; Brooks, C.L. Flexible CDOCKER: Development and application of a pseudo-explicit structure-based docking method within CHARMM. J. Comput. Chem. 2016, 37, 753–762. [Google Scholar] [CrossRef]
  40. Shen, M.-Y.; Sali, A. Statistical potential for assessment and prediction of protein structures. Protein Sci. 2006, 15, 2507–2524. [Google Scholar] [CrossRef]
  41. John, B.; Sali, A. Comparative protein structure modeling by iterative alignment, model building and model assessment. Nucleic Acids Res. 2003, 31, 3982–3992. [Google Scholar] [CrossRef] [PubMed]
  42. Pei, J.; Lai, L. Binding Site Detection and Druggability Prediction of Protein Targets for Structure- Based Drug Design. Curr. Pharm. Des. 2013, 19, 2326–2333. [Google Scholar] [CrossRef]
  43. Brooks, B.R.; Bruccoleri, R.E.; Olafson, B.D.; States, D.J.; Swaminathan, S.; Karplus, M. CHARMM: A program for macro-molecular energy, minimization, and dynamics calculations. J. Comput. Chem. 1983, 4, 187–217. [Google Scholar] [CrossRef]
  44. Bownik, A.; Sokołowska, N.; Ślaska, B. Effects of apomorphine, a dopamine agonist, on Daphnia magna: Imaging of swimming track density as a novel tool in the assessment of swimming activity. Sci. Total. Environ. 2018, 635, 249–258. [Google Scholar] [CrossRef]
  45. Mendes, J.; Pereira, J.; Pereira, T. Variability of Heart Rate in Athletes and Non Athletes. Eur. J. Public Heal. 2019, 29. [Google Scholar] [CrossRef]
  46. Chaswal, M.; Kapoor, R.; Batra, A.; Verma, S.; Yadav, B.S. Heart rate variability and cardiovascular reflex tests for assessment of autonomic functions in preeclampsia. Int. J. Hypertens. 2018, 2018, 8163824. [Google Scholar] [CrossRef] [PubMed]
  47. Hsiao, C.-D.; Wu, H.-H.; Malhotra, N.; Liu, Y.-C.; Wu, Y.-H.; Lin, Y.-N.; Saputra, F.; Santoso, F.; Chen, K.H.-C. Expression and purification of recombinant ghk tripeptides are able to protect against acute cardiotoxicity from exposure to water-borne-copper in zebrafish. Biomolecules 2020, 10, 1202. [Google Scholar] [CrossRef] [PubMed]
  48. Cikes, M.; Solomon, S.D. Beyond ejection fraction: An integrative approach for assessment of cardiac structure and function in heart failure. Eur. Heart J. 2015, 37, 1642–1650. [Google Scholar] [CrossRef] [PubMed]
  49. van Dalen, E.C.; van den Brug, M.; Caron, H.N.; Kremer, L.C. Anthracycline-induced cardiotoxicity: Comparison of recommendations for monitoring cardiac function during therapy in paediatric oncology trials. Eur. J. Cancer 2006, 42, 3199–3205. [Google Scholar] [CrossRef] [PubMed]
  50. Chengode, S. Left ventricular global systolic function assessment by echocardiography. Ann. Card. Anaesth. 2016, 19, 26–34. [Google Scholar] [CrossRef]
  51. Tyrer, P.J.; Lader, M.H. Response to Propranolol and Diazepam in Somatic and Psychic Anxiety. BMJ 1974, 2, 14–16. [Google Scholar] [CrossRef]
  52. Gebauer, D.L.; Pagnussat, N.; Piato, Â.L.; Schaefer, I.C.; Bonan, C.D.; Lara, D.R. Effects of anxiolytics in zebrafish: Similarities and differences between benzodiazepines, buspirone and ethanol. Pharmacol. Biochem. Behav. 2011, 99, 480–486. [Google Scholar] [CrossRef]
  53. Finn, J.; Hui, M.; Li, V.W.T.; Lorenzi, V.; de la Paz, N.; Cheng, S.H.; Lai-Chan, L.; Schlenk, D. Effects of propranolol on heart rate and development in Japanese medaka (Oryzias latipes) and zebrafish (Danio rerio). Aquat. Toxicol. 2012, 122, 214–221. [Google Scholar] [CrossRef]
  54. Huggett, D.; Brooks, B.; Peterson, B.; Foran, C.; Schlenk, D. Toxicity of select beta adrenergic receptor-blocking pharmaceuticals (b-blockers) on aquatic organisms. Arch. Environ. Contam. Toxicol. 2002, 43, 229–235. [Google Scholar] [CrossRef]
  55. Sweetman, S.; Martindale, I. The complete drug reference pharmaceutical press. J. Med. Libr. Assoc. 2002, 100, 75–76. [Google Scholar]
  56. Jozefowski, S.J.; Plytycz, B. Characterization of β-adrenergic receptors in fish and amphibian lymphoid organs. Dev. Comp. Immunol. 1998, 22, 587–603. [Google Scholar] [CrossRef]
  57. Fraysse, B.; Mons, R.; Garric, J. Development of a zebrafish 4-day embryo-larval bioassay to assess toxicity of chemicals. Ecotoxicol. Environ. Saf. 2006, 63, 253–267. [Google Scholar] [CrossRef] [PubMed]
  58. Owen, S.F.; Giltrow, E.; Huggett, D.B.; Hutchinson, T.H.; Saye, J.; Winter, M.J.; Sumpter, J.P. Comparative physiology, pharmacology and toxicology of β-blockers: Mammals versus fish. Aquat. Toxicol. 2007, 82, 145–162. [Google Scholar] [CrossRef] [PubMed]
  59. Ternes, T.A. Occurrence of drugs in German sewage treatment plants and rivers. Water Res. 1998, 32, 3245–3260. [Google Scholar] [CrossRef]
  60. Santos, L.; Araujo, A.; Fachini, A.; Pena, A.; Delerue-Matos, C.; Montenegro, M. Ecotoxicological aspects related to the presence of pharmaceuticals in the aquatic environment. J. Hazard. Mater. 2010, 175, 45–95. [Google Scholar] [CrossRef] [PubMed]
  61. Black, J.W.; Crowther, A.F.; Shanks, R.G.; Smith, L.H.; Dornhorst, A.C. A new adrenergic betareceptor antagonist. Lancet 1964, 1, 1080–1081. [Google Scholar] [CrossRef]
  62. Margiotta-Casaluci, L.; Owen, S.F.; Rand-Weaver, M.; Winter, M.J. Testing the translational power of the zebrafish: An interspecies analysis of responses to cardiovascular drugs. Front. Pharmacol. 2019, 10, 893. [Google Scholar] [CrossRef]
  63. Marchant-Forde, J.; Lay, D., Jr.; Pajor, E.; Richert, B.; Schinckel, A. The effects of ractopamine on the behavior and physiology of finishing pigs. J. Anim. Sci. 2003, 81, 416–422. [Google Scholar] [CrossRef] [PubMed]
  64. Lonare, S.; Sole, S.; Umap, S.J. Exposure to ractopamine induces behavioural and reproductive alterations in zebrafish (danio rerio). Toxicol. Int. 2018, 25, 31–42. [Google Scholar]
  65. Lopes, E.; De Sousa, R.V.; Zangeronimo, M.; Pereira, A.N.D.J.; Coelho, M.D.R.; Ferreira, M.S.D.S.; Lima, R.R.; Marcondes, F.; Napimoga, M.H.; Pereira, L.J. Metabolic and behavioral effects of ractopamine at continuous low levels in rats under stress. Braz. Arch. Biol. Technol. 2015, 58, 406–413. [Google Scholar] [CrossRef]
  66. Srinivasan, A. Propranolol: A 50-year historical perspective. Ann. Indian Acad. Neurol. 2019, 22, 21–26. [Google Scholar] [CrossRef] [PubMed]
  67. Prichard, B.N.C.; Gillam, P.M.C. Use of Propranolol (Inderal) in Treatment of Hypertension. BMJ 1964, 2, 725–727. [Google Scholar] [CrossRef]
  68. Hayes, P.E.; Schulz, S.C. Beta-blockers in anxiety disorders. J. Affect. Disord. 1987, 13, 119–130. [Google Scholar] [CrossRef]
  69. Frese, D.A.; Reinhardt, C.D.; Bartle, S.J.; Rethorst, D.N.; Bawa, B.; Thomason, J.D.; Loneragan, G.H.; Thomson, D.U. Effect of ractopamine hydrochloride and zilpaterol hydrochloride on cardiac electrophysiologic and hematologic variables in finishing steers. J. Am. Veter. Med Assoc. 2016, 249, 668–677. [Google Scholar] [CrossRef]
  70. Schwerte, T.; Prem, C.; Mairösl, A.; Pelster, B. Development of the sympatho-vagal balance in the cardiovascular system in zebrafish (Danio rerio) characterized by power spectrum and classical signal analysis. J. Exp. Biol. 2006, 209, 1093–1100. [Google Scholar] [CrossRef]
  71. Brown, D.; Ryan, K.; Daniel, Z.; Mareko, M.; Talbot, R.; Moreton, J.; Giles, T.C.B.; Emes, R.; Hodgman, C.; Parr, T.; et al. The Beta-adrenergic agonist, Ractopamine, increases skeletal muscle expression of Asparagine Synthetase as part of an integrated stress response gene program. Sci. Rep. 2018, 8, 15915. [Google Scholar] [CrossRef] [PubMed]
  72. Kim, H.M.; Suman, S.; Li, S.; Nair, M.N.; Beach, C.M.; Edenburn, B.M.; Boler, D.D.; Dilger, A.C.; Felix, T.L. Ractopamine Influences Muscle Proteome Profile of Postmortem Beef Longissimus Lumborum. Meat Muscle Biol. 2017, 1, 140. [Google Scholar] [CrossRef]
  73. Barrett, A.; Cullum, V.A. The biological properties of the optical isomers of propranolol and their effects on cardiac arrhythmias. Br. J. Pharmacol. 1968, 34, 43–55. [Google Scholar] [CrossRef]
  74. Ruuskanen, J.O.; Xhaard, H.; Marjamäki, A.; Salaneck, E.; Salminen, T.; Yan, Y.-L.; Postlethwait, J.H.; Johnson, M.S.; Larhammar, D.; Scheinin, M. Identification of Duplicated Fourth α2-Adrenergic Receptor Subtype by Cloning and Mapping of Five Receptor Genes in Zebrafish. Mol. Biol. Evol. 2004, 21, 14–28. [Google Scholar] [CrossRef]
  75. Wang, Z.; Nishimura, Y.; Shimada, Y.; Umemoto, N.; Hirano, M.; Zang, L.; Oka, T.; Sakamoto, C.; Kuroyanagi, J.; Tanaka, T. Zebrafish β-adrenergic receptor mRNA expression and control of pigmentation. Gene 2009, 446, 18–27. [Google Scholar] [CrossRef]
  76. Johnston, S.; Staines, D.; Klein, A.; Marshall-Gradisnik, S. A targeted genome association study examining transient receptor potential ion channels, acetylcholine receptors, and adrenergic receptors in Chronic Fatigue Syndrome/Myalgic Encephalomyelitis. BMC Med. Genet. 2016, 17, 79. [Google Scholar] [CrossRef] [PubMed]
  77. Shappell, N.; Feil, V.; Smith, D.; Larsen, G.; McFarland, D. Response of c2c12 mouse and turkey skeletal muscle cells to the β-adrenergic agonist ractopamine. J. Anim. Sci. 2000, 78, 699–708. [Google Scholar] [CrossRef] [PubMed]
  78. Maltin, C.A.; Delday, M.I.; Hay, S.M.; Smith, F.G.; Reeds, P.J. Propranolol apparently separates the physical and compositional characteristics of muscle growth induced by clenbuterol. Biosci. Rep. 1987, 7, 51–57. [Google Scholar] [CrossRef] [PubMed]
  79. Futaki, S.; Goto, Y.; Ohgoshi, Y.; Yaku, H.; Kawaguchi, O.; Suga, H. Denopamine (beta1-selective adrenergic receptor agonist) and isoproterenol (non-selective beta-adrenergic receptor agonist) equally increase heart rate and myocardial oxygen consumption in dog heart. Jpn. Circ. J. 1991, 55, 972–982. [Google Scholar] [CrossRef] [PubMed]
  80. Kossack, M.; Hein, S.; Juergensen, L.; Siragusa, M.; Benz, A.; Katus, H.A.; Most, P.; Hassel, D. Induction of cardiac dysfunction in developing and adult zebrafish by chronic isoproterenol stimulation. J. Mol. Cell. Cardiol. 2017, 108, 95–105. [Google Scholar] [CrossRef]
Cells 10 02449 g001 550
Figure 1. The RAC effect on zebrafish photomotor activities. (A) Total distance traveled of zebrafish larvae after 24-h incubation to different doses of RAC in (A) light cycle and (B) dark cycle. Total rotation movement of zebrafish larvae after 24-h incubation to different doses of RAC in (C) light cycle and (D) dark cycle. Total burst movement of zebrafish larvae after exposed to different doses of RAC for 24 h in (E) light cycle and (F) dark cycle. The data were shown as the Mean ± SEM and analyzed by Kruskal–Wallis test followed with Dunn’s multiple comparisons test. Different letters a, b, c, d, e, and f above columns indicate significant statistical differences with p < 0.05 (n = 96 for all groups). Black bar: control, red bar: RAC 0.1 ppm, green bar: RAC 1 ppm, purple bar: RAC 2 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Figure 1. The RAC effect on zebrafish photomotor activities. (A) Total distance traveled of zebrafish larvae after 24-h incubation to different doses of RAC in (A) light cycle and (B) dark cycle. Total rotation movement of zebrafish larvae after 24-h incubation to different doses of RAC in (C) light cycle and (D) dark cycle. Total burst movement of zebrafish larvae after exposed to different doses of RAC for 24 h in (E) light cycle and (F) dark cycle. The data were shown as the Mean ± SEM and analyzed by Kruskal–Wallis test followed with Dunn’s multiple comparisons test. Different letters a, b, c, d, e, and f above columns indicate significant statistical differences with p < 0.05 (n = 96 for all groups). Black bar: control, red bar: RAC 0.1 ppm, green bar: RAC 1 ppm, purple bar: RAC 2 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Cells 10 02449 g001
Cells 10 02449 g002 550
Figure 2. Oxygen consumption of zebrafish larvae at 5 dpf after 24-h incubation in different doses of RAC. (A) Experiment setting used to oxygen consumption measurement by using commercial instrument of Loligo Microplate Respirometry System. (B) The time chronology of dissolved oxygen level measured in the microplate. (C) The relative oxygen consumption for zebrafish embryos treated with different dose of RAC. The data are shown as the means ± SEM and analyzed by Two-way ANOVA either with Kruskal–Wallis One-way ANOVA test with Dunn’s correction as post-hoc multiple comparison test (n = 43–46). Different letters a and b above columns indicate significant statistical differences with p < 0.05. Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Figure 2. Oxygen consumption of zebrafish larvae at 5 dpf after 24-h incubation in different doses of RAC. (A) Experiment setting used to oxygen consumption measurement by using commercial instrument of Loligo Microplate Respirometry System. (B) The time chronology of dissolved oxygen level measured in the microplate. (C) The relative oxygen consumption for zebrafish embryos treated with different dose of RAC. The data are shown as the means ± SEM and analyzed by Two-way ANOVA either with Kruskal–Wallis One-way ANOVA test with Dunn’s correction as post-hoc multiple comparison test (n = 43–46). Different letters a and b above columns indicate significant statistical differences with p < 0.05. Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Cells 10 02449 g002
Cells 10 02449 g003 550
Figure 3. Cardiac physiology analysis of zebrafish larvae at 72 hpf after exposed to different doses of RAC. Several important parameters like (A) heart rate in the ventricle, (B) heart rate in the atrium, (C) stroke volume, (D) cardiac output, (E) ejection fraction and (F) shortening fraction were measured and compared. The data are expressed as the means ± SEM and were analyzed by one-way ANOVA test followed with Fisher’s LSD test. Different letters a, b, and c above columns indicate significant statistical differences with p < 0.05 (n = 27 for all groups). Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Figure 3. Cardiac physiology analysis of zebrafish larvae at 72 hpf after exposed to different doses of RAC. Several important parameters like (A) heart rate in the ventricle, (B) heart rate in the atrium, (C) stroke volume, (D) cardiac output, (E) ejection fraction and (F) shortening fraction were measured and compared. The data are expressed as the means ± SEM and were analyzed by one-way ANOVA test followed with Fisher’s LSD test. Different letters a, b, and c above columns indicate significant statistical differences with p < 0.05 (n = 27 for all groups). Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Cells 10 02449 g003
Cells 10 02449 g004 550
Figure 4. Cardiac rhythm analysis of zebrafish larvae at 72 hpf after 24-h exposure to different doses of RAC. Several important parameters like (A) SD1 of Poincare plot in atrium, (B) SD2 of Poincare plot in atrium, (C) SD1 of Poincare plot in ventricle, (D) SD2 of Poincare plot in ventricle, (E) A–V Interval and (F) V–A Interval were measured and compared. The data are expressed as the means ± SEM and were analyzed by Kruskal–Wallis test followed with Dunn’s multiple comparisons test (n = 27). Same letter ‘a’ above columns indicate no significant statistical differences with p > 0.05. Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Figure 4. Cardiac rhythm analysis of zebrafish larvae at 72 hpf after 24-h exposure to different doses of RAC. Several important parameters like (A) SD1 of Poincare plot in atrium, (B) SD2 of Poincare plot in atrium, (C) SD1 of Poincare plot in ventricle, (D) SD2 of Poincare plot in ventricle, (E) A–V Interval and (F) V–A Interval were measured and compared. The data are expressed as the means ± SEM and were analyzed by Kruskal–Wallis test followed with Dunn’s multiple comparisons test (n = 27). Same letter ‘a’ above columns indicate no significant statistical differences with p > 0.05. Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Cells 10 02449 g004
Cells 10 02449 g005 550
Figure 5. Evaluation of blood flow velocity in zebrafish larvae after exposed to RAC. (A) Schematic diagram showing the position for image capture to calculate blood flow velocity. (B) Magnification of red circle in (A) to show the position of dorsal aorta and post cardinal vein. (C) The time chronology of blood flow rate monitoring in the dorsal aorta of zebrafish embryos. (D) Average and (E) maximal blood flow velocity in the dorsal aorta of zebrafish larvae aged at 72-h post fertilization (hpf) after 24 h incubations in RAC from 48 hpf onwards. The data are shown as the means ± SD and were analyzed by one-way ANOVA test followed with Fisher’s LSD test. Different letters a and b above columns indicate significant statistical differences with p < 0.05 (n = 21 for control, n = 20 for all RAC groups). Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Figure 5. Evaluation of blood flow velocity in zebrafish larvae after exposed to RAC. (A) Schematic diagram showing the position for image capture to calculate blood flow velocity. (B) Magnification of red circle in (A) to show the position of dorsal aorta and post cardinal vein. (C) The time chronology of blood flow rate monitoring in the dorsal aorta of zebrafish embryos. (D) Average and (E) maximal blood flow velocity in the dorsal aorta of zebrafish larvae aged at 72-h post fertilization (hpf) after 24 h incubations in RAC from 48 hpf onwards. The data are shown as the means ± SD and were analyzed by one-way ANOVA test followed with Fisher’s LSD test. Different letters a and b above columns indicate significant statistical differences with p < 0.05 (n = 21 for control, n = 20 for all RAC groups). Black bar: control, red bar: RAC 0.1 ppm, orange bar: RAC 4 ppm, and blue bar: RAC 8 ppm.
Cells 10 02449 g005
Cells 10 02449 g006 550
Figure 6. Molecular docking of RAC and zebrafish endogenous β-adrenergic receptor. (A) The identified binding pocket in the middle of the hollow cylinder of the seven transmembrane domain sub-types homology modeling structures of β-adrenergic receptor (adra1aa as an example). Three-dimensional (3D) (B) and 2D (C) illustrations of interactions between RAC and endogenous zebrafish β-adrenergic receptor showing the three hydrogen bond formations between RAC and adra1aa at the positions Glu101, Ser418, and Asp120 (highlighted by the green and blue dotted line).
Figure 6. Molecular docking of RAC and zebrafish endogenous β-adrenergic receptor. (A) The identified binding pocket in the middle of the hollow cylinder of the seven transmembrane domain sub-types homology modeling structures of β-adrenergic receptor (adra1aa as an example). Three-dimensional (3D) (B) and 2D (C) illustrations of interactions between RAC and endogenous zebrafish β-adrenergic receptor showing the three hydrogen bond formations between RAC and adra1aa at the positions Glu101, Ser418, and Asp120 (highlighted by the green and blue dotted line).
Cells 10 02449 g006
Cells 10 02449 g007 550
Figure 7. Locomotory activity evaluation after co-incubation of different doses of RAC and PROP in zebrafish larvae after 24-h incubation. Totally, six endpoints of (A) Light cycle distance, (B) Dark cycle distance, (C) Light cycle rotation, (D) Dark cycle rotation, (E) Light cycle burst, and (F) Dark cycle burst were measured and compared. The statistics were demonstrated as the mean ± SEM and analyzed by Kruskal–Wallis test followed with Dunn’s multiple comparisons test. Different letters a, b, c, and d above columns indicate significant statistical differences with p < 0.05 (n = 48). Abbreviations: R8 (RAC 8 ppm), R8P1 (RAC 8 ppm and PROP 1 ppm), and R8P4 (RAC 8 ppm and PROP 4 ppm). Black bar: control, blue bar: R8, cream bar: R8P1, and pink bar: R8P4.
Figure 7. Locomotory activity evaluation after co-incubation of different doses of RAC and PROP in zebrafish larvae after 24-h incubation. Totally, six endpoints of (A) Light cycle distance, (B) Dark cycle distance, (C) Light cycle rotation, (D) Dark cycle rotation, (E) Light cycle burst, and (F) Dark cycle burst were measured and compared. The statistics were demonstrated as the mean ± SEM and analyzed by Kruskal–Wallis test followed with Dunn’s multiple comparisons test. Different letters a, b, c, and d above columns indicate significant statistical differences with p < 0.05 (n = 48). Abbreviations: R8 (RAC 8 ppm), R8P1 (RAC 8 ppm and PROP 1 ppm), and R8P4 (RAC 8 ppm and PROP 4 ppm). Black bar: control, blue bar: R8, cream bar: R8P1, and pink bar: R8P4.
Cells 10 02449 g007
Cells 10 02449 g008 550
Figure 8. Cardiac physiology alterations of zebrafish larvae at 72 hpf after co-incubation of RAC and PROP at different doses. Several important endpoints like: (A) heart rate in the ventricle, (B) heart rate in the atrium, (C) stroke volume, (D) cardiac output, (E) ejection fraction and (F) shortening fraction (G) A–V interval, (H) V–A interval, (I) SD of Poincare plot in 1 atrium, (J) SD 2 of Poincare plot in atrium, (K) SD 1 of Poincare plot in ventricle, and (L) SD 2 of Poincare plot in ventricle were measured and compared. The data are expressed as the means ± SEM and were analyzed either by one-way ANOVA followed with Fisher’s LSD test or Kruskal–Wallis test with Dunnett’s correction for multiple comparison test. Different letters a, b, and c above columns indicate significant statistical differences with p < 0.05 (n = 20–22). Abbreviations: R8 (RAC 8 ppm), R8P1 (RAC 8 ppm and PROP 1 ppm), and R8P4 (RAC 8 ppm and PROP 4 ppm). Black bar: control, blue bar: R8, cream bar: R8P1, and pink bar: R8P4.
Figure 8. Cardiac physiology alterations of zebrafish larvae at 72 hpf after co-incubation of RAC and PROP at different doses. Several important endpoints like: (A) heart rate in the ventricle, (B) heart rate in the atrium, (C) stroke volume, (D) cardiac output, (E) ejection fraction and (F) shortening fraction (G) A–V interval, (H) V–A interval, (I) SD of Poincare plot in 1 atrium, (J) SD 2 of Poincare plot in atrium, (K) SD 1 of Poincare plot in ventricle, and (L) SD 2 of Poincare plot in ventricle were measured and compared. The data are expressed as the means ± SEM and were analyzed either by one-way ANOVA followed with Fisher’s LSD test or Kruskal–Wallis test with Dunnett’s correction for multiple comparison test. Different letters a, b, and c above columns indicate significant statistical differences with p < 0.05 (n = 20–22). Abbreviations: R8 (RAC 8 ppm), R8P1 (RAC 8 ppm and PROP 1 ppm), and R8P4 (RAC 8 ppm and PROP 4 ppm). Black bar: control, blue bar: R8, cream bar: R8P1, and pink bar: R8P4.
Cells 10 02449 g008
Cells 10 02449 g009 550
Figure 9. Evaluation of blood flow velocity and oxygen consumption in zebrafish larvae after co-incubating RAC and PROP. (A) Maximal and (B) average blood flow velocity in the dorsal aorta of zebrafish larvae aged at 72-h post fertilization (hpf) after 24 h incubations in RAC–PROP from 48 hpf onwards. The data are shown as the means ± SEM and were analyzed by One-way ANOVA test with Fisher’s LSD test (n Control = 41, n R4 = 40, n R4P1 = 20, and n R4P4 = 20). (C) Oxygen consumption of zebrafish larvae at 5 dpf after 24-h incubation in different doses of RAC and RAC-PROP combination. The data are shown as the means ± SD and analyzed by One-way ANOVA test with Dunn’s correction as post-hoc multiple comparison test (n = 42–46). Different letters a, b, and c above columns indicate significant statistical differences with p < 0.05. Abbreviations: R4 (RAC 4 ppm), R4P1 (RAC 4 ppm, PROP 1 ppm), R4P4 (RAC 4 ppm, PROP 4 ppm), R8 (RAC 8 ppm), R8P1 (RAC 8 ppm, PROP 1 ppm), and R8P4 (RAC 8 ppm, PROP 4 ppm). Black bar: control, orange bar: R4, silver bar: R4P1, gold bar: R4P4, blue bar: R8, cream bar: R8P1, and pink bar: R8P4.
Figure 9. Evaluation of blood flow velocity and oxygen consumption in zebrafish larvae after co-incubating RAC and PROP. (A) Maximal and (B) average blood flow velocity in the dorsal aorta of zebrafish larvae aged at 72-h post fertilization (hpf) after 24 h incubations in RAC–PROP from 48 hpf onwards. The data are shown as the means ± SEM and were analyzed by One-way ANOVA test with Fisher’s LSD test (n Control = 41, n R4 = 40, n R4P1 = 20, and n R4P4 = 20). (C) Oxygen consumption of zebrafish larvae at 5 dpf after 24-h incubation in different doses of RAC and RAC-PROP combination. The data are shown as the means ± SD and analyzed by One-way ANOVA test with Dunn’s correction as post-hoc multiple comparison test (n = 42–46). Different letters a, b, and c above columns indicate significant statistical differences with p < 0.05. Abbreviations: R4 (RAC 4 ppm), R4P1 (RAC 4 ppm, PROP 1 ppm), R4P4 (RAC 4 ppm, PROP 4 ppm), R8 (RAC 8 ppm), R8P1 (RAC 8 ppm, PROP 1 ppm), and R8P4 (RAC 8 ppm, PROP 4 ppm). Black bar: control, orange bar: R4, silver bar: R4P1, gold bar: R4P4, blue bar: R8, cream bar: R8P1, and pink bar: R8P4.
Cells 10 02449 g009
Table
Table 1. Flexible molecular docking score for each homology modeling structure of endogenous β-adrenergic receptor subtype. The zebrafish subtype, adra1aa, binding to RAC, gains the highest CDOCKER scores, 47.63, while the human adra2a (6kux) gains 44.15. The zebrafish subtype, adra1aa, binding to PROP with a chiral center, gains the highest CDOCKER scores for S form, 27.33, while the human adra2a (6kux) gains 21.83 and 20.62, in R form and S form, respectively. PROP overall has lower binding affinity than RAC.
Table 1. Flexible molecular docking score for each homology modeling structure of endogenous β-adrenergic receptor subtype. The zebrafish subtype, adra1aa, binding to RAC, gains the highest CDOCKER scores, 47.63, while the human adra2a (6kux) gains 44.15. The zebrafish subtype, adra1aa, binding to PROP with a chiral center, gains the highest CDOCKER scores for S form, 27.33, while the human adra2a (6kux) gains 21.83 and 20.62, in R form and S form, respectively. PROP overall has lower binding affinity than RAC.
ZFIN Gene IDGene NameChromosomal PositionCDOCKER Score
Ractopamine		CDOCKER Score
Propranolol (R)		CDOCKER Score
Propranolol (S)
ZDB-GENE-030131-2831adra1aa847.6325.3725.97
ZDB-GENE-060503-384adra1ab1044.7722.9822.40
ZDB-GENE-120510-1adra1ba2142.8522.5122.41
ZDB-GENE-041114-51adra1bb1444.8325.5827.33
ZDB-GENE-090312-203adra1d142.5023.1723.37
ZDB-GENE-021010-1adra2a2243.6821.5323.25
ZDB-GENE-021010-2adra2b846.9124.9324.48
ZDB-GENE-021010-3adra2c141.9024.2023.57
ZDB-GENE-021010-4adra2da1444.2724.4222.44
ZDB-GENE-021010-5adra2db2142.7425.1722.50
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abbas, K.; Saputra, F.; Suryanto, M.E.; Lai, Y.-H.; Huang, J.-C.; Yu, W.-H.; Chen, K.H.-C.; Lin, Y.-T.; Hsiao, C.-D. Evaluation of Effects of Ractopamine on Cardiovascular, Respiratory, and Locomotory Physiology in Animal Model Zebrafish Larvae. Cells 2021, 10, 2449. https://doi.org/10.3390/cells10092449

AMA Style

Abbas K, Saputra F, Suryanto ME, Lai Y-H, Huang J-C, Yu W-H, Chen KH-C, Lin Y-T, Hsiao C-D. Evaluation of Effects of Ractopamine on Cardiovascular, Respiratory, and Locomotory Physiology in Animal Model Zebrafish Larvae. Cells. 2021; 10(9):2449. https://doi.org/10.3390/cells10092449

Chicago/Turabian Style

Abbas, Kumail, Ferry Saputra, Michael Edbert Suryanto, Yu-Heng Lai, Jong-Chin Huang, Wen-Hao Yu, Kelvin H.-C. Chen, Ying-Ting Lin, and Chung-Der Hsiao. 2021. "Evaluation of Effects of Ractopamine on Cardiovascular, Respiratory, and Locomotory Physiology in Animal Model Zebrafish Larvae" Cells 10, no. 9: 2449. https://doi.org/10.3390/cells10092449

APA Style

Abbas, K., Saputra, F., Suryanto, M. E., Lai, Y. -H., Huang, J. -C., Yu, W. -H., Chen, K. H. -C., Lin, Y. -T., & Hsiao, C. -D. (2021). Evaluation of Effects of Ractopamine on Cardiovascular, Respiratory, and Locomotory Physiology in Animal Model Zebrafish Larvae. Cells, 10(9), 2449. https://doi.org/10.3390/cells10092449

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop