Original Research ARTICLE
Identification of the antiepileptic racetam binding site in the synaptic vesicle protein 2A by molecular dynamics and docking simulations
- 1Laboratorio de Modelado Molecular y Diseño de fármacos, Departamento de Bioquímica de la Escuela Superior de Medicina del Instituto Politécnico Nacional, México City, Mexico
- 2Laboratorio de Nutrición Experimental, Laboratorio de Oncología Experimental and Laboratorio de Neuroquímica, Instituto Nacional de Pediatría, México City, Mexico
- 3División de Estudios de Posgrado, Facultad de Química, Universidad Autónoma de Querétaro, Santiago de Querétaro, Mexico
Synaptic vesicle protein 2A (SV2A) is an integral membrane protein necessary for the proper function of the central nervous system and is associated to the physiopathology of epilepsy. SV2A is the molecular target of the anti-epileptic drug levetiracetam and its racetam analogs. The racetam binding site in SV2A and the non-covalent interactions between racetams and SV2A are currently unknown; therefore, an in silico study was performed to explore these issues. Since SV2A has not been structurally characterized with X-ray crystallography or nuclear magnetic resonance, a three-dimensional (3D) model was built. The model was refined by performing a molecular dynamics simulation (MDS) and the interactions of SV2A with the racetams were determined by docking studies. A reliable 3D model of SV2A was obtained; it reached structural equilibrium during the last 15 ns of the MDS (50 ns) with remaining structural motions in the N-terminus and long cytoplasmic loop. The docking studies revealed that hydrophobic interactions and hydrogen bonds participate importantly in ligand recognition within the binding site. Residues T456, S665, W666, D670 and L689 were important for racetam binding within the trans-membrane hydrophilic core of SV2A. Identifying the racetam binding site within SV2A should facilitate the synthesis of suitable radio-ligands to study treatment response and possibly epilepsy progression.
Epilepsy is the most common chronic brain disorder that affects people of all ages. More than 50 million people worldwide have epilepsy (WHO, 2014). It is characterized by recurrent seizures, which may vary from a brief lapse of attention or muscle jerks, to severe and prolonged convulsions (WHO, 2014). In most patients with epilepsy, seizures respond to available medications. However, a significant number of patients -especially in the setting of medical-intractable epilepsies- may experience different degrees of memory or cognitive impairment, behavioral abnormalities or psychiatric symptoms, which limit their daily functioning. As a result, in many patients, epilepsy may resemble a neurodegenerative disease (Ono and Galanopoulou, 2012). Epileptic seizures and/or epileptogenesis may functionally alter brain regions involved in cognitive processing, contributing to the progressive nature of epilepsy; additionally neurodegenerative cellular mechanisms may also participate. To cure epilepsy, both epileptogenesis and the associated neuro-degeneration have to be stopped and, if possible, reversed. This will require early detection through biomarkers that can reliably predict disease progression (Ono and Galanopoulou, 2012).
Synaptic vesicle protein 2A is an integral membrane protein present on all synaptic vesicles, being nearly ubiquitous in the CNS (Lynch et al., 2004). The existing evidence suggests that SV2A is critical for the proper function of the CNS. It regulates vesicle exocytosis by modulating the concentration of pre-synaptic calcium (Janz et al., 1999; Chang and Sudhof, 2009; Wan et al., 2010); maybe by binding to synaptotagmin 1 (Schivell et al., 1996, 2005; Yao et al., 2010), a calcium sensor protein implicated in synaptic vesicle exocytosis. In addition, it has been proposed that SV2A may control the release of neurotransmitters as a gel matrix or could transport ions, such as chloride or calcium (Bajjalieh et al., 1994; Janz et al., 1999; Reigada et al., 2003; Mendoza-Torreblanca et al., 2013). The SV2A protein has been associated to the physiopathology of epilepsy; newborn mice lacking SV2A experience severe seizures and die within 3 weeks, suggesting multiple neural alterations (Crowder et al., 1999; Janz et al., 1999). Moreover, SV2A protein expression is significantly reduced in brain tissue obtained from epileptic patients and in rats during epileptogenesis, suggesting that decreased SV2A contributes to the progression of epilepsy (Feng et al., 2009; Toering et al., 2009; van Vliet et al., 2009). SV2A is the molecular target of the second-generation antiepileptic drug LEV and its structural analogs BRIV, SEL and UCB-30889 (Gillard et al., 2003; Lynch et al., 2004; Frycia et al., 2010). Although LEV may acutely modulate ion channels and other targets; chronically, LEV may lead to decreased transmitter release predominantly by binding to SV2A (Surges et al., 2008; Crepeau and Treiman, 2010; Lyseng-Williamson, 2011; Deshpande and Delorenzo, 2014).
Knowing the in vivo expression pattern of SV2A in epilepsy patients may be essential for evaluating increased epileptogenicity and thus disease progression (van Vliet et al., 2009). Moreover, SV2A expression could also be useful as a biomarker for treatment response (de Groot et al., 2011). For this purpose, sensitive neuroimaging methods such as PET may enhance our ability to detect SV2A expression (Ono and Galanopoulou, 2012). Because PET requires the use of a radio-ligand to label the protein target, as a first approximation to synthesize a radio-ligand with a high affinity and specificity, we investigated the specific binding site(s) for racetams in SV2A and the residues involved in their interaction. Therefore, an in silico study was performed to identify the binding site for LEV and other racetams within SV2A. Since SV2A has not been structurally characterized with X-ray crystallography or nuclear magnetic resonance, it was necessary to generate and validate a 3D model of this protein. This model was refined by a MDS and docking studies were performed to decipher the interactions between the racetams and SV2A.
Materials and Methods
The SV2A protein sequence was retrieved from the National Center for Biotechnology Information (NP_476558.2) and UniProt (Q02563). There was no difference between these databases; the protein sequence was thus used as an input for the I-Tasser server (Roy et al., 2010); which provided a cluster of five 3D models for this protein. The Rampage server was used to obtain the Ramachandran plot, which shows the φ and ψ torsion angles for all protein residues. Additionally, the CPHmodels-3.2 and 3D-JIGSAWv3.0 servers were used to confirm the reproducibility of the SV2A 3D model provided by I-Tasser.
System Preparation for Molecular Dynamics Simulation
The model of the SV2A protein was inserted into a pre-equilibrated membrane consisting of POPC molecules. The resulting system was then placed in a hexagonal prism shaped box with its symmetry axis, which was perpendicular to the plane of the bilayer membrane, in the z direction. The system was then solvated using simple point charge model water molecules. The geometry and size of the simulation box were carefully selected to reduce to a minimum the number of water molecules and thus decrease the computation time. In order to remove solvent molecules accidentally introduced in the hydrophobic region of the lipid bilayer, 19 Na+ ions were added to neutralize the whole system.
The simulation box contained one molecule of SV2A, 342 POPC lipids, 47,671 water molecules and 19 ions, resulting in a total of 172,298 atoms. The OPLS-AA force field was used, it has been used to simulate and describe several TM proteins with good results (Tieleman et al., 2006). The simulation was done by applying the half-𝜀 double-pair list method to ensure compatibility between the Berger united atom parameters, used for lipids (Lindahl and Edholm, 2000; Chakrabarti et al., 2010), with the OPLS-AA force field employed for protein and ions (Jorgensen et al., 1996; Kaminski et al., 2001).
Molecular Dynamics Simulation
The system in its entirety was minimized for 500 steps using the steepest descent optimization algorithm prior to generating the MDS trajectory. Solvent–protein and lipid–protein contacts were optimized by a simulation lasting 5 ns, with harmonic constraints in the heavy atoms of the protein (all, except for hydrogen). The force constants of these positional restraints were sequentially reduced. Afterward, an unrestrained trajectory lasting 50 ns was generated. The simulation was performed using a 2 fs time step. The pressure was set at 1 bar using a semi-isotropic Parrinello and Rahman (1981) barostat (Nosé and Klein, 1983) to allow independent modification of the box dimensions in the XY plane and in the Z dimension. The temperature was kept constant at 310 K by a Nosé (1984) and Hoover (1985) thermostat. The linear constraint solver algorithm was employed to remove bond vibrations (Hess et al., 1997). The particle mesh Ewald method (Darden et al., 1993; Essmann et al., 1995) coupled to periodic boundary conditions was used to simulate long-range electrostatic interactions using a direct space cutoff of 1.2 nm and a grid spacing of 0.15 nm. Because periodic boundary conditions can produce periodic artifacts, especially if combined with Ewald methods (Essmann et al., 1995), the minimum distance between protein molecules in adjacent boxes was calculated as a function of time. The simulation box size was adequate for the system, based on the cutoff used for electrostatic and short range contributions to the potential function. The van der Waals interactions were computed using periodic boundary conditions coupled to a spherical cutoff of 1.2 nm. The MDS was performed using GROMACS 4.6.1 (Hess et al., 2008). Additionally, a replicate MDS was carried out in order to confirm the results of the first MDS.
Docking and MDS simulations were combined by screening different SV2A conformers. LEV, BRIV, SEL and UCB-30889 were the selected racetam ligands. The criteria to include these ligands were as follows: LEV specifically binds to SV2A (Lynch et al., 2004) and is the lead anti-convulsive racetam. BRIV, SEL and UCB-30889 also bind to SV2A, but with higher affinity than LEV (Gillard et al., 2003, 2011; Matagne et al., 2009). BRIV and SEL have been the two leading candidates to replace LEV, since they were more potent and effective in animals model studies (Kaminski et al., 2012). Finally, UCB-30889 has been extensively used to characterize the binding properties of SV2A (Gillard et al., 2003; Lambeng et al., 2005; Shi et al., 2011). Interactions between the racetam ligands and nine SV2A conformers corresponding to every 2 ns during the last 16 ns of the MDS trajectory (when the system reached equilibrium) were simulated. In order to obtain the minimum energy conformation for each ligand structure, they were pre-optimized with AM1 (a semi-empirical method) and then optimized using the B3LYP/6-31G(d) basis set, both implemented in Gaussian 03 (Frisch et al., 1998). Dockings were performed using AutoDock 4.0.1 and AutoDock Tools 1.5.6 (Morris et al., 2008).
Before starting the docking simulations, hydrogen atoms were added to the polar atoms (considering a 7.4 pH value) and the Kollman charges were assigned for all atoms in the receptor. All rotating bonds, torsional degrees of freedom, atomic partial charges and non-polar hydrogens of the ligands were assigned. Then the ligands were docked onto the SV2A protein using a 70 Å × 70 Å × 70 Å grid box which covered all putative residues involved in racetam recognition according to the literature (Shi et al., 2011). A grid spacing of 0.375 Å was used under the hybrid Lamarckian Genetic Algorithm with an initial population of 100 randomly placed individuals and 1 × 107 as the maximum number of energy evaluations. All other docking parameters remained at their default settings. The resulting docked orientations within a root mean square deviation (RMSD) of 0.5 Å were clustered together. The lowest energy cluster for each compound, returned by AutoDock, was used for further analysis. The script and files used were prepared with AutoDock Tools 1.5.6 and this software was also used to visualize all the ligand–protein complexes together with VMD 1.9.1 (Humphrey et al., 1996). All the artwork was created using Gnuplot 4.6, VMD 1.8, Pymol 1.5 and Photoshop CS5 Extended 12.
After sending the SV2A protein sequence (FASTA format) to the I-Tasser server, it returned five possible 3D models with 12 α-helical TM segments, two long loops and a long cytoplasmic N-terminus. The C-score was between -2.84 to -3.02, the TM-score was 0.39 ± 0.13 and the RMSD was 15.4 ± 3.4 Å. The Ramachandran plot showed φ and ψ angles of 94.7–96.2 % in the favored-allowed regions of the SV2A models. In four of the models, the TM region was invaded importantly by the N-terminus. They were discarded for this reason and the model with an N-terminus separate from the TM region was selected (model 3; Figure 1). Moreover, the CPHmodels-3.2 and 3D-JIGSAWv3.0 servers were used to build 3D models of SV2A to confirm reproducibility. Both servers generated a structure quite similar to the I-Tasser model, predicting 12 TM regions.
FIGURE 1. 3D SV2A model. The I-Tasser server yielded a model of SV2A with 742 amino acids distributed in 12 TM domains consisting of α-helices, an undefined intra-luminal domain and two cytoplasmic domains: a relatively large N-terminus with a synaptotagmin 1 binding site (arrow) plus an ATP binding motif (arrow head) and a large loop between TM domains 6–7 with another ATP binding motif (two arrows heads).
Molecular Dynamics Simulations
The selected SV2A model was inserted into a POPC membrane and a MDS was performed to obtain a highly reliable 3D structure. Figure 2 shows the simulation system after structural minimization (Figure 2A) and after 50 ns of MDS (Figure 2B). Notice the rearrangement of the SV2A protein, which indicates that conformational changes occur in the extra and intra-membrane domains, being greatest in loop structure components. The MDS allowed SV2A to re-accommodate within the membrane region, producing a more reliable system to perform docking studies of known ligands.
FIGURE 2. Views of the simulation system, (A) after minimization and (B) after the MDS. Notice that SV2A rearranges within the lipid bilayer.
Figure 3 shows the RMSD values for α-carbons of the whole backbone (Figure 3A) or α-carbons of the TM backbone helical domains (Figure 3B). RMSD values typically help to visualize global structural changes to determine protein stability. The lower RMSD values (0.21 ± 0.04 nm) of the TM helical domains (Figure 3B), compared to the whole backbone (Figure 3A), indicate strong conformational changes in extra-membrane regions. The RMSD slope became modest at ∼20 ns and tended to zero beyond ∼30–35 ns. The RMSD values suggest that the system reaches equilibrium, as confirmed by the Rg values (Figure 4).
FIGURE 3. Root mean square deviation values for (A) the whole backbone and (B) the backbone of the TM region as a function of time. The gray bands correspond to the SD and the line within them represents the average value.
FIGURE 4. Radius of gyration values for (A) the whole backbone and (B) the backbone of the TM region as a function of time. The gray bands correspond to the SD and the line within them represents the average value.
The Rg of the whole protein backbone did not show a clear trend throughout the 50 ns trajectory (Figure 4A); it oscillated around an average value of 3.43 ± 0.02 nm. However, the variation was much lower when only the backbone atoms of the TM helical domain were considered (Figure 4B), although the Rg values for this region tended to increase slightly. This could be due to the fact that SV2A is an integral membrane protein with fairly stable TM domains and flexible extra-membrane portions.
The number of intra-molecular (Figure 5A), protein–water (Figure 5B) and protein–lipid (Figure 5C) H-bonds was determined as a function of time, based on geometrical criteria (donor–acceptor distance cutoff <3.5 Å and donor-H-acceptor angle <30°). Changes in inter-helical hydrogen bonding are associated with the conformational dynamics of membrane proteins (Bondar and White, 2012). Most H-bonds corresponded to protein–solvent interactions, while only a few were formed between SV2A and the lipid heads. There was significant variation in the number of intra-molecular H-bonds that may correspond to a conformational change; however, this number reached an equilibrium value during the last part of the trajectory.
FIGURE 5. Number of H-bonds as a function of time, (A) protein intra-molecular, (B) protein–water and (C) protein–lipid. The gray bands correspond to the SD and the line within them represents the average value.
To obtain more structural details about the behavior of the protein, the RMSF of the α-carbon positions throughout the whole trajectory were calculated. These values indicated that the most flexible parts of SV2A were the loops located outside the membrane; i.e., the N-terminus, the large cytoplasmic loop between the sixth-seventh TM regions and the large intra-luminal loop between the seventh-eighth TM domains (residues 476–594; Figure 6).
FIGURE 6. Graphic RMSF representations. (A) RMSF of the atomic positions and (B) representation of protein flexibility as a gradient from low to high fluctuations.
On the other hand, protein solvent contacts were studied using SAS, which is a measure of the protein surface exposed to the solvent. The total SAS and the contributions of hydrophobic or hydrophilic groups to the SAS were determined as a function of time (Figure 7A). No significant changes were observed along the trajectory. The SAS of the protein exposed to water decreased slightly (Figure 7B), whereas the SAS of the protein exposed to lipids increased a little (Figure 7C).
FIGURE 7. Evolution of the total, hydrophobic and hydrophilic protein SAS; (A) exposed to water plus lipids, (B) exposed to water or (C) exposed to lipids as a function of time.
Solvent accessible surface area by residues exposed to the membrane hydrophobic core helps to evidence the localization of several amino acids. Figure 8 shows the total area exposed to the membrane by residues. Some residues in the loops, N-terminus and ATP binding motif interacted with the membrane. The TM helical domains can be determined, since the highest SAS values correspond to residues totally exposed to the membrane hydrophobic core (Figure 8).
FIGURE 8. Solvent accessible surface area by residues exposed to the membrane inner core. The highest values correspond to amino acids totally exposed to membrane hydrophobic core.
In principal components analysis, the α-carbon matrix of SV2A shows a positive correlation between atomic motions of the large extra-vesicular loop and N-terminal segment (Figure 9A). It revealed that only 2–3 components are relevant for the trajectory. The projection of the trajectory onto the two principal eigenvectors shows a structure that reached equilibrium in the last portion of the trajectory; illustrated by denser yellow–red regions (Figure 9B).
FIGURE 9. PCA analysis. (A) Variance–covariance matrix of the SV2A α-carbon atoms for the 50 ns of the MDS. (B) Projection of the 50 ns trajectory onto the two principal eigenvectors obtained from the diagonalization of variance–covariance matrix; time along the trajectory is represented by a denser gradient.
The results of the replicate MDS were, in general, in agreement with the data obtained in the first MDS (data not shown). Although details varied, all MDS parameters behaved in a similar manner. In particular, the TM domains arrived at a conformational equilibrium; whereas the extra-membrane components remained mobile.
The values for the ΔG of binding between SV2A and the racetams varied between -5.26 and -8.82 Kcal/mol throughout the last 16 ns of the MDS trajectory (Table 1). The results showed the same affinity pattern in all evaluated conformations: UCB30889 > BRIV > SEL > LEV; which is in agreement with experimental SV2A binding data (Table 2; Noyer et al., 1995; Gillard et al., 2003, 2011; Bennett et al., 2007).
TABLE 1. Free energy of binding and residues involved in the interaction between the racetam derivatives and nine SV2A conformers during the MD trajectory (corresponding to every 2 ns during the last 16 ns), bold residues interacted most frequently.
TABLE 2. Mean ΔG of binding during the MD trajectory and experimental SV2A binding data (see text for references).
Ligand binding in dockings showed a clear preference for the following residues: T456, S665, W666, D670 and L689. These residues are within TM domains 7, 10 and 11. According to the SAS analysis these amino acids show values of zero or very close, which means that they are totally exposed to the membrane hydrophilic core (i.e., the pore or channel). Contacts between SV2A and the racetams were dominated by hydrophobic interactions and H-bonds. Every racetam showed a particular pattern of interactions when bound to SV2A; even if they shared some residues, the interaction was not the same in all cases. Figure 10 shows schematic representations of the interactions between the last SV2A conformer (50 ns) of the MDS and LEV (Figure 10A), BRIV (Figure 10B), SEL (Figure 10C) or UCB30889 (Figure 10D).
FIGURE 10. Putative racetam binding site within the last SV2A conformer of the MDS. 3D models and schematic representations of the interactions between SV2A and (A) LEV, (B) BRIV, (C) SEL or (D) UCB30889 are shown. Hydrophobic interactions and H-bonds participate importantly in ligand recognition within the binding site (amino acids in circles).
Trans-membrane proteins are difficult to crystallize; thus computational techniques such as protein folding using homology modeling and MDS have been very useful tools to obtain information about their structural properties. The I-Tasser server was used to generate a full-length 3D model of SV2A with adequate structural predictions and acceptable quantitative assessments (Zhang and Skolnick, 2004; Zhang, 2008; Nugent and Jones, 2012). Modeling the extra-membrane loops was much more challenging than the TM regions, primarily due to sparse contacts or poor predictions of them; additionally, inherent loop flexibility imposed extra difficulties. Nevertheless, I-Tasser modeled the N-terminus and the two long loops quite similar to the secondary structure predicted by UniProt (2014). It is clear, however, that more sophisticated strategies will be required to overcome a number of current limitations in this area (Nugent and Jones, 2012).
Although the obtained structure was acceptable, it was further refined by running a MDS. It showed the time dependent behavior of SV2A and provided detailed information on its fluctuations and conformational changes. The RMSD results showed that SV2A reached equilibrium in the last 15 ns of the MDS and this was in agreement with the Rg values. This was also true for the number of H-bonds, since they remained constant during the last part of the MDS. The results showed that the TM domains reached equilibrium and were likely to be correctly modeled; whereas the regions outside the membrane still remained mobile (also shown by the RMSF values). This is a general occurrence for TM proteins (Nugent and Jones, 2012) and it has been reported for other proteins, such as Bcl-2 (Ilizaliturri-Flores et al., 2013).
Extra-membrane domains, such as the long intra-luminal loops and N-terminus of SV2A, have post-translational modifications and interactions with other molecules, which might stabilize the protein conformation in natural conditions (Schivell et al., 1996, 2005; Nowack et al., 2010). The lack of post-translational changes and contacts with other molecules (which were not modeled due to their added complexity) may contribute to the extensive loop motion. Moreover, the extra-membrane regions of SV2A were evocative of intrinsically disordered domains within other proteins, such as Bcl-2 (Rautureau et al., 2010). All of the above, plus difficulties associated to proteins embedded in the lipid bilayer, may explain why SV2A has not been crystallized (Hipolito et al., 2014). The relevance of post-translational modifications on the conformational stability could be explored in future works with new, longer MD simulations. It is important to note that the racetam binding site is not part of the extra-membrane regions and is within the pore, a region in equilibrium for at least the last 15 ns of the MDS; thus, our model proves to be of acceptable quality.
Synaptic vesicle protein 2A underwent amino acids rearrangements, mostly in the loop structure components. Figure 2 and the SAS analysis showed an increase in the protein portion exposed to lipids; suggesting that amino acids initially exposed to the aqueous solvent were moving into the lipid bilayer. It is likely that they required a more hydrophobic environment or more favorable interactions with the lipid heads. In principal components analysis, the variance–covariance matrix of atomic motions showed which atoms moved together in a concerted way. It confirmed the conformational change that occurred in the large extra-vesicular loop and N-terminus, which were the most difficult regions to model and also the most flexible parts of the protein. The projection of the trajectory onto the two principal eigenvectors indicated that the protein experienced serious structural evolution during the first 15 ns of the MDS; however, its structure reached equilibrium in the last part of the trajectory.
We performed docking studies to explore the recognition pattern of racetams in SV2A. Docking studies can be performed on protein structures obtained from homology modeling, which maintains the backbone conformation (Tovchigrechko et al., 2002). However, the SV2A 3D structure obtained from I-Tasser was refined by performing a MDS; the advantage was that the side chain residues were energetically accommodated to their neighbors, resulting in several different conformers. Since the docking procedure lacks protein flexibility, we used the protein conformers from the MDS to obtain more reliable data, as has been published elsewhere (Sixto-Lopez et al., 2014). The MDS and docking simulations were thus combined, since using only one conformation may omit several structural properties and result in non-biological dockings (Ramirez-Duran et al., 2013). Our docking studies suggested a consensus binding site for racetam ligands within SV2A constituted by five residues: T456, S665, W666, D670 and L689. Tryptophan 666 directly engages the racetams, mostly by electrostatic interactions. The SV2A protein with a mutated W666 did not restore synaptic depression in neurons lacking SV2A and SV2B; showing that this residue is vital for SV2A neurotransmitter function (Nowack et al., 2010). Additionally, binding of UCB30889 to SV2A is altered when this amino acid is mutated, indicating that it may participate in ligand recognition (Shi et al., 2011).
The SV2A protein has been proposed to be a transporter owing to its high degree of homology with other mayor facilitator super-family (MFS) carriers (Bajjalieh et al., 1992; Feany et al., 1992; Gingrich et al., 1992) and its compact funnel-shaped structure with a visible indentation in the center indicative of a pore opening (Lynch et al., 2008). Tryptophan 666 is conserved in all SV2 isoforms and is homologous to a tryptophan in the 10th TM domain of MFS transporters which is vital for their activity (Nowack et al., 2010; Shi et al., 2011). Thus, it may provide the necessary hydrophobic milieu for transport of an endogenous substrate or LEV, in addition to its involvement in racetam binding. Residues Y461 and Y462 also frequently interacted with the racetams (Table 1); both are homologous to functional amino acids of MFS transporters (Shi et al., 2011). Additionally, Y461 is homologous to a residue important for chloride carrier activity (Bostick and Berkowitz, 2004).
It has been shown that 14 TM residues alter binding of UCB-30889 to SV2A when mutated (Shi et al., 2011). The present results coincide in the importance of W666 for ligand binding. The other four residues reported here were not mutated in the previous study. Y462 was also repeatedly found in this study; however, the 12 other amino acids described in the prior work were not regularly detected. These discrepancies may be due to the following: (1) SV2A was transiently expressed in COS-7 cells which are far removed from neurons, (2) the true functionality of SV2A mutants was unclear since they were not likely to exit the endoplasmic reticulum or Golgi (Nowack et al., 2010), (3) the residues that altered binding to UCB-30889 were not necessarily part of the ligand pocket; they may affect other aspects, such as changing the native conformation or stability of the protein. This, in turn, may change important aspects, such racetam access, retention or expulsion. These phenomena are not modeled by docking studies, which only predict ligand affinity for a region of the protein. A MDS of SV2A together with the racetams would be useful to further refine the interactions of SV2A and its racetam ligands.
The 3D model generated by I-Tasser and the structure obtained after the MDS show that the TM portions of SV2A may be accessed by the racetams. As previously suggested, LEV and related racetams likely enter into recycling vesicles, gain access to SV2A from its luminal side and bind in its TM segments (Meehan et al., 2011, 2012; Mendoza-Torreblanca et al., 2013). Once bound to SV2A, racetams inhibit its function; reducing the ready releasable vesicular pool and synaptic transmission, preferentially during long high frequency activity (Xu and Bajjalieh, 2001).
Modeling of the complete SV2A protein was challenging due to its large size (742 amino acids), its long N-terminus and its two long loops. However, a valid 3D model was obtained and the 50 ns of MDS were enough to reach structural equilibrium of the protein. The racetam binding site in SV2A was identified by docking studies employing different SV2A snapshots to explore the protein motion involved in ligand recognition; it showed that five residues participate in ligand recognition. Given the biological and pharmacological importance of SV2A, knowing its conformation is essential to direct further experimental work in order to elucidate its function and attempt the discovery of new ligands. Additionally, knowing the racetam binding site within SV2A should facilitate synthesizing suitable radio-ligands for PET studies with a high affinity and specificity. This, in turn, should allow the in vivo evaluation of treatment response and possibly disease progression.
All authors contributed to the written manuscript. All authors have given approval to the final version of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to MD.USe Innovations for their collaboration with the SV2A MDS, Sara Navarrete Hernández for her administrative assistance and Julio C. López Dávila for his support in the implementation and use of GROMACS. The research was supported by the following grants: Consejo Nacional de Ciencia y Tecnología 180069 and 132353; Instituto Nacional de Pediatría 111/2012; Instituto de Ciencia y Tecnología del Distrito Federal PIRIVE09-9; Secretaria de Investigación y Posgrado 20150322; Estimulo al Desempeño de los Investigadores and Comisión de Operación y Fomento de Actividades Académicas del Instituto Politécnico Nacional; Programa Iberoamericano de Ciencia y Tecnología para el Desarrollo.
AM1, Austin Model 1; BRIV, brivaracetam; CNS, central nervous system; H-bonds, hydrogen bonds; ΔG, free energy; LEV, levetiracetam; MFS, mayor facilitator family; MDS, molecular dynamics simulation; OPLS-AA, optimized potentials for liquid simulations all atoms; PET, positron emission tomography; POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine; Rg, radius of gyration; RMSD, root mean square deviation; RMSF, root mean square fluctuations; SAS, solvent accessible surface area; SD, standard deviation; SEL, seletracetam; SV2A, synaptic vesicle protein 2A; 3D, three-dimensional; TM, trans-membrane; Uniprot, universal protein resource; VMD, visual molecular dynamics.
Bajjalieh, S. M., Peterson, K., Shinghal, R., and Scheller, R. H. (1992). SV2, a brain synaptic vesicle protein homologous to bacterial transporters. Science 257, 1271–1273. doi: 10.1126/science.1519064
Bostick, D. L., and Berkowitz, M. L. (2004). Exterior site occupancy infers chloride-induced proton gating in a prokaryotic homolog of the ClC chloride channel. Biophys. J. 87, 1686–1696. doi: 10.1529/biophysj.104.042465
Chakrabarti, N., Neale, C., Payandeh, J., Pai, E. F., and Pomes, R. (2010). An iris-like mechanism of pore dilation in the CorA magnesium transport system. Biophys. J. 98, 784–792. doi: 10.1016/j.bpj.2009.11.009
Crowder, K. M., Gunther, J. M., Jones, T. A., Hale, B. D., Zhang, H. Z., Peterson, M. R., et al. (1999). Abnormal neurotransmission in mice lacking synaptic vesicle protein 2A (SV2A). Proc. Natl. Acad. Sci. U.S.A. 96, 15268–15273. doi: 10.1073/pnas.96.26.15268
de Groot, M., Aronica, E., Heimans, J. J., and Reijneveld, J. C. (2011). Synaptic vesicle protein 2A predicts response to levetiracetam in patients with glioma. Neurology 77, 532–539. doi: 10.1212/WNL.0b013e318228c110
Feng, G., Xiao, F., Lu, Y., Huang, Z., Yuan, J., Xiao, Z., et al. (2009). Down-regulation synaptic vesicle protein 2A in the anterior temporal neocortex of patients with intractable epilepsy. J. Mol. Neurosci. 39, 354–359. doi: 10.1007/s12031-009-9288-2
Frycia, A., Starck, J. P., Jadot, S., Lallemand, B., Leclercq, K., Lo Brutto, P., et al. (2010). Discovery of indolone acetamides as novel SV2A ligands with improved potency toward seizure suppression. ChemMedChem 5, 200–205. doi: 10.1002/cmdc.200900436
Gillard, M., Fuks, B., Leclercq, K., and Matagne, A. (2011). Binding characteristics of brivaracetam, a selective, high affinity SV2A ligand in rat, mouse and human brain: relationship to anti-convulsant properties. Eur. J. Pharmacol. 664, 36–44. doi: 10.1016/j.ejphar.2011.04.064
Gillard, M., Fuks, B., Michel, P., Vertongen, P., Massingham, R., and Chatelain, P. (2003). Binding characteristics of [3H]ucb 30889 to levetiracetam binding sites in rat brain. Eur. J. Pharmacol. 478, 1–9. doi: 10.1016/j.ejphar.2003.08.032
Gingrich, J. A., Andersen, P. H., Tiberi, M., El Mestikawy, S., Jorgensen, P. N., Fremeau, R. T., et al. (1992). Identification, characterization, and molecular cloning of a novel transporter-like protein localized to the central nervous system. FEBS Lett. 312, 115–122. doi: 10.1016/0014-5793(92)80917-6
Hess, B., Bekker, H. H. J., Berendsen, C., and Fraaije, J. G. E. M. (1997). LINCS: a linear constraint solver for molecular simulations. J. Comp. Chem. 18, 1463–1472. doi: 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H
Hess, B., Kutzner, C., Van Der Spoel, D., and Lindahl, E. (2008). GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theor. Comput. 4, 435–447. doi: 10.1021/ct700301q
Hipolito, C. J., Bashiruddin, N. K., and Suga, H. (2014). Protein cocrystallization molecules originating from in vitro selected macrocyclic peptides. Curr. Opin. Struct. Biol. 26, 24–31. doi: 10.1016/j.sbi.2014.03.001
Ilizaliturri-Flores, I., Correa-Basurto, J., Benitez-Cardoza, C. G., and Zamorano-Carrillo, A. (2013). A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations. J. Biomol. Struct. Dyn. 32, 1707–1719. doi: 10.1080/07391102.2013.833858
Janz, R., Goda, Y., Geppert, M., Missler, M., and Sudhof, T. C. (1999). SV2A and SV2B function as redundant Ca2+ regulators in neurotransmitter release. Neuron 24, 1003–1016. doi: 10.1016/S0896-6273(00)81046-6
Jorgensen, W. L., Maxwell, D. S., and Tirado-Rives, J. (1996). Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 118, 11225–11236. doi: 10.1021/ja9621760
Kaminski, G., Friesner, R. A., Tirado-Rives, J., and Jorgensen, W. L. (2001). Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides. J. Phys. Chem. B 105, 6474–6648. doi: 10.1021/jp003919d
Kaminski, R. M., Gillard, M., and Klitgaard, H. (2012). “Targeting SV2A for Discovery of Antiepileptic Drugs,” in Jasper’s Basic Mechanisms of the Epilepsies, 4th Edn, eds J. L. Noebels, M. Avoli, M. A. Rogawski, R. W. Olsen, and A. V. Delgado-Escueta (Bethesda, MD: National Center for Biotechnology Information). Available at: http://www.ncbi.nlm.nih.gov/books/NBK98183/
Lambeng, N., Gillard, M., Vertongen, P., Fuks, B., and Chatelain, P. (2005). Characterization of [(3)H]ucb 30889 binding to synaptic vesicle protein 2A in the rat spinal cord. Eur. J. Pharmacol. 520, 70–76. doi: 10.1016/j.ejphar.2005.07.029
Lindahl, E., and Edholm, O. (2000). Mesoscopic undulations and thickness fluctuations in lipid bilayers from molecular dynamics simulations. Biophys. J. 79, 426–433. doi: 10.1016/S0006-3495(00)76304-1
Lynch, B. A., Lambeng, N., Nocka, K., Kensel-Hammes, P., Bajjalieh, S. M., Matagne, A., et al. (2004). The synaptic vesicle protein SV2A is the binding site for the antiepileptic drug levetiracetam. Proc. Natl. Acad. Sci. U.S.A. 101, 9861–9866. doi: 10.1073/pnas.0308208101
Lynch, B. A., Matagne, A., Brannstrom, A., Von Euler, A., Jansson, M., Hauzenberger, E., et al. (2008). Visualization of SV2A conformations in situ by the use of Protein Tomography. Biochem. Biophys. Res. Commun. 375, 491–495. doi: 10.1016/j.bbrc.2008.07.145
Matagne, A., Margineanu, D. G., Potschka, H., Loscher, W., Michel, P., Kenda, B., et al. (2009). Profile of the new pyrrolidone derivative seletracetam (ucb 44212) in animal models of epilepsy. Eur. J. Pharmacol. 614, 30–37. doi: 10.1016/j.ejphar.2009.04.024
Meehan, A. L., Yang, X., Mcadams, B. D., Yuan, L., and Rothman, S. M. (2011). A new mechanism for antiepileptic drug action: vesicular entry may mediate the effects of levetiracetam. J. Neurophysiol. 106, 1227–1239. doi: 10.1152/jn.00279.2011
Mendoza-Torreblanca, J. G., Vanoye-Carlo, A., Phillips-Farfan, B. V., Carmona-Aparicio, L., and Gomez-Lira, G. (2013). Synaptic vesicle protein 2A: basic facts and role in synaptic function. Eur. J. Neurosci. 38, 3529–3539. doi: 10.1111/ejn.12360
Nowack, A., Yao, J., Custer, K. L., and Bajjalieh, S. M. (2010). SV2 regulates neurotransmitter release via multiple mechanisms. Am. J. Physiol. Cell Physiol. 299, C960–C967. doi: 10.1152/ajpcell.00259.2010
Noyer, M., Gillard, M., Matagne, A., Henichart, J. P., and Wulfert, E. (1995). The novel antiepileptic drug levetiracetam (ucb L059) appears to act via a specific binding site in CNS membranes. Eur. J. Pharmacol. 286, 137–146. doi: 10.1016/0014-2999(95)00436-O
Nugent, T., and Jones, D. T. (2012). Accurate de novo structure prediction of large transmembrane protein domains using fragment-assembly and correlated mutation analysis. Proc. Natl. Acad. Sci. U.S.A. 109, E1540–E1547. doi: 10.1073/pnas.1120036109
Ramirez-Duran, L. A., Rosales-Hernandez, M. C., Hernandez-Rodriguez, M., Mendieta-Wejebe, J. E., Trujillo-Ferrara, J., and Correa-Basurto, J. (2013). Mapping myeloperoxidase to identify its promiscuity properties using docking and molecular dynamics simulations. Curr. Pharm. Des. 19, 2204–2215. doi: 10.2174/1381612811319120008
Reigada, D., Diez-Perez, I., Gorostiza, P., Verdaguer, A., Gomez De Aranda, I., Pineda, O., et al. (2003). Control of neurotransmitter release by an internal gel matrix in synaptic vesicles. Proc. Natl. Acad. Sci. U.S.A. 100, 3485–3490. doi: 10.1073/pnas.0336914100
Schivell, A. E., Batchelor, R. H., and Bajjalieh, S. M. (1996). Isoform-specific, calcium-regulated interaction of the synaptic vesicle proteins SV2 and synaptotagmin. J. Biol. Chem. 271, 27770–27775. doi: 10.1074/jbc.271.44.27770
Schivell, A. E., Mochida, S., Kensel-Hammes, P., Custer, K. L., and Bajjalieh, S. M. (2005). SV2A and SV2C contain a unique synaptotagmin-binding site. Mol. Cell. Neurosci. 29, 56–64. doi: 10.1016/j.mcn.2004.12.011
Shi, J., Anderson, D., Lynch, B. A., Castaigne, J. G., Foerch, P., and Lebon, F. (2011). Combining modelling and mutagenesis studies of synaptic vesicle protein 2A to identify a series of residues involved in racetam binding. Biochem. Soc. Trans. 39, 1341–1347. doi: 10.1042/BST0391341
Sixto-Lopez, Y., Gomez-Vidal, J. A., and Correa-Basurto, J. (2014). Exploring the potential binding sites of some known HDAC inhibitors on some HDAC8 conformers by docking studies. Appl. Biochem. Biotechnol. 173, 1907–1926. doi: 10.1007/s12010-014-0976-1
Surges, R., Volynski, K. E., and Walker, M. C. (2008). Is levetiracetam different from other antiepileptic drugs? Levetiracetam and its cellular mechanism of action in epilepsy revisited. Ther. Adv. Neurol. Disord. 1, 13–24. doi: 10.1177/1756285608094212
Tieleman, D. P., Maccallum, J. L., Ash, W. L., Kandt, C., Xu, Z., and Monticelli, L. (2006). Membrane protein simulations with a united-atom lipid and all-atom protein model: lipid-protein interactions, side chain transfer free energies and model proteins. J. Phys. Condens. Matter. 18, S1221–S1234. doi: 10.1088/0953-8984/18/28/S07
Toering, S. T., Boer, K., De Groot, M., Troost, D., Heimans, J. J., Spliet, W. G., et al. (2009). Expression patterns of synaptic vesicle protein 2A in focal cortical dysplasia and TSC-cortical tubers. Epilepsia 50, 1409–1418. doi: 10.1111/j.1528-1167.2008.01955.x
UniProt, C. (2014). Uniprot Consortium/SV2A/Q02563. Available at: http://www.uniprot.org.
van Vliet, E. A., Aronica, E., Redeker, S., Boer, K., and Gorter, J. A. (2009). Decreased expression of synaptic vesicle protein 2A, the binding site for levetiracetam, during epileptogenesis and chronic epilepsy. Epilepsia 50, 422–433. doi: 10.1111/j.1528-1167.2008.01727
Wan, Q. F., Zhou, Z. Y., Thakur, P., Vila, A., Sherry, D. M., Janz, R., et al. (2010). SV2 acts via presynaptic calcium to regulate neurotransmitter release. Neuron 66, 884–895. doi: 10.1016/j.neuron.2010.05.010
WHO. (2014). World Health Organization. Available at: http://www.who.int/mediacentre/factsheets/fs999/en/
Yao, J., Nowack, A., Kensel-Hammes, P., Gardner, R. G., and Bajjalieh, S. M. (2010). Cotrafficking of SV2 and synaptotagmin at the synapse. J. Neurosci. 30, 5569–5578. doi: 10.1523/JNEUROSCI.4781-09.2010
Keywords: SV2A, epilepsy, levetiracetam, brivaracetam, seletracetam
Citation: Correa-Basurto J, Cuevas-Hernández RI, Phillips-Farfán BV, Martínez-Archundia M, Romo-Mancillas A, Ramírez-Salinas GL, Pérez-González OA, Trujillo-Ferrara J and Mendoza-Torreblanca JG (2015) Identification of the antiepileptic racetam binding site in the synaptic vesicle protein 2A by molecular dynamics and docking simulations. Front. Cell. Neurosci. 9:125. doi: 10.3389/fncel.2015.00125
Received: 06 November 2014; Accepted: 17 March 2015;
Published online: 10 April 2015
Edited by:Victoria Campos-Peña, Instituto Nacional de Neurologia y Neurocirugia, Mexico
Reviewed by:Benjamín Florán, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Mexico
Alan Talevi, National University of La Plata, Argentina
Copyright © 2015 Correa-Basurto, Cuevas-Hernández, Phillips-Farfán, Martínez-Archundia, Romo-Mancillas, Ramírez-Salinas, Pérez-González, Trujillo-Ferrara and Mendoza-Torreblanca. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Julieta G. Mendoza-Torreblanca, Laboratorio de Neuroquímica, Instituto Nacional de Pediatría, Insurgentes Sur No. 3700-C, Colonia Insurgentes Cuicuilco, Delegación Coyoacán, México City 04530, Mexico firstname.lastname@example.org