Alpha7 nicotinic acetylcholine receptor agonists: prediction of their binding affinity through a molecular mechanics Poisson-Boltzmann surface area approach

Terahertz Photons Improve Cognitive Functions in Posttraumatic Stress Disorder

Posttraumatic stress disorder (PTSD) is a serious psychosis leading to cognitive impairment. To restore cognitive functions for patients, the main treatments are based on medication or rehabilitation training but with limited effectiveness and strong side effects. Here, we demonstrate a new treatment approach for PTSD by using terahertz (THz) photons stimulating the hippocampal CA3 subregion. We verified that this method can nonthermally restore cognitive function in PTSD rats in vivo. After THz photon irradiation, the PTSD rats' recognitive index improved by about 10% in a novel object recognition test, the PTSD rats' accuracy improved by about 100% in a shuttler box test, the PTSD rats' numbers to identify target box was about 5 times lower in a Barnes maze test, and the rate of staying in new arm increased by approximately 40% in a Y-maze test. Further experimental studies found that THz photon (34.5 THz) irradiation could improve the expression of NR2B (increased by nearly 40%) and phosphorylated NR2B (increased by about 50%). In addition, molecular dynamics simulations showed that THz photons at a frequency of 34.5 THz are mainly absorbed by the pocket of glutamate receptors rather than by glutamate molecules. Moreover, the binding between glutamate receptors and glutamate molecules was increased by THz photons. This study offers a nondrug, nonthermal approach to regulate the binding between the excitatory neurotransmitter (glutamate) and NR2B. By increasing synaptic plasticity, it effectively improves the cognitive function of animals with PTSD, providing a promising treatment strategy for NR2B-related cognitive disorders.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors

Changeux et al. recently suggested that the SARS-CoV-2 spike (S) protein may interact with nicotinic acetylcholine receptors (nAChRs). Such interactions may be involved in pathology and infectivity. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the S protein, to investigate this 'nicotinic hypothesis'. We examine the binding of the Y674-R685 loop of the S protein to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results indicate that Y674-R685 has affinity for nAChRs and the region responsible for binding contains the PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. In particular, R682 has a key role in the stabilisation of the complexes as it forms interactions with loops A, B and C in the receptor's binding pocket. The conformational behaviour of the bound Y674-R685 region is highly dependent on the receptor subtype, adopting extended conformations in the α4β2 and α7 complexes and more compact ones when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of the simulations of the full-length S protein show that the Y674-R685 region is accessible for binding, and suggest a potential binding orientation of the S protein with nAChRs.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

In Silico Modeling of the α7 Nicotinic Acetylcholine Receptor: New Pharmacological Challenges Associated with Multiple Modes of Signaling

The α7 nicotinic acetylcholine receptor is a homopentameric ion-channel of the Cys-loop superfamily characterized by its low probability of opening, high calcium permeability, and rapid desensitization. The α7 receptor has been targeted for the treatment of the cognitive symptoms of schizophrenia, depression, and Alzheimer's disease, but it is also involved in inflammatory modulation as a part of the cholinergic anti-inflammatory pathway. Despite its functional importance, in silico studies of the α7 receptor cannot produce a general model explaining the structural features of receptor activation, nor predict the mode of action for various ligand classes. Two particular problems in modeling the α7 nAChR are the absence of a high-resolution structure and the presence of five potentially nonequivalent orthosteric ligand binding sites. There is wide variability regarding the templates used for homology modeling, types of ligands investigated, simulation methods, and simulation times. However, a systematic survey focusing on the methodological similarities and differences in modeling α7 has not been done. In this work, we make a critical analysis of the modeling literature of α7 nAChR by comparing the findings of computational studies with each other and with experimental studies under the main topics of structural studies, ligand binding studies, and comparisons with other nAChR. In light of our findings, we also summarize current problems in the field and make suggestions for future studies concerning modeling of the α7 receptor.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

Lupin Peptides Modulate the Protein-Protein Interaction of PCSK9 with the Low Density Lipoprotein Receptor in HepG2 Cells

Proprotein convertase subtilisin/kexin type 9 (PCSK9) has been recently identified as a new useful target for hypercholesterolemia treatment. This work demonstrates that natural peptides, deriving from the hydrolysis of lupin protein and absorbable at intestinal level, are able to inhibit the protein-protein interaction between PCSK9 and the low density lipoprotein receptor (LDLR). In order to sort out the best potential inhibitors among these peptides, a refined in silico model of the PCSK9/LDLR interaction was developed. Docking, molecular dynamics (MD) simulations and peptide binding energy estimations, by MM-GBSA approach, permitted to select the two best candidates among tested peptides that were synthesized and evaluated for their inhibitory activity. The most active was P5 that induced a concentration dependent inhibition of the PCSK9-LDLR binding, with an IC₅₀ value equal to 1.6 ± 0.33 μM. Tested at a 10 μM concentration, this peptide increased by 66 ± 21.4% the ability of HepG2 cells to take up LDL from the extracellular environment.

