Mechanism of intermolecular hydroacylation of vinylsilanes catalyzed by a rhodium(I) olefin complex: a DFT study

Density functional theory (DFT) was used to investigate the Rh(I)-catalyzed intermolecular hydroacylation of vinylsilane with benzaldehyde. All intermediates and transition states were optimized completely at the B3LYP/6-31G(d,p) level (LANL2DZ(f) for Rh). Calculations indicated that Rh(I)-catalyzed intermolecular hydroacylation is exergonic, and the total free energy released is −110 kJ mol⁻¹. Rh(I)-catalyzed intermolecular hydroacylation mainly involves the active catalyst CA2, rhodium–alkene–benzaldehyde complex M1, rhodium–alkene–hydrogen–acyl complex M2, rhodium–alkyl–acyl complex M3, rhodium–alkyl–carbonyl–phenyl complex M4, rhodium–acyl–phenyl complex M5, and rhodium–ketone complex M6. The reaction pathway CA2 + R2 → M1b → T1b → M2b → T2b1 → M3b1 → T4b → M4b → T5b → M5b → T6b → M6b → P2 is the most favorable among all reaction channels of Rh(I)-catalyzed intermolecular hydroacylation. The reductive elimination reaction is the rate-determining step for this pathway, and the dominant product predicted theoretically is the linear ketone, which is consistent with Brookhart’s experiments. Solvation has a significant effect, and it greatly decreases the free energies of all species. The use of the ligand Cp′ (Cp′ = C₅Me₄CF₃) decreased the free energies in general, and in this case the rate-determining step was again the reductive elimination reaction.     

A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances

Abstract

In this study for the first time we have revealed and investigated in details 123 different prototropic tautomers of the most stable conformer of the quercetin molecule using quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of QM theory. We have found that in the most energetically favorable prototropic tautomer mobile hydrogen atoms are localized at the О3, О3′, О4′, О5, and О7 exocyclic oxygen atoms. Molecular tautomers are in the range of the Gibbs free energies from 0.0 to 69.8?kcal·mol^?1, while zwitterionic ones – from 30.1 до 172.8?kcal·mol^?1 at normal conditions. It was also reliably established that the weakest point causing the decyclization of the molecule is its C ring – this reaction is launched by the transition of the proton from the C8H group to the endocyclic O1 oxygen atom. All prototropic tautomers, except two cases, are joined by the intramolecular cooperative specific interactions (from 1 to 5) – H-bonds and attractive van der Waals contacts, which have been revealed and characterized by QTAIM analysis.

Communicated by Ramaswamy H. Sarma     

Theoretical investigation of tautomerism in N-hydroxy amidines

A DFT study with QST3 approach method is used to calculate kinetic, thermodynamic, spectral and structural data of tautomers and transition state structures of some N-hydroxy amidines. All tautomers and transition states are optimized at the B3LYP/6-311++g** and B3LYP/aug-cc-pvtz level, with good agreement in energetic result with energies obtained from CBS-QB3, a complete basis set composite energy method. The result shows that the tautomer a (amide oxime) is more stable than the tautomer b (imino hydroxylamine) as is reported in the literature. In addition, our finding shows that, the energy difference between two tautomers is only in about 4–10 kcal/mol but the barrier energy found in traversing each tautomer to another one is in the range of 33–71 kcal/mol. Therefore, it is impossible to convert these two tautomers to each other at room temperature. Additionally, transition state theory is applied to estimate the barrier energy and reaction rate constants of the hydrogen exchange between tautomers in presence of 1–3 molecules of water. The computed activation barrier shows us that the barrier energy of solvent assisted tautomerism is about 9–20 kcal/mol and lower than simple tautomerism and this water-assisted tautomerism is much faster than simple tautomerism, especially with the assisting two molecules of water.     

Fast cyclic electron transport around photosystem I in leaves under far-red light: a proton-uncoupled pathway?

