Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates.

In this work, the effect of the range of dispersive interactions in determining the three-phase coexistence line of the CO2 and CH4 hydrates has been studied. In particular, the temperature (T3) at which solid hydrate, water, and liquid CO2/gas CH4 coexist has been determined through molecular dynamics simulations using different cutoff values (from 0.9 to 1.6 nm) for dispersive interactions. The T3 of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures T. If the hydrate melts, then T > T3. Conversely, if the hydrate grows, then T < T3. The effect of the cutoff distance on the dissociation temperature has been analyzed at three different pressures for CO2 hydrate: 100, 400, and 1000 bar. Then, we have changed the guest and studied the effect of the cutoff distance on the dissociation temperature of the CH4 hydrate at 400 bar. Moreover, the effect of long-range corrections for dispersive interactions has been analyzed by running simulations with homo- and inhomogeneous corrections and a cutoff value of 0.9 nm. The results obtained in this work highlight that the cutoff distance for the dispersive interactions affects the stability conditions of these hydrates. This effect is enhanced when the pressure is decreased, displacing the T3 about 2-4 K depending on the system and the pressure.

Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate.

Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.

Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate.

In this work, the effects of finite size on the determination of the three-phase coexistence temperature (T3) of the carbon dioxide (CO2) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved (hydrate-aqueous solution-liquid CO2) are placed together in the same simulation box. By varying the number of molecules of each phase, it is possible to analyze the effect of simulation size and stoichiometry on the T3 determination. In this work, we have determined the T3 value at 8 different pressures (from 100 to 6000 bar) and using 6 different simulation boxes with different numbers of molecules and sizes. In two of these configurations, the ratio of the number of water and CO2 molecules in the aqueous solution and the liquid CO2 phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO2 in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO2 available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above T3, overestimating the T3 value obtained from direct coexistence simulations. The simulation results obtained for the CO2 hydrate confirm the sensitivity of T3 depending on the size and composition of the system, explaining the discrepancies observed in the original work by Míguez et al. [J. Chem Phys. 142, 124505 (2015)]. Non-stoichiometric configurations with larger unit cells show a convergence of T3 values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations.

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log10 (J/(m3 s)) = -7(5).

Intestinal variation of serotonin, melatonin, and digestive enzymes activities along food passage time through GIT in Salmo salar fed with supplemented diets with tryptophan and melatonin.

In teleosts, peripheral serotonin (5-HT) and melatonin (MEL) are synthesised in the gastrointestinal tract (GIT) and regulate secretion and motility processes. Their production is regulated by diet and the passage of food through the GIT. This study aimed to evaluate how intestinal 5-HT, melatonin, and the activity of digestive enzymes varied with food passage time through GIT in Atlantic salmon (Salmo salar). We fed fish diets supplemented with tryptophan and melatonin (L-Trp 2.5% and MEL 0.01%) and measured the activity of digestive enzymes (amylase, lipase, and total protease) in the pyloric caeca, midgut, and hindgut at different times after feeding. 5-HT levels increased in all GIT portions and diets at 120 min post-intake and were highest in the pyloric caeca. Intestinal enzymatic activity was varied with diet, post-intake time and in different intestinal portions. In conclusion, food passage time directly affects GIT 5-HT secretion and digestive enzyme activity in S. salar, and diet composition regulates S. salar GIT function.

Computational study of the interplay between intermolecular interactions and CO2 orientations in type I hydrates.

Carbon dioxide (CO2) molecules show a rich orientation landscape when they are enclathrated in type I hydrates. Previous studies have described experimentally their preferential orientations, and some theoretical works have explained, but only partially, these experimental results. In the present paper, we use classical molecular dynamics and electronic density functional theory to advance in the theoretical description of CO2 orientations within type I hydrates. Our results are fully compatible with those previously reported, both theoretical and experimental, the geometric shape of the cavities in hydrate being, and therefore, the steric constraints, responsible for some (but not all) preferential angles. In addition, our calculations also show that guest-guest interactions in neighbouring cages are a key factor to explain the remaining experimental angles. Besides the implication concerning equation of state hydrate modeling approximations, the conclusion is that these guest-guest interactions should not be neglected, contrary to the usual practice.

On interfacial properties of tetrahydrofuran: Atomistic and coarse-grained models from molecular dynamics simulation.