Allosteric Modulation of Alpha7 Nicotinic Receptors: Mechanistic Insight through Metadynamics and Essential Dynamics

Increasing attention has recently been devoted to allosteric modulators, as they can provide inherent advantages over classic receptor agonists. In the field of nicotinic receptors (nAChRs), the main advantage is that allosteric modulators can trigger pharmacological responses, limiting receptor desensitization. Most of the known allosteric ligands are "positive allosteric modulators" (PAMs), which increase both sensitivity to receptor agonists and current amplitude. Intriguingly, some allosteric modulators are also able to activate the α7 receptor (α7-nAChR) even in the absence of orthosteric agonists. These compounds have been named "ago-allosteric modulators" and GAT107 has been studied in depth because of its unique mechanism of action. We here investigate by molecular dynamics simulations, metadynamics, and essential dynamics the activation mechanism of α7-nAChR, in the presence of different nicotinic modulators. We determine the free energy profiles associated with the closed-to-open motion of the loop C, and we highlight mechanistic differences observed in the presence of different modulators. In particular, we demonstrate that GAT107 triggers conformational motions and cross-talk similar to those observed when the α7-nACh receptor is in complex with both an agonist and an allosteric modulator.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

Nicotine-modulated subunit stoichiometry affects stability and trafficking of α3β4 nicotinic receptor

Heteromeric nAChRs are pentameric cation channels, composed of combinations of two or three α and three or two β subunits, which play key physiological roles in the central and peripheral nervous systems. The prototypical agonist nicotine acts intracellularly to upregulate many nAChR subtypes, a phenomenon that is thought to contribute to the nicotine dependence of cigarette smokers. The α3β4 subtype has recently been genetically linked to nicotine dependence and lung cancer; however, the mode of action of nicotine on this receptor subtype has been incompletely investigated. Here, using transfected mammalian cells as model system, we characterized the response of the human α3β4 receptor subtype to nicotine and the mechanism of action of the drug. Nicotine, when present at 1 mm concentration, elicited a ∼5-fold increase of cell surface α3β4 and showed a more modest upregulatory effect also at concentrations as low as 10 μM. Upregulation was obtained if nicotine was present during, but not after, pentamer assembly and was caused by increased stability and trafficking of receptors assembled in the presence of the drug. Experimental determinations as well as computational studies of subunit stoichiometry showed that nicotine favors assembly of pentamers with (α3)2(β4)3 stoichiometry; these are less prone than (α3)3(β4)2 receptors to proteasomal degradation and, because of the presence in the β subunit of an endoplasmic reticulum export motif, more efficiently transported to the plasma membrane. Our findings uncover a novel mechanism of nicotine-induced α3β4 nAChR upregulation that may be relevant also for other nAChR subtypes.

Mechanism of falcipain-2 inhibition by α,β-unsaturated benzo[1,4]diazepin-2-one methyl ester

Falcipain-2 (FP-2) is a papain-family cysteine protease of Plasmodium falciparum whose primary function is to degrade the host red cell hemoglobin, within the food vacuole, in order to provide free amino acids for parasite protein synthesis. Additionally it promotes host cell rupture by cleaving the skeletal proteins of the erythrocyte membrane. Therefore, the inhibition of FP-2 represents a promising target in the search of novel anti-malarial drugs. A potent FP-2 inhibitor, characterized by the presence in its structure of the 1,4-benzodiazepine scaffold and an α,β-unsaturated methyl ester moiety capable to react with the Cys42 thiol group located in the active site of FP-2, has been recently reported in literature. In order to study in depth the inhibition mechanism triggered by this interesting compound, we carried out, through ONIOM hybrid calculations, a computational investigation of the processes occurring when the inhibitor targets the enzyme and eventually leads to an irreversible covalent Michael adduct. Each step of the reaction mechanism has been accurately characterized and a detailed description of each possible intermediate and transition state along the pathway has been reported.

Enantioselective drug–protein interaction between mexiletine and plasma protein

Objectives

This study examined the interaction of mexiletine enantiomers with human plasma, human serum albumin (HSA), and human α1-acid glycoprotein (hAGP), and characterized the binding modes of mexiletine enantiomers with hAGP in the molecular level.

Methods

Enantiomer separation of mexiletine was performed using precolumn derivatization chiral HPLC. The ultrafiltration technique was used to separate the free mexiletine in plasma matrix. Molecular dynamics simulations and free energy calculations were assessed using molecular mechanics and the generalized Born surface area method.

