Use of calculated cation-pi binding energies to predict relative strengths of nicotinic acetylcholine receptor agonists

300-fold higher neuro- and immunotoxicity from low-redox transformation of carbamazepine

Current challenges in (eco)toxicology are in understanding the transformation of (reactive) substances, and how transformation affects toxic modes of action. Empirical assessment of transformation products of, practically an infinite number of substances, via experimentation, is impossible. Predicting transformation products for (benchmarking) compounds from conditions, facilitates risk analyses. This study applied calculus to predict transformation products of an important environmental and medicinal/toxicological marker, carbamazepine. As radicals are ubiquitous in humans and the environment, we looked into radical-mediated transformations of carbamazepine as a benchmark. We calculated proportions of their speciation states as function of redox conditions, which we took as pH and O2 concentration, describing transformation via covalent and ionic interactions. Formation of ring-contracted products with neuro-immunological activity is thermodynamically favored under anaerobic conditions and at low pH. Experimentally observed product distributions and toxicities reflect that pattern. Our predictive method may support toxicity predictions for other substances and conditions 'similar' to the current case study via interpolation. This paves the way for a more coherent, effective and easier risk assessment of transformation products.

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.

Purinergic control of lysenin's transport and voltage-gating properties

Lysenin, a pore-forming protein extracted from the coelomic fluid of the earthworm Eisenia foetida, manifests cytolytic activity by inserting large conductance pores in host membranes containing sphingomyelin. In the present study, we found that adenosine phosphates control the biological activity of lysenin channels inserted into planar lipid membranes with respect to their macroscopic conductance and voltage-induced gating. Addition of ATP, ADP, or AMP decreased the macroscopic conductance of lysenin channels in a concentration-dependent manner, with ATP being the most potent inhibitor and AMP the least. ATP removal from the bulk solutions by buffer exchange quickly reinstated the macroscopic conductance and demonstrated reversibility. Single-channel experiments pointed to an inhibition mechanism that most probably relies on electrostatic binding and partial occlusion of the channel-conducting pathway, rather than ligand gating induced by the highly charged phosphates. The Hill analysis of the changes in macroscopic conduction as a function of the inhibitor concentration suggested cooperative binding as descriptive of the inhibition process. Ionic screening significantly reduced the ATP inhibitory efficacy, in support of the electrostatic binding hypothesis. In addition to conductance modulation, purinergic control over the biological activity of lysenin channels has also been observed to manifest as changes of the voltage-induced gating profile. Our analysis strongly suggests that not only the inhibitor's charge but also its ability to adopt a folded conformation may explain the differences in the observed influence of ATP, ADP, and AMP on lysenin's biological activity.

Energy for wild-type acetylcholine receptor channel gating from different choline derivatives

Agonists, including the neurotransmitter acetylcholine (ACh), bind at two sites in the neuromuscular ACh receptor channel (AChR) to promote a reversible, global change in protein conformation that regulates the flow of ions across the muscle cell membrane. In the synaptic cleft, ACh is hydrolyzed to acetate and choline. Replacement of the transmitter's ester acetyl group with a hydroxyl (ACh→choline) results in a + 1.8 kcal/mol reduction in the energy for gating generated by each agonist molecule from a low- to high-affinity change of the transmitter binding site (ΔG(B)). To understand the distinct actions of structurally related agonist molecules, we measured ΔG(B) for 10 related choline derivatives. Replacing the hydroxyl group of choline with different substituents, such as hydrogen, chloride, methyl, or amine, increased the energy for gating (i.e., it made ΔG(B) more negative relative to choline). Extending the ethyl hydroxide tail of choline to propyl and butyl hydroxide also increased this energy. Our findings reveal the amount of energy that is available for the AChR conformational change provided by different, structurally related agonists. We speculate that a hydrogen bond between the choline hydroxyl and the backbone carbonyl of αW149 positions this agonist's quaternary ammonium group so as to reduce the cation-π interaction between this moiety and the aromatic groups at the binding site.

Insights into intramolecular Trp and His side-chain orientation and stereospecific π interactions surrounding metal centers: an investigation using protein metal-site mimicry in solution

Metal-binding scaffolds incorporating a Trp/His-paired epitope are instrumental in giving novel insights into the physicochemical basis of functional and mechanistic versatility conferred by the Trp-His interplay at a metal site. Herein, by coupling biometal site mimicry and (1)H and (13)C NMR spectroscopy experiments, modular constructs EDTA-(L-Trp, L-His) (EWH; EDTA=ethylenediamino tetraacetic acid) and DTPA-(L-Trp, L-His) (DWH; DTPA=diethylenetriamino pentaacetic acid) were employed to dissect the static and transient physicochemical properties of hydrophobic/hydrophilic aromatic interactive modes surrounding biometal centers. The binding feature and identities of the stoichiometric metal-bound complexes in solution were investigated by using (1)H and (13)C NMR spectroscopy, which facilitated a cross-validation of the carboxylate, amide oxygen, and tertiary amino groups as the primary ligands and indole as the secondary ligand, with the imidazole (Im) N3 nitrogen being weakly bound to metals such as Ca(2+) owing to a multivalency effect. Surrounding the metal centers, the stereospecific orientation of aromatic rings in the diastereoisomerism is interpreted with the Ca(2+)-EWH complex. With respect to perturbed Trp side-chain rotamer heterogeneity, drastically restricted Trp side-chain flexibility and thus a dynamically constrained rotamer interconversion due to π interactions is evident from the site-selective (13)C NMR spectroscopic signal broadening of the Trp indolyl C3 atom. Furthermore, effects of Trp side-chain fluctuation on indole/Im orientation were the subject of a 2D NMR spectroscopy study by using the Ca(2+)-bound state; a C-H2(indolyl)/C-H5(Im(+)) connectivity observed in the NOESY spectra captured direct evidence that the N-H1 of the Ca(2+)-Im(+) unit interacted with the pyrrole ring of the indole unit in Ca(2+)-bound EWH but not in DWH, which is assignable to a moderately static, anomalous, T-shaped, interplanar π(+)-π stacking alignment. Nevertheless, a comparative (13)C NMR spectroscopy study of the two homologous scaffolds revealed that the overall response of the indole unit arises predominantly from global attractions between the indole ring and the entire positively charged first coordination sphere. The study thus demonstrates the coordination-sphere/geometry dependence of the Trp/His side-chain interplay, and established that π interactions allow (13)C NMR spectroscopy to offer a new window for investigating Trp rotamer heterogeneity near metal-binding centers.