A Synthetic Combinatorial Strategy for Developing α-Conotoxin Analogs as Potent α7 Nicotinic Acetylcholine Receptor Antagonists*

α-Conotoxins are peptide neurotoxins isolated from venomous cone snails that display exquisite selectivity for different subtypes of nicotinic acetylcholine receptors (nAChR). They are valuable research tools that have profound implications in the discovery of new drugs for a myriad of neuropharmacological conditions. They are characterized by a conserved two-disulfide bond framework, which gives rise to two intervening loops of extensively mutated amino acids that determine their selectivity for different nAChR subtypes. We have used a multistep synthetic combinatorial approach using α-conotoxin ImI to develop potent and selective α7 nAChR antagonists. A positional scan synthetic combinatorial library was constructed based on the three residues of the n-loop of α-conotoxin ImI to give a total of 10,648 possible combinations that were screened for functional activity in an α7 nAChR Fluo-4/Ca2+ assay, allowing amino acids that confer antagonistic activity for this receptor to be identified. A second series of individual α-conotoxin analogs based on the combinations of defined active amino acid residues from positional scan synthetic combinatorial library screening data were synthesized. Several analogs exhibited significantly improved antagonist activity for the α7 nAChR compared with WT-ImI. Binding interactions between the analogs and the α7 nAChR were explored using a homology model of the amino-terminal domain based on a crystal structure of an acetylcholine-binding protein. Finally, a third series of refined analogs was synthesized based on modeling studies, which led to several analogs with refined pharmacological properties. Of the 96 individual α-conotoxin analogs synthesized, three displayed ≥10-fold increases in antagonist potency compared with WT-ImI.