Key findings

Significant differences in enantioselective binding to human plasma were observed (R > S). The hAGP–mexiletine binding profile exhibited similar enantioselectivity (R > S) to that in human plasma, whereas HSA–mexiletine interaction was S > R at pH 7.4. Moreover, the results of comparative studies indicated that mexiletine had the highest binding affinity for F1-S, a variant of hAGP. Based on the computational studies, residues such as Arg90, Leu79, Ser89 and Phe89 showed an energy difference of more than −0.35 kcal/mol between the enantiomers.

Conclusions

hAGP may be one of the key proteins leading to the enantioselective protein bindings of mexiletine in human plasma (R > S). The residues Arg90, Leu79, Ser89 and Phe89 of hAGP may have important roles in the observed enantioselectivity.

Application of conformational clustering in protein-ligand docking

Protein-Ligand docking is a powerful technique routinely employed in structure-based drug design. Despite many reported success stories, docking is not always able to provide an accurate and easily interpretable prediction of the structure of the bound complex formed by a small organic molecule and a pharmacologically relevant target. Cluster analysis can represent a versatile and readily available postprocessing tool to be employed in combination with protein-ligand docking to simplify the evaluation of the results and help to overcome present limitations of docking protocols.

Identification of a negative allosteric site on human α4β2 and α3β4 neuronal nicotinic acetylcholine receptors

Acetylcholine-based neurotransmission is regulated by cationic, ligand-gated ion channels called nicotinic acetylcholine receptors (nAChRs). These receptors have been linked to numerous neurological diseases and disorders such as Alzheimer's disease, Parkinson's disease, and nicotine addiction. Recently, a class of compounds has been discovered that antagonize nAChR function in an allosteric fashion. Models of human α4β2 and α3β4 nicotinic acetylcholine receptor (nAChR) extracellular domains have been developed to computationally explore the binding of these compounds, including the dynamics and free energy changes associated with ligand binding. Through a blind docking study to multiple receptor conformations, the models were used to determine a putative binding mode for the negative allosteric modulators. This mode, in close proximity to the agonist binding site, is presented in addition to a hypothetical mode of antagonism that involves obstruction of C loop closure. Molecular dynamics simulations and MM-PBSA free energy of binding calculations were used as computational validation of the predicted binding mode, while functional assays on wild-type and mutated receptors provided experimental support. Based on the proposed binding mode, two residues on the β2 subunit were independently mutated to the corresponding residues found on the β4 subunit. The T58K mutation resulted in an eight-fold decrease in the potency of KAB-18, a compound that exhibits preferential antagonism for human α4β2 over α3β4 nAChRs, while the F118L mutation resulted in a loss of inhibitory activity for KAB-18 at concentrations up to 100 µM. These results demonstrate the selectivity of KAB-18 for human α4β2 nAChRs and validate the methods used for identifying the nAChR modulator binding site. Exploitation of this site may lead to the development of more potent and subtype-selective nAChR antagonists which may be used in the treatment of a number of neurological diseases and disorders.

Design, synthesis, and pharmacological characterization of novel spirocyclic quinuclidinyl-Δ2-isoxazoline derivatives as potent and selective agonists of α7 nicotinic acetylcholine receptors

A set of racemic spirocyclic quinuclidinyl-Δ(2)-isoxazoline derivatives was synthesized using a 1,3-dipolar cycloaddition-based approach. Target compounds were assayed for binding affinity toward rat neuronal homomeric (α7) and heteromeric (α4β2) nicotinic acetylcholine receptors. Δ(2) -Isoxazolines 3 a (3-Br), 6 a (3-OMe), 5 a (3-Ph), 8 a (3-OnPr), and 4 a (3-Me) were the ligands with the highest affinity for the α7 subtype (K(i) values equal to 13.5, 14.2, 25.0, 71.6, and 96.2 nM, respectively), and showed excellent α7 versus α4β2 subtype selectivity. These compounds, tested in electrophysiological experiments against human α7 and α4β2 receptors stably expressed in cell lines, behaved as partial α7 agonists with varying levels of potency. The two enantiomers of (±)-3-methoxy-1-oxa-2,7-diaza-7,10-ethanospiro[4.5]dec-2-ene sesquifumarate 6 a were prepared using (+)-dibenzoyl-L- or (-)-dibenzoyl-D-tartaric acid as resolving agents. Enantiomer (R)-(-)-6 a was found to be the eutomer, with K(i) values of 4.6 and 48.7 nM against rat and human α7 receptors, respectively.

Different interaction between the agonist JN403 and the competitive antagonist methyllycaconitine with the human alpha7 nicotinic acetylcholine receptor