Fast cyclic electron transport (CET) around photosystem I (PS I) was observed in sunflower (Helianthus annuus L.) leaves under intense far-red light (FRL) of up to 200 μmol quanta m⁻² s⁻¹. The electron transport rate (ETR) through PS I was found from the FRL-dark transmittance change at 810 and 950 nm, which was deconvoluted into redox states and pool sizes of P700, plastocyanin (PC) and cytochrome f (Cyt f). PC and P700 were in redox equilibrium with K _e = 35 (ΔE _m = 90 mV). PS II ETR was based on O₂ evolution. CET [(PS I ETR) − (PS II ETR)] increased to 50–70 μmol e⁻ m⁻² s⁻¹ when linear electron transport (LET) under FRL was limited to 5 μmol e⁻ m⁻² s⁻¹ in a gas phase containing 20–40 μmol CO₂ mol⁻¹ and 20 μmol O₂ mol⁻¹. Under these conditions, pulse-saturated fluorescence yield F _m was non-photochemically quenched; however, F _m was similarly quenched when LET was driven by low green or white light, which energetically precluded the possibility for active CET. We suggest that under FRL, CET is rather not coupled to transmembrane proton translocation than the CET-coupled protons are short-circuited via proton channels regulated to open at high ΔpH. A kinetic analysis of CET electron donors and acceptors suggests the CET pathway is that of the reversed Q-cycle: Fd → (FNR) → Cyt c_n → Cyt b_h → Cyt b_l → Rieske FeS → Cyt f → PC → P700 →→ Fd. CET is activated when PQH₂ oxidation is opposed by high ΔpH, and ferredoxin (Fd) is reduced due to low availability of e⁻ acceptors. The physiological significance of CET may be photoprotective, as CET may be regarded as a mechanism of energy dissipation under stress conditions.     

Molecular modeling of dissociative and non-dissociative chemisorption of nitrosamine on close-ended and open-ended pristine and Stone-Wales defective (5,5) armchair single-walled carbon nanotubes

The nitrosamine adsorbed on close-ended and open-ended pristine and Stone-Wales defective (5,5) armchair single-walled carbon nanotubes (SWCNTs) was studied using the B3LYP/6-31G(d) method. Structure optimization of all possible adsorption configurations based on the combination of two nitrosamine (amino- and imino-) isomers and four types of nanotubes was carried out. The most stable configuration for the nitrosamine adsorbed on the (5,5) armchair SWCNTs was found to be dissociative chemisorption. The adsorption energies of the most stable structures of the adsorption complexes of close-ended and open-ended pristine SWCNTs with the imino isomer of nitrosamine were −127.15 and −137.14 kcal mol⁻¹, respectively.     

Electrochemical bromochlorination of peracetylated glycals

Peracetylated glycals—3,4,5-tri-O-acetyl-d-glucal (1a), 3,4,5-tri-O-acetyl-d-galactal (1b) and 3,4-di-O-acetyl-6-deoxy-l-glucal (1c)—have been bromochlorinated by a suitable halogenating agent, generated electrochemically from a mixture of bromides and chlorides in dichloromethane. The reaction was performed in two ways: (i) by a constant current electrolysis (2 F mol⁻¹) of bromides and substrates in a milieu containing an excess of chlorides (Br^?/1/Cl^? = 1:1:6.8) and (ii) by anodic generation of free chlorine from chlorides (2 F mol⁻¹) and subsequent addition of bromides and substrates in a ratio Br^?/1 = 1:1. The corresponding 2-bromo-2-deoxy-glycopyranosyl chlorides were obtained in high yields.     

Base Pairing at the Stem-loop Junction in the SL1 Kissing Complex of HIV-1 RNA: A Thermodynamic Study Probed by Molecular Dynamics Simulation

Abstract

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1_Lai genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal·mol^?1, ?11 kcal·mol^?1and ?0.037 kcal·mol^?1-K^?1 at 300° K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.     

Biotransformation of arenobufagin and cinobufotalin by Alternaria alternata

Abstract

The biotransformation of arenobufagin (1) and cinobufotalin (2), isolated from the natural medicine Chan Su, by Alternaria alternata AS 3.4578 was carried out. Incubation of 1 and 2 afforded six metabolites: 3-oxo-arenobufagin (1a), ψ-bufarenogin (1b), 3-oxo-ψ-bufarenogin (1c), 3-oxo-Δ⁴-derivative of cinobufotalin (2a), 3-oxo-cinobufotalin (2b) and 12β-hydroxycinobufotalin (2c). Among them, metabolites 1a, 1c and 2c are new compounds and their structures were characterized on the basis of their spectroscopic data (NMR, MS and IR). Compounds 1, 2, 1b, 2a and 2b were evaluated for their cytotoxicity against HepG2 and MCF-7 human cancer cells, and all of them showed significant inhibitory activities.     