␣-Conotoxins are peptide neurotoxins isolated from venomous cone snails that display exquisite selectivity for different subtypes of nicotinic acetylcholine receptors (nAChR). They are valuable research tools that have profound implications in the discovery of new drugs for a myriad of neuropharmacological conditions. They are characterized by a conserved two-disulfide bond framework, which gives rise to two intervening loops of extensively mutated amino acids that determine their selectivity for different nAChR subtypes. We have used a multistep synthetic combinatorial approach using ␣-conotoxin ImI to develop potent and selective ␣ 7 nAChR antagonists. A positional scan synthetic combinatorial library was constructed based on the three residues of the n-loop of ␣-conotoxin ImI to give a total of 10,648 possible combinations that were screened for functional activity in an ␣ 7 nAChR Fluo-4/Ca 2؉ assay, allowing amino acids that confer antagonistic activity for this receptor to be identified. A second series of individual ␣-conotoxin analogs based on the combinations of defined active amino acid residues from positional scan synthetic combinatorial library screening data were synthesized. Several analogs exhibited significantly improved antagonist activity for the ␣ 7 nAChR compared with WT-ImI. Binding interactions between the analogs and the ␣ 7 nAChR were explored using a homology model of the amino-terminal domain based on a crystal structure of an acetylcholine-binding protein. Finally, a third series of refined analogs was synthesized based on modeling studies, which led to several analogs with refined pharmacological properties. Of the 96 individual ␣-conotoxin analogs synthesized, three displayed >10-fold increases in antagonist potency compared with WT-ImI.
Venomous marine cone snails continue to provide researchers with a valuable source of ligands able to pharmacologically distinguish between various classes of ion channels and receptors (1). Their venom is used to swiftly immobilize their fast moving prey, a result of the singular and combined actions of between 50 and 200 individual peptide neurotoxin components known as conotoxins. Conotoxins are gene-encoded products that are expressed in venom glands as natural combinatorial peptide libraries (2). It is estimated that across the 500 species that make up entire genus of cone snails, there are more than 100,000 individual conotoxins with unique neuropharmacological properties (3).
Nicotinic acetylcholine receptors (nAChRs) 4 are ligandgated ion channels that mediate the fast excitatory cholinergenic neurotransmission in the central and peripheral nervous systems (4). They are involved in cognitive function, pain pathways, and nicotine reward mechanisms (5). They are pentameric complexes composed of five closely related subunits arranged around a central cation conducting pore. The subunits are divided into the ␣ 1 , ␤ 1 , ␥, and ␦/⑀ subunits forming the muscle-type nAChRs and the ␣ 2 -␣ 10 and ␤ 2 -␤ 4 subunits forming neuronal nAChRs (6). Each nAChR subtype is characterized by distinct pharmacological properties and is involved in a range of physiological and pathophysiological functions. For example, ␣ 7 nAChRs are believed to play an important role in learning and memory; hence they have been proposed to be involved in diseases such as schizophrenia and Alzheimer disease (5). The large number of physiologically relevant nAChR subtypes and the numerous diverse physiological functions regulated by them underlines the importance of developing ligands that can specifically target different subtypes of these receptors.
␣-Conotoxins are selective antagonists for different subtypes of nAChRs (7). They typically contain 12-20 amino acid residues and consist of two highly conserved disulfide bonds connected via a [1][2][3], [2][3][4] (globular) arrangement, giving rise to a two-loop cysteine framework, denoted m and n, respectively ( Fig. 1). With the exception of the conserved cysteines, each position in the ␣-conotoxin sequence is extensively mutated. The variability of amino acids and functionalities contained within the loops has been shown to determine the potency and subtype selectivity profile for the nAChRs (8 -16). The published x-ray co-crystal structures of several ␣-conotoxins in complex with acetylcholine-binding proteins (AChBPs) have been reported, providing structural insight into ␣-conotoxin/ nAChR binding at the molecular level (17)(18)(19)(20). AChBPs are a class of water-soluble proteins that display significant sequence homology with the ligand-binding domain of ␣ 7 nAChRs. These co-crystal structures have been used to yield several new ␣-conotoxin analogs characterized by even more potent antagonist activities and different selectivity profiles than the native toxins (14).
The ␣4/3-conotoxins are a rich source of functionally diverse inhibitors of all ␣-subunit-containing nAChRs (21). Several examples, including ␣-conotoxins ImI, ImII, and RgIA, that target the ␣ 7 and ␣ 9 ␣ 10 nAChRs are well characterized, and more are continually being discovered (22). Interestingly, ImI and RgIA differ by only three residues and share an SDPR m-loop sequence, with variability occurring in the n-loop (Fig. 1). However, ImI and RgIA are selective for the ␣ 7 and ␣ 9 ␣ 10 nAChRs, respectively. Furthermore, the SDPR m-loop motif is also found in the ␣4/7-conotoxins EpI and Vc1.1, and it has been suggested that the ␣4/3-conotoxin n-loop residues may be important for refining nAChR selectivity (21).
Given that conotoxins are traditionally discovered from cone snail venom extracts that are themselves the natural combinatorial mixtures of bioactive peptides, it occurred to us that synthetic mixture-based combinatorial libraries of conotoxin would be of value in structure-activity relationship studies. With their highly conserved cysteine frameworks and intervening loops of hypermutated amino acids that gives rise to rigid three-dimensional conformations, conotoxins represent attractive combinatorial templates. Positional scanning synthetic combinatorial libraries (PS-SCL) are mixture-based libraries that provide a rapid means to acquire information regarding all possible variable positions within a chemical framework (23)(24)(25). They contain both defined and mixture positions of diversity, enabling the most active functionalities for each position to be identified directly from initial screening data, allowing the design of individual compounds (24, 26 -28).
In this study, we have developed novel ligands for the ␣ 7 nAChR using a multistep synthetic combinatorial approach through modification of the three residues forming the n-loop of ␣-conotoxin ImI. We have also employed homology modeling and docking studies using the ImI-AChBP complex to rationalize observed differences in antagonist potency and to refine the pharmacological properties of these analogs.
Crude samples were individually oxidized by gently agitating in a solution of 0.1 M ammonium bicarbonate, 50% isopropyl alcohol, pH 8.2, for 2 days at room temperature in an open vessel. PS-SCL mixtures were injected onto a preparative C 18 RP-HPLC column and eluted from isocratic 62% acetonitrile followed by lyophilization to remove associated buffer salts and residual scavengers. Individual analogs were eluted from a linear gradient of water/acetonitrile (2-62% for 30 min) containing 0.1% trifluoroacetic acid and fractionated accordingly. Regioselective synthesis of compound 36 (globular and ribbon isomers) was achieved by stirring the reduced/S-acetomidomethyl-protected peptide in 0.1 M NH 4 HCO 3 for 2 days, and isolating by preparative RP-HPLC. The partially oxidized peptides (ϳ10 mg) were dissolved in 80% methanol (25 ml), and 0.1 M HCl was added (1 ml), followed by 0.1 M I 2 in methanol (2.5 ml). The solution was stirred for 10 min before quenching with 0.1 M sodium thiosulfate. The fully oxidized peptides were finally isolated using preparative RP-HPLC.
NMR Spectroscopy-NMR data were recorded on a sample dissolved in 90% H 2 O, 10% D 2 O at pH 3.5 using a Bruker ARX 600 MHz spectrometer. Two-dimensional total correlation spectroscopy and nuclear Overhauser effect spectroscopy NMR were acquired for all compounds, and additional DQF-COSY and E-COSY experiments were recorded for compound 36. For compound 36, peak intensities from a nuclear Overhauser effect spectroscopy spectrum with a mixing time of 200 ms at 285 K were used to determine distance restraints. Backbone dihedral restraints were determined from 3 J HN-H␣ coupling constants, and angles were restrained to Ϫ60 Ϯ 30°for 3 J HN-H␣ Ͻ 5.8 Hz (Cys 2 , Cys 3 , Asp 5 , Tyr 10 , and Cys 12 ), Ϫ120 Ϯ 30°for 3 J HN-H␣ Ͼ 8.0 Hz (Abu11), and Ϫ100 Ϯ 80°when H ␣iϪ1 Ϫ HN i is significantly stronger than H ␣iϪ1 Ϫ HN i (Ser 4 and Leu 9 ). Intra-residue nuclear Overhauser effect and 3 J HN-H␤ coupling patterns were used in assigning 1 angle restraints of some side chains (Cys 2 , Cys 3 , Ser 4 , Arg 7 , and Cys 8 to 180 Ϯ 30°and Cys 12 to Ϫ60 Ϯ 30°).
Initial structures were generated using Cyana software (32), and the final structure calculations were performed using a simulated annealing protocol with CNS (33) as described previously (34). Fifty structures were calculated, and the 20 structures with the lowest overall energies were retained for analysis. Structures were visualized using MOLMOL (35) and analyzed with PROMOTIF (36) and PROCHECK_NMR (37).
Pharmacology-␣-Conotoxin analogs were functionally characterized at the ␣ 7 nAChR in the fluorescence-based Fluo-4/Ca 2ϩ assay using a stable ␣ 7 -GH3 cell line (38) as described previously (14). Cell lines were split into poly-D-lysine-coated black 96-well clear bottom plates (BD Biosciences), following a 64 -72-h incubation. On the day of the assay, the culture medium was aspirated, and the cells were incubated in 50 l of assay buffer (Hanks' buffered saline solution containing 20 mM HEPES, 1 mM CaCl 2 , 1 mM MgCl 2 , and 2.5 mM probenecid, pH 7.4) supplemented with 6 mM Fluo-4/AM (Molecular Probes, Eugene, OR) at 37°C for 1 h. The buffer was aspirated, and the cells were washed once with 100 l of assay buffer, followed by addition of 100 l of assay buffer supplemented with 100 M genistein and various concentrations of the ␣-conotoxin analogs. Following 30 min of incubation at 37°C in a 5% CO 2 incubator, the 96-well plate was assayed in a NOVOstar TM microplate reader (BMG Labtechnologies, Offenburg, Germany) measuring emission at 520 nm caused by excitation at 485 nm before and up to 60 s after addition of 33 l of agonist solution in assay buffer. The compounds were characterized in duplicate at least three times using EC 80 -EC 90 concentrations of acetylcholine as agonist.
Selected ␣-conotoxin analogs were also functionally characterized at a HEK293 cell line stably expressing the rat ␣ 3 ␤ 4 nAChR and at a HEK293T cell line stably expressing the mouse ␣ 4 ␤ 2 nAChR (39,40). The characterization at the ␣ 3 ␤ 4 -HEK293 cell line was performed in the Fluo-4/Ca 2ϩ assay. This assay was performed essentially as described for the ␣ 7 -GH3 cell line, except that genistein was not included in the assay buffer and EC 80 -EC 90 concentrations of epibatidine, not acetylcholine, were used as agonist. The characterization at the ␣ 4 ␤ 2 -HEK293 cell line was performed in the Fluorescent Imaging Plate Reader Membrane Potential Blue TM assay (Molecular Devices) as described previously, using EC 80 -EC 90 concentrations of epibatidine as agonist (14).
Homology Modeling and Molecular Dynamics-A homology model of the extracellular amino-terminal domain of human nAChR ␣ 7 was built using the Aplysia AChBP x-ray crystal structure bound with ImI (Protein Data Bank code 2BYP (19)). Sequence alignment of nAChR and AChBP was performed through ClustalW (41) and Prime; a protein structure prediction suite from Schrödinger LLC was used for model building (42). The aligned sequence of query and template showed 26% homology and 67% similarity. A homology structure was built using this aligned template taking into account the effects of solvent and the bound conotoxin through the various algorithms implemented in Prime. Missing query residues that did not match or align well with the template sequence were built using an ab initio procedure (43). Water molecules were retained from the template, and all the steric clashes were refined through minimization using OPLS2005 force field. The final model showed a backbone root mean square deviation of 0.49 Å with the template structure. Using this model, two separate docking models were generated using high and low affinity binding analogs from the second series of compounds. In both of these models, the pentameric nAChR structure was retained, superimposed from the template structure, and mutated according to the second iteration library ␣ 7 nAChR functional assays. In model 1, ImI was mutated to Leu-Trp-Abu (compound 28, high antagonist potency analog), and in model 2 ImI was mutated to Ala-Phe-Arg (compound 5, low antagonist potency analog).
To further refine the docked complexes, short 500-ps Molecular Dynamics simulations were performed using the NAMD2 program (44) with the CHARMM 29b2 force field (45). The initial docked and refined structures from the previous step were taken as starting points for Molecular Dynamics simulations. The system setup procedure was initiated by adding hydrogen atoms and a box of TIP3 water molecules (solvation), such that there was at minimum 13.0 Å of water between the surface of the protein and the edge of the simulation box using the Solvate plug-in of the Visual Molecular Dynamics program (46). Any added bulk water molecules within 2.5 Å of the protein were excluded. To maintain the electrical neutrality of the system, appropriate numbers of ions (15 Na ϩ ) were added using the Autoionize plug-in in Visual Molecular Dynamics programs, which were initially set at least 7.0 Å away from the surface of the protein.

Construction of PS-SCL
An initial PS-SCL based on the three residues in the ␣-conotoxin ImI n-loop was constructed to identify amino acid substitutions in this loop that confer activity for the ␣ 7 nAChR. Three sub-libraries were prepared, where O n is a single defined position, and X is an equimolar mixture of 22 natural and nonnatural L-amino acids (see Table 1). A total of 66 mixtures, each containing 484 compounds, were obtained in three sub-libraries consisting of a total of 10,648 possible individual conotoxins. Cysteine was omitted from the O n and X-positions to avoid interfering with the conserved disulfide bond framework during oxidative folding. Methionine was also omitted from the library to avoid complications associated with the formation of its corresponding methionine sulfoxide. As such, isosteric Abu and norleucine were included on the O n -and X-positions, together with Nva to complete the series of side chains containing hydrophobic alkyl groups. Additionally, hydroxyproline, a commonly occurring post-translational modification found in several conotoxins, was also incorporated into this library.
The PS-SCL was assembled using manual solid-phase peptide synthesis with Boc chemistry using polypropylene tea bags (29). X-positions were coupled using the portioning mixing strategy, which greatly simplifies library assembly compared with other approaches such as divide-couple-recombine or split syntheses (25). However, competitive coupling rates can be predominantly influenced by the relative nucleophilicity and steric hindrance of the incoming activated amino acid. Therefore, predetermined ratios were used by adjusting the relative concentration of each amino acid in the X-position (31).
Following assembly of the PS-SCL, each sample was cleaved using a two-step low-high HF procedure. Each sample was oxidized as a crude mixture in aqueous buffer consisting of 0.1 M NH 4 HCO 3 , pH 8.2, containing 50% isopropyl alcohol, which has been reported as the optimum condition for generating the native disulfide bond isomer of WT-ImI (47), and subsequently isolated from the oxidation buffer by elution with 62% isocratic acetonitrile on preparative RP-HPLC.

Functional Characterization of Mixture-based Libraries
The PS-SCL mixtures were screened for inhibitory ␣ 7 nAChR signaling in ␣ 7 -GH3 cells using the Fluo-4/Ca 2ϩ assay (Table 1 and Fig. 2) (38). We and others have previously found that that the ␣ 7 nAChR in this and other fluorescence-based functional assays displays functional properties that are in agreement with reported electrophysiological recordings (38,48).
Screening and deconvolution of the PS-SCL library revealed information regarding amino acid residues in positions O 9 , O 10 , and O 11 in WT-ImI for optimal binding to the ␣ 7 nAChR. In the case of the O 9 -position, the most active amino acid was defined as Nva, since this mixture exhibited a similar IC 50 value to WT-ImI (IC 50 WT-ImI /IC 50 PS-SCL ratio 0.82). Introduction of Leu and Ile residues in position O 9 also resulted in mixtures with high antagonist potencies, displaying IC 50 WT-ImI /IC 50 PS-SCL ratios of 0.23 and 0.21, respectively. Notably, each of these three mixtures exhibited significantly higher antagonist potencies than the mixtures containing the Ala residue also present in this position in WT-ImI (IC 50 WT-ImI /IC 50 PS-SCL ratio 0.10). In contrast, the most active amino acid at the O 10 -position was the Trp residue also found in WT-ImI (IC 50 WT-ImI /IC 50 PS-SCL ratio 0.43). Other residues introduced at the O 10 -position, which gave rise to mixtures with high antagonist potencies, were the other aromatic side chain bearing amino acids Phe and Tyr (IC 50 WT-ImI /IC 50 PS-SCL ratios 0.25 and 0.23, respectively). Finally, in the O 11 -position, the PS-SCL mixture with the Arg residue, also found in this position in WT-ImI, exhibited weaker antagonist activity (IC 50 WT-ImI /IC 50 PS-SCL ratio 0.10) compared with mixtures with His (ratio 0.30) and Trp (ratio 0.25), followed by mixtures with Abu and Tyr (both ratios of 0.16). To this end, deconvolution of the PS-SCL, including native WT-ImI residues, involved choosing Ala, Nva, and Leu residues in the O 9 -position, Trp, Phe, and Tyr residues in the O 10 -position, and Abu, Arg, His, and Trp residues in the O 11 -position.

Construction and Functional Characterization of a Second Iteration Library of Individual Analogs
Combinations of the residues giving rise to the most potent PS-SCL mixtures were used to construct a second series library consisting of a total of 36 individual analogs (see Table 2). Although oxidation using 0.1 M NH 4 HCO 3 , 50% isopropyl alco- hol yielded one major isomer as determined by LC/MS, several analogs were accompanied by a second minor peak with identical mass (ϳ10% of the total yield of oxidized conotoxin). This was presumably the ribbon isomer (i.e. [1][2][3][4], [2][3] disulfide bond connectivity) and in all cases was readily separated using preparative RP-HPLC. The individual analogs were tested for activity at the ␣ 7 nAChR in the Fluo-4/Ca 2ϩ assay ( Table 2 and Fig. 3a). Individual analogs are defined in this study as being very high (Ͼ10fold increase compared with that of WT-ImI), high (6 -10-fold increase), medium (0 -6-fold increase), low (10-fold decrease) antagonist potency, or inactive (Ͼ10-fold decrease).
For single amino acid substitutions in WT-ImI, nearly all analogs containing an Ala residue at O 9 displayed low antagonist potencies at the ␣ 7 nAChR (Table 2 and Fig. 3a), with compounds 5, 8, 11, and 12 displaying no significant activity. Compounds 2 and 10 were the only analogs containing the native Ala at position O 9 that exhibited similar antagonist potency to WT-ImI. Although a considerable number of the compounds containing Nva and Leu at O 9 clearly showed increased activity, the analogs with single substitutions of Ala for Nva or Leu in WT-ImI (compounds 13 and 25) both displayed IC 50 values comparable with WT-ImI. The analog containing single amino acid substitutions at O 10 with Tyr (compound 9) exhibited a 5-fold lower antagonist potency compared with WT-ImI, whereas substitution with Phe (compound 5) was inactive. At O 11 , substitution of Arg 11 for His (compound 2) was the only single amino acid substitution in WT-ImI with similar antagonist potency as the wild type peptide, as substitutions with Trp (compound 3) or Abu (compound 4) at this position yielded analogs with significantly lower antagonist potencies.
The analogs that displayed the most profound increases in antagonist potencies contained multiple substitutions of amino acids in the n-loop. The majority of the analogs containing Nva or Leu at O 9 displayed significantly increased antagonist potencies relative to WT-ImI. Furthermore, the most potent analogs containing Leu at O 9 also contained Abu at O 11 (compounds 28, 32, and 36). However, most compounds containing Nva at O 9 displayed medium antagonist potency compared with Leu-containing analogs. Although it seems to be of little importance whether the O 10 -position in the analogs was a Trp or a Tyr residue for their activities at the ␣ 7 nAChR, analogs containing these two amino acids demonstrated increased antagonist potencies compared with analogs containing Phe in this position. Compounds 14,16,22,28,32, and 36 all demonstrated notably higher potency than WT-ImI, with compound 36 exhibiting 10 times more antagonist potency than the native ␣-conotoxin. Importantly, compound 14 (Nva 9 -Trp 10 -His 11 ), which was predicted to be the most active analog based on the initial PS-SCL screening data, was included in the group of analogs exhibiting high antagonist potency.

Structural Characterization of Selected ImI Analogs
Structural variations resulting from hypermutation of the ImI n-loop for a selection of the most potent analogs were examined by NMR spectroscopy. NMR spectra were obtained for compounds 14, 16, 22, 24, 28, 32, and 36. NMR data were recorded at 285 K and showed single spin systems for each residue for all peptides. The spectra were assigned, and chemical shift analysis indicated that the H␣ shift values of all peptides were comparable with those of the published values for WT-ImI (49), with only small local variations around the sites of mutation (Fig. 4a). Therefore, overall the chemical shifts of the analogs examined by NMR indicated the peptides were similar in structure to the globular isomer of WT-ImI.
The disulfide bond connectivity of compound 36 was unambiguously confirmed using a regioselective disulfide bond forming strategy. Both the globular and ribbon isomers were synthesized utilizing the orthogonal S-Acm and S-4-methylbenzyl combination in a two-step oxidation procedure. Co-injection of both isomers with the sample obtained from random oxidation confirmed its identity as the globular disulfide bond isomer.
The full three-dimensional structure of compound 36 was determined. Sufficient nuclear Overhauser effects were observed to calculate a three-dimensional structure using a simulated annealing protocol in CNS. A set of 50 structures was calculated with 51 distance restraints consisting of 38 sequential, 9 medium range, and 4 long range restraints. Eight angle restraints were included and 6 1 angle restraints. The 20 lowest energy structures were chosen to represent the structure of compound 36, and structural statistics for this NMR ensemble are provided in supplemental Table 1, and the structural data have been deposited in the Biological Magnetic Resonance Bank. A well defined structure was obtained as shown by the root mean square deviation over residues 1-12 of 0.48 Ϯ 0.28 Å for the backbone atoms and 0.99 Ϯ 0.35 Å for heavy atoms. Analysis of the structures with PROMOTIF (36) identified a 3 10 helix between residues 2 and 4, which is illustrated by the ribbon depiction in Fig. 4b, and a type I turn from residues 5 to 8.
A comparison of the structure of compound 36 with native ImI (Protein Data Bank code 1IM1) (51), shown in Fig. 4c, shows the high degree of similarity between the peptide backbones of the two molecules. The root mean square deviation of the backbone atoms of the mean structures from the two ensembles is only 0.82, although most of this difference is around the substituted residues in the n-loop. Furthermore, a comparison of the NMR solution structure of compound 36 with that of the known ImI/AChBP crystal structure showed that overall the backbone architecture is very similar with backbone root mean square deviation of only ϳ1 Å (18). The 21 conformations of NMR structure showed difference only in the dynamics of side chains of residues Arg 7 and Tyr 10 .

Homology Modeling and Docking Studies
In an effort to further refine the activity and to gain additional insights into the binding interactions of the ␣-conotoxin ImI analogs with the ␣ 7 nAChR, computational modeling and simulations were used for selected analogs. A homology model of the amino-terminal binding domain of ␣ 7 nAChR was built using the AChBP/ImI x-ray crystal structure (Protein Data Bank code 2BYP) as a template (18). Two docking models were obtained using high and low antagonist potency analogs from the second iteration library. The pentameric nAChR structure was retained, superimposed from the template structure, and for each model the co-crystallized ␣-conotoxin was mutated according to the second iteration library ␣ 7 nAChR functional assays. In Model 1, ImI was mutated to Leu 9 -Trp 10 -Abu 11 (compound 28, high antagonist potency), and in Model 2 ImI was mutated to Ala 9 -Phe 10 -Arg 11 (compound 5, inactive). Further refinement and analysis of both models were made on the basis of 500 ps of molecular dynamics simulations. Although a complete computation-based explanation of the structure-activity relationship studies of the 36 compounds in the second series of individual compounds is beyond the scope of this paper, the focused analysis on the final models of the three mutated residues revealed the following trends that were valuable in designing a third iteration library (Fig. 5).

TABLE 2 Functional properties of the second series of individual ␣-conotoxin analogs at the ␣ 7 -GH3 cell line in the Fluo-4/Ca 2؉ assay
The assay was performed as described under "Experimental Procedures," and the data are the means of 3-5 individual experiments performed in duplicate. EC 80 -EC 90 concentrations of acetylcholine were used for the experiments. Fold differences were calculated as follows: (IC 50 WT-ImI/IC 50 analog) Ϫ1 for analogs with increased antagonistic potency compared with WT-ImI (IC 50 analog/IC 50 WT-ImI)ϩ1 for analogs with decreased antagonistic potency compared with WT-ImI. Examples of concentration-inhibitions curves for selected high antagonist potency analogs are shown in supplemental Fig. S1. O 9 -position (Fig. 5a)-The Leu 9 residue in compound 28 was observed to participate in minor hydrophobic interactions with the ␣ 7 receptor, but overall, the residue was in a pocket lined with the following hydrophobic residues: Leu 118 (A) and Trp 54 (A) on one side and Met 159 (A) and Ser 33 (A) on the other. Two polar residues, Gln 56 (A) and Gln 116 (A), also contributed to this pocket. All of these residues were part of only one subunit (shown as subunit A with residues in the text and in cyan in Fig.  5a) of the ␣ 7 nAChR. Given the size of this pocket, it was only partially occupied by Leu 9 from the ␣-conotoxin. From these analyses, we hypothesized that mutations at this position for residues that are either too small and/or with less hydrophobic side chain will further reduce the hydrophobic contacts thus resulting in a loss of affinity/antagonist potency. O 10 -position (Fig. 5b)-This residue is situated in an outside pocket that upon visual inspection appears amenable for accommodating aromatic residues. Aromatic substitutions display cationinteraction with Arg 78 (A) in some of the conformations during the Molecular Dynamics simulation. Overall, the receptor residues surrounding this residue originated from both the subunits of nAChR, as shown in cyan and magenta in Fig. 5b.
O 11 -position (Fig. 5c)-This residue is outside the binding pocket and appears not to participate in any major interactions during the simulation of complex with compound 28.  The structure of compound 36, obtained from NMR analysis, was also modeled into the binding site of homology model of the ␣ 7 receptor (Fig. 6). The analysis of this complex identified several amino acid residues from the receptor that contribute to the ␣-conotoxin-binding site. As established previously from the known crystal and modeled structures, the binding site of the ␣-conotoxin is surrounded by charged aromatic and hydrophobic residues. Notably, as shown in Fig. 6, these residues are as follows: acidic, Asp 163  Water molecules modeled from the crystal structure of AChBP, and optimized along with docked compound 36, appear to play an important role in the binding. These water molecules primarily act as direct or indirect mediators between the binding site residues and compound 36. These interactions cannot be seen in Fig.  6, where only the mutated residues of the ␣-conotoxin are shown (the various interactions of these water molecules are shown in the Molecular Operating Environment interaction map in supplemental Fig. S2).

Functional Characterization of a Third Iteration Library of Individual Compounds
In an effort to further refine the ␣ 7 nAChR activities of the most potent ␣-conotoxin ImI analogs from the second series library, a third series of individual compounds was generated.  Two sub-libraries were synthesized to probe the O 9 -and O 10positions, respectively, based on docking studies of selected analogs into the ␣ 7 nAChR homology model and functional data for second series library analogs at the ␣ 7 nAChR. The first sub-library consisted of a series of analogs probing the O 10position and incorporated six non-natural aromatic amino acids, including 4-nitrophenylalanine (Nph), napth-1-ylalanine (Nal), O-methyltyrosine (Mty), 4-aminophenylalanine (Aph), 2,6-dimethyltyrosine (Dmt), 3-pyridinylalanine (Pyr), as well the nonaromatic hydrophobic amino acids Ile and cyclohexylalanine (Cha). These were used in combination with Leu or Nva residues in the O 9 -position, and His, Abu, or Trp residues in the O 11 -position to give a total of 48 individual analogs. The second sub-library probed the O 9 -position with a series of hydrophobic amino acids including Ile, Val, Cha, Nle, Phe, and Nal residues in combination with Tyr or Trp in the O 10 -position and a residue His in the O 11 -position to give 12 additional analogs. In total, 60 additional individual analogs were synthesized in the third iteration library, and their antagonist potencies were determined at the ␣ 7 nAChR in the Fluo-4/Ca 2ϩ assay (Table 3 and Fig. 3b).
Although a majority of the analogs in the O 9 sub-library displayed medium antagonist potencies, none of these compounds surpassed corresponding analogs containing Nva or Leu at this position. However, profound increases in antagonistic potency were observed for several analogs in the O 10 sub-library. Compounds 49 (Nva 9 -Dmt 10 -His 11 ), 56 (Leu 9 -Aph 10 -Abu 11 ), and 81 (Nva 9 -Dmt 10 -Trp 11 ) all exhibited very high antagonist potencies, with the most active analog, compound 56, exhibiting a more than 14-fold increase compared with WT-ImI. Furthermore, there is trend in antagonistic potency of analogs bearing a strongly activating substituent on the aromatic ring, with several analogs containing Tyr-, Dmt-and Aph-in the O 10 -position exhibiting high or very high antagonist potency. On the other hand, analogs containing Mty, Nph, or Pyr in the O 10 -position all exhibited similar antagonist potencies to WT-ImI. Analogs containing Nal in the O 10 -position displayed significantly increased antagonist potencies across all samples and displayed comparable antagonist potency to corresponding analogs containing Trp. On the other hand, the nonaromatic amino acid Cha and Ile O 10 analogs generally exhibited lower antagonistic potencies. Five of six compounds with Cha at O 10 (compounds 43, 51, 59, 67, and 75) showed low antagonist potency. In general, these observations are in agreement with the trend derived in the modeling where it was concluded that an aromatic residue at O 10 is beneficial for the activity.
The pairwise structure-activity relationships for both the second and third iteration libraries of individual analogs were examined to determine the optimal combinations of amino

Functional properties of the third series of individual ␣-conotoxin analogs at the ␣ 7 -GH3 cell line in the Fluo-4/Ca 2؉ assay
The assay was performed as described under "Experimental Procedures," and the data are the means of 3-5 individual experiments performed in duplicate. EC 80 -EC 90 concentrations of acetylcholine were used for the experiments. Fold differences were calculated as described in Table 2. Examples of concentration-inhibition curves for selected high and very high antagonist potency analogs are shown in supplemental Fig. S1. acid functionalities for the synthesized analogs (Fig. 7). For the O 9 -O 10 combination (Fig. 7a), the Nva 9 -Dmt 10 combination led to two of the three very high antagonist potency analogs. Similarly, the Leu 9 -Aph 10 , Nva 9 -Aph 10 , and Leu 9 -Nal 10 combinations lead the same numbers of analogs displaying high to very high antagonist potencies. Other combinations for O 9 -O 10 that yielded high antagonist potency include Nva 9 -Trp 10 , Nva 9 -Tyr 10 , Leu 9 -Trp 10 , and Leu 9 -Tyr 10 . For the O 9 -O 11 combination (Fig. 7b), all Nva 9 -Trp 11 analogs displayed increased antagonist potency compared with WT-ImI, with four compounds exhibiting high or very high potency and 11 displaying medium antagonist potencies. The Nva 9 -His 11 and Leu 9 -Abu 11 combinations also yielded four compounds with high or very high potency, with one and two low antagonist potency analogs, respectively. For the O 10 -O 11 combination (Fig. 7c), Aph 10 -Abu 11 , Dmt 10 -His 11 , and Dmt 10 -Trp 11 containing analogs each yielded one compound with very high and one compound with medium antagonist potency. Analogs containing the Aph 10 -Trp 11 , Nal 10 -Trp 11 , and Trp 10 -Abu 11 combinations each yielded two high antagonist potency analogs. The most dramatic change in activity is observed for the analogs containing Tyr 10 -Abu 11 (Fig. 7c) that yield one compound with high antagonist potency (compound 36, Leu 9 ; IC 50 ϭ 0.44 M) and one compound regarded as inactive (compound 12, Ala 9 ; IC 50 ϳ 50 M).
Further Functional Characterization of Selected ImI Analogs-The high and very high antagonist potency analogs identified in the second and third iteration libraries were screened for activities at ␣ 4 ␤ 2 and ␣ 3 ␤ 4 nAChR subtypes to assess their selectivity for the ␣ 7 nAChR toward these two other major neuronal nAChR subtypes. This was achieved using cell lines stably expressing the rat ␣ 3 ␤ 4 nAChR and the mouse ␣ 4 ␤ 2 nAChR using the Fluo-4/Ca 2ϩ assay and the fluorescent imaging plate reader membrane potential blue assay, respectively (39,40). No significant antagonist activity was observed for any of these analogs at these two subtypes at concentrations up to 300 M, except for analog 81, which only displayed moderate inhibition (ϳ50%) at a concentration of 300 M (supplemental Table S2). This indicates that all very high and high antagonist potency ImI analogs are characterized by high degrees of selectivity for the ␣ 7 nAChR over ␣ 4 ␤ 2 and ␣ 3 ␤ 4 .

DISCUSSION
To date, a majority of structure-activity relationship studies on ␣-conotoxins as nAChR ligands have largely focused on single amino acid substitutions to study the roles of specific positions within the m-and n-loop regions (8 -16). In this study, we have utilized a three-step synthetic combinatorial strategy to study a specific region (i.e. the n-loop) of ␣-conotoxin ImI to develop novel analogs with improved antagonist properties for the ␣ 7 nAChR.
Synthesis of Mixture-based ␣-Conotoxin Libraries-Although other methods have previously been reported for the high throughput synthesis of conopeptides (52,53), our approach has involved parallel Boc-SPPS using tea bags, twostep low-high HF cleavage with a multiple cleavage manifold and oxidative formation of disulfide bonds. To our knowledge, this is the first study that has applied synthetic mixture-based combinatorial libraries specifically to conotoxins, which has proven to be a time-and cost-effective approach to generating large amounts of information regarding specific residues in a conotoxin sequence.
It is known that mature conotoxin sequences contain all the necessary information to achieve correct in vitro oxidative folding (54). Even though our synthetic approach utilized oxidation conditions optimized for the synthesis of WT-ImI, individual analogs were often accompanied by the presence of other nonnative disulfide bond isomers. Moreover, the formation of disulfide bond isomers in PS-SCL mixtures cannot be determined directly. As such, misfolded conotoxins would be expected to be present in the PS-SCL mixtures. Given that correct conotoxin folding is sequence-driven, misfolded analogs would be expected to be inactive and would not contribute to any activity in the initial PS-SCL screening assay. However, the non-native ribbon isomer of ␣-conotoxin AuIB displays a 10-fold increased activity at the ␣ 3 ␤ 4 nAChR compared with the native ␣-conotoxin (50). In this regard, highly active misfolded isomers in mixtures would be expected to increase the level of diversity of conotoxin PS-SCLs. These can be identified when they are individually synthesized.
Structure-Activity Relationships of the ␣-Conotoxin ImI Analogs at the ␣ 7 nAChR-Through screening and deconvolution of a PS-SCL via a fluorescence-based functional assay at the ␣ 7 nAChR, specific amino acid residue substitutions in the n-loop were identified as promising candidates for the development of more potent ␣-conotoxin ImI analog antagonists.
Of the total 96 individual analogs synthesized in this study, three analogs were found to exhibit very high potency, 11 high potency, 58 medium potency, and 20 low potency, and four analogs were found to be inactive as ␣ 7 nAChR antagonists. Furthermore, a selection of the most potent ␣ 7 nAChR antagonists identified in this series was found to be highly selective for this subtype over the ␣ 3 ␤ 4 and ␣ 4 ␤ 2 subtypes. The most significant increases in antagonist potency were observed for analogs containing the Nva 9 -Dmt 10 -His 11 , Leu 9 -Aph 10 -Abu 11 , and Nva 9 -Dmt 10 -Trp 11 combinations (compounds 49, 56, and 81), which exhibited ϳ12-, 14-, and 10-fold increases, respectively, when compared with WT-ImI. To our knowledge, these compounds represent the most potent ␣-conotoxin ImI analogs reported to date.
One of the most notable observations arising from this study is the complexity of the contribution of the n-loop of ␣-conotoxin ImI to its binding to ␣ 7 nAChR, where no single substitution resulted in increased antagonist potency. Rather, substitutions acted synergistically to achieve significant increases in antagonist potency, underlining the efficacy of our combinatorial approach at multiple positions. Although this suggests that the functionalities of the n-loop residues in native toxin are quite optimized for its interactions with the receptors, the antagonist activity of it can be significantly enhanced by combinations of two or three residue substitutions in the loop.
Substitutions of Ala with Nva or Leu residues in the O 9 -position of the peptide were found to be optimal for ␣ 7 nAChR activity in the PS-SCL mixtures. This constituted the rationale for the focus on these residues as O 9 substitutions in the individual analogs in the second and third iteration libraries.
Molecular modeling of selected compounds from the second iteration library revealed that small amino acids such as Ala at O 9 might lack sufficient hydrophobic interactions within the binding sub-pocket. Actually, only Nva stood out in the PS-SCL analysis, as the functional data for the Leu 9 analog mixture did not differ significantly from many other PS-SCL O 9 mixtures ( Table 1). The superiority of Nva over Leu in this position translated from the PS-SCL into the individually synthesized analogs.
The presence of an aromatic residue in the O 10 -position of ␣-conotoxin ImI, as in the case of the native Trp, was observed to be crucial for optimal receptor binding. This observation is supported by the fact that analogs containing Cha and Ile in the O 10 -position typically display lower antagonist potencies than their corresponding Trp, Nph, Nal, Mty, Aph, Dmt, and Pyr analogs. Deviations from this general trend do exist in the series, however, as the Ile 10 -containing analogs 52 and 84 are quite potent as ␣ 7 antagonists (Table 3).
Another observation regarding the O 10 residue is that large aromatic side chain functionalities such as the native Trp and Nal in this position generally increase the antagonistic potency of the analogs. This observation might be associated with cationinteraction between the aromatic residue at O 10 and Arg 78 (A) in the binding site. Significantly, analogs with aromatic residues bearing strongly activating moieties that increase electron density into the aromatic ring (i.e. Tyr, -OH; Dmt, -OH, and -CH 3 ; and Aph, -NH 2 substituents) in the O 10 -position also yield considerable increases in antagonistic potency, indicating stronger cation-interactions with Arg 78 (A) in the receptor-binding site. As such, each of the very high antagonist potency analogs (compounds 49, 56, and 81) possessed either Dmt or Aph at position O 10 . Conversely, Mty, Nph, and Pyr containing analogs each exhibited no significant difference in antagonistic potency, which is possibly a result of a deactivated aromatic ring that weakens the interaction with Arg 78 (A).
In general, substitutions in the O 11 -position in the n-loop had minor effects on antagonistic potency. Furthermore, it can be confirmed by comparisons of the IC 50 values of the various groups of three analogs containing the same residues in the O 9and O 10 -positions and having Trp, His, or Abu residues in the O 11 -position. Finally, Abu, His, and Trp as O 11 residues were all found in the four very high antagonist potency analogs. This can be explained in terms of the O 11 -position being oriented away from the binding site and not taking part in any significant interactions with the receptor per se. Thus, the functional data of the analogs are in agreement with the observations from the molecular modeling. However, it should be stressed that the O 11 -position was not probed to the same extent in the third series of compounds as the O 9 -and O 10 -positions.
One of the most striking outcomes from this study is that the most potent individual compound predicted from the screening data (Nva 9 -Trp 10 -His 11 , compound 14) was among the most potent analogs in the second series of compounds. This result underlines the efficacy of PS-SCL mixtures for identifying analogs with amino acid combinations giving rise to improved functional activities for a given pharmacological target.
Additional probing of other regions of ␣-conotoxin ImI could also lead to higher potency antagonists. For example, we recently reported that substitution of Pro 6 with 5-(R)-phenylproline in ␣-conotoxin ImI induces -stacking interactions with another underutilized binding pocket (14). As such, we speculate that substituting Pro 6 with 5-(R)-phenylproline in the very high potency antagonists arising from this study would result in even greater increases in antagonistic potency at the ␣ 7 nAChR.
Future Perspectives-Given the importance of the ␣ 7 nAChR in memory, cognition, and in overall mental health (5,6), the development of potent and selective ligands for this receptor would be of major significance for researchers studying its role in neuropathological conditions such as schizophrenia, Alzheimer, and Parkinson disease. Although this study has focused on the development of ␣-conotoxin ImI analogs with increased antagonist activity and subtype selectivity at the ␣ 7 nAChR, it is expected that screening of the ␣-conotoxin ImI n-loop PS-SCL for activity at the ␣ 9 ␣ 10 nAChR will identify a different set of active amino acid functionalities, allowing the design of individual analogs with enhanced activities for this receptor. Finally, the construction of PS-SCLs based on other classes of conotoxins seems to be a promising strategy to develop new subtype-selective ligands targeting other nAChR subtypes, just as other disulfide-rich venom peptides could be used as scaffolds in PS-SCL-based design for other therapeutically relevant targets.