We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombic intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.

Molecular dynamics simulation of CO2 hydrates: Prediction of three phase coexistence line.

The three phase equilibrium line (hydrate-liquid water-liquid carbon dioxide) has been estimated for the water + carbon dioxide binary mixture using molecular dynamics simulation and the direct coexistence technique. Both molecules have been represented using rigid nonpolarizable models. TIP4P/2005 and TIP4P/Ice were used for the case of water, while carbon dioxide was considered as a three center linear molecule with the parameterizations of MSM, EPM2, TraPPE, and ZD. The influence of the initial guest occupancy fraction on the hydrate stability has been analyzed first in order to determine the optimal starting configuration for the simulations, paying attention to the influence of the two different cells existing in the sI hydrate structure. The three phase coexistence temperature was then determined for a pressure range from 2 to 500 MPa. The qualitative shape of the equilibrium curve estimated is correct, including the high pressure temperature maximum that determines the hydrate re-entrant behaviour. However, in order to obtain quantitative agreement with experimental results, a positive deviation from the classical Lorentz-Berthelot combining rules must be considered.

Influence of the long-range corrections on the interfacial properties of molecular models using Monte Carlo simulation.

We analyze the influence of the long-range corrections, due to the dispersive term of the intermolecular potential energy, on the surface tension using direct simulation of the vapour-liquid interface of different molecular models. Although several calculation methods have been proposed recently to compute the fluid-fluid interfacial properties, the truncation of the intermolecular potential or the use of the tail corrections represents a contribution relevant from a quantitative perspective. In this work, a simplified model for methane, namely a spherical Lennard-Jones intermolecular potential, has been considered first, and afterwards other models including rigid non polarizable structures with both Lennard-Jones sites and point electric charges, representing some of the most popular models to describe water (namely the original TIP4P model, and the TIP4P/Ew and TIP4P/2005 versions), and carbon dioxide (MSM, EPM2, TraPPE, and ZD models) have been studied. Our results show that for all cases tested, including those in which the electrostatic interactions may be predominant, an incomplete account of the long-range corrections produces a systematic underestimation of the computed interfacial tension.

On interfacial tension calculation from the test-area methodology in the grand canonical ensemble.

We propose the extension of the test-area methodology, originally proposed to evaluate the surface tension of planar fluid-fluid interfaces along a computer simulation in the canonical ensemble, to deal with the solid-fluid interfacial tension of systems adsorbed on slitlike pores using the grand canonical ensemble. In order to check the adequacy of the proposed extension, we apply the method for determining the density profiles and interfacial tension of spherical molecules adsorbed in slitlike pore with different pore sizes and solid-fluid dispersive energy parameters along the same simulation. We also calculate the solid-fluid interfacial tension using the original test-area method in the canonical ensemble. Agreement between the results obtained from both methods indicate that both methods are fully equivalent. The advantage of the new methodology is that allows to calculate simultaneously the density profiles and the amount of molecules adsorbed onto a slitlike pore, as well as the solid-fluid interfacial tension. This ensures that the chemical potential at which all properties are evaluated during the simulation is exactly the same since simulations can be performed in the grand canonical ensemble, mimicking the conditions at which the adsorption experiments are most usually carried out in the laboratory.

Vapor-liquid interfacial properties of rigid-linear Lennard-Jones chains.

We have obtained the interfacial properties of short rigid-linear chains formed from tangentially bonded Lennard-Jones monomeric units from direct simulation of the vapour-liquid interface. The full long-range tails of the potential are accounted for by means of an improved version of the inhomogeneous long-range corrections of Janecek [J. Phys. Chem. B 110, 6264-6269 (2006)] proposed recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2009)] valid for spherical as well as for rigid and flexible molecular systems. Three different model systems comprising of 3, 4, and 5 monomers per molecule are considered. The simulations are performed in the canonical ensemble, and the vapor-liquid interfacial tension is evaluated using the test-area method. In addition to the surface tension, we also obtain density profiles, coexistence densities, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the chain length and rigidity on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapor-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The surface tension has been scaled by critical properties and represented as a function of the difference between coexistence densities relative to the critical density.

Glucosensing capacity of rainbow trout telencephalon.