A new acylated flavonol diglycoside from <Emphasis Type="Italic">Atriplex littoralis</Emphasis>

A new acetylated flavonol glycoside: patuletin 3-O-[5′″-O-feruloyl-β-D-apiofuransyl (1′″→2′′)-β-D-glucopyranoside] (2), together with a known patuletin 3-O-β-D-glucopyranoside (1) were isolated from the aerial part of Artiplex littoralis L. (Chenopodiacease). Their structures were elcidated by acid hydrolysis and spectroscopic methods including UV, ¹H, ¹³C NMR and ESI-MS for both compounds, additionally 2D-NMR, HSQC, HMBC experiments were performed for 2.     

A biophysical and computational study unraveling the molecular interaction mechanism of a new Janus kinase inhibitor Tofacitinib with bovine serum albumin

The interaction of a recently certified kinase inhibitor Tofacitinib (TFB) with bovine serum albumin (BSA) has been studied, by spectroscopic and molecular docking studies. Spectrofluorimetric measurements at 3 different temperatures (288, 298, and 310 K) showed that TFB quench the intrinsic fluorescence of BSA upon forming a nonfluorescent complex. The intrinsic fluorescence data showed that TFB binds to BSA with binding constant (K _b) of approximately 10⁴M⁻¹, affirming a significant affinity of TFB with BSA. The decrease in Stern‐Volmer quenching constant with increasing temperature exhibited the static mechanism of quenching. Negative value of ΔG (−6.94 ± 0.32 kcal·mol⁻¹), ΔH (−7.87 ± 0.52 kcal·mol⁻¹), and ΔS (−3.14 ± 0.42 cal·mol⁻¹·K⁻¹) at all 3 temperatures declared the reaction between BSA and TFB to be spontaneous and exothermic. Far‐UV circular dichroism spectroscopy results demonstrated an increase in helical content of BSA in the presence of TFB. Moreover, dynamic light scattering measurements showed that TFB resulted into a decrease in the hydrodynamic radii (from 3.6 ± 0.053 to 2.9 ± 0.02 nm) of BSA. Molecular docking studies confirmed that TFB binds near site II on BSA, hydrogen bonding, and hydrophobic interaction were involved in the BSA‐TFB complex formation. The present study characterizing the BSA‐TFB interaction could be significant towards gaining an insight into the drug pharmacokinetics and pharmacodynamics and also in the direction of rational drug designing with better competence, against emerging immune‐mediated diseases, ie, alopecia and rheumatoid arthritis.     

DPT tautomerization of the long A∙A* Watson-Crick base pair formed by the amino and imino tautomers of adenine: combined QM and QTAIM investigation

Combining quantum-mechanical (QM) calculations with quantum theory of atoms in molecules (QTAIM) and using the methodology of sweeps of the energetic, electron-topological, geometric and polar parameters, which describe the course of the tautomerization along the intrinsic reaction coordinate (IRC), we showed for the first time that the biologically important A?A* base pair (C_s symmetry) formed by the amino and imino tautomers of adenine (A) tautomerizes via asynchronous concerted double proton transfer (DPT) through a transition state (TS), which is the A⁺?A^? zwitterion with the separated charge, with C_s symmetry. The nine key points, which can be considered as electron-topological “fingerprints” of the asynchronous concerted A?A*?A*?A tautomerization process via the DPT, were detected and completely investigated along the IRC of the A?A*?A*?A tautomerization. Based on the sweeps of the H-bond energies, it was found that intermolecular antiparallel N6Н?N6 (7.01 kcal mol^?1) and N1H?N1 (6.88 kcal mol^?1) H-bonds are significantly cooperative and mutually reinforce each other. It was shown for the first time that the A?A*?A*?A tautomerization is assisted by the third C2H?HC2 dihydrogen bond (DHB), which, in contrast to the two others N6H?N6 and N1H?N1 H-bonds, exists within the IRC range from ?2.92 to 2.92 Å. The DHB cooperatively strengthens, reaching its maximum energy 0.42 kcal mol^?1 at IRC?=??0.52 Å and minimum energy 0.25 kcal mol^?1 at IRC?=??2.92 Å, and is accompanied by strengthening of the two other aforementioned classical H-bonds. We established that the C2H?HC2 DHB completely satisfies the electron-topological criteria for H-bonding, in particular Bader’s and all eight “two-molecule” Koch and Popelier’s criteria. The positive value of the Grunenberg’s compliance constant (5.203 Å/mdyn) at the TS_A?A*?A*?A proves that the C2H?HC2 DHB is a stabilizing interaction. NBO analysis predicts transfer of charge from σ(C2–H) bonding orbital to σ*(H–C2) anti-bonding orbital; at this point, the stabilization energy E⁽²⁾ is equal to 0.19 kcal mol^?1 at the TS_A?A*?A*?A.     