The interaction of the agonist JN403 with the human (h) alpha7 nicotinic acetylcholine receptor (AChR) was compared to that for the competitive antagonist methyllycaconitine (MLA). The receptor selectivity of JN403 was studied on the halpha7, halpha3beta4, and halpha4beta2 AChRs. The results established that the cationic center and the hydrophobic group found in JN430 and MLA are important for the interaction with the AChRs. MLA preincubation inhibits JN403-induced Ca(2+) influx in GH3-halpha7 cells with a potency 160-fold higher than that when MLA is co-injected with JN403. The most probable explanation, based on our dynamics results, is that MLA (more specifically the 3-methyl-2,5-dioxopyrrole ring and the B-D rings) stabilizes the resting conformational state. The order of receptor specificity for JN403 is as follows: halpha7 > halpha3beta4 ( approximately 40-fold) > halpha4beta2 ( approximately 500-fold). This specificity is based on a larger number of hydrogen bonds between the carbamate group (another pharmacophore) of JN403 and the halpha7 sites, the electrostatic repulsion between the positively charged residues around the halpha3beta4 sites and the cationic center of JN403, fewer hydrogen bonds for the interaction of JN403 with the halpha3beta4 AChR, and an unfavorable van der Waals interaction between JN403 and the alpha4-beta2 interface. The higher receptor specificity for JN403 could be important for the treatment of alpha7-related disorders, including dementias, pain-related ailments, depression, anxiety, and wound healing.

How to obtain statistically converged MM/GBSA results

The molecular mechanics/generalized Born surface area (MM/GBSA) method has been investigated with the aim of achieving a statistical precision of 1 kJ/mol for the results. We studied the binding of seven biotin analogues to avidin, taking advantage of the fact that the protein is a tetramer with four independent binding sites, which should give the same estimated binding affinities. We show that it is not enough to use a single long simulation (10 ns), because the standard error of such a calculation underestimates the difference between the four binding sites. Instead, it is better to run several independent simulations and average the results. With such an approach, we obtain the same results for the four binding sites, and any desired precision can be obtained by running a proper number of simulations. We discuss how the simulations should be performed to optimize the use of computer time. The correlation time between the MM/GBSA energies is approximately 5 ps and an equilibration time of 100 ps is needed. For MM/GBSA, we recommend a sampling time of 20-200 ps for each separate simulation, depending on the protein. With 200 ps production time, 5-50 separate simulations are required to reach a statistical precision of 1 kJ/mol (800-8000 energy calculations or 1.5-15 ns total simulation time per ligand) for the seven avidin ligands. This is an order of magnitude more than what is normally used, but such a number of simulations is needed to obtain statistically valid results for the MM/GBSA method.

Insights into docking and scoring neuronal alpha4beta2 nicotinic receptor agonists using molecular dynamics simulations and QM/MM calculations

A combined quantum mechanical (QM)-polarized docking and molecular dynamics approach to study the binding mode and to predict the binding affinity of ligands acting at the alpha4beta2-nAChR is presented. The results obtained in this study indicate that the quantum mechanical/molecular mechanics docking protocol well describes the charge-driven interactions occurring in the binding of nicotinic agonists, and it is able to represent the polarization effects on the ligand exerted by the surrounding atoms of the receptor at the binding site. This makes it possible to properly score agonists of alpha4beta2-nAChR and to reproduce the experimental binding affinity data with good accuracy, within a mean error of 2.2 kcal/mol. Moreover, applying the QM-polarized docking to an ensemble of nAChR conformations obtained from MD simulations enabled us to accurately capture nAChR-ligand induced-fit effects.

5-(2-Pyrrolidinyl)oxazolidinones and 2-(2-pyrrolidinyl)benzodioxanes: synthesis of all the stereoisomers and alpha4beta2 nicotinic affinity

The four stereoisomers of 2-oxazolidinone 5-substituted with 1-methyl-2-pyrrolidinyl (1), of 1,4-benzodioxane 2-substituted with the same residue (2) and of the nor-methyl analogue of this latter (2a) were synthesized as candidate nicotinoids. Of the 12 compounds, two N-methylated pyrrolidinyl-benzodioxane stereoisomers, namely those with the same relative configuration at the pyrrolidine stereocentre as (S)-nicotine, bind at alpha4beta2 nicotinic acetylcholine receptor with submicromolar affinity. Consistently with the biological data, docking analysis enlightens significant differences in binding site interactions not only between 1 and 2, but also between 2 and 2a and between the stereoisomers of 2 accounting for the critical role played, in the case of the pyrrolidinyl-benzodioxanes, by the chirality of both the stereolabile and stereostable stereogenic atoms, namely the protonated tertiary nitrogen and the two asymmetric carbons.