1236
Vol. 62, No. 12
arena (70cm diameter) with 30-cm high walls. Brightness at B3LYP/6–31++g(d,p) level.
the level of the maze arms and open-field centre was 60lx.
Preparation of the Receptors The structures of the
Procedure The ETM is an anxiety model that evokes two transporters of serotonin (SERT), dopamine (DAT) and
defensive responses in the same rat, namely, inhibitory avoid- norepinephrine (NET) modelled in Apo form were kindly
ance latency from the closed arm (baseline, avoidance 2 and provided by Ravna et al.26) The active sites were identified
3) and one-way escape latency from the open arm (escape 1, 2 through the overlapping of the three-dimensional structure
and 3), which have been related to generalized anxiety disor- of the LeuT submitted to crystallography with (R)-fluoxetine
der and panic disorder, respectively. Locomotion was assessed (PDB: 3GWV)42) in the three modelled transporters. Ligand
in the open field following each drug treatment as a control (R)-fluoxetine was added to the transporter structures by geo-
for nonspecific motor effects.
metric docking. LeuT from Aquifex aeolicus (pdb: 3GWV)
After 21d of treatment with 8-PN (10mg/kg), fluoxetine was chosen as a geometric docking mold for fluoxetine be-
(15mg/kg) or the vehicle were administered by gavage and be- cause this protein belongs to the Sodium Neurotransmitter
havioural tests were performed. The doses of fluoxetine were Symporter Family (SNF), the same as the transporters SERT,
chosen based on their effects in previous studies performed in DAT and NET (protein family id: PF00209), according to
the ETM.32) The doses of 8-PN were chosen based on in vivo evaluation performed by server Pfam.43) Furthermore, pdb
study using the adult rat model33) and in a previous study per- 3GWV is the only structural mold linked to (R)-fluoxetine that
formed in the ETM with a similar compound.34) On the 20th was available.
day of treatment, the animals were gently handled for 5min
(R)-Fluoxetine-receptor complexes were submitted to 60000
and pre-exposed to one of the open arms of the ETM for energy minimization steps by conjugated gradient for the re-
30min. The open arm exit was closed with a wooden barrier moval of any steric hindrance with the NAMD2 program.44)
mounted between the central area of the maze and the proxi- The objective of docking (R)-fluoxetine to the three receptors
mal end of the arm. It has been shown that such pre-exposure was to serve as a docking guide for other ligands which had
makes the escape task more sensitive to antipanic drugs, as it unknown bonding mechanisms.
shortens the withdrawal latencies from the open arm during
Ligand Docking The docking protocol was validated by
redocking assays of (R)-fluoxetine in the three minimized
the test.35) The ETM test was performed 24h later.
The ETM test started with the inhibitory avoidance task. complexes. The protocol was considered valid when the rmsd
Each animal was placed at the distal end of the enclosed calculated from the overlapping of the best pose onto the
ETM arm facing the intersection. The time that the rats took ligand was smaller than an average of 0.5Å in four assays of
to leave this arm with four paws was recorded (baseline la- each complex, in order to avoid false positive results.
tency). This measurement was repeated in two subsequent
The docking simulations for the standards (R)-fluoxetine,
trials (avoidance 1 and 2) at 30-s intervals. Thirty second after (S,S)-reboxetine and (S)-bupropion and for enantiomers
the last avoidance trial, the rats were placed at the end of the (R)-8-PN and (S)-8-PN were performed with the AutoDock
open arm to which they had been previously exposed and the 4.2.3 program implemented at the interface PyRx 0.9,45) ap-
latency to leave this arm with four paws was recorded in three plying the redocking-validated protocol (hybrid Lamarckian
consecutive trials (escape 1, 2 and 3) at 30-s intervals. A cut- Genetic Algorithm).
off time of 300s was established for avoidance and the escape
The energy evaluation grid was chosen according to the
latency. Thirty second after the ETM test, the animals were crystallographic structure of (R)-fluoxetine for each receptor
placed inside the circular arena for 5min to evaluate locomo- and was centred on coordinates (X=27.980, Y=20.251 and
tion. The total distance travelled was recorded with a video Z=24.870) in SERT, (X=28.301, Y=21.269 and Z=20.420) in
tracking system (Ethovision; Noldus, Holland) for analysis.
DAT and (X=25.639, Y=20.555 and Z=20.010) in NET, with
The results observed in the ETM and circular arena were grid points in the x, y and z axes set to 50×50×50 and sepa-
submitted to one-way ANOVA. When appropriate, the Dun- rated by 0.375Å. The initial population size and maximum
can post hoc test was used. The significance level was set at number of energy evaluations were set to 10. The docked
p<0.05. The Statistica Six Sigma software was used for statis- results within an rmsd of 2.0Å were clustered and the final
tical analysis.
results of each ligand were selected considering both the
Molecular Docking: Ligand Preparation Standards embedded empirical binding free energy evaluation and the
(R)-fluoxetine, (S,S)-reboxetine, (S)-bupropion (Fig. 4) were clustering analysis.
selected according to literature data36–39) and confirmed as
Energy Minimization The lowest energy poses obtained
the actually desired transporters using the SEA search tool.40) by docking the ligands (S,S)-reboxetine and (S)-bupropion
The three-dimensional structures of these compounds were and enantiomers (R)-8-PN and (S)-8-PN were exported, in-
obtained from the ZINC database (codes ZINC01530638, corporated into the receptor structures and submitted again
ZINC00002284 and ZINC00020228, respectively). It is note- to 60000 energy minimization steps with program NAMD2.
worthy that at pH 7.0, all standards have positively charged In the energy minimization procedures, the field force
nitrogen; the protonated version of the structures were used in CHARMM C35b2–C36a2 was adopted for the proteins, while
the molecular docking calculations.
for the ligands, it was generated in the same format by the
The enantiomers structures (R)-8-PN and (S)-8-PN were SwissParam server.46)
generated with the Gaussian 09 program.41) The stable ge-
Energy minimization was simulated with the complexes
ometries of the compounds were obtained by calculating immersed in a box with water at least 10Å from the protein
the potential energy surface (PES) through the HF/3-21G outermost surface. Either Na+ or Cl− counter ions were added
level of theory. The most stable geometries were optimized in appropriate amounts to neutralize the system charges. The
by density functional theory (DFT) calculations with the temperature and the pressure were adjusted to 300K and 1