Design,synthesis, molecular docking and biological evaluation of thiophen-2-iminothiazolidine derivatives for use against Trypanosoma cruzi

In this study, we designed and synthesized a series of thiophen-2-iminothiazolidine derivatives from thiophen-2-thioureic with good anti-Trypanosoma cruzi activity. Several of the final compounds displayed remarkable trypanocidal activity. The ability of the new compounds to inhibit the activity of the enzyme cruzain, the major cysteine protease of T. cruzi, was also explored. The compounds 3b, 4b, 8b and 8c were the most active derivatives against amastigote form, with significant IC₅₀ values between 9.7 and 6.03 μM. The 8c derivative showed the highest potency against cruzain (IC₅₀ = 2.4 μM). Molecular docking study showed that this compound can interact with subsites S1 and S2 simultaneously, and the negative values for the theoretical energy binding (E_b = −7.39 kcal·mol⁻¹) indicates interaction (via dipole–dipole) between the hybridized sulfur sp³ atom at the thiazolidine ring and Gly66. Finally, the results suggest that the thiophen-2-iminothiazolidines synthesized are important lead compounds for the continuing battle against Chagas disease.     

Structure of the O-antigen and characterization of the O-antigen gene cluster of Escherichia coli O108 containing 5,7-diacetamido-3,5,7,9-tetradeoxy-l-glycero-d-galacto-non-2-ulosonic (8-epilegionaminic) acid