To assess the hypothesis of glucosensing systems present in fish telencephalon, we first demonstrated in rainbow trout, by in situ hybridisation, the presence of glucokinase (GK). Then, we assessed the response of glucosensing markers in rainbow trout telencephalon 6 hours after i.c.v. treatment with glucose or 2-deoxyglucose (inducing glucoprivation). We evaluated the response of parameters related to the mechanisms dependent on GK, liver X receptor (LXR), mitochondrial activity, sweet taste receptor and sodium-glucose linked transporter 1 (SGLT-1). We also assessed mRNA abundance of neuropeptides involved in the metabolic control of food intake (agouti-related protein, neuropeptide Y, pro-opiomelanocortin, and cocaine- and amphetamine-related transcript), as well as the abundance and phosphorylation status of proteins possibly involved in linking glucosensing with neuropeptide expression, such as protein kinase B (AkT), AMP-activated protein kinase (AMPK), mechanistic target of rapamycin and cAMP response element-binding protein (CREB). The responses obtained support the presence in the telencephalon of a glucosensing mechanism based on GK and maybe one based on LXR, although they do not support the presence of mechanisms dependent on mitochondrial activity and SGLT-1. The mechanism based on sweet taste receptor responded to glucose but in a converse way to that characterised previously in the hypothalamus. In general, systems responded only to glucose but not to glucoprivation. Neuropeptides did not respond to glucose or glucoprivation. By contrast, the presence of glucose activates Akt and inhibits AMPK, CREB and forkhead box01. This is the first study in any vertebrate species in which the response to glucose of putative glucosensing mechanisms is demonstrated in the telencephalon. Their role might relate to processes other than homeostatic control of food intake, such as the hedonic and reward system.

Daily and photoperiodic melatonin binding changes in the suprachiasmatic nuclei, paraventricular thalamic nuclei, and pars tuberalis of the female Siberian hamster (Phodopus sungorus).

Using quantitative autoradiography, 2-(125)I-melatonin binding was investigated throughout the light:dark cycle in the suprachiasmatic nuclei (SCN), paraventricular nuclei (PVT), and pars tuberalis (PT) of adult female Siberian hamsters kept for 10 weeks in either long or short photoperiods (LP or SP, respectively). Plasma melatonin concentrations were measured by radioimmunoassay, and the sexual status of the animals was established by visual inspection of vaginal smears and by weighing uteri after sacrifice. The SCN displayed neither daily nor photoperiod-dependent variations in specific binding. Melatonin receptors in these nuclei would be regulated neither by plasma melatonin nor by the light:dark cycle or sexual steroids. By contrast, melatonin receptor density in the PT displayed both strong daily (maximal values during the first half of the light period and minimal values during the night) and photoperiod-dependent (maximal values in LP) variations. These variations dependent on changes in the maximal binding (Bmax) without differences in the dissociation constant (Kd). Daily melatonin receptor densities in the PT of LP- and SP-exposed animals might be regulated by the dark:light transition but not by melatonin. Daily profiles of 2-(125)I-melatonin-specific binding in the PT were independent of photoperiod. Factors underlying the photoperiod-related variations presently are unknown. Concerning the PVT, weak variations in specific binding were detected in SP only when time points were grouped according to the light or dark periods. It is not yet possible to conclude whether they have any physiological relevance. These results show clearly that the regulation of melatonin receptors varies among structures (SCN, PVT, and PT) in the Siberian hamster and is also totally different from that found in the rat.

Understanding the Phase Behavior of Tetrahydrofuran + Carbon Dioxide, + Methane, and + Water Binary Mixtures from the SAFT-VR Approach.

The high-pressure phase diagrams of the tetrahydrofuran(1) + carbon dioxide(2), + methane(2), and + water(2) mixtures are examined using the SAFT-VR approach. Carbon dioxide molecule is modeled as two spherical segments tangentially bonded, water is modeled as a spherical segment with four associating sites to represent the hydrogen bonding, methane is represented as an isolated sphere, and tetrahydrofuran is represented as a chain of m tangentially bonded spherical segments. Dispersive interactions are modeled using the square-well intermolecular potential. In addition, two different molecular model mixtures are developed to take into account the subtle balance between water-tetrahydrofuran hydrogen-bonding interactions. The polar and quadrupolar interactions present in water, tetrahydrofuran, and carbon dioxide are treated in an effective way via square-well potentials of variable range. The optimized intermolecular parameters are taken from the works of Giner et al. (Fluid Phase Equil. 2007, 255, 200), Galindo and Blas (J. Phys. Chem. B 2002, 106, 4503), Patel et al. (Ind. Eng. Chem. Res. 2003, 42, 3809), and Clark et al. (Mol. Phys. 2006, 104, 3561) for tetrahydrofuran, carbon dioxide, methane, and water, respectively. The phase diagrams of the binary mixtures exhibit different types of phase behavior according to the classification of van Konynenburg and Scott, ranging from types I, III, and VI phase behavior for the tetrahydrofuran(1) + carbon dioxide(2), + methane(2), and + water(2) binary mixtures, respectively. This last type is characterized by the presence of a Bancroft point, positive azeotropy, and the so-called closed-loop curves that represent regions of liquid-liquid immiscibility in the phase diagram. The system exhibits lower critical solution temperatures (LCSTs), which denote the lower limit of immiscibility together with upper critical solution temperatures (UCSTs). This behavior is explained in terms of competition between the incompatibility with the alkyl parts of the tetrahydrofuran ring chain and the hydrogen bonding between water and the ether group. A minimum number of unlike interaction parameters are fitted to give the optimal representation of the most representative features of the binary phase diagrams. In the particular case of tetrahydrofuran(1) + water(2), two sets of intermolecular potential model parameters are proposed to describe accurately either the hypercritical point associated with the closed-loop liquid-liquid immiscibility region or the location of the mixture lower- and upper-critical end-points. The theory is not only able to predict the type of phase behavior of each mixture, but also provides a reasonably good description of the global phase behavior whenever experimental data are available.

Gradation of the stress response in rainbow trout exposed to stressors of different severity: the role of brain serotonergic and dopaminergic systems.

After an intense acute stressor, fish develop a metabolic and behavioural response that usually lasts for several hours. Brain monoaminergic systems, particularly the serotonergic system, appear to play a key role in the central regulation of the stress response. However, the influence of stressor severity on brain monoaminergic systems and on the induced stress responses is yet poorly understood. We hypothesise that serotonergic system could have a direct role in the integration of sensory information during stressor exposure and in the organisation of the subsequent integrated stress response. According to our hypothesis, a low stressor intensity would induce a low response of brain serotonergic system and therefore stress responses of low magnitude and duration. To test this hypothesis, we exposed fish to handling disturbance for 5 s, 15 s or 3 min. We sampled fish at 0 (controls), 3, 15, 45 and 240 min after the start of the stress protocol. Brain levels of serotonin, dopamine and their respective main oxidative metabolites were quantified, along with plasma levels of stress markers (catecholamines, cortisol, glucose and lactate). Regarding stress markers, the 5-s and 15-s stress protocols induced similar and relatively low elevations in all parameters assessed. As expected, the 3-min protocol induced responses of a higher intensity and duration in all plasma parameters. Interestingly, the alterations of brain monoaminergic systems did not follow the same trend. The three stress protocols induced increases in the serotonergic activity in all brain regions analysed (hypothalamus, telencephalon and medulla oblongata), independently of the duration of the handling disturbance, whereas the effects on the dopaminergic system were minor and brain region-dependent. These data suggest that the brain serotonergic system, although likely involved in the recognition of the stressor stimuli, is not the only actor determining the magnitude and duration of the acute stress response in trout.

Effects of insulin treatment on the response to oleate and octanoate of food intake and fatty acid-sensing systems in rainbow trout.

We hypothesized that food intake and the response of fatty acid (FA)-sensing systems in hypothalamus, liver, and Brockmann bodies of rainbow trout to raised levels of oleate (OL) or octanoate (OCT) is modified by insulin treatment. To assess this hypothesis, 15 fish per group received intraperitoneally 10-mL/kg injection of saline solution alone (control), or containing insulin (2-mg bovine insulin/kg body mass), OL (300 µg/kg), OCT (300 µg/kg), insulin + OL, or insulin + OCT to be sampled 6 h later to assess parameters related to FA sensing. Our results suggest that the modulatory role of insulin on the responses of hypothalamic FA-sensing systems to changes in circulating levels of OL or OCT was of minor importance in contrast to the mammalian model. However, this is in contrast with the effects observed in another experiment assessing changes in food intake after similar treatments because insulin treatment enhanced the anorectic effects of FA alone, and the effect was especially relevant (P < 0.001) for OCT, in contrast with the mammalian model where this FA is not inducing an anorectic response. In liver and Brockmann bodies, insulin treatment enhanced the responses to OL or OCT treatment in parameters related to FA sensing. Therefore, we provide for the first time in fish, and in a non-mammalian vertebrate, evidence for the modulation of FA-sensing systems by insulin.

Peptidergic modulation of serotonin release from cultured rat pinealocytes.

This paper describes the effects of beta-adrenergic and peptidergic inputs on serotonin (5-HT) synthesis, outflow and metabolism into melatonin in cultured dissociated rat pinealocytes. The spontaneous outflow of 5-HT from pinealocytes was high as demonstrated by the elevated levels of extracellular 5-HT accumulated in the medium (about 5 ng/h/70,000 pineal cells). The beta-adrenergic agonist isoproterenol (ISO) used at concentrations up to 10(-6) M induced a moderate (+20-40%) increase in intra- and extracellular 5-HT levels together with a large release of melatonin. At a higher ISO stimulation (10(-5) M), the intra- and extracellular levels of 5-HT were significantly (-25-30%) reduced whereas melatonin secretion was dramatically increased. This is interpreted as a large 5-HT mobilization for melatonin synthesis and release, consequently reducing both the intracellular pool and outflow of 5-HT. The peptides vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase activating peptide (PACAP) up to 10(-7) M induced always a moderate (+20-30%) increase in intra- and extracellular levels of 5-HT. However, the use of nM concentrations of VIP or PACAP together with 10(-6) M ISO induced a decrease in 5-HT outflow (-25-30%) and a dramatic increase in melatonin secretion as did 10(-5) M ISO alone. Neuropeptide Y (NPY) is another pineal peptide which induced a stimulation of 5-HT outflow (+30-40%) although its effect on melatonin release was marginal. The above results are discussed in term of the multineuronal regulation of the synthetic and secretory activities of the rat pineal gland.

Changes in pineal indoleamines in rats after single melatonin injections: evidence for a diurnal sensitivity to melatonin.

We recently determined that melatonin stimulated serotonin (5-HT) secretion from rat pineal glands by increasing 5-HT release from the pinealocytes (microM melatonin concentrations) and by inhibiting 5-HT uptake in the pineal sympathetic nerve endings (mM melatonin concentrations). The present study investigated whether a single melatonin injection could alter the content of indoleamines in the rat pineal gland, as well as its possible dependence on the daytime of administration. Melatonin (150 micrograms/kg) was i.p. injected at 8 time points (11.00 h, 14.00 h, 17.00 h, 20.00 h, 23.00 h, 02.00 h, 05.00 h and 08.00 h) to rats kept in 12:12 h light:dark cycle (lights on at 07.00 h). Melatonin injections in the afternoon (17:00 h) and late in the nighttime (02.00 h and 05.00 h) decreased pineal 5-HT content 90 min later. The levels of 5-hydroxyindoleacetic acid (5-HIAA) were also decreased 90 min after the melatonin treatment at 14.00 h, 17.00 h and 02.00 h. The effect of melatonin on 5-HT content was a long-fasting effect (still evident after 180 min) only when injected at 02.00 h, whereas 5-HIAA levels were found to be decreased 180 min after melatonin treatment at 14.00 h and 23.00 h. No changes in these compounds were detected 240 min after melatonin treatment. Moreover, melatonin did not change 5-hydroxytryptophan levels at any of the daytime points studied. By contrast, 90 min after the injection of melatonin at 20.00 h, an increased content of pineal N-acetylserotonin was observed. This effect of melatonin could be mediated through a phase alteration of the pineal N-acetyltransferase activity rhythm by acting on the suprachiasmatic clock, although a direct melatonin effect on the pineal rhythmic function cannot be excluded. The effects of the hormone on 5-HT and 5-HIAA contents agree with previous findings on the inhibitory effect of pharmacological doses of melatonin on pineal 5-HT uptake, which presumably would result in a decreased intraneuronal content of 5-HT and its acid metabolite. These data point to an acute regulatory action of exogenous melatonin on the pineal melatonin synthesis pathway which seems to be limited to two daytime phases; the afternoon-early evening period and the second half of the night.

Evidence for a regulatory role of melatonin on serotonin release and uptake in the pineal gland.

Melatonin has been proposed to exert some regulatory actions within the pineal gland itself. The present study examined the effect of melatonin on the release of serotonin (5-HT) and 5-hydroxyindoleacetic acid (5-HIAA) from rat pineal glands by using an in vitro perifusion system. Melatonin induced a concentration-dependent stimulatory effect on 5-HT secretion from 10(-6) M to 10(-3) M. Maximal effects were obtained with melatonin 10(-3) M and concentrations lower than 10(-6) M were without effect. The secretion of 5-HIAA was inhibited by melatonin 10(-3) and 10(-4) M, but it was increased when pineals were incubated with 10(-5) and 10(-6) M of melatonin. The indoleamine secretion was also studied on peripherally denervated rat pineal glands. Basal output of 5-HT from these glands was increased when compared with those from control rats. In contrast, the secretion of 5-HIAA was strongly reduced after removal of the sympathetic input to the pineal gland. Melatonin 10(-3) M failed to stimulate 5-HT release from denervated pineal glands, although it inhibited 5-HIAA secretion. In contrast, melatonin 10(-5) M enhanced 5-HT release without altering 5-HIAA output. Fluoxetine, a 5-HT uptake inhibitor, produced similar effects than mM concentrations of melatonin on the indoleamine secretion from control pineal glands, but it had no effect on glands taken from peripherally denervated rats. These data suggest that mM concentrations of the pineal hormone are able to stimulate 5-HT release from the pinealocyte, while mM concentrations of melatonin increase extracellular 5-HT by inhibiting its reuptake in the adrenergic nerve endings. These findings are discussed in relation to the possible role of melatonin regulating the intra- and extracellular availability of 5-HT in the pineal gland and its significance as an autocrine factor.

Changes in serotonin level and turnover in discrete hypothalamic nuclei after pinealectomy and melatonin administration to rats.

The influence of the pineal gland on the hypothalamic serotonergic function was examined by studying the effects of long-term pinealectomy (1 month) and melatonin replacement (500 micrograms/kg; 10 days) on serotonin (5-HT) and 5-hydroxyindoleacetic acid (5-HIAA) content as well as on the in vivo 5-HT synthesis rate in discrete hypothalamic nuclei. Pinealectomy was followed by a significant decrease of 5-HT content in the anterior hypothalamic nuclei (AHN) and the ventromedial hypothalamic nuclei (VMHN), and also in 5-HIAA content in lateral (LPON) and medial preoptic nuclei (MPON). The 5-HT synthesis rate, estimated from the accumulation of 5-hydroxytryptophan after blockade of the 1-amino acid decarboxylase activity, were also decreased in the AHN and the paraventricular hypothalamic nuclei (PVHN) of pinealectomized rats. In contrast, an enhanced 5-HT synthesis rate and basal 5-HIAA content were found in the suprachiasmatic nuclei (SCN) after pinealectomy. Daily treatment with melatonin for 10 days reversed most of the effects induced by pinealectomy. Thus, melatonin increased the levels of 5-HT in the AHN and VMHN, and slightly increased the 5-HIAA content in preoptic nuclei. In addition, melatonin increased the 5-HT synthesis rate in the AHN and VMHN, but also in the MPON, VMHN and dorsomedial hypothalamic nuclei (DMHN) where pinealectomy had no effect. By contrast, melatonin treatment did not affect SCN 5-HT synthesis rate, although it decreased 5-HIAA levels. The results demonstrate that melatonin is able to stimulate 5-HT metabolism in most of the hypothalamic areas, but inhibits SCN 5-HT function. Some of the effects of melatonin seems to be exerted by modulating the synthesis of the amine, although melatonin likely also interacts with other regulatory processes of 5-HT function (i.e. release/uptake). The well defined presence of melatonin receptors in the rat SCN, and its absence in other hypothalamic structures, suggest that this may be the mechanism mediating the differential response to endogenous melatonin. Moreover, the larger effect of exogenous melatonin in relation to pinealectomy suggests the presence of melatonin unespecific effects possibly owing to supraphysiological doses. The present findings may be relevant for the mode of action of melatonin and its implication in several endocrine and behavioral functions mediated by serotonergic neurons.