On mild acid degradation of the lipopolysaccharide of Escherichia coli O108, the O-polysaccharide was isolated and studied by sugar analysis and one- and two-dimensional ¹H- and ¹³C-NMR spectroscopy. The polysaccharide was found to contain an unusual higher sugar, 5,7-diacetamido-3,5,7,9-tetradeoxy-l-glycero-d-galacto-non-2-ulosonic acid (di-N-acetyl-8-epilegionaminic acid, 8eLeg5Ac7Ac). The following structure of the tetrasaccharide repeating unit of the polysac-charide was established: →4)-α-8eLegp5Ac7Ac-(2→6)-α-D-Galp-(1→3)-α-L-FucpNAc-(1→3)-α-D-GlcpNAc-(1→. Functions of the E. coli O108 antigen biosynthetic genes, including seven putative genes for synthesis of 8eLeg5Ac7Ac, were assigned by sequencing the O-antigen gene cluster along with comparison with gene databases and known biosynthetic pathways for related nonulosonic acids.     

Consequence of one-electron oxidation and one-electron reduction for aniline

Quantum-chemical calculations were performed for all possible isomers of neutral aniline and its redox forms, and intramolecular proton-transfer (prototropy) accompanied by π-electron delocalization was analyzed. One-electron oxidation (PhNH₂ – e → [PhNH₂]^+•) has no important effect on tautomeric preferences. The enamine tautomer is preferred for oxidized aniline similarly as for the neutral molecule. Dramatical changes take place when proceeding from neutral to reduced aniline. One-electron reduction (PhNH₂ + e → [PhNH₂]^-•) favors the imine tautomer. Independently on the state of oxidation, π- and n-electrons are more delocalized for the enamine than imine tautomers. The change of the tautomeric preferences for reduced aniline may partially explain the origin of the CH tautomers for reduced nucleobases (cytosine, adenine, and guanine).     

Molecular modeling study of isoindolines as L-type Ca2+ channel blockers by docking calculations

Two series of isoindolines 1(a–g) and 2(a–g) were found by docking calculations to be possible L-type Ca²⁺ channel (LCC) blockers. The theoretical 3-D model of the outer vestibule and the selective filter of the LCC was provided by Professor Lipkind; this model consists of transmembrane segments S5 and S6 and P-loops contributed by each of four repeats (I, II, III, and IV) of Ca_v 1.2. Therefore, two well-known LCC blockers, nifedipine 3 and (R)-ethosuccinimide 4 were also evaluated, and their binding sites on the LCC were identified and compared with those obtained for 1(a–g) and 2(a–g). Analysis of the results shows that the target compounds tested probably could be LCC blockers, since they interact with or near the glutamic acid residues Glu393, Glu736, Glu1145 and Glu1446 (the EEEE locus), which belong to the LCC selectivity region. The ∆G values for all of the Ca²⁺ channel ligands are between−10.78 and −3.67 (kcal mol⁻¹), showing that LCC-1b, -1e and -1f complexes are more stable than the other compounds tested. Therefore, theoretically calculated dissociation constants K _d (μM) were obtained for all compounds. Comparing these values reveals that compounds 1b (0.0244 μM), 1e (0.0176 μM) and 1f (0.0125 μM) exhibit more affinity for the LCC than the other compounds. This screening shows that the two series of isoindolines probably could act as LCC blockers.     

Teratocarcinoma transplantation rejection loci: Genetic localization of the Gt-1 locus on chromosome 17 and the expression of alternate alleles

A series of congenic mice on the BALB/c genetic background have been employed to localize a teratocarcinoma transplantation rejection locus, Gt-1, to the K side of the H-2 locus on chromosome 17. Previous studies have placed the Gt-1 ^sv allele about 8 centimorgans away from the H-2 ^b or H-2 ^bv1 locus. Teratocarcinomas derived from 129/sv mice, Gt-1 ^sv (H-2K ^bv1/H-2D ^bv1), are rejected by BALB/c (H-2K ^d/H-2D^d) and BALB-G mice (H-2K ^d/H2D-^b, but form tumors in BALB-B (H-2K ^b/H2D ^b) and BALB/5R5 mice (H-2K ^b/H2D ^d). In the reciprocal tumor-rejection test, a BALB/c teratocarcinoma was rejected by immunized BALB·B mice, but formed tumors in the immunized isogenic BALB/c mouse. These studies demonstrate the reciprocal expression of two Gt-1 alleles, one Gt-1 ^c, in BALB/c mice, and the other, Gt-1 ^sv, in the congenic BALB·B mice. Shedlovsky and co-workers have placed the Gt-1 locus in a similar location on the K side of the H-2 locus on chromosome 17.     

Host-determined T cell fine specificity for self-H-2 in radiation bone-marrow chimeras of standard C57BL/6 (H-2b) mutantHz1 (H-2ba), and F1 mice

Reciprocal radiation bone-marrow chimeras were produced between the standard C57BL/6 (=B6) and the mutant B6.C-H-2 ^ba (=Hz1) strain. When infected with vaccinia virus, these chimeras, as well as an (Hz1 × B6)= Hz1 chimera, produced cytotoxic cells that killed vaccinia-infected H-2K^kH-2D^b target cells but failed to kill virus-infected H-2K^bH-2D^d cells. Virus-infected (Hz1 × B6)F₁ B6 chimeras, however, killed both types of target. These experiments demonstrate strict T-cell specificity capable of differentiating between two molecules that apparently differ by a single amino acid substitution.     

Inhibition by light of CO2 evolution from dark respiration: Comparison of two gas exchange methods

Two approaches to determine the fraction (μ) of mitochondrial respiration sustained during illumination by measuring CO₂ gas exchange are compared. In single leaves, the respiration rate in the light (`day respiration' rate R<sub>d</sub>) is determined as the ordinate of the intersection point of A–c<sub>i</sub> curves at various photon flux densities and compared with the CO<sub>2</sub> evolution rate in darkness (`night respiration' rate R_n). Alternatively, using leaves with varying values of CO₂ compensation concentration (Γ), intracellular resistance (r_i) and R_n, an average number for μ can be derived from the linear regression between Γ and the product r_iċR_n. Both methods also result in a number c* for that intercellular CO₂ concentration at which net CO₂ uptake rate is equal to –R_d. c* is an approximate value of the photocompensation point Γ* (Γ in the absence of mitochondrial respiration), which is related to the CO₂/O₂ specificity factor of Rubisco S_c/o. The presuppositions and limitations for application of both approaches are discussed. In leaves of Nicotiana tabacum, at 22 °C, single leaf measurements resulted in mean values of μ = 0.71 and c_* = 34 μmol mol⁻¹. At the photosynthetically active photon flux density of 960 μmol quanta m⁻² s⁻¹, nearly the same numbers were derived from the linear relationship between Γ and r_iċR_n. c* and R_d determined by single leaf measurements varied between 31 and 41 μmol mol⁻¹ and between 0.37 and 1.22 μmol m⁻² s⁻¹, respectively. A highly significant negative correlation between c* and R_d was found. From the regression equation we obtained estimates for Γ* (39 μmol mol⁻¹), S_c/o (96.5 mol mol⁻¹) and the mesophyll CO₂ transfer resistance (7.0 mol⁻¹ m² s). This revised version was published online in June 2006 with corrections to the Cover Date.     

Microbial transformation of ginsenoside Rb<Subscript>1</Subscript> by <Emphasis Type="Italic">Acremonium strictum</Emphasis>

Preparative-scale fermentation of ginsenoside Rb₁ (1) with Acremonium strictum AS 3.2058 gave three new compounds, 12β-hydroxydammar-3-one-20 (S)-O-β-d-glucopyranoside (7), 12β, 25-dihydroxydammar-(E)-20(22)-ene-3-O-β-d-glucopyranosyl-(1→2)-β-d-glucopyranoside (8), and 12β, 20 (R), 25-trihydroxydammar-3-O-β-d-glucopyranosyl-(1→2)-β-d-glucopyranoside (9), along with five known compounds, ginsenoside Rd (2), gypenoside XVII (3), ginsenoside Rg₃ (4), ginsenoside F₂ (5), and compound K (6). The structural elucidation of these metabolites was based primarily on one- and two-dimensional nuclear magnetic resonance and high-resolution electron spray ionization mass spectra analyses. Among these compounds, 2–6 are also the metabolites of ginsenoside Rb₁ in mammals. This result demonstrated that microbial culture parallels mammalian metabolism; therefore, A. strictum might be a useful tool for generating mammalian metabolites of related analogs of ginsenosides for complete structural identification and for further use in pharmaceutical research in this series of compounds. In addition, the biotransformation kinetics was also investigated.     

Homology modeling and dynamics of the extracellular domain of rat and human neuronal nicotinic acetylcholine receptor subtypes α4β2 and α7

In recent years, it has become clear that the neuronal nicotinic acetylcholine receptor (nAChR) is a valid target in the treatment of a variety of diseases, including Alzheimer’s disease, anxiety, and nicotine addiction. As with most membrane proteins, information on the three-dimensional (3D) structure of nAChR is limited to data from electron microscopy, at a resolution that makes the application of structure-based design approaches to develop specific ligands difficult. Based on a high-resolution crystal structure of AChBP, homology models of the extracellular domain of the neuronal rat and human nAChR subtypes α4β2 and α7 (the subtypes most abundant in brain) were built, and their stability assessed with molecular dynamics (MD). All models built showed conformational stability over time, confirming the quality of the starting 3D model. Lipophilicity and electrostatic potential studies performed on the rat and human α4β2 and α7 nicotinic models were compared to AChBP, revealing the importance of the hydrophobic aromatic pocket and the critical role of the α-subunit Trp—the homolog of AChBP-Trp 143—for ligand binding. The models presented provide a valuable framework for the structure-based design of specific α4β2 nAChR subtype ligands aimed at improving therapeutic and diagnostic applications. Figure Electrostatic surface potential of the binding site cavity of the neuronal nicotinic acetylcholine receptor (nAChR). Nicotinic models performed with the MOLCAD program: a rat α7, b rat α4β2, c human α7, d human α4β2. All residues labeled are part of the α7 (a,c) or α4 (b,d) subunit with the exception of Phe 117, which belongs to subunit β2 (d). Violet Very negative, blue negative, yellow neutral, red